Abstract
When the evolution of discrete time quantum walk is carried out for particles, the ramble state is prone to error due to the influence of system noise. A multiparticle quantum walk error correction algorithm based on the two-lattice Bose–Hubbard model is proposed in this study. First, two point Bose–Hubbard models are constructed according to the local Euclidean generator, and it is proved that the two elements in the model can be replaced arbitrarily. Second, the relationship between the transition intensity and entanglement degree of the particles in the model is obtained by using the Bethe hypothesis method. Third, the position of the quantum lattice is coded and the quantum state exchange gate is constructed. Finally, the state replacement of quantum walk on the lattice point is carried out by switching the walker to the lattice point of quantum error correction code, and the replacement is carried out again. The entanglement of quantum particles in the double-lattice Bose–Hubbard model is simulated numerically. When the ratio of the interaction between particles and the transition intensity of particles is close to 0, the entanglement operation of quantum particles in the model can be realized by using this algorithm. According to the properties of the Bose–Hubbard model, quantum walking error correction can be realized after particle entanglement. This study introduces the popular restnet network as a training model, which increases the decoding speed of the error correction circuit by about 33%. More importantly, the lower threshold limit of the convolutional neural network (CNN) decoder is increased from 0.0058 under the traditional minimum weight perfect matching (MWPM) to 0.0085, which realizes the stable progress of quantum walk with high fault tolerance rate.
1 Introduction
A quantum walk is one of the effective methods to realize quantum computation [1,2]. However, because of the randomness of quantum walk, it is difficult to correct errors when the particles are affected by system noise during particle evolution. The physical system used in this study is a homogeneous system, which extends the two-lattice Bose–Hubbard model to the Bose–Hubbard model, and the arbitrary replacement of two particles does not lead to a new quantum state of the system [3–6].
A quantum walk algorithm can be implemented on many graph structures, and its function is very significant [7–14]. Also, quantum walks have important applications in quantum cryptography and quantum communication [15,16]. Noise research in the quantum walk algorithm is also one of the key points, and general quantum walk calculation based on the Bose–Hubbard model has been proposed [17,18]. On this basis, the multiparticle quantum walk error correction algorithm based on the two-lattice Bose–Hubbard model is further studied. As early as 2003, Shapira et al. used computer simulation to study the influence of unitary noise on one-dimensional quantum walk and showed the variation of probability distribution of quantum walk under the influence of noise [19]. Also, specific noise types are studied, such as the non-Markov continuous time quantum walk algorithm with dynamic noise [20] and the quantum Bernoulli central limit theorem in the quantum walk algorithm [21]. In 2015, Ambainis et al. improved the potential barrier using the Grover algorithm, improved the amplitude of quantum walking coin state, and reduced the impact of noise during particle quantum walking [21,23]. In 2016, Wang et al. proposed to construct time-varying quantum walking with infinite degrees of freedom by using quantum Bernoulli noise [24]. In 2018, Du et al. constructed quantum gates using the quantum walk algorithm under noise, and the algorithm increased the number of addresses on the graph or the ratio of jump strength to potential could improve the coherence time, thus inhibiting the decoherence [25]. In 2019, Claudia et al. address the use of quantum walks as a quantum probe to characterize defects and perturbations occurring in complex, classical, and quantum networks [26]. These algorithms are designed to reduce the influence of noise on particles during the quantum walk.
Similar to the aforementioned algorithm, a multiparticle quantum walk error correction algorithm is proposed to reduce the influence of noise in the quantum walk. The two-lattice Bose–Hubbard model is constructed for the multiparticle quantum ramble, and it is proved that the quantum state does not change after the displacement of any two particles in the model. Using the encoding method of quantum states mentioned in the study in reference [17], entanglement of quantum states in two lattices is generated by controlling the ratio of transition intensity and interaction between particles [27]. Controlling the ratio of model parameters to generate entanglement is easier than other methods [15,28]. The convolutional neural network decoder has been introduced in detail in the study mentioned in reference [29–33]. In this study, the threshold of the quantum walk error correction circuit has been greatly improved by using the convolution operation of the decoder and the improvement of the training speed. The evolution operator of quantum walk is improved and quantum state error correction is realized. Finally, the advantages and disadvantages of the algorithm are analyzed. Figure 1 shows the core idea of the work.
FIGURE 1
The sections of this article are organized in the following manner. In Section 2, we briefly introduce the background knowledge of discrete quantum walks. In Section 3, the two-lattice Bose–Hubbard model is constructed. In Section 4, the position of the quantum lattice is coded and the quantum state exchange gate is constructed. In Section 5, an analysis of error correction performance is performed, and Section 6concludes the study.
2 Discrete quantum walk
A quantum walk is a quantized model for a classical random walk. The discrete time models discussed here include the walker position state and the “coin” state. The position n of the walker is a vector of infinite dimensional Hilbert space , and the basis vector of Hilbert space is {|n⟩: n ∈ Z}, which is called the computational basis for position space.
The evolution of the walk depends on the state of a quantum coin. Suppose the walker is on position |n⟩, after the quantum coin is flipped “heads”, the walker goes to position |n + 1⟩ in the next step. If the coin is flipped “tails”, the walker will go to position |n − 1⟩. Now, the Hilbert space of the whole ramble system iswhere is a two-dimensional Hilbert space of the quantum coin (different walkers have different throwing operators), and the basis vector of Hilbert space is {|0⟩, |1⟩}. The aforementioned ground state is also the computational basis of the quantum coin space, and the computational basis for space is {|i⟩|n⟩, i ∈ {0, 1}, n ∈ Z}.
P is the coin operator that determines the direction of the quantum walk. The operator for the walker to move from position |n⟩ to position |n + 1⟩ or position |n − 1⟩ is called the transition operator S.We calculated the representation of the transition operator S under the computational basis of Hilbert space :
3 Two-lattice Bose–Hubbard model
It is assumed that the graph structure of quantum walk is G, and the vertices V in the graph represent lattice points and satisfy the local Euclidean symmetry [4]. The two-lattice Bose–Hubbard model is used to describe the discrete time quantum walk in multiparticle interaction on the graph and to correct the error of quantum walk.
To construct the two-lattice Bose–Hubbard model, we first introduce the boson creation (annihilation) operator, , where j = 1, 2, …, n, which satisfies the commutation relation:
Then we get the Hamiltonian of the two-lattice Bose–Hubbard model:where t is the real parameter to describe the particle jumping intensity (in this work t = 1), ϵj is the local potential, and is the interaction between particles at the lattice point:where is the interaction between two particles, is the interaction between three particles, U1 and U2 are the strength of the interaction between particles, and .
In Eq. 6, the constraint of the two sets {j1, j2, …, jk} and is that the sum extends over all lattice points, and the two sets are not the same. Also, the operator operation condition is satisfied:where is the number operator at the lattice point j. Operator satisfy the commutation relation:
According to the Hamiltonian of the model, the model can make multiple particles transition at the same time without limiting the transition of two adjacent lattice points. In the study mentioned in reference [17], the limit condition of Hilbert space of n boson is calculated in the 2n dimensional computational Hilbert space , and the location is encoded by subset and , along with an isomorphism , where and are interchangeable. The quantum state in the model system provides a computer state encoding on the n-qubits, . A Hamiltonian was obtained by coupling with states space (=) outside of the computational space, which generated entanglement in . It is proved that any initial state can be mapped to and can be returned to recoded as a valid computed state.
4 SWAP gate operator construction
Assume two lattices V0 and V1 (Figure 2), where V0 is the location of the quantum walk and V1 is the error-correcting coded location. The corresponding two lattices in V0 and V1 represent a qubit whose calculated ground state is |0⟩ and |1⟩, respectively. The corresponding qubits in the lattice must be entangled, and we show how quantum lattices are encoded and how qubits generate entanglement. Suppose the initial n-qubit state |Φ⟩ satisfies the following condition:
FIGURE 2
The function of the four calculated ground-state mapping relations and the creation (annihilation) operators for a 10 dimensional space is
The Hamiltonian of the two-lattice Bose–Hubbard model is reduced towhen the interaction of two particles and the local potential ϵj are considered only. In Eq.6, N is the particle number, and the subscript j is the lattice position, and the parameter tij is
In the two-lattice Bose–Hubbard model, the six lattice points are coded according to the Yang diagram, so the lattice coding sequence is as follows:
According to the aforementioned Eqs 11, 12, 13, we can get the energy matrix E6 by matrix diagonalization. From the energy matrix E6, the eigenvalues of lattice particles and the ground state of the model can be obtained when determining the values of U/t, which is helpful for the establishment of the two-lattice Bose–Hubbard model. Moreover, the entangled degree calculation system mentioned in the study in reference [4] is used to calculate the entangled degree between the ground-state particles:where N is the number of particles, M is the number of lattices, and (ϕ)i is the reduced density matrix at the ith lattice point. The next step is to find the relationship between the specific value U/t and the degree η of entanglement. The lattice energy matrix of the seven lattice points is shown below.
According to the eigenenergy, the Bethe hypothesis method can be obtained, as follows:where − t∑m,nC(m, n)/F(m, n) = 1, m and n (m ≥ n) are the number in the lattice points, when m = n and C(m, n) = 1 and when m ≠ n and C(m, n) = 2. According to Eq. 16, the calculation formula of entanglement is improved:We can get the degree of entanglement between the two lattice points to 1 at U/t → 0, as shown in Figure 3. The SWAP gate is obtained by limiting [4] to U/t ≈ 0:where a and b need to be evaluated according to the degree of entanglement η to ensure the degree of entanglement between two particles.
FIGURE 3
5 Error correction and analysis
5.1 Quantum walk error correction
Because the generator of local Euclidean symmetry is used to construct the two-lattice Bose–Hubbard model, it is necessary to modify the quantum ramble operator. Assuming that the generator in the model (n lattices) is , the modification of quantum walk transfer operator and coin operator is as follows:
Because it is currently in a two-dimensional physical system, the aforementioned formula only takes two generative elements to construct the quantum walk evolution operator.
In the previous section, we showed that the quantum states of two corresponding lattice points can be entangled by controlling the ratio of the parameters U/t. Thus, according to the Properties of the Bose–Hubbard model, when particle qubits wander to lattice point α, it can be exchanged with another lattice point α′, which becomes entangled with the particle qubits through the exchange gate and corrects the error.
In the space of N-qubits, a single-bit quantum measurement will appear at each lattice point when the evolution of the particles is in the discrete time quantum walk. The measurement operator is:where, |v⟩i,0 is the encoded quantum state in the first lattice in Figure 2, which is the state of quantum walk; |v⟩i,1 is the quantum state in the second lattice, which is used to assist quantum error correction. Also, the lattice locations of |v⟩i,0 and |v⟩i,1 are entangled in the model. The quantum state |v⟩i,0 has superposition state in the process of quantum walk:where ψ0,n and ψ1,n are the amplitude, |ξ0⟩ and |ξ1⟩ are the coin state, and |n⟩ is the position state. In the process of quantum walk, there are mainly two kinds of errors, which are phase inversion and bit inversion. However, the amplitude in Eq. 23 is generated by the superposition state, so the detection and correction of walker errors are mainly aimed at bit inversion errors. The quantum walk error correction circuit is shown in Figure 4. In order to suppress the error generated in the quantum walk, the noise error based on the double-lattice Bose–Hubbard model in the quantum walk is reduced by reducing the circuit gate overhead. Taking depolarization noise as an example, the data map is obtained through the training of the decoder. The decoder uses the convolutional neural network [27–29] in the current hot machine learning algorithm as the dominant algorithm to decode the quantum information in the line. The specific decoding training model is shown in Figure 5. By training with different restnet layers, the accuracy and speed of training are improved; the speed is 1/3 higher than before and the accuracy reaches 99.82%.
FIGURE 4
FIGURE 5
5.2 Algorithm analysis
At present, there are many research achievements in quantum walking. For example, the quantum state operator based on the Bose–Hubbard model is given in the study in reference [17], including single-bit operation, entanglement operation, and permutation operation. This study mainly studies the error correction of system noise during a quantum walk.
The algorithm is based on the two-lattice Bose–Hubbard model, and the evolution of the quantum walk algorithm is carried out. In the first section, it is proved that no new quantum state will be generated after the particles are replaced in the model [17]. Based on the model, this algorithm can correct the quantum walk error by using quantum error correction coding. In the second section, according to the Bethe hypothesis method mentioned in reference [4], we can get the relationship between energy eigenvalues, particle transition intensity, particle interaction, and entanglement (Eq. (17)). Therefore, two lattice points can be entangled by controlling the value of U/t, and the step of generating entanglement can be omitted compared with the study mentioned in reference [17].
5.3 Error correction performance analysis
We extracted non-local regularities from noise and performed transfer learning in various tasks. Applying this advantage to the cost of qubits passing through quantum gates can reduce the cost of qubits. The qubits contain auxiliary qubits in the synthetic measurement process, and the logic overhead is the cost of the auxiliary qubits during the synthetic measurement process. Threshold is an effective way to characterize fault tolerance performance. Specifically, applying quantum error correction for efficient quantum computing can suppress the logical error rate to arbitrarily low levels when the physical error rate of the qubits is below a certain threshold. The CNN used in this study is used for experiments. The number of simulated qubits varies from 3.8 × 103 to 1.5 × 105 bits. The specific training data quantification is shown in Table 1. Comparing the overhead of MWPM and CNN optimization with thresholds of 0.0058 and 0.0085. Figure 6A shows that with the increase of the number of qubits both the MWPM overhead and the CNN overhead increase, but the increase in CNN-optimized overhead is significantly lower than that of the MWPM overhead. Figure 6B shows that as the number of qubits increases, the optimization overhead of CNN also gets much lower than the original one; although when the number of qubits is 2.0 × 108, the overhead of the optimized CNN is slightly higher than the original one. At the same time, comparing the results under different thresholds in Figure 6, it can be found that the larger threshold, the greater the overhead of quantum circuit gate.
TABLE 1
| Trainable dataset | Steps | Accuracy (%) | |
|---|---|---|---|
| MWPM | 1.48 × 105 | 4.8 × 104 | 75.388 |
| RestNet7 | 1.37 × 103 | 3.1 × 103 | 84.256 |
| RestNet14 | 2.75 × 103 | 2.9 × 103 | 88.792 |
| RestNet21 | 3.79 × 103 | 1.7 × 103 | 96.753 |
Unified quantization of network layers under different decoders.
FIGURE 6
6 Conclusion
In this study, a multiparticle quantum walk error correction algorithm based on the two-lattice Bose–Hubbard model is proposed. This algorithm controls the proportion of the interaction between particles and the transition intensity of particles to realize the entanglement of quantum state in two lattices and then corrects the quantum walking error by using the invariant property of model particle replacement. Under the condition of the threshold of 0.0085, the restnet network layer is used as a training model, and the error noise of quantum walk is reduced by decoding the error correction circuit model, so as to achieve a more stable walk circuit. Compared with the traditional quantum walk error correction, the threshold limit is increased from 0.0058 under the traditional MWPM to 0.0085, and the speed is increased by a full 1/3. However, when the number of particles or the size of the system exceeds a certain threshold, it is impossible to accurately control the transition strength of particles and the interaction between particles, thus destroying the entanglement operation. The development of quantum convolutional neural network is relatively mature at present, and it is the focus of the next research, in preparation for further improving the fault tolerance performance.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
S-MW: conceptualization, methodology, and software; Y-JQ: data curation and writing—original draft preparation; H-WW: visualization, investigation, and supervision; ZC: writing—reviewing and editing; H-YM: conceptualization, funding acquisition, and resources.
Funding
This project was supported by the National Natural Science Foundation of China (Grant No. 61772295), Natural Science Foundation of Shandong Province, China (Grant Nos. ZR2021MF049, ZR2019YQ01), and Project of Shandong Provincial Natural Science Foundation Joint Fund Application (ZR202108020011).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
AharonovYDavidovichLZaguryN. Quantum random walks. Phys Rev A (Coll Park) (1993) 48:1687–90. 10.1103/PhysRevA.48.1687
2.
KempeJ. Quantum random walks: An introductory overview. Contemp Phys (2003) 44:307–27. 10.1080/00107151031000110776
3.
AndersonMHEnsherJRMatthewsMRWiemanCECornellEA. Observation of bose-einstein condensation in a dilute atomic vapor. Science (1995) 269:198–201. 10.1126/science.269.5221.198
4.
ZwierleinMWStanCASchunckCHRaupachSMGuptaSHadzibabicZet alObservation of bose-einstein condensation of molecules. Phys Rev Lett (2003) 91:250401. 10.1103/PhysRevLett.91.250401
5.
KennedyCJBurtonWCChungWCKetterleW. Observation of bose–einstein condensation in a strong synthetic magnetic field. Nat Phys (2015) 11:859–64. 10.1038/nphys3421
6.
AvelineDCWilliamsJRElliottERDutenhofferCKelloggJRKohelJMet alObservation of bose–einstein condensates in an earth-orbiting research lab. Nature (2020) 582:193–7. 10.1038/s41586-020-2346-1
7.
GhosalADebP. Quantum walks over a square lattice. Phys Rev A (Coll Park) (2018) 98:032104. 10.1103/PhysRevA.98.032104
8.
WongTG. Isolated vertices in continuous-time quantum walks on dynamic graphs. Phys Rev A (2019) 100:062325. 10.1103/PhysRevA.100.062325
9.
SzigetiBEHomaGZimborásZBarankaiN. Short-time behavior of continuous-time quantum walks on graphs. Phys Rev A (Coll Park) (2019) 100:062320. 10.1103/PhysRevA.100.062320
10.
CaoWYangYLiDDongJZhouYShiW. Quantum state transfer on unsymmetrical graphs via discrete-time quantum walk. Mod Phys Lett A (2019) 34:1950317. 10.1142/S0217732319503176
11.
ZhanH. An infinite family of circulant graphs with perfect state transfer in discrete quantum walks. Quan Inf Process (2019) 18:369–26. 10.1007/s11128-019-2483-3
12.
FengYShiRShiJZhouJGuoY. Arbitrated quantum signature scheme with quantum walk-based teleportation. Quan Inf Process (2019) 18:154–21. 10.1007/s11128-019-2270-1
13.
RhodesMLWongTG. Quantum walk search on the complete bipartite graph. Phys Rev A (Coll Park) (2019) 99:032301. 10.1103/PhysRevA.99.032301
14.
QiangXWangYXueSGeRChenLLiuYet alImplementing graph-theoretic quantum algorithms on a silicon photonic quantum walk processor. Sci Adv (2021) 7:eabb8375. 10.1126/sciadv.abb8375
15.
WangYLouXFanZWangSHuangG. Verifiable multi-dimensional (t,n) threshold quantum secret sharing based on quantum walk. Int J Theor Phys (2022) 61:1–17. 10.1007/s10773-022-05009-w
16.
LouXWangSRenSZanHXuX. Quantum identity authentication scheme based on quantum walks on graphs with ibm quantum cloud platform. Int J Theor Phys (Dordr) (2022) 61:40–15. 10.1007/s10773-022-04986-2
17.
UnderwoodMSFederDL. Bose-hubbard model for universal quantum-walk-based computation. Phys Rev A (2012) 85:052314. 10.1103/PhysRevA.85.052314
18.
YeTian-YuGengMao-JieXuTian-JieChenYing. Efficient semiquantum key distribution based on single photons in both polarization and spatial-mode degrees of freedom. Quantum Information Processing, 2022, 21(4) 123.
19.
ShapiraDBihamOBrackenAHackettM. One-dimensional quantum walk with unitary noise. Phys Rev A (Coll Park) (2003) 68:062315. 10.1103/PhysRevA.68.062315
20.
BenedettiCBuscemiFBordonePParisMG. Dynamics of quantum correlations in colored-noise environments. Phys Rev A (Coll Park) (2013) 87:052328. 10.1103/PhysRevA.87.052328
21.
WangCWangCTangYRenS. Quantum walk in terms of quantum Bernoulli noise and quantum central limit theorem for quantum Bernoulli noise. Adv Math Phys (2018) 2018:1–9. 10.1155/2018/2507265
22.
AmbainisAWongTG. Correcting for potential barriers in quantum walk search. arXiv preprint (2015). 10.48550/arXiv.1505.02035
23.
YeTian-YuLiHong-KunHuJia-Li. Semi-quantum key distribution with single photons in both polarization and spatial-mode degrees of freedom. International Journal of Theoretical Physics, 2020, 59(9) 2807–2815.
24.
WangCYeX. Quantum walk in terms of quantum Bernoulli noises. Quan Inf Process (2016) 15:1897–908. 10.1007/s11128-016-1259-2
25.
DuYLuLLiY. A rout to protect quantum gates constructed via quantum walks from noises. Sci Rep (2018) 8:7117–1. 10.1038/s41598-018-25550-1
26.
BenedettiCRossiMAParisMG. Continuous-time quantum walks on dynamical percolation graphs. EPL (Europhysics Letters) (2019) 124:60001. 10.1209/0295-5075/124/60001
27.
CaiWMaYWangWZouCSunL. Bosonic quantum error correction codes in superconducting quantum circuits. Fundam Res (2021) 1:50–67. 10.1016/j.fmre.2020.12.006
28.
QiuTLiHXieM. Coherent generation and manipulation of stationary light pulses encoded in degrees of freedom of polarization and orbital angular momentum. Phys Rev A (Coll Park) (2019) 100:013844. 10.1103/PhysRevA.100.013844
29.
ChenHZhangYKalraMKLinFChenYLiaoPet alLow-dose ct with a residual encoder-decoder convolutional neural network. IEEE Trans Med Imaging (2017) 36:2524–35. 10.1109/TMI.2017.2715284
30.
XieTGrossmanJC. Crystal graph convolutional neural networks for an accurate and interpretable prediction of material properties. Phys Rev Lett (2018) 120:145301. 10.1103/PhysRevLett.120.145301
31.
LiJZhangHChenJZ. Structural prediction and inverse design by a strongly correlated neural network. Phys Rev Lett (2019) 123:108002. 10.1103/PhysRevLett.123.108002
32.
MaskaraNKubicaAJochym-O’ConnorT. Advantages of versatile neural-network decoding for topological codes. Phys Rev A (Coll Park) (2019) 99:052351. 10.1103/PhysRevA.99.052351
33.
VaronaSMartin-DelgadoMA. Determination of the semion code threshold using neural decoders. Phys Rev A (2020) 102:032411. 10.1103/PhysRevA.102.032411
Summary
Keywords
quantum error correction, multiparticle quantum walk, Bethe hypothesis, Bose–Hubbard model, threshold
Citation
Wang S-M, Qu Y-J, Wang H-W, Chen Z and Ma H-Y (2022) Multiparticle quantum walk–based error correction algorithm with two-lattice Bose–Hubbard model. Front. Phys. 10:1016009. doi: 10.3389/fphy.2022.1016009
Received
10 August 2022
Accepted
24 August 2022
Published
20 September 2022
Volume
10 - 2022
Edited by
Tianyu Ye, Zhejiang Gongshang University, China
Reviewed by
Hao Cao, Anhui Science and Technology University, China
Liyun Hu, Jiangxi Normal University, China
Xiaoping Lou, Hunan Normal University, China
Updates
Copyright
© 2022 Wang, Qu, Wang, Chen and Ma.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hong-Yang Ma, hongyang_ma@aliyun.com
This article was submitted to Quantum Engineering and Technology, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.