Improving Vector Evaluated Particle Swarm Optimisation Using Multiple Nondominated Leaders

The vector evaluated particle swarm optimisation (VEPSO) algorithm was previously improved by incorporating nondominated solutions for solving multiobjective optimisation problems. However, the obtained solutions did not converge close to the Pareto front and also did not distribute evenly over the Pareto front. Therefore, in this study, the concept of multiple nondominated leaders is incorporated to further improve the VEPSO algorithm. Hence, multiple nondominated solutions that are best at a respective objective function are used to guide particles in finding optimal solutions. The improved VEPSO is measured by the number of nondominated solutions found, generational distance, spread, and hypervolume. The results from the conducted experiments show that the proposed VEPSO significantly improved the existing VEPSO algorithms.


Introduction
Multiobjective optimisation (MOO) problems involve the simultaneous minimisation/maximisation of multiple objective functions, which usually conflict with each other. Due to the conflict between objective functions, a single solution could not satisfy all objective functions. Hence, MOO problem usually results in a set of tradeoffs or nondominated solutions. The vector evaluated particle swarm optimisation (VEPSO) [1] algorithm has been widely used to solve MOO problems [2][3][4][5][6][7]. As an example, VEPSO algorithm has been implemented in solving DNA sequence problem by minimising four objective functions, namely, measure , similarity, continuity, and hairpin, and two constraints, namely, melting temperature and GC content [7]. Compared to DNA sequence design using binary particle swarm optimization which produces single set of DNA sequences [8], VEPSO is able to generate several sets of good DNA sequences which fulfil the four objective functions and two constraints.
The VEPSO algorithm is adapted from the vector evaluated genetic algorithm (VEGA) [9], in which each swarm optimises one objective function by using the best solution from another swarm as a guidance. However, the VEPSO suffers from performance drawback. Therefore, it is improved by redefining the selection of the guidance from nondominated solution, known as VEPSOnds [10]. Although VEPSOnds has shown better performance than conventional VEPSO, the VEPSOnds suffers from weak performance in terms of lacking solution distributions and convergence to the true Pareto front.
Other than VEPSOnds, there are various MOO algorithms which used nondominated solution to guide particle in finding the optimum solutions for MOO problem. For example, in Multiobjective particle swarm optimisation (MOPSO) algorithm [11,12], all nondominated solutions are separated into groups according to their location in the objective space. A guiding solution for each particle is then randomly selected from the group containing the fewest solutions. Besides, in nondominated sorting PSO (NSPSO) algorithm [13], which uses the main mechanism of the nondominated sorting fenetic algorithm-II [14], each particle is guided by a nondominated solution that is randomly 2 The Scientific World Journal selected using the niche count and the nearest neighbour density estimator. A nondominated solution is selected based on binary tournament selection for the purpose of guiding the other particles in the optimised MOPSO (OMOPSO) algorithm [15]. Additionally, Abido [16] introduces the use of two nondominated solutions, which are called the local set and the global set. The guide is selected based on the nearest distance in objective space between each particle and each member of the nondominated solution of both sets.
Noticeably, most particle swarm optimisation-(PSO-) based MOO algorithms, including conventional VEPSO and VEPSOnds, only use one solution as the particle guide. In particular, in VEPSOnds, particles from a swarm will be guided by the nondominated solution which has the best fitness at one objective function. Thus, the particles may guide the searching with limited information about the other objective functions during the optimisation process. Therefore, VEPSOnds can be further improved by using more than one nondominated solution as particle guide. In this context, this improved VEPSO algorithm will use the best solution from all swarms as guidance during the optimisation process.
The next section of this paper explains the particle swarm optimisation (PSO), the conventional VEPSO, VEP-SOnds algorithm, and the proposed VEPSO algorithms. The following section presents the experimental work and the description of the benchmark test problems and performance measures and the discussion of the results. The final section concludes the proposed technique and discusses few possible future works.

