This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Necessary and sufficient condition for constructing
a single qudit insertion/deletion code
and its decoding algorithm

Taro Shibayama Kaijo Junior and Senior High School, 3-6-1 Okubo, Shinjuku, Tokyo, Japan, 169-0072. E-mail: shibayama@kaijo.ed.jp
Abstract

This paper shows that Knill-Laflamme condition, known as a necessary and sufficient condition for quantum error-correction, can be applied to quantum errors where the number of particles changes before and after the error. This fact shows that correctabilities of single deletion errors and single insertion errors are equivalent. By applying Knill-Laflamme condition, we generalize the previously known correction conditions for single insertion and deletion errors to necessary and sufficient level. By giving an example that satisfies this condition, we construct a new single qudit insertion/deletion code and explain its decoding algorithm.

1 Introduction

Quantum error-correcting codes play an important role in quantum information theory and have been actively studied since the 1990s[22, 3, 5]. The necessary and sufficient condition for quantum error-correction (KL condition) given by Knill-Laflamme in 1997[9] is extremely useful. Quantum errors in KL condition are mainly assumed to be those in which the number of particles does not change before and after the error, such as errors represented by unitary matrices. Recently, however, quantum errors where the number of particles changes, such as quantum insertion/deletion errors, have been attracting attention. Since 2020, several examples of quantum insertion/deletion error-correcting codes have been reported[13, 7, 17, 16, 19, 6, 12], and applications of quantum deletion codes have also been reported recently[1].

Classical insertion/deletion error-correcting codes were first given by Levenstein in 1966[11] and have been actively studied in recent years since then[8, 2, 15, 23, 4]. The most important property of classical codes is the equivalence between the error-correctability of insertion and deletion errors[11, 8]. Although it has been an open problem whether the equivalence holds in quantum theory, it has recently been shown that KL condition for deletions is equivalent to that for insertions of separable states[21]. In this paper, we first prove that KL condition can be used as a necessary and sufficient condition for error-correction even when the number of particles changes before and after the error.

By using KL conditions, the single deletion error-correction condition by Nakayama-Hagiwara[14] and the single insertion error-correction condition by Shibayama-Hagiwara[20] are improved to reach the necessary and sufficient level. These are also greatly generalized in that they give construction conditions not only for binary codes but also for non-binary codes. We also construct a new single qudit insertion/deletion code by giving an example that satisfies these conditions, and explain its decoding algorithm in detail.

This paper is organized as follows. Section 2 describes the variables and notations used in this paper. Section 3 explains that KL condition is also available for errors that change the number of particles. In Section 4, we define single deletion and single insertion errors and show that their correctabilities are equivalent. Section 5 gives necessary and sufficient conditions for the correction of single insertion/deletion errors. In Section 6, we give an example of the code and explain its decoding algorithm. Finally, Section 7 summarizes.

2 Preliminaries

The symbols and notations defined in this section will be used throughout this paper. Let nn be a positive integer and [n]:={1,2,,n}[n]:=\{1,2,\dots,n\}. We denote by l\mathbb{C}^{l} the ll-dimensional vector space over a complex field \mathbb{C} for an integer l2l\geq 2. Let |0,|1,|2,,|l1|0\rangle,|1\rangle,|2\rangle,\dots,|l-1\rangle be the standard orthogonal basis of l\mathbb{C}^{l} and :={0,1,2,,l1}\mathcal{L}:=\{0,1,2,\dots,l-1\}. Set |𝒙:=|x1|x2|xnln|\bm{x}\rangle:=|x_{1}\rangle\otimes|x_{2}\rangle\otimes\dots\otimes|x_{n}\rangle\in\mathbb{C}^{l\otimes n} for an nn-tuple 𝒙=x1x2xnn\bm{x}=x_{1}x_{2}\dots x_{n}\in\mathcal{L}^{n}. Here \otimes is the tensor product operation. Let 𝒙|:=|𝒙\langle\bm{x}|:=|\bm{x}\rangle^{\dagger} denote the conjugate transpose of |𝒙|\bm{x}\rangle. A positive semi-definite Hermitian matrix of trace 11 is called a density matrix. We denote by S(ln)S(\mathbb{C}^{l\otimes n}) the set of all density matrices of order lnl^{n}. An element of S(ln)S(\mathbb{C}^{l\otimes n}) is called an nn-qudit quantum state. In this paper, we also use a complex vector |ϕln|\phi\rangle\in\mathbb{C}^{l\otimes n} for representing a pure state |ϕϕ|S(ln)|\phi\rangle\langle\phi|\in S(\mathbb{C}^{l\otimes n}).

For vectors |0L,|1L,|2L,,|l1Lln|0_{L}\rangle,|1_{L}\rangle,|2_{L}\rangle,\dots,|l-1_{L}\rangle\in\mathbb{C}^{l\otimes n} of length 11 orthogonal to each other, the image of a linear map Enc{\rm Enc} that satisfies

|ψ\displaystyle|\psi\rangle :=iαi|il,\displaystyle:=\sum_{i\in\mathcal{L}}\alpha_{i}|i\rangle\in\mathbb{C}^{l},
|Ψ\displaystyle|\Psi\rangle :=iαi|iLln,\displaystyle:=\sum_{i\in\mathcal{L}}\alpha_{i}|i_{L}\rangle\in\mathbb{C}^{l\otimes n},

and Enc(|ψ)=|Ψ{\rm Enc}(|\psi\rangle)=|\Psi\rangle is called a quantum code, and the map Enc{\rm Enc} is called the encoder for the code. Here, vectors |0L,|1L,|2L,,|l1Lln|0_{L}\rangle,|1_{L}\rangle,|2_{L}\rangle,\dots,|l-1_{L}\rangle\in\mathbb{C}^{l\otimes n} are collectively called logical codewords.

Definition 2.1 (Quantum error).

For positive integers n,nn,n^{\prime}, take a set 𝒜={Aa}\mathcal{A}=\{A_{a}\} of lnl^{n^{\prime}}-by-lnl^{n} matrices such that aAaAa\sum_{a}A_{a}^{\dagger}A_{a} is the identity operator on S(ln)S(\mathbb{C}^{l\otimes n}). The map E:S(ln)S(ln)E:S(\mathbb{C}^{l\otimes n})\rightarrow S(\mathbb{C}^{l\otimes n^{\prime}}) defined by

E(ρ):=aAaρAa\displaystyle E(\rho):=\sum_{a}A_{a}\rho A_{a}^{\dagger}

is called a quantum error. Here, 𝒜\mathcal{A} is called a Kraus set for the error EE and an element of 𝒜\mathcal{A} is called a Kraus operator.

If n=nn=n^{\prime} and a linear map EE is a quantum channel, that is, completely positive and trace-preserving, then EE is a quantum error[10]. Definition 2.1 extends the ordinary definition of quantum error to include the case of nnn\neq n^{\prime}. This allows the case where the number of particles changes before and after an error, such as quantum insertion or quantum deletion, to be considered a quantum error.

If |ψS(l)|\psi\rangle\in S(\mathbb{C}^{l}) can be obtained by the operations allowed by quantum mechanics on E(Enc(|ψ))S(ln)E({\rm Enc}(|\psi\rangle))\in S(\mathbb{C}^{l\otimes n^{\prime}}), then the quantum code CC is correctable for the quantum error EE. In this paper, the quantum code CC is called 𝒜\mathcal{A}-correcting when |00|ψS(ln)|0\dots 0\rangle\otimes|\psi\rangle\in S(\mathbb{C}^{l\otimes n^{\prime}}) can be obtained by a finite number of measurements and unitary transformations. When CC is an 𝒜\mathcal{A}-correcting code, the original state |ψS(l)|\psi\rangle\in S(\mathbb{C}^{l}) can be obtained by deleting the 11st through the (n1)(n^{\prime}-1)th particles. The decoding in this paper uses the measurement described by the measurement operator ={Mk}\mathcal{M}=\{M_{k}\} that satisfies the completeness relation kMkMk=I\sum_{k}M_{k}^{\dagger}M_{k}=I. If we perform a measurement \mathcal{M} under a state ρS(ln)\rho\in S(\mathbb{C}^{l\otimes n^{\prime}}), the probability to get an outcome kk is p(k)=Tr(MkMkρ)p(k)={\rm Tr}(M_{k}^{\dagger}M_{k}\rho) and the after measurement state is MkρMk/p(k)M_{k}\rho M_{k}^{\dagger}/p(k). Here, we denote the sum of the diagonal elements of the square matrix MM by Tr(M){\rm Tr}(M) and the identity operator by II.

