<?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.2018.92011</article-id><article-id pub-id-type="publisher-id">JMP-81875</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>
 
 
  Junction Surface Recombination Concept as Applied to Silicon Solar Cell Maximum Power Point Determination Using Matlab/Simulink: Effect of Temperature
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Bakary</surname><given-names>Dit Dembo Sylla</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>Ibrahima</surname><given-names>Ly</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>Ousmane</surname><given-names>Sow</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Babou</surname><given-names>Dione</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>Youssou</surname><given-names>Traore</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>Grégoire</surname><given-names>Sissoko</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Electromecanical department, Polytechnic School of Thiès, Thiès, Senegal</addr-line></aff><aff id="aff3"><addr-line>University Institute of Technology, University of Thiès, Thiès, Senegal</addr-line></aff><aff id="aff1"><addr-line>Laboratory of Semiconductors and Solar Energy, Physics Department, Faculty of Science and Technology, University Cheikh Anta Diop, Dakar, Senegal</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>gsissoko@yahoo.com(GS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>18</day><month>01</month><year>2018</year></pub-date><volume>09</volume><issue>02</issue><fpage>172</fpage><lpage>188</lpage><history><date date-type="received"><day>15,</day>	<month>December</month>	<year>2017</year></date><date date-type="rev-recd"><day>16,</day>	<month>January</month>	<year>2018</year>	</date><date date-type="accepted"><day>19,</day>	<month>January</month>	<year>2018</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>
 
 
  In this work, we study the method for determining the maximum of the minority carrier recombination velocity at the junction Sf
  <sub>max</sub>, corresponding to the maximum power delivered by the photovoltaic generator. For this, we study the temperature influence on the behavior of the front white biased solar cell in steady state. By solving the continuity equation of excess minority carrier in the base, we have established the expressions of the photocurrent density, the recombination velocity on the back side of the base Sb, and the photovoltage. The photocurrent density and the photovoltage are plotted as a function of Sf, called, minority carrier recombination velocity at the junction surface, for different temperature values. The illuminated I-V characteristic curves of the solar cell are then derived. To better characterize the solar cell, we study the electrical power delivered by the base of the solar cell to the external charge circuit as either junction surface recombination velocity or photovoltage dependent. From the output power versus junction surface recombination velocity Sf, we have deduced an eigenvalue equation depending on junction recombination velocity. This equation allows to obtain the maximum junction recombination velocity Sf
  <sub>max</sub> corresponding to the maximum power delivered by the photovoltaic generator, throughout simulink model. Finally, we deduce the conversion efficiency of the solar cell.
 
</p></abstract><kwd-group><kwd>Silicon Solar Cell-Junction Surface Recombination Velocity-Maximum Power</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Many techniques exist and provide photovoltaic generators to operate at maximum points of their characteristics [<xref ref-type="bibr" rid="scirp.81875-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref2">2</xref>] . The method used in this work enables us to determine the maximum power point based on electronic transport parameters of the solar cell. Contrary to the most commonly used MPPT technical, which takes into account the controllers particularly adapted for managing a non-linear source by forcing the generator to work at its maximum power point (MPP) [<xref ref-type="bibr" rid="scirp.81875-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref5">5</xref>] .</p><p>This study aims to determination of the maximum power point supplied by a photovoltaic generator in static regime under the effect of temperature. For this, we give through the continuity equation, the expressions of minority carrier generation rate, and the density of the excess minority carrier in the base [<xref ref-type="bibr" rid="scirp.81875-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref7">7</xref>] .</p><p>The expressions of the photocurrent density, the photovoltage, the recombination velocity of the minority carrier at the back surface Sb and the electrical power, all depending on temperature are produced and graphically represented as function of minority carrier recombination velocity at the junction [<xref ref-type="bibr" rid="scirp.81875-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref10">10</xref>] .</p><p>The I(Sf)-V(Sf) characteristic curves of the photocurrent density as function of photovoltage, the power as a function of minority carrier recombination velocity at the junction or photovoltage dependent, have also graphically represented [<xref ref-type="bibr" rid="scirp.81875-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref13">13</xref>] .</p><p>A transcendental equation giving the minority carrier recombination velocity to the maximum power points Sf<sub>max</sub> is determined. It is graphically represented by matlab/Simulink as a minority carrier recombination velocity at the junction Sf (Time (secs)). Finally, the graphical results are compared with those of the model [<xref ref-type="bibr" rid="scirp.81875-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref16">16</xref>] .</p></sec><sec id="s2"><title>2. Theory</title><p>Consider a crystalline silicon solar cell (n<sup>+</sup>-p-p<sup>+</sup>). Its structure is illustrated in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>z is the depth in the base of the solar cell measured from the junction emitter- base (z = 0) to the back surface (z = H). H is the base thickness.</p><p>For a given illumination of the solar cell, different processes take place in the base, it is the absorption of the light, the minority carrier generation, bulk and surfaces recombination and diffusion. These processes can be simplified to the one dimensional continuity equation with boundary conditions:</p><p>∂ 2 δ ( z , T ) ∂ z 2 − δ ( z , T ) L 2 ( T ) = − 1 D ( T ) ⋅ G ( z ) (1)</p><p>δ ( z , T ) represents the excess minority carrier density in the solar cell base at position z.</p><p>With</p><p>L ( T ) = τ ⋅ D ( T ) (2)</p><p>L(T) is the minority carrier diffusion length in the base and is function of the temperature. It also represents the average distance traveled during the lifetime (τ) by the minority carrier before their recombination.</p><p>D(T) is the electron diffusion coefficient in the base given by the well-known Einstein relation [<xref ref-type="bibr" rid="scirp.81875-ref18">18</xref>] , temperature dependent, given as:</p><p>D ( T ) = μ ( T ) K b q T (3)</p><p>μ(T) is the mobility coefficient for electron and depends on the temperature [<xref ref-type="bibr" rid="scirp.81875-ref19">19</xref>] , its expression is given by the following equation:</p><p>μ ( T ) = 1.43 &#215; 10 9 T − 2.42   cm 2 ⋅ V − 1 ⋅ s − 1 (4)</p><p>K<sub>b</sub> is the Boltzmann constant, q is the elementary charge of an electron.</p><p>G(z) is the minority carriers generation rate at position z in the base, given by [<xref ref-type="bibr" rid="scirp.81875-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref21">21</xref>] .</p><p>G ( z ) = ∑ i = 1 3 a i e − b i ⋅ z (5)</p><p>The coefficient a<sub>i</sub> and b<sub>i</sub> are obtained from tabulated values of the radiation in A.M 1.5 condition [<xref ref-type="bibr" rid="scirp.81875-ref22">22</xref>] .</p><p>The excess minority carrier density is obtained from differential Equation (1) resolution and is given by.</p><p>δ ( z , T ) = A ⋅ cosh ( z L ( T ) ) + B ⋅ sinh ( z L ( T ) ) − ∑ i = 1 3 K i ⋅ e − b i ⋅ z (6)</p><p>The coefficients A and B are obtained with the boundary conditions at the emitter-base junction and at the back surface of the cell [<xref ref-type="bibr" rid="scirp.81875-ref22">22</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref23">23</xref>] :</p><p>i) At the junction: emitter-base (z = 0)</p><p>∂ δ ( z , T ) ∂ z | z = 0 = S f D ( T ) δ ( z = 0 , T ) | z = 0 (7)</p><p>ii) At the back side (z = H)</p><p>∂ δ ( z , T ) ∂ z | z = H = − S b D ( T ) δ ( z = H , T ) | z = H (8)</p><p>Sf is the excess minority carrier recombination velocity at the junction emitter-base, it also characterizes the operating point of the solar cell [<xref ref-type="bibr" rid="scirp.81875-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref25">25</xref>] on the illuminated current-voltage characteristic.</p><p>Sb is the excess minority carrier recombination velocity on the back of the base [<xref ref-type="bibr" rid="scirp.81875-ref25">25</xref>] . Its expression is obtained from the derivative (<xref ref-type="fig" rid="fig2">Figure 2</xref>) of the photocurrent density for large Sf values [<xref ref-type="bibr" rid="scirp.81875-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref20">20</xref>] .</p><p>[ ∂ J p h ∂ S f ] = 0 (9)</p><p>From the Equation (9), Sb is then expressed as dependent on the following parameters L(T), D(T), μ(T), b<sub>i</sub> and H:</p><p>S b ( T ) = D ( T ) L ( T ) ⋅ ∑ i = 1 3 L ( T ) ⋅ b i ( e b i ⋅ H − cosh ( H L ( T ) ) ) − sinh ( H L ( T ) ) L ( T ) ⋅ b i ⋅ sinh ( H L ( T ) ) + cosh ( H L ( T ) ) − e b i ⋅ H (10)</p><sec id="s2_1"><title>2.1. Photocurrent Density</title><p>The expression of the photocurrent density is obtained from the density of the minority charge carriers in the base. It’s given by the following relation:</p><p>J p h ( S f , T ) = q ⋅ D ( T ) ⋅ [ ∂ δ ( z , S f , T ) ∂ z ] z = 0 (11)</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref> represents the photocurrent density profile as function of the junction</p><p>recombination velocity for different temperature values.</p><p>On this curve, the photocurrent density is very low in the vicinity of the open-circuit (Sf &lt; 2 &#180; 10<sup>2</sup> cm/s), then increases very rapidly with the junction minority carrier recombination velocity to reach an asymptotic value at large Sf values (Sf &gt; 5 &#180; 10<sup>5</sup> cm/s) corresponding to solar cell short-circuit photocurrent density. The short circuit current density decreases as the temperature increases. This decreasing due to the disordered movement of the charge carriers because the increasing of temperature causes a thermal agitation.</p></sec><sec id="s2_2"><title>2.2. Photovoltage</title><p>The expression of the photovoltage at the terminal of the solar cell, when the latter is subjected to a multispectral illumination, is obtained by the Boltzmann relation [<xref ref-type="bibr" rid="scirp.81875-ref26">26</xref>] .</p><p>V p h ( S f , T ) = V T ⋅ log [ N b n i 2 ( T ) δ ( S f , T ) + 1 ] (12)</p><p>V<sub>T</sub> is the thermal voltage, it is given as follows:</p><p>V T = K b q T (13)</p><p>T is the absolute temperature, it’s included between 300 - 350 K.</p><p>N<sub>b</sub> is the acceptor atom doping rate in the base.</p><p>n i 2 ( T ) = A ⋅ T 3 ⋅ exp ( − E g K b ⋅ T ) (14)</p><p>n i ( T ) is the law of conservation (generation rate must be equal to recombination rates of charge carriers n = p = n<sub>i</sub>) [<xref ref-type="bibr" rid="scirp.81875-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref28">28</xref>] .</p><p>E<sub>g</sub> is the energy gap, it corresponds to the difference between the energy of the conduction band Ec and the valence band E<sub>g</sub>. E<sub>g</sub> = 1.12 &#180; 1.6 &#180; 10<sup>−19</sup> J,</p><p>A is a specific constant of the material, A = 3.87 &#180; 10<sup>16</sup> cm<sup>−3</sup>∙K<sup>−3/2</sup> [<xref ref-type="bibr" rid="scirp.81875-ref18">18</xref>] .</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> represents the photovoltage as function of junction recombination velocity for different temperature values.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> shows that, at low junction minority carrier recombination velocity values, the photovoltage is maximal and constant, and corresponds to the open-circuit voltage (vicinity of the open-circuit). There is a blockage and storage of minority carriers at the junction. The phototvoltage decreases for large junction minority carrier recombination velocity values and it becomes low in the vicinity of the short-circuit (large Sf values), thus the minority carriers crossed over the junction and participate to the production of photocurrent. The increasing of temperature leads to a slight increase in the creation of electron - hole pairs.</p></sec><sec id="s2_3"><title>2.3. Study of the I(Sf)-V(Sf) Characteristic</title><p>The profile of the I(Sf)-V(Sf) characteristic for different temperature values is represented on <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><p>We notice that, the short-circuit photocurrent density decreases when the temperature increases and the open circuit voltage increases with increasing of temperature.</p></sec><sec id="s2_4"><title>2.4. Study of the Power and the Maximum Power Point</title><sec id="s2_4_1"><title>2.4.1. Study of the Power</title><p>The equivalent electrical circuit of an actual solar cell under illumination is represented by <xref ref-type="fig" rid="fig5">Figure 5</xref>. This circuit treats the solar cell as an ideal current generator which sells a photocurrent density depending on the illumination, connected in parallel with a diode and a shunt resistance R<sub>sh</sub> and in series with a series resistance Rs.</p><p>The Ohm law applied to <xref ref-type="fig" rid="fig5">Figure 5</xref>, allows to extract the solar cell base electrical power issued to the external charge circuit as following (15):</p><p>P ( S f , T ) = V p h ( S f , T ) ⋅ I ( S f , T ) (15)</p><p>Applying the first law of Kirchhoff to the circuit of <xref ref-type="fig" rid="fig5">Figure 5</xref>, the current density supplied to the external charge circuit is given by the following relation [<xref ref-type="bibr" rid="scirp.81875-ref16">16</xref>] .</p><p>I ( S f , T ) = I p h ( S f , T ) − I d ( S f , T ) − I S h ( S f , T ) (16)</p><p>I<sub>d</sub> is the diode current density, its expression is given by the following relation:</p><p>I d ( S f , T ) = q ⋅ S f 0 ⋅ ( n i ( T ) ) 2 N b ⋅ exp ( V p h ( S f , T ) V T − 1 ) (17)</p><p>I<sub>sh</sub> is the shunt current density, in the case of an ideal current generator, (R<sub>s</sub><sub>h</sub> tends to infinity; I<sub>Sh</sub> = 0).</p><p>Sf<sub>0</sub> is the intrinsic junction recombination velocity associated with the losses of charge carriers induced by shunt resistance. It characterizes the good quality of the solar cell [<xref ref-type="bibr" rid="scirp.81875-ref29">29</xref>] .</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref> represent the variations of the electrical power in function of the recombination velocity of the charge carriers at the junction and the photovoltage for different temperature values.</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> shows that, the power increases with the junction surface recombination velocity until a value Sf<sub>max</sub> then it decreases. The maximum power point varies with the temperature. Increasing photocurrent density is due to a decreasing of the energy larger band gap E<sub>g</sub>. At the same time, we observe an increasing of the diode current density resulting a decreasing values of temperature.</p></sec><sec id="s2_4_2"><title>2.4.2. Maximum Power Point and Efficiency</title><p>The maximum power point of a generator photovoltaic corresponds to the</p><p>photocurrent density-photovoltage couple generating the maximum electrical power of this solar cell.</p><p>Four essential data can be used to determine the photocurrent density-photovoltage characteristic of the photovoltaic generator [<xref ref-type="bibr" rid="scirp.81875-ref17">17</xref>] .</p><p>・ The photocurrent density of short-circuit noted J<sub>sc</sub></p><p>・ The open-circuit voltage noted V<sub>oc</sub></p><p>・ The maximal photocurrent density noted J<sub>ph</sub><sub>max</sub></p><p>・ The maximal photovoltage noted V<sub>ph</sub><sub>max</sub></p><p>The product of the maximal photocurrent density and the maximal photovoltage J<sub>ph</sub><sub>max</sub> &#215; V<sub>ph</sub><sub>max</sub> gives a maximal power P<sub>max</sub>.</p><p>We determine the maximal minority carrier recombination velocity at the junction corresponding to the maximal power point [<xref ref-type="bibr" rid="scirp.81875-ref30">30</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref31">31</xref>] . Let:</p><p>∂ P ∂ S f = 0 (18)</p><p>Let Sf<sub>max</sub> denote this maximal junction recombination velocity corresponding to the maximal power point. It depends on:</p><p>i) electronic transport parameters (L, μ, D, Sf, Sb, τ, n<sub>i</sub>, N<sub>b</sub>) in the solar cell that are related to physical phenomena,</p><p>ii) geometrical one (H), i.e. in the 1D model</p><p>iii) material absorption coefficients (b<sub>i</sub>).</p><p>All these parameters are taking into account for the determination of the maximum power based on electron diffusion coefficient variation with temperature and the solar cell thickness H [<xref ref-type="bibr" rid="scirp.81875-ref18">18</xref>] .</p><p>From the Equation (18), the calculation gives the transcendental equation dependent of the junction minority carrier recombination velocity.</p><p>and the temperature. It’s given by the following expressions:</p><p>With</p><p>M ( S f , T ) = 1 S f max L [ 1 − S f max L Y 1 D + S f max L ] (19)</p><p>And</p><p>N ( S f , T ) = [ Γ max ( 0 , T ) ( Γ max ( 0 , T ) + n i 2 N b ) ⋅ ( S f max ⋅ L + Y 1 ⋅ D ) ] * [ 1 log ( N b ⋅ Γ max ( 0 , T ) n i 2 + 1 ) ] (20)</p><p>Γ max ( 0 , T ) is the density of the minority carrier on maximum power point, its expression is given by the following relation:</p><p>Γ max ( 0 , T ) = K ⋅ D ⋅ [ Y 2 + Y 1 − b i ⋅ L S f max ⋅ L + Y 1 ⋅ D ] (21)</p><p>With</p><p>K = n ⋅ a i ⋅ L 2 ⋅ cos ( θ ) D ⋅ ( L 2 ⋅ b i 2 − 1 ) (22)</p><p>Y 1 = D / L ⋅ sinh ( H / L ) + S b ⋅ cosh ( H / L ) D / L ⋅ cosh ( H / L ) + S b ⋅ sinh ( H / L ) (23)</p><p>Y 2 = ( D ⋅ b i − S b ) ⋅ exp ( − b i ⋅ H ) D / L ⋅ cosh ( H / L ) + S b ⋅ sinh ( H / L ) (24)</p><p>For the transcendental equation modelling a simulink model is used [<xref ref-type="bibr" rid="scirp.81875-ref32">32</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref33">33</xref>] . This model is composed of an inputs block that includes a system allowing to define the set simulation type, simulation parameters, and preferences (powergui), the ramp block, which output a ramp signal starting at the specified time with the start time and initial output estimated at zero. To file block, allows the variable may be created as a matlab time series, an array, or a matlab structure. We used matlab time series which is used for any data type. It also allows defining the decimation and the sample time (−1 for inherited). The constant block containing the different values of temperature. The Sf<sub>max</sub> block containing the expression of the transcendental equation. Finally, the output block composed of a creator bus, this block creates a bus signal from its inputs, it also allows defining the inherit bus signal names from input ports, the number of inputs, the signal in the bus and the output data type. The time scope block, this block defined the number of input ports (two in this model), input processing: elements as channels (sample based), time span of Sf<sub>max</sub> block (10 seconds), the time units in metric (based on time span), the time displays ofset at zero, show time axis-labels, the minimum and maximum limit of y-axis and y-axis label.</p><p>The model of the transcendental equation is represented in <xref ref-type="fig" rid="fig8">Figure 8</xref> as follows:</p><p>The graphical resolution of this model provides the values of Sf<sub>max</sub> defined by the intercept point of two curves represented on <xref ref-type="fig" rid="fig9">Figure 9</xref>. To each operating point generating the maximum electrical power delivered by the photovoltaic generator corresponds to a Sf<sub>max</sub> value.</p><p><xref ref-type="fig" rid="fig9">Figure 9</xref> shows the increasing of the Sf<sub>max</sub> values then the temperature increases. Reflecting the increasing of the maximal power as the temperature increases.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The numerical values of Sf<sub>max</sub> corresponding to the maximal power point for different temperature values</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Temperature (K)</th><th align="center" valign="middle" >Points of intersection for each given temperature values p (cm/s)</th><th align="center" valign="middle" >Sf<sub>max</sub> (cm/s)</th></tr></thead><tr><td align="center" valign="middle" >308 K</td><td align="center" valign="middle" >3.396</td><td align="center" valign="middle" >2.443 &#215; 10<sup>3</sup></td></tr><tr><td align="center" valign="middle" >318 K</td><td align="center" valign="middle" >3.883</td><td align="center" valign="middle" >7.620 &#215; 10<sup>3</sup></td></tr><tr><td align="center" valign="middle" >328 K</td><td align="center" valign="middle" >4.370</td><td align="center" valign="middle" >23.659 &#215; 10<sup>3</sup></td></tr><tr><td align="center" valign="middle" >338 K</td><td align="center" valign="middle" >4.913</td><td align="center" valign="middle" >80.909 &#215; 10<sup>3</sup></td></tr><tr><td align="center" valign="middle" >349 K</td><td align="center" valign="middle" >5.501</td><td align="center" valign="middle" >31.686 &#215; 10<sup>4</sup></td></tr></tbody></table></table-wrap><p>We observe the intercept points on the figure corresponding to the Sf<sub>max</sub> values. These Sf<sub>max</sub> values correspond to an operating condition of the solar cell at the maximum power point.</p><p>The results obtained by the model on <xref ref-type="fig" rid="fig9">Figure 9</xref> corresponding the numerical values of Sf<sub>max</sub> for each maximal power point are given in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>The influence of Sf<sub>max</sub> on the temperature can be plotted from a theoretical model of Sf<sub>max</sub> given by the following expression (25):</p><p>S f max = b ⋅ e a T (25)</p><p>The resolving of the right equation from <xref ref-type="fig" rid="fig1">Figure 1</xref>0 provides the coefficients a and b values follows as: a = 0.052 cm∙s<sup>−1</sup>/K, this value represents the slope corresponding to the vertical variation versus the horizontal variation of the right and b = 3.216 cm/s corresponding a temperature of 305 K, it is determined from the y intercepts and it is homogeneous to a maximal surface recombination velocity Sf<sub>max</sub>. Sf<sub>max</sub> increases with temperature on a low positive slope.</p></sec><sec id="s2_4_3"><title>2.4.3. The Efficiency</title><p>The conversion efficiency of the solar cell is a ratio between the maximal power supplied by the solar cell and the absorbed incident light power, writing as follows:</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> The photocurrent density I<sub>max</sub> values, the photovoltage V<sub>max</sub>, the power Pmax and the efficiency Ƞmax corresponding to the maximal power point for different temperature values</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Temperature (K)</th><th align="center" valign="middle" >308 K</th><th align="center" valign="middle" >318 K</th><th align="center" valign="middle" >328 K</th><th align="center" valign="middle" >338 K</th><th align="center" valign="middle" >349 K</th></tr></thead><tr><td align="center" valign="middle" >I<sub>max</sub> (mA/cm<sup>2</sup>)</td><td align="center" valign="middle" >0.03236</td><td align="center" valign="middle" >0.03132</td><td align="center" valign="middle" >0.03033</td><td align="center" valign="middle" >0.02935</td><td align="center" valign="middle" >0.02836</td></tr><tr><td align="center" valign="middle" >I<sub>d</sub> (Sf<sub>max</sub>) (mA)</td><td align="center" valign="middle" >3.023 &#215; 10<sup>−4</sup></td><td align="center" valign="middle" >2.923 &#215; 10<sup>−4</sup></td><td align="center" valign="middle" >2.83 &#215; 10<sup>−4</sup></td><td align="center" valign="middle" >2.741 &#215; 10<sup>−4</sup></td><td align="center" valign="middle" >2.65 &#215; 10<sup>−4</sup></td></tr><tr><td align="center" valign="middle" >V<sub>max</sub> (mV)</td><td align="center" valign="middle" >0.4891</td><td align="center" valign="middle" >0.5113</td><td align="center" valign="middle" >0.5314</td><td align="center" valign="middle" >0.5514</td><td align="center" valign="middle" >0.5713</td></tr><tr><td align="center" valign="middle" >P<sub>max</sub> (W/cm<sup>2</sup>)</td><td align="center" valign="middle" >0.015679</td><td align="center" valign="middle" >0.015864</td><td align="center" valign="middle" >0.015967</td><td align="center" valign="middle" >0.016032</td><td align="center" valign="middle" >0.016051</td></tr><tr><td align="center" valign="middle" >η<sub>max</sub> (%)</td><td align="center" valign="middle" >15.679</td><td align="center" valign="middle" >15.864</td><td align="center" valign="middle" >15.967</td><td align="center" valign="middle" >16.032</td><td align="center" valign="middle" >16.051</td></tr></tbody></table></table-wrap><p>η = J max ⋅ V max P incident (26)</p><p>P<sub>incident</sub> is the absorbed incident light power by the solar cell, with P<sub>incident</sub> = 100 mW/cm<sup>2</sup> in the standards AM 1.5 condition [<xref ref-type="bibr" rid="scirp.81875-ref22">22</xref>] .</p><p>For representation of efficiency, we have deduced on the characteristic curve of the photocurrent density as function of photovoltage according to <xref ref-type="fig" rid="fig4">Figure 4</xref>, the graphical values corresponding to the maximal power point, the maximal photocurrent density, the maximal photovoltage, the maximal power and the maximal conversion efficiency of the solar cell for different temperature values. These results are given in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>Figures 11-13 represent the maximal photocurrent density I<sub>max</sub>, the maximal photovoltage V<sub>max</sub> and the maximal conversion efficiency η<sub>max</sub> of the solar cell as function of temperature.</p><p>The resolution of line equation y = γT + χ, allows to obtain the line coefficient γ and χ from <xref ref-type="fig" rid="fig1">Figure 1</xref>1, with γ = −9.756 &#215; 10<sup>−5</sup> mA∙cm<sup>−2</sup>/K and χ = 0.03227 mA/cm<sup>−2</sup> corresponding at the temperature of 305 K. On this curve, we notice a decreasing of the maximal photocurrent I<sub>max</sub> when the temperature increases,</p><p>that explains a negative slope of the line. Therefore the decreasing of the photocurrent density as the temperature increases is realized.</p><p>The graphical representation of the maximal photovoltage versus temperature is given by <xref ref-type="fig" rid="fig1">Figure 1</xref>2.</p><p>In <xref ref-type="fig" rid="fig1">Figure 1</xref>2, the resolution of line equation provide the coefficient values γ = 0.002 mV/K and χ = 0.4842 mV for a temperature of 305 K. We notice a low positive slope, it means that the photovoltage is maximal and is equal to the open-circuit voltage. The increasing of the temperature implies a minor increasing of the creation electron-hole pairs.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>3 represents the variation of the maximal photovoltaic efficiency conversion as function of temperature.</p><p>These results shows that, the maximal junction recombination velocity Sf<sub>max</sub>, the maximal photovoltage and the conversion efficiency photovoltaic increase when the temperature increases, contrary the maximal photocurrent density which decreases with increasing of temperature. This shows the variation of the maximum power point as the temperature increases [<xref ref-type="bibr" rid="scirp.81875-ref34">34</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref35">35</xref>] . The method in this work shows us, how to optimize the electrical power supplied by the solar cell, by use of Sf<sub>max</sub> as a research function of the maximal power point. The method is mostly based on the variation of electron diffusion coefficient in the base D, the diffusion length L, the mobility coefficient of the electron μ, the intrinsic concentration of minority carriers in the base (ni) and the recombination velocity of the charge carriers at the back side Sb all in function of temperature with the possibility to vary the thickness of the base, the electron lifetime in the base and the electron doping rate Nb. Its strength is that, its depend on the electronic structuration parameters based on the physical mechanism of solar cell, contrary the converters and different tracking systems forcing the solar cell operating at the maximum power point (MPP) based on the macroscopic parameter of solar cell [<xref ref-type="bibr" rid="scirp.81875-ref36">36</xref>] [<xref ref-type="bibr" rid="scirp.81875-ref37">37</xref>] .</p></sec></sec></sec><sec id="s3"><title>3. Conclusions</title><p>In this work, after resolution the expression of the density of the minority charge carriers in excess in the base, the photocurrent density and the photovoltage, the I-V characteristic is proposed. This study shows us a decreasing of the short- circuit photocurrent and an increasing of the open-circuit photovoltage when the temperature increases. The decreasing of the short-circuit photocurrent manifests by a decreasing of the density of the excess minority carrier which crosses the junction when the temperature increases, which leads a increasing of the open-circuit photovoltage.</p><p>From the i-v characteristic of <xref ref-type="fig" rid="fig4">Figure 4</xref>, we studied the delivered electrical power by the base of the solar cell as the function of the density of minority charge carriers at the junction Sf and the photovoltage (<xref ref-type="fig" rid="fig6">Figure 6</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>). We note that the power increases with the recombination velocity of the charge carriers at the junction Sf and of the photovoltage until a maximal value which represents the maximum power then it decreases and cancel of a values corresponding of the open-circuit photovoltage (large Sf).</p><p>A transcendental equation allowing to obtain the maximal recombination velocity of the charge carriers Sf<sub>max</sub> corresponding to the maximal power point of the solar cell is determined as the function of temperature. It depends only on phenomelogical and geometrical parameters of the solar cell (L, μ, D, Sf, Sb, τ, ni, Nb and H) and the material absorptions coefficients bi, in the determination of the maximum power based on variation of the electron diffusion coefficient in the base as a function of temperature and the thickness H of the solar cell [<xref ref-type="bibr" rid="scirp.81875-ref18">18</xref>] .</p><p>The transcendental equation is modeled from a simulink model (<xref ref-type="fig" rid="fig8">Figure 8</xref>), then graphically represented in <xref ref-type="fig" rid="fig9">Figure 9</xref>, and allows to extract the numerical values of Sf<sub>max</sub> for different values of the temperature noted in <xref ref-type="table" rid="table1">Table 1</xref>. The results in <xref ref-type="table" rid="table1">Table 1</xref> show that the recombination velocity of the charge carriers at the junction corresponding to the maximal power point increases, and varies with the temperature (<xref ref-type="fig" rid="fig9">Figure 9</xref>). Finally, we studied the conversion efficiency of the solar cell from the graphical values extracted of the I-V characteristic curve of <xref ref-type="fig" rid="fig4">Figure 4</xref>, then plotted I<sub>max</sub>, V<sub>max</sub> and η<sub>max</sub> (Figures 11-13) as the function of the temperature compared to the Sf<sub>max</sub> values of <xref ref-type="fig" rid="fig1">Figure 1</xref>0 obtained from the simulink model (<xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="fig" rid="fig9">Figure 9</xref>).</p></sec><sec id="s4"><title>Cite this paper</title><p>Sylla, B.D.D., Ly, I., Sow, O., Dione, B., Traore, Y. and Sissoko, G. (2018) Junction Surface Recombination Concept as Applied to Silicon Solar Cell Maximum Power Point Determination Using Matlab/Simulink: Effect of Temperature. Journal of Modern Physics, 9, 172-188. https://doi.org/10.4236/jmp.2018.91011</p></sec></body><back><ref-list><title>References</title><ref id="scirp.81875-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Benmoussa, W.C., Amara, S. and Zerga, A. (2007) Optimisation du rendement d’une photopile. Revue des Energies Renouvelables ICRESD-07 Tlemcen, 301-306.</mixed-citation></ref><ref id="scirp.81875-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Gueye, S., Diallo, H.L., Ndaye, M., Dione, M.M. and Sissoko, G. (2013) International Journal of Emerging Technology and Advanced Engineering (IJTEA), 3, 1-9.</mixed-citation></ref><ref id="scirp.81875-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Barro, F.I., Mbodji, S., Ndiaye, M., Ba, E. and Sissoko, G. (2008) Influence of Grains Size and Grains Boundaries Recombination on the Space-Charge Layer Thickness z of Emitter-Base Junction’s n+-p-p+ Solar Cell. Proceedings of 23rd European Photovoltaic Solar Energy Conference and Exhibition, 604-607.</mixed-citation></ref><ref id="scirp.81875-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Ndoye, S., Ndiaye, M., Diao, A., Dione, M.M., Diarisso, D., Bama, A.O.N., Ly, I., Sow, G., Maiga, A.S., Foulani, A., Barro, F.I. and Sissoko, G. (2010) Modelling and Simuling the Powering System of a Base Transmitter Station with a Standalone Photovoltaic Generator. Proceedings of 25th European Photovoltaic Solar Energy Conference and Exhibition, 5208-5211.</mixed-citation></ref><ref id="scirp.81875-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Diarisso, D., Ly, I., Sow, G., Sow, O., Gaye, I., Barro, F.I. and Sissoko, G. (2002) Research Journal of Applied Sciences, Engineering and Technology, 4, 3740-3745.</mixed-citation></ref><ref id="scirp.81875-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Barro, F.I., Maiga, A.S., Wereme, A. and Sissoko, G. (2010) Phys. Chem. News, 56, 76-84.</mixed-citation></ref><ref id="scirp.81875-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Sissoko, G., Correa, A., Nanema, E., Diarra, M.N., Ndiaye, A.L. and Adj, M. (1998) Recombination Parameters Measurement in Silicon Double Sided Surface Field Cell. Proceeding of the World Renewable Energy Congress, 20-25 September 1998, 1856-1859.</mixed-citation></ref><ref id="scirp.81875-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Gaye, I., Corréa, A., Ba, B., Ndiaye, A.L., Nanéma, E., Ba, A.B.B., Adj, M. and Sissoko, G. (1996) Renewable Energy, 3, 1598-1601.</mixed-citation></ref><ref id="scirp.81875-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Ly, I., Wade, M., Diallo, H.L., EL Moujtaba, M.A.O., Lemtabott, O.H., Mbodji, S., Diasse, O., Ndiaye, A., Gaye, I., Barro, F.I., Wereme, A. and Sissoko, G. (2011) Irradiation Effect on the Electrical Parameters of a Bifacial Silicon Solar Cell under Multispectral Illumination. Proceedings of 26th European Photovoltaic Solar Energy Conference and Exhibition, 785-788.</mixed-citation></ref><ref id="scirp.81875-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Sow, O., Zerbo, I., Mbodji, S., Ngom, M.I., Diouf, M.S. and Sissoko, G. (2012) International Journal of Sciences and Technology, 1, 230-246.</mixed-citation></ref><ref id="scirp.81875-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Corréa, A., Gaye, I., Ba, B., Ndiaye, A.L. and Sissoko, G. (1994) Renewable Energy, 5, 166-168.</mixed-citation></ref><ref id="scirp.81875-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Madougou, S., Made, F., Boukary, M.S. and Sissoko, G. (2007) Advanced Materials Research, 18-19, 303-312.</mixed-citation></ref><ref id="scirp.81875-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Mbodji, S., Ly, I., Diallo, H.L., Dione, M.M., Diasse, O. and Sissoko, G. (2012) Research Journal of Applied Sciences, Engineering and Technology, 4, 1-7.</mixed-citation></ref><ref id="scirp.81875-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Zoungrana, M., Dieng, B., Lemrabott, O.H., Touré, F., Ould EL Moujtaba, M.A., Sow, M.L. and Sissoko, G. (2012) Research Journal of Applied Sciences, Engineering and Technology, 4, 2967-2972.</mixed-citation></ref><ref id="scirp.81875-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Sissoko, G., Nanéma, E., Corréa,A., Biteye, P.M., Adj, M. and Ndiaye, A.L. (1996) Renewable Energy, 3, 1848-1851.</mixed-citation></ref><ref id="scirp.81875-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Mbodji, S., Ly, I., Diallo, H.L., Dione, M.M., Diasseand, O. and Sissoko, G. (2012) Research Journal of Applied Sciences, Engineering and Technology, 4, 1-7.</mixed-citation></ref><ref id="scirp.81875-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Nam, L., Rodot, M., Nijs, J., Ghannam, M. and Coppye, J. (1992) Journal de Physique III, EDP Sciences, 2, 1305-1316.</mixed-citation></ref><ref id="scirp.81875-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Diouf, M., Gaye, S., Thiam, I., Fall, A., M, M.F., Ly, I. and Sissoko, G. (2014) Current Trends in Technology &amp; Sciences, 3, 372-375. http://www.ctts.in/</mixed-citation></ref><ref id="scirp.81875-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Kunst, M. and Sanders, A. (1992) Semiconductor Science and Technology, 7, 51-59.</mixed-citation></ref><ref id="scirp.81875-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Zerbo, I., Koalaga, Z., Barro, F.I., Zougmore, F., Ndiaye, A.L., Diao, A. and Sissoko, G. (2004) Journal des Sciences, 4, 42-46.</mixed-citation></ref><ref id="scirp.81875-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Furlan, J. and Amon, S. (1985) Solid-State Electronics, 28, 1241-1243. https://doi.org/10.1016/0038-1101(85)90048-6</mixed-citation></ref><ref id="scirp.81875-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Mohammad, S.N. (1987) Journal of Applied Physics, 28, 767-772.</mixed-citation></ref><ref id="scirp.81875-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Sissoko, G., Museruka, C., Corréa, A., Gaye, I. and Ndiaye, A.L. (1996) Renewable Energy, 3, 1487-1490.</mixed-citation></ref><ref id="scirp.81875-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Sissoko, G., Sivoththanam, S., Rodot, M. and Mialhe, P. (1992) Constant Illumination-Induced Open Circuit Voltage Decay (CIOCVD) Method, as Applied to High Efficiency Si Solar Cells for Bulk and Back Surface Characterization. 11th European Photovoltaic Solar Energy Conference and Exhibition, Poster 1B, Montreux, 12-16 October 1992, 352-354.</mixed-citation></ref><ref id="scirp.81875-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Diallo, H.L., Maiga, A.S., Wereme, A. and Sissoko, G. (2008) The European Physical Journal Applied Physics, 42, 203-211. https://doi.org/10.1051/epjap:2008085</mixed-citation></ref><ref id="scirp.81875-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Thurmond, C.D. (1975) Journal of the Electrochemical Society, 122, 133-141. https://doi.org/10.1149/1.2134410</mixed-citation></ref><ref id="scirp.81875-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">El-Adawi, M.K. and Al-Nuaim, I.A. (2002) Vacum, 64, 33-36. https://doi.org/10.1016/S0042-207X(01)00370-0</mixed-citation></ref><ref id="scirp.81875-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Peykov, P. and Aceves, M. (2004) Superficies y Vacío, 17, 29-31.</mixed-citation></ref><ref id="scirp.81875-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Thiam, N.D., Diao, A., Ndiaye, M., Dieng, A., Thiam, A., Sarr, M., Maiga, A.S. and Sissoko, G. (2012) Research Journal of Applied Sciences, Engineering and Technology, 4, 4646-4655.</mixed-citation></ref><ref id="scirp.81875-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Ly Diallo, H., Wade, M., Ly, I., NDiaye, M., Dieng, B., Lemrabott, O.H., Ma&amp;iuml;ga, A.S. and Sissoko, G. (2002) Research Journal of Applied Sciences, Engineering and Technology, 4, 1672-1676.</mixed-citation></ref><ref id="scirp.81875-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Ly Diallo, H., Dieng, B., Ly, I., Dione, M.M., Ndiaye, M., Lemrabott, O.H., Bako, Z.N., Wereme, A. and Sissoko, G. (2002) Research Journal of Applied Sciences, Engineering and Technology, 4, 2626-2631.</mixed-citation></ref><ref id="scirp.81875-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Dione, B., Sow, O., Wade, M., Ibrahima, L.Y., Mbodji, S. and Sissoko, G. (2016) Experimental Processus for Acquisition Automatic Features of I-V Properties and Temperature of the Solar Panel by Changing the Operating Point. Scientific Research Publishing Inc., 3985-4000.</mixed-citation></ref><ref id="scirp.81875-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Na, W., Carley, T., Ketcham, L., Zimmer, B. and Chen, P. (2016) Journal of Power and Energy Engineering, 4, 61-76. https://doi.org/10.4236/jpee.2016.49006</mixed-citation></ref><ref id="scirp.81875-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Diouf, M.S., Ly, I., Wade, M., Diatta, I., Traore, Y., Ndiaye, M. and Sissoko, G. (2016) Journal of Scientific and Engineering Research, 3, 289-297.</mixed-citation></ref><ref id="scirp.81875-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Sy, K.M., Diène, A., Tamba, S., Diouf, M.S., Diatta, I., Dièye, M., Traoré, Y. and Sissoko, G. (2016) Journal of Scientific and Engineering Research, 3, 433-445.</mixed-citation></ref><ref id="scirp.81875-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Kechar, E., Azzag, E. and Toua&amp;iuml;bia, I. (2015) International Journal of Scientific Research &amp; Engineering Technology, 3, 71-77.</mixed-citation></ref><ref id="scirp.81875-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Abbes, H., Abid, H., Loukil, K., Toumi, A. and Abid, M. (2013) International Journal of Control, Energy and Electrical Engineering, 17, 1-6.</mixed-citation></ref></ref-list></back></article>