A Chaotic Improved Artiﬁcial Bee Colony for Parameter Estimation of Photovoltaic Cells

: The search for new energy resources is a crucial task nowadays. Research on the use of solar energy is growing every year. The aim is the design of devices that can produce a considerable amount of energy using the Sun’s radiation. The modeling of solar cells (SCs) is based on the estimation of the intrinsic parameters of electrical circuits that simulate their behavior based on the current vs. voltage characteristics. The problem of SC design is deﬁned by highly nonlinear and multimodal objective functions. Most of the algorithms proposed to ﬁnd the best solutions become trapped into local solutions. This paper introduces the Chaotic Improved Artiﬁcial Bee Colony (CIABC) algorithm for the estimation of SC parameters. It combines the use of chaotic maps instead random variables with the search capabilities of the Artiﬁcial Bee Colony approach. CIABC has also been modiﬁed to avoid the generation of new random solutions, preserving the information of previous iterations. In comparison with similar optimization methods, CIABC is able to ﬁnd the global solution of complex and multimodal objective functions. Experimental results and comparisons prove that the proposed technique can design SCs, even with the presence of noise.


Introduction
The world is constantly changing, and societies demand energy to continue growing and living.New resources are necessary in order to avoid an energetic crisis.This fact is particularly true since fossil fuels have been overexploited for decades.In this sense, it is a need to find and explore new energy sources that can maintain a balance between price and cleanness.Solar energy (also called photovoltaic (PV) energy (PVE)) has attracted the attention of the scientific community because it is present all over the world.PV modules are applied to transform solar radiation into electrical energy, and in the last decades, the use of such modules has increased based on the features of PVE.The main advantages of using PVE instead other sources are that it is emission-free, available all over the world and easy to install [1].However, since PVE use is focused on domestic purposes, the installation cost is expensive and also the cost of maintenance is higher.These situations occur primarily because the technology is not completely developed.Another cause is the outdoors environment that directly affects the solar modules making necessary their frequent replacement [2].In the same context, the efficiency of photovoltaic modules (or PV cells) depends on environmental factors such as temperature or radiation that cannot be controlled [3][4][5].This situation can be solved by including energy storage systems that increase the cost of the PV system.
Based on the drawbacks described above it is necessary to generate new alternative methods that help to increase the use of PVE.For example, to improve the performance of SCs requires an accurate method for design and modeling them.This situation has attracted the attention of researchers that are looking for new approaches that contribute to the design of efficient PV modules.This task is not trivial, but it is helpful for modeling, testing, control and simulation of PV systems [6][7][8][9].
The design process of SCs requires the definition of a relationship between the current (I) and voltage (V) considering the internal parameters of the cells.A mathematical model is then used to generate a representation of the elements that conforms this device, in practical terms, here is simulated the desired output of current vs. voltage (I-V).In the related literature, two approaches are used for to generate this output, (1) the single diode (SD) model and (2) the double diode (DD) model.Such methods involve electrical circuits that define the PV modules.The SD has five parameters, and the DD has seven parameters [1], their values are unknown, and their proper calibration defines the performance of the SC.In this context, it is necessary to estimate the diode saturation current, the series resistance, and the diode ideality factor for both circuits.Considering this fact, the biggest problem could be summarized as developing a mechanism able to find the best configuration of the parameters that approximate the results to the experimental data from the real SC [1].
The identification of the best parameter of PV cells can be formulated as an optimization problem.Considering the Root Mean Squared Error (RMSE) as an objective function is possible to define an algorithm to search the optimal values for the SD or DD model.The aim of using RMSE is to reduce the difference between the output of the mathematical model and an experimental dataset.It is important to notice that the experimental data is commonly obtained from measurements that involve imprecision that is reflected as noise.As a result, the optimization problem is established in a multimodal search space that contains several suboptimal solutions increasing its complexity [10,11].
To affront the problems of solar cells design, the use of classical deterministic methodologies has been proposed.Some of the most important approaches are:
Moreover, other interesting approaches like the one proposed in [15] consider the diode models as dynamic systems to estimate the best values.The use of this kind of methods implies some drawbacks, for example, the application of differentiability and convexity [16].Moreover, the initialization of the candidate solutions affects their convergence, and they can be trapped in suboptimal solutions [16].
On the other hand, stochastic techniques are a good alternative to overcome the disadvantages that the deterministic methods present.Heuristic and metaheuristic algorithms are part of stochastic approaches.They are robust techniques able to explore complex search spaces and accurately find the best solutions considering simple initial conditions [1,17,18].Some examples of these algorithms applied to modelling SCs are the following:
These algorithms are a good alternative to finding the best parameter configuration of SC diode models.However, based on the No-Free-Lunch (NFL) theorem a single metaheuristic algorithm cannot be used to solve all the optimization problems [24].In other words, depending on the complexity of the search space, the search techniques could fall into local optima.This fact affects the performance and convergence in the iterative process.However, even with these limits, metaheuristic approaches have higher probabilities of obtaining a global solution in the SC design in comparison to deterministic methods.
The Artificial Bee Colony (ABC) is a metaheuristic algorithm that mimics the behavior of honey bees [25].Different studies about the performance of the ABC algorithm show that it is able to obtain optimal solutions to various real-life problems [26][27][28][29][30].The main advantage of ABC is the absence of a local search strategy in the iterative process.Under the ABC perspective, there exist food source positions, nectar amount and different kinds of honeybees.The food sources are positions in the search space; each food source has a nectar amount value that determines its quality.The nectar amount is defined by the fitness function.ABC employs three different operators that are capable of avoiding local optimal values in complex problems.Each operator represents one of the various kinds of bees that are used on the ABC.They are the worker bees, the onlooker bees, and the scout bees.The bees represent the entire optimization process that involves the exploration and exploitation of the search space.However, like other similar approaches, ABC has some parameters that must be carefully defined to obtain a good performance [31].The main three parameters of ABC are the number of food sources, the limit value to determine the abandoned solutions and the number of foraging cycles.Such parameters are not random values, and they are selected depending on the implementation and problem.Considering this fact, the correct setting of these parameters is a complex task that is performed by the designer.In addition, the ABC applies the exploitation process in its search space by using onlookers and worker bees.The onlooker bee phase starts after finishing the worker bee phase depending on a probability value.In this context, the exploitation process is performed depending on this value, and it changes in each foraging cycle.In our study, we found that the selection of the probability value is a challenge and should be done carefully.Therefore, this paper uses a chaotic map to control this value.
On the other hand, non-linear dynamic systems define chaos as the behavior of a complex system, where small changes in the starting conditions can lead to very significant changes over time.Such changes can be random and unpredictable.The properties of chaotic systems have been applied to several optimization techniques to improve the accuracy of the optimization algorithm or to escape from local minima.The use of chaos instead random signals in metaheuristics considerably improve their performance [32].In most of the cases the chaotic versions of these optimization methods, increase the diversity of the solution and the capability to avoid local solutions in the search space.
This paper introduces the Chaotic Improved Artificial Bee Colony (CIABC) for the problem of solar cell design.In this context this article has two main goals, (1) generate an enhanced version of the ABC (CIABC); (2) apply CIABC for the parameter estimation of solar cells and photovoltaic modules.The standard ABC has been previously used for parameter estimation of SC [21].However, considering that SC modeling is a complex optimization problem and according to the NFL.The use of ABC does not represent the most accurate solution.CIABC employs chaotic maps instead of the random values on the onlooker bee step of the optimization process.Moreover, another improvement is introduced in the scout bee phase.It permits the use of the best solution so far to generate new elements of the population, instead of the use of random values.Based on such modifications the CIABC perseveres the information of the best option at each iteration and also enhance its performance converge using chaos.In this sense, the aim of this paper is to present an alternative improved method that accurately estimates the parameters that define the output of SC and photovoltaic (PV) panels.In this implementation, a selected dataset of measurements and the CIABC is applied over the diode models to minimize the Root Mean Squared Error (RMSE) [33,34].The RMSE is the fitness function, and for this problem, it determines if the values of the SC or the PV panels designed by CIABC is close to the values from the dataset.
For experimental purposes two datasets are used in this paper.The first one is for single SC, and it was extracted for the datasheet provided by the manufacturer.The second dataset corresponds to experimental information extracted from real PV modules [35].The aim of using two information sources is to provide evidence of the capabilities of the CIABC for benchmark and real data.In addition, the use of the second dataset provides evidence that the proposed algorithm can work with models that include a different amount of diodes.The property of the standard ABC to work with multi-dimensional search spaces is also contained in the CIABC.An increment on the number of diodes in the model is reflected in the number of dimensions of the optimization problem.On the other hand, to show the proficiency of the CIABC, several experiments are performed and comparisons with similar algorithms selected from the state-of-the-art.Some statistical comparisons have also been conducted in order to verify the efficacy of the improvements included in the ABC.In addition, the computational effort of the CIABC is compared with the standard version of ABC.Experimental evidence indicates that CIABC is practically immune to the sensitivity generated by noisy conditions and it has a high performance regarding accuracy, robustness and preserving low computational requirements.
The remainder of the paper is organized as follows: Section 2 describes the Preliminaries, and the problem of solar cell modeling, the standard ABC, and the chaotic map concept are introduced.In Section 3 the proposed chaotic improved ABC is introduced.Section 4 discusses the results of the proposed algorithm.The conclusions and the future works are examined in Section 5.