3 Correctable condition for quantum errors with non-square Kraus operators

This paper extends the recovery superoperator given in Ref.[9] to apply to non-square Kraus operators.

Definition 3.1 (Recovery superoperator).

The set ={Rr}\mathcal{R}=\{R_{r}\} of square matrices of order lnl^{n^{\prime}} such that the following two conditions are satisfied is called a recovery superoperator of (C,𝒜)(C,\mathcal{A}).

  • For any ii\in\mathcal{L}, RrR_{r}\in\mathcal{R}, Aa𝒜A_{a}\in\mathcal{A}, there exists λr,a\lambda_{r,a}\in\mathbb{C} such that

    RrAa|iL=λr,a|00i.\displaystyle R_{r}A_{a}|i_{L}\rangle=\lambda_{r,a}|0\dots 0i\rangle.
  • \mathcal{R} satisfies the completeness relation rRrRr=I\sum_{r}R_{r}^{\dagger}R_{r}=I.

In Theorem 3.2, the case n=nn=n^{\prime} is a well-known fact, and the 3rd condition in the theorem is known as KL condition[9]. Here we extend it and claim that it also holds for the case of non-square Kraus operators.

Theorem 3.2.

For a quantum code CC with logical codewords |0L,|1L,|2L,,|l1Lln|0_{L}\rangle,|1_{L}\rangle,|2_{L}\rangle,\dots,|l-1_{L}\rangle\in\mathbb{C}^{l\otimes n} and a quantum error E:S(ln)S(ln)E:S(\mathbb{C}^{l\otimes n})\rightarrow S(\mathbb{C}^{l\otimes n^{\prime}}) with the Kraus set 𝒜\mathcal{A}, the following three conditions are equivalent to each other:

  • 1.{\it 1.}

    The quantum code CC is 𝒜\mathcal{A}-correcting.

  • 2.{\it 2.}

    There exists a recovery superoperator of (C,𝒜)(C,\mathcal{A}).

  • 3.{\it 3.}

    (KL condition) For any Aa,Ab𝒜A_{a},A_{b}\in\mathcal{A}, and any i,ji,j\in\mathcal{L}, there exists μa,b\mu_{a,b}\in\mathbb{C} such that

    iL|AaAb|jL=μa,bδi,j,\displaystyle\langle i_{L}|A_{a}^{\dagger}A_{b}|j_{L}\rangle=\mu_{a,b}\delta_{i,j},

    where δi,j\delta_{i,j} denotes the Kronecker delta function.

In the case of non-square Kraus operators, it may not be possible to correct to the post-encoding state because the number of particles has changed, but it is possible to correct to the pre-encoding state. The concept of the proof of Theorem 3.2 is based on the one by Knill-Laflamme[9], but needs some modifications in this sense.

Proof.

This proof consists of three steps.

(Step 1) First, we prove that 12{\it 1}\Rightarrow{\it 2}.

Assume that the outcomes obtained in the decoding process are in turn k1,k2,,kmk_{1},k_{2},\dots,k_{m}, and the operators corresponding to them are in turn

U0,1,U1,2,U2,,m,Um,\displaystyle U_{0},\mathcal{M}^{1},U_{1},\mathcal{M}^{2},U_{2},\dots,\mathcal{M}^{m},U_{m},

where, U0,U1,U2,,UmU_{0},U_{1},U_{2},\dots,U_{m} are unitary operators and 1={Mk11},2={Mk22},,m={Mkmm}\mathcal{M}^{1}=\{M^{1}_{k_{1}}\},\mathcal{M}^{2}=\{M^{2}_{k_{2}}\},\dots,\mathcal{M}^{m}=\{M^{m}_{k_{m}}\} are measurement operators. Since the product of unitary matrices is a unitary matrix and the identity matrix is a unitary matrix, we may assume that the measurement and unitary transformation are performed alternately as described above. Note that the operator used may vary depending on the outcomes obtained each time. In this case, define R𝒌=UmMkmmU2Mk22U1Mk11U0R_{\bm{k}}=U_{m}M_{k_{m}}^{m}\dots U_{2}M_{k_{2}}^{2}U_{1}M_{k_{1}}^{1}U_{0} with 𝒌=k1k2km\bm{k}=k_{1}k_{2}\dots k_{m}, we can show that ={R𝒌}\mathcal{R}=\{R_{\bm{k}}\} is a recovery superoperator. Since

R𝒌(aAa(iαi|iL)(iαi¯iL|)Aa)R𝒌\displaystyle R_{\bm{k}}\left(\sum_{a}A_{a}\left(\sum_{i\in\mathcal{L}}\alpha_{i}|i_{L}\rangle\right)\left(\sum_{i\in\mathcal{L}}\overline{\alpha_{i}}\langle i_{L}|\right)A_{a}^{\dagger}\right)R_{\bm{k}}^{\dagger}
=j=1mpj(kj)(iαi|00i)(iαi¯00i|)\displaystyle=\prod_{j=1}^{m}p_{j}(k_{j})\left(\sum_{i\in\mathcal{L}}\alpha_{i}|0\dots 0i\rangle\right)\left(\sum_{i\in\mathcal{L}}\overline{\alpha_{i}}\langle 0\dots 0i|\right)

for any |ψ=iαi|il|\psi\rangle=\sum_{i\in\mathcal{L}}\alpha_{i}|i\rangle\in\mathbb{C}^{l}, we have

R𝒌Aa|iL=j=1mpj(kj)|00i,\displaystyle R_{\bm{k}}A_{a}|i_{L}\rangle=\sqrt{\prod_{j=1}^{m}p_{j}(k_{j})}|0\dots 0i\rangle,

where pj(kj)p_{j}(k_{j}) is the probability of obtaining the outcome kjk_{j} in the measurement j\mathcal{M}^{j} for 1jm1\leq j\leq m. We can also check that 𝒌R𝒌R𝒌=I\sum_{\bm{k}}R_{\bm{k}}^{\dagger}R_{\bm{k}}=I by considering that MkMkMk=I\sum_{M_{k}\in\mathcal{M}}M_{k}^{\dagger}M_{k}=I for the same measurement \mathcal{M} and UU=IU^{\dagger}U=I for any unitary matrix UU.

(Step 2) Next, we prove that 23{\it 2}\Rightarrow{\it 3}.

For any Aa,Ab𝒜A_{a},A_{b}\in\mathcal{A}, and any i,ji,j\in\mathcal{L}, we obtain

iL|AaAb|jL\displaystyle\langle i_{L}|A_{a}^{\dagger}A_{b}|j_{L}\rangle =iL|Aa(rRrRr)Ab|jL\displaystyle=\langle i_{L}|A_{a}^{\dagger}\left({\sum_{r}R_{r}^{\dagger}R_{r}}\right)A_{b}|j_{L}\rangle
=riL|AaRrRrAb|jL\displaystyle=\sum_{r}\langle i_{L}|A_{a}^{\dagger}R_{r}^{\dagger}R_{r}A_{b}|j_{L}\rangle
=rλr,a¯λr,b00i|00j\displaystyle=\sum_{r}\overline{\lambda_{r,a}}\lambda_{r,b}\langle 0\dots 0i|0\dots 0j\rangle
=(rλr,a¯λr,b)δi,j.\displaystyle=\left(\sum_{r}\overline{\lambda_{r,a}}\lambda_{r,b}\right)\delta_{i,j}.

Thus, KL condition is satisfied.

(Step 3) Finally, we prove that 31{\it 3}\Rightarrow{\it 1}.

Fix ii\in\mathcal{L} and let VilnV^{i}\subset\mathbb{C}^{l^{n^{\prime}}} be the vector space spanned by Aa|iLA_{a}|i_{L}\rangle for all Kraus operators Aa𝒜A_{a}\in\mathcal{A}. We define the basis {|u1i,|u2i,,|udi}\{|u_{1}^{i}\rangle,|u_{2}^{i}\rangle,\dots,|u_{d}^{i}\rangle\} of ViV^{i} by applying the Gram-Schmidt orthonormalization as follows. We define

|u~1i\displaystyle|\tilde{u}_{1}^{i}\rangle :=A1|iL,\displaystyle:=A_{1}|i_{L}\rangle,
|u~ki\displaystyle|\tilde{u}_{k}^{i}\rangle :=Ak|iLp=1k1u~pi|Ak|iLu~pi|u~pi|u~pi\displaystyle:=A_{k}|i_{L}\rangle-\sum_{p=1}^{k-1}\frac{\langle\tilde{u}_{p}^{i}|A_{k}|i_{L}\rangle}{\langle\tilde{u}_{p}^{i}|\tilde{u}_{p}^{i}\rangle}|\tilde{u}_{p}^{i}\rangle

