Location Management in LTE Networks using Multi-Objective Particle Swarm Optimization

Long-term evolution (LTE) and LTE-advance (LTE-A) are widely used efficient network technologies serving billions of users, since they are featured with high spectrum efficiency, less latency, and higher bandwidth. Despite remarkable advantages offered by these technologies, signaling overhead remains a major issue in accessing the network. In particular, the load of signaling is mainly attributed to location management. This paper proposes an efficient approach for minimizing the total signaling overhead of location management in LTE networks using multi-objective particle swarm optimization (MOPSO). Tracking area update (TAU) and paging are considered to be the main elements of the signaling overhead of optimal location management in LTE. In addition, the total inter-list handover contributes significantly to the total signaling overhead. However, the total signaling cost of TAU and paging is adversely related to the total inter-list handover. Two cost functions should be minimized, the first is the total signaling cost of TAU and paging and the second is the total signaling overhead. The trade-off between these two objectives can be circumvented by MOPSO, which alleviates the total signaling overhead. A set of non-dominated solutions on the Pareto-optimal front is defined and the best compromise solution. The proposed algorithm results feasible compromise solution, minimizing the signaling overhead and the consumption of the power battery of a user. The efficacy and the robustness of the proposed algorithm have been proven using large scale environment problem illustrative example. The location management in LTE networks using MOPSO best compromise solution has been compared to a mixed integer non-linear programming (MINLP) algorithm. Location management mobility management entity MME pooling clustering SON Distributed Centralized pooling scheme fuzzy implementation setup LP-CPLEX

to higher requirements on networks, since signaling overhead directly affects the cellular performance. High levels of signaling overhead result in increased computational cost and more power consumption [1]. This could take a form of high load and data traffic in the core network which harms the cellular network reliability. The development in signaling overhead may tackle the problem of data traffic by 50% [4]. In addition, LTE and LTE-A have a flat IP-based heterogeneous architecture that increases the signaling overhead [1]. Accordingly, the research focus was shifted from signaling overhead to handling greater amount of data and increasing speed connectivity. In fact, several aspects could lead to signaling overhead, however, it can be primarily attributed to user mobility.
Signaling overhead is triggered by transmission or reception of data between mobile networks and user equipment (UEs). The performance of signaling overhead can be evaluated by two factors, namely, paging and tracking area update (TAU). From one side, paging is normally actuated by mobility management entity (MME) placing an idle-user within the tracking area (TA) of the system. TA is commonly termed routing area or location area in other applications. TA is a virtual area, which combines a group of cells located within a particular area [1]. The main drawback of TA is the service interruption which occurs when a particular cell gets reassigned to a different TA. In addition, the edge users located between two or more tracking areas (TAs) are subject to ping-pong effect. Tracking area list (TAL) is the state-of-the-art extension of TA proposed by the 3rd Generation Partnership Project (3GPP). Several TAs can be combined to form TAL. Thus, the structure of TA and TAL is nearly similar. However, TAL is characterized by greater flexibility and diversity of cells. MME defines the location of a user via his/her most recent registered location. Accordingly, MME is updated by the last TAL. From the other side, TAU is initiated by the user to the MME, this procedure includes reporting an update. This update governs transferring the user from a cell in one TAL to another TAL. TAL has a prominent role of minimizing the total signaling overhead generated by TAU and paging. Thus, the design of TAL could allow to adapt various TAs and reduce the ping-pong effect. Once the mobile user switches into an idle mode, the signals generated by TAU and paging are diminished through TA or, in other words, via TAL. Both TAU and paging are employed to track the location of the user, and subsequently give a continuous update to the evolved packet core (EPC). Consequently, the total signaling cost of TAU and paging have a major signaling impact on EPC, allowing them to play a vital role for the location management between UEs. An important issue to consider is that the inter-list handover resides on three factors: UEs, probability of the user to move between cells, and the cost of inter-MME. A careful look at the aforementioned discussion reveals a bottleneck effect between the total signaling cost of TAU and paging from one side, and the total cost of the inter-list handover from the other side [5,6]. Thus, the location management addressed in this work is NP-hard, and locating approximate solution is NP-hard. Therefore, this work handles two main factors of optimal location management in LTE, the first factor is the total signaling cost of TAU and paging, and the second factor is the total inter-list handover. The ultimate objective is to minimize the total signaling overhead which is a result of paging, TAU and inter-list handover. The problem is tackled by optimizing TAL, which is successfully achieved by minimizing the total signaling cost of TAU and paging, defined as the first objective, and the total inter-list handover -the second objective, consequently attenuating the total signaling overhead. It should be remarked that TAU and paging are inversely related to the total inter-list handover. In an effort to circumvent the abovementioned challenge, the problem can be divided into twolayers. The two-layered approach can be carried out by mixed integer non-linear programming (MINLP) [7] which divides the problem into two sub-problems. The main disadvantage of this approach using MINLP is that one objective is favored over the other. By consequence, the final solution obtained by MINLP confines the solution on the Pareto-optimal front. Unfortunately, the solution obtained will not guarantee the best compromise solution.
Thus, the main focus is two optimize the two objectives simultaneously. According to the conflict between the aforementioned two objectives, the trade-off can be achieved through a best compromise among a set of non-dominated solutions on the Pareto-optimal front using multi-objective particle swarm optimization (MOPSO) [8]. In this regard, the signaling overhead problem, which is mainly attributed to the total signaling of paging, TAU and inter-list handover is studied for an arbitrary network. For demonstration and numerical results, a large scale network is considered. Numerical results show the robustness for the proposed algorithm of finding the best compromise solution by solving the trade-off between objective one, which is the total signaling cost of TAU and paging, and objective two, which is the total interlist handover. Also, the results guarantee the minimization of the signaling overhead problem to lower levels. In addition, the proposed algorithm significantly reduces the power consumption of the user's device. The contributions of this study can be summarized as follows: 1) The location management problem has been reformulated and approached as a multi-objective optimization problem.
2) The problem is solved as a minimization problem considering the two above-mentioned objectives using MOPSO. 3) A fuzzy based mechanism to extract the best compromise solution has been implemented. 4) The superiority of the solution obtained by MOPSO over MINLP has been proven.
The rest of the paper is organized as follows: Related work and an overview of mobility management techniques is presented in Section II. Section III defines the problem, states the system model and formulates the problem. The description of MOPSO, the related flow chart, and MOPSO implementation is presented in Section IV. Numerical results and the performance of the proposed deployment algorithm compared to MINLP are illustrated in Section V. Section VI includes the summary of the work and the concluding remarks.

