Multiparameter Adaptive Optimisation of MSE Osseous Expansion Position

The accuracy of the implantation position of the MSE osseous expansion anchorage implant is a key issue in the treatment of osseous expansion, which also suﬀers from the drawback of falling into prematureness too early when the standard QPSO algorithm is used for its multiobjective optimisation. In this study, an adaptive improved QPSO algorithm is proposed to address the above problems. Firstly, a partitioned retrieval strategy is used to divide the population into an auxiliary class group and a main iterative population, and the respective search iteration intervals of the populations are assigned, thus optimising the initialisation mechanism of the standard algorithm, and then, the pheromone mechanism in the ACO algorithm is introduced to make the particles in the QPSO algorithm carry pheromones, and the particles determine their direction of travel by sensing the pheromone concentration in each path, thus improving the search ability of the particles. Experimental simulation results show that the improved strategy proposed in this study eﬀectively improves the population diversity of the standard QPSO algorithm, avoids the algorithm from entering local optimum, and has good application in MSE osseous expansion position optimisation.


Introduction
e MSE (maxillary skeletal expander), also known as microimplant brace-assisted rapid maxillary expansion, is characterised by its reliance on four BMK anchorage implants in the maxilla, which act directly on the midpalatal suture with the aid of the anchorage implants to open the midpalatal suture and widen the upper arch by drawing the bone, effectively avoiding the tilt compensation associated with traditional expansion devices with the lack of sufficient skeletal effect [1][2][3]. However, the implant position of the anchorage implants is currently determined by the surgeon's own experience, so the rapid and accurate acquisition of the optimal implant position has become an urgent problem.
Studies on the optimisation of MSE osseous expansion position have been less reported in the domestic and international literature, mostly on the analysis of the effect of MSE osseous expansion. Paredes et al. [4] evaluated the differences in bone, alveolar bone, and tooth composition induced by the osseous expander in the middle of the microimplant-supported surface, and Moon et al. [5] studied the molar inclination and the surrounding alveolar bone changes after MSE osseous expansion by CBCT. Abedini et al. [6] analysed and studied the changes in the human facial soft tissues after MSE osseous expansion. Hartono et al. [7] performed a finite element analysis of the palatal alveolar bone, palatal suture, zygomatic suture, microbone, and the surrounding bone of the first molar after MSE osseous expansion implantation. Cantarella et al. [8] improved the implant position of MSE based on biomechanics. Giudice et al. [9] optimised the implantation parameters of four BMK anchorage implants to improve their stability.
Optimisation of the MSE osseous expansion position is essentially a multiobjective optimisation process, and the QPSO (quantum-behaved particle swarm optimisation) algorithm is a common multiobjective optimisation algorithm currently available. QPSO algorithm has effective applications in various fields, such as parameter selection in industrial manufacturing [10,11], image segmentation [12], and prediction problems in systems such as SLAM [13,14].
Kumar et al. [15] proposed an improved QPSO algorithm based on constrained optimisation to improve the multiobjective optimisation performance. Guo et al. [16] used a multilevel perturbation strategy to optimise the QPSO algorithm to improve the convergence performance and search range. Zhang et al. [17] proposed a quantum-behaved particle swarm algorithm based on generalised space transform search, which solved the problem of premature convergence of the standard algorithm. Grotti et al. [18] proposed an improved QPSO algorithm based on multiobjective archiving, which optimised the multiobjective optimisation effect. Dong et al. [19] used BP neural networks to improve the QPSO algorithm, solving the problems of slow convergence and easy oscillation. Chen et al. [20] proposed a hybrid quantum algorithm based on a local optimisation strategy and an improved optimised rotation angle. Although the abovementioned improved algorithms have better multiobjective optimisation effects, they have certain limitations in the application of MSE osseous expansion position optimisation.
In this study, an adaptive improved QPSO algorithm is proposed to address the phenomenon of early falling into prematureness exhibited by the standard QPSO algorithm in the application of MSE osseous expansion position optimisation. e initialisation mechanism is optimised using a partitioned retrieval strategy, and the search capability of the standard algorithm is improved using a pheromone mechanism, thus reducing the probability of early falling into prematureness of the standard QPSO algorithm.

e Parametric Mapping Model for the MSE Osseous
Expansion Position. In order to obtain the optimal position of the MSE osseous expansion, the width of the midpalatal suture at the cusp (SOC) x 1 , the maxillary basal bone width (WMB) x 2 , the width of the midpalatal suture at the molar (SOM) x 3 , the nasal width (NW) x 4 , the maxillary alveolar bone width (MAW) x 5 , the intermolar width of palatal root (IWR) x 6 , the molar buccal cusp width (IWC) x 7 , the inclination of the right molar x 8 , the inclination of the left molar x 9 , and the molar palatal cusp angle x 10 were first selected as variable factors to construct a parametric mapping model of the MSE osseous expansion position, as shown in Figure 1. e methods for obtaining the variable factors are shown in Table 1.
On the basis of the above variable factors, a dataset of the optimal MSE osseous expansion positions was collected for each case, containing the three-dimensional spatial positions A, B, C, and D of the four implants, as shown in Figure 2.
Let the set of parameters for the optimal MSE osseous expansion position be f(x); then, the parametric mapping model for the MSE osseous expansion position can be expressed as where a i (i � 1, 2, . . . , 10) is the parameter variation factor.

