Red-tailed hawk algorithm for numerical optimization and real-world problems

This study suggests a new nature-inspired metaheuristic optimization algorithm called the red-tailed hawk algorithm (RTH). As a predator, the red-tailed hawk has a hunting strategy from detecting the prey until the swoop stage. There are three stages during the hunting process. In the high soaring stage, the red-tailed hawk explores the search space and determines the area with the prey location. In the low soaring stage, the red-tailed moves inside the selected area around the prey to choose the best position for the hunt. Then, the red-tailed swings and hits its target in the stooping and swooping stages. The proposed algorithm mimics the prey-hunting method of the red-tailed hawk for solving real-world optimization problems. The performance of the proposed RTH algorithm has been evaluated on three classes of problems. The first class includes three specific kinds of optimization problems: 22 standard benchmark functions, including unimodal, multimodal, and fixed-dimensional multimodal functions, IEEE Congress on Evolutionary Computation 2020 (CEC2020), and IEEE CEC2022. The proposed algorithm is compared with eight recent algorithms to confirm its contribution to solving these problems. The considered algorithms are Farmland Fertility Optimizer (FO), African Vultures Optimization Algorithm (AVOA), Mountain Gazelle Optimizer (MGO), Gorilla Troops Optimizer (GTO), COOT algorithm, Hunger Games Search (HGS), Aquila Optimizer (AO), and Harris Hawks optimization (HHO). The results are compared regarding the accuracy, robustness, and convergence speed. The second class includes seven real-world engineering problems that will be considered to investigate the RTH performance compared to other published results profoundly. Finally, the proton exchange membrane fuel cell (PEMFC) extraction parameters will be performed to evaluate the algorithm with a complex problem. The proposed algorithm will be compared with several published papers to approve its performance. The ultimate results for each class confirm the ability of the proposed RTH algorithm to provide higher performance for most cases. For the first class, the RTH mostly got the optimal solutions for most functions with faster convergence speed. The RTH provided better performance for the second and third classes when resolving the real word engineering problems or extracting the PEMFC parameters.

