Generation of multipartite entanglement between spin-1 particles with bifurcation-based quantum annealing.

Quantum annealing is a way to solve a combinational optimization problem where quantum fluctuation is induced by transverse fields. Recently, a bifurcation-based quantum annealing with spin-1 particles was suggested as another mechanism to implement the quantum annealing. In the bifurcation-based quantum annealing, each spin is initially prepared in [Formula: see text], let this state evolve by a time-dependent Hamiltonian in an adiabatic way, and we find a state spanned by [Formula: see text] at the end of the evolution. Here, we propose a scheme to generate multipartite entanglement, namely GHZ states, between spin-1 particles by using the bifurcation-based quantum annealing. We gradually decrease the detuning of the spin-1 particles while we adiabatically change the amplitude of the external driving fields. Due to the dipole-dipole interactions between the spin-1 particles, we can prepare the GHZ state after performing this protocol. We discuss possible implementations of our scheme by using nitrogen vacancy centers in diamond.

Introduction

Quantum annealing (QA) is a technique for solving combinational optimization problems[1-3]. The solution of the combinational optimization problems is embedded in the ground state of the Ising Hamiltonian[4], which is called the problem (or target) Hamiltonian. We use the transverse magnetic fields to induce quantum fluctuation, and this Hamiltonian is called the driving Hamiltonian. After preparing a ground state of the driving Hamiltonian, we gradually decrease the amplitude of the transverse driving fields while we slowly increase the strength of the Ising Hamiltonian. If the dynamics is adiabatic, the ground state of the problem Hamiltonian can be prepared[5]. Previous studies mainly focus on the use of two-level systems for QA[6-9].

The other mechanisms using bifurcation were proposed to induce the quantum fluctuations for QA. It is known that a parametrically driven Kerr nonlinear oscillator (KPO) shows the bifurcation[10]. A quantum superposition of two distinct states of the KPO can be generated by using quantum adiabatic evolution through its bifurcation point. Moreover, we can use this system as a qubit for a gate type-quantum computer[11]. Previous researches reveal that we can use the KPO for QA to find a ground state of Ising Hamiltonians[12,13].

Recently, Takahashi shows that we can use spin-1 particles for the bifurcation-based QA[14]. For non-interacting spin-1 systems, the initial state is , and degenerate states are prepared at the end of the evolution, which is similar to the bifurcation mechanism of the KPO. On the other hand, for interacting spin-1 systems, the problem Hamiltonian is encoded in a subspace spanned by . Each spin-1 particle is initially prepared in , and adiabatic changes of the Hamiltonian including the coupling between the spin-1 particles provide a ground state of the problem Hamiltonian[14].

Here, we propose a scheme to generate the GHZ states between spin-1 particles by using the bifurcation-based QA. Suppose that there are dipole-dipole interactions between the spin-1 particles. By choosing suitable parameters, the GHZ states have the lowest energy. This means that, starting from a trivial ground state of with longitudinal fields, we adiabatically change the Hamiltonian, and we can obtain the GHZ states where we add external transversal fields in the middle of the dynamics. Importantly, due to the degeneracy of the ground states of the target Hamiltonian, the energy gap between the ground state and excited states becomes small during QA. However, we show that the total Hamiltonian commutes with a parity operator, and this symmetry can suppress the non-adiabatic transitions during QA. Although this kind of the symmetry protected mechanism was discussed in the conventional QA[15-19], we firstly utilize the symmetry protected mechanism for the bifurcation-based QA. Moreover, as a possible implementation, we discuss the use of nitrogen vacancy (NV) centers in diamond, and they are spin-1 particles that are candidates to realize quantum information processing.

The paper is structured as follow. In “Quantum annealing” section, we review the conventional QA and bifurcation-based QA to find a ground state of the Ising Hamiltonian. In “The nitrogen vacancy centers in diamond” section, we review the NV ceners in diamond. In “Generation of the GHZ states with bifurcation-based quantum annealing” section, we introduce our scheme to generate the GHZ states with the bifurcation-based QA. In “Numerical simulations to generate GHZ states with bifurcation-based quantum annealing” section, we perform numerical simulations to evaluate the performance of our scheme. In “Conclusion” section, we summarize our results.

Quantum annealing

Conventional quantum annealing with spin-1/2 particles

Here, we review the conventional QA with spin-1/2 particles[1-3]. The main aim of QA is to prepare a ground state of the following Ising-type Hamiltonian.where L denotes the number of spins, denotes a longitudinal field at the j-th spin, and denotes the coupling strength between the i-th spin and j-th spin. We also use a driver Hamiltonian to induce the quantum fluctuation as follows.where denotes transverse fields. The total Hamiltonian is described as follows.where T denotes the time to implement QA. In QA, we prepare a ground state of , and let this state evolve by the total Hamiltonian. It is known that, as long as an adiabatic condition is satisfied, we can obtain a ground state of the total Hamiltonian.

Bifurcation-based quantum annealing with spin-1 particles