QPSO-Based Parametric Multiobjective Optimisation.
In this study, QPSO is used for multiobjective optimisation of the parametric mapping model based on quantum behaviour. Unlike the PSO algorithm, the QPSO algorithm no longer depends on the position and velocity of the individuals for its current population state during the update and iteration process. In quantum space, because of the particle and fluctuating nature of individual particles, the quantum behaviour of particles cannot be described by definite values of particle position and particle velocity as in the PSO algorithm, so instead of particle position and velocity, wave functions are used. Since the intensity of the wave function at a point in quantum space is proportional to the probability density of the particle's appearance, it is shown in equation (2) as where X is the position vector of the particle in the quantum space, Q is the probability density function, and ψ(X, t) is the wave function, while the probability density function should satisfy the normalization condition: e control of particle motion is achieved by solving Schrödinger's equation to obtain the probability density function and probability distribution function for the appearance of a particle at a point in space, as shown in equations (4) and (5):  Mathematical Problems in Engineering where p is the particle attractor, L is the potential well δ length, and X denotes the particle position. In order to solve the change from quantum search space to classical space, Monte Carlo simulation is used to solve the equation for the position update of the particle in classical space. e particle is moved inside a potential well centred on p, and its position is determined by where u is a random number uniformly distributed on the interval (0, 1). In a D-dimensional search space, let the particle population consist of N particles with a potential problem solution; then, the position of the i particle at the t iteration can be expressed as where the range of i is 1 ∼ N.
In the QPSO algorithm, the individual best position and the population best position are calculated as When the particle i is in t + 1 iteration, then the position of the particle can be expressed according to equation (6) as where L i d (t) is then expressed by where α is the constriction-broadening factor and C d (t) is the average optimal position, which is the centroid of all particles' own optimal positions and can be determined from us, the equation for the evolution of the position update of the particle is obtained as . (12) Combined with the parametric mapping model of the MSE osseous expansion position, the set of parameters for its optimal MSE osseous expansion position is f(x). e objective of the multiobjective optimisation is to find the optimal parameter change factor so that the deviation of the Width across the buccal cusp of the right and left first premolar x 3 Distance between a point on each side of the midpalatal suture corresponding to the bilateral maxillary first molars x 4 Maximum width of the nasal part at the first maxillary molar level on both sides x 5 Distance between the lowest points of the alveolar bone on the buccal side of the first maxillary molar level, bilaterally x 6 Distance between the points of the palatal roots on each side of the first maxillary molar level, bilaterally x 7 Distance between the points of the buccal cusps on both sides at the first maxillary molar level bilaterally x 8 Medial angle between the line of the buccal cusp of the right molar and the cusp of the palatal root and the palatal plane x 9 Medial angle between the line of the buccal cusp of the left molar and the cusp of the palatal root and the palatal plane Mathematical Problems in Engineering 3 MSE osseous expansion position is minimised; therefore, the specific process of its multiobjective optimisation is as follows: Step 1: input the parameter variation factor into the algorithm and initialize the population, randomly generating N particles Step 2: calculate the average optimal position of the population particles Step 3: calculate the current individual fitness function value of each particle, and update it if it is better than or less than the individual fitness function value of the particle in the last iteration Step 4: calculate the global optimal position in the population Step 5: for each dimension of the particles in the population individuals, calculate the position of a random point Step 6: perform the evolution of the position update according to equation (12) Step 7: if it is satisfied, then end the update and iteration to obtain the optimal parameter change factor

Algorithm Simulation Test.
Simulation experiments were conducted on the QPSO-based parametric multiobjective optimisation strategy, and the optimal position error obtained by the algorithm was set as the evaluation value, and three sets of experimental data were used for the optimisation search test, as shown in Table 2. e test results of the QPSO algorithm are shown in Figure 3.
Based on the above experimental results, it can be seen that the QPSO algorithm shows excellent search performance in the early stage of the algorithm, but falls into local optimum in the middle stage.