• Parallelism is utilized in population-based algorithms such as SA and HA.Different individuals are sent out at once to complete a task, and the results are compared.According to the comparison, additional ideas are utilized to assess the population evolution or to create new individuals.• Acceptance is utilized in three cases: 1. Admit temporary solutions that weaken objective function due to search space expansion, 2. Handling the cost function's constraints.There are two ways to cope with the constraints.The first way excludes any solutions if they match the start conditions.The second procedure is used if any solution can be allocated a numerical value.In this case, all solutions must participate, and the initial conditions may match infeasible solutions.3. Adding restrictions to allowed solutions that improve the optimal solution by, at minimum, the limiting level.When comparing values produced from previous calculations, this strategy aids in avoiding numerical issues.• Elitism: the best solution must be kept from one iteration to the next in repeating population-based algo- rithms.The elitism principle is utilized to accomplish this by retaining the individual who outlined the best solution and utilizing it as a reference for the following iteration or upgrading it if a better solution is located.The notion of elitism may also be applied to many individuals, directing an élite group of them to the next iteration.
• Selection is a probability-based method that generates new individuals from a pool of available ones.This approach may incorporate weights into the probabilistic selection process, in which random individuals are chosen to produce new ones.• Decay or Reinforcement: Decay allows for more initial freedom, followed by incremental flexibility constraints.
This approach is based on a decremental updating factor (< 1) used at each iteration.Reinforcement is employed similarly in some cases by applying an incremental factor higher (> 1).• Immunity is gained through finding characteristics of solutions that lead to appropriate settings.Immunity prefers solutions with characteristics similar to those attributes.• Self-Adaptation is a process that allows the parameters of algorithms to be updated in response to the progress of optimization.
A novel high-efficiency SA optimization algorithm is suggested in this study to compete with existing optimizers.The basic concept of the suggested optimizer is based on the hunting skills of Red-tailed Hawks, one of the most intelligent birds of prey.The benefits of evolutionary and swarm approaches have been included in the RTH algorithm's design to outperform the existing optimization algorithms.The proposed RTH contains three stages: high soaring, low soaring, and stooping and swooping stages.The high soaring stage is quite similar to evolutionary approaches' search behavior.The RTH algorithm relays and gathers all search points, beginning with the original position and progressing to the best position.The second and third phases imitate SA's behavior while moving to the best position, using the information of the preceding site for each point.
Despite the availability of metaheuristic algorithms, there is a continuous need for developing and creating new algorithms that may provide better performance for such problems than the actual ones.The NFL theory confirms this need for new optimizers.To this end, this paper suggests this novel algorithm.The main contribution of this paper is to propose a novel nature-inspired metaheuristic optimization algorithm for better solving optimization issues.The proposed algorithm benefit from the unique hunting strategy of the red-tailed hawk and applies it to tackle various optimization problems.The suggested RTH algorithm's performance has been examined on three optimization problem classes.The first class contains three types of optimization functions: • Twenty-three standard benchmark functions: unimodal, multimodal, and fixed-dimensional multimodal functions.
Inspiration and behavior during hunting.The red-tailed hawk (Buteo jamaicensis) is a bird of prey that breeds over much of North America, from Alaska's interior and northern Canada to Panama and the West Indies.The red-tailed hawk lives in various environments and elevations, such as deserts, grasslands, forests, agricultural fields, and cities.The red-tailed hawk is a predatory, carnivorous eater.Almost every little animal they come upon may be seen as possible prey 49 .Small mammals, such as rodents, are their most common prey, although they also eat birds, fish, reptiles, invertebrates, and amphibians.Prey varies greatly depending on geographical and seasonal availability; however, rodents comprise 85 percent of a hawk's diet 50 .
The red-tailed hawk soars with its wings in a mild dihedral, flapping as little as possible to save energy.Unlike other hawks, the red-tails can fly for long distances thanks to this feature.Because soaring is the most efficient flying mode for these hawks, it is utilized more frequently 51 .It moves between 32 and 64 km/h (20-40 mph)  when soaring or flapping its wings.The large wings allow the red-tailed to reach 190 km/h (120 mph) when plunging 52 .Red-tailed hawks can fly fast and powerfully while repeatedly diving at perceived threats during nest protection 53 .As illustrated in Fig. 1, the red-tailed has three types of flying: • High soaring (Fig. 1a): It flies highly with its wings in a mild dihedral, flapping as little as possible to save energy to explore the selected area.• Low soaring (Fig. 1b): after selecting the target position, the red-tailed fly with a low soaring in a spiral move- ment around the prey.This movement allows it to detect the best location and time to hit the target.• Stooping and swooping (Fig. 1c): after selecting the best location and moment in the previous step, the red-tailed swooped its prey by stooping and raising its acceleration (from 32-64 to 190 km/h) in a curved direction.
Mathematical model.The suggested red-tailed hawk (RTH) algorithm mimics the red-tailed hawk's hunting behavior.The actions taken at each hunt stage are presented and modeled.This algorithm includes three stages, high soaring, low soaring, and stooping and swooping.
1 High soaring: the red-tailed hawk will soar far into the sky, looking for the best location in terms of food availability.Figure 2 illustrates the behavior of the red-tailed hawks during the high soaring stage, and Eq. ( 1) represents the mathematical model of this stage : where X(t) represents the red-tailed hawk position at the iteration t, X best is the best-obtained position, X mean is the positions' mean, Levy represents the levy flight distribution function that can be calculated according to Eq. ( 2), and TF(t) denotes the transition factor function that can be calculated according to Eq. ( 3).
(1)  where T max represents the max number of iterations. 2 Low soaring: the hawk surrounds the prey by flying much lower to the ground in a spiral line.This stage is illustrated in Fig. 3, and its model can be expressed as follows: where x and y denote direction coordinates which can be calculated as follows where R 0 represents the initial value of the radius [0.5-3],A denotes the angel gain [5-15], rand is a random gain [0-1], and r is a control gain [1, 2].These parameters help the hawk fly around the prey with spiral movements, as explained in Fig. 4. 3 Stooping and Swooping: In this stage, the hawk suddenly stoops and attacks the prey from the best-obtained position in the low soaring stage.Figure 2 explains the behavior of the red-tailed hawks during this stage.This stage can be modeled as follows: where each step size can be calculated as follows ( 2) where α and G are the acceleration and the gravity factors, respectively, they can be defined as follows: where α represents the hawk's acceleration that increases with the increase of t to enhance the convergence speed, and G is the gravity effect that decreases to reduce the exploitation diversity when the hawk is much near the prey.This phase is explained in Fig. 5.  Figure 6 presents the evolution of the solutions of a 2-dimensional problem as a function of the iterations.During the first iterations, the hawk explored the search space (exploration phase).It took significant steps to detect the prey position.The length of these steps decreases as the hawk is near the prey position.The hawk spent a small number of iterations in the exploration due to the hybrid evolutionary-swarm updating mechanism.It avoided attempting on the local minim thanks to the high soaring stage based on the Levy function.The low soaring enhances the convergence speed in the exploitation phase, where the stooping and the swooping stages strengthen the RTH algorithm's accuracy.Figure 7 presents the hunting behavior for 2-dimensional and 3-dimensional problems.The hawk explores the whole search space from these figures to avoid attempting local