Photovoltaic Models
In the process of solar cell design, it is crucial to define a mathematical model that is used to estimate the internal parameters of the SC.Commonly electronic circuits help to set this model; in this sense, the single diode (SD) and double diode (DD) method are widely employed for describe SCs [36].This section introduces both SD and DD that are also adapted to be considered as an optimization problem.

Single Diode Model
The single diode model is presented in Figure 1.The equivalent circuit contains one diode that shunts the photo-generated current source I ph .This diode is configured as a rectifier and to model its non-physical ideality; the SD model considers an extra parameter [1,17,37].The SD is commonly used in the related literature for PV modeling.This circuit is very easy to implement and has only five parameters to estimate.From Figure 1, the current of the solar cell is obtained using Equation (1): In Equation ( 1), I t is the terminal current, I ph the photo-generated current, I sd is the diode current.Meanwhile I sh is the shunt resistor current.To obtain a more accurate model of the PV cells, it is used the Shockley diode equation.Therefore Equation ( 1) is rewritten as shown in Equation ( 2): The internal parameters of the diode are included in Equation ( 2), where I sd is the diode saturation current.V t is the terminal voltage whereas the series and shunt resistances are represented by R s and R sh respectively.n is a non-physical ideality factor.Some parameters are also considered for the Shockley diode equation, some of them are constants extensively used in semiconductors physics.The magnitude of charge on an electron q = 1.602 × 10 −19 C (coulombs), the Boltzmann constant k = 1.380 × 10 −23 (J/ • K) and T that is the cell temperature ( • K).In Equation ( 2), the parameters to be estimated are R s , R sh , I ph , I sd , and n.The estimation or identification of such values is reflected in the performance of SC, for that reason, this task is critical in PV systems.