II. RELATED WORK
Location management is utilized by various technologies and is popular among many researchers [9,10]. Actually, location management has motivated scholars to investigate the signaling overhead problem. TAU and paging are essential parts of location management [9,11]. Number of stereotype techniques aimed to attenuate the signaling overhead, for instance (velocity-based, timer-based, and movement-based [11][12][13]). The control-plane components, for instance, service gateway (SG) or MME, have been employed to regulate the signaling load [14][15][16]. Two schemes of MME pooling, namely the centralized and distributed MME pooling schemes have been investigated measuring the load in signaling from user mobility [14]. However, the study in [14] did not examine the structure of TA and TAL. As introduced in [15] a pool of SGs assists the TAs in terminating the intermittent connections once UE exits one TA and enters another. The minimization of relocation frequency has been studied through relocating the data by SGs to support high-mobility users [16]. The shortcoming of the approach proposed in [15,16] is that it targets only non-idle users, which could weaken the service quality through SG re-allocation.
The signaling overhead problem is mainly attributed to TAL assignment. The TAL assignment research can be divided into two categories: static and dynamic TAL configuration. TAL design using conventional optimization techniques has been introduced through case studies using a linear programming (LP) CPLEX (LP-CPLEX) optimizer [17,18]. The results in [17,18] showed significantly better results in terms of signaling overhead compared to the static configuration of TAL. The solution proposed in [17,18] has been extended in a new case study with less complex protocol of LP-CPLEX optimizer [19]. Recently, centralized and distributed MME pooling scheme has been employed to attempt addressing the signaling overhead problem [6]. The work in [6] used MINLP optimizer to solve the trade-off between the total signaling cost of TAU and paging and the total inter-list handover. It was found that the centralized MME pooling scheme outperforms the distributed MME pooling scheme in terms of minimum total signaling overhead. It is important to recall that the location management in LTE networks is NP-hard and finding appropriate solution is NP-hard. Therefore, a critical drawback of the study in [6] is that the algorithm used favored one objective over the other one, thus the solution obtained is far from being the best compromise solution.
Particle swarm optimization (PSO) is a global search technique which considers the natural behavior of bird swarming [20]. It has been used effectively to solve: the problem of resource allocation and maximization of the system throughput while keeping the minimum user rate requirement [21], sensor deployment in designing sensor networks [22], near optimal deployment of sensors with Voronoi diagram evaluating the cost function [23], and optimal tuning of fuzzy feedback filter [24]. PSO algorithm introduced in [20,25] is applicable for single objective maximization problems. The multi-objective version of PSO is termed MOPSO. Due to the fact that a single optimal solution is no longer available in multiobjective context, MOPSO is concerned with selecting the best compromise solution from a set of non-dominated solutions [8,[26][27][28]. Moreover, evolutionary techniques have been efficient of solving several recent problems, operational aircraft maintenance routing [29][30][31], directional steering [32,33], and optimization of membership functions of a fuzzy logic controller utilized in control applications [24,28,34,35]. In the next section, the problem of location management in LTE networks is formulated to be subsequently solved by MOPSO.

III. PROBLEM FORMULATION
This study aims to minimize the two main contributing factors to the total signaling overhead problem, namely, the total signaling cost of TAU and paging and the total interlist handover. This can be accomplished by defining the optimal deployment of overlapping TALs. The total signaling overhead problem can be tackled by one of the two popular schemes: centralized pooling scheme [5,14] or distributed pooling scheme [14]. Both schemes associate the distribution of TAL with MME. The main difference between centralized scheme and distributed one is the method of calculating TAU. According to the result of comparison between the two schemes, distributed scheme requires extra signal load from MME relocation, thus, it can be stated that centralized scheme is more efficient than the distributed scheme [5,6]. Therefore, the problem considered in this section utilizes a centralized MME pooling scheme [5,14]. Both schemes are illustrated in Figure 1, where the upper portion depicts the centralized MME pooling scheme, while the lower portion shows the distributed one.

A. Important Notations
Table I presents math and network notation that will be used throughout the paper. Table II lists notation related to the multi-objective particle swarm optimization algorithm. For any N ∈ N, n ∈ [1, N ] indicates that n is an integer number between 1 and N . z ∈ {0, 1} denotes a binary number or more simple z ∈ Z. Throughout the paper,

B. Preliminaries and Design Hypothesis
The cells are allocated within every TAL through the centralized MME pooling scheme. Every TA refers to a single cell, hence, TA and cell will be used interchangeably. The centralized MME pooling scheme is depicted in the upper portion of Figure 1. The model includes two layers of assignments.   : Position of parameter j within particle i, ∀i, j ∈ N v i,j : Velocity of parameter j within particle i, ∀i, j ∈ N jth parameter of ith best global particle, ∀i, j ∈ N P : Number of parameters within 1 particle, P ∈ N The first layer considers the Cells-to-TALs/MMEs assignment associated with the core network. The second layer presents the TALs/MMEs to UEs assignment, which relates to every cell within the system. The first layer of the system evaluates the UEs mobility pattern among the cells, then it allocates the cells in TALs/MMEs with the goal of attenuating the total signaling overhead caused by TAU and paging. In the second layer, a number of TALs is positioned at a portion of UEs termed σ l k . The selection of σ l k defines the usage portion of every single TAL/MME per cell which leads to interference of TALs through the cells. This technique guarantees a broader assortment of cells in the list.
1) Cell/TA-to-TAL/MME Assignment: The allocation of cell/TA-to-TAL/MME in the centralized pooling scheme is the main objective of the presented model. Every TAL refers to a certain one-to-one MME basis. The allocation of Cellsto-TAL/MME is given through a binary decision variable C l k,n ∈ Z for all l, k, n ∈ N, k, n ∈ [1, N ] and l ∈ [1, L]. In order to proceed with the model and problem formulation, it is necessary to present the following example. Example 1. As an illustrative example of C l k,n ∈ Z allocation in the search space, define C l k,n to be a cell-to-TAL assignment. For one list (l = 1), define the maximum number of cells/TAs to ben = 3 with N = 4 MMEs or lists. Accordingly, the total number of possible combinations in matrixL 1 = C l is 3×4, where C l k,n is a binary decision created in the system. One can specify the decision variables in list 1 as Equation (1) in Example 1 shows thatL 1 is symmetric with Also, it demonstrates that number of nonzero rows/columns has to be 3. In addition, number of ones in any given nonzero row/column has to be 3 with one zero row/column. Therefore, a general cell-to-TAL assignment with L lists and N MMEs or lists can be expressed bȳ The assignment of the maximum number of cells within a list has to satisfy the following expression where MX denotes maximum number of TAs in a list. Also, according to the above-mentioned discussion and Example 1, the cell/TA assignment in every TAL takes the following form C l k,n = C l n,k , ∀l ∈ [1, L], and k, n, k = n ∈ [1, N ] (4) 2) TAL-to-UE: According to the fact that any TAL associated with MME could be assigned to more than one cell, the assigned cell has not to violate the MME load of any nonzero n column, which could be achieved by where σ l k denotes the percentage of usage of list l associated with cell k.

