A multi-objective optimisation approach with improved pareto-optimal solutions to enhance economic and environmental dispatch in power systems

This work implements the recently developed nth state Markovian jumping particle swarm optimisation (PSO) algorithm with local search (NS-MJPSOloc) awareness method to address the economic/environmental dispatch (EED) problem. The proposed approach, known as the Non-dominated Sorting Multi-objective PSO with Local Best (NS-MJPSOloc), aims to enhance the performance of the PSO algorithm in multi-objective optimisation problems. This is achieved by redefining the concept of best local candidates within the search space of multi-objective optimisation. The NS-MJPSOloc algorithm uses an evolutionary factor-based mechanism to identify the optimum compromise solution, a Markov chain state jumping technique to control the Pareto-optimal set size, and a neighbourhood’s topology (such as a ring or a star) to determine its size. Economic dispatch refers to the systematic allocation of available power resources in order to fulfill all relevant limitations and effectively meet the demand for electricity at the lowest possible operating cost. As a result of heightened public consciousness regarding environmental pollution and the implementation of clean air amendments, nations worldwide have compelled utilities to adapt their operational practises in order to comply with environmental regulations. The (NS-MJPSOloc) approach has been utilised for resolving the EED problem, including cost and emission objectives that are not commensurable. The findings illustrate the efficacy of the suggested (NS-MJPSOloc) approach in producing a collection of Pareto-optimal solutions that are evenly dispersed within a single iteration. The comparison of several approaches reveals the higher performance of the suggested (NS-MJPSOloc) in terms of the diversity of the Pareto-optimal solutions achieved. In addition, a measure of solution quality based on Pareto optimality has been incorporated. The findings validate the effectiveness of the proposed (NS-MJPSOloc) approach in addressing the multi-objective EED issue and generating a trade-off solution that is both optimal and of high quality. We observed that our approach can reduce \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼6.4% of fuel costs and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼9.1% of computational time in comparison to the classical PSO technique. Furthermore, our method can reduce \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼9.4% of the emissions measured in tons per hour as compared to the PSO approach.

In the power industry, recent research has been driven to focus on issues related to reviewing industrial design or operations in order to lower pollution and emissions to the environment.These ecological effects usually appear from thermal power stations 1-3 due to the growing public awareness of environmental protection.This is due to the fact that these power stations consume significant amount of fuel and, most recently, the world has seen rising fuel costs along with its environmental impacts.The emission dispatching option, which aims to reduce both emissions and fuel costs, is a desirable short-term alternative.This strategy has attracted a lot of interest recently [4][5][6][7][8] because it only needs a slight modification to the fundamental economic dispatch to account for emissions.
In power systems, one of the most important optimisation problems is known as Economic Environmental Dispatch (EED) that is sometimes referred to as Economic Emission Dispatch (EED).The EED's primary objective is to ascertain which configurations of components within a power generation system result in the most efficient generation of power.The catch is, however, that the solution needs to be both economically viable and environmentally friendly in order to be considered acceptable.The significance of the EED problem has been growing substantially as people all over the world become more aware of the need to preserve the natural environment.The goal is crystal clear: to simultaneously cut down on the overall expense of fuel and the pollution that it causes to our environment.The PSO has become an increasingly popular method for addressing the EED problem over the course of the past several years.This is mostly attributable to the fact that it is a straightforward method that is both effective and good at locating global optimal solutions to similar optimisation problems.
The ELD problem of power systems has been successfully resolved by using PSO techniques as discussed in section Related work.In the existing literature, PSO algorithms iteratively modify the parameters of a swarm of particles to converge toward the ideal solution by maximizing the distribution of power generation across numerous units while taking into account several restrictions including fuel cost, power demand, and generator limits.To meet the real-time requirements for dynamic power system operation, however, more developments are still required to improve the performance of PSO methods for ELD, particularly in handling larger and more complex power systems, integrating renewable energy sources, taking uncertainties into account, and increasing computational efficiency.Furthermore, we believe that there still exists a gap to concentrate on creating advanced PSO-based strategies that combine PSO's advantages with other optimisation methods to get around traditional PSO algorithms' drawbacks and produce more durable and dependable solutions for the ELD issue in contemporary power systems.
This paper proposes and employs a novel nth state Markovian jumping PSO algorithm with a local search (NS-MJPSOloc) method to solve the economic/environmental dispatch problem.Subsequently, the newly developed algorithm implies the theory of local search capability.Using this capability, the problem search space having multiple optima is thoroughly explored.It is known that the canonical versions of the PSO algorithms are based on global search.Similarly to other multi-objective evolutionary algorithms, an evolutionary factor-based mechanism is used to identify the optimum compromise solution, and a Markov chain state jumping technique is used to control the Pareto-optimal set size.The results of several runs on the common IEEE topology test system are compared to other methods described in the literature.The efficacy and potential of the proposed (NS-MJPSOloc) approach are shown to solve the multi-objective EED problem.The major contributions of this research are as follows: • we propose and employ a novel nth state Markovian jumping PSO algorithm with local search (NS-MJPSO- loc) method to solve the economic/environmental dispatch problem; • the proposed algorithm uses an evolutionary factor-based mechanism to identify the optimum compromise solution; • a Markov chain state jumping technique is used to control the Pareto-optimal set size along with a neighbour- hood's topology (such as a ring or a star) to determine its size; and • the algorithm is implemented for the economic dispatch problem in the domain of power systems.The experimental outcomes of the proposed NS-MJPSOloc approach has been verified on IEEE 30 Bus and 15-unit Systems.
The rest of the paper is organized as follows.A summary of the related work is offered in Section Related work.
The problem statement along with the constraints is explained in section Problem statements.The optimisation problem is formulated in section Problem formulation.The concept of the multi-objective optimisation is elaborated in section The concept of multi-objective optimisation.The proposed multi-objective optimisation algorithm is explained in section The proposed NS-MJPSOloc algorithm.Performance evaluation of the proposed algorithm is discussed in section Performance evaluation.The experimental setup is explained in section Experimental setup and evaluation metrics are given in section Evaluation metrics.The obtained results and findings are illustrated in section Results and discussion.Finally, the concluding remarks along with future research directions are summarized in section Conclusions and future work.

