Quantum Genetic Algorithm Based on Qutrits and Its Application

Typical approaches to designing quantum genetic algorithms are based on a concept of a qubit, a two-level quantum system. But many-valued quantum logic is more perspective from the point of view of the computational power. This paper proposes a quantum genetic algorithm based on a three-level quantum system in order to accelerate evolutionary process. Simulation using a set of standard test functions proves that the given algorithm is more effective and precise than the conventional quantum genetic algorithm.


Introduction
Evolutionary algorithms are one of the most widely used methods for solving the global optimization or search problem [1].In its classical implementation, the smallest unit of information is a bit, a structure which may be in one of the following states: 0 and 1.It is important that a classical bit does not allow any intermediate states.
Quantum genetic algorithm (QGA) is a new evolutionary algorithm based on a combination of quantum computation and conventional genetic algorithms technology [2,3].This algorithm can be applied to the same set of problems the conventional genetic algorithm is used for, but it allows significantly accelerating the evolutionary process through quantum parallelization and entanglement of the quantum state.Probabilistic mechanism of the quantum computations combined with the evolutionary algorithm provides a global search for a solution with a rapid convergence and a small value of population.These algorithms have demonstrated their effectiveness for solving combinatorial and functional optimization problems, mechanical engineering optimization problems, image processing, and many others [4][5][6][7][8].
QGA main concepts were first proposed by Narayanan and Moore [9].The basic unit of information used in a quantum computations is a qubit, a quantum system, which may be in the |0⟩ basis state or the |1⟩ basis state.Quantum nature of the qubit lies in the superposition principle, under which the qubit may be in any state which is a linear combination of basis states: with a normalization constraint: As a consequence of the superposition principle, qubit state space is incomparably greater than state space of a classical bit.Information, stored in  0 2 and  1 2 amplitudes, is the quantum part of the information.Practically,  0 2 g1  1 2 are numbers that specify the probability amplitudes of a qubit being in the states |0⟩ and |1⟩, respectively.The outcome of quantum measurement is a qubit in one of the basis states.However, it is important that the result of the measurement is not deterministic, like it is in the classical calculation, but probabilistic.
QGA uses a matrix representation for practical implementation: ) .

Mathematical Problems in Engineering
Qubit (1) state in this representation can be defined as A structured set containing  qubits represents a quantum chromosome.State vector of such chromosome is a superposition state of 2  basis states of the register | 1 ,  2 , . . .,   ⟩,  1 ,  2 , . . .,   = {0, 1}: All the information about a qubits system is represented by state vector |⟩.The only thing that can be done to such a system is transforming the initial state vector |⟩ into a new state |  ⟩ by applying a quantum gate.So quantum genetic algorithm is a state transition process from some starting state to an ending state using a quantum gate algorithm.The information about the problem solution is only contained by the vector of the system ending state |  ⟩.
The main disadvantage appearing during quantum state rotation is the need to use a table to look up the rotating angle, which limits the universality of the search.Fixed rotating angle has a negative impact on the search speed, which is the reason for QGA sometimes being implemented as an adaptive process for determining the rotating angle.To enhance the local searching ability and to get out of the local optimal solutions, GDA can be extended with conventional operators traditionally used in the classical genetic algorithm, such as quantum mutation operation and quantum disaster operation [10].

Quantum Genetic Algorithm
Based on Qutrits (QGA_3) Qudit is a structure with several, more than two, states, which can also be used to encode a chromosome.A qudit is basically a unit of quantum information, which may be in any of  states or in any superposition of those.Ternary quantum logic is the simplest type of the manyvalued logic.The basic unit of memory is called a qutrit.It has three basis states, |0⟩, |1⟩, and |2⟩.The state of the qutrit can be represented as a superposition in the form of a linear combination: with a normalization constraint: Matrix representation: For the practical implementation, qutrit (6) state can be represented as A system containing  qutrits has 3  basis states (as opposed to 2  states for binary logic).When switching from binary to many-valued logic, we get an exponential increase in the number of basis states, which leads to an increase in the algorithm performance for the same search accuracy.
Quantum chromosome length is determined by the desired search accuracy , search area [ min ,  max ], and the number of quantum system basis states: where  is a number of quantum system basis states.With  = 2, search accuracy  = 10 −6 , and search area [−1, 1], the length of the quantum chromosome has to be at least 21 qubits, while with  = 3 only 14 qutrits are required.

