Research Article | Open Access

Moustafa Y. Ghannam, Husain A. Kamal, "Modeling Surface Recombination at the p-Type Interface via Dangling Bond Amphoteric Centers", *Advances in Condensed Matter Physics*, vol. 2014, Article ID 857907, 9 pages, 2014. https://doi.org/10.1155/2014/857907

# Modeling Surface Recombination at the p-Type Interface via Dangling Bond Amphoteric Centers

**Academic Editor:**Mohindar S. Seehra

#### Abstract

An integral model is proposed for recombination at the silicon/silicon dioxide (Si/SiO_{2}) interface of thermally oxidized p-type silicon via *P*_{b} amphoteric centers associated with surface dangling bonds. The proposed model is a surface adaptation of a model developed for bulk recombination in amorphous silicon based on Sah-Shockley statistics which is more appropriate for amphoteric center recombination than classical Shockley-Read-Hall statistics. It is found that the surface recombination via amphoteric centers having capture cross-sections larger for charged centers than for neutral centers is distinguished from Shockley-Read-Hall recombination by exhibiting two peaks rather than one peak when plotted versus surface potential. Expressions are derived for the surface potentials at which the peaks occur. Such a finding provides a firm and plausible interpretation for the double peak surface recombination current measured in gated diodes or gated transistors. Successful fitting is possible between the results of the model and reported experimental curves showing two peaks for surface recombination velocity versus surface potential. On the other hand, if charged and neutral center capture cross-sections are equal or close to equal, surface recombination via amphoteric centers follows the same trend as Shockley-Read-Hall recombination and both models lead to comparable surface recombination velocities.

#### 1. Introduction

Recombination at the Si/SiO_{2} interface has been classically modeled using the two-charge-state Shockley-Read-Hall (SRH) statistics [1] based on Fermi-Dirac occupation probability. Defects at the Si/SiO_{2} interface, however, have been identified by Electron Paramagnetic and Electron Spin Resonance (EPR, ESR) as dangling bonds, referred to as centers, precisely located at the Si/SiO_{2} interface [2–4], and ascertained to be amphoteric [3–6]. Recombination via amphoteric centers is better described by Sah-Shockley recombination statistics [7] which has been used to model bulk recombination due to dangling bond defects in amorphous silicon (a-Si) [8] taking into account the multicharge nature of the centers and the energy correlation between them. Recently, models for recombination via amphoteric dangling bond defects at the a-Si/c-Si interface of heterojunction silicon solar cells using Sah-Shockley statistics have been developed [9, 10], and their application to recombination at the Si/SiO_{2} interface has been proposed [10].

In the present work a model for surface recombination via amphoteric dangling bond centers at the Si/SiO_{2} interface is developed based on the available model for bulk recombination in a-Si [8]. The surface recombination velocity is calculated using the proposed model and using SRH classical recombination adapted to the surface. In both models, recombination centers are assumed to form a continuum distributed throughout the energy gap and carrier dependence on surface potential at the thermally oxidized Si/SiO_{2} interface is taken into account. Comparison between the results of both amphoteric center and SRH models is presented. Finally, by fitting previously reported experimental data with the results of the proposed amphoteric center model, an attempt is made to give a convincing explanation for the double peak surface recombination velocity extracted from recombination current in gated diodes and gated transistors.

#### 2. Surface Recombination Physical Parameters and Models

Let us assume that minority electrons are injected into a p-type silicon region through a pn junction at , travel through the region, and reach its oxide passivated surface at where they recombine with holes. The band diagram at the p-type oxide passivated surface sketched in Figure 1 shows a band bending that extends over a shallow surface space charge layer in the c-Si substrate which results in an electric field extending throughout the layer (surface field layer) and in a surface potential . The surface potential is positive if the bending is downwards (surface depletion) and negative if the bending is upward (accumulation).

The surface carrier concentrations depend on the surface potential such that at low level injectionwhere is the p-type doping concentration and is the excess minority electron concentration at the inner edge of the surface space charge region associated with a split in the Fermi Energy. In (1c) bulk recombination in the very thin surface space charge (field) region is neglected which is usually a valid approximation. Hence, the surface recombination rate can be defined at the inner edge of the surface space charge region and is related to the effective surface recombination velocity through

##### 2.1. Shockley-Read-Hall (SRH) Surface Recombination Model