Nomenclature
The list of abbreviations shown in Table 1 and the list of mathematical notations shown in Table 3 are used in the rest of the paper.We believe that these tables will help all readers to quickly understand all the mathematical formulas mentioned in this paper.

Related work
In 2014 Han et al. 9 delved deep into the environmental and economic dispatch of a micro-grid, encompassing diverse energy sources like photovoltaic generation, wind turbines, and more.The study underscored the efficacy of an improved linearly decreasing weight PSO algorithm, emphasizing its theoretical and practical feasibility.In 2016 Tlijani et al. 10 presented an extended version of the conventional DEED (dynamic economic environmental dispatch), aiming to mitigate ramp rate violations across consecutive dispatch periods.The goal was clear: to consistently meet periodic load demands.In 2017 11 illuminated the potential of a multi-objective PSO algorithm, leveraging both the Pareto criterion and fuzzy logic, to address environmental pollution in economic dispatch.In 2018 12 elegantly formulated the power dispatch challenge as a dual-objective optimisation problem.The mission was dual-pronged: simultaneous minimization of fuel cost and emissions.Considering its concealed nature, false data injection attacks (FDIAs) have attracted a lot of attention in the field of cyber-physical power systems (CPPS).Improving CPPS cybersecurity requires an understanding of likely attacker actions.Nonetheless, the majority of FDIA models now in use frequently concentrate on the implications of attacks or the effects of attackers alone.In response, a unique multi-objective stealthy FDIA strategy is presented in 13 within the framework of an AC grid model.In order to maximize the impact of the attack and minimize tainted measurements while keeping stealth, the suggested attack model is presented as a multi-objective optimisation problem.Additionally, in order to improve the attack vector's generation efficiency, a novel representation mechanism is proposed to characterize the positions and parameters of injected states 13 .
In 2019 14 proposed a refined version of PSO to tackle the EED conundrum of thermal electric power units.The innovative Space Reduction strategy was employed to pinpoint the Pareto optimal solution within the designated search space.In 2020 15 integrated DE with Quantum PSO (QPSO) to address the short-term EED challenge of microgrids.In another study, Mehrpour et al. 16 focused on the dynamic load and emission dispatch in daily cycles, especially considering the potential impacts of renewable energy sources.In 2022 17 showcased the Perfectly Convergent PSO (PCPSO) for addressing combined economic and multiple emissions dispatch challenges.The study meticulously considered the ramifications of various pollutants, employing cubic functions with seven price penalty factors.The history of EED when viewed through the perspective of these research articles, presents a picture of continuous innovation and development.Microgrids and other forms of renewable energy are only two examples of how this industry has expanded its scope in pursuit of solutions that are good for the economy while also being friendly to the environment.The application and refining of PSO have been the constant thread spinning throughout this exploration, establishing its strength as a solid tool for EED difficulties.Environmental/economic dispatch (EED) problems have been resolved in a variety of ways [4][5][6][7][8][18][19][20][21][22] . In genral, there are three methods to resolve the EED issue.The first approach involves considering the level of emissions as a constraint with a tolerable limit 4 .However, demonstrating the relationship between cost and emissions in this formulation is quite difficult.
The second method addresses the emission as a separate objective beside the traditional cost objective [5][6][7][8] .However, the EED problem was simplified to a single objective function by linearly combining the two objectives or only taking into account a single objective at each stage of optimisation.inevitably, this method finds marginally non-dominated solutions and necessitates more runs than the required number of Pareto-optimal solutions.In power systems, the Economic Emission Dispatch (EED) problem is a well-known constrained multi-objective optimisation problem.It strives to meet a variety of operational requirements while concurrently minimizing expenses and emissions.Even though a lot of solutions have been devised to deal with it, however, this problem remains difficult and challenging because of the unpredictable and inconsistent nature of renewable energy sources (RES) like wind and solar, in particular, when they are integrated into the system.The authors in 23 present a novel constrained multi-objective optimisation technique called CMOPEO-EED, with the goal of improving EED performance in the presence of renewable power generation.
The third method treats fuel cost and emission as dual objectives at the same time.For the EED issue, multiobjective search heuristic and fuzzy membership-based optimisation approaches have been explored 18,19 , and 20 , whereas, algorithmic approaches do not offer a logical framework for guiding the search to the Pareto-optimal front, and it is extremely hard to expand these methods to include other objectives.These methods need a lot of computational power and, subsequently, take a lot of time.Numerous non-dominated alternatives can be observed in a single run using strategies based on multi-objective genetic algorithms, as described in 21,22 .Premature convergence is a problem with genetic algorithm-based methods, and the method described in 21 requires a lot of computational work because of the ranking process that takes place throughout the fitness assignment phase.
The PSO approach provides an adaptable and diverse strategy to enhance and evolve global and local exploration capabilities, despite genetic algorithms and heuristic approaches.Compared to the genetic algorithm, it typically produces faster convergence rates 24 .PSO has been applied with remarkable popularity in the past decade to a variety of power system problems, such as the economic power dispatch issue 25,26 , and 27 .It has been shown and documented that PSO has the potential to deal with non-smooth and non-convex economic power dispatch problems 26,27 .subsequently, the fuel cost was the only factor taken into account for optimisation when the problem was defined as an ordinary dispatch problem.
In order to obtain an edge of optimal solutions, it is predominantly necessary to redefine global and local best persons when switching from a uni-objective to a multi-objective PSO.There is a set of non-dominated solutions rather than an absolute global best in multi-objective PSO.Additionally, there could not be a single local best individual for every swarm particle.In a multi-objective domain, selecting the global and local best for steering the swarm particles turns into a challenging problem.
Numerous real-world problems can be formulated as multi-objective optimisation problems (MOPs), in which it is necessary to simultaneously optimize several, frequently conflicting objectives.Finding a set of solutions that are unable to be enhanced in one area without compromising another is the aim of addressing a mixed optimal problem (MOP).MOEAs, or multi-objective evolutionary algorithms, have become a popular and successful method for handling MOPs 28 .As MOEAs may produce roughly optimum solutions in a single run and do not require specific assumptions like continuity or differentiability, they are used in the majority of the related works.Furthermore, these methods are based on randomised search algorithms that draw inspiration from Darwin's theory of natural selection.Although MOEAs have obvious benefits, it is vital to remember that they require a large number of objective function evaluations, which could make them unfeasible for some applications requiring a lot of computational power.
The computational complexity of evaluating the objective functions and the flexibility of the input parameters are the two main determinants of the computational cost of MOEAs at each iteration 29,30 .Since many MOPs in the first instance are unable to be represented analytically, evaluations must be done through laborious simulations.Though there have been several attempts to reduce the MOEA execution time by utilizing the knowledge collected throughout the search process, these approaches usually lead to a consensus regarding the accuracy level of the final findings.Applications using high-dimensional spaces, including multi-objective programming (MOP) problems with four or more objectives or those with hundreds or thousands of decision variables, fall under the second instance.Large-scale MOPs and many-objective optimisation problems (MaOPs) are the terms used to describe these kinds of MOPs 31 .They considerably increase MOEAs' runtime 32 .Furthermore, large population size is needed for some applications in order to improve accuracy or cover more ground in the search space [33][34][35] .While the majority of MOEAs function within an expected polynomial time for the size of their population, storage capacity limitations can provide difficulties 36,37 .
Recent reports indicate that PSO variants including NS-MJPSOloc implementation require only a small amount of work to resolve power system issues.In order to handle the EED problem with thermal dispatch and various renewable energy sources, Wang and Singh 38 proposed a fuzzified MOPSO.The method offers a fuzzification process for choosing the world's top person while considering the world's top as an entire area rather than simply a single spot.On the other hand, each particle is maintained with a single local best solution.This will reduce search efficiency and is against the multi-objective optimisation principle.By breaking down the original optimisation problem into smaller problems, Kitamura et al. 39 presented a modified MOPSO to optimize an energy management system.This method, however, faces serious limitations when there is a significant interaction between the constraints in several sub-problems.In their MOPSO-based solution to the congestion management issue, Hazra and Sinha 40 demonstrated how to reduce both costs and congestion at the same time.The sigma method 41 is used in this strategy to identify the ideal local particle guide.However, the use of the www.nature.com/scientificreports/sigma values heightens the PSO's already considerable selection pressure.In some situations, such as multifrontal difficulties, this may lead to early convergence.In order to calculate the contributions of generators to the transmission system, a vector-evaluated PSO (VEPSO) was suggested and studied 42 .Depending on each objective independently, VEPSO selects portions of the future generation from the older generation.However, choosing people who excel in one area without considering the other areas raises the issue of eliminating those with average performance, who might be very helpful in finding compromise solutions 43 .In 44,45 , the authors have proposed solutions for enhancements in evolutionary algorithms and optimization techniques 44 outlines a strategy for improving the effectiveness of SAEAs (Surrogate-assisted evolutionary algorithms) utilizing unevaluated solutions.A crucial element of MOEA, the offspring generator, has received little attention and lacks a design concept.In order to overcome this problem, regularity evolution (RE), an offspring generation model for MOEAs, is introduced in 45 .The authors in 46 proposed an adaptive algorithm for control system optimizations to address the issue of selecting optimal starting values for connection weight parameters in MPIDNN (multivariable PID neural network).The authors suggested a constrained optimization problem and an adaptive population extremal optimization-based MPIDNN method called PEO-MPIDNN, which minimizes exponential time and system errors.
In 47,48 , the authors have proposed novel approaches in multi-objective optimization and algorithm efficiency within machine learning environments.For online sequential learning machines, in 47 a multi-objective model selection approach is suggested to enhance the target output error, control quality, and channel equalization.To build channel equalization models and accomplish classification selection and equalization, the proposed method makes use of feedback compensation and adaptive equalization control.The authors in 48 suggest ASDNSGA-II, an improved fast NSGA-II based on an adaptive crossover technique and unique congestion strategy.The proposed approach balances the convergence and variety of the decision and object spaces, hence improving the selection strategy. 49discusses the application of federated learning (FL) and blockchain technology in IIoT.To lower energy usage and application latency, the study focuses on FL Aware Multi-Objective Modeling in Decentralized Microservices Assisted IIoT Systems.To optimize workload allocation and application delay, the Blockchain-Enabled FL Algorithm Framework (DLEBAF) is designed with three strategies: (i) deadline-efficient task sequencing and scheduling (DETS), (ii) latency-efficient task scheduling (LETS), and (iii) energy-efficient task scheduling (EETS).Table 2 shows the advantages and disadvantages of the existing works.We believe the data in Table 2 will help our reader to quickly identify what is missing in the current literature and what is further needed to improve the efficiency of EED problems in power systems.

