Shohei Watabe1,2, Yuya Seki3, Shiro Kawabata3. 1. Department of Physics, Faculty of Science Division I, Tokyo University of Science, Shinjuku, Tokyo, 162-8601, Japan. shoheiwatabe@rs.tus.ac.jp. 2. Nanoelectronics Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki, 305-8568, Japan. shoheiwatabe@rs.tus.ac.jp. 3. Nanoelectronics Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki, 305-8568, Japan.
Abstract
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 Λ-type system, which may show the higher success probability than the conventional spin-1/2 model.
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 Λ-type system, which may show the higher success probability than the conventional spin-1/2 model.
Quantum annealing is an interesting approach for finding the optimal solution of combinatorial optimization problems by using the quantum effect[1-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], several hardwares have been investigated and developed[9-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-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-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) eigenenergies, and (iii) non-uniform effects of the spin-interaction 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-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-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.
Quantum Wajnflasz–Pick Model
A conventional quantum annealing consists of the spin-1/2 model, where the time dependent Hamiltonian is given by[1]where are a problem and driver Hamiltonian, respectively, and is the time scaled by the annealing time T. The problem Hamiltonian with the number of spins N, which encodes the desired optimal solution, has a non-trivial ground state. In contrast, the driver Hamiltonian has a trivial ground state, where the driver Hamiltonian must not be commutable with the problem Hamiltonian . A problem Hamiltonian and driver Hamiltonian are typically given bywhere are the Pauli matrices, is the coupling strength between spins, is the local longitudinal magnetic field, and is the local transverse magnetic field. The time-dependent total Hamiltonian gradually changes from the driver Hamiltonian to the problem Hamiltonian . 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 in Eqs. (2) and (3) with the spin matrices of the quantum Wajnflasz-Pick model . The spin operator is given by[30]where is the number of the degeneracy of the upper (lower) states. The spin-operator in the driver Hamiltonian is given bywhere is a matrix with the off-diagonal term , given by
Figure 1
Schematic setup of an interacting degenerate two-level system, called the quantum Wajnflasz–Pick model.
Schematic setup of an interacting degenerate two-level system, called the quantum Wajnflasz–Pick model.Here, is a parameter of the internal transition between the degenerated upper/lower states. The matrix is the matrix, where all the elements is unity, which gives the transition between the upper and lower states. The constant is the normalization factor, where the maximum eigenvalue is normalized to be +1, so as to be equal to the maximum eigenvalue of .In the following, for the consistency to the earlier work[30], we consider a uniform transverse field , and also take the parameter of the internal transition to be real with . In the case where , we havewith , 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 , the magnetic fields , and the annealing time . 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.
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 equationIn this method, the time-evolution of the wave function is calculated by using the Cayley’s form[47]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 , which is consistent with the earlier work[30]. For example, in the case where and for , the time-dependence of the ground state population of the problem Hamiltonian, given by , clearly shows that this quantity in the quantum Wajnflasz–Pick model is greater than that in the conventional spin-1/2 model (Panels (a) and (b) in Fig. 2). Here, is the ground state of the problem Hamiltonian, and is the wave function obtained from the time-dependent Schrödinger equation. In the case where and for , on the other hand, the ground state population of the problem Hamiltonian in the quantum Wajnflasz–Pick model is less than that in the spin-1/2 model (Panels (c) and (d) in Fig. 2).
Figure 2
Population of the ground state of the problem Hamiltonian in the quantum Wajnflasz–Pick model and that of the conventional spin-1/2 model, where . The scaled time is given by . We employed the number of spins both in the quantum Wajnflasz–Pick model and in the spin-1/2 model. We used the parameters , , and .
Population of the ground state of the problem Hamiltonian in the quantum Wajnflasz–Pick model and that of the conventional spin-1/2 model, where . The scaled time is given by . We employed the number of spins both in the quantum Wajnflasz–Pick model and in the spin-1/2 model. We used the parameters , , and .Compare the success probability of the quantum Wajnflasz–Pick model, , with that of the conventional spin-1/2 model denoted as , where is the final state obtained from the time-dependent 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 , the energy gap is known to close exponentially and the first-order phase transition emerges in the absence of the longitudinal 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.
Figure 3
Success probability of a quantum Wajnflasz–Pick model scaled by that of the conventional spin-1/2 model as a function of longitudinal magnetic field h and the coupling strength between degenerate internal states. The parameters are the same as those used in Fig. 2.
Success probability of a quantum Wajnflasz–Pick model scaled by that of the conventional spin-1/2 model as a function of longitudinal magnetic field h and the coupling strength between degenerate internal states. The parameters are the same as those used in Fig. 2.In the case where , 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 , the success probability of the quantum Wajnflasz–Pick model is almost equal to that of the case where , 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).
Figure 4
Success probability of a quantum Wajnflasz–Pick model scaled by that of spin-1/2 model . We consider the following degeneracy case: in Panels (a) and in Panels (b). We used , , and .
Success probability of a quantum Wajnflasz–Pick model scaled by that of spin-1/2 model . We consider the following degeneracy case: in Panels (a) and in Panels (b). We used , , and .
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 , 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 , 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 . 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.
Figure 5
Eigenenergies of the instantaneous Hamiltonian in the quantum Wajnflasz–Pick model (blue) and those in the reduced spin-1/2 model (black) as a function of . The parameters are the same as those in Fig. 2.
Eigenenergies of the instantaneous Hamiltonian in the quantum Wajnflasz–Pick model (blue) and those in the reduced spin-1/2 model (black) as a function of . The parameters are the same as those in Fig. 2.For example, in the case where , the single-spin Hamiltonian in the quantum Wajnflasz–Pick model is decomposable, where the irreducible representation is given byfor arbitrary values of s, by using the unitary operatorwhere , and . 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 formSince the initial ground state of the single-spin Hamiltonian is given by in the original quantum Wajnflasz–Pick model, this state can be mapped to . It indicates that the initial ground state can be also projected to the Hilbert space of the reduced Hamiltonian .This reduction of the single-spin problem in the case where 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 formwherewith
Figure 6
Schematics of an original quantum Wajnflasz–Pick model (a) and its reduced model (b).
Schematics of an original quantum Wajnflasz–Pick model (a) and its reduced model (b).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 in the reduced Hamiltonian is the same as that of the original Wajnflasz–Pick model. The effective longitudinal magnetic field in the reduced Hamiltonian also reaches the same value as that of the original Wajnflasz–Pick model at the end of the annealing: . 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 is a positive real number such that the maximum eigenvalue of is unity, and we take . Then, always holds during the annealing time . In the case where the longitudinal magnetic field is very large, , the effect of the original longitudinal magnetic field is dominant compared with the effective additional term except at the very early stage of the annealing . 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 .In the case where the original longitudinal magnetic field is not large, the effective additional longitudinal magnetic field cannot be neglected compared with . When the effective additional field is in the same direction as the original longitudinal field, the total effective longitudinal magnetic field 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 , 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 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 , 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.
Figure 7
Excited state energies measured from the ground state energy of the instantaneous Hamiltonian in the reduced model (red) and those in the conventional spin-1/2 model (blue). The parameters are the same as those in Fig. 2.
Excited state energies measured from the ground state energy of the instantaneous Hamiltonian in the reduced model (red) and those in the conventional spin-1/2 model (blue). The parameters are the same as those in Fig. 2.Behavior of success probability is also explained by the reference of the annealing time[49]whereHere, and 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 much larger than . Let be an instantaneous reference time of the annealing. The maximum value of this time 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 is in the same direction as the original longitudinal field (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 (Fig. 3). The maximum value of in the effective Wajnflasz–Pick model has larger values than that of the spin-1/2 model, where the effective additional field is in the opposite direction to the original longitudinal field (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 (Fig. 3).
Figure 8
Instantaneous reference annealing time as a function of the scaled time . The parameters are the same as those in Fig. 2.
Instantaneous reference annealing time as a function of the scaled time . The parameters are the same as those in Fig. 2.In order to perform the scaling analysis of the minimum energy gap , we consider the -spin model in the absence of the longitudinal magnetic field:where the transverse magnetic field is homogeneous. Replacement of with provides the conventional -spin model, where the first order phase transition emerges, and the minimum energy gap is known to close exponentially as increases in the case where 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 can be given byup to the constant energy shift, where , , and . By using the commutation relation and by following the standard argument of the angular momentum, where the total spin conserves, the Hilbert space can be spanned by states , where and with . The diagonal elements of this Hamiltonian is given by , and the off-diagonal elements are . Since the ground state of this model is given by the case , we diagonalize the matrix of the reduced Hamiltonian. We compare the minimum energy gap of this model reduced from the quantum Wajnflasz–Pick model with that of the conventional -spin model composed of the spin-1/2 system (Eq. (24) with and ). 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: .
Figure 9
Minimum energy gap as a function of the number of spins on a linear-log scale (a log-log scale in the inset). The minimum energy gap is obtained from the exact diagonalization of the ferromagnetic -spin model with . We have used and .
Minimum energy gap as a function of the number of spins on a linear-log scale (a log-log scale in the inset). The minimum energy gap is obtained from the exact diagonalization of the ferromagnetic -spin model with . We have used and .
Random Coupling
In the random spin-spin coupling case, where 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 - plane, where the higher success probability may be obtained than the conventional spin-1/2 model. They are almost ranged from 0.02 to 0.15 in the second and forth orthants in the - 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.
Figure 10
Averaged success probability of a quantum Wajnflasz–Pick model scaled by that of spin-1/2 model in a randomly generated coupling strength case. We employed the coupling strength randomly generated from the gaussian distribution function, where the mean is zero and the variance is 1/N. We used , and . The success probabilities and are averaged values of 100 samplings in each data point.
Averaged success probability of a quantum Wajnflasz–Pick model scaled by that of spin-1/2 model in a randomly generated coupling strength case. We employed the coupling strength randomly generated from the gaussian distribution function, where the mean is zero and the variance is 1/N. We used , and . The success probabilities and are averaged values of 100 samplings in each data point.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 in addition to the random interactions . The success probabilities and are almost equal in the weak internal state coupling case ( in Fig. 11). In the strong internal state coupling case ( 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.
Figure 11
Success probability of the quantum Wajnflasz–Pick model vs. success probability of the conventional spin-1/2 model. We take 1000 samples of problem hamiltonian with the random coupling strength as well as the random longitudinal magnetic field , both of which are generated from the standard Gaussian distribution. The means of and are zeros, and the variances are 1/N and 1/, respectively. For each problem set, we consider four cases and ±1. We used , , and .
Success probability of the quantum Wajnflasz–Pick model vs. success probability of the conventional spin-1/2 model. We take 1000 samples of problem hamiltonian with the random coupling strength as well as the random longitudinal magnetic field , both of which are generated from the standard Gaussian distribution. The means of and are zeros, and the variances are 1/N and 1/, respectively. For each problem set, we consider four cases and ±1. We used , , and .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 .
Discussion
In the case where , 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 byIndeed, after we interchange elements of second and third rows in the spin matrices defined in Eq. (9) in the case where , as well as we interchange elements of second and third columns, simultaneously, we find the following mapswhere we have introduced quadrupolar operators[51,52]and () is the real (imaginary) part of . Since and hold, we find that is the operator of the conserved quantity in the case where the parameter is a real number. The coupling of is mapped into the interaction , which is a kind of the biquadratic interaction with respect to the spin. In short, the interacting quantum Wajnflasz–Pick model with can be mapped into the spin-1 model with an artificial biquadratic interaction. In particular, in the case where , there is the hidden symmetry related to , 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 and at an arbitrary time s, which can be proven in the case where the parameter is a real number and the condition holds. In Supplementary Information, we show that the Hamiltonian of the interacting quantum Wajnflasz–Pick model with arbitrary 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 -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 line of 87Rb[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 . 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 bywhere we take , and the coefficient is a normalization factor so as the maximum eigenvalues of are unity. The Hamiltonian of the quantum annealing with the Λ-type system is given by Eqs. (1), (4) and (5), where 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 , which corresponds to a degenerate two-level system, the unitary transformationcan map the Hamiltonian to the following block diagonal form:
Figure 12
Success probability of the Λ-type system in the - plane, compared with that of the conventional spin-1/2 model . The number of spin both in the Λ-type system and in the spin-1/2 model is . We used the parameter sets , and .
Success probability of the Λ-type system in the - plane, compared with that of the conventional spin-1/2 model . The number of spin both in the Λ-type system and in the spin-1/2 model is . We used the parameter sets , and .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 . 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.
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.
Authors: R Barends; A Shabani; L Lamata; J Kelly; A Mezzacapo; U Las Heras; R Babbush; A G Fowler; B Campbell; Yu Chen; Z Chen; B Chiaro; A Dunsworth; E Jeffrey; E Lucero; A Megrant; J Y Mutus; M Neeley; C Neill; P J J O'Malley; C Quintana; P Roushan; D Sank; A Vainsencher; J Wenner; T C White; E Solano; H Neven; John M Martinis Journal: Nature Date: 2016-06-09 Impact factor: 49.962