Surface recombination has been traditionally treated using Shockley-Read-Hall (SRH) statistics which uses Fermi-Dirac occupation probability function (occupancy), , for recombination centers at energy . At equilibrium is given by which states that the probability of occupation of a recombination center at energy is smaller than 1/2 if is larger than the Fermi Energy , greater than 1/2 if , and exactly equal 1/2 when .

The SRH surface recombination rate is expressed by an expression similar to the SRH bulk recombination rate with replacing the bulk electron and hole concentrations and by the surface potential dependent surface carrier concentrations and . Unlike previous reports the interface recombination centers in the present work do not occur at a single discrete energy level with a density (/cm^{2}) but form a continuum in the energy bandgap with a state density (/cm^{2}/eV). Consequently, the SRH surface recombination rate is to be determined from the integration over all energies in the bandgap from the valence band edge up to the conduction band edge
Using (2), the effective surface recombination velocity can then be obtained from
where is the thermal velocity (*≈*10^{7} cm/s), is the intrinsic carrier concentration, is the intrinsic Fermi Energy (midgap energy level), is Boltzmann constant, is the absolute temperature, and and represent the electron and hole capture cross-sections of the centers.

##### 2.2. Surface Recombination Amphoteric Dangling Bond Model

It has been argued that the two-state charge Fermi-Dirac occupation function has some limitations in a-Si abundantly hosting amphoteric dangling bonds [8, 11] and that multicharge Sah-Shockley correlated electron statistics [7] is more appropriate for this case. This may also apply to recombination at the Si/SiO_{2} where defects are mainly due to dangling bond amphoteric centers. Hence we propose to use the amphoteric center recombination model developed for bulk recombination in a-Si after proper adaptation to the surface as done for SRH surface recombination.

Since the amphoteric recombination model is not very well known to many readers, the main formulation aspects detailed in [8] are repeated in the appendix for convenience. Besides the flaw distribution density , the model parameters include electron capture cross-sections of positive and neutral centers and and hole capture cross-sections of negative and neutral centers and . The total surface recombination rate via a continuum of amphoteric centers in the energy gap is obtained by integrating the recombination rate at a single level over all the energies from to : and the effective surface recombination velocity follows from (2):

#### 3. Numerical Results

A MATLAB code is developed and used for calculating the effective surface recombination velocity and , according to (3c) and (4b), at the interface of an oxide passivated p-type Si surface having a doping concentration of /cm^{3}. The recombination centers, although initially having a shape distribution [5], are assumed to be uniformly distributed in the energy gap as it happens to be after annealing in hydrogen or forming gas [12–14]. After this technological step the midgap center density drops from the 10^{11}/cm^{2}/eV range to the 10^{10}/cm^{2}/eV and the correlation energy between the centers is in the range 0 to 0.5 eV. The capture cross-sections of amphoteric centers are reported to be in the range 10^{−16} cm^{2} [15, 16] and it is widely accepted that can be up to 100 times larger than [17]. Values for , , , and close to these typical values are used in the calculations carried out here.

As detailed in [8] and repeated in the appendix, the correlation energy between amphoteric centers influences the occupation probability which affects the recombination rate and velocity. Therefore two extreme values will be considered: and eV. In the first case and according to , the equilibrium occupation probability of centers considered to be positively charged is high for centers having , equals 25% for centers with , and vanishes for centers with more than below . The probability of occupation for centers considered to be negatively charged is high for centers having , equals 25% for centers with , and vanishes for centers with more than above . Finally, the probability of occupation for centers considered to be neutral is symmetric around , equals 50% for , and vanishes when is more than below or above . For the extreme case eV the center of symmetry of the occupation probability is moved from to , that is, to 0.2 eV below . In that case, the probability of occupation of positive centers is 1/3 at , is higher than 1/3 for centers having , and vanishes for more than below . The probability of occupation of negative centers is equal to for centers with * *eV, is higher for centers with lower energies, and vanishes for centers with more than above . Finally, the probability of occupation of neutral centers is 2/3 for and * *eV and almost one between * *eV and * *eV.

##### 3.1. Effective Surface Recombination Velocity versus Excess Carrier Concentration