Problem statements
In order to solve the EED problem, the fuel cost and emission objective functions should be minimized while adhering to a number of equality and inequality constraints.In section Problem objectives, we describe the objectives of EED problem.In section Problem constraints, various constraints are explained.Finally, the EED problem is mathematically described in section Problem formulation.Table 3 describes a list of all mathematical notations and their brief description.We believe that this table will help all readers to quickly understand all the mathematical formulas.
The problem in question is typically formulated as follows.

Problem objectives
We consider two objectives for the EED problem in our optimisation problem i.e. fuel cost and carbon emissions.The mathematical foundations of both objectives are described in subsequent subsections.

Smooth cost minimization
The goal of the ELD problem is to generate electricity while satisfying all system constraints at the lowest possible cost per unit (fuel cost in US dollars).The cost of each generator is calculated using the quadratic function in the smooth or simplified ELD problems.The equality and inequality constraints are the fundamental limitations of the simplified cost functions.The total cost of fuel in US dollars per hour, F(P), is represented as: The objective function is a minimisation problem and is given by: where C represents the total generation cost and F i represents the objective function of generator i.Furthermore, N is the total number of generators, a i , b i , and c i are the ith generator's cost coefficients, and P Gi is the generator's real power output of the ith generator.P Gi , which is described as follows.
The index of a particle is represented by n, while the number of units or generators is represented by D. For example, P i1 denotes the power produced by unit 1.The dimension of the problem in this function is [Population times Number of systems units].However, on the basis of instances, appropriate constraints are taken into account.