Results and discussion
RTH's performance was evaluated using several test functions.Three classes of test functions are used in this study: • Standard test functions 54 , • CEC 2020 benchmark test functions 55 , • CEC 2022 benchmark test functions 56 .

Standard test functions.
Unimodal, multimodal, and fixed-dimension multimodal test functions are used to evaluate the RTH performance 54 .Unimodal test functions (F1-F7) are used to put an algorithm's exploitation abilities to the test, whereas multimodal test functions (F8-F13) are used to experiment with the algorithm's exploration performance.The fixed dimension test functions (TF14-TF23) demonstrate the algorithm exploration capability 57 .These functions are provided in Table A1 (appendix).For all the tests, the population size for each algorithm has been set at 30 with max iterations of 1000.This study ran RTH 30 times and reported the findings, including the average and standard deviation values of the best-so-far solutions obtained in each run to provide statistically significant results.The test is carried out for six other algorithms to underline the effectiveness and superiority of RTH to that of different methods: Particle Swarm Optimization (PSO) 37 as the best-known optimization method, Salp Swarm Algorithm (SSA) 38 as one of the most used algorithms, Hunger Games Search (HGS) 20 , COOT algorithm 19 , Artificial Eco-system Optimizer (AEO) 4 , and Aquila Optimizer (AO) 21 as recently developed metaheuristics.Table 2 summarizes the parameter setting for the RTH algorithm, where the parameters are set by the try-and-error method.The other algorithms are based on their default parameters.Table 3 shows the results of RTH and the other used algorithms on these standard test functions using the best, average, and standard deviation (StD) values.The obtained results in Table 3 and Fig. 9 demonstrated that RTH has superior exploration capabilities to other methods for five of six functions in terms of finding the optimal solution with the lowest number of iterations.These functions approve the ability of the RTH to handle complicated problems.Fixed-dimensional multimodal function (F 14-23) similar to the multimodal test functions, these functions include a lot of local optima with more complexity due to the minimization of the negative cost function.Hence, they are suitable for evaluating the accuracy of the algorithm.As illustrated in Fig. 10 and reported in Table 3, the proposed RTH has reached the global optimum for all the considered functions.However, its robustness, represented by the standard deviation (StD), has been reduced to the increased complexity.RTH's exploration is due to its different optimization phases and the acceleration and gravity effects.

