Implications of Gravitational Particle CreationView this Special Issue
Research Article | Open Access
Reissner-Nordström Black Holes Statistical Ensembles and First-Order Thermodynamic Phase Transition
We apply Debbasch proposal to obtain mean metric of coarse graining (statistical ensemble) of quantum perturbed Reissner-Nordstöm black hole Then we seek its thermodynamic phase transition behavior. Our calculations predict first-order phase transition which can take Bose-Einstein’s condensation behavior.
Every observation in any arbitrary system is necessarily finite which deals with a finite number of measured quantities with a finite precision. A given system is therefore generally susceptible of different, equally valid descriptions and building the bridges between those different descriptions is the task of statistical physics (see introduction in  for more discussion). Nonlinearity property of Einstein’s metric equations causes their averaging to be nontrivial. Various possible ways of averaging the geometry of space time have already been proposed by [2–8], but none of them seems fully satisfactory (see section 7 in  for full discussion). Debbasch used an alternative way to averaging Einstein’s metric equation in . To do so, he chose a general framework where the mean metric still obeys the equations of general theory of relativity. In his approach averaging and/or coarse graining a gravitational field changes the matter content of space time called “apparent matter” which in cosmological context is related to the dark energy (see [9–12]). So general relativity mean field theory can propose a physical meaning for unknown cosmological dark energy/matter via the “apparent matter”. In the Debbasch approach, statistical ensemble of metric is ensembles of histories and not ensembles of states. This is different basically with ordinary statistical mechanics of classical and/or quantum particles. From the latter point of view, it has been known for a long time that black holes in asymptotically flat space times do not admit stable equilibrium states in the canonical ensemble (see introduction in ). But from the former point of view the Debbasch gives in  general proposal to obtain a mean field theory for the general theory of relativity. In his model members of the ensembles will be labeled by the symbol where is an arbitrary probability space . To each , there are corresponding metric tensor , compatible connections , and the Einstein metric equation (see  and section 2 in ). All members of the ensemble correspond to the same macroscopic history of the space time manifold, in particular to a given same mean metric and corresponding mean connection . As application of his model Debbasch and coworkers considered statistical ensemble of Schwarzschild black holes as nonvacuum solutions of mean Einstein metric equation by using Kerr-Schild coordinates They calculated nonvanishing temperature of mean metric where single Schwarzschild black hole is well known which has nonvanishing temperature as a vacuum solution of the Einstein equation. They discussed their results with special emphasis on their connections with the context of astrophysical observations . Extreme RNBH with has vanishing temperature (see next section) and regular Kerr-Schild coordinates are not applicable to obtain mean metric similar to the Schwarzschild one because the coarse graining space time turns out not to be a black hole . Hence Chevalier and Debbasch used analytic continuation of the Kerr-Schild coordinates as to obtain mean metric of extreme classical black hole in . According to the Debbasch approach we are free to choose types of coarse graining and/or ensemble space to obtain mean metric of the space times ensemble under consideration. We should point that topology of ensemble space times must be similar to topology of their mean metric (see ) which restrict us to choose an analytic continuation of Kerr-Schild coordinates for extreme RNBH. In short, with Debbasch proposal the averaging process does not change topology between ensemble of the curved space times and the corresponding mean space time. Precisely, the averaging process modifies the horizon radius and changes the energy-momentum tensor of space time but not total energy or mass of the black holes ensemble. Really the averaging process just redistributes without any change in the total mass which means that the total energy of the black holes dose not change by the coarse graining proposal.
Similar to study of thermodynamic behavior of single RNBH , we seek thermodynamic aspect of mean metric of nonextreme RNBHs ensemble in this work, by applying the Debbasch approach to evaluate the mean and/or coarse graining metric. Organization of the paper is as follows.
In Section 2, we calculate mean metric of ensemble of RNBHs. In Section 3 we obtain locations of mean metric horizons. In Section 4 we calculate interior and exterior horizons entropy, temperature, heat capacity, Gibbs free energy, and pressure of RNBHs mean metric. In Section 5 we calculate interior and exterior horizons luminosity and corresponding mass loss equation of quantum perturbed RN mean metric. Section 6 denotes concluding remark and discussion.
2. RNBHs Ensemble and Mean Metric
Exterior metric tensor of a single charged, nonrotating, spherically symmetric body is given byThis is metric solution of Einstein-Maxwell equation and is called RNBH in which and are corresponding ADM mass and electric charge defined in units where Equating for arbitrary spherically symmetric hypersurface , one can obtain apparent (exterior) horizon radius as and Cauchy (interior) horizon radius as which appear only for One can obtain mass independent relation between and as With particular choice (called extreme and/or Lukewarm RNBH) these horizons coincide as Clearly the RNBH metric solution (1) leads to Schwarzschild one by setting for which we will have and . Temperature of a single RNBH can be obtained for interior and exterior horizons as  which reduce to a zero value for extreme (Lukewarm) RNBH because of They show positive (negative) temperature for exterior (interior) horizons. Negative temperatures of systems have physical meaning and happen under particular conditions. More authors studied conditions where the physical systems are taken to have negative temperatures. See  for temperatures of interior and exterior horizons of Kerr-Newman black hole. One can see [17–19] for negative temperature of nongravitational systems. In the nature, materials are obtained which have interesting properties like negative refraction index and reversibility of the Doppler’s effect, and so the phase and group velocity (velocity of energy propagation) have opposite singes. In these systems temperature will have negative values (see  and references therein). Such systems are called dual system (left-handed) of direct counterpart (right-handed conventional materials). Absolute temperature is usually bounded to be positive but its violation is shown in  by Braun et al. They showed, under special conditions, however negative temperatures where high energy states are more occupied than low energy states. Such states have been demonstrated in localized systems with finite, discrete spectra. They used the Bose-Hubbard Hamiltonian and obtained attractively interacting ensemble of ultra-cold bosons at negative temperature which are stable against collapse for arbitrary atom number. Furman et al. studied in  behavior of quantum discord of dipole-dipole interacting spins in an external magnetic field in the whole temperature range They obtained that negative temperatures, which are introduced to describe inversions in the population in a finite level system, provide more favorable conditions for emergence of quantum correlations including entanglement. At negative temperature the correlations become more intense and discord exists between remove spins being in separated states. According to the documentation and looking to diagrams of the present work, one can be convinced that a quantum perturbed mean metric of coarse graining RNBHs will be exhibited finally with a first-order phase transition and Bose-Einstein condensation state microscopically. According to the Debbasch approach  ensemble of the nonextreme RNBHs is collection of coarse graining RNBHs indexed by a 3-dimensional real parameter where is the three balls of radius as follows:The metric solution (1) is convenient to be rewritten with Kerr-Schild coordinates by transformingas follows (see [10–12]):whereand is the Euclidean norm of the vector . It should be pointed that all metric solutions of Einstein’s field equation will have simple form by using Kerr-Schild coordinates. They are decomposed into the well-known flat Minkowski background metric and null vector fields as where and is a scalar function (see  and references therein). Now, we must choose a probability measure. Hence we follow the assumption presented in  and choose uniform probability measure in which is probability density of this measure with respect to Lebesgue measure as with Applying the Kerr-Schild radial coordinate (in case of extreme RNBH where we must use analytic continuation of the Kerr-Schild coordinates as (see discussion given in the introduction)), we extend single RNBH metric (4) to obtain metric of coarse graining and/or statistical ensemble of RNBHs as follows:where and Using perturbation series expansion method and averaging the metric (6) against we obtain mean metric of (6) such that (see  for details of calculations)where , ,andIt is simple to show that the mean metric (7) reduces to a single RNBH metric (4) by setting . We can rewrite the mean metric (7) in the static frame by defining the Schwarzschild coordinates. To do so, we first choose a suitable local frame with coordinates asandwhereIn the latter case the mean metric (7) readswhere we definedandWe now seek location of mean metric horizons.
3. Horizons Location for Mean Metric
One can obtain event horizon location of the mean metric (15) by solving and location of apparent (interior and exterior) horizons by solving null condition which leads to the equation such thatThe above equation has not exactly analytic solution for but for small we can use perturbation series expansion to evaluate the event horizon location. To do so we first define for which the horizon equation (18) can be written as The latter equation has a real solution as for We know that for a single RN black hole and so the condition reads for which horizon of the ensemble of statistical RN black holes is not destructed by raising if we want to apply perturbation series expansion method to obtain asymptotically behavior of the event horizon solution versus the parameters Thus we must obtain perturbation series expansion form of the event horizon but for as follows. Insertingand solving (18) as order by order, we obtainwhere and denote apparent exterior and Cauchy (interior) horizon radiuses of the mean metric (7), respectively. Inserting (9) and (19), one can obtain perturbation series expansion of (12) which up to terms in order of becomeswhere we definedArea equation of apparent horizon hypersurface of the spherically symmetric static mean metric (15) is defined by which up to terms in order of readswhere we definedAccording to Bekenstein-Hawking entropy theorem we have the result that given by (23) will be entropy function of exterior (interior) horizon of the mean metric (15). Black holes containing multiple horizons have several corresponding temperatures. Such a black hole will be in-equilibrium thermally throughout the space time where the temperature has a gradient between the horizons. Thermal equilibrium is possible only if horizon radiuses and so the corresponding temperatures become equal (see, for instance, [22, 23]). The latter situations happen for an extreme RNBH where and so We now calculate thermodynamic characteristics of interior and exterior horizons of the nonextreme mean metric of RNBHs statistical ensemble.
4. Mean Metric Thermodynamics
In the next section we will consider massless, chargeless quantum scalar field effects on luminosity of the quantum perturbed coarse graining RNBHs where its electric charge becomes invariant quantity. Hence it is useful to define dimensionless black hole mass and ensemble factor in what follows. In the latter case exterior horizon entropy of mean metric (15) can be obtained up to terms in order of as follows:and its interior horizon entropy becomeswhere andDiagrams of entropies (25) and (26) are plotted versus in Figure 4. They show that for a single RNBH () in limits but for an ensemble of RNBHs for which we use , they reach infinity In fact for physical systems the entropy itself must be positive function but its variations may reach some negative values. Hence we define difference between interior horizon entropy and exterior horizon entropy asand total entropy such as follows:Diagrams of and are plotted in Figure 3. Fortunately these diagrams show that, for a single RNBH where , we will have by decreasing and but for ensemble of RNBHs with we have while . Hence and should be considered as physical entropies of coarse graining RNBHs. Decrease of entropy causes some negative temperatures (see Figure 2) in thermodynamic systems containing bounded energy levels. In the latter case there is a critical temperature for which the system exhibits a phase transition reaching Bose-Einstein condensation state microscopically. In thermodynamics, increase of entropy means an increase of disorder or randomness in natural systems. It measures heat transfer of the system for which heat flows naturally from a warmer to a cooler substance. Decrease of entropy means an increase of orderliness or organization of microstates of a system. To do so the substance of a system must loose heat in the transfer process. Individual systems can experience negative entropy, but overall, natural processes in the universe trend toward positive entropy. Negative entropy was first introduced for living things by Ervin Schrödinger in 1944 as the reverse concept of entropy, to describe the order that can emerge from chaos . The heat generated by computations in the information theory is other applications for negative entropy concept (see [25–28] for more discussions). However we consider and to be physical entropies of RNBHs statistical ensemble containing two horizons which is in accord with positivity condition of the Bekenstein-Hawking entropy theorem. Our coarse graining RNBHs can be considered as a two-level thermodynamical system with upper bound finite energy because it has two dual (interior and exterior) horizons. We now calculate exterior (interior) horizon temperature of the RNBHs mean metric (15) as follows:Their diagrams are plotted against in Figure 2 for . For we see that has some negative (positive) values and their sign is changed when . We also plotted diagram for versus in Figure 2. They show that for reaching zero value at for While , for after that to obtain a finite positive maximum value. This maximum has smaller value for with respect to situations where we choose In ordinary statistical physics, negative temperatures are taken into account when the system has upper bound (maximum finite) energy for which entropy is continuously increasing but the energy and temperature decrease and vice versa. In the latter case the system reaches Bose-Einstein condensation state microscopically. Energy upper bound of our system is its total mass for which we have Regarding quantum matter effects on mean metric we will show in Section 5 that mass of mean metric decreases finally as (see Figure 1). Bose-Einstein condensation state needs a phase transition which happens when sign of heat capacity is changed. Hence we now calculate interior and exterior horizon of mean metric heat capacity which up to terms in order of , at constant electric charge and ensemble radius , becomeTheir diagrams are plotted against in Figure 5. They show that sign of is changed at for but sign of is changed at for We plot also diagrams of versus in Figure 5. They show a changing of sign for when and but not for In case we see for but its absolute value exhibits a minimum value. When we see which decreases monotonically to negative infinite value for Changing of sign of exterior horizon heat capacity means that a phase transition happens when the quantum perturbed RNBHs ensemble reaches its stable state with minimum mass To determine order kind of this phase transition we should study behavior of the corresponding Gibbs free energy as follows.
Exterior and interior horizon Gibbs free energies are defined bywhere entropy is given by (24) and electric potential is defined byInserting , and (30), the above Gibbs energy equation readsin which we haveWe plot diagrams of the above equations against in Figure 6. They show that has minimum zero value at but raises to by decreasing for In case , we see when Furthermore we plot diagrams of versus in Figure 6. We see when for but for . decreases to a positive minimum value by decreasing and then reaches positive infinite value. The latter behavior shows changing the sign of first derivative of when decreases and/or which means that the phase transition is first order.
One of other suitable quantities which should be calculated is pressure of black hole microparticles which coincide with the interior horizon as follows. If a quantum particle is collapsed inside of the interior (exterior) horizon then its de Broglie wave length must be at least We use de Broglie quantization condition on quantum particles as where is Planck constant and is momentum of in-falling quantum particles inside of the horizons. In Plank units where we can writein which is difference of momentum of quantum particles which move from exterior horizon to interior horizon. For they move for durations We now use the latter assumptions to rewrite Newton’s second law as is dimensionless force which affects interior horizon surface. When the system becomes stable mechanically, then must be balanced by the electric force of the system defined by Spherically symmetric condition of the system causes choosing some radial motions for quantum particles located inside of the statistical ensemble of RNBHs. However one can define pressure of moving charged quantum particle on the interior horizon aswhich by inserting (22) and using some simple calculations readsWe plot diagram of the above pressure in Figure 8. They show that in case for all values of Diagrams show that is vanishing when Also we plot diagram for versus It shows for In case where one can see when for but The latter results predict dark matter behavior of the interior horizon matter counterpart where for positive mass there is some “negative” pressure. How can mass of mean metric RNBHs decrease? Dynamically this is possible if we consider corrections of quantum matter field interacting with the mean metric of RNBHs as follows. This makes the mean metric of RNBHs unstable quantum mechanically. In the next section we assume interaction of the mean metric of RNBHs statistical ensemble with massless, chargeless quantum scalar field for which will be invariant of the system and so there is not any electromagnetic radiation. In other words there will be only mass interaction between quantum scalar field and ensemble of the RNBHs. They reduce usually to the well-known Hawking thermal radiation of the quantum perturbed mean metric which is causing mass loss of the mean RNBHs. For such a quantum mechanically unstable mean metric we now calculate its luminosity, mass loss process, and switching off effect.
5. Mean Quantum RNBH Mass Loss
We applied massless, chargeless quantum scalar field Hawking thermal radiation effects on single quantum unstable RNBH and calculated time dependence mass loss function in . We obtained that the evaporating quantum perturbed RNBH exhibits switching off effect before its mass disappears completely. It should be pointed that electric charge of the black hole is invariant of the system because there is no electromagnetic interaction between its electric charge and chargeless quantum matter scalar field. Thus mass of the RNBH decreases to reach nonvanishing remnant stable mini Lukewarm black hole with . In other words its luminosity is eliminated while its mass is not eliminated completely (see figures 9, 10, and 11 given in ). Here we study mass loss and switching off effect of quantum perturbed mean metric (15). This is a dynamical approach to describe that how mean metric of RNBHs statistical ensemble exhibits a phase transition leading to a possible Bose-Einstein condensation state microscopically. Line element of the evaporating mean metric (15) can be written near the exterior horizon as Vaida form (see, for instance, ):with the associated stress energy tensorwhere is advance Eddington-Finkelstein coordinates system. Subscript denotes the word , and denotes expectation value of quantum matter scalar field stress tensor operator evaluated in its vacuum state. The black hole luminosity is defined by the following equation from point of view of distant observer located in Applying (40) and (41), (42) becomeswhere negative sign describes inward flux of negative energy across the horizon. This causes the mean metric horizon of RNBHs statistical ensemble to shrink. In the latter case quantum particles of matter content of the black hole are in high energy state and so one can assume that the quantum black hole behaves as a black body radiation of which luminosity is defined by well-known Stefan-Boltzman law as follows:where is surface area of the black body and is its temperature. is Stefan-Boltzman coupling constant of which dimensions become as in units . If (44) satisfies (43), then we can obtain mass loss equation of the mean metric of RNBHs statistical ensemble such thatwhere the normalization constant depends linearity on the number of massless, chargeless quantum matter fields and will control the rate of evaporation. Inserting (27) one can show that the luminosity (44) for RNBHs mean metric (15) becomeswhere and should be inserted from (25) and (30), respectively. Applying (19), (20), (21), (22), (27), and some simple calculations, we can show that the mean mass loss equation (45) for RNBHs mean metric (15) readswhere we used (19), (20), (21), (22), (27), , , and to calculate which up to terms in order of become given in (47) is integral constant for which evaporating mean mass of RNBHs statistical ensemble reaches its final value as Also we defined dimensionless advance Eddington-Finkelstein time coordinate as follows:When exterior horizon of quantum evaporating RNBHs mean metric reduces to scale of its interior horizon as then one can use similar equations for luminosity and mass loss equations (46) and (47) for interior horizon as follows:whereDiagrams of the luminosities (46) and (50) and the evaporating mean RNBHs mass loss equations (47) and (51) are plotted versus mass parameter in Figure 7. They show that evaporating quantum unstable mean mass of RNBHs final state reaches remnant stable cold mini Lukewarm RNBH with final mass where its causal singularity is still covered by its shrunken horizon and its luminosity reaches zero value. We see that invariant conditions on the black hole electric charge causes the Penrose cosmic censorship hypothesis to be valid while the black hole metric is evaporated where the casual singularity of mean metric (15) defined by is still covered by their smallest scale horizons hypersurface with no naked singularity.
6. Summary and Discussion
According to the Debbasch approach we calculated mean metric of RNBHs statistical ensemble to obtain locations of interior and exterior horizons. We calculated corresponding entropy, temperature, heat capacity, Gibbs free energy, and pressure. At last section of the paper we considered interaction of massless, chargeless quantum scalar matter field on quantum perturbed mean metric of coarse graining RNBHs. Our mathematical calculations predict evaporation of the mean metric which reduces to a remnant stable mini black hole metric with nonvanishing mass. Before the evaporation reaches its final state, the mean metric exhibits a first-order phase transition and Bose-Einstein condensation state happens microscopically. Our results approve outputs of the published work  qualitatively in which the author studied thermodynamic behavior of a single RN black hole.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
- F. Debbasch, “What is a mean gravitational field?” The European Physical Journal B - Condensed Matter, vol. 37, no. 2, pp. 257–269, 2003.
- T. Buchert, “On average properties of inhomogeneous fluids in general relativity: dust cosmologies,” General Relativity and Gravitation, vol. 32, no. 1, pp. 105–125, 2000.
- T. Buchert, “On average properties of inhomogeneous fluids in general relativity: perfect fluid cosmologies,” General Relativity and Gravitation, vol. 33, no. 8, pp. 1381–1405, 2001.
- T. Futamase, “A new description for a realistic inhomogeneous universe in general relativity,” Progress of Theoretical and Experimental Physics, vol. 86, no. 2, pp. 389–399, 1991.
- T. Futamase, “General relativistic description of a realistic inhomogeneous universe,” Progress of Theoretical and Experimental Physics, vol. 89, no. 3, pp. 581–597, 1993.
- T. Futamase, “Averaging of a locally inhomogeneous realistic universe,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 53, no. 2, pp. 681–689, 1996.
- M. Kasai, “Construction of inhomogeneous universes which are Friedmann-Lematre-Robertson-Walker on average,” Physical Review Letters, vol. 69, no. 16, pp. 2330–2332, 1992.
- R. M. Zalaletdinov, “Averaging problem in general relativity, macroscopic gravity and using Einstein’s equations in cosmology,” Bulletin of the Astronomical Society of India, vol. 25, p. 401, 1997.
- C. Chevalier and F. Debbasch, “Is matter an emergent property of Space time?” 2010, https://www.semanticscholar.org/paper/Is-Matter-an-emergent-property-of-Space-Time-Chevalier-Debbasch/170660f108a1eaaa463144dcdeeec9d10f1c6a26.
- F. Debbasch and Y. Ollivier, “Observing a Schwarzschild black hole with finite precision,” Astronomy & Astrophysics , vol. 433, no. 2, pp. 397–404, 2005.
- C. Chevalier and F. Debbasch, “Thermal statistical ensembles of classical extreme black holes,” Physica A: Statistical Mechanics and Its Applications, vol. 388, no. 5, pp. 628–638, 2009.
- C. Chevalier, M. Bustamante, and F. Debbasch, “Thermal statistical ensembles of black holes,” Physica A: Statistical Mechanics and Its Applications, vol. 376, pp. 293–307, 2007.
- S. W. Hawking and D. N. Page, “Thermodynamics of black holes in anti-de Sitter space,” Communications in Mathematical Physics, vol. 87, no. 4, pp. 577–588, 1983.
- G. Grimmett and D. Stiraker, Probability and Random Process, Oxford University Press, Oxford, UK, 2nd edition, 1994.
- H. Ghaffarnejad, “Classical and quantum Reissner-Nordström black hole thermodynamics and first order phase transition,” Astrophysics and Space Science, vol. 361, article 7, 2016.
- L. Bo and L. Wen-Biao, “Negative temperature of inner horizon and planck absolute entropy of a Kerr-Newman black hole,” Communications in Theoretical Physics, vol. 53, no. 1, pp. 83–86, 2010.
- J. C. Flores and L. P. Chilla, “Theoretical thermodynamics connections between dual (left-handed) and direct (right-handed) system: entropy, temperatue, pressure and heat capacity,” Physca B, vol. 476, p. 88, 2015.
- S. Braun, J. P. Ronzheimer, M. Schreiber et al., “Negative absolute temperature for motional degrees of freedom,” Science, vol. 339, no. 6115, pp. 52–55, 2013.
- G. B. Furman, S. D. Goren, V. M. Meerovich, and V. L. Sokolovsky, “Quantum correlations at negative absolute temperatures,” Quantum Information Processing, vol. 13, no. 12, pp. 2759–2768, 2014.
- H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Herlt, Exact Solutions of Einstein's Equations, Cambridge University Press, Cambridge, UK, 2nd edition, 2009.
- H. Ghaffarnejad and M. Farsam, “Statistical ensembles of Reissner-Nordstrom black holes and thermodynamical first order phase transition,” https://arxiv.org/abs/1603.08408.
- H. Nariai, “On some static solutions of Einstein's gravitational field equations in a spherically symmetric case,” Science Reports of Tohoku University, vol. 34, pp. 160–167, 1950.
- R. Bousso and S. W. Hawking, “(Anti-)evaporation of Schwarzschild-de SITter black holes,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 57, no. 4, pp. 2436–2442, 1998.
- E. Schrödinger, What Is Life - The Physical Aspect of the Living Cell, Cambridge University Press, 1944.
- L. del Rio, J. Aberg, R. Renner, O. Dahlsten, and V. Vedral, “The thermodynamic meaning of negative entropy,” Nautre, vol. 474, p. 61, 2011.
- S. P. Mahulikar and H. Herwig, “Exact thermodynamic principles for dynamic order existence and evolution in chaos,” Chaos, Solitons & Fractals, vol. 41, no. 4, pp. 1939–1948, 2009.
- E. Oliva, “Entropy and negentropy: applications in game theory,” in Dynamics, Games and Science, J.-P. Bourguignon, R. Jeltsch, A. A. Pinto, and M. Viana, Eds., vol. 1 of CIM Series in Mathematical Sciences, pp. 513–525, Springer, Cham, Switzerland, 2015.
- L. F. del Castillo and P. V. Cruz, “Thermodynamical formulation of living systems and their evolution,” Journal of Modern Physics, vol. 2, p. 379, 2012.
- P. Hajicek and W. Israel, “What, no black hole evaporation?” Physics Letters A, vol. 80, no. 1, pp. 9-10, 1980.
Copyright © 2019 Hossein Ghaffarnejad and Mohammad Farsam. 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.