Carbon emission minimization
In order to develop a mathematical model for emission reduction in a power system, the challenges including power generation, emission factors, and various constraints are being considered.A generalized mathematical model for minimizing emissions is provided here: Decision variables: In the emission reduction model, we define the following decision variables.P Gi : Power generated by the ith generator.E i : Emission factor caused by (per unit generated power) by the ith generator.Hence, the total sum of emissions produced by the entire committed unit (power system) is formulated as follows: (4) Solving economic dispatch problem using Cuckoo Search Algorithm Testing the method on western Algerian electrical power system The work lacks an investigation of the computational complexity of the CSA as well as any discussion on its convergence features.Furthermore, the paper fails to address the CSA's ability to effectively manage uncertainties or fluctuations within the power system, a critical aspect for practical implementation.The research does not address the constraints or potential disadvantages associated with the MOSSA approach, such as its computing complexity or susceptibility to variations in problem size and complexity.Furthermore, the analysis conducted utilising the MOSSA technique lacks comprehensiveness in terms of evaluating the statistical significance and robustness of the acquired data.www.nature.com/scientificreports/Emission pollutants that are produced by fossil-fuelled thermal units, comprising sulphur oxides SO x and nitro- gen oxides NO x , could be analysed independently.However, in order to provide an illustration, the aggregate ton/h emission E(P i ) among these particles is equal to Eq. ( 5): as described in 5 and 8 .
whereas, α, β, γ , δ and ǫ are various parameters describing different emission coefficients for each generator in the power system.

Problem constraints
We consider three constraints for the EED problem in our multi-objective optimisation problem i.e. power generation (lower and upper limits), power stability (production meets demand), and security constraints.The mathematical foundations of these constraints are described in subsequent subsections.

Power generation constraints
Every single generator's realistic power output is constrained by its upper and lower limits for reliable operation in the manner described below: The generated power at time t must be equal to the demand from the total loads side as given below: (5)  www.nature.com/scientificreports/