The effective surface recombination velocity is calculated using SRH and amphoteric center recombination models and plotted as a function of the excess electron concentration at the edge of the surface field region. Since nonnegligible positive fixed charges are present in the passivating oxide layer, the p-type surface would be depleted which creates a surface field region, as shown in Figure 1, and results in band bending at the surface and a positive surface potential. Therefore calculations of are carried out at flat band condition (surface potential ) and for a positive surface potential and V. The results in Figure 2(a) represent the case where the correlation energy is equal to zero, is uniform and equal to 10^{10}/cm^{2}/eV throughout the energy gap, and all capture cross-sections are equal to 10^{−16 }cm^{2}. It is always assumed in the present analysis that the capture cross-sections of the SRH centers and of neutral amphoteric centers are equal such that and . The impact of asymmetry in other capture cross-sections (, , ) and of positive correlation energy will be discussed in later sections.

**(a)**

**(b)**

By inspecting Figure 2(a) it appears that and follow the same trend versus excess electron concentration at the surface, . The value of is larger than but stays in the same order of magnitude. Both values get closer as the surface potential is increased and maybecome larger at higher surface potential or higher correlation energy as will be seen in later sections. For several orders of magnitude in the lower range exhibits its highest value that stays constant, then starts to decay gradually in the higher range, and continues to decay to reach very small values close to unity at high values of . Such a behavior can be predicted from (3c) since, for a constant value of the surface potential, for several orders of magnitude in the low range the term present in the denominator is negligible compared to the term which results in a constant and high value of . As increases may become significantly boosted such that it starts to affect the value of the denominator of (3c) which results in a gradual decrease in the value of as depicted in Figure 2(a). At positive surface potential 0.1 and 0.2 V, which are typical values at the p-type surface due to the presence of positive fixed oxide charges in the SiO_{2} passivating layer, the value of is exponentially boosted and that of exponentially reduced as predicted by (1a) and (1b) which shifts the onset of degradation to smaller values of *.* Moreover, due to the reduced value of the value of the denominator becomes smaller which boosts the maximum constant value of at low .

###### 3.1.1. Impact of the Correlation Energy between Amphoteric Centers

A major property of recombination via amphoteric centers is that the recombination centers are energy correlated. The correlation energy affects the occupation probability and hence would certainly affect the recombination rate and recombination velocity. In the previous section the correlation energy was assumed to be zero which would logically mean a reduced recombination activity. As depicted in Figure 2(a), the surface recombination velocity via amphoteric centers indeed increases with increasing correlation energy. For a correlation energy eV close to the maximum value, the surface recombination velocity displayed in Figure 2(b) is more than doubled compared to its values plotted in Figure 2(a) for and exceeds the value of at V and 0.2 V. On the other hand, does not seem to affect the trend nor the ranges characterizing the dependence of on .

###### 3.1.2. Asymmetric Electron and Hole Capture Cross-Sections

It has been reported that the capture cross-section of electrons may be up to 100 times larger than the capture cross-sections of holes [17]. Therefore is calculated assuming that the electron capture cross-sections of SRH centers and of the neutral amphoteric centers are equal () and are 100 times larger than their hole capture cross-sections such that . Equality is also maintained between capture cross-sections of charged and neutral center ( and ). The impact of the asymmetry seems to affect only the value of since both and are boosted by approximately a factor of 100 compared to the results in Figure 2(a) which reflects the ratio between and but not the general trend it has versus , as depicted in Figure 3.

###### 3.1.3. Asymmetric Carrier Capture Cross-Sections of Neutral and Charged Centers

Another fundamental difference between SRH statistics and Sah-Shockley statistics is the possibility in the latter to have different capture cross sections for neutral and charged centers. This would lead to a discrepancy between and which was not found when capture cross sections of neutral and charged states were assumed to be equal. Indeed, a significant discrepancy is found between and when and are not equal to one, as depicted in Figure 4(a) obtained for V and eV. These results are obtained assuming equal capture cross-sections of electrons and holes of SRH and neutral centers . When > 1 () the constant value of in the low range is significantly enhanced and the range itself is strongly reduced which is the result of a boost in the product strongly affecting the electron capture rate. On the other hand, is boosted in the high voltage range when () and the range is extended to higher as a result of a boost in the product strongly affecting the hole capture rate. In the presence of both asymmetries and ( and ) exhibits a double step behavior resulting from the superposition of the low voltage and high voltage trends displayed in Figure 4(a).

