Enhancing quantum annealing performance by a degenerate two-level system

Quantum annealing is an innovative idea and method for avoiding the increase of the calculation cost of the combinatorial optimization problem. Since the combinatorial optimization problems are ubiquitous, quantum annealing machine with high efficiency and scalability will give an immeasurable impact on many fields. However, the conventional quantum annealing machine may not have a high success probability for finding the solution because the energy gap closes exponentially as a function of the system size. To propose an idea for finding high success probability is one of the most important issues. Here we show that a degenerate two-level system provides the higher success probability than the conventional spin-1/2 model in a weak longitudinal magnetic field region. The physics behind this is that the quantum annealing in this model can be reduced into that in the spin-1/2 model, where the effective longitudinal magnetic field may open the energy gap, which suppresses the Landau--Zener tunneling providing leakage of the ground state. We also present the success probability of the $\Lambda$-type system, which may show the higher success probability than the conventional spin-1/2 model.


I. INTRODUCTION
Quantum annealing is an interesting approach for finding the optimal solution of combinatorial optimization problems by using the quantum effect [1][2][3][4]. The combinatorial optimization problems are ubiquitous in the real social world, therefore the spread of quantum annealing machine with high efficiency and high scalability will give impacts and benefits on many fields, such as an industry including drug design [5], financial portfolio problem [6], and traffic flow optimization [7]. After the commercialization of superconducting quantum annealing machine by D-Wave Systems inc. [8], hardware has been investigated and developed [9][10][11][12][13].
However, there are bottlenecks for implementing scalable quantum annealing machine; for the conventional and scalable quantum annealing machine may not have a high success probability for finding the solution of a combinatorial optimization problem because of the emergence of the first order phase transition, where the energy gap between ground state and the first exited state closes exponentially as a function of the system size [2]. In this case, it necessitates an exponentially long annealing time for finding the solution of the problem [14][15][16]. In the case of the second oder phase transition, on the other hand, an annealing time for finding the solution may scales polynomially as a function of the system size [17].
To propose an idea for finding high success probability is one of the most important and challenging issue in the field of quantum annealing. One of the approaches for obtaining the high success probability is to engineer the scheduling function for the driving Hamiltonian and the problem Hamiltonian, such as a monotonically increasing scheduling function satisfying the local adiabatic condition [18], the reverse quantum annealing [19] implemented in D-wave 2000Q [20], inhomogeneous sweeping out of local transverse magnetic fields [21,22], and a diabatic pulse application [23]. Another is to add an artificial additional Hamiltonian for suppressing the emergence of the excitations with avoiding the slowing down of annealing time, which is called shortcuts to adiabaticity by the counterdiabatic driving [24][25][26][27], and to add an additional Hamiltonian for avoiding the first order phase transition [17,28,29]. In this paper, we study the possibility of other approach: to employ a variant spin, such as a qudit, in the quantum annealing architecture.
Recently, two of the authors have studied the quantum phase transition in a degenerate two-level spin system, called the quantum Wajnflasz-Pick model, where an internal spin state is coupled to all the same energy internal states with a single coupling strength, and to all the different energy internal states with the other single coupling strength [30]. In the earlier study, this model is found to show a several kinds of phase transition while annealing; single or double first-order phase transitions as well as a single second-order phase transition, depending on an internal state coupling parameter [30], which suggests that the quantum annealing of this model may be controlled by an internal state tuning parameter. However, the study is based on the static statistical approach using the mean-field theory, because only the order of the phase transition has been interested in. Therefore, the enhancement of the success probability for quantum annealing based on degenerate two-level systems is not clear yet. Furthermore, they employed a fully-connected uniform interacting system, and it is unclear whether their idea works that a double (or even-number of) first-order phase transition while annealing would bring the system back into the ground state at the end of the annealing, where the even number of the Landau-Zener tunneling may happen with respect to the ground state.
In the present paper, we clarify the success probability of the quantum annealing in the quantum Wajnflasz-Pick model, focusing on (i) the Schrödinger dynamics, (ii) arXiv:1906.11393v4 [cond-mat.stat-mech] 23 Dec 2019 eigenenergies, and (iii) non-uniform effects of the spininteraction as well as the longitudinal magnetic field. We find that the quantum Wajnflasz-Pick model is more efficient than the conventional spin-1/2 model in the weak longitudinal magnetic field region as well as in the strong coupling region between degenerate states. We also find that the quantum Wajnflasz-Pick model is reducible into a spin-1/2 model, where effect of the transverse magnetic field in the original Hamiltonian emerges in the reduced Hamiltonian not only as the effective transverse magnetic field but also as the effective longitudinal magnetic field. As a result, this model may provide the higher success probability in the case where the effective longitudinal magnetic field opens the energy gap between the ground state and the first excited state. We also evaluate the success probability in another variant spin, a Λ-type system [31][32][33][34][35][36][37][38][39][40], which has three internal levels. This model also shows the higher success probability than the conventional spin-1/2 model in the weak magnetic field region.
A multilevel system is ubiquitous, which can be seen, for example, in degenerate two-level systems in atoms [41,42], Λ-type atoms [31,32,34], Λ-, V -, Θ-and ∆-type systems in the superconducting circuits [33,[35][36][37][38][39][40]43] as well as Λ-type systems in the nitrogen-vacancy centre in diamond [44]. We hope that insights of our results in the degenerate two-level system and knowledge of their reduced Hamiltonian inspire and promote further study as well as future engineering of quantum annealing.