Power stability constraints
The power stability constraint means that the total amount of generated power P Gi must be equal to the total amount of demand P D plus the actual transmission line losses P loss .Hence, mathematically it is given by: In practice, there are numerous ways to determine transmission line losses, including the power flow and the B matrix technique.Another approach has been used in our implementation, and it entails solving the load flow problem with equality bounds on both reactive and real power at each bus in the way that is described below.
where i refers to a particular bus index, designating the bus at where power is produced (P Gi ) , consumed (P Di ) , or where the magnitude of the voltage (V i ) and phase angle (i) are obtained.However, j indicates an additional bus index that designates a different bus in the network.The equations take into account the contributions from many buses (from j = 1 to j = n) , as indicated by the summation ( ) over j, where n represents the total number of buses in the system.Subsequently, P G i indicates the power generated by source i, P Di is the power consumed by load i, V i is the voltage magnitude at bus i, G ij is the conductance between buses i and j.Similarly, B ij determines the susceptance between buses i and j.Eventually, θ i and θ j illustrates the voltage phase angles at buses i and j, respectively.Thus, the actual power transmission losses can be measured with the following equation.
In this equation, the power loss can be calculated by combining the inputs of every element (indexed by k).Moreover, g k is used to denote the electrical conductance within the kth line connecting bus i to bus j, and N represents the total number of transmission lines.

Security constraints
Security and integrity constraints can be expressed empirically to define the limits or requirements that must be met for safe and reliable operation.Hence, the transmission line loading S is constrained by its upper limit for secure operation as follows: It is important to operate any generator within its minimum and maximum capacity.This should be noted that going above the capacity limit will compromise the system's security, reliability, and dependability.

Problem formulation
This problem ought to be mathematically represented as a non-linear constrained MOPs by combining both constraints and objectives simultaneously.

Minimize(P G ) in terms of:
The equality constraint g is dependent on the variables P G , F(P G ) , and E(P G ) .This condition makes sure that specific requirements are fulfilled.However, the parameter h is represented as the inequality constraint that is further dependent on P, which ensures that some specific requirements are met.

The concept of multi-objective optimisation
In the real-world, simultaneous optimisation of multiple objective functions occurs in many problems.These functions usually have incommensurate, predominantly opposing, and contradictory objectives.Considering such competing objective functions, multi-objective optimisation produces a set of optimal solutions instead of just one.Numerous solutions are optimal because no one can be said to be superior to any other with regard to all objective functions.Pareto-optimal solutions are exactly these optimal approaches that have been referred to, as discussed in 59 .
The two competing solutions, x 1 and x 2 , associated with a MOP, can either dominate the other or else, depend- ing on the situation.If the subsequent two conditions are met, a solution x 1 to a minimising problem dominates x 2 regardless of compromising clarity: (8) In case, any of the preceding conditions are compromised, the outcome x 1 fails to dominate over the solution x 2 .The non-dominated solution is x 1 if it dominates the solution x 2 .The Pareto-optimal set is defined as the set of solutions that are non-dominated across the whole search space.The Pareto optimal front is another name for this set.

General overview
In recent years, study on evolutionary approaches illustrates that population-based algorithms are well-suited for solving multi-objective optimisation problems.They can also be effectively used to overcome many of the limitations of traditional single objective strategies, including their sensitivity to the Pareto-optimal front's shape and their requirement of numerous runs to find various Pareto-optimal solutions.
Traditionally, the main goal of a MOP algorithm is to maintain population diversity in the set of Pareto optimal solutions, besides steering the search towards the Pareto-optimal front.PSO has been inspired as a promising social heuristic method with an adaptable and diverse strategy to improve and modify the capabilities of both local and global exploration in recent years.In contrast, considering that there is not a definite global best in multi-objective PSO, only a set of non-dominated solutions, transforming classical single objective-based to a multi-objective-based PSO necessitates redefining global and local best candidates.Additionally, there could not be a single local best solution with every swarm particle.In a multiobjective space, selecting the global and local best for steering the population turns into a challenging problem.
The challenges of evolving a multi-objective variant of the standard PSO are resolved by the proposed approach, which involves a process for choosing between the global and local best candidates.It is important to note that the suggested MOPSO technique has been applied with remarkable success to a number of challenging standard test problems in the field of multi-objective optimisation 60 .