IEEE CEC 2020.
The IEEE Congress on Evolutionary Computation (CEC) offers yearly benchmark functions to test and evaluate optimization algorithms' ability to identify optimal solutions (minimization problems).This collection of functions is commonly recognized by the year as standard functions established by the IEEE Congress on Evolutionary Computation, for example, CEC 2020 and CEC 2022.These test functions include unimodal, basic multimodal, hybrid, and composition functions.First, the proposed RTH will be compared with the considered algorithms using the ten functions of the IEEE CEC 2020 55 .The characteristics of CEC 2020 functions are presented in Table A2.Their 3D maps are shown in Fig. 11 for the unimodal function (F1), Fig. 12 for the basic functions (F2-F4), and Fig. 13 for the composition functions (F8-F10).
All the algorithms are run 30 times for each function.The results for these functions for D = 15 are presented in Table 4 and.These include the best, the main, and the standard deviation (StD).The fitness evolutions during the optimization process for each function are presented in Fig. 14.
The three last rows show the scores of each algorithm.It can be seen from the L1 row that the proposed RTH can find the optimal solution for five functions from ten, followed by MGO by three functions, then HGS, AVOA, and GTO by two functions, where the other algorithms achieved only one optimal solution.From the penultimate row, the proposed RTG can accomplish the best mean results for five problems, followed by MGO by four, then HGS by two.The last row presents the StD results that confirm the algorithm's robustness.Based on these results, the proposed RTH provided the best results in this criterion four times.The evolution of the average cost function is presented in Fig. 13.From this figure, the proposed RTH has fast convergence speed compared to the other algorithms.This can be explained by the unique exploitation and exploitation strategies of the proposed RTH algorithm.
The functions of the CEC 2020 benchmark will be tested with increased search space dimension where D is equal to 20.This allows us to investigate the performance with more complex cases.The results are provided in Table 5.
It can be seen from the L1 row that the proposed RTH can find the optimal solution for six functions, followed by MGO by three, then HGS, AVOA, and GTO by two functions.From the penultimate row, the proposed RTG has got the best mean results for five problems, followed by MGO by three, then HGS by two.From the last row, the StD results of the RTH is the bet by five times, which confirms the algorithm's robustness.Based on these results, the proposed RTH provided the best results in this criterion four times.The evolution of the average cost function is presented in Fig. 15.Similar to the previous case, the proposed RTH has fast convergence speed compared to the other algorithms.
ANOVA is the abbreviation for analysis of variance.ANOVA is frequently used to evaluate equality across different means by comparing variance between groups with variation within groups.This test will be used to assess the results between the optimizers for both D = 15 and D = 20 cases.The best outcome for each case is marked in red.
For the unimodal function (function 01), the ANOVA ranking is presented in Fig. 16.From this figure; the RTH generates the best solutions in terms of accuracy with the lowest variation range.On the other hand, the FO and the HHO cannot find any solution.Hybrid functions (fnc 05-07) ANOVA results are presented in Fig. 18.The results of function 06 are not included because of the exact similarity of the results.These figures confirm the outperformance of the proposed RTH for the three functions compared to the other considered algorithms.The MGO and GTO algorithms performed the second best after the RTH.
Composition functions (fnc 08-10) ANOVA results are presented in Fig. 19.For function 08; the MGO provided the best performance where the proposed RTH for the best result, but its variation range is higher due to some ultimate results far from the best.The performance of the algorithms is much closer to each other except for the FO and the HHO for function 09, with light superiority for the MGO.The RTH provides more stable performance and higher accuracy for function 10, followed by the MGO.IEEE CEC 2022.In this case, the proposed RTH will be compared with the considered algorithms using the 12 functions of the IEEE CEC 2022 56 .The characteristics of CEC 2022 functions are presented in Table A3.Their 3D maps are illustrated in Fig. 20 for the unimodal function (F1), Fig. 21 for the basic functions (F2-F5), and Fig. 22 for the composition functions (F9-F12).
Similar to the CEC 2020 cases, all the algorithms are run 30 times for each function for D = 10 and D = 20.The obtained results for these functions for D = 21 are presented in Table 6.
From the L1 row, the proposed RTH can find the optimal solution for eight from twelve functions, followed by GTO by six functions, then the HGS by five functions, then AVOA by four.This can approve its ability to get the optimal solution and escape from the local solutions.In comparison, FO, AO, and HHO cannot provide any best solution.From the penultimate row, the MGO achieved the best mean results by six functions, followed by the proposed RTH by five functions, then the HGS and GTO by two.This approves the accuracy of the proposed RTH.FO, COOT, AO, HHO, and AVOA algorithms didn't get any mean best results for all the tests.From the last row, the proposed RTH provided the best results by five functions, followed by the proposed MGO by four times.This can approve its robustness.Figure 23 presents the fitness evolution for each function of the CEC 2022 with D = 10.These figures support the findings of the CEC 2020, where the proposed RTH gives an excellent convergence speed.
For D = 20 case, the results are provided in www.nature.com/scientificreports/From the L1, the proposed RTH can find the optimal solution for six from twelve functions, followed by GTO by only three functions, then the HGS by two functions.The increased complexity can explain the decreased number of optimal solutions due to increased search space dimensions.These results approve the RTH's ability to get the optimal solution better and escape from the local solutions even in complicated problems.From L2, the MGO achieved the best mean results by five functions, followed by the proposed RTH by four.This approves the accuracy of the proposed RTH compared to the MGO algorithm.From L3, the proposed RTH provided the best results by five functions, followed by the proposed MGO by three times.This can approve its robustness.Compared to D = 10, the RTH can generally handle problems that include optimization variables better.Figure 24 presents the fitness evolution for each function of the CEC 2022 with D = 20.These figures approve the findings of the CEC 2020 and CEC 2022 with D = 10, where the proposed RTH presents a faster convergence speed.
Figure 25 shows the ANOVA ranking results for unimodal functions for both cases.From this figure, the RTH and the MGO generate the best solutions in terms of accuracy with the lowest variation range for D = 10.However, when the search space dimensions increased, the performance of the MGO decreased where the RTH performance was constant.
The basic functions (fnc 02-05) results are presented in Fig. 26.For fnc 02 for both cases (D = 15 and D = 20), all the optimizers except for the FO and the HHO got the optimal solution where the MGO and RTH provided the lower variation range for both cases.The results of fnc 03 show that the proposed RTH cannot resolve this problem effectively.The MGO and the HGS provide the best performance for this case.For fnc 04, all the optimizers provided near performance with a slight advantage to the MGO.The same comments can be extracted for fnc 05.
Hybrid functions (fnc 06-08) ANOVA results are presented in Fig. 27.The results of the fnc 06 confirm the outperformance of the proposed RTH.All the optimizers provide similar performance for the other two functions with an advantage to the MGO, HGS, and RTH algorithms.
Composition functions (fnc 09-12) ANOVA results are presented in Fig. 28.For function 09, the MGO and the RTH provided the best performance for both cases.For the other functions, the performance of the  algorithms is much closer to each other, with light superiority for the FO and the HGS for fnc 10, the GTO and the RTH for fnc 11, and the MGO and the COOT for fnc 12.
These mathematical test functions are used to elaborate the performance of the proposed RTH algorithm compared to other new and robust algorithms.The statistical results are summarized in Fig. 29.
From the results obtained, the RTH algorithm can provide very competitive performance for solving various optimization problems.The accuracy and stability of this proposed algorithm have been tested compared to the prementioned algorithms.The convergence speed has been checked and confirmed.After mathematically approving its performance, the proposed RTH algorithm will be tested with some published papers for several engineering tests.