for k2k\geq 2 to obtain vectors |u~1i,|u~2i,,|u~diVi|\tilde{u}_{1}^{i}\rangle,|\tilde{u}_{2}^{i}\rangle,\dots,|\tilde{u}_{d}^{i}\rangle\in V^{i} that are orthogonal to each other. Note that even if the Kraus set 𝒜\mathcal{A} is an infinite set, the above inductive operation must finish because the dimension dd of ViV^{i} is finite. In addition, even if the Kraus set 𝒜\mathcal{A} is uncountable, the orthogonal vectors |u~1i,|u~2i,,|u~diVi|\tilde{u}_{1}^{i}\rangle,|\tilde{u}_{2}^{i}\rangle,\dots,|\tilde{u}_{d}^{i}\rangle\in V^{i} can be obtained by selecting the Kraus operators A1,A2,,AdA_{1},A_{2},\dots,A_{d} in any order such that A1|iL,A2|iL,,Ad|iLA_{1}|i_{L}\rangle,A_{2}|i_{L}\rangle,\dots,A_{d}|i_{L}\rangle are linearly independent. By adjusting the lengths of all the obtained orthogonal vectors |u~1i,|u~2i,,|u~diVi|\tilde{u}_{1}^{i}\rangle,|\tilde{u}_{2}^{i}\rangle,\dots,|\tilde{u}_{d}^{i}\rangle\in V^{i} to 11, the basis {|u1i,|u2i,,|udi}\{|u_{1}^{i}\rangle,|u_{2}^{i}\rangle,\dots,|u_{d}^{i}\rangle\} of ViV^{i} is obtained. Note here that dd is independent of ii.

From the definition of the basis {|u1i,|u2i,,|udi}\{|u_{1}^{i}\rangle,|u_{2}^{i}\rangle,\dots,|u_{d}^{i}\rangle\}, we can represent

|uki=ck,1A1|iL+ck,2A2|iL++ck,kAk|iL\displaystyle|u_{k}^{i}\rangle=c_{k,1}A_{1}|i_{L}\rangle+c_{k,2}A_{2}|i_{L}\rangle+\dots+c_{k,k}A_{k}|i_{L}\rangle (1)

for k[d]k\in[d]. Note that ck,1,ck,2,,ck,kc_{k,1},c_{k,2},\dots,c_{k,k}\in\mathbb{C} do not depend on ii. This fact can be proved by mathematical induction as follows. From u~1i|u~1i=iL|A1A1|iL=μ1,1\langle\tilde{u}_{1}^{i}|\tilde{u}_{1}^{i}\rangle=\langle i_{L}|A_{1}^{\dagger}A_{1}|i_{L}\rangle=\mu_{1,1} for any ii\in\mathcal{L},

|u1i=|u~1i|u~1i=A1|iLμ1,1\displaystyle|u_{1}^{i}\rangle=\frac{|\tilde{u}_{1}^{i}\rangle}{\||\tilde{u}_{1}^{i}\rangle\|}=\frac{A_{1}|i_{L}\rangle}{\sqrt{\mu_{1,1}}}

holds and μ1,1\mu_{1,1} does not depend on ii, thus it holds when k=1k=1. When kmk\leq m for an integer m[d1]m\in[d-1], assume that

|u~ki=ck,1A1|iL+ck,2A2|iL++ck,kAk|iL\displaystyle|\tilde{u}_{k}^{i}\rangle=c^{\prime}_{k,1}A_{1}|i_{L}\rangle+c^{\prime}_{k,2}A_{2}|i_{L}\rangle+\dots+c^{\prime}_{k,k}A_{k}|i_{L}\rangle

and ck,1,ck,2,,ck,kc^{\prime}_{k,1},c^{\prime}_{k,2},\dots,c^{\prime}_{k,k}\in\mathbb{C} are all independent of ii. Then,

u~ki|u~ki\displaystyle\langle\tilde{u}_{k}^{i}|\tilde{u}_{k}^{i}\rangle =p=1kq=1kck,p¯ck,qiL|ApAq|iL\displaystyle=\sum_{p=1}^{k}\sum_{q=1}^{k}\overline{c^{\prime}_{k,p}}c^{\prime}_{k,q}\langle i_{L}|A_{p}^{\dagger}A_{q}|i_{L}\rangle
=p=1kq=1kck,p¯ck,qμp,q\displaystyle=\sum_{p=1}^{k}\sum_{q=1}^{k}\overline{c^{\prime}_{k,p}}c^{\prime}_{k,q}\mu_{p,q}

and

u~ki|Am+1|iL\displaystyle\langle\tilde{u}_{k}^{i}|A_{m+1}|i_{L}\rangle =p=1kck,p¯iL|ApAm+1|iL\displaystyle=\sum_{p=1}^{k}\overline{c^{\prime}_{k,p}}\langle i_{L}|A_{p}^{\dagger}A_{m+1}|i_{L}\rangle
=p=1kck,p¯μp,m+1\displaystyle=\sum_{p=1}^{k}\overline{c^{\prime}_{k,p}}\mu_{p,m+1}

are both independent of ii. Therefore, every coefficient of

|u~m+1i\displaystyle|\tilde{u}_{m+1}^{i}\rangle =Am+1|iLp=1mu~pi|Am+1|iLu~pi|u~pi|u~pi\displaystyle=A_{m+1}|i_{L}\rangle-\sum_{p=1}^{m}\frac{\langle\tilde{u}_{p}^{i}|A_{m+1}|i_{L}\rangle}{\langle\tilde{u}_{p}^{i}|\tilde{u}_{p}^{i}\rangle}|\tilde{u}_{p}^{i}\rangle
=Am+1|iLp=1m(u~pi|Am+1|iLu~pi|u~piq=1pcp,qAq|iL)\displaystyle=A_{m+1}|i_{L}\rangle-\sum_{p=1}^{m}\left(\frac{\langle\tilde{u}_{p}^{i}|A_{m+1}|i_{L}\rangle}{\langle\tilde{u}_{p}^{i}|\tilde{u}_{p}^{i}\rangle}\sum_{q=1}^{p}c^{\prime}_{p,q}A_{q}|i_{L}\rangle\right)

does not depend on ii. Hence, all coefficients of |um+1i|u_{m+1}^{i}\rangle are also independent of ii, and it is shown that it holds for k=m+1k=m+1.

Define Mk:=i|ukiuki|M_{k}:=\sum_{i\in\mathcal{L}}|u_{k}^{i}\rangle\langle u_{k}^{i}| for k[d]k\in[d]. For the space VlnV\subset\mathbb{C}^{l\otimes n^{\prime}} with the basis {|ukik[d],i}\{|u_{k}^{i}\rangle\mid k\in[d],i\in\mathcal{L}\}, choose a basis {|ej}\{|e_{j}\rangle\} of its orthogonal complementary space VV^{\perp} and define M:=j|ejej|M_{\emptyset}:=\sum_{j}|e_{j}\rangle\langle e_{j}|. Then :={M}{Mkk[d]}\mathcal{M}:=\{M_{\emptyset}\}\cup\{M_{k}\mid k\in[d]\} is a set of measurement operators because the completeness relation

MM+k[d]MkMk=I\displaystyle M_{\emptyset}^{\dagger}M_{\emptyset}+\sum_{k\in[d]}M_{k}^{\dagger}M_{k}=I

is satisfied. Here, II is the identity matrix of order lnl^{n^{\prime}}.

The measurement \mathcal{M} is performed on the state after the quantum error EE of the encoded state |ΨS(ln)|\Psi\rangle\in S(\mathbb{C}^{l\otimes n}). Then, the probability to get the outcome k[d]k\in[d] is

p(k)\displaystyle p(k) =Tr(MkMk(aAa|ΨΨ|Aa))\displaystyle={\rm Tr}\left(M_{k}^{\dagger}M_{k}\left(\sum_{a}A_{a}|\Psi\rangle\langle\Psi|A_{a}^{\dagger}\right)\right)
=𝒙na𝒙|MkMkAa|ΨΨ|Aa|𝒙\displaystyle=\sum_{\bm{x}\in\mathcal{L}^{n^{\prime}}}\sum_{a}\langle\bm{x}|M_{k}^{\dagger}M_{k}A_{a}|\Psi\rangle\langle\Psi|A_{a}^{\dagger}|\bm{x}\rangle
=aΨ|Aa(𝒙n|𝒙𝒙|)MkMkAa|Ψ\displaystyle=\sum_{a}\langle\Psi|A_{a}^{\dagger}\left(\sum_{\bm{x}\in\mathcal{L}^{n^{\prime}}}|\bm{x}\rangle\langle\bm{x}|\right)M_{k}^{\dagger}M_{k}A_{a}|\Psi\rangle
=aΨ|AaMkMkAa|Ψ.\displaystyle=\sum_{a}\langle\Psi|A_{a}^{\dagger}M_{k}^{\dagger}M_{k}A_{a}|\Psi\rangle.