Double Diode Model
Considering the double diode (DD) model equivalent circuit presented in Figure 2, the two diodes are used to shunt the photo-generated current source I ph .The first diode is configured as a rectifier.Meanwhile, the second diode represents the recombination current and other non-idealities of PV cells [1].From Figure 2, the current of the solar cell is obtained using a modified version of Equation ( 1), defined as follows: The parameters of Equation ( 3) are similar to Equation (1), the main difference is that I d1 and I d2 correspond to the first and second diode currents, respectively.In Equation ( 4) an accurate approximation of the DD obtained using the Shockley diode equation to include the internal parameter of the diode is presented: where I sd1 and I sd2 are the diffusion and saturation current, respectively.Meanwhile, the diffusion and recombination diode ideality factors are represented using n 1 and n 2 .The rest of the parameters have the same definition as in Equation ( 2).In this context, in the DD model, seven unknown parameters should be properly estimated, such elements are R s , R sh , I ph , I sd1 , I sd2 , n 1 and n 2 .
On the other hand, Table 1 shows the ranges of SD and DD parameters.The values of I sd1 and I sd2 in the DD are the same that I sd in the SD model.This situation also occurs for the intrinsic diode parameters.In this paper, the values of Table 1 are selected due that they are extensively used in the related literature [13,[38][39][40].However, a good definition of the limits of each variable is desired in order to obtain a solution with real physical meaning [41].Laudani et al. proposed an interesting analysis of the maximum and minimum values that have a physical meaning [41].However, since the aim of this paper is to address the parameter estimation as an optimization problem, the values of Table 1 possess some tolerances that permit us to create a feasible search space.

Solar Cells Design as an Optimization Problem
In the mathematical definition of both the SD and DD circuits described previously, there exist a different number of parameters that should be estimated.Using Equations ( 2) and ( 4) is possible to define an objective function that measures the quality of the estimated set of parameters (candidate solution).This equation defines if a set of parameters produce an accurate approximation between the output of the model and the measurements from the real SC.The error function for SD has then defined as: and for the DD model we have: From Equations ( 5) and ( 6) the values of V t and I t are measurements from the real CS, meanwhile, x is a vector with the parameters of the model, and it is defined as: x = R s , R sh , I ph , I sd , n for the single diode circuit and x = R s , R sh , I ph , I sd1 , I sd2 , n 1 , n 2 for the double diode circuit.The terms f SD and f DD are used to evaluate the grade of similarity of the current values computed using the estimated parameters (x) and a model with the real values defined by I t .The optimization process defined for this problem requires the minimization of the difference between the output that is modified at each iteration the values estimated by x.In this context using a dataset of N E elements, the selected objective function is the Root Mean Square Error (RMSE), and it is defined in Equation (7): From Equation ( 7), M is used to select the model DD o SD.The dataset used in the estimation could be extracted from a commercial PV cell provided by the manufactured in a datasheet or generated by experimental measurements.Here is important to mention that in most of the cases, the data is not accurate and contains a certain degree of noise.Such imprecisions are reflected in the search space where the parameters and the objective function are defined.The affectations are multi-modality and noisy features that also implies bad performance for the used search strategies [10,11,42].

Standard Artificial Bee Colony
The Artificial Bee Colony (ABC) is a swarm intelligence metaheuristic algorithm inspired by foraging behavior of real honey bees' colony in nature; it was firstly proposed by Karaboga in 2005 [25].This algorithm contains three groups of bees; the first one is the worker bees.This kind of bees searches for new sources of food around other food locations in their memory.The information about this process is passed to a second group that includes the onlooker bees.The bees of this group are used to select a food source depending on the information provided by the first group.The third consists of scout bees which randomly search for a food source.The entire process of the standard ABC is described step by step in Algorithm 1, and an explanation of its operators is also provided in this subsection.Firstly, the ABC algorithm creates a random population of N solutions that represents the first group (employed bees) x i ∈ R d , i = 1, 2, . . ., N .Each new solution v i can be created based on x i as follows [30]: where x k is a neighbor employed bee of x i , ϕ ij ∈ [−1, 1] and it is generated randomly.The fitness function value for f (x i ) and f (v i ) is calculated for x i and v i respectively; then if f (x i ) > f (v i ) the solution x i is cleared from the memory of the first group and v i is added.The fitness function value f (x i ) which is obtained from the worker bee group is passed to the second group (onlooker bees).Then the roulette wheel selection method is applied to choose the x i that has a higher probability of fitness function P i that is determined as: Each one of the onlooker bees is updating its solution by the same method used by the worker bees.The onlooker bees also check the new solution and old one to determine whether the old solution will be cleared from the memory or not.If there is no improvement in the solutions after a specific number of repetitions, these solutions are rejected; then the scout bee group searching for a new solution to update x j as: where x ij is a parameter to be optimized for the i-th worker bee, x max j and x min j are the upper and lower bounds for x ij respectively, and δ is a random number.After a new solution x ij is created, then it becomes a worker bee.where x ij is a parameter to be optimized for the i-th worker bee, x max j and x min j are the upper and lower bounds for x ij respectively, and δ is a random number.After a new solution x ij is created, then it becomes a worker bee.

Chaotic Maps
A chaotic map is defined as a method that generates non-repetitive numbers which have some properties such as stochastically intrinsic and showing irregular conduct that is also sensitive to the initial values [43].The behavior of a non-linear system is changed widely whenever the initial values present a small difference.According to these properties, the population variety can be maintained, improve the performance of determining the global optimum and escape from local optima [44,45].The chaos can be determined as a discrete-time dynamical system defined as: where t max is the dimension of the map (the number of iterations in this study).Meanwhile, f cp i k is the function that generates the chaotic model and it is described using one of the maps presented in Equation (12).Such maps are widely used in the related literature [44,45].They are many types of chaotic maps such as Sinusoidal, Logistic, Singer or Tent.In this paper, the Tent map is selected, since its performance is better than other maps according to our study and this will be illustrated in Section 4.3.The Tent map is then defined as:

Chaotic Improved Artificial Bee Colony
The proposed algorithm Chaotic Improved Artificial Bee Colony (CIABC) is introduced in this section.The improvement of the standard ABC can be explained in two steps defined as follows: (1) The use of a chaotic map.The Tent chaotic map is combined with ABC to improve the "onlooker bee" phase, which works to select a food source depending on the information provided by the first group as shown in Equation ( 9).The Tent map is used instead of the random number which is applied to start this phase.This modification makes ABC inherit the strengths of chaos such as the ergodic and non-repetition properties.This will help in improving the exploitation process as well as reduce the computational time that may occur if the random value is less than the probability for a long time.(2) Updating the solution if there is no improvement occurs with the onlooker bee.When the trial counter exceeds the limit of improving the solution in the "onlooker bee" phase, the standard ABC applies Equation (10) to generate a new solution randomly by the "scout bee" group.Whereas, the CIABC uses the best solution obtained so far to update the "scout bee" group's solution instead of Equation (10), that makes the CIABC searches around the best solution rather than random solutions.
Figure 3 illustrates the entire CIABC sequence.It starts by defining the standard ABC parameters and initializing the population.Also, it generates the matrix of a chaotic map.Then the main cycle of the algorithm starts by dealing with each of worker bees, after that, the probabilities of food sources are calculated.The onlooker bees cycle is started by checking the probability with the chaotic map (CM) value and if it greater than CM the position and fitness values will be updated.Then, if there is no improvement in the solutions after a specific number of repetitions, these solutions are rejected, and scout bee group generates a new solution using the best solution obtained.Finally, this cycle is repeated until the satisfied condition is met.

Results and Discussion
This section presents the result obtained after applying the CIABC to the problem of SC design.The experiments performed are divided in two: (1) parameter estimation of a single solar cell using and (2) parameter estimation of PV modules.In order to test the performance of the proposed method, for a single SC, a set of experimental I-V data has been used [21].Such data is applied to determine the parameter of the photovoltaic model.For experimental purposes is considered a commercial silicon solar cell (from the R.T.C. Company, Paris, France) with a diameter of 57 mm under standard test conditions (STC).For the STC the solar cell (SC) works at 1 Sun (1000 W/m 2 ) with a temperature T = 33 • C.
Under the optimization context, the search domain for the SC parameters according to the related literature [21] is presented in Table 1.On the other hand, Table 3 shows the values of the dataset; it consists of 23 samples that are widely used in several approaches for the design of photovoltaic [19,21].The CIABC method is configured using a number of food sources N = 200.Meanwhile, the maximum number of iterations is set to 10,000.The experiments were performed using "Matlab 2014b" under Windows 10, 64 bit on an Intel Core2 Duo CPU with 4 GB RAM.

Performance Measures
The performance of the proposed modified version of ABC is evaluated using different metrics; they are listed in Table 2 [21].

Metric Formula
Relative Error R err Mean Bias Error (MBE) From Table 2, I tm is the measured value, I te the estimated value, N E = 37 (the number of all experiments).The min(I te ) and max(I te ) are the minimum and maximum values (respectively) of I te over N E .

Experimental Results of CIABC
Tables 3 and 4, as well as Figures 4 and 5 illustrate the results of the CIABC method in determining the unknown parameters of SC using the single diode (SD) and the double diode (DD) models.Table 3 displays the I te and R err for SD and DD models after estimating the parameters of them, in addition to the original data.Table 4 shows the estimated values of the SD parameters which have been predicted by the CIABC model.These values are also compared with those obtained by the other selected algorithms.Comparisons illustrate that the RMSE value of the CIABC (9.8602 × 10 −4 ) is better than the original ABC and the MABC (9.862 × 10 −4 and 9.861 × 10 −4 respectively), that indicates the CIABC increased the performance of another version of ABC.Also, the RMSE of CIABC equals the values obtained by STLBO and CSO.In addition, CIABC outperformed all other algorithms; however, the worst value is 0.00139 for CPSO.Moreover, several methods are used to determine the parameters of SD model and it uses does not depend on the swarm.For example, in [48] the authors proposed a method which used the reduced forms.From reference [48], the results of the second case study (since it used the same dataset) that is listed in Table 10     Table 5 illustrates the values and the corresponding RMSE obtained for the DD model parameters by the CIABC and eight different algorithms.In this table, the RMSE value of the CIABC (9.8262 × 10 −4 ) is smaller than those obtained by standard ABC and MABC which proves that the CIABC is indeed still better than other ABC versions in both of SD and DD models.In addition, the CIABC obtains the best performance compared with all algorithms except STLBO and CSO which yielded 9.8248 × 10 −4 and 9.8252 × 10 −4 respectively.Furthermore, in Table 6 the results of several performance metrics are provided to analyze the accuracy of CIABC algorithm for both SD and DD models.These metrics are the mean and the standard deviation (STD) of RMSE values over 40 runs; as well as, the NRMSEMAE, NMAE, MBE and NMBE.These results prove that the CIABC algorithm is consistency and has a high efficiency to determine unknown parameters of SD and DD models.
Furthermore, the convergence performance of CIABC for SD and DD models are illustrated in Figure 4. From this figure can be observed that the CIABC, in both models, has fast convergence; whereby it achieved an optimal value before 2000 iterations.In the same context, the proposed CIABC requires less amount of time to achieve the optimal solutions, for a single run in SD it takes 1300 s and for DD s. Figure 5 shows the current vs. voltage and power vs. voltage at different five temperatures (33 • C, 50 • C, 75 • C and 100 • C) for the DD and SD models.It demonstrates that when the solutions yielded by CIABC are adopted, both solar cell models can accurately represent the characteristics of the solar cell.It also indicates that CIABC obtained an accurate approximation for I tm .However, the temperatures affect the current values in the CIABC and also in the power values.

Sensitivity Analysis
In this section the performance of Tent map is explained against different seven chaotic maps such as Chebyshev, Circle, Gauss/mouse, Logistic, Piecewise, Sine and Sinusoidal [49].Table 7 shows the mean and standard deviation (STD) of the RMSE values obtained using CIABC based on eight chaotic maps.From this table, is possible to conclude that Tent map has the lower mean RMSE (also the STD) compared with other maps.Also, the circle map and Piecewise map is in the second rank for dimension 5 and dimension 7, respectively.

CIABC on PV Panels
In this experimental section, another two real datasets are used to further investigate the performance of the proposed CIABC approach, where the CIABC is used to identify the parameters of PV panels that namely, a polycrystalline and monocrystalline one [35].In this experiment, the parameters of SC are modified as R sh = R p and I sh = I p in order to differentiate them from the single cell experiments.
In general, the PV panel M solar cells which are interconnected in series, as well as the current of this panel (I) equals to I M .Assuming that the same level of photon flux is received by the M cells, then the same current and voltage are generated from these cells [35].The most of the related literature the SD model is used for each element of the PV panel.According to this fact, to calculate the I of the entire PV panel, the I-V relation introduced in Equation ( 2) can be written as: Based on the values of the current (I SC ) (when the cells are short-circuited) and the voltage (V ) (when the cells are open-circuited), the I-V data is computed from each panel.The results of the proposed CIABC method are compared against two other methods called ABC and the results given in [35].

The Polycrystalline PV Panel Description
For the experiments a commercial model STM6-120/36 solar panel manufactured by Schutten Solar (Nanjing, Jiangsu, China) is considered which consists of 36 polycrystalline cells (of size 156 mm × 156 mm) in series form.For this solar panel different measures are taken that represent the real data.This panel is formed by 22 points at a temperature of 55 • C; I SC = 7.48 A, V OC = 19.21V, V M = 14.93 V and I M = 6.83 A. The entire dataset is provided in the second and third columns of Table 8.On the other hand, Tables 8 and 9 illustrates the calculated current and the identified parameters of the PV panel.From these tables, it can be observed that the results of [35] are better than the ABC method.However, the proposed CIABC gives the best results.Also, we can find that the proposed method is faster than the other algorithm which takes less time to reach the stable value.In other words, the modification performed over the CIABC does not affect the computational effort of its operators.Moreover, Figure 6 shows the estimated current and power of the three methods (ABC, CIABC, and [35]) the measured data against the voltage.From this figure, it can be seen that the curve of the proposed CIABC method is the closest curve to the curve of the measurement data.Tables 10 and 11 give the comparison results of the CIABC with the ABC and reference [35] algorithms according to the estimated current with it absolute value for error (|Error|) and the estimated parameters, respectively.From these tables, it can be observed that the CIABC has better values in terms of |Error|, MAE and RMSE.Also, the ABC is better than the results of reference [35].Also, regarding the time consumption, the CIABC need a smaller amount of time to estimate the parameters than the other algorithms.
The estimated current and the estimated power, at a temperature of 51 • C, are plotted against the voltage as in Figure 7. From this figure, it can be observed that the ABC and CIABC give better results than reference [35].

Conclusions
This paper has proposed a new method to solve the photovoltaic cell design problem by estimating the parameters of solar cells.The proposed algorithm is based on improving the Artificial Bee Colony algorithm in two stages.Firstly, the chaotic Tent map is combined with the "onlooker bee" phase instead of the random number which used to start this phase.This improvement makes ABC inherit the strength of chaos such as the ergodic and non-repetition properties and saves the computational time that may be needed if the random value is less than the probability for a long time.Secondly, instead of generating a random population using the scout bees if no improvement is made, the best solution obtained is used; that makes the ABC searche around the best value rather than random solutions.
The proposed algorithm is tested by using a dataset of experimental values.The CIABC results have been compared with eight different algorithms used in the state of the art; these algorithms are STLBO, GOTLBO, ABC, IGHS, ABSO, SA, CSO and MABC.The performance of the proposed algorithm has been evaluated based on its robustness and accuracy.The results show that the performance of the CIABC algorithm is close to the results of STLBO and CSO for SD and DD models, whereas, it is better than all other algorithms.Moreover, the effects of chaotic Tent map are compared with other different chaotic maps and its performance is better than other maps in both the SD and DD models.
In addition, the proposed CIABC was used to identify the parameters of PV panels based on polycrystalline material, and it gives results better than two other methods.In future work, we will implement other interesting models for solar cells.Also, some other modifications will be introduced, extending the use of this optimization algorithm for renewable energy problems.

Figure 1 .
Figure 1.Single diode model of a solar cell.

Figure 2 .
Figure 2. Double diode model of a solar cell.

4 .
Convergence curve of CIABC for unknown parameter identification of (a) SD and (b) DD models.

Figure 5 .
Figure 5.For the SD model: (a) Measured voltage vs. CABC-power at different temperatures; (b) Measured voltage vs. CABC computed current for different temperatures; DD model: (c) Measured voltage vs. CABC-power at different temperatures; (d) Measured voltage vs. CABC computed current for different temperatures.

Figure 6 .
Figure 6.The estimated (a) current and (b) power by the proposed CIABC method based on polycrystalline cells.

Figure 7 .
Figure 7.The Estimated (a) Current and (b) Power by the Proposed CIABC method based on polycrystalline.

Table 1 .
The range of the solar cell parameters.

While c < max iteration do
5. Generate new solutions for the worker bees by using Equation (8). 6. Evaluate worker bees, then perform the greedy selection.7.Compute the probability values by Equation (9).8. t = 0, and compute the probability 9.

While t<number of onlooker bees 10. If probability of food source > rand value 11
. Generate new solutions for the onlooker bees and evaluate them.12. Perform the greedy selection process for the onlookers.

Table 2 .
The measures used to calculate the performance.
are 1.1388 × 10 −2 , 8.8437 × 10 −4 , 8.9605 × 10 −3 and 7.7301 × 10 −4 for 2.A, 2.B, 2.C and 2.D respectively.By comparing these results with the obtained using CIABC from in Table4, it can be seen that the CIABC is better than 2.A and 2.C while it is less accurate than two other cases.

Table 3 .
Terminal (V tm I tm ) measurements and relative error values for double and single diode models.

Table 4 .
Estimated values of SD model and the corresponding RMSEusing different algorithms.

Table 5 .
Estimated values of DD model and the corresponding RMSE using different algorithms.

Table 6 .
The performance metrics for SD and DD models over 40 runs.

Table 7 .
Comparison between Tent map and seven chaotic maps for both dimensions.

Table 8 .
The calculated current and its Absolute Error (|Error|) results.

Table 9 .
The results of the parameters estimation and their accuracy for a polycrystalline cell solar panel.

Table 10 .
The current and the (|Error|) for monocrystalline PV panels.

Table 11 .
The results of the parameters estimation and their accuracy for a monocrystalline solar panel