C. Total Signaling of TAU and paging
The signaling cost of TAU and paging is a fundamental consideration when placing an idle-user in a cellular network [11,13]. In fact, the total signaling overhead created by TAU is a result of a user movement from a cell in one TAL to a different TAL. Whereas, the total signaling overhead created by paging is due to a necessity to locate a user via the core network. Since the total signaling overhead cost is inversely related to the size of TAL, there is negative correlation between the total signaling created by TAU and paging, and the size of TAL. This is the case due to the fact that the greater the number of cells within a certain list the lower is the probability for a user to move between lists.

D. Problem Description and Cost Functions
This subsection presents the centralized MME pooling problem [5,6,14]. The objective is to minimize the signaling overhead, which can be accomplished by minimizing two objective functions, namely, the total signaling cost of TAU and paging, and the total inter-list handover. The total signaling cost of TAU of cell k for any n nonzero column is equivalent to with UE k being the total number of UEs served by cell k, U C being the TAU cost of UE moving from one list to another, H C being the MME cost relocation over handover process, and O l M defining whether list l belongs to MME M or not. The total cost of paging of cell k for any n nonzero column can be evaluated by where G a denotes the arrival rate of paging and G C refers to the paging cost of equipment of a particular user. Recall the constraints in Equations (2) C l k,n = C l n,k , ∀l ∈ [1, L], k, n, k = n ∈ [1, N ] where Equations (11), (12), and (13) refer to boundary constraints. Accordingly, the total signaling cost of TAU and paging of the ith element is defined by and the first objective function is given by Hence, the total signaling overhead introduced by TAU and paging could be minimized through the minimization of J 1 in Equation (15). On the other side, the total inter-list handover of the kth element for any n nonzero column is with the second objective function being defined by such that the minimization of total signaling overhead could be achieved by diminishing the total inter-MME reallocation cost, which is defined by J 2 in Equation (17). The fluid flow model is popular and frequently employed to resemble the mobility behavior of a user in a system [5]. The fluid flow model presents the UEs traffic outflow rate of an enclosed area. In this specific case enclosed area refers to a single cell. For a given cell k, let P M denote its perimeter, ρ k refer to UE density of the cell, and v be the UE average velocity. Hence, the average of cell crossings with respect to time is defined by implying that the cells are hexagonal-shaped with length L H , which means that P M = 6L H . According to the above-mentioned discussion and the problem statement, the two objectives J 1 and J 2 are inversely related. A weighting scheme can be applied to the two conflicting objectives [24]. However, there is no guarantee that the obtained solution would be the best compromise solution [28].