Core concepts and notions
The brief descriptions, definitions, and terminologies of the proposed NS-MJPSOloc algorithm are given as follows: • Particle, current position: X i (t) is a candidate solution, where d is the total number of optimized parameters, and is represented by an d-dimensional vector.X i,d (t) = [x 1,d (t), ..., x n,d (t)] , where X i,d (t) represents the position of the particle ith with respect to the dimension d, or the value of the dimension parameter d in the candidate solution ith, describes the particle ith at time t.• Swarm size, population: S(t), represents a distinct set of n particles in time t, whereas, S(t) = [x 1 (t), ..., x n (t)] • Velocity vector, V i (t) : The parameter identified as velocity adjusts how each particle V i,d (t) moves in the d-dimensional search space.In order to identify optimal or near-optimum solutions, it coordinates the ways to exploit and explore the swarm in the search space.
, where the parameter t represents time.• N states Markov Jumping: It is a mathematical illustration of a series of occurrences or states where the next state is solely dependent upon the current state.Furthermore, it relates to a particular type of Markov chain with a set of states represented by N. • Evolutionary factor, E f : In order to automatically adjust the inertia weight and acceleration coefficients, an evolutionary factor was developed that determines four specific evolutionary stages such as convergence, exploitation, exploration, and jumping out in each generation [61][62][63][64][65] .The evolutionary factor is able to consider data about population distribution.In this paper, we define four states using the evolutionary factor.The following expression represents the mean distance of each particle in the whole swarm: wherein Eq. ( 17), S represents population size and D represents the dimension of the search space and Eq. (18),P dg denotes the global best particle among the P d (i) , P d(min) and P d(max) are particles with minimum and maximum distances, respectively.• Inertia weight, ω(t) : The control parameter known as inertia weight, or ω(t) , serves to determine how sig- nificantly the preceding velocities influence the current velocity.As a result, it influences the balance between the global and local exploring capacities.A large inertia weight is recommended for the early stages of the search process to improve global exploration, however, a smaller inertia weight is proposed for the later stages to improve local exploration.It is assumed that E f is significantly large in the jump state and small in the convergence state.The evolutionary factor E f and the inertia weight ω share several characteristics.As a result, the mapping ω(E f ) is defined as follows: (15) ∀i ∈ {1, 2, . . ., N obj } : • Neighbourhood K i , Lbest: A subgroup of neighbouring particles, labelled K i , this group is known as the neighbourhood of particle i.Based on the neighbourhood's topology (such as a ring or a star), its size can be determined.Within its defined neighbourhood search space, K i indicates the indexes of neighbouring particles for each particle i. Considering the Lbest strategy, the velocity update model for particle i at time t + 1 can be expressed as follows: Whereas, the Lbest i (t) indicates the particle's best position with the best fitness value observed in its imme- diate neighbourhood.

Computational stream
In the computational steam of the proposed NS-MJPSOloc algorithm, we have n number of particles with D dimensional parameters and neighbourhood K i .It can be described in the following Algorithm 1, Algorithm 2, and Algorithm 3.
Initialize particles' positions, velocities, personal bests, and local bests while termination condition is not met do for each particle i do Update particle's velocity: Update particle's position: ( made to encourage particles to move more quickly in the direction of the best particle overall.A greater c 2 accelerates convergence by amplifying the effects of the particle's individual best position and the global best position.However, a smaller c 1 minimises the influence of the particle's previous velocity, enabling it to more successfully explore new regions.The values for c 1 and c 2 in the aforementioned case are set to c 1 (δ(4)) = 1.8 and c 2 (δ(4)) = 2.2, respectively.These parameters are selected to allow escape from local optima and to encourage quick progress towards the particle that performs best globally.
• Exploration state In the exploration state, the focus is on analysing new optima while avoiding getting stuck in existing ones.In this phase, c 1 is used with a greater value, whereas c 2 is used with a relatively smaller value.By favouring the past velocity of the particle c 1 over the impact of the global best position c 2 , these changes aim to promote individual exploration.The values for c 1 and c 2 in the precise instance are set at c 1 (δ(3)) = 2.2 and c 2 (δ(3)) = 1.8, respectively.By letting the particle rely more on its own velocity to explore new areas of the search space, these values are set to encourage individual exploration.To establish a balance between exploration and exploitation, acceleration coefficients are frequently dynamically adjusted throughout the optimisation process.• Exploitation state The main purpose of the exploitation phase is to employ the local knowledge that the particles have while looking for the best solutions in the potential region.Typically, this state comes after the exploration state and before the convergence state.In this state, c 1 serves a substantially bigger value, whereas c 2 is used with a relatively smaller value.With these adjustments, the particle's own best position (local infor- mation) will be given more weight, but the impact of the global best position will still be taken into account.
The values for c 1 and c 2 in the precise instance are set at c 1 (δ(2)) = 2.1 and c 2 (δ(2)) = 1.9, respectively.These values were selected to strike a compromise between using each particle's local information and investigating the prospective region as a whole.• Convergence state Finally, in the convergence state, the swarm has a tendency to form dense clusters and become close to the overall best solution.The proposed neighbourhood's topology is shown in Fig. 2.However, there is a chance of early convergence, in which the particles may become trapped in less-than-ideal solutions and stop further exploring the search space.The acceleration coefficients c 1 and c 2 are modified to address this problem and maintain search diversity.The values for c 1 and c 2 in the exact scenario are set at c 1 (δ(1)) = 2 and c 2 (δ(1)) = 2, respectively, in the convergence stage.These parameters have been selected to promote exploration and preserve search diversity within the swarm, while also pointing the particles in the direction of the present global area.The particles strike a compromise between exploration and exploitation by setting c 1 and c 2 to the same value.This method enables the swarm to carry on searching and maybe find better solutions, preventing premature convergence.