**(a)**

**(b)**

##### 3.2. Effective Surface Recombination Velocity versus Surface Potential

Due to the exponential dependence of the surface carrier concentrations and on the surface potential , the surface recombination activity is also a strong function of . Such a dependence is very well established and has been investigated thoroughly in previous work [18–20] based on SRH recombination statistics. To our knowledge this work is the first rigorous analysis that treats surface recombination via amphoteric centers at the Si/SiO_{2} interface, and hence the dependence of on the surface potential needs to be investigated and compared to the classical SRH dependence. Typical results for the effective surface recombination velocity and versus surface potential are plotted in Figure 5. A single peak in the SRH recombination rate and SRH surface recombination velocity dependence on is predicted by (3b) and (3c), respectively. This peak occurs when , which, when using (1b) and (1a), would occur at a surface potential given by
For K, , and the given typical values /cm^{3} and /cm^{3} and according to (5a), would peak at V which is confirmed in Figure 5. Further increase in surface potential results in a decrease of due to the dominance of in the denominator of (3c).

On the other hand, our analysis proves that versus surface potential is distinguished from by the appearance of two peaks, one to the left and one to the right of the SRH peak, as depicted in Figure 5. The formation of these peaks is strongly dependent on the ratios and since, like SRH, they correlate to the relationship between the terms involving the products and . In this case, however, due to the presence of neutral, positively charged, and negatively charged centers, amphoteric center recombination involves four products instead of two, namely, , , , and . Assuming the competition between and at large values of leads to a first peak when occurring at given by
which, for example, for results in the low surface potential peak occurring at V which is confirmed in Figure 5. On the other hand, the competition between and at large values of leads to a first peak when occurring at given by
which, again for , results in the high surface potential peak occurring at V which is also confirmed in Figure 5. Note that, when or is equal to one, the magnitude of the corresponding peak is small and tends to be close to the SRH value, as explained in Section 3.1, and occurs exactly at the same surface potential as predicted by (5b) and (5c) when . Note that the two peak behavior is pronounced only if both and are large and comparable. On the other hand, if or are asymmetric, the peak associated with the larger is maintained while the peak associated with the smaller disappears or shows as a hump or a shoulder at the peak position determined by (5b) and (5c). It is clear that according to (5b) and (5c), the whole plot in Figure 5 may be shifted to the right or to the left by changing the excess electron concentration , and the peak positions may be stretched away from each other or brought closer to each other by modifying the ratios and . The two-peak behavior and the peak-hump behavior have been previously reported for surface recombination velocity deduced from measurements of the recombination current at the Si/SiO_{2} interface in gated diodes and gated transistors (e.g., [18, 21, 22]). The presence of mobile ions in the oxide layer [21], nonuniformities at the Si/SiO_{2} interface [18], interface states having different capture cross-sections [22] are previous interpretations for the two-peak behavior which in our opinion still need to be justified.

Using the proposed amphoteric center recombination model a successful fitting is demonstrated in Figure 6 for an experimental double peak surface recombination—versus plot previously reported in the literature for a p-type surface doping concentration /cm^{3} and excess electron concentration /cm^{3} [22]. It is possible to obtain such fitting using several combinations with reasonable values for the parameters , , , , and and therefore there is no need to insist on a specific combination. It is impossible, however, to fit these data using SRH recombination unless mathematical oriented assumptions with no physical justifications are proposed such as splitting the centers into groups with different capture cross-sections. On the other hand, the present analysis firmly attributes the two-peak surface recombination velocity to surface recombination via amphoteric centers with larger capture cross-sections of charged centers than of neutral centers.

#### 4. Conclusions