Engineering optimization problems
As mathematics and engineering real-world complex examples, the proposed RTH algorithm will be used to solve seven real-world complex engineering problems.
(a) Optimal design of I-Shaped beam: this issue, which tries to reduce the vertical deflection of the beam, is another common engineering optimization problem 34 .It satisfies the cross-sectional area and stress restrictions simultaneously under specific loads, as explained in Fig. 30.This optimization problem can be expressed as follows: where the optimization variables are: the width of the flange (x 1 ), the height of the section (x 2 ), the thickness of the web (x 3 ), and the thickness of the flange (x 4 ).
From Table 8, the proposed RTH obtains the best solutions for solving this optimization problem compared to other published algorithms.In addition, the PFOA provided similar performance to the proposed RTH, but the RTH achieved the stop criteria after 300 iterations, proving better convergence speed and accuracy.( 9) 10 ≤ x 1 ≤ 50, 10 ≤ x 2 ≤ 80, 0.9 ≤ x 3 ≤ 5, 0.9 ≤ x 4 ≤ 5 www.nature.com/scientificreports/(b) Three-bar truss design: In this example, a 3-bar planar truss is considered 34 , as represented in Fig. 31.The objective function and its constraints are presented as follows Table 9 provides the statistical results of these algorithms.It is clear that the best results are obtained using the proposed RTH algorithm and the AHA optimizer.However, the required number of iterations is only 670, (10)    www.nature.com/scientificreports/much lower than the other reported algorithms, including the AHA, which requires 15,000 iterations to achieve the optimal value.
(c) Design of a tubular column: this problem illustrates the structure of a uniform column with a tubular cross-section that can support a compressive load at the lowest possible cost 34 , as explained in Fig. 32.
The two design variables in this issue are the mean diameter of the column (x 1 ) and the thickness of the tube (x 2 ).The column is constructed from a material having a yield stress of y and an elasticity modulus of E. The following equation represents this optimization problem: www.nature.com/scientificreports/This problem has already been tackled using several algorithms.The best results of these algorithms, as well as the suggested RTH, are shown in Table 10.Based on these results, the proposed RTH achieved the best results with the lowest number of iterations (300).( 11)  71 .As shown in Fig. 33, the weight of the speed reducer must be reduced under 11 restrictions in this optimization problem.Face width (x 1 ), the module of teeth (x 2 ), number of teeth in the pinion (x 3 ), length of the first shaft between bearings (x 4 ), length of the second shaft between bearings (x 5 ), the diameter of first shafts (x 6 ), and diameter of the second shafts (x 7 ) are the seven variables in this issue.This problem can be modeled as     The results of RTH and other published algorithms are reported in Table 11.Among these algorithms, the RTH has the lowest number of iterations, while its metrics in terms of best, mean, and StD are better than those other algorithms.
(e) Piston lever: the essential goal of this problem is to identify the piston components (x 1 ), (x 2 ), (x 3 ), and (x 4 ) by reducing the oil volume when the piston's lever is raised from 0° to 45°3 4 .This problem is presented in Fig. 34, and it can be modeled as follows: (12) f (x) = 0.7854x 1 x 2 2 3.3333x x 5 ≤ 8.3, 2.9 ≤ x 6 ≤ 3.9, 5 ≤ x 7 ≤ 5.5     www.nature.com/scientificreports/According to the obtained results in Table 12, the SNS performs better than the proposed RTH algorithm in terms of mean and StD.However, achieving the stop criteria took many iterations (5000), whereas the RTH requires only 130 iterations.f)Corrugated design: This task seeks to reduce the weight of a corrugated bulkhead in a chemical tanker 66 .The optimization variables are width (x 1 ), depth (x 2 ), length (x 3 ), and plate thickness (x 4 ).Its mathematical model of this problem can be expressed as follows: (13)       www.nature.com/scientificreports/ The statistical results of this problem are provided in Table 13.Although AOA achieved the best results, these results were only obtained after 40,000 iterations, unlike the proposed RTH for the stop value after 129 iterations, ( 14)   which approves its solving speed.g)Design of tension/compression spring: tension/compression spring design challenge aims to decrease the weight of a tension/compression spring, as detailed in 78 and Fig. 35.Minimum deflection, shear stress, surge frequency, outside diameter limitations, and design factors all play a role in this problem.The optimization variables include the mean coil diameter (x 1 ), the wire diameter (x 2 ), and the number of active coils (x 3 ).This problem can be expressed as follows:

Proton Exchange membrane fuel cell parameters' identification
Fuel cells are recent energy-generation devices that produce efficient power by directly using fuel (mainly hydrogen) and oxygen while emitting no pollutants or noise (water, heat, and electricity).Because of its benefits and wide range of applications, proton exchange membrane fuel cells (PEMFCs) that employ polymers as electrolytes are becoming more prominent 85,86 .Fuel cell modeling is a crucial stage that enables the investigation and improvement of its performance 87 .The PEMFC model is built from nonlinear differential equations that explain internal chemical processes.These equations contain various empirical parameters that must be accurately identified to improve model accuracy.As a result, identifying these parameters is essential for creating an accurate model.

PEMFC model.
The PEMFC output voltage can be provided as follows 88 (15)  where T represents the operating temperature, P H2 and P O2 are the partial pressures of the hydrogen and the oxygen, respectively.They can be calculated as follows: where R ha and R hc represent the vapor humilities of both anode and cathode, P a and P c present the inlet pressures of both anode and cathode (atm), A is the electrode surface (cm 2 ), i is the current of the FC (A), and P H2O represents the water vapor saturation pressure(atm).
(      The activation voltage losses (V act ) can be obtained as where ξ 1 , ξ 2 , ξ 3 , ξ 4 denote semi-empirical parameters; C O2 is the concentration of oxygen at the cathode's surface (mol.cm−3 ).It can be calculated as follows: The ohmic losses (V ohm ) can be obtained as where R c is the resistance of the connectors, and R m is the resistance of the membrane.R m can be calculated as where l and A m represent the membrane thickness (cm) and surface (cm 2 ), respectively, and ρ m represents the membrane-specific resistivity (ohm × cm).ρ m can be obtained as where λ is the membrane material's water content.
Objective function.The sum square error (SSE) between the measured (V data ) and model output data (V data ) will be used as an objective function.The objective function can be constructed as follows: where N represents the data size, and x represents a vector containing seven unknown parameters.
The FC data are compared with those generated by the model, and the fitness value is calculated based on the error between them.The model is developed in MATLAB script, and the data are loaded from the Excel sheet file.The identification is an iterative process that updates the candidate solutions at each iteration by sending them to the MATLAB script that includes the FC model and simulates it after that, and then generates the fitness value.This process repeated until the last iteration.
Results.The suggested RTH algorithm will be used to extract the seven unknown characteristics of three PEM fuel cells: the NedStack PS6, the BCS 500W, and the SR-12 500W.Table 15 87 provided the accurate values of the parameters, testing operating conditions, and measurement data for the tested PEMFC types.The upper and lower limits of the empirical parameters are presented in Table 16.The results are compared to those published for the same FC types.
Table 17 shows the comparative findings for the BCS 500 W type, Table 18 for the NedStack PS6, and Table 19 for the SR-12 PEM 500 W. The comparison findings are primarily based on sum square error (SSE).
Based on the findings reported in Table 17, the suggested method-based RTH performed the best of all the reported techniques in this comparison to extract the parameters of the BCS 500W.The proposed RTH algorithms obtained the minimal fitness function (SSE) by 1.14 × 10 −2 .On the other hand, the mean fitness values were similar to the best value.The RTH provides the lowest STD values compared to the different algorithms (9.37 × 10 −8 ), indicating its robustness.Figure 36 depicts the experimental and estimated voltage and power curves of BCS 500W using the suggested RTH algorithm.The estimated voltage and power curves match the experimental curves.These curves demonstrate the proposed algorithm's accuracy in deriving the best BCS 500W parameters.
Similar to the BDC 500W case, the suggested method-based RTH performed the best of all the cited methods in this comparison to extract the parameters of the NedStack PS6.The proposed RTH best result is 2.058, similar to the mean values.The RTH provides the lowest StD values compared to the different algorithms (1.79 × 10 −6 ), approving its robustness.Figure 37 depicts the experimental and estimated voltage and power curves of NedStack PS6 using the suggested RTH algorithm.The estimated voltage and power curves match the experimental curves.
Similar to the previous cases, the suggested method-based RTH performed the best of all the cited methods in this comparison to extract the parameters of the SR 500W.The proposed RTH best result is 0.5607, similar to the mean values.The RTH also provides the lowest StD value by 1.79 × 10 −6 , approving its robustness.Figure 38 shwos the experimental and estimated voltage and power curves of SR 500W using the suggested RTH algorithm.
The obtained results for the standard test functions evaluated the performance of the proposed RTH algorithm for both exploitation and exploitation phases compared to other MAs.From these results, the proposed RTH has excellent exploitation and exploration abilities.In addition, its convergence speed has been approved.Then, the proposed RTH efficiently has been evaluated for several real-world applications.Its performance has been compared to several published algorithms for each problem.The achieved results confirm its ability to handle these problems.Finally, a PEMFC parameters extraction has been performed using the proposed RTH algorithm.The results have been compared to other well-known algorithms to approve its performance.The results achieved were excellent.Hence, its performance has been elaborated and approved for various types of optimization problems.

Conclusion and future works
This paper has proposed a novel metaheuristic optimization algorithm named the red-tailed hawk (RTH) algorithm to solve various optimization tasks and problems.The proposed RTH is inspired by the red-tailed hawk's hinting behaviors of a predatory bird.A mathematical model has been developed to replicate the behavior of red-tailed hawks.The main contribution of this paper is to propose a new optimizer that has high robustness and fast convergence speed when solving various optimization problems.RTH's performance was firstly evaluated

Figure 1 .
Figure 1.Behavior of red-tailed hawk during hunting.

Figure 2 .
Figure 2. Behavior of red-tailed hawk during high soaring stage.

Figure 3 .
Figure 3. Behavior of red-tailed hawk during low soaring stage.

( 8 )Figure 4 .Figure 5 .
Figure 4. Evolution of the direction coordinates as a function of iterations.

Figure 6 .
Figure 6.Agents' evolution as a function of the iterations for a 2-dimensional problem.

Figure 8 .
Figure 8. Convergence curves of the unimodal functions.

Figure 36 .
Figure 36.Voltage and power curves of BCS 500W using the proposed RTH algorithm.
(k) − V model (k, x)) 2 (26) x = ζ 1 ζ 2 ζ 3 ζ 4 R c b . Due to its unique updating mechanism gets closer to the optimal solution without falling into local ones with a fast convergence rate.To better understand the characteristics of the proposed RTH and its basics, the RTH concepts are presented in Table1.