Since uki|Aa|iL\langle u_{k}^{i}|A_{a}|i_{L}\rangle does not depend on ii from Equation (1), if we set this value to βk,a\beta_{k,a}\in\mathbb{C}, then

uki|Aa|jL={βk,ai=j0ij\displaystyle\langle u_{k}^{i}|A_{a}|j_{L}\rangle=\begin{cases}\beta_{k,a}&i=j\\ 0&i\neq j\end{cases}

holds. Therefore, we obtain

MkAa|Ψ\displaystyle M_{k}A_{a}|\Psi\rangle =i|ukiuki|Aajαj|jL\displaystyle=\sum_{i\in\mathcal{L}}|u_{k}^{i}\rangle\langle u_{k}^{i}|A_{a}\sum_{j\in\mathcal{L}}\alpha_{j}|j_{L}\rangle
=ijαj|ukiuki|Aa|jL\displaystyle=\sum_{i\in\mathcal{L}}\sum_{j\in\mathcal{L}}\alpha_{j}|u_{k}^{i}\rangle\langle u_{k}^{i}|A_{a}|j_{L}\rangle
=iαiβk,a|uki\displaystyle=\sum_{i\in\mathcal{L}}\alpha_{i}\beta_{k,a}|u_{k}^{i}\rangle
=βk,aiαi|uki.\displaystyle=\beta_{k,a}\sum_{i\in\mathcal{L}}\alpha_{i}|u_{k}^{i}\rangle. (2)

Hence, we have

p(k)\displaystyle p(k) =aβk,a¯βk,a(iαi¯uki|)(iαi|uki)\displaystyle=\sum_{a}\overline{\beta_{k,a}}\beta_{k,a}\left(\sum_{i\in\mathcal{L}}\overline{\alpha_{i}}\langle u_{k}^{i}|\right)\left(\sum_{i\in\mathcal{L}}\alpha_{i}|u_{k}^{i}\rangle\right)
=a|βk,a|2.\displaystyle=\sum_{a}|\beta_{k,a}|^{2}. (3)

The state after the measurement \mathcal{M} when the outcome kk is obtained is iαi|ukiS(ln)\sum_{i\in\mathcal{L}}\alpha_{i}|u_{k}^{i}\rangle\in S(\mathbb{C}^{l\otimes n^{\prime}}), since

Mk(aAa|ΨΨ|Aa)Mk\displaystyle M_{k}\left(\sum_{a}A_{a}|\Psi\rangle\langle\Psi|A_{a}^{\dagger}\right)M_{k}^{\dagger} =aMkAa|ΨΨ|AaMk\displaystyle=\sum_{a}M_{k}A_{a}|\Psi\rangle\langle\Psi|A_{a}^{\dagger}M_{k}^{\dagger}
=aβk,aβk,a¯(iαi|uki)(iαi¯uki|)\displaystyle=\sum_{a}\beta_{k,a}\overline{\beta_{k,a}}\left(\sum_{i\in\mathcal{L}}\alpha_{i}|u_{k}^{i}\rangle\right)\left(\sum_{i\in\mathcal{L}}\overline{\alpha_{i}}\langle u_{k}^{i}|\right)
=p(k)(iαi|uki)(iαi¯uki|)\displaystyle=p(k)\left(\sum_{i\in\mathcal{L}}\alpha_{i}|u_{k}^{i}\rangle\right)\left(\sum_{i\in\mathcal{L}}\overline{\alpha_{i}}\langle u_{k}^{i}|\right)

from Equation (2). On the other hand, from the definition of MM_{\emptyset}, the probability of obtaining the outcome \emptyset is 0.

For each k[d]k\in[d], take one unitary matrix UkU_{k} such that Uk|uki=|00iU_{k}|u_{k}^{i}\rangle=|0\dots 0i\rangle for any ii\in\mathcal{L}. When the outcome kk is obtained, applying the unitary operator UkU_{k} to the state after the measurement \mathcal{M}, we obtain

Uk(iαi|uki)=iαi|00i=|00|ψ.\displaystyle U_{k}\left(\sum_{i\in\mathcal{L}}\alpha_{i}|u_{k}^{i}\rangle\right)=\sum_{i\in\mathcal{L}}\alpha_{i}|0\dots 0i\rangle=|0\dots 0\rangle\otimes|\psi\rangle.

Therefore, the code CC is 𝒜\mathcal{A}-correcting. The original state |ψ=iαi|iS(l)|\psi\rangle=\sum_{i\in\mathcal{L}}\alpha_{i}|i\rangle\in S(\mathbb{C}^{l}) can be obtained by deleting the 11st through the (n1)(n^{\prime}-1)th particles and error-correction is completed. ∎

4 Equivalence of quantum single deletion and single insertion error-correctabilities

As typical errors that change the number of particles, this section defines single deletion errors and single insertion errors using Kraus operators.

Let m0m\geq 0 be an integer and let |ϕl|\phi\rangle\in\mathbb{C}^{l} with ϕ|ϕ=1\langle\phi|\phi\rangle=1. For an integer p[m+1]p\in[m+1], we define a lml^{m}-by-lm+1l^{m+1} matrix Dp,|ϕmD_{p,|\phi\rangle}^{m} and a lm+1l^{m+1}-by-lml^{m} matrix Ip,|ϕmI_{p,|\phi\rangle}^{m} as

Dp,|ϕm\displaystyle D_{p,|\phi\rangle}^{m} :=𝕀l𝕀l(p1)timesϕ|𝕀l𝕀l(mp+1)times,\displaystyle:=\underbrace{\mathbb{I}_{l}\otimes\dots\otimes\mathbb{I}_{l}}_{(p-1)\,{\rm times}}\otimes\langle\phi|\otimes\underbrace{\mathbb{I}_{l}\otimes\dots\otimes\mathbb{I}_{l}}_{(m-p+1)\,{\rm times}},
Ip,|ϕm\displaystyle I_{p,|\phi\rangle}^{m} :=𝕀l𝕀l(p1)times|ϕ𝕀l𝕀l(mp+1)times.\displaystyle:=\underbrace{\mathbb{I}_{l}\otimes\dots\otimes\mathbb{I}_{l}}_{(p-1)\,{\rm times}}\otimes|\phi\rangle\otimes\underbrace{\mathbb{I}_{l}\otimes\dots\otimes\mathbb{I}_{l}}_{(m-p+1)\,{\rm times}}.

Here, 𝕀l\mathbb{I}_{l} denotes the identity matrix of order ll.

Definition 4.1 (Single deletion error).

For a quantum state ρS(ln)\rho\in S(\mathbb{C}^{l\otimes n}), we define a quantum single deletion error as a map Edel:S(ln)S(l(n1))E^{del}:S(\mathbb{C}^{l\otimes n})\rightarrow S(\mathbb{C}^{l\otimes(n-1)}) expressed as

Edel(ρ):=p[n](p(p)bDp,|bn1ρDp,|bn1),\displaystyle E^{del}(\rho):=\sum_{p\in[n]}\left(p^{-}(p)\sum_{b\in\mathcal{L}}D_{p,|b\rangle}^{n-1}\rho{D_{p,|b\rangle}^{n-1}}^{\dagger}\right),

with p[n]p(p)=1\sum_{p\in[n]}p^{-}(p)=1 for a non-negative-valued function pp^{-}. I.e., the Kraus set for the single deletion error EdelE^{del} is

𝒟1:={p(p)Dp,|bn1|p[n],b}.\displaystyle\mathcal{D}_{1}:=\left\{\sqrt{p^{-}(p)}D_{p,|b\rangle}^{n-1}\,\middle|\,p\in[n],b\in\mathcal{L}\right\}. (4)

Although quantum single deletion error is often defined as a partial trace at an unknown position, it can also be defined as in Definition 4.1 from the discussion in Ref. [21].

Definition 4.2 (Single insertion error).

Suppose that a single qudit σS(l)\sigma\in S(\mathbb{C}^{l}) is represented as σ=bpb|ϕbϕb|\sigma=\sum_{b\in\mathcal{L}}p_{b}|\phi_{b}\rangle\langle\phi_{b}| in the form of a spectral decomposition. For a quantum state ρS(ln)\rho\in S(\mathbb{C}^{l\otimes n}), we define a quantum single insertion error of σ\sigma as a map Eσins:S(ln)S(l(n+1))E^{ins}_{\sigma}:S(\mathbb{C}^{l\otimes n})\rightarrow S(\mathbb{C}^{l\otimes(n+1)}) expressed as

Eσins(ρ):=p[n+1](p+(p)bpbIp,U|bnρIp,U|bn),\displaystyle E^{ins}_{\sigma}(\rho):=\sum_{p\in[n+1]}\left(p^{+}(p)\sum_{b\in\mathcal{L}}p_{b}I_{p,U|b\rangle}^{n}\rho{I_{p,U|b\rangle}^{n}}^{\dagger}\right),

with p[n+1]p+(p)=1\sum_{p\in[n+1]}p^{+}(p)=1 for a non-negative-valued function p+p^{+}. Here, UU is a unitary matrix such that U|b=|ϕbU|b\rangle=|\phi_{b}\rangle for every bb\in\mathcal{L}. I.e., the Kraus set for the single insertion error EσinsE^{ins}_{\sigma} is

1:={p+(p)pbIp,U|bn|p[n+1],b}.\displaystyle\mathcal{I}_{1}:=\left\{\sqrt{p^{+}(p)p_{b}}I_{p,U|b\rangle}^{n}\,\middle|\,p\in[n+1],b\in\mathcal{L}\right\}. (5)

Quantum single insertion error is often defined as the insertion of a quantum state into an unknown position, but it is known that if the state before the error is pure, it can be expressed in a Kraus representation as in Definition 4.2[20]. Throughout this paper, the discussion is based on the assumption that the quantum state before the error is pure. This means that the codeword is pure, which is a weak assumption.

It is pointed out that the following theorem presented by the author in 2022 does not show the equivalence of error-correctability[20].

Theorem 4.3.

[Theorem III.4,[20]] KL condition for single deletion error is equivalent to that for single insertion error.

From Theorem 3.2, reported for the first time in this paper, the equivalence of error-correctability is also shown.

Theorem 4.4.

It is equivalent for quantum code to be 𝒟1\mathcal{D}_{1}-correcting and 1\mathcal{I}_{1}-correcting.

5 A necessary and sufficient condition for single qudit insertion/deletion error-correction

This section discusses the code CC defined by the following logical codewords. Take AinA_{i}\subset\mathcal{L}^{n} for every ii\in\mathcal{L} and define

|iL:=1|Ai|𝒂Ai|𝒂.\displaystyle|i_{L}\rangle:=\frac{1}{\sqrt{|A_{i}|}}\sum_{{\bm{a}}\in A_{i}}|\bm{a}\rangle.

A necessary and sufficient condition for the correction of single deletion errors in the code CC is as follows:

Theorem 5.1.

It is a necessary and sufficient condition for the code CC to be 𝒟1\mathcal{D}_{1}-correcting that the following two conditions are satisfied simultaneously.

  • (C1del: ratio condition): For each p1,p2[n]p_{1},p_{2}\in[n] and b1,b2b_{1},b_{2}\in\mathcal{L}, the following values are constant regardless of ii\in\mathcal{L}:

    |Δp1,b1(Ai)Δp2,b2(Ai)|/|Ai|.\displaystyle{|\Delta^{-}_{p_{1},b_{1}}(A_{i})\cap\Delta^{-}_{p_{2},b_{2}}(A_{i})|}\,/\,{|A_{i}|}.
  • (C2del: distance condition): For any p1,p2[n]p_{1},p_{2}\in[n] and any b1,b2,i,jb_{1},b_{2},i,j\in\mathcal{L}, if iji\neq j, then

    |Δp1,b1(Ai)Δp2,b2(Aj)|=0.\displaystyle|\Delta^{-}_{p_{1},b_{1}}(A_{i})\cap\Delta^{-}_{p_{2},b_{2}}(A_{j})|=0.

We refer to (C1del) and (C2del) collectively as (Cdel). Here,

Δp,b(A):={a1ap1ap+1ana1ap1bap+1anA}\displaystyle\Delta_{p,b}^{-}(A):=\{a_{1}\dots a_{p-1}a_{p+1}\dots a_{n}\mid a_{1}\dots a_{p-1}ba_{p+1}\dots a_{n}\in A\}

for a non-empty set AnA\subset\mathcal{L}^{n} and p[n]p\in[n], bb\in\mathcal{L}.

Proof.

For p[n]p\in[n] and b,ib,i\in\mathcal{L}, we have

Dp,|bn1|iL\displaystyle D_{p,|b\rangle}^{n-1}|i_{L}\rangle =1|Ai|𝒂AiDp,|bn1|𝒂=1|Ai|𝒂~Δp,b(Ai)|𝒂~.\displaystyle=\frac{1}{\sqrt{|A_{i}|}}\sum_{\bm{a}\in A_{i}}D_{p,|b\rangle}^{n-1}|\bm{a}\rangle=\frac{1}{\sqrt{|A_{i}|}}\sum_{\tilde{\bm{a}}\in\Delta^{-}_{p,b}(A_{i})}|\tilde{\bm{a}}\rangle.

Hence, for p1,p2[n]p_{1},p_{2}\in[n] and b1,b2,i,jb_{1},b_{2},i,j\in\mathcal{L}, we obtain

iL|Dp1,|b1n1Dp2,|b2n1|jL\displaystyle\langle i_{L}|{D_{p_{1},|b_{1}\rangle}^{n-1}}^{\dagger}D_{p_{2},|b_{2}\rangle}^{n-1}|j_{L}\rangle =|Δp1,b1(Ai)Δp2,b2(Aj)||Ai||Aj|.\displaystyle=\frac{|\Delta^{-}_{p_{1},b_{1}}(A_{i})\cap\Delta^{-}_{p_{2},b_{2}}(A_{j})|}{\sqrt{|A_{i}||A_{j}|}}.

Therefore, from Theorem 3.2, the condition (Cdel) is a necessary and sufficient condition for 𝒟1\mathcal{D}_{1}-correcting. ∎

Theorem 5.1 extends the code by Nakayama-Hagiwara[14] with l=2l=2 to the case where ll is any positive integer, and further generalizes it to be necessary and sufficient level of error-correction. Note that the decoder shown in Ref.[14] is consequently the same as the decoder shown in the proof of Theorem 3.2 in this paper.

A necessary and sufficient condition for the correction of single insertion errors in the code CC is as follows:

Theorem 5.2.

It is a necessary and sufficient condition for the code CC to be 1\mathcal{I}_{1}-correcting that the following two conditions are satisfied simultaneously.

  • (C1ins: ratio condition): For each p1,p2[n+1]p_{1},p_{2}\in[n+1] and b1,b2b_{1},b_{2}\in\mathcal{L}, the following values are constant regardless of ii\in\mathcal{L}:

    |Δp1,b1+(Ai)Δp2,b2+(Ai)|/|Ai|.\displaystyle|\Delta^{+}_{p_{1},b_{1}}(A_{i})\cap\Delta^{+}_{p_{2},b_{2}}(A_{i})|\,/\,|A_{i}|.
  • (C2ins: distance condition): For any p1,p2[n+1]p_{1},p_{2}\in[n+1] and any b1,b2,i,jb_{1},b_{2},i,j\in\mathcal{L}, if iji\neq j, then

    |Δp1,b1+(Ai)Δp2,b2+(Aj)|=0.\displaystyle|\Delta^{+}_{p_{1},b_{1}}(A_{i})\cap\Delta^{+}_{p_{2},b_{2}}(A_{j})|=0.

We refer to (C1ins) and (C2ins) collectively as (Cins). Here,

Δp,b+(A)\displaystyle\Delta_{p,b}^{+}(A)\coloneqq {a1ap1bapana1ap1apanA}\displaystyle\{a_{1}\dots a_{p-1}ba_{p}\dots a_{n}\mid a_{1}\dots a_{p-1}a_{p}\dots a_{n}\in A\}

for a non-empty set AnA\subset\mathcal{L}^{n} and p[n+1]p\in[n+1], bb\in\mathcal{L}.

Proof.

For p[n+1]p\in[n+1] and b,ib,i\in\mathcal{L}, we have

Ip,|bn|iL\displaystyle I_{p,|b\rangle}^{n}|i_{L}\rangle =1|Ai|𝒂AiIp,|bn|𝒂=1|Ai|𝒂~Δp,b+(Ai)|𝒂~.\displaystyle=\frac{1}{\sqrt{|A_{i}|}}\sum_{\bm{a}\in A_{i}}I_{p,|b\rangle}^{n}|\bm{a}\rangle=\frac{1}{\sqrt{|A_{i}|}}\sum_{\tilde{\bm{a}}\in\Delta^{+}_{p,b}(A_{i})}|\tilde{\bm{a}}\rangle.

Hence, for p1,p2[n+1]p_{1},p_{2}\in[n+1] and b1,b2,i,jb_{1},b_{2},i,j\in\mathcal{L}, by noting that there exists a unitary matrix Up1,p2U_{p_{1},p_{2}} such that Ip1,U|b1n=Up1,p2Ip1,|b1nI_{p_{1},U|b_{1}\rangle}^{n}=U_{p_{1},p_{2}}I_{p_{1},|b_{1}\rangle}^{n} and Ip2,U|b2n=Up1,p2Ip2,|b2nI_{p_{2},U|b_{2}\rangle}^{n}=U_{p_{1},p_{2}}I_{p_{2},|b_{2}\rangle}^{n}, we obtain

iL|Ip1,U|b1nIp2,U|b2n|jL\displaystyle\langle i_{L}|{I_{p_{1},U|b_{1}\rangle}^{n}}^{\dagger}I_{p_{2},U|b_{2}\rangle}^{n}|j_{L}\rangle =iL|Ip1,|b1nIp2,|b2n|jL\displaystyle=\langle i_{L}|{I_{p_{1},|b_{1}\rangle}^{n}}^{\dagger}I_{p_{2},|b_{2}\rangle}^{n}|j_{L}\rangle
=|Δp1,b1+(Ai)Δp2,b2+(Aj)||Ai||Aj|.\displaystyle=\frac{|\Delta^{+}_{p_{1},b_{1}}(A_{i})\cap\Delta^{+}_{p_{2},b_{2}}(A_{j})|}{\sqrt{|A_{i}||A_{j}|}}.

Therefore, from Theorem 3.2, the condition (Cins) is a necessary and sufficient condition for 1\mathcal{I}_{1}-correcting. ∎

Theorem 5.2 extends ll to any positive integer for codes with l=2l=2 in Ref.[20] and greatly generalizes the error-correction condition to necessary and sufficient level. Note that although the 4-qubit code satisfies the condition (Cins), the decoder given in the proof of Theorem 3.2 in this paper is not the one by Hagiwara[6], but corresponds to the one given in Ref.[18].

Theorem 5.3.

Conditions (Cdel) and (Cins) are equivalent.

Theorem 5.3 follows immediately from Theorem 5.1, Theorem 5.2 and Theorem 4.4, but can also be proved classically as follows:

Classical proof of Theorem 5.3.

For integers p1,p2[n]p_{1},p_{2}\in[n], assume that p1p2p_{1}\leq p_{2} without loss of generality. The sequence 𝒙Δp1,b1(A)Δp2,b2(A)\bm{x}\in\Delta_{p_{1},b_{1}}^{-}(A)\cap\Delta_{p_{2},b_{2}}^{-}(A) and the sequence 𝒙~Δp1,b1+(A)Δp2+1,b2+(A)\tilde{\bm{x}}\in\Delta_{p_{1},b_{1}}^{+}(A)\cap\Delta_{p_{2}+1,b_{2}}^{+}(A) correspond one-to-one, where

𝒙\displaystyle\bm{x} =x1xp11xp1xp21xp2+1xn1,\displaystyle=x_{1}\dots x_{p_{1}-1}x_{p_{1}}\dots x_{p_{2}-1}x_{p_{2}+1}\dots x_{n-1},
𝒙~\displaystyle\tilde{\bm{x}} =x1xp11b1xp1xp21b2xp2+1xn1.\displaystyle=x_{1}\dots x_{p_{1}-1}b_{1}x_{p_{1}}\dots x_{p_{2}-1}b_{2}x_{p_{2}+1}\dots x_{n-1}.

Therefore, we obtain

|Δp1,b1(A)Δp2,b2(A)|=|Δp1,b1+(A)Δp2+1,b2+(A)|.\displaystyle|\Delta_{p_{1},b_{1}}^{-}(A)\cap\Delta_{p_{2},b_{2}}^{-}(A)|=|\Delta_{p_{1},b_{1}}^{+}(A)\cap\Delta_{p_{2}+1,b_{2}}^{+}(A)|.

On the other hand, we have

|Δp,b1+(A)Δp,b2+(A)|={|A|b1=b20b1b2\displaystyle|\Delta_{p,b_{1}}^{+}(A)\cap\Delta_{p,b_{2}}^{+}(A)|=\begin{cases}|A|&b_{1}=b_{2}\\ 0&b_{1}\neq b_{2}\end{cases}

for any p[n+1]p\in[n+1]. From the above, it is shown that (C1del) and (C1ins) are equivalent.

From the equivalence of classical insertion codes and classical deletion codes, the equivalence of (C2del) and (C2ins) is obvious. ∎

6 Example of insertion/deletion qudit code and its decoding algorithm

By defining the sets A0,A1A_{0},A_{1}, and A2A_{2} as follows, we can construct an example of a 66-qudit single insertion/deletion error-correcting code with l=3l=3:

A0\displaystyle A_{0} ={001122,112200,220011},\displaystyle=\{001122,112200,220011\},
A1\displaystyle A_{1} ={002211,110022,221100},\displaystyle=\{002211,110022,221100\},
A2\displaystyle A_{2} ={001100,112211,220022}.\displaystyle=\{001100,112211,220022\}.

We can easily check that these sets satisfy the condition (Cdel) by Table 1.

Table 1: Δp,b(Ai)\Delta^{-}_{p,b}(A_{i}) for p[6]p\in[6] and b,i{0,1,2}b,i\in\{0,1,2\}
Δp,b(A0)\Delta^{-}_{p,b}(A_{0}) A0={001122,112200,220011}A_{0}=\{001122,112200,220011\}
b=0b=0 b=1b=1 b=2b=2
p=1p=1 {01122}\{01122\} {12200}\{12200\} {20011}\{20011\}
p=2p=2 {01122}\{01122\} {12200}\{12200\} {20011}\{20011\}
p=3p=3 {22011}\{22011\} {00122}\{00122\} {11200}\{11200\}
p=4p=4 {22011}\{22011\} {00122}\{00122\} {11200}\{11200\}
p=5p=5 {11220}\{11220\} {22001}\{22001\} {00112}\{00112\}
p=6p=6 {11220}\{11220\} {22001}\{22001\} {00112}\{00112\}
Δp,b(A1)\Delta^{-}_{p,b}(A_{1}) A1={002211,110022,221100}A_{1}=\{002211,110022,221100\}
b=0b=0 b=1b=1 b=2b=2
p=1p=1 {02211}\{02211\} {10022}\{10022\} {21100}\{21100\}
p=2p=2 {02211}\{02211\} {10022}\{10022\} {21100}\{21100\}
p=3p=3 {11022}\{11022\} {22100}\{22100\} {00211}\{00211\}
p=4p=4 {11022}\{11022\} {22100}\{22100\} {00211}\{00211\}
p=5p=5 {22110}\{22110\} {00221}\{00221\} {11002}\{11002\}
p=6p=6 {22110}\{22110\} {00221}\{00221\} {11002}\{11002\}
Δp,b(A2)\Delta^{-}_{p,b}(A_{2}) A2={001100,112211,220022}A_{2}=\{001100,112211,220022\}
b=0b=0 b=1b=1 b=2b=2
p=1p=1 {01100}\{01100\} {12211}\{12211\} {20022}\{20022\}
p=2p=2 {01100}\{01100\} {12211}\{12211\} {20022}\{20022\}
p=3p=3 {22022}\{22022\} {00100}\{00100\} {11211}\{11211\}
p=4p=4 {22022}\{22022\} {00100}\{00100\} {11211}\{11211\}
p=5p=5 {00110}\{00110\} {11221}\{11221\} {22002}\{22002\}
p=6p=6 {00110}\{00110\} {11221}\{11221\} {22002}\{22002\}

By Theorem 5.3, these sets also satisfy the condition (Cins). Therefore, the quantum code whose logical codewords are

|0L\displaystyle|0_{L}\rangle :=13(|001122+|112200+|220011),\displaystyle:=\frac{1}{\sqrt{3}}(|001122\rangle+|112200\rangle+|220011\rangle),
|1L\displaystyle|1_{L}\rangle :=13(|002211+|110022+|221100),\displaystyle:=\frac{1}{\sqrt{3}}(|002211\rangle+|110022\rangle+|221100\rangle),
|2L\displaystyle|2_{L}\rangle :=13(|001100+|112211+|220022),\displaystyle:=\frac{1}{\sqrt{3}}(|001100\rangle+|112211\rangle+|220022\rangle),

is 𝒟1\mathcal{D}_{1}-correcting and 1\mathcal{I}_{1}-correcting. Using the decoding algorithm presented under Theorem 3.2, we can correct any single deletion error or any single insertion error as follows.

6.1 Decoding algorithm for single deletion errors

Since n=6n=6, from Equation (4) there are 18 Kraus operators for the single deletion error, which are labeled newly as AaA_{a}, where a{1,2,,18}a\in\{1,2,\dots,18\}. The dimension of the vector space ViV^{i} spanned by Aa|iLA_{a}|i_{L}\rangle for all Kraus operators AaA_{a} is 99, and its basis {|u1i,|u2i,,|u9i}\{|u_{1}^{i}\rangle,|u_{2}^{i}\rangle,\dots,|u_{9}^{i}\rangle\} is defined by applying the Gram-Schmidt orthonormalization. For example, for i=0i=0, we have the following results. Compare with Table 1.

|u10=|01122,\displaystyle|u_{1}^{0}\rangle=|01122\rangle,
|u20=|22011,\displaystyle|u_{2}^{0}\rangle=|22011\rangle,
|u30=|11220,\displaystyle|u_{3}^{0}\rangle=|11220\rangle,
|u40=|12200,\displaystyle|u_{4}^{0}\rangle=|12200\rangle,
|u50=|00122,\displaystyle|u_{5}^{0}\rangle=|00122\rangle,
|u60=|22001,\displaystyle|u_{6}^{0}\rangle=|22001\rangle,
|u70=|20011,\displaystyle|u_{7}^{0}\rangle=|20011\rangle,
|u80=|11200,\displaystyle|u_{8}^{0}\rangle=|11200\rangle,
|u90=|00112.\displaystyle|u_{9}^{0}\rangle=|00112\rangle.

When the measurement \mathcal{M} is performed on the quantum state Edel(α0|0L+α1|1L+α2|2L)S(35)E^{del}(\alpha_{0}|0_{L}\rangle+\alpha_{1}|1_{L}\rangle+\alpha_{2}|2_{L}\rangle)\in S(\mathbb{C}^{3\otimes 5}), the probability of obtaining the outcome k[9]k\in[9] is

p(k)={13(p(1)+p(2))k=1,4,713(p(3)+p(4))k=2,5,813(p(5)+p(6))k=3,6,9\displaystyle p(k)=\begin{cases}\displaystyle\frac{1}{3}(p^{-}(1)+p^{-}(2))&k=1,4,7\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle\frac{1}{3}(p^{-}(3)+p^{-}(4))&k=2,5,8\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle\frac{1}{3}(p^{-}(5)+p^{-}(6))&k=3,6,9\end{cases}

from Equation (3). Then, the unitary operator UkU_{k} corresponding to each outcome k[9]k\in[9] is applied and the 11st through 44th particles are deleted. Thus, the original quantum state α0|0+α1|1+α2|2)S(3)\alpha_{0}|0\rangle+\alpha_{1}|1\rangle+\alpha_{2}|2\rangle)\in S(\mathbb{C}^{3}) can be obtained.

6.2 Decoding algorithm for single insertion errors

From Equation (5) there are 2121 Kraus operators for the single insertion error, which are labeled newly as AaA_{a}, where a{1,2,,21}a\in\{1,2,\dots,21\}. The dimension of the vector space ViV^{i} spanned by Aa|iLA_{a}|i_{L}\rangle for all Kraus operators AaA_{a} is 2121, and its basis {|u1i,|u2i,,|u21i}\{|u_{1}^{i}\rangle,|u_{2}^{i}\rangle,\dots,|u_{21}^{i}\rangle\} is defined by applying the Gram-Schmidt orthonormalization. We have

|u7j+1i\displaystyle|u_{7j+1}^{i}\rangle =\displaystyle= I1,|j6|iL,\displaystyle\!\!\!\!I_{1,|j\rangle}^{6}|i_{L}\rangle,
|u7j+2i\displaystyle|u_{7j+2}^{i}\rangle =\displaystyle= 98(13I1,|j6|iL+I2,|j6|iL),\displaystyle\!\!\!\!\sqrt{\frac{9}{8}}\bigg{(}-\frac{1}{3}I_{1,|j\rangle}^{6}|i_{L}\rangle+I_{2,|j\rangle}^{6}|i_{L}\rangle\bigg{)},
|u7j+3i\displaystyle|u_{7j+3}^{i}\rangle =\displaystyle= 65(14I1,|j6|iL14I2,|j6|iL+I3,|j6|iL),\displaystyle\!\!\!\!\sqrt{\frac{6}{5}}\bigg{(}-\frac{1}{4}I_{1,|j\rangle}^{6}|i_{L}\rangle-\frac{1}{4}I_{2,|j\rangle}^{6}|i_{L}\rangle+I_{3,|j\rangle}^{6}|i_{L}\rangle\bigg{)},
|u7j+4i\displaystyle|u_{7j+4}^{i}\rangle =\displaystyle= 1513(110I1,|j6|iL+110I2,|j6|iL25I3,|j6|iL+I4,|j6|iL),\displaystyle\!\!\!\!\sqrt{\frac{15}{13}}\bigg{(}\frac{1}{10}I_{1,|j\rangle}^{6}|i_{L}\rangle+\frac{1}{10}I_{2,|j\rangle}^{6}|i_{L}\rangle-\frac{2}{5}I_{3,|j\rangle}^{6}|i_{L}\rangle+I_{4,|j\rangle}^{6}|i_{L}\rangle\bigg{)},
|u7j+5i\displaystyle|u_{7j+5}^{i}\rangle =\displaystyle= 3932(113I1,|j6|iL+113I2,|j6|iL413I3,|j6|iL313I4,|j6|iL+I5,|j6|iL),\displaystyle\!\!\!\!\sqrt{\frac{39}{32}}\bigg{(}\frac{1}{13}I_{1,|j\rangle}^{6}|i_{L}\rangle+\frac{1}{13}I_{2,|j\rangle}^{6}|i_{L}\rangle-\frac{4}{13}I_{3,|j\rangle}^{6}|i_{L}\rangle-\frac{3}{13}I_{4,|j\rangle}^{6}|i_{L}\rangle+I_{5,|j\rangle}^{6}|i_{L}\rangle\bigg{)},
|u7j+6i\displaystyle|u_{7j+6}^{i}\rangle =\displaystyle= 9683(132I1,|j6|iL132I2,|j6|iL+18I3,|j6|iL+332I4,|j6|iL\displaystyle\!\!\!\!\sqrt{\frac{96}{83}}\bigg{(}-\frac{1}{32}I_{1,|j\rangle}^{6}|i_{L}\rangle-\frac{1}{32}I_{2,|j\rangle}^{6}|i_{L}\rangle+\frac{1}{8}I_{3,|j\rangle}^{6}|i_{L}\rangle+\frac{3}{32}I_{4,|j\rangle}^{6}|i_{L}\rangle
1332I5,|j6|iL+I6,|j6|iL),\displaystyle\!\!\!\!-\frac{13}{32}I_{5,|j\rangle}^{6}|i_{L}\rangle+I_{6,|j\rangle}^{6}|i_{L}\rangle\bigg{)},
|u7j+7i\displaystyle|u_{7j+7}^{i}\rangle =\displaystyle= 8368(283I1,|j6|iL283I2,|j6|iL+883I3,|j6|iL+683I4,|j6|iL\displaystyle\!\!\!\!\sqrt{\frac{83}{68}}\bigg{(}-\frac{2}{83}I_{1,|j\rangle}^{6}|i_{L}\rangle-\frac{2}{83}I_{2,|j\rangle}^{6}|i_{L}\rangle+\frac{8}{83}I_{3,|j\rangle}^{6}|i_{L}\rangle+\frac{6}{83}I_{4,|j\rangle}^{6}|i_{L}\rangle
2683I5,|j6|iL1983I6,|j6|iL+I7,|j6|iL),\displaystyle\!\!\!\!-\frac{26}{83}I_{5,|j\rangle}^{6}|i_{L}\rangle-\frac{19}{83}I_{6,|j\rangle}^{6}|i_{L}\rangle+I_{7,|j\rangle}^{6}|i_{L}\rangle\bigg{)},

for i,j{0,1,2}i,j\in\{0,1,2\}. Unlike the deletion case, since Aa|iLA_{a}|i_{L}\rangle are not all orthogonal to each other, the result is complicated.

When the measurement \mathcal{M} is performed on the quantum state Eins(α0|0L+α1|1L+α2|2L)S(37)E^{ins}(\alpha_{0}|0_{L}\rangle+\alpha_{1}|1_{L}\rangle+\alpha_{2}|2_{L}\rangle)\in S(\mathbb{C}^{3\otimes 7}), the probability of obtaining the outcome k[21]k\in[21] is

p(k)={pj(p+(1)+19p+(2)+19p+(3))k=7j+1pj(89p+(2)+118p+(3))k=7j+2pj(56p+(3)+215p+(4)+215p+(5))k=7j+3pj(1315p+(4)+365p+(5))k=7j+4pj(3239p+(5)+1396p+(6)+1396p+(7))k=7j+5pj(8396p+(6)+3617968p+(7))k=7j+6pj(6883p+(7))k=7j+7\displaystyle p(k)=\begin{cases}\displaystyle p_{j}\left(p^{+}(1)+\frac{1}{9}p^{+}(2)+\frac{1}{9}p^{+}(3)\right)&k=7j+1\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle p_{j}\left(\frac{8}{9}p^{+}(2)+\frac{1}{18}p^{+}(3)\right)&k=7j+2\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle p_{j}\left(\frac{5}{6}p^{+}(3)+\frac{2}{15}p^{+}(4)+\frac{2}{15}p^{+}(5)\right)&k=7j+3\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle p_{j}\left(\frac{13}{15}p^{+}(4)+\frac{3}{65}p^{+}(5)\right)&k=7j+4\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle p_{j}\left(\frac{32}{39}p^{+}(5)+\frac{13}{96}p^{+}(6)+\frac{13}{96}p^{+}(7)\right)&k=7j+5\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle p_{j}\left(\frac{83}{96}p^{+}(6)+\frac{361}{7968}p^{+}(7)\right)&k=7j+6\vskip 6.0pt plus 2.0pt minus 2.0pt\\ \displaystyle p_{j}\left(\frac{68}{83}p^{+}(7)\right)&k=7j+7\end{cases}

for j{0,1,2}j\in\{0,1,2\} from Equation (3). Then, the unitary operator UkU_{k} corresponding to each outcome k[21]k\in[21] is applied and the 11st through 66th particles are deleted. Thus, the original quantum state α0|0+α1|1+α2|2)S(3)\alpha_{0}|0\rangle+\alpha_{1}|1\rangle+\alpha_{2}|2\rangle)\in S(\mathbb{C}^{3}) can be obtained.