E. Overview of model decomposition using MINLP
The problem formulation of the model mentioned above is NP-hard, thus, the problem can be divided into two sub-problems. One sub-problem may consider the cell-to-TAL/MME assignment which is allocated periodically to minimize the inter-list handover rate of UEs from one cell to another. This sub-problem mainly includes Equations (9), (10), (16) and (17). The second sub-problem could assume the set of TALs/MMEs given and, thereby, using the given data one can find the optimum usage ratio for each TAL/MME. The other sub-problem mainly includes Equations (11), (13), (6), (7), (14) and (15). For a full description of the previously mentioned two sub-problems using MINLP visit [6,7]. A scrutiny look at the final solution obtained by MINLP, one can find that the solution is confined on the Paretooptimal front. However, one objective is strictly favored the optimum solution over the other one, which is a significant weakness. Thus, the best compromise solution cannot be achieved using this method. Consequently, a random artificial search technique such as optimization algorithms based on the social behavior of animals could help in approaching a set of solutions on the Pareto-optimal of the NP-hard problem [36]. In that case, the two above-mentioned sub-problems would be solved simultaneously through the optimization algorithm and the best compromise solution would be defined.

IV. MULTI-OBJECTIVE PARTICLE SWARM OPTIMIZATION A. Particle Swarm Optimization
Particle Swarm Optimization (PSO) is considered to be an evolutionary heuristic technique. This technique imitates the cognitive and social behavior of animals in their natural environment such as fish schooling and bird swarming [20,25]. In the very beginning, PSO is initiated by releasing an initial random population of particles into the space. In this paper, the population size is fixed. Every particle within the space is a potential solution or, in other words, a potential candidate. It is worth mentioning that particle and candidate will be used interchangeably. Each single particle has a set of parameters such that the set of parameters swarms in the space irregularly in a multi-dimensional space. The swarming of particles is aimed at determining an optimal solution. The position of a particle represents a possible solution including the optimal solution. The velocity of every particle has a major impact of heading the best potential solution or alternatively, the best candidate position. Additionally, the subsequent position and velocity of each particle depend on: current velocity, current position, current position of the neighboring particles, and best position of the neighboring particles. Let P ∈ N and N ∈ N denote number of parameters to be optimized and population size, respectively. For any ith candidate solution within the space, define the candidate solution (position) at iteration t ∈ N by X i = [x i,1 (t) , . . . , x i,P (t)] ∈ R P and the velocity of every candidate by V i = [v i,1 (t) , . . . , v i,P (t)] ∈ R P for all t = 1, 2, . . . , T and i = 1, 2, . . . , N . Accordingly, the velocity update of the jth parameter within the ith swarming particle is given by [20,25] where α (t) denotes an inertial factor to be continuously reduced with every iteration such that α (t) = 0.99α (t − 1), rand 1 and rand 2 are real random numbers between 0 and 1. k 1 and k 2 refer to weighting factors associated with personal and social influence, respectively. x i,j denotes a best local particle within the previous generation while x i,j refers to the best global particle within all previous generations, for all i = 1, 2, . . . , N and j = 1, 2, . . . , P. It should be noted that the iterative solution starts at t = 1 continues until the total number of iterations set by the user is reached. The position update of a parameter within the swarming particle associated with Equation (19) can be expressed by [20,25] The candidate solution X i is initialized between the minimum and maximum boundaries of the space, termed x min and x max , respectively, such that x i,j (t) ∈ x j min , x j max . On the other side, the velocity with N intv being the number of intervals.