Table 3 .
standard test functions results.The best results are marked in bold.

Table 4 .
CEC 2020 benchmark test function results with D = 15.The best results are marked in bold.

Table 5 .
CEC 2020 benchmark test function results with D = 15.The best results are marked in bold.

Table 6 .
CEC 2022 benchmark test function results for D = 10.The best results are marked in bold.

Table 7 .
CEC 2022 benchmark test function results for D = 10.The best results are marked in bold.

Table 8 .
Comparison of RTH results with the other algorithms for the I-shaped beam problem.

Table 9 .
Comparison of RTH results with the other used algorithms for the three-bar truss design problem.The best results are marked in bold.

Table 10 .
RTH results and the other used algorithms for tubular column designing problem.The best results are marked in bold.

Table 11 .
Comparison of RTH results with other algorithms for speed reducer designing problems.

Table 12 .
Comparison of RTH results with the other used algorithms for the piston lever problem.The best results are marked in bold.

Table 13 .
Comparison of RTH results with the other used algorithms for the corrugated bulkhead design problem.The best results are marked in bold.

Table 14
compares the RTH statistical results with other MAs.The RTH algorithm solves this problem in only 130 iterations, much less than the other considered algorithms with better results.

Table 14 .
Comparison of RTH results with the other used algorithms for tension/compression spring design problem.The best results are marked in bold.V act is the activation voltage losses, V ohm is the ohmic voltage losses, V con represents the concentration voltage losses, and E Nernest denotes the thermodynamic potential voltage (the Nernst voltage).E Nernest can be calculated as follows: where

Table 15 .
The characteristics of the considered PEMFCs.

Table 17 .
The extracted BSC 500W FC parameters.The best results are marked in bold.

Table 18 .
The extracted NedStack PS6 FC parameters.The best results are marked in bold.

Table 19 .
The extracted SR-12 500W FC parameters.The best results are marked in bold.