Multiobjective Optimization
where ⃗ = { 1 , 2 , . . . , } is the decision variable vector which represents the possible solution, is the number of objectives, and ∈ R → R is the objective function. { , ℎ } ∈ R → R are the inequality and equality constraint function, respectively. The Pareto optimality concept is defined as follows.  ⃗ ( → ). Nondominated solution is also known as Pareto optimal solution. Definition 3. The set of nondominated solutions of a MOO problem is known as Pareto optimal set, P.
Definition 4. The set of objective vectors with respect to P is known as the Pareto front, PF = { ⃗ ( ⃗ ) ∈ R | ⃗ ∈ P}. PF for a two-objective problem is illustrated as the black circles in Figure 1.
The motivation of MOO is to find as many nondominated solutions as possible according to the objective functions and constraints. However, it is possible to have different solutions which map to the same fitness value in objective space. Therefore, it will be more challenging to find more nondominated solutions.

Original Particle Swarm Optimisation Algorithm.
Particle swarm optimisation (PSO) is a population-based stochastic optimisation algorithm introduced by Kennedy and Eberhart [17]. This algorithm finds an optimal solution using a method inspired by the social behaviour of birds flocking and fish schooling. In the PSO algorithm, an individual is known as a particle, and it holds the possible solution to the optimisation problem, given its position. A particle explores the search space, looking for a better solution with respect to the objective functions defined by the optimisation problem. The search process requires the particle to compare its current position with the best positions that it and the whole swarm have found, so that all particles collaborate with each other.
The PSO algorithm is shown in Algorithm 1. Consider a minimisation problem in which a swarm of particles are flying around in an -dimensional search space, each with a position ( = 1, 2, . . . , ; = 1, 2, . . . , ) representing the possible solution. At initialization stage, all particles are randomly positioned in the search space with random velocity, V ( ). Subsequently, the objective fitness → ( ) of each particle is evaluated based on the objective function for ( ). After that, the particle's best position, Best ( ), is set to its initial position. Additionally, the swarm's best position, Best( ), is the best Best ( ) among all particles, as in (2), where is the swarm of particles In the search process, the algorithm will iterate until the maximum number of iterations is reached. Within an iteration, the velocity and position of each particle are updated using (3) and (4), respectively, where is the constriction factor and is the inertia weight. The 1 and 2 are both random numbers ranging from zero to one. The 1 and 2 are the cognitive and social constants, respectively, which control the attraction of the Best ( ) and Best( ). Then, the ⃗ ( ) for each particle is evaluated again. After updating the fitness, the new position of particle is compared with Best ( ), and the more optimal of the two is saved as Best ( ). Next, the Best( ) is updated as well with the best among all Best ( ), as in (2). When the search process ended, the Best( ) will then represent the best solution found for the problem by this algorithm.

Vector Evaluated Particle Swarm Optimisation Algorithm.
The VEPSO algorithm, introduced by Parsopóulos and Vrahatis [1], uses the multiswarms concept from the VEGA algorithm [9]. Each swarm optimises one objective function using the Best( ) from another swarm. In the VEPSO algorithm, the Best ( ) which has the best fitness with respect to the th objective is the Best( ) for the th swarm, as in (5) The flow of the VEPSO algorithm is given as in Algorithm 2. For problem with objective functions, VEPSO algorithm is similar to that of the PSO but some processes are repeated for all -swarm and nondominated solutions are recorded in an archive. However, the velocity update is reformulated and it is given in (6). Note that the particles in the th swarm will fly using Best ( ) where is defined by (7). The sharing of Best( ) between swarms is illustrated in Figure 2: The nondominated solutions are recorded in an archive after the objective functions are evaluated. In the recording process, the fitness → ( ) of each particle is compared to all others, before it is compared to the nondominated solutions in the archive, using the Pareto optimality criterion, so that the archive only contains nondominated solutions. At the end of the computation, all nondominated solutions are the possible solutions to the problem.