A model is proposed for surface recombination via amphoteric defects at the Si/SiO_{2} interface of thermally oxidized p-type silicon surface. The model is an adaptation to the surface of the model developed for the bulk recombination in amorphous silicon based on Sah-Shockley multicharge statistics for energy correlated amphoteric dangling bond defects. The results indicate that surface recombination via amphoteric centers behaves in a very similar way to SRH recombination when capture cross-sections of charged and neutral centers are equal or close to equal. On the other hand, if electron or hole capture cross-sections are larger for charged centers than for neutral centers the surface recombination velocity is enhanced at small or at high excess electron concentration. If both asymmetries are present, the superposition of both trends leads to a double step behavior of versus excess electron concentration. In addition, in the presence of such asymmetries the surface recombination velocity exhibits two peaks versus surface potential, which relate to a larger electron capture cross-section of positively charged centers and to a larger hole capture cross-section for negatively charged centers than for neutral centers. Expressions are derived for the surface potentials at which the two peaks occur. Recombination via amphoteric centers with asymmetric capture cross-sections for charged and neutral centers is a firm and plausible interpretation for the double peak surface recombination velocity extracted from measurements of surface recombination current in gated diode and gated transistor versus gate voltage.

#### Appendix

#### A. Expression for the Dangling Bond Amphoteric Surface Recombination Rate

The model proposed here is an extension to the model suggested for recombination via dangling bonds in the bulk of amorphous silicon [8] but applied to the recombination via dangling bond at the Si/SiO_{2} interface. Centers may be neutral , positive , or negative and correlated with a positive correlation energy as depicted in Figure 7. Transitions are allowed between dangling bonds and valence or conduction bands such that
The electron capture rates and represent electron capture processes associated with capture cross-sections and , and the rates and represent electron emission processes. The rates and represent hole capture processes associated with capture cross-sections and . The rates and represent hole emission processes.

In the dark, three independent conservation equations can be written.

##### A.1. The Steady State Electron Recombination Rate

One has where and represent, respectively, the occupation probabilities and electron capture coefficients associated with a neutral or positive center.

##### A.2. The Steady State Hole Recombination Rate

One has where and are the occupation probability and hole capture coefficients associated with a neutral or negative center, respectively.

##### A.3. The Steady State Rate of Change of

One has The emission rates are obtained from the principle of detailed balance stating that each charge state of the dangling bonds must be in equilibrium with the band states such that

###### A.3.1. Resulting in the Emission Probabilities

One has The occupation probabilities , , and at equilibrium in set (A.6) are obtained from the solution of set (A.2) through set (A.5) at equilibrium resulting in Then, by equating in set (A.2) and in set (A.3), the occupation probabilities at nonequilibrium are obtained and are given by The recombination rate via centers at is calculated by substituting set (A.8) in set (A.2) or in set (A.3) to get

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is supported by Kuwait University Research Administration Grant EE03/12. The authors are indebted to Eng. Mohamed Abdelazim Koura for helping with MATLAB.

#### References