B. Pareto-optimal front and clustering
For a multi-objective optimization problem, or to be more specific, two-objective optimization problem, the trade-off between any two solutions can follow two possible scenarios: either one dominates the other or none dominates the other. For instance, if X 1 dominates X 2 according to the first objective, while X 2 dominates X 1 based on the second objective, both solutions are considered non-dominated and are located on the Pareto-optimal front. All the solution positioned on the Pareto-optimal front are combined into one set. Number of solutions on the Pareto-optimal front could be extremely large. Clustering, however, can reduce the number of non-dominated solutions within a large set while keeping the desired characteristics of the trade-off. In MOPSO, the local set of nondominated solutions is denoted by S with S i = X i (1). After which, non-dominated solutions are to be added to form a large set S i . Clustering is performed based on the distance between two pairs. The pairs separated by large distance are retained while minimal distance pairs are combined into one cluster [37]. Therefore, the local non-dominated set can be resized to satisfy the user settings, for instance S i , has the size S L . The non-dominated global set is denoted by S and includes all the non-dominated solutions from the first to the last generation. In the same spirit, clustering is performed to keep pairs within large distance and combine pairs within small distance. In that case, the non-dominated global set will not exceed a predefined size, for example S has the size S G . It should be remarked that X i and X k belong to the sets S i and S , respectively. The steps of the clustering algorithm can be found in [8,26].
Individual distances between all the members of the local set S i , and the members of the global set S are measured with respect to the objective space. Let X i and X i denote members of S i and S , respectively. For the case when X i and X i give minimum distances, they are to be selected as local best and global best of the ith particle, respectively.