The Improved VEPSO Algorithm by Incorporating Nondominated Solutions.
In the search process of conventional VEPSO, as in Figure 3(a), particles from a swarm are optimised using the Best ( ) from another swarm that has the best fitness at the objective function optimised by the other swarm. However, based on the velocity update of conventional VEPSO in (5)  unless there is a Best ( ) that has better fitness than that at the -objective. Consequently, in a two-objective MOO problem, the Best 1 ( ) of the first swarm is not updated when particle in the first swarm has found a solution with equal fitness at the first objective and better fitness at the second objective. Thus, particles from the second swarm will be guided toward the Best 1 ( ).
Due to this limitation, Lim et al. [10] have introduced an improved VEPSO algorithm by incorporating nondominated solutions (VEPSOnds). In VEPSOnds, as specified by (8), the Best ( ) is still the solution with best fitness at -objective function but is selected from the set of nondominated solutions and not from all Best ( ) of the -swarm where is a nondominated solution and P is the set of nondominated solutions in the archive.
This improvement is illustrated in Figure 3(b) where the Best ( ) is always the best solution with respect to -objective function because the other objective functions are considered as well. Hence, particles from the second swarm can converge faster towards the Best 1 ( ), which is a nondominated solution. As a result, better quality of Pareto front is obtained. From an algorithm perspective, the VEPSOnds is similar to the conventional VEPSO except that (5) in Algorithm 2 is replaced with (8).

The Improved VEPSO Using Multiple nondominated
Leader. Based on the results of VEPSOnds [10], this algorithm suffers weak performance in obtaining solutions that has a weak diversity performance where the solution distributions along the Pareto front are not well distributed. Besides, in comparison to other state-of-the-art MOO algorithm, the VEPSOnds also has a problem in convergence where the obtained solution is far distant from the Pareto front. This weak performance could possibly be caused by the fact that particles in each swarm are guided by one Best ( ) only so the obtained solutions do not well diverse to the other objective functions.
Thus, the use of nondominated solutions to enhance the VEPSO algorithm can be further improved by the use of multileader concept in this work. According to (6), which is the velocity equation of the VEPSO, the particles of a swarm are guided by its Best( ) and another swarm's Best( ). For example, as shown in Figure 4(a), the particles from the second swarm optimise the second objective function using Best 1 ( ) only, which may not be the solution that has the best fitness with respect to the second objective function. Thus, this original mechanism of VEPSO may limit the convergence rate of the algorithm. Therefore, an improved VEPSO algorithm is proposed by including Best 2 ( ) as additional guidance to optimise both objective functions, as shown in Figure 4 Hence, the general velocity equation of this improved VEPSO is formulated as in (9) where for each , 2 , and 2 are independent constant and random values, respectively. In addition, from (9), as compared to the improved VEPSO at previous section, the particles will search toward the nondominated solutions which located at different end of the Pareto front. Therefore, The Scientific World Journal the diversity performance of the algorithm is expected to be better as the search area is wider, rather than a single point. Because the improved VEPSO algorithm uses multiple nondominated solutions as particle guides, or leaders, this algorithm is called VEPSO using multiple nondominated leaders (VEPSOml). Also, a polynomial mutation mechanism from NSGA-II [14] is used to modify particle positions at some probability. By mutating the position of some particles out of the locally optimal solution, this mechanism broadens the search for a globally optimal solution. In this study, the position of one out of every fifteen particles is mutated in the algorithm. Therefore, the complete VEPSO algorithm using multiple nondominated leaders is shown in Algorithm 3.