For example, if σ=12|00|+13|11|+16|22|S(3)\sigma=\frac{1}{2}|0\rangle\langle 0|+\frac{1}{3}|1\rangle\langle 1|+\frac{1}{6}|2\rangle\langle 2|\in S(\mathbb{C}^{3}) is inserted in the 44th position, that is, if we consider the single insertion error where p+(4)=1p^{+}(4)=1, p+(p)=0p^{+}(p)=0 for any p4p\neq 4, and p1=1/2p_{1}=1/2, p2=1/3p_{2}=1/3, p3=1/6p_{3}=1/6, the error-correction process is as shown in Figure 1.

Refer to caption
Figure 1: Error-correction process for the 6-qudit insertion code when σ=12|00|+13|11|+16|22|S(3)\sigma=\frac{1}{2}|0\rangle\langle 0|+\frac{1}{3}|1\rangle\langle 1|+\frac{1}{6}|2\rangle\langle 2|\in S(\mathbb{C}^{3}) is inserted in the 4th position

7 Conclusion

In this paper, we proved that KL condition [9] can be used as a necessary and sufficient condition for correcting quantum errors with changing number of particles. By using KL condition, we showed that the error-correctabilities of single deletion and single insertion are equivalent. Furthermore, the error-correctability conditions for single deletion given by Nakayama-Hagiwara[14] and for single insertion given by Shibayama-Hagiwara[20] were improved to necessary and sufficient level. We also constructed a new single qudit insertion/deletion code by giving an example that satisfies the conditions, and provided a decoding algorithm for the code.