C. Extraction of the best compromise solution
Once the Pareto-optimal front is available, one solution should be selected as the best compromise solution. The fuzzy approach is implemented to identify the best compromise solution. Define M as a number of non-dominated objectives on the Pareto-optimal front and N Obj as a number of objectives. Let the set of objective functions on the Paretooptimal front beJ c = J c 1 ,J  (22) for m = 1, 2, . . . , M , and c = 1, 2, . . . , N Obj . In our case there are two objectives, thus N Obj = 2. Next, the best compromise solution is obtained bȳ Given the solution ofμ m as in Equation (23) for m = 1, 2, . . . , M , the best compromise solution isμ m that has the maximum value. Figure 2 illustrates the flow chart of MOPSO computational algorithm.

D. Implementation of the MOPSO algorithm
For simplicity, the implementation process of the MOPSO algorithm can be subdivided into two stages: Stage I: Recall Cell/TA-to-TAL/MME Assignment in Section III, specifically Equations (2), (3) and (4). There are N ∈ N cells in the system withn ∈ [1, N ) cells being distributed within the list andn < N : such thatL ∈ Z LN ×N denotes all lists and the associated superscript refers to the list number. Thus, the number of the cell-to-TAL combination is given by L×N . It should be noted thatL l is a symmetric matrix for l = 1, 2, . . . , L, such that L l = L l ∈ Z N ×N , with being the transpose of the matrix. There aren nonzero rows withinL l and the rest of rows (N −n) are equal to zero, while the dimension of each row is 1 × N . Let L l l ∈ Z 1×N be a nonzero row, such that l is an integer between 1 and N ( l ∈ N and l ∈ [1, N]). In fact,n numbers within L l l are equal to 1 and the rest are zeros. Also,n nonzero columns within the square matrix L l are identical, for l = 1, 2, . . . , L. Accordingly, number of variables to be optimized and allocated within one list is N . In the L lists in Equation (24), there are L × N binary parameters to be optimized and subsequently allocated such that Stage II: Consider TAL-to-UE in Section III, Equation (5). Each TAL is associated with a specific MME and their relationship can be expressed as N k=1 σ l k · C l k,n = 1, l = 1, 2, . . . , L with σ l = σ l 1 σ l 2 · · · σ l N ∈ R 1×N for l = 1, 2, . . . , 10. Therefore, in the L lists in Equation (26), there are another L × N real parameters to be optimized. The decision variable σ l has to match the L l l , for instance In order to obtain σ l , define the following auxiliary variablē σ l = σ l 1σ l 2 · · ·σ l N ∈ R 1×N withσ l k ∈ [0, 1] for all l = 1, 2, . . . , 10 and k = 1, 2, . . . , N . Hence, we optimize for X 2 such that After optimizing X 2 , σ l can be defined as follows where C 1 l ,k is a cell associated with the nonzero row L l l . Next, to satisfy the equality in Equation (26), one has This concludes the implementation process of the MOPSO algorithm. Therefore, Stage I and Stage II summarize the optimization process and number of parameters to be optimized. It is worth re-mentioning that X 1 and X 2 in Stage I and Stage II, respectively, would be solved by MOPSO algorithm simultaneously. In total, there are 2LN parameters to be optimized. Let X i denote the vector of parameters to be optimized with dimension 2LN × 1 for all i = 1, 2, . . . , N . One can find the complete vector (X i ) to be The total number of parameters to be optimized in one particle is P = 2LN with X 1 ∈ Z N ×1 and X 2 ∈ R N ×1 . A set of random particles is initiated such that every single particle has P = 2LN parameters, and there are N particles within one iteration/generation. Therefore, the number of objective functions obtained from one generation is 2LN . Velocity and position of the ith particle are defied according to Equations (19) and (20), respectively. All the remaining steps are detailed above.