Initialisation Mechanism Based on Partitioned Search.
Without degrading the convergence performance of the algorithm, a new initialisation strategy is proposed in this study in order to improve the quality of each particle and the diversity at the late stage of the population. Based on two types of populations, the auxiliary population and the main iterative population, the auxiliary population is defined as € X i , satisfying 1 ≤ i ≤ m, and m is the number of auxiliary populations. Meanwhile, the main population is defined as _ X. In order to create the main iterative population _ X, each auxiliary population € X i is given a random initialisation. As a result, such random initialisation throughout the search space leads to the generation of diverse populations in uncertain regions with high quality for each particle in the population. e auxiliary populations are updated and moved according to the basic method of the QPSO algorithm, and then, the globally optimal individuals of each auxiliary population € X i are added to the main population _ X. e basic steps for creating the main population _ X are shown in Figure 4. e search space is then partitioned into a number of zones and the auxiliary populations are initialised in the search intervals of these partitions. e partitioning strategy follows a better coverage of the optimised area, resulting in a good view of each auxiliary search interval. e whole partition [min , max] is divided into several iterative search regions with a total number of max β auxiliary intervals, and a start node of the partitioning process is defined as the centre of the search region.
To define the boundaries of each partitioned region, calculate the step size and centre of each partitioned region according to equations (13) and (14): e total number of each partitioned region is then solved according to equations (15) and (16): where zone β is the current number of populations, i.e., zone β � 1, 2, . . . , max β . e partitioned search space is depicted in Figure 5.
Different partitioned intervals have different sizes and the size of the auxiliary population € X i is defined as € X is , which is based on the size of the corresponding search interval. e size of the partition interval follows the principle that the smaller the search space, the fewer individuals in the search space, and conversely, the larger the search interval, the more individuals in the search space. en, the size of the partition interval is shown as where zone β is the number of particles in each auxiliary partition and H β is a fixed value, which in this study is taken as 2.

Pheromone Mechanism-Based Search Capability
Optimisation. In order to further improve the search capability of the QPSO algorithm, the pheromone mechanism in the ACO algorithm is introduced in this study, so that the particles in the QPSO algorithm carry pheromones. e state of a particle can be represented by the wave function ψ(Y):

4
Mathematical Problems in Engineering where each particle can release pheromones,;τ ij (t) indicates the content of pheromones on the path (i, j) of the particle at the moment t and sets to start with the same content of pheromones on the path. e particle decides the direction of its travel by sensing the pheromone concentration on each path, and the next position shift is influenced by both the  amount of heuristic information η ij (t) and residual information τ ij (t) on each path. θ denotes the relative importance of the residual pheromone and c denotes the relative importance of the expected value.
In this study, we assume that an update of the pheromone on the path is performed after all particles of the population have cycled once, and the update formula is as follows: where ρ is the residual coefficient of the pheromone on the path. To avoid infinite accumulation of information, ρ takes values in the range [0, 1), while 1 − ρ denotes the decay coefficient of the information. Δτ ij (t) is the pheromone increment on the path (i, j) and the initial moment is Δτ ij (t) � 0. Δτ k ij (t) is the amount of information that the particle k at the moment t left in the path (i, j).
When the particle k passes through the path (i, j) in this cycle, then we have in the other cases, Δτ k ij � 0, where Q is a constant and L k represents the total length of the path travelled by the kth particle in this cycle.
In this study, we use the average of all the particles' individual optima mbest, which is the average best position, i.e., where the characteristic length L is where α is the constriction-broadening factor, which usually works better when decreasing linearly from 1.0 to 0.5, so we have where t is the current number of evolutionary generations and T is the number of terminating evolutionary generations. e evolution equation for its particles becomes where p i,j (t) is the centre of the potential well of the ith particle in the jth dimension:

Performance Simulations of the Improved QPSO
Algorithm. In order to verify the performance of the improved QPSO algorithm, simulation tests were carried out. e optimal position error obtained by the algorithm was set as the evaluation value, and three sets of experimental data were used for the optimisation test, as shown in Table 3.
e results of the optimisation search tests of the improved QPSO algorithm are shown in Figure 6.
From the above simulation results, it can be seen that the improved QPSO algorithm proposed in this study, due to the optimisation of the initialisation mechanism and the search capability, improves the population diversity and avoids the algorithm from entering into local optimum.

Application of the MSE Osseous Expansion Position.
In order to test the application performance of the improved QPSO algorithm proposed in this study, the detection parameters of the cases were optimised to obtain the best MSE osseous expansion position, and the detection parameters of the cases are shown in Table 4.    erefore, the error of the improved QPSO algorithm in finding the optimal MSE osseous expansion positions for this case is 0.0975, which is much smaller than that of the standard QPSO algorithm of 0.3959, reflecting the superiority of the improved QPSO algorithm proposed in this study.

Summary
e optimal MSE osseous expansion arch position is obtained by the computer intelligence algorithm, and a mathematical model is established to produce implant guides using 3D printing technology, which helps to improve the success rate of MSE treatment and reduce its risk. In this study, an adaptive improved QPSO algorithm is proposed to address the problem that the multiobjective optimisation process of the standard QPSO algorithm for the parametric mapping model of MSE osseous expansion position is prone to local optimum, and the experimental simulation results show that the improved QPSO algorithm proposed in this study has better multiobjective optimisation performance and is well applied in the optimisation of MSE osseous expansion position.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.