References

  • [1] Yuta Aoki and Mitsuru Tada. A quantum secret sharing scheme using a 4-qubit code. IEICE Tech. Rep., IT2023-38, in Japanese, 123(338):45–50, 2024.
  • [2] Tilo Buschmann and Leonid V. Bystrykh. Levenshtein error-correcting barcodes for multiplexed dna sequencing. BMC Bioinformatics, 14(1), 2013.
  • [3] A. Robert Calderbank and Peter W. Shor. Good quantum error-correcting codes exist. Physical Review A, 54(2):1098–1105, 1996.
  • [4] Roni Con, Amir Shpilka, and Itzhak Tamo. Explicit and efficient constructions of linear codes against adversarial insertions and deletions. IEEE Transactions on Information Theory, 68(10):6516–6526, 2022.
  • [5] Daniel Gottesman. Stabilizer codes and quantum error correction. arXiv preprint arXiv:9705052, 1997.
  • [6] Manabu Hagiwara. The four qubits deletion code is the first quantum insertion code. IEICE Communications Express, 10(5):243–247, 2021.
  • [7] Manabu Hagiwara and Ayumu Nakayama. A four-qubits code that is a quantum deletion error-correcting code with the optimal length. In The Proceedings of 2020 IEEE International Symposium on Information Theory (ISIT2020), pages 1870–1874, 2020.
  • [8] Albertus S. J. Helberg and Hendrik C. Ferreira. On multiple insertion/deletion correcting codes. IEEE Transactions on Information Theory, 48(1):305–308, 2002.
  • [9] Emanuel Knill and Raymond Laflamme. Theory of quantum error-correcting codes. Physical Review A, 55(2):900–911, 1997.
  • [10] K. Kraus, A. Böhm, J.D. Dollard, and W.H. Wootters. States, Effects, and Operations: Fundamental Notions of Quantum Theory. Lecture Notes in Physics. Springer Berlin Heidelberg, 1983.
  • [11] Vladimir I Levenshtein. Binary codes capable of correcting deletions, insertions, and reversals. Soviet physics doklady, 10:707–710, 1966.
  • [12] Ken Nakamura and Takayuki Nozaki. Decoding algorithm correcting single-insertion plus single-deletion for non-binary quantum codes. arXiv preprint arXiv:2409.10924, 2024.
  • [13] Ayumu Nakayama and Manabu Hagiwara. The first quantum error-correcting code for single deletion errors. IEICE Communications Express, 9(4):100–104, 2020.
  • [14] Ayumu Nakayama and Manabu Hagiwara. Single quantum deletion error-correcting codes. In The Proceedings of 2020 International Symposium on Information Theory and Its Applications (ISITA2020), pages 329–333, 2020.
  • [15] Takayuki Nozaki. Bounded single insertion/deletion correcting codes. In The Proceedings of 2019 IEEE International Symposium on Information Theory (ISIT2019), pages 2379–2383, 2019.
  • [16] Yingkai Ouyang. Permutation-invariant quantum coding for quantum deletion channels. In The Proceedings of 2021 IEEE International Symposium on Information Theory (ISIT2021), pages 1499–1503, 2021.
  • [17] Taro Shibayama. Construction of single quantum deletion codes via combinatorial conditions and adjacency matrices. Quantum Information Processing, 20(9), 2021.
  • [18] Taro Shibayama. A generalization of the four qubits single insertion error-correcting code. IEICE Communications Express, 13(4):126–129, 2024.
  • [19] Taro Shibayama and Manabu Hagiwara. Permutation-invariant quantum codes for deletion errors. In The Proceedings of 2021 IEEE International Symposium on Information Theory (ISIT2021), pages 1493–1498, 2021.
  • [20] Taro Shibayama and Manabu Hagiwara. Equivalence of quantum single insertion and single deletion error-correctabilities, and construction of codes and decoders. In The Proceedings of 2022 IEEE International Symposium on Information Theory (ISIT2022), pages 2957–2962, 2022.
  • [21] Taro Shibayama and Yingkai Ouyang. The equivalence between correctability of deletions and insertions of separable states in quantum codes. In The Proceedings of 2021 IEEE Information Theory Workshop (ITW2021), pages 1–6, 2021.
  • [22] Peter W. Shor. Scheme for reducing decoherence in quantum computer memory. Physical Review A, 52:R2493–R2496, 1995.
  • [23] Evagoras Stylianou, Lorenz Welter, Rawad Bitar, Antonia Wachter-Zeh, and Eitan Yaakobi. Equivalence of insertion/deletion correcting codes for d-dimensional arrays. In The Proceedings of 2022 IEEE International Symposium on Information Theory (ISIT2022), pages 814–819, 2022.