Qutrit Encoding.
Matrix representation of a quantum chromosome is a structure which consists of  qutrits.
Here {  ,   ,   } determine the state of th qutrit, 1nd  qutrits determine one individual in the population.The initial state of the qutrits holds no information about the system state, so the easiest way to initialize the base population is by setting all probability state amplitudes to be equal to one another [2].This means that at the end of the initialization each qutrit is in the state 2.2.Observation of Genes.Classical information about the problem's solution is located in the ending state vector |⟩  and can be retrieved as a result of quantum observation.Qutrit in one of the basis states, which is obtained as a result of quantum observation, is a classical representation of a quantum chromosome.Based on the approach described in [11], the following pseudocode for quantum -qutrit chromosome state observation can be suggested (Algorithm 1).The product of its ← 1 ( 7) else (8)  ← 2 (9) end if (10)  , respectively.The algorithm leads to the quantum chromosome being transformed into its classical representation in the ternary numbering system.
3. Quantum Rotating Gates.All information about the problem and the algorithm for solving it is contained in the quantum gate, so its algorithm is the key issue of any quantum genetic algorithm construction.
The pseudocode for the genetic operator updating the qutrit state (quantum gate of the system evolution) as a result of probability amplitudes   ,   , and   modification is presented in Algorithm 2.
The algorithm's concept can be illustrated by the image given below.The vertical lines set the probability amplitudes , , and .If the classical representation of the population best individual has, for example, a 2 on the position , all the amplitudes, except for , are multiplied by the factor . Its value lies in a range between 0 and 1 and is determined empirically.As a result, only the probability amplitude , which corresponds to the population best individual on the previous evolution stage, will increase.At the same time, probability amplitudes  and  will decrease.
b ∈ {0 1 2} It is important, that the described quantum gate algorithm eliminates the necessity to use a lookup table, which exists in the traditional QGA.
Initialized population ( = 0) includes all the possible solutions with the same probability.This means that the QGA process starts with a random search.As the evolution progresses in time, the distribution undergoes changes (Figure 1).
The mechanisms for the local search with  > 100 are illustrated in Figure 2.They demonstrate the role of  parameter.Its fixed value has a bad impact on the convergence speed.So, it is desirable to use an adaptive quantum gate operator.The listed data shows that it can influence the individual  in three possible ways: (1) The fitness of the individual is improved, but it does not become the best; ( end for (7) end if (8) end if (9) end for Algorithm 3: Quantum disaster.
(2) The fitness of the individual is improved, and it becomes the best; (3) The fitness of the individual is worsened.
This means that the bigger the number of individuals located in the area of the best individual , the more effective the process of the local convergence.
To sum up, QGA starts with the global search and switches to the local search automatically due to the change in the structural characteristics of the population individuals' distribution.

Quantum Disaster Operation.
Quantum population evolution may cause its falling into local optimal solution.This question is particularly important when the search area is large (big number of parameters) or when the fitness function has a complex topology around the global optimum.
Taking into consideration the small number of the quantum population, which is usually not bigger than  = 10 individuals, getting out of the local optimal solution can be achieved by a quantum disaster operation [12].It can be summarized as initialization of some population individuals, except for the best, with an initial state.Such operation "virtually" enlarges the size of the population and extends the area for optimal value search.
In the simplest case, quantum disaster operation can be implemented as Algorithm 3.

Simulation Test
QGA 3 based on a qutrits system is implemented using C++ programming language.The tests have been performed on the Intel Celeron CPU G1840 2.80 GHz, 4.0 GB RAM, Windows 7 Professional environment.
The majority of the quantum algorithms use similar primary ideas and approaches, which are based on the qubit representation of the population.For example, some algorithms based on the quantum computation ideas are presented in [13].The studies of the quantum algorithms and the possibilities of their modification and improvement are carried out in [14,15].Other potential approaches also include the option of the transition to qutrits representation of the population and can be considered as directions for further research.Considering this, the only algorithms compared by effectiveness and performance were QGA and the suggested QGA 3.
The analysis and performance evaluation have been performed on a number of test functions [16]: The function has one global minimum solution: Coldstein Price Function The function has one global minimum solution: Schaffer Function The function has one global minimum solution: Monopole and Six-Peak Camelback Function The function has one global maximum solution: Dual-Pole and Six-Peak Camelback Function The function has two global minimum solutions: The function has one global maximum solution: Ackley Function The function has one global minimum solution: Rastrigin Function −5.12 ≤   ≤ 5.12,  = 1, 2. The function has one global minimum solution: Optimization and Results.For QGA and QGA 3 simulation testing the following typical parameters were used: population size  = 10, total generations of iteration  = 500, and precision  = 1 ⋅ 10 −6 .
The Choice of  Parameter.The value of parameter  determines the evolution of the quantum system and has a key role in the algorithm.It determines the balance between the global (quantum observation operation) and local (quantum gate) search for the found solution specification.
To choose the optimal  value we have investigated its influence on the behavior of the main parameters which characterize the quantum system evolution, such as the average fitness of the best individual and the standard deviation of function value over 1000 algorithm runs.Its characteristic behavior for the example of  8 function optimization is presented in Figures 3 and 4. The quantum disaster operation was not applied in this case.
The research shows that if  < 0.95, the determinative operation in the population evolution is the quantum observation.Global search and the best individual's approximation to the optimal value take place because of its probabilistic nature.At the same time, the quantum rotation angle is greater than 15 ∘ and is too big to perform a local search.
As  value increases, the rotation angle approaches its optimal value ∼ 4 ∘ .It corresponds to  ∼ 0.996 ÷ 0.998: with this parameter value the fitness function and standard deviation reach the minimum value (see Figures 3 and 4).
Quantum system behavior for  > 0.998 is critically dependent on its evolution time (see Figure 5).Indeed, with  = 250 system simply does not have enough time to relax to the ending state, which corresponds to the fitness function optimal value.Enlarging the value of evolution time parameter  = 1000 practically simply causes the  algorithm running time to increase, with no optimal value improvement.
The role of the quantum disaster operation can be illustrated by Figure 6.The results of the  8 function optimization are given for parameters count  = 6.Taking into consideration the  8 complex topology around the optimum, the size of the population  = 50.
Based on the given analysis it has been established that for QGA 3 without the quantum disaster operation  ≈ 0.998 is the optimal value.Adding the quantum disaster operation to the QGA 3 algorithm allows producing the best result with  ≈ 0.975.
The optimization results of the  1 - 8 test functions with a standard set of input parameters are listed in Table 1.In Table 1,  opt denotes the function value of optimum;  is the average value of the best individual fitness function; sd is standard deviation of fitness function;  is average elapsed time need for the algorithm to run.
A search for  7 and  8 functions global minimum has been carried out separately, using parameter count  = 6 and a standard set of other algorithm parameters (see Table 2).The results reveal that QGA 3 is more, but not enough, effective, especially in regard to function  8 .
The result can be significantly improved by enlarging the population size at least to  = 50 and adding a quantum  disaster operation, with the other parameters remaining unchanged.The results of experiments after the described improvements to the algorithm are listed in Table 3.

Conclusions
In this paper a new quantum genetic algorithm QGA 3 has been proposed.The algorithm is based on the idea of using qutrits for the representation of the population individuals.A simple algorithm, which does not require a lookup table, has been used for constructing the evolution operator which is responsible for the quantum system state updates.
The experiments have been performed using a typical set of test functions.On the basis of the experimental results we have concluded that QGA 3 allows decreasing the computational time while being more productive and having better convergence, compared to the conventional QGA.
We have also obtained that, in a case with a lot of parameters to be optimized, the QGA 3 efficacy can be improved by adding a quantum disaster operation.In the same time, the increase in algorithm operating time does not exceed 4%.
Further promising investigations include the research of -valued quantum logic units, qudits and their possible applications for representation of the population individuals.This is expected to be a subject of separate studies.

Figure 1 :Figure 2 :
Figure 1: Average distribution of the population individuals  over the search area on different stages of the evolution using an example of  4 function optimization for 100 runs.

Figure 3 :
Figure 3: Average population fitness as a function of  parameter.Evolution time  = 500.

Figure 4 :
Figure 4: Standard deviation of the population fitness as a function of  parameter.Evolution time  = 500.

Figure 5 :
Figure 5: Average fitness of the best individual as a function of  parameter for different evolution time parameter values.

Figure 6 :
Figure 6: Quantum disaster operation influence on the average population fitness as a function of  parameter:  -QGA 3 without the disaster operation; × -QGA 3 with the disaster operation.

Table 3 :
The influence of the quantum disaster operation on QGA 3 efficacy.