II. QUANTUM WAJNFLASZ-PICK MODEL
A conventional quantum annealing consists of the spin-1/2 model, where the time dependent Hamiltonian is given by [ whereĤ z,x are a problem and driver Hamiltonian, respectively, and s ≡ t/T is the time t ∈ [0, T ] scaled by the annealing time T . The problem HamiltonianĤ z with the number of spins N , which encodes the desired optimal solution, has a non-trivial ground state. In contrast, the driver HamiltonianĤ x has a trivial ground state, where the driver HamiltonianĤ x must not be commutable with the problem HamiltonianĤ z . A problem Hamiltonian and driver Hamiltonian are typically given bŷ whereσ x,z are the Pauli matrices, J ij is the coupling strength between spins, h z i is the local longitudinal magnetic field, and h x i is the local transverse magnetic field.
The time-dependent total HamiltonianĤ(s) gradually changes from the driver HamiltonianĤ x to the problem HamiltonianĤ z . If the Hamiltonian changes sufficiently slowly, the quantum adiabatic theorem guarantees that the initial quantum ground state follows the instantaneous ground state of the total Hamiltonian [45]. We can thus finally obtain a non-trivial ground state of the problem Hamiltonian starting from the trivial ground state of the deriver Hamiltonian making use of the Schrödinger dynamics.
The quantum Wajnflasz-Pick model is a quantum version of the Wajnflasz-Pick model [46], which can describe one of the interacting degenerate two-level systems. In the language of the quantum annealing, the problem Hamiltonian and the driver Hamiltonian are respectively given by [30] (Schematic picture of this model is shown in Fig. 1.) The Hamiltonian of this model can be simply obtained by replacing the Pauli matricesσ x,z in Eqs. (2) and (3) with the spin matrices of the quantum Wajnflasz-Pick model τ x,z . The spin operatorτ z is given by [30] where g u(l) is the number of the degeneracy of the upper (lower) states. The spin-operatorτ x in the driver Hamiltonian is given bŷ where A(l) is a (l × l) matrix with the off-diagonal term ω, given by Here, ω is a parameter of the internal transition between the degenerated upper/lower states. The matrix 1(m, n) is the (m × n) matrix, where all the elements is unity, which gives the transition between the upper and lower states. The constant c is the normalization factor, where the maximum eigenvalue is normalized to be +1, so as to be equal to the maximum eigenvalue ofτ z . In the following, for the consistency to the earlier work [30], we consider a uniform transverse field h x i ≡ 1, and also take the parameter of the internal transition to be real ω = ω * with ω > −1. In the case where (g u , g l ) = (2, 1), we havê with c = (ω + √ 8 + ω 2 )/2, which is a kind of the ∆-type system [38].
In this paper, we employ the common sets of parameters in both quantum Wajnflasz-Pick model and the conventional spin-1/2 model, including the coupling strength J ij , the magnetic fields h z,x i , and the annealing time T . By using these parameters, we can obtain the same spin configuration (+1 or −1) in the ground state of the problem Hamiltonian. We thus compare efficiency of these models from the success probability.

III. SCHRÖDINGER DYNAMICS
In order to numerically calculate the success probability of the quantum annealing, we employ the Crank-Nicholson method [47] for solving the Schrödinger equation In this method, the time-evolution of the wavefunction is calculated by using the Cayley's form [47] |Ψ(t + ∆t) = 1 − iĤ∆t/2 Although the inverse matrix is needed, this method conserves the norm of the wave-function and is second-order accurate in time [47]. We first consider the fully connected model, where the spin-spin coupling is ferromagnetic and the longitudinal magnetic field is uniform h z i ≡ h, which is consistent with the earlier work [30]. For example, in the Compare the success probability of the quantum Wajnflasz-Pick model, P ≡ | Ψ(T )|Ψ 0 (T ) | 2 , with that of the conventional spin-1/2 model denoted as P 1/2 , where |Ψ(T ) is the final state obtained from the timedependent Schrödinger equation. In almost all regions in the ω-h plane, efficiencies of both models are almost the same, where the ratio of the success probability of the quantum Wajnflasz-Pick model and that of the conventional spin-1/2 model is almost unity (Fig. 3). On the other hand, in the regime of the weak longitudinal magnetic field h, we can find higher or lower efficiency regions in the quantum Wajnflasz-Pick model, compared with the spin-1/2 model. In the spin glass model, a non-trivial state may emerge in the weak longitudinal magnetic field limit [48]. In a p-spin model where p = 3, 5, 7, · · · , the energy gap is known to close exponentially and the firstorder phase transition emerges in the absence of the lon- gitudinal magnetic field [15]. In this sense, it is of interest that the quantum Wajnflasz-Pick model may provide the high efficiency in the weak longitudinal magnetic field region.
In the case where (g u , g l ) = (2, 2), where the numbers of upper and lower states are equal, the success probability of the quantum Wajnflasz-Pick model is almost equal to that of the conventional spin-1/2 model (Panel (a) in Fig. 4). In the case where (g u , g l ) = (3, 2), the success probability of the quantum Wajnflasz-Pick model is almost equal to that of the case where (g u , g l ) = (2, 1), where the differences between the number of the upper states and that of lower states are the same in both cases ( Fig. 3 and Panel (b) in Fig. 4).

IV. EIGENVALUES
Eigenvalue spectrum of the instantaneous Hamiltonian may help to understand these higher or lower success probabilities of the quantum Wajnflasz-Pick model than that of spin-1/2 model, although eigenvalues of the instantaneous Hamiltonian shows tangled spaghetti structures (Fig. 5). For example, in the case where (ω, h) = (0.8, −0.1), the energy gap between the ground state and the first excited state clearly closes once, which causes the low success probability (Panel (c) in Fig. 5). In the case where (ω, h) = (0.8, 0.02), the ground state and the first excited state are finally merged at the annealing time, where the degeneracy would cause the high success probability (Panel (a) in Fig. 5). However, according to the following discussion, it will be found that the latter explanation would not be correct in the case where (ω, h) = (0.8, 0.02). From panels (b) and (d) in Fig. 5, many crossings of eigenvalues are found to emerge. It suggests that there are no matrix elements in some states, and we may find symmetry behind the present quantum Wajnflasz-Pick model, where the Hamiltonian would be block diagonalized by a unitary operatorÛ . Since the energy spectrum of the original quantum Wajnflasz-Pick model shows very complicated behavior, it would be better to find out the reason of the high/low success probability from the reduced Hamiltonian, which are truly relevant for the efficiency of the quantum annealing. For example, in the case where (g u , g l ) = (2, 1), the single-spin Hamiltonian in the quantum Wajnflasz-Pick model is decomposable, where the irreducible representation is given bŷ for arbitrary values of s, by using the unitary operator where As a result, we may reduce a quantum annealing problem in the single-spin quantum Wajnflasz-Pick model into that of the spin-1/2 model, the Hamiltonian of which is given in the form Since the initial ground state of the single-spin Hamiltonian is given by |Ψ(s = 0) ∝ (c/2, c/2, 1) T in the original quantum Wajnflasz-Pick model, this state can be mapped toÛ |Ψ(s = 0) ∝ (c/ √ 2, 0, 1) T . It indicates that the initial ground stateÛ |Ψ(s = 0) can be also projected to the Hilbert space of the reduced Hamiltonian H(s).
This reduction of the single-spin problem in the case where (g u , g l ) = (2, 1) can be generalized to an interacting N -spin problem (Fig. 6). A quantum annealing problem of the original quantum Wajnflasz-Pick model is reduced into that of the spin-1/2 model, given in the formĤ where with As in the single-spin case, the initial ground state of the original N -spin quantum Wajnflasz-Pick model can be also projected to the Hilbert space of the reduced Hamiltonian (15). The coupling J ij in the reduced Hamiltonian is the same as that of the original Wajnflasz-Pick model. The effective longitudinal magnetic field h z eff,i in the reduced Hamiltonian also reaches the same value as that of the original Wajnflasz-Pick model at the end of the annealing: h z eff,i (s = 1) = h z i . Eigenvalues of the reduced spin-1/2 model exactly trace eigenvalues in the original Wajnflasz-Pick model (Fig. 5). The time-dependence of the ground state population of the problem Hamiltonian is confirmed to show the completely same behavior between the reduced model and the original model. This effective model clearly explains behavior of success probability of the quantum Wajnflasz-Pick model shown in Fig. 3. Note that the coefficient c is a positive real number such that the maximum eigenvalue of τ x is unity, and we take h x i = 1. Then, h i (s) ≥ 0 always holds during the annealing time 0 ≤ s ≤ 1. In the case where the longitudinal magnetic field h z i is very large, |h z i | |ω|h i (0), the effect of the original longitudinal magnetic field h z i is dominant compared with the effective additional term ωh i (s) except at the very early stage of the annealing s |ωh i (0)/h z i |. In this case, the problem Hamiltonian in the reduced model is almost the same as that in the conventional spin-1/2 model in Eq. (2). As a result, the success probability of the quantum Wajnflasz-Pick model is almost the same as that of the conventional spin-1/2 model, which provides P P 1/2 .
In the case where the original longitudinal magnetic field h z i is not large, the effective additional longitudinal magnetic field ωh i (s) cannot be neglected compared with h z i . When the effective additional field is in the same direction as the original longitudinal field, the total effective longitudinal magnetic field h z eff,i (s) is enhanced, which opens the energy gap between the ground state and the first excited state (Panels (a) and (b) in Fig. 7). This region is given by the condition ωh z i > 0, which is consistent with the result shown in Fig. 3. As a result, the success probability of the quantum Wajnflasz-Pick model become superior to that of the conventional spin-1/2 model. When the effective additional field is in the opposite direction to the original longitudinal field, the total effective longitudinal magnetic field h z eff,i (s) is diminished, which closes the energy gap between the ground state and the first excited state (Panels (c) and (d) in Fig. 7). This region is given by the condition ωh z i < 0, which is consistent with the result shown in Fig. 3. As a result, the success probability of the quantum Wajnflasz-Pick model become inferior to that of the conventional spin-1/2 model.
Behavior of success probability is also explained by the reference of the annealing time [49] T where Here, |Ψ 0(1) (s) and E 0(1) (s) are the wave functions and eigenenergies of the ground (first-excited) state with respect to the instantaneous Hamiltonian, respectively. Annealing machine needs the annealing time T much larger than T . Let T * ≡ b(s)/∆ 2 (s) be an instantaneous reference time of the annealing. The maximum value of this time T * in the reduced Wajnflasz-Pick model given in (15) is suppressed compared with that of the conventional spin-1/2 model, where the effective additional field ωh i (s) is in the same direction as the original longitudinal field h z i (Panels (a) and (b) in Fig. 8). It is consistent with the case where the quantum Wajnflasz-Pick model is more efficient than the conventional spin-1/2 model in the region where ωh z i > 0 (Fig. 3). The maximum value of T * in the effective Wajnflasz-Pick model has larger values than that of the spin-1/2 model, where the effective additional field ωh i (s) is in the opposite direction to the original longitudinal field h z i (Panels (c) and (d) in Fig. 8). It is consistent with the case where the quantum Wajnflasz-Pick model is less efficient than the conventional spin-1/2 model in the region where ωh z i < 0 (Fig. 3).
In order to perform the scaling analysis of the minimum energy gap ∆ min ≡ min[E 1 (s) − E 0 (s)], we consider the p-spin model in the absence of the longitudinal magnetic field: where the transverse magnetic field is homogeneous. Replacement of τ x,y i with σ x,y i provides the conventional p-spin model, where the first order phase transition emerges, and the minimum energy gap is known to close exponentially as N increases in the case where p is odd [15]. After mapping to the subspace spanned by the spin-1/2 model, the reduced Hamiltonian of the quantum Wajnflasz-Pick model with (g u , g l ) = (2, 1) can be reduced tô up to the constant energy shift, where Γ z ≡ ωh x /(2c),  Figure 9 clearly shows that the minimum energy gap closes exponentially in the conventional spin-1/2 model, and the gap closes polynomially in the model reduced from the quantum Wajnflasz-Pick model. This polynomial gap closing originates from the emergence of the effective longitudinal magnetic field in the reduced model: Γ z = ωh x /(2c) = 0.

V. RANDOM COUPLING
In the random spin-spin coupling case, where J ij are randomly generated by the gaussian distribution function [50] the density plot of the mean-value of the success probability is similar to the uniform coupling case. The maximum (minimum) value of the success probability is, however, suppressed (increased) compared with the uniform coupling case (Fig. 10). The variances of the success probability of the quantum Wajnflasz-Pick model are almost ranged from 0.03 to 0.06 in the first and third orthants in the ω-h plane, where the higher success probability may be obtained than the conventional spin-1/2 model. They Averaged success probability of a quantum Wajnflasz-Pick model P scaled by that of spin-1/2 model P 1/2 in a randomly generated coupling strength case. We employed the coupling strength Jij randomly generated from the gaussian distribution function, where the mean is zero and the variance is 1/N . We used N = 4, h x i = 1 and T = 10. The success probabilities P and P 1/2 are averaged values of 100 samplings in each data point. are almost ranged from 0.02 to 0.15 in the second and forth orthants in the ω-h plane, where the lower success probability may be obtained. In the spin-1/2 model, the variance of the success probability is almost within the range from 0.03 to 0.06 in all the orthants.
The discussion above is in the case for a uniform longitudinal magnetic field. In the following, we discuss the case of random longitudinal magnetic fields h z i in addition to the random interactions J ij . The success probabilities P and P 1/2 are almost equal in the weak internal state coupling case (ω = ±0.1 in Fig. 11). In the strong internal state coupling case (ω = ±1 in Fig. 11), the distribution is broaden. Although we can find cases where the conventional spin-1/2 model is superior to the quantum Wajnflasz-Pick model, we can also find many cases where the quantum Wajnflasz-Pick model is superior to the conventional spin-1/2 model, where the success probability is close to the unity compared with the conventional spin-1/2 model.
In these random coupling cases, it may not be definitely concluded that the quantum Wajnflasz-Pick model is always more efficient than the conventional spin-1/2 model. The variance is relatively large, and there are cases where the quantum Wajnflasz-Pick model is inferior to the conventional spin-1/2 model (Fig. 11). However, we can find many cases where the quantum Wajnflasz-Pick model is possibly more efficient than the conventional spin-1/2 model. In the quantum Wajnflasz-Pick model and its reduced model, we have chances to find a better solution of the combinatorial optimization problem. In real annealing machines, we can extract a better solution after performing many sampling experiments by tuning ω.

VI. DISCUSSION
In the case where (g u , g l ) = (2, 1), the spin matrix in the quantum Wajnflasz-Pick model is represented by a (3 × 3)-matrix, which suggests that the quantum Wajnflasz-Pick model in this case may be mapped into the model represented by the spin-1 matrices given bŷ Indeed, after we interchange elements of second and third rows in the spin matrices defined in Eq. (9) in the case where (g u , g l ) = (2, 1), as well as we interchange elements of second and third columns, simultaneously, we find the following mapŝ where we have introduced quadrupolar operators [51, 52] and ω ( ω) is the real (imaginary) part of ω. Since [q z , (Ŝ x ) 2 ] = 0 and [q x , (Ŝ x ) 2 ] = i( ω/c)Ŝ x hold, we find that (Ŝ x ) 2 is the operator of the conserved quantity in the case where the parameter ω is a real number. The coupling ofτ z iτ z j( =i) is mapped into the interactionq z iq z j( =i) , which is a kind of the biquadratic interaction with respect to the spin. In short, the interacting quantum Wajnflasz-Pick model with (g u , g l ) = (2, 1) can be mapped into the spin-1 model with an artificial biquadratic interaction. In particular, in the case where ω ∈ R, there is the hidden symmetry related to (Ŝ x ) 2 , which indicates that the quantum Wajnflasz-Pick model is reducible in this case.
It is general that an interacting quantum Wajnflasz-Pick model is reducible to the conventional spin-1/2 model. It holds for an arbitrary number of the degeneracy (g u , g l ) and at an arbitrary time s, which can be proven in the case where the parameter ω is a real number and the condition ω > −1 holds. In Supplementary  Information, we show that the Hamiltonian of the interacting quantum Wajnflasz-Pick model with arbitrary (g u , g l ) can be projected to the spin-1/2 model, and the initial ground state in the original quantum Wajnflasz-Pick Hamiltonian is also projected to the reduced Hilbert space. It indicates that the quantum annealing in the quantum Wajnflasz-Pick model can be always described by the reduced Hamiltonian.
As shown in Supplementary Information, this projection holds not only in the 2-body interacting quantum Wajnflasz-Pick model, but also in the N -body interacting model. It indicates that if the quantum Wajnflasz-Pick model is embedded into the Lechner-Hauke-Zoller (LHZ) architecture [53,54], it can be also projected into the LHZ architecture composed of the spin-1/2 model, where the effective additional magnetic fields may emerge. The present quantum Wajnflasz-Pick model is a degenerate two-level system in the presence of the transverse magnetic field. The possibility of the implementation of the degenerate two-level system has been discussed for the D 2 line of 87 Rb [41,42]. The quantum Wajnflasz-Pick model is also similar to the ∆-type cyclic artificial atom in the superconducting circuit [38,43]. In the ∆-type artificial atom, the population is controllable by making use of the amplitudes and/or phases of microwave pulses, where the amplitudes alone controls the population in the conventional three-level system (Λ-type system) [43]. However, the ∆-type system in the superconducting circuit is not an exactly degenerate two-level system. With this regard, it may be difficult to directly implement our model in the ∆-type cyclic artificial atom in the superconducting circuit. Actually, it may be feasible to employ the spin-1/2 model with the scheduling function inspired by the quantum Wajnflasz-Pick model, in the case where the Schrödinger dynamics without the dissipation holds.
The quantum Wajnflasz-Pick model is one of the qudit models, which is a kind of the artificial ∆-type system [38,43] in the case where (g u , g l ) = (2, 1). The question naturally arises whether the Λ-type system also shows the higher success probability than the conventional spin-1/2 model. The spin matrix of the Λ-type system we employ here is given bŷ where we take |ε| ≤ 1, and the coefficient c ≡ √ 1 + κ 2 is a normalization factor so as the maximum eigenvalues ofτ x,z are unity. The Hamiltonian of the quantum annealing with the Λ-type system is given by Eqs. (1), (4), and (5), whereτ x,z are replaced with those given in (32). The success probability in the Λ-type system is found to be higher than that in the conventional spin-1/2 model, in the case where ε is small in the weak longitudinal magnetic field region, which is similar to the case of the quantum Wajnflasz-Pick model (Panels (a) and (b) in Fig. 12). When is large, on the other hand, the success probability is drastically suppressed (Panel (c) in Fig. 12). In the case of a single Λ-spin system with = 0, which corresponds to a degenerate two-level system, the unitary transformation can map the HamiltonianĤ(s) = −sh zτ z − (1 − s)h xτ x to the following block diagonal form: As a result, after exchanging the first and second columns and also the first and second rows, we may reduce a quantum annealing problem in this Λ-spin model into that of the spin-1/2 model, the Hamiltonian of which is given byĤ(s) = −sh z σ z /2 − (1 − s)h x σ x − sh z /2. Although the Λ-type system may provide the higher success probability than the conventional spin-1/2 model, the effect of dark states (never employed states) on the quantum annealing in the general Λ-spin case and its reduction to the spin-1/2 model in the many-spin system would be important issues for future study.
To summarize, we have demonstrated that qudit models, such as the quantum Wajnflasz-Pick model as well as the Λ-type system, may provide the higher success probability than the conventional spin-1/2 model in the weak magnetic field region. We have analytically shown that the quantum Wajnflasz-Pick model can be reduced into the spin-1/2 model, where effect of the transverse magnetic field in the original Hamiltonian emerges as the effective additional longitudinal magnetic field in the reduced Hamiltonian, which possibly opens the energy gap between the ground state and the first excited state in the reduced Hamiltonian. Since qubits have experimental advantages for the manipulation, the direct implementation of the reduced spin-1/2 model may be convenient for the quantum annealing. On the other hand, the reduction to the subspace in terms of the spin-1/2 model is useful only in the case where we focus on the Schrödinger dynamics. If we consider the dissipation as a realistic system, the transition between the subspaces emerges. The efficiency of the quantum annealing in this system is open for further study.

VII. CONCLUSIONS
We studied the performance of the quantum annealing constructed by one of the degenerate two-level systems, called the quantum Wajnflasz-Pick model. This model shows the higher success probability than the conventional spin-1/2 model in the region where the longitudinal magnetic field is weak. The physics behind this is that the quantum annealing of this model can be reduced into that of the spin-1/2 model, where the effective longitudinal magnetic field in the reduced Hamiltonian may open the energy gap between the ground state and the first excited state, which gives rise to the suppression of the Landau-Zener transition. The reduction of the quantum Wajnflasz-Pick model to the spin-1/2 model is general at an arbitrary time as well as in an arbitrary number of degeneracies. We also demonstrated that the Λ-type system also shows the higher success probability than the conventional spin-1/2 model in the weak magnetic field regions. We hope that studying quantum annealing with variant spins, and utilizing the insight of their reduced model will promote further development of high performance quantum annealer.

ACKNOWLEDGEMENT
We thank R. van Bijnen, W. Lechner, Y. Matsuzaki, T. Ishikawa, T. Yamamoto, and T. Nikuni for fruitful discussions and comments. Two of the authors (S.W. and S.K.) were supported by Nanotech CUPAL, Japan Science and Technology Agency (JST). Y.S. and S.K. were supported by the New Energy and Industrial Technology Development Organization (NEDO), Japan.

SUPPLEMENTAL INFORMATION
We have shown in (15) that the quantum Wajnflasz-Pick model with the N spins with (g u , g l ) = (2, 1) can be reduced into the spin-1/2 model. In the single-spin case, in particular, if we take the Hamiltonian (1) projected into the Hilbert space spanned by |φ u,l is represented bŷ It can be clearly found that this representation is exactly the same form as the block matrix shown in (12).
In this Supplementary Information, by generalizing this idea, we show that the quantum annealing problem of the quantum Wajnflasz-Pick model with the N -spins with arbitrary number of degeneracy g u,l can be reduced into the spin-1/2 model. First, consider the Hamiltonian of the quantum Wajnflasz-Pick model can be expanded in the Maclaurin series. Equation (39) is a generalization of (1) with (4) and (5). We here introduce eigenstates of a single spinτ z i at a site i as where |u k i for k = 1, 2, · · · , g u and |l k i for k = 1, 2, · · · , g l are eigenstates whose eigenvalues ofτ z i is +1 and −1, respectively.
In the following, we first prove a lemma I: if the parameter ω is a real number, the Hamiltonian of the N -spin quantum Wajnflasz-Pick model can be decomposed into two partsĤ whereP ≡ N i=1P i withP i = |φ u i φ u | i + |φ l i φ l | i is a projection operator. A local projection operatorP i is spanned by two bases |φ u,l i , where It indicates that the Hilbert space of the quantum Wajnflasz-Pick model can be reduced to a subspace spanned by |φ u,l i . We also prove a lemma II: if the parameter ω is a real number and the condition ω > −1 holds, the ground state of the initial HamiltonianĤ x with h x i > 0 is an element of the Hilbert space spanned by |φ u,l i . According to these two lemmas I and II, in the case where ω ∈ R and ω > −1, the quantum annealing problem in the quantum Wajnflasz-Pick model is represented as a model where the local spin has two states-the spin-1/2 model. We first consider the lemma I. SinceP 2 =P , a necessary and sufficient condition providing (42) is given by The condition (45) can be reduced into P ,Ĥ = s P ,Ĥ z + (1 − s) P ,Ĥ x = 0.
Here, we will easily prove that [P ,Ĥ z ] = 0, for the projection operatorP is composed of bases that diagonalizê H z . Indeed, sinceĤ z is composed ofτ z i , what we need to show is A spectral representation of the single-spin operatorτ z i can be represented by eigenstates ofτ z i , given in the form The projection operatorP i is also composed of eigenstates ofτ z i . We can thus immediately conclude that (47) holds, which provides [P ,Ĥ z ] = 0. The result (47) can also provide the following representation The remain we need to prove is The necessary and sufficient condition for (50) for arbitrary h x j is because we can expand a term in (50) in the following way In order to show (51), it is convenient to introduce the spectral representation of the single-spin operatorτ x i , given byτ Here, a constant c is the normalization factor such that the spectral norm ofτ x i is to be unity. We can also representτ x i aŝ Here, ω and ω in (56) are the real and imaginary parts of ω, respectively. By using the representation (56), we can obtain the following result: From the result (60), we find that (51) holds in the case where ω is a real number: ω = 0, which provideŝ As a result, (45) is found to be hold, and the Hamiltonian H are reducible and block diagonalizable, independent of the time s as well as a structure of the spin coupling. Finally, we will prove the lemma II, where the ground state of the initial HamiltonianĤ x with h x i > 0 is an element of the Hilbert space spanned by |φ u,l i , if the conditions ω ∈ R and ω > −1 hold. In the following, we assume that the parameter ω is a real number. Since the initial driver HamiltonianĤ x is a sum of a single-site spin operator −h x iτ x i with h x i > 0, let |Ψ(s = 0) ≡ N i=1 |ψ 0 i be the initial ground state, where |ψ 0 i is the ground state of −h x iτ x i with h x i > 0. In order to prove the lemma II,