Experimental setup
In this paper, the goal was to effectively solve nonlinear constrained optimisation problems using the NS-MJP-SOloc algorithm.In order to do this, an additional method of neighbourhood awareness (as shown in Fig. 2) was added to the proposed algorithm to evaluate the viability of potential solutions at each stage of the search.This process is used in the study to guarantee that the non-dominated solutions are both optimal and compliant with the set of constraints.The proposed NS-MJPSOloc algorithm can successfully handle a challenging optimisation problem with several objectives given various constraints.Table 4 describes various experimental parameter settings for all compared algorithms.A desktop PC Corei5 with 3.30GHz, 8GB RAM, Windows 10 Enterprise was used to carry out the evaluation of the proposed NS-MJPSOloc approach.Furthermore, the MATLAB R2018b application is used for algorithm development, data analysis, visualization and production of results.A collection of parameters from Table 4 with generators data in Table 5 was used during the optimisation runs.
The maximum number of generations was established at 1000, the number of particles was decided to be 100 and the tests were repeated 30 times.The results shown in various tables are averaged over these multiple runs.The datasets used in the experiments were taken from previous published sources 63 and online open source code repositries (https:// github.com/P-N-Sugan than).Furthermore, the power system, bus unit, and other datasets generated and/or analysed during the current study are publicly available in the Github repository, and can be accessed at (https:// github.com/ evgen ytsyd enov/ ieee1 18_ power_ flow_ data).Moreover, various IEEE bus systems for power systems used within the experimental work are publicly available online.
A threshold of 25 solutions was placed in place to keep the Pareto-optimal set at a suitable size.The neighbourhood-best set has a maximum size of 10 solutions as well.In order to manage and regulate the size of these www.nature.com/scientificreports/solution sets, a neighbourhood strategy is used if the number of non-dominated Pareto-optimal solutions in either the global best set or the neighbourhood-best set exceeds their respective boundaries.

Evaluation metrics
The proposed algorithm is compared with other closest rivals using several performance evaluation metrics such as fuel costs (measured in US dollars per hour), electricity/power generation through each generator (MW), carbon emission (measured in tons per hour), and computational time (measured in seconds).

Results and discussion
The proposed NS-MJPSOloc technique was used for the first time on the typical IEEE 30-bus 6-generator test system to assess its efficacy.This system is frequently referenced in the literature, and numerous strategies have already been tried on it with promising results.Figure 3 shows the single-line diagram of the IEEE test system, and references 5,8 provide in-depth information about the system.Table 5 lists the values of the fuel costs and the emission coefficients related to the generators.
Three cases with different levels of complexity were taken into account in order to demonstrate the efficacy of the proposed algorithm.These examples were chosen to show how well the algorithm performed under various conditions and tasks.• With P Loss : In this scenario, transmission losses P loss are considered along with power balancing and genera- tion capacity constraints.• Without P Loss : Here transmission losses P loss are ignored, while power balancing and generation capacity constraints are taken into consideration.• All constraints: In this case, all relevant constraints were considered as described in the problem constraints, i.e., Section Problem constraints.
The aim of this study is to provide a smart solution to the ELD problem with line flow and emission limits.Considering the IEEE 30-bus system, the EED problem in power systems is investigated and evaluated with several plausible assumptions.Furthermore, to achieve the EED schedule with the least amount of generation and cost  63 .
of the generating units, this paper uses a variety of PSO variants, such as GA, EP, PSO, and DE.The performance of a newly developed PSO variant i.e., NS-MJPSOloc is also investigated and compared with other methods.The IEEE 30-bus system's generating characteristics are listed in Table 5.
The study rendered employing a variety of intelligent algorithms mentioned earlier, and the results of the transitional cost analysis for the IEEE 30-bus system are shown in Table 6.The results shown that our proposed NS-MJPSOloc technique can approximately 2.3%, 3.0%, 2.5%, 2.4%, and 0.4% fuel costs per hour as compared to GA, EP, DE, PSO, and NS-MJPSO techniques, respectively.Tables 7, 8, 10 also show the IEEE 30-bus system's convergence criteria, given that line flow constraints are taken into consideration.Table 10 compares the minimal total production costs attained utilizing the above algorithms for a demand of 283.4 MW.
The data summarized in Table 6 shows that, when compared to the overall minimal production costs obtained by using (EP) technique, the costs obtained in the research are noticeably higher.The (EP) algorithm requires more generations to reach convergence than the other techniques (GA, PSO, and DE).This shows that finding the best solution for (EP) might require additional iterations.
The (GA) needs additional solution time compared to the other methods evaluated in the research.The parameters that the aforementioned methods such as (GA, EP, PSO, and DE) apply at various times are detailed Scaling factor, transition prob ( β, φ) -0.008/0.0016-1.9/1.9 0.9/0.90.9/0.9Time taken (s) for the IEEE www.nature.com/scientificreports/ in Table 7.We observed that the proposed algorithm could save significant computational times as compared to other closest rivals.For example, for a 15-unit bus system, NS-MJPSOloc runs approximately 49%, 12.7%, 9.1%, 7.8%, and 1.9% more efficiently than GA, EP, DE, PSO, and NS-MJPSO algorithms, respectively.Similarly, these values were noted as 66.2%, 33.5%, 17.7%, 40.6%, and 26.3% for a 30-unit bus system.The effectiveness and efficiency of each method in solving the ELD problem are greatly influenced by these parameters.Figures 4  and 5 show the convergence rates of the proposed NS-MJPSOloc's and DE's evolutionary processes throughout a 6-unit bus system, respectively.Furthermore, Figs. 6 and 7 show the results of the NS-MJPSOloc's and DE's optimal dispatch load analysis for the 6-unit power system, respectively.Moreover, the study shows that the (EP) technique works better than other algorithms in terms of reaching a lower overall production cost for the specified demand.To achieve convergence, though, more generations must be produced.The (GA), in contrast, takes longer to solve problems than the other methods.The parameter settings of each method are presumably revealed in Table 4, which explains why each algorithm performed differently in the article.
The generator attributes for the 15-unit system are presented in Table 8.Furthermore, Table 9 provides an overview of the comparison of ELD with line flow constraints using different intelligent techniques for a 15-unit system.We observed that our approach can reduce ∼ 6.4% of fuel costs in comparison to the classical PSO technique.Furthermore, approximately 0.3%, 0.61%, 4.7%, and 0.07% fuel costs can be saved by the proposed   10 represents the actual power generation output of the 15-unit system and the transition cost of the 15-unit system, determined by a variety of intelligent techniques.Table 7 demonstrates the convergence conditions for the 15-unit system concerning emission and line flow constraints using GA, EP, PSO, and DE.A breakdown of the minimum total production costs for demand of 2630 MW using smart techniques like GA, EP, PSO, and DE is shown in Table 10.A description of the cost estimation is given in the Smooth cost minimization section.As  reported in Table 10, our method can reduce approximately 9.4% of the emissions measured in tons per hour as compared to the PSO approach.
The results in Table 7 to Table 8 show that for both the IEEE 30-bus system and the 15-unit systems, the PSO algorithm exhibits rapid convergence and requires less time.Although GA requires more time for convergence, it outperforms other intelligent techniques for the 15-unit system when the total minimum generation cost is taken into account.Table 11 shows the solutions for cost and emission optimized using the proposed multiobjective algorithm.

Conclusions and future work
In this paper, the environmental/economic power dispatch (EED) optimisation problem was addressed using a newly developed, neighbourhood-aware, n-states Markovian jumping PSO (NS-MJPSOloc) algorithm.The (NS-MJPSOloc) algorithm is a variant of the conventional PSO approach that was created specifically to handle MOPs.The fuel cost and environmental impact were two competing objectives, taken into account when formulating the EED problem.The results derived show the significant potential and effectiveness of the proposed method (NS-MJPSOloc) in solving the multi-objective EED problem.The method also demonstrated a significant capacity to generate a variety of evenly distributed solutions within the non-dominated set.While comparing the simulation results, it was obvious that the (NS-MJPSOloc) method outperformed other variants of the PSO algorithm in the diversity and quality of the Pareto-optimal solutions achieved.These findings indicate that the proposed method holds promise in enhancing the optimisation process and facilitating better decision-making in power dispatch optimisation, considering both economic and environmental aspects.The evaluation of the proposed approach demonstrated that it can reduce ∼6.4% of fuel costs and ∼9.1% of computational time, and ∼9.4% of the emissions measured in tons per hour, in comparison to the classical PSO technique.
The following are limitations and potential future directions for this research.www.nature.com/scientificreports/The primary limitation of the proposed version is the complexity is increased exponentially as compared to the simplified framework of the PSO algorithm.The (NS-MJPSOloc) algorithm could potentially be further improved and optimized to increase its effectiveness and efficiency in resolving multi-objective EED problems.Robust mechanisms for the adjustment of various parameters, other neighbourhood topologies, dynamic adaptive mechanisms, constraints handling, multi-objective optimisation metrics, and merging with other algorithms.By following these research directions, the (NS-MJPSOloc) algorithm can be further enhanced, and its adaptability to different scenarios can be improved.
Implementing the proposed approach in a real-world power dispatch environment would provide useful information about its usefulness and efficacy.Furthermore, the incorporation of renewable energy sources into power networks expands the EED problem formulation to consider these sources and their irregular role.This would overcome the hassle created by the integration of renewable energy sources and make the optimisation procedure more analogous to contemporary power systems.Ultimately, these efforts can contribute to more effective and reliable solutions for multi-objective EED optimisation problems, advancing sustainable energy management practices in power systems.In the future, we will extend this work to an effective energy management system within the domain of smart grids. https://doi.org/10.1038/s41598-024-62904-4 jt = P D load,t

Figure 1 .
Figure 1.Illustration of the the jumping-out state.

Figure 2 .
Figure 2. The neighbourhood's topology using four different structures.

Figure 4 .
Figure 4. Convergence of the proposed algorithm over a 6-unit system.

Figure 5 .
Figure 5. Convergence of the DE evolutionary process throughout a 6-unit system.

Figure 6 .
Figure 6.Results of NS-MJPSOloc optimal dispatch load analysis for the 6-unit power system.

Figure 7 .
Figure 7. Using DE, the determined optimal dispatch load for a 6-unit power system.

Table 1 .
List of Abbreviations.

Table 2 .
Advantages and disadvantages of existing works.The research does not address the issue of scalability in relation to the algorithm's performance when applied to larger power systems that encompass a greater number of thermal power plants.The provided analysis lacks a comprehensive examination of the algorithm's computational efficiency, specifically in terms of the duration necessary for convergence to the best solution.

Table 3 .
List of mathematical symbols and notations.

Table 4 .
Parameter settings for all the compared algorithms.

Table 6 .
Method Comparison with closest rivals.

Table 7 .
Algorithm Performance on IEEE 30 Bus and 15-unit Systems.

Table 9 .
Comparison of ELD with line flow constraints using different intelligent techniques for the 15-unit system.

Table 10 .
The optimal solutions for cost and emission are optimized separately (in terms of fuel costs and carbon emissions).

Table 11 .
Solutions for cost and emission optimized using multi-objective.