Let us review a bifurcation-based quantum annealing with spin-1 particles[14]. We consider the following driving Hamiltonianwhere , , , and , . We slowly change C(t) from a positive large value to a negative large value while A(t) has a finite but a small value in the middle of QA. The problem Hamiltonian is given asand the total Hamiltonian is given as

We set , and the ground state of the total Hamiltonian at is . By letting this state evolve by the total Hamiltonian, we obtain the ground state of the problem Hamiltonian as long as the adiabatic condition is satisfied.

The nitrogen vacancy centers in diamond

We review the Hamiltonian of the NV centers in diamond. The NV center is a spin-1 patricle, and there is a dipole-dipole interaction between the NV centers. The Hamiltonian is described as followswhere denotes a zero-field splitting at the j-th spin, denotes a strain at the j-th spin, denotes the flip-flop interaction between the j-th spin and k-th spin, and denotes the Ising interaction between the j-th spin and k-th spin. It is worth mentining that we can change the values of () by changing the temperature (amplitude of the applying electric fields)[20-24].

The NV center is a promising candidate to realize quantum information processing. The NV center can be coupled with magnetic fields, electric fields, and temperature, and pressure[22,25,26]. We can polarize the NV centers by illuminating a green laser, and also we can readout the spin state by using the photoluminescence from the NV centers[25-27]. Moreover, the NV center has a long coherence time such as a few milliseconds[28-30]. The NV center can be coherently coupled with an optical photon[31]. These properties are prerequisite for the NV centers to be candidates for the quantum sensing[25,32-34], a quantum memory for a superconducting qubit[35-38], quantum communication[39], and distributed quantum computation[40]. Especially, NV centers could be used to realize an entanglement-enhanced quantum sensing with the GHZ states[41-48] or could be used for a quantum network with encoding where the GHZ states are resource to construct an error correcting code[49-51].

Generation of the GHZ states with bifurcation-based quantum annealing

We explain our scheme to generate a GHZ state between spin-1 particles with the bifurcation-based QA. The GHZ states are defined as . The schematic is shown in Fig. 1a. We consider to apply our scheme with the NV centers in diamond. Importantly, in an experiment, it is difficult to have a negative value of . Although we can slightly change the value of by changing the temperature, the value of is as large as GHZ, and there is no experiment to change the value of the zero field splitting to the negative values, which requires a frequency shift of a few GHZ. To overcome this problem, we adopt an idea of a spin-lock QA where the system driven by microwave fields is in a rotating frame[52-54]. The advantage of this scheme is that the detuning between the resonant frequency of the spins and the microwave frequency plays an role of the longitudinal fields, and we can easily set the negative detuning by setting a suitable value of the microwave frequency. When the NV centers are arranged in a one dimensional chain and microwave driving field are applied along x direction, the Hamiltonian is described as follows.where () denotes the amplitude (frequency) of the microwave driving at the j-th NV center. The dipole-dipole interactions decrease by where denotes the distance between the spins. For example, is satisfied. In a rotating frame defined by , we obtainwhere we define and we use a rotating wave approximation (RWA). In the real experiments, we can easily change the frequency of the microwave driving while the dynamical control of the zero-field splitting is difficult. So we assume that is constant while we change during QA. Throughout of our paper, we set the following.where we set . At (), the ground state of the Hamiltonian after the RWA is approximately described by for . On the other hand, at , degenerate ground states of the Hamiltonian after the RWA are described by for and . We plot an energy diagram of the Hamiltonian (9) in Fig. 1b, and we confirm that the energy gap between the ground state and first excited state becomes smaller as the time t approaches to T.

Figure 1

(a) Schematic of our proposal. Spin-1 particles with zero-field splitting are arranged in a one-dimensional chain. There are dipole-dipole interactions between the spin-1 particles. We globally apply time-independent microwave fields to implement the bifurcation-based quantum annealing. (b) We plot eigenenergies against the time with the unit of milliseconds where we adopt the Hamiltonian in the Eq.(9), which is derived by using a rotating wave approximation. The energy gap between the ground state and first excited state becomes smaller as t approaches to T. We set the parameters as , kHz, kHz, kHz, kHz, kHz, kHz, , ms.

Importantly, the Hamiltonian in the Eq. (9) commutes with a parity operator of , and we have while we have . Therefore, by preparing a state of , the adiabatic change in the Hamiltonian allows us to create the state of where non-adiabatic transitions between and are prohibited due to the difference of the symmetry.

Numerical simulations to generate GHZ states with bifurcation-based quantum annealing

To evaluate the performance of our scheme, we perform numerical simulations to plot the fidelity between the target GHZ state and the state after QA. Here, we adopt the Hamiltonian in the Eq. (8). To consider the decoherence, we use the following GKSL master equation[55,56]where denotes a decoherence rate and denotes a lindblad operator at the j-site. Throughout of this paper, we use , which corresponds to magnetic field noise that is typical for the NV centers[57-60]. We define a fidelity as .