Performance
Measure. MOO algorithms face difficulty in converging to and distributing the nondominated solutions over the true Pareto front, PF . Hence, the algorithm performance is measured by the quality of the obtained Pareto front, PF . Several performance measures are used for comparison to highlight any improvement in the proposed algorithm.
The number of solutions (NS) measured will calculate the total number of nondominated solutions found by an algorithm. The best algorithm, by this measure, gives the most nondominated solutions. A more advanced measure uses the generalized distance (GD) [18], which is a popular measure of convergence [14]. This performance measure first evaluates the average distance between the true Pareto front and the one obtained by the algorithm. Equation (10) is used to compute the average distance, with a smaller value corresponding to a better performance. Then, the minimum distance of a nondominated solution from the true Pareto front is calculated using (11) In addition, SP [14] is a commonly used measure of the diversity performance, or the distribution of nondominated solutions [14] is used. Equations (12), (13), and (14) evaluate the diversity performance, as measured by SP. The and are the Euclidean distances between the boundary solution and the nondominated solutions returned by the algorithm and the true Pareto front, respectively. The Euclidean distance between two solutions can be calculated using (13). Thus, SP actually measures the average distance of one solution and of the next solution to all nondominated solutions returned by the algorithm as well as two boundary solutions in the true Pareto front. Hence, it is desirable that the Pareto front returned by the algorithm produces a small SP: Additionally, the hypervolume (HV) [19] measures the area (in a two-objective problem) or the volume between a reference point, and the Pareto front with respect to the nondominated solutions obtained by the algorithm, as illustrated in Figure 5. Thus, it is desirable that the Pareto front returned by the algorithm produces a large HV.