- W. Shockley and W. T. Read, “Statistics of the recombinations of holes and electrons,”
*Physical Review*, vol. 87, no. 5, pp. 835–842, 1952. View at: Publisher Site | Google Scholar - Y. Nishi, “Study of silicon-silicon dioxide structure by electron spin resonance I,”
*Japanese Journal of Applied Physics*, vol. 10, pp. 52–62, 1971. View at: Google Scholar - P. M. Lenahan and P. V. Dressendorfer, “Effect of bias on radiation-induced paramagnetic defects at the silicon-silicon dioxide interface,”
*Applied Physics Letters*, vol. 41, no. 6, pp. 542–544, 1982. View at: Publisher Site | Google Scholar - E. H. Poindexter and P. J. Caplan, “Characterization of Si/SiO
_{2}interface defects by electron spin resonance,”*Progress in Surface Science*, vol. 14, no. 3, pp. 201–294, 1983. View at: Google Scholar - E. H. Poindexter, G. J. Gerardi, M.-E. Rueckel, P. J. Caplan, N. M. Johnson, and D. K. Biegelsen, “Electronic traps and Pb centers at the Si/SiO
_{2}interface: band-gap energy distribution,”*Journal of Applied Physics*, vol. 56, no. 10, pp. 2844–2849, 1984. View at: Publisher Site | Google Scholar - P. M. Lenahan and P. V. Dressendorfer, “An electron spin resonance study of radiation-induced electrically active paramagnetic centers at the Si/SiO
_{2}interface,”*Journal of Applied Physics*, vol. 54, no. 3, pp. 1457–1460, 1983. View at: Publisher Site | Google Scholar - C.-T. Sah and W. Shockley, “Electron-hole recombination statistics in semiconductors through flaws with many charge conditions,”
*Physical Review*, vol. 109, no. 4, pp. 1103–1115, 1958. View at: Publisher Site | Google Scholar - F. Vaillant and D. Jousse, “Recombination at dangling bonds and steady-state photoconductivity in a-Si:H,”
*Physical Review B*, vol. 34, no. 6, pp. 4088–4098, 1986. View at: Publisher Site | Google Scholar - C. Leendertz, R. Stangl, T. F. Schulze, M. Schmidt, and L. Korte, “A recombination model for a-Si:H/c-Si heterostructures,”
*Physica Status Solidi C*, vol. 7, no. 3-4, pp. 1005–1010, 2010. View at: Publisher Site | Google Scholar - S. Olibet, E. Vallat-Sauvain, L. Fesquet et al., “Properties of interfaces in amorphous/crystalline silicon heterojunctions,”
*Physica Status Solidi A*, vol. 207, no. 3, pp. 651–656, 2010. View at: Publisher Site | Google Scholar - S. Steingrube, R. Brendel, and P. P. Altermatt, “Limits to model amphoteric defect recombination via SRH statistics,”
*Physica Status Solidi A*, vol. 209, no. 2, pp. 390–400, 2012. View at: Publisher Site | Google Scholar - W. Füssel, M. Schmidt, H. Angermann, G. Mende, and H. Flietner, “Defects at the Si/SiO
_{2}interface: their nature and behaviour in technological processes and stress,”*Nuclear Instruments and Methods in Physics Research A*, vol. 377, no. 2-3, pp. 177–183, 1996. View at: Publisher Site | Google Scholar - R. R. Razouk and B. E. Deal, “Dependence of interface state density on silicon thermal oxidation process variables,”
*Journal of the Electrochemical Society*, vol. 126, no. 9, pp. 1573–1581, 1979. View at: Google Scholar - P. J. Caplan, E. H. Poindexter, B. E. Deal, and R. R. Razouk, “ESR centers, interface states, and oxide fixed charge in thermally oxidized silicon wafers,”
*Journal of Applied Physics*, vol. 50, no. 9, pp. 5847–5854, 1979. View at: Publisher Site | Google Scholar - P. M. Lenahan, N. A. Bohna, and J. P. Campbell, “Radiation-induced interface traps in MOS devices: capture cross section and density of states of P
_{b1}silicon dangling bond centers,”*IEEE Transactions on Nuclear Science*, vol. 49, no. 6, pp. 2708–2712, 2002. View at: Publisher Site | Google Scholar - D. Sands, K. M. Brunson, and M. H. Tayarani-Najaran, “Measured intrinsic defect density throughout the entire band gap at the <100> Si/SiO
_{2}interface,”*Semiconductor Science and Technology*, vol. 7, no. 8, pp. 1091–1096, 1992. View at: Publisher Site | Google Scholar - S. W. Glunz, D. Biro, S. Rein, and W. Warta, “Field-effect passivation of the SiO
_{2}-Si interface,”*Journal of Applied Physics*, vol. 86, no. 1, pp. 683–691, 1999. View at: Google Scholar - M. Y. Ghannam, R. P. Mertens, R. F. De Keersmaecker, and R. J. Van Overstraeten, “Electrical characterization of the Boron-doped Si-SiO
_{2}interface,”*IEEE Transactions on Electron Devices*, vol. 32, no. 7, pp. 1264–1271, 1985. View at: Google Scholar - M. Y. Ghannam and R. P. Mertens, “Surface recombination current with a nonideality factor greater than 2,”
*Electron device letters*, vol. 10, no. 6, pp. 242–244, 1989. View at: Google Scholar - M. W. Hillen and J. Holsbrink, “The base current recombining at the oxidized silicon surface,”
*Solid State Electronics*, vol. 26, no. 5, pp. 453–463, 1983. View at: Google Scholar - M. W. Hillen, G. Greeuw, and J. F. Verweij, “On the mobility of potassium ions in SiO
_{2},”*Journal of Applied Physics*, vol. 50, no. 7, pp. 4834–4837, 1979. View at: Publisher Site | Google Scholar - E. Yablonovitch, R. Swanson, W. Eades, and B. Weinberger, “Electron-hole recombination at the Si-SiO
_{2}interface,”*Applied Physics Letters*, vol. 48, no. 3, pp. 245–247, 1986. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2014 Moustafa Y. Ghannam and Husain A. Kamal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.