We plot the fidelities against t for without decoherence in Fig. 2a. When there is no strain, the fidelity is more than 0.999, and this means that the adiabatic condition is reasonably satisfied. When we add the effect of the strain, the fidelity becomes as small as 0.979 (0.925) for () kHz, as shown in Fig. 2a. This comes from the fact that a ground state of the Hamiltonian with the strain is not the GHZ state. To obtain a high-fidelity GHZ state among the NV centers, it is crucial to suppress the effect of the strain by applying suitable amount of the electric fields. In the real experiment, we have GHz. However, the computational cost becomes expensive when is much larger than the other parameters. Therefore, throughout of this paper, we set MHz. Since we confirm that the dynamics does not significantly change even when we increase and around this parameter range, we believe that our numerical simulations are still useful to predict the experimental results for GHz.

Figure 2

(a) We plot a fidelity against the time with the unit of milliseconds where we adopt the Hamiltonian without a rotating wave approximation. We set the parameters as kHz, kHz, , kHz, kHz, kHz, , ms, , and MHz. (b) We plot a fidelity against the time under the effect of decoherence with kHz. We use the same parameters as that used in (a) except the decoherence rate.

Also, we plot the fidelity under the effect of decoherence against t for in Fig. 2b. Compared with the fidelity by using the unitary dynamics plotted in Fig. 2a, the fidelity becomes smaller as expected. However, the fidelity is still around 0.9, and so these results show that we can generate the GHZ states even under noisy environments.

Importantly, there was an experimental demonstration to generate an entanglement between two NV centers[61]. However, the previous scheme requires a complicated pulse sequence, and the necessary number of the pulse operations increases as the number of NV centers increases. Moreover, the NV centers should be individually controlled by using frequency selectivity. On the other hand, our protocol just requires global applications of the microwave pulses without individual adressing of the NV centers, which would be beneficial to generate a GHZ states with more than two NV centers.

Moreover, we plot the fidelities against t for with and without decoherence, as shown in Fig. 3a,b, respectively. When we consider the unitary dynamics, the fidelities with are comparable with those with , as shown in Fig. 3a. This means that, for , the adiabatic conditions are reasonably satisfied. With decoherence, the fidelities becomes worse than those without decoherence. However, as shown in Fig. 3b, the fidelities are still around 0.9. Again, these results show the practicality of our scheme.

Figure 3

(a) We plot a fidelity against the time with the unit of milliseconds where we adopt the Hamiltonian without a rotating wave approximation. We set the parameters as kHz, kHz, , , kHz, kHz, , kHz, , ms, , and MHz. (b) We plot a fidelity against the time under the effect of decoherence with kHz. We use the same parameters as that used in (a) except the decoherence rate.

Finally, we plot the fidelity agaisnt the number of qubits. For this purpose, let us firstly show the fidelity with in Fig. 4a. To save the computational cost, we adopt the Hamiltonian with a rotating wave approximation. This result shows that the fidelity is as high as 0.89. Then, we plot the fidelity against the number of qubits in Fig. 4b. Roughly speaking, the fidelity decreases by 0.03 as we add another NV center.

Figure 4

(a) We plot a fidelity against the time with the unit of milliseconds where we adopt the Hamiltonian with a rotating wave approximation. We set the parameters as kHz, kHz, kHz, kHz, , , kHz, , ms, and, kHz. (b) We plot the fidelity agains the number of qubits where we adopt the Hamiltonian with a rotating wave approximation. We choose kHz for . For and , we use the same parameters as that used in Figs. 2 and 3 except the strain. For , we use the same parameters as that used in (a).

Conclusion

In conclusion, we propose a scheme to generate GHZ states between spin-1 particles by using bifurcation-based QA. Suppose that there are dipole-dipole couplings between the spin-1 particles. After each spin-1 particle is prepared in , we slowly turn on the microwave driving, and we finally turn off the the microwave driving in an adiabatic way. We show that adiabatic changes in frequency and amplitude of the microwave driving fields provide a GHZ states after QA. Although the energy gap between the ground state and first excited state becomes nearly degenerate when we turn off the microwave driving fields, we show that a symmerty of the Hamiltonian protects the state from the non-adiabatic transitions. Our scheme could be useful for possible applications to quantum information processing by using nitrogen vacancy centers in diamond. Our method could be extended to the preparation of other entangled states such as spin cat states (that are considered to be generalized GHZ states[62-64]) or W states[65]. There are four key points to implement our scheme. First, we need to design a problem Hamiltonian whose ground state is the target entangled state. Second, the Hamiltonian should commute with a symmetry operator (such as a parity operator). Third, the ground state of the driver Hamiltonian should belong to the same sector of the symmetry operator as the ground state of the problem Hamiltonian. Fourth, the first excited state of the problem Hamiltonian should belong to a different sector of the symmetry operator as the ground state of the problem Hamiltonian. If these conditions are satisfied, our method can be used to prepare other entangled states. Detailed investigation whether these conditions are satisfied for other entanglement is beyond the cope of the paper, and therefore we leave this as a future work.

The datasets generated and/or analysed during the current study are not publicly available due internal rules of the AIST but are available from the corresponding author on reasonable request.