Test Problems.
Because different features in MOO problems are responsible for decreasing the likelihood of obtaining Pareto front with good convergence and diversity, the standard test functions with well-defined true Pareto fronts are important for testing optimisation algorithms. Five test functions from Zitzler et al. [20] (ZDT) are used here for this reason. The ZDT test problems have two objectives and are formulated with one feature in each problem. ZDT5 is not used because it is binary coded, whereas this work focuses on real-value problems. During testing, the GD, SP, and HV measure require the true Pareto front for the ZDT test problems, the standard database generated by jMetal (http://jmetal.sourceforge.net/problems.html) is used for this purpose. Additionally, all test problems used here are set up as recommended by [20].   VEPSOml is conducted without the use of polynomial mutation as to clarify that the polynomial mutation is not the sole reason for any performance improvement. Thus, this experiment compares the conventional VEPSO and the VEPSOnds without mutation against two different variations of VEPSOml: VEPSOml1 is the VEPSOml without mutation and VEPSOml2 is the VEPSOml with mutation, respectively.
All improved VEPSO algorithms are compared to the conventional VEPSO algorithm. Hence, similar parameters are used for all experimented algorithms which are listed in Table 1. In addition, the archive size is set to 100 solutions and is controlled by removing the nondominated solutions with the smallest crowding distance [14]. Each test problem is simulated for 100 runs on each algorithm to obtain statistical results for a fair comparison because the convergence and diversity performance varies in each run. Table 2 lists the performance of each algorithm on the ZDT1 test problem. In the NS measure, the number of nondominated solutions significantly increases in all improved algorithms. Under the GD measure, VEPSOnds performs approximately 10 times better than conventional VEPSO. However, under the same measure, VEPSOml1 shows a more dramatic improvement, performing approximately 100 times better than VEPSO, as the concept of multiple nondominated leaders shows its benefit in finding more accurate solutions. Additionally, when the polynomial mutation is included, as in VEPSOml2, the GD performance improved much better at about 600% as compared to the conventional VEPSO. Under the SP measure, VEPSOnds also gives a significant improvement in performance. Meanwhile, the VEPSOml1 and VEPSOml2 show significant improvement in diversity performance as compared to the VEPSOnds. This shows the significance of using more than one nondominated solution which diversify the search toward the nondominated solutions at different end. The above mentioned improvements are supported by the higher HV measures when compared to the conventional VEPSO, which indicates that they return better Pareto fronts. Figure 6 shows plots of the nondominated solutions with the best GD measure returned by each algorithm tested on ZDT1. From the first plot, it is clear that the nondominated solutions obtained by VEPSO are far away from the true Pareto front, which explains the poor performance of this algorithm for this test problem. In addition, the nondominated solutions are distributed unevenly, and so VEPSO has a larger SP value. Meanwhile for all the improved VEPSO algorithms, their nondominated solutions fall very close to the true Pareto front. However, VEPSOnds produces a distribution of nondominated solutions that contain empty spaces along the true Pareto front, which results in higher SP value as compared to the other improved VEPSO algorithms. Table 3 lists the performance of the algorithms on the ZDT2 test problem. The average number of nondominated solutions found by VEPSOnds1 slightly improves over the number found by VEPSO, but VEPSOnds2, VEPSOml1, and VEPSOml2 greatly improve over VEPSO by this same than the reference point, . On the other hand, all improved algorithms are able to create a hypervolume, especially the VEPSOml2 which produce the largest hypervolume. Figure 7 displays the nondominated solutions, plotted for each the best GD measure obtained for each algorithm using the ZDT2 test problem. The first plot shows that VEPSO returns nondominated solutions that are far from the true Pareto font and poorly distributed. Although VEPSOnds and VEPSOml1 return a low GD measure, the number of nondominated solutions is found to have low value, which is clearly displayed in the second and third plots, respectively, of Figure 6. In fact, there is only one nondominated solution found by both algorithms which falls exactly on the true Pareto front and yields a GD value of zero. On the other hand, the fourth plot of Figure 6 shows that VEPSOml2 returns the nondominated solutions that converge nicely and are well distributed over the true Pareto front. Besides, the nondominated solutions found by VEPSOml2 distributed evenly which yield a good SP value. Table 4 lists the performance of the algorithms on the ZDT3 test problem. All improved VEPSO algorithms are able to find more nondominated solutions than the conventional VEPSO algorithm. In addition, the performances of the improved VEPSO algorithms, with respect to convergence, improve on conventional VEPSO, while VEPSOml2 shows the greater improvement. However, by the SP measure, the VEPSOnds algorithm performs worse than the conventional VEPSO algorithm. However, although the SP value of the  conventional VEPSO algorithm is better, the superior convergence of the VEPSOnds algorithm maintains its performance advantage. In contrast, both improved VEPSO algorithm using multiple nondominated leaders show better SP measure than the conventional VEPSO, which strengthen the hypothesis that using multiple nondominated leaders will improve diversity performance. In addition, the HV value of the conventional VEPSO algorithm is smaller than of all improved algorithms which suggest that the improved algorithms have better performance. Figure 8 displays the nondominated solutions, plotted for the best GD measure obtained for each algorithm using the ZDT3 test problem. The nondominated solutions returned by the conventional VEPSO algorithm were distributed equally but not well converged with respect to the true Pareto front.
On the other hand, the nondominated solutions from all improved VEPSO algorithms are well converged with respect to the true Pareto front. However, the nondominated solutions returned by the VEPSOnds algorithm are denser at the upper left of the Pareto front, which causes the increase in its SP value. In contrast, the nondominated solutions obtained by both VEPSOml algorithms are equally distributed over the Pareto front and yield better SP value. Table 5 lists the performance of the algorithms on the ZDT4 test problem. The average number of nondominated solutions obtained by VEPSO is relatively low, while all improved VEPSO algorithms found most of the solutions. In this test, the conventional VEPSO algorithm produced a very large GD value due to the multimodality feature in the test problem, and so the improved VEPSO algorithms  clearly performed better in this respect. However, the diversity performance of nondominated solutions returned by conventional VEPSO is small compared to the VEPSOnds algorithm. Once again, the use of multiple nondominated leaders in VEPSO algorithms could diversify the search and result in better diversity performance. Additionally, all algorithms produce a hypervolume from the reference point, and all improved algorithms return larger HV values than the conventional algorithm. Figure 9 displays the nondominated solutions, plotted for the best GD measure obtained for each algorithm using the ZDT4 test problem. The first plot shows that VEPSO converges to the Pareto front but only manages to obtain a single nondominated solution. The VEPSOnds algorithm not only converges to the Pareto front but also returns a diverse set of nondominated solutions. On the other hands, both VEPSOml also returned the nondominated solutions with good convergence but they are not well distributed as compared to the VEPSOnds, in this case. Thus, the VEPSOnds shows better HV value as compared to the VEPSOml. Table 6 lists the performance of the algorithms on the ZDT6 test problem. All algorithms find a similar number of nondominated solutions. In the GD measure, all algorithms are capable of returning the nondominated solutions that converge well to the Pareto front. On the other hand, both VEPSOml1 and VEPSOml2 algorithms outperform the conventional VEPSO and VEPSOnds algorithm in the GD measure. In addition, the SP and HV values for each algorithm are similar. However, the VEPSOml2 algorithm shows superiority in getting the minimum SP value and average HV value. As can be predicted from the similar quantitative performance of the algorithms on ZDT6, the plot of nondominated solutions returned by each algorithm is very similar, especially in convergence performance, as shown in Figure 10. The plots do show that VEPSO has slightly less diversity compared to VEPSOnds and VEPSOml2 because of some small gaps in coverage along the middle of the Pareto front. On the other hand, the VEPSOml1 shows weak distribution of nondominated solutions over the Pareto front. In contrast, the nondominated solutions found by VEPSOnds and VEPSOml2 completely cover the true Pareto front and are spaced out equally.
As seen from the results of all the test problems, the VEPSO algorithms using multiple nondominated leaders shows more improvement in terms of convergence and diversity of the nondominated solutions found than the VEP-SOnds. The additional leader, specifically the nondominated solution with respect to the objective function optimised by a swarm, not only guides the particles to optimise the objective function with respect to the swarm. It also increased the search area because all leaders used to guide the particles are located at the different end of the Pareto front.

Analysis of the Number of Particles.
This experiment analysed the performance of the VEPSOml2 algorithm with various numbers of particles. Similar parameters from the previous experiment were used except for the total number of particles as it is equally divided into two swarms; the total number of particles was varied to be 10, 30, 50, 100, 300, 500, and 1000. Figure 11 shows plots of the performance measures for each benchmark problem against the total number of particles.
The VEPSOml2 algorithm performance improved as the number of particles increased. The performance of the VEPSOml2 algorithm was sufficient when there were 100 particles computed for 250 iterations, which corresponds to 25000 function evaluations. However, the performance of the VEPSOml2 algorithm exhibited better results when the total number of particles was increased. Unfortunately, when the number of particles is increased, the algorithm requires more computational effort to solve the problem.

Analysis of the Number of Iterations.
This experiment investigated the performance of VEPSOml2 for various numbers of iterations. The number of iterations was fixed to be 10, 30, 50, 100, 300, 500, 1000, 3000, 5000, and 10 000. Meanwhile, the other parameters were kept the same as in the previous experiment except that the number of particles, which were divided equally among swarms, was fixed to 100 divided equally between all swarms. Figure 12 plots the performance measures for each benchmark problem against the number of iterations.
As expected, the performance of VEPSOml2 is improved when the number of iterations was increased. When 100 particles were used, the VEPSOml2 algorithm started to yield acceptable results when there were 500 iterations, which is equivalent to 50000 function evaluations. However, if computational cost is not critical, the VEPSOml2 algorithm could use 3000 iterations because the performance saturated after this value.

Benchmarking with the State-of-the-Art Multiobjective Optimisation Algorithms.
For benchmarking, the VEP-SOml2 algorithm was compared to four other state-ofthe-art MOO algorithms: nondominated sorting genetic algorithm-II (NSGA-II) [14], strength Pareto evolutionary algorithm 2 (SPEA2) [21], archive-based hybrid scatter search (AbYSS) [22], and the speed-constrained multiobjective PSO (SMPSO) algorithm [23]. All algorithms only computed 25000 function evaluations, and the archive size was set to 100 for fair comparison. The population size for NSGA-II was set to 100 for optimisation. The Simulated Binary Crossover (SBX) operator was used with crossover probability = 0.9. The polynomial mutation [24] operator was also used with mutation probability = 1/ . Meanwhile, the distribution indices for both operators were set to = = 20. The parameters in SPEA2 were set the same as in NSGA-II. The The Scientific World Journal population size for AbYSS was set to 20 and the pairwise combination parameters 1 and 2 were both set to 10. In addition, the polynomial mutation parameters in AbYSS were also set similarly as in NSGA-II and SPEA2. algorithms, its minimum GD was still the smallest among them. However, the SP was, on average, better than NSGA-II. Interestingly, the HV measure of VEPSOml2 was as good as those of the other algorithms. Table 8 presents the performance measure of the algorithms for the ZDT2 problem. The VEPSOml2 was sufficiently competitive at obtaining a reasonable number of solutions. In the GD measure, on average, the VEPSOml2

18
The Scientific World Journal  algorithm was as good as the other algorithms, but SMPSO had greater performance. Surprisingly, the VEPSOml2 algorithm was able to obtain a better minimum GD measure than the SMPSO algorithm. Additionally, the SP measure of the VEPSOml2 algorithm was better than those of the other algorithms except SMPSO. All algorithms had similar HV values, but VEPSOml2 yielded the best HV performance. The performance measures for the ZDT3 problem for all algorithms are listed in Table 9. Both SMPSO and VEPSOml2 were unable to obtain the maximum number of solutions consistently for all 100 runs but still yielded solutions within a reasonable range. Noticeably, the average GD measure for VEPSOml2 was the largest among all algorithms. However, the diversity for VEPSOml2 was similar to that of the others. Moreover, although the HV value of VEPSOml2 was the smallest, it still yielded a very large HV. Table 10 presents the performance measures for the algorithms for the ZDT4 problem. VEPSOml2 faced great challenges from the multiple local optima featured in this problem, where it cause the algorithm to obtain a very small number of solutions. Additionally, the convergence and diversity of VEPSOml2 were bad, as indicated by the very large GD and SP values. As expected, the HV performance was also very poor because the multiple local optima feature is one of the natural weaknesses of PSO-based algorithms [26,27].  Finally, the performance measures for the ZDT6 problem for all algorithms are listed in Table 11. VEPSOml2 algorithm was inconsistent in obtaining the maximum number of solutions. Moreover, the convergence and diversity measures for VEPSOml2 were significantly larger than those for the other algorithms. However, the VEPSOml2 algorithm was able to obtain the minimum GD value. Additionally, the HV performance for VEPSOml2 was relatively weak, on average, but its maximum HV value was the largest of all the algorithms.
An overall performance comparison for state-of-theart algorithms against VEPSOml2 was investigated in this experiment. In some cases, the VEPSOml2 algorithm yielded better results than some of the other algorithms.

Conclusions
Most PSO-based MOO algorithms, including conventional VEPSO and VEPSOnds, only use one solution as the particle guide. Thus VEPSOml is proposed in this study where the particles are guided by multiple nondominated solutions while retaining the unique information shared between swarms that are inherent in conventional VEPSO. Five ZDT test problems were used to investigate the performance of the improved VEPSO algorithm based on the measures of the number of nondominated solutions found, the generational distance, the spread, and the hypervolume. The proposed VEPSOml algorithm obtained a higher-quality Pareto front as compared to conventional VEPSO and VEP-SOnds. The VEPSOml2 algorithm that included polynomial mutation has exhibited further improvement for most of the performance measures.
Using more nondominated solutions as particle guides yielded faster convergence performance improvements, especially for the ZDT1, ZDT2, and ZDT3 test problems. The use of more than one leader reduced the risk of trapping at local Pareto front. In future, the success of using two leaders motivates the investigation of the use of more than two leaders during the optimisation process.