E. Computational flow using MOPSO
The computational flow of the proposed location management of LTE networks using MOPSO technique can be briefly summarized as follows: Step 1 (Initialization): initiate t = 1 and generate N random particles X i (1) ∈ R P ∀i = 1, 2, . . . , N where X i (1) = Evaluate the objective function of every particle. For every particle, generate a local best set S i (1) with X i (1) = X i (1) and S i (1) = X i (1) for all i = 1, 2, . . . , N . Search for non-dominated solutions and generate a non-dominated global set S (1) where the closest member of S (1) to X i (1) is marked as the global best X i (1) of the ith particle. Set the initial value of the inertia factor α (1).
Step 2 (iteration update): update the iteration number to t = t + 1.
Step 3 (Inertia factor update): update the inertia factor.
Step 4 (Velocity update): update the velocity V i (t) as defined in Equation (19) given the particle X i (t − 1), the local best X i (t − 1), and the global best X i (t − 1) for all i = 1, 2, . . . , N .
Step 5 (Position update): update the position X i (t) in accordance with Equation (20) given the particle X i (t − 1) and the velocity V i (t) for all i = 1, 2, . . . , N .
Step 6 (Non-dominated local set update): The updated position X i (t) is added to the set S i (t) for all i = 1, 2, . . . , N . Keep all non-dominated solutions and eliminate all dominated solutions within the local best set S i (t). If the size of S i (t) exceeds the predefined value S L , apply clustering and reduce the size to S L .
Step 7 (Non-dominated global set update): Set S (t) = S (t − 1). If any of the non-dominated solutions within S i (t) are non-dominated by S (t) and are not its members, add it to the set S (t). If the size of S (t) exceeds the predefined limit S G , apply clustering to reduce the size to S G .
Step 8 (Local best and global best update): The individual distances between members of S i and members of S (t) are evaluated with respect to the space of the objective function. The pair that gives minimum distance is selected as local best and global best of the ith particle, respectively.
Step 9 (Stopping criteria): If the maximum number of iterations defined by the user is reached, stop the program. Otherwise, go to Step 2.

V. PERFORMANCE EVALUATION AND SIMULATIONS
Optimization of the LTE networks location management using MOPSO is evaluated and compared to MINLP using MATLAB . In order to achieve reliable assessment, both algorithms, MOPSO and MINLP, have been implemented at different speed levels. The speed range has been subdivided into four levels: very slow speed stands for (0-8 m/s); slow speed indicated (8-16 m/s); normal speed should be interpreted as (16-25 m/s); and high speed implies speeds in the (25-33 m/s) range. The environment is assumed to have 30 cells with an average of 100 users are evenly distributed throughout the systems of cells. The users are allocated in a regular manner through the cells within the system. The system consists of 10 lists with 0.05 rate of paging. Clearly, the number of lists in the system equals the total number of MME. Table III outlines the full details of the model parameters associated with LTE networks.

A. MOPSO Set-up
The aim of the simulation section is to examine the convergence of MOPSO in comparison with MINLP. Recall Stage I subsection of the Cell/TA-to-TAL/MME Assignment (Section IV), more specifically, Equations (24), (3) and (4). As mentioned above, there are N = 30 cells in the system withn = 16 cells being distributed within the list. Thus, L l ∈ Z 30×30 denotes a complete list. Define L l l ∈ Z 30×1 to be a nonzero row, such that l is any integer between 1 and 30. L l l includes 16 cells with the value of 1 while the rest of the cells have the value of 0. For L = 10 lists in Equation (24), number of parameters to be optimized is 300 such that Each TAL is related to MME and the relation can be expressed by 30 k=1 σ l k · C l k,n = 1, l = 1, 2, . . . , 10 with σ l = σ l 1 σ l 2 · · · σ l 30 ∈ R 1×30 for l = 1, 2, . . . , 10. A new variableσ l is introduced such thatσ l = σ l 1σ l 2 · · ·σ l 30 ∈ R 1×30 . For 10 lists in Equation (26), one has 300 real parameters to be optimized, such that therefore, the decision σ l is defined by Therefore, there are 600 parameters to be optimized. Define vector X i as the total number of parameters to be optimized such that And the problem of optimization is defined as a minimization problem using MOPSO. The setting parameters of MOPSO are listed in Table IV with T = 400 being the total number of iterations. Also, it should be noted that the setting of parameters in Table IV were selected as a result of a set of trials. The inertia factor α (t) is calculated by α (t) = 0.99α (t − 1) and initiated at α (1) = 1.

B. MOPSO Total Signaling Cost
Using MOPSO, the total signaling of TAU and paging (J 1 ) and the total inter-list handover (J 2 ) are presented and compared to MINLP at four different speeds mentioned above.
The optimization algorithms, MOPSO and MINLP have been implemented and compared to assess: 1) The minimum value of the two objective functions J 1 and J 2 obtained by each algorithm. 2) The best compromise solution by MOPSO and the MINLP solution.
3) The totaling signaling overhead and average battery power consumption.
In order to evaluate the robustness of MOPSO algorithm in terms of its ability to approach the optimal solution, the algorithm has been implemented four times at each of the four speed levels considering different random initialization conditions of the X vector. The speed ranges examined are the following: [0,8], [8,16], [16,25] and [25,33] m/sec. For clarity, three out of the four trials have been presented. Additionally, the following color notation is used: red color demonstrates a best compromise solution, blue illustrates a solution on the Pareto-optimal front and magenta represents a candidate solution as well as search history obtained using MOPSO. Also, black refers to a solution obtained by MINLP. An illustration of the search history of a single trial of the speed [0,8] (m/sec) within the two-dimensional space is depicted in Figure 3 presenting all candidate solutions, Pareto-optimal front and the best compromise solution. Figure 4 demonstrate the Pareto-optimal front in blue color and the best compromise solutions in red color obtained by MOPSO plotted against MINLP in black color at speed ranges [8,16]. In spite of the fact that the solution generated by MINLP is located on the Pareto-optimal front as depicted in Figure 4, MINLP favored the objective function J 2 over J 1 . MOPSO, on the contrary, captured a set of solutions on the Pareto-optimal front with the trade-off represented by the best compromise solution.
The best compromise solution creates a balance between the inter-list handover which resides on UEs, probability of the user to travel between cells, and the cost of inter-MME, from one side, and the total signaling cost of TAU and paging which impact the EPC significantly by influencing the location management of UEs, on the other side. The significance of the proposed algorithm can be observed in terms of total signaling overhead and battery power consumption. The consumption estimation of the triggered TAU signal in a regular smartphone is considered 10 mW [5]. Figure 5 illustrates the total signaling overhead of MOPSO best compromise solution average values versus MINLP. MOPSO exhibited lower total signaling overhead than MINLP at various speed ranges as presented in Figure 5. As Figure 6 confirms that MOPSO has a significant advantage over MINLP with respect to battery power consumption. The average values of MOPSO best compromise solutions associated with battery power consumption versus MINLP are depicted in Figure 6. The overall quality of MOPSO convergence very close to the optimal solution is presented in terms of mean and standard deviation (STD). Table V lists the best compromise solutions of non-dominated objective functions recorded on the Paretooptimal front considering four trials at various speed ranges. Also, Table V outlines the mean and STD of the four trials at every speed range. The accuracy of convergence can be    evaluated by the percentage of STD with respect to the means, which is known as a relative standard deviation (RSD). RSD of J 1 is 0.7%, 0.7%, 1.3% and 1.26% and RSD of J 2 is 1.4%, 2.8%, 1.2% and 0.68% at the four speed ranges. Indeed, the RSD of J 1 and J 2 illustrates the remarkable performance of MOPSO when applied to the LTE networks location management. Likewise, the total power consumption has been evaluated in Table VI in terms of the mean and STD of the results obtained by MOPSO algorithm versus MINLP at different speeds. It can be noticed that MOPSO generated lower values of total power consumption than MINLP.

VI. CONCLUSION
Cellular networks are essential as they offer a variety of connectivity solutions. LTE and LTE-A provide reliable and fast wireless network service. Nonetheless, the total signaling overhead has been a concern, in particular with the rapid growth of more advanced cellular phones and phone applications. In fact, location management which allocate the idle user is directly related to the total signaling overhead. The total signaling overhead can be assessed by two elements, the total signaling cost of TAU and paging and the total inter-list handover. However, these two elements are adversely related. This study proposes a new formulation of the total signaling overhead as a true multi-objective problem where both conflicting objectives are optimized simultaneously. This novel approach is based on the multi-objective particle swarm optimization (MOPSO) which allows to minimize the total