Quantum LDPC Codes Based on Cocyclic Block Matrices

Motivated by a family of binary cocyclic block matrices over GF(2), we proposed a construction method to gain the stabilizer of long-length quantum error-correction codes (QECCs). Stabilizer quantum codes (SQCs) can be obtained by the different rows of the yielded circulant permutation matrices; hence, the quantum codes have the virtue of a fast construction algorithm. The recursive relation of a block matrix is employed in the proposed approach, so that the generator matrix of quantum cocyclic codes with long length can be constructed easily. Furthermore, the obtained quantum codes have the low-density advantage of there being no 4-cycles in the Tanner graph.


Introduction
The theoretical and experimental research on virtual quantum secure communication has gradually matured.However, due to the constraints of single photon sources, single photon detectors and the high cost of quantum devices, its practical application still faces many technical challenges.Extending the error correction coding technology in classical communication to the field of quantum information, the breakthrough in quantum error correction coding technology will undoubtedly become an important link in the practical application of quantum secure communication in the future, although quantum errorcorrecting codes (QECCs) are of only nominal significance at present in both quantum key distribution (QKD) and quantum secure direct communication (QSDC) [1,2].Therefore, QECCs play theoretically a vital role in quantum information fields owing to their effects in fighting decoherence in quantum computing large-scale data [3][4][5][6][7][8][9][10].The appeared errors from quantum channel noise and possible eavesdropper's interference can be corrected by good QECCs.These codes are useful for building quantum fault-tolerant circuits.In general, an important class of quantum codes called stabilizer codes is firstly proposed by Gottesman [3].Currently, in quantum coding theory, the stabilizer formalism of QECCs and its special form called Calderbank-Shor-Steane (CSS) are its main construction methods [4][5][6].In quantum code construction theory, how to design the generator matrix resort to block matrix is the key.However, the error correction ability of short quantum codes is limited.The construction method of long-length quantum codes is becoming more important.In order to achieve high error correction ability, low-density parity check (LDPC) codes can achieve the maximum rate of reaching the Shannon limit with an error probability approaching zero when the codeword length tends to infinity.Since Gallager first proposed LDPC codes in the 1960s [11], two decades later, it caught the considerable interest of many researchers again [12][13][14][15].Subsequently, with the rapid development of quantum information, many research achievements in its quantum versions have been made in this field [16][17][18][19][20][21][22][23].However, compared to the in-depth study of classical LDPC codes, the corresponding quantum LDPC codes are still relatively less studied.One main obstacle is the difficulty of obtaining a more effective algorithm of iterative coding.With the application of large-scale data in quantum information fields, it becomes more important to encode massive data by making use of LDPC codes.
In quantum coding theory, the main problem for constructing low-density quantum code is how to obtain the generator for the stabilizer.Based on the current research results in this field, we will further study the encoding methods of quantum long-codes stabilizers.In this paper, a class of quantum stabilizer codes spirited with a family of cocyclic matrices is proposed [24,25].The involved cocyclic block matrices are based on the Jacket matrix, of which the main property is that its inverse matrix can be gained by block-wise inverse [26,27].According to the properties, its physical decoder circuit can be realized relatively simply.Meanwhile, the constructed quantum codes also have no 4-cycles in their Tanner graphs, so their parity check matrices have sparse characteristics.
The remainder of the paper is arranged as follows.In Section 2, we provide some necessary background and notations for QECCs.Then, in Sections 3 and 4, in terms of circulant permutation matrices, we investigate the construction of quantum codes generated from cocyclic block matrices which are applied to construct a generator matrix of quantum codes.Furthermore, its poverties of low density are also analyzed in this section.Finally, some concluding remarks are drawn in Section 4.

Preliminaries
In this section, we briefly follow here the relevant terms and notation for quantum error-correcting codes and its stabilizer, which are used to be constructed in the paper.

General Construction Methods of QECCs
In terms of the depolarizing channel, there are four basic operators performing on a single qubit, i.e., σ 0 , σ x , σ z and σ y = σ x σ z , that can form a set P whose eigenvalues are ±1 and ±i; here, σ x , σ z are the Pauli transformations.The stabilizer formalism of quantum codes can provide a compact description of codewords, especially for larger codes.Generally, a stabilizer S of quantum code C[[N, k, d]] is a nontrivial subgroup of the N-operator group here, any element in E * = {σ 0 , σ x , σ z , σ y } is a Pauli matrix.It shows that the set P N consists of all tensor products of Pauli matrices on N qubits again with multiplicative factors, which is a Hilbert codeword space of N binary codes.Therefore, for a stabilizer code C of encoding k logical qubits into N physical qubits, there are 2 N−k independent elements in 2 k -dimension code space C(S) associated with S of form To obtain the generator matrix of a stabilizer, two constraints should be met.One is that any pair of elements in stabilizer S commute, and the second condition implies that all elements of S are Hermitian, i.e., they have an eigenvalue of '±1'.It means that S is specified by N − k commuting and independent generators, {M i : 1 ≤ i ≤ N − k}.Correspondingly, if stabilizer quantum code is to be formed, N − k generators of stabilizer S should be firstly gained, which are described by matrix where According to the property of commuting generators, the elements of its row vector satisfy that the symplectic inner product is equal to zero.It means that an arbitrary operation γ i in the N-qubit depolarizing channel belongs to G N , which may be uniquely shown by where x i s , z i s ∈ {0, 1} for 1 ≤ s ≤ N, and τ is a factor among the eigenvalues.Without considering the factor τ, based on single-qibit operator γ i , we denote a 2N-dimensional concatenated vector as For any two concatenated row vectors γ i = (x i |z i ) and γ v = (x j |z j ) in generator matrix G N , the symplectic inner product is defined by Thereby, two operations γ i and γ j satisfy the commuting condition if The vector γ i in Equation ( 5) is generally used to seek the generators of the stabilizer.

Error of Quantum Error-Correction Code and Bound
All errors that occur in different quantum codes can be labeled with string vector over field F 2 when the states are transmitted through quantum channels.Single qubit error correction for quantum codes includes three classes, i.e, flip error, phase error or phase-flip error, which can be described with three corresponding Pauli operators: σ x , σ z and σ y , respectively.Assuming a single error qubit described as flip error X(a) and phase error Z(a) for a ∈ F 2 , a state |x ∈ C acted by error is denoted as Correspondingly, the errors on N qubits can be denoted as vector e = σ X x σ Z z , for X = (x 1 , x 2 , . . ., x N ) and Z = (z 1 , z 2 , . . ., z N ) ∈ F N 2 .Hence, it acts on an N-qubit basis state |Q = (q 1 , q 2 , . . ., q N ) over F N Furthermore, the quantum Hamming bound of a general quantum code C[[N, k, d]] is required to satisfy the following inequality condition [28] t which may correct up to t = [(d − 1)/2] quantum error bits.

Quantum Codes Based on Cocyclic Quasi-Jacket Block Matrices
Motivated by the center weight Hadamard matrices [26], a family of orthogonal matrices called Jacket transform was proposed [27], which is widely employed in signal processing, encoding, mobile communication, etc. [29][30][31][32].Especially, Hadamard and DFT matrices belong to this kind of Jacket matrix.

Definitions and Terms
Based on the general definition of a square matrix, in this section, we presented a construction method of quantum codes over binary matrices.As the following, for obtaining the stabilizer generators of quantum codes over the complex number field with a bigger size, we consider a class of binary cocyclic Jacket matrices based on a circular permutation matrix to construct its concatenated generator matrix, i.e., here, G x and G z are generator matrices of the encoding bit and phase, respectively.Since a cocyclic approach was introduced in [24], some interesting codes come from useful cocyclic matrices, such as the Hadamard matrix that is the matrix representation of a cocycle [25].

Definition 1.
If set G is a finite group of order v and A is an Abelian group of order w, respectively, a mapping: ϕ: G × G → A is called a cocycle while it satisfies the following particular cocycle equation From the definition, it is easy to show that ϕ(g, According to the cocycle mapping, if the rows and columns of a square matrix can be indexed by the elements of G in terms of some distribution, the matrix is called a cocyclic matrix, which is described as where the map ϕ(g, h) denotes a position (g, h) of the entry in matrix M ϕ .If the cocycle ϕ is symmetric, M ϕ will be a symmetric matrix.Furthermore, on the basis of the cocylic matrices, some cocyclic codes are primitively derived from Hadamard matrices or from related block designs [33].In this case, a v × v Hadamard matrix H with entries {±1} such that HH T = vI v is used to construct codes with length v = 4 s for integer s ≥ 1.If 1s and −1 s are replaced by 1 s and 0 s, respectively, the Hadamard matrix becomes a binary one.
A main superiority of cocyclic codes is that the long-length codes may be easily gained with a circulant manner; hence, it may be applicable to the generalized classical case and further quantum codes.
where a is the normalized constant, we call matrix J N a block Jacket matrix.Obviously, the Hadamard matrix belongs to the defined Jacket matrix.Furthermore, if there is a permutation matrix T , and its inverse matrix can be described as J −1 N = T • [a ij ] N /a, matrix J N is a quasi-Jacket block matrix [34].This class of matrices may be generalized by a simple binary transformation.Matrix decomposition with a Kronecker product of identity matrices can successively design a large length of codes that resorts to lower-order coefficient matrices.Meanwhile, this factorization decomposition generates a family of block circular sparse Jacket matrices which have good performance of LDPC codes.In the applications, the derived circulant permutation matrices with a Jacket pattern could lead to a very simple encoding algorithm.

Cocyclic Matrix Derived from Jacket Block Matrices
On the basis of the above definitions, we shall apply the index mapping to construct generator matrices G x and G z based on cocyclic Jacket matrices with order N = p m .Firstly, we form an Abelian group by factoring a quasi-Jacket block matrix, of which elements will be used to construct a stabilizer.
are two cocyclic Jacket matrices, their Kronecker product A u ⊗ B v will be also a cocyclic Jacket matrix with size uv [35].
Proof.Because A u and B v are Jacket matrices, in terms of the definition of Jacket matrix, it has where c a , c b are two constants.As a result, their Kronecker product is i.e., A u ⊗ B v is a Jacket matrix.
On the other hand, suppose that the row and column index orders of the two matrices as where s = r, c, and g a s1 (or g b s1 ) and g τ c1 (or g b c1 ) refer to the first row and column indices of matrix A u (or B v ), respectively.Hence, the row and column index orders of matrix Correspondingly, define Because both A u and B v are cocyclic matrices, ϕ A and ϕ B satisfy the condition shown in Equation (12).As a result, it is verified that ϕ AB also meets it with the index orders shown in Equation (18).It means that A u ⊗ B v also is a cocyclic matrix.
More generally, by making use of the mathematical methods of induction, it is easily verified that a large-length matrix with order N = p 1 p 2 • • • p n is still a Jacket matrix and cocylic matrix, while Jacket matrix with a Kronecker product also achieves the definitions of Jacket matrix and cocylic matrix.
In the following, we further illustrate that the generator matrix may be derived from a recursive cocyclic matrix.
Theorem 2. A recursive p-order cocyclic matrix J p m for a integer m can be composed of a family of cocyclic matrices where for that, J p is a cocyclic matrix and I p is an identity matrix.
Proof.We use the mathematical induction method to prove the result.In fact, while n = 1, the conclusion is obviously valid, i.e., J p = I p 0 ⊗ J p ⊗ I p 0 .Since the hypothesis is true for integer n = m, we consider that it holds for n = m + 1. Namely, the recursive equation holds.According to the property of the Kronecker product, it obeys the following formula for any square matrices A, B, C and D. Hence, It means that Equation ( 24) is equal to Equation ( 26), i.e., the case of n = m + 1 also holds.
In terms of the above conclusion, based on its decomposition factors, we can form a factor group Resorting to the factor group A, we consider constructing the stabilizer with two isomorphic classes of groups.
Lemma 1.Take factor group A in Equation ( 23) and its isomorphic factor group B = {B 1 N , B 2 N , • • • , B m N } with matrix J p .In terms of the property of the Kronecker product, if matrices J p and J p satisfy the commuting condition, i.e., their symplectic inner product is equal to zero, it is easy to check that any pair of factors A i N and B i N (1 ≤ i ≤ m) in the two groups also satisfies the commuting condition (7).

It means that any pair of A i
N and B i N can be taken as the generator matrices G x and G z of quantum codes, respectively, when two matrices J p and J p commute.Resorting to the binary matrixes, we consider building the generator matrix of quantum stabilizer codes.On the one hand, we define a class of fundamental matrices α and β derived from binary Hadamard matrix H as follows, respectively: It is easy to check that the introduced matrices have the following properties Furthermore, define a permutation matrix with 2q order as Based on the above fact, with simple computation, we can obtain a 4-order matrix and its inverse matrix which is a quasi-Jacket block matrix in terms of the definition.
Similarly, we may also obtain the 16-order quasi-Jacket block matrix as and its inverse matrix According to its definition, the matrix obviously is a quasi-Jacket block matrix.
Generally, by resorting to the recursive relationships based on J 4 defined by (30), we can obtain the quasi-Jacket block matrices with larger size 4 m by using the Kronecker product, i.e., which may be used to generate the generator matrices G x (or G z ).
On the other hand, to obtain the matrix G z (or G x ), we consider another class of matrix called the permutation matrix for an even number p ≥ 4. According to the definition, it is obviously a finite cyclic group P = {I, P, . . ., P p−1 } formed by a family of matrices P i , 0 ≤ i ≤ p − 1 is an Abelian group based on the matrix-multiplication operation, and I = P 0 is its unit element.Here, given that the number p is an even number and h is an odd number, respectively, a sequence of matrices whose exponents are odd can form a subset P o = {P, P 3 , . . ., P p−1 } of Abelian group P.
Based on the generated two classes of matrices, i.e., the cocyclic Jacket matrix and permutation matrix, we consider obtaining the stabilizer for quantum codes.Let G x with a size of N = 4 m be obtained by the cocyclic Jacket matrix described before.In the following, we consider how to construct another matrix G z .Firstly, we take two 4-order basic permutation matrices P 4 in Equation (34) with parameters h = 1, 3 respectively, i.e., {ω 2 , ω 4 }, such that Then, it provides that P N = P 4 ⊗ P N−1 , and it takes G z = P 4 ⊗ P N/4 .As a result, the concatenated generator matrix G is built.
In the following, we use the mathematical induction to deduce that two matrices G x and G z in the obtained concatenated generator matrix G satisfy self-orthogonal and commuting conditions, i.e., the obtained generator matrix should meet the conditions of self-orthogonal and commuting, so that it may be used to construct a stabilizer of quantum codes.In fact, it is obvious that matrix G z satisfies the self-orthogonal, so we firstly consider that the case of m = 1 for G x meets the self-orthogonal constrain.Namely, in the following matrix any two rows in G x and G z are self-orthogonal and commuting.Assuming it meets the case of m − 1, we consider the result for m.In fact, if any two rows in concatenated block submatrix [α β] meet the conditions, it can also do so in [β α].Without a loss of generality, here, we only consider the rows in submatrix [α β].By making use of the Kronecker product of any two row vectors in G x 4 , it is easy to check that they meet the orthogonal condition.And then, the symplectic inner product in Equation ( 6) is proved as following.Define complement operators 1 = 0 and 0 = 1, according to number theory; then, it is possible to easily check that matrices P = ( pij ) with parameters h = 1, 3 commute.Namely, any two row vectors in concatenated matrices ( P|P) satisfy the commuting condition in Equation (7).In fact, it can be checked that the elements in the Abelian group P commute when their cycle indexes both are odd or even numbers.The matrices satisfy isomorphic J ∼ = P when its corresponding orthogonal matrix also satisfies isomorphic J ∼ = P.The complement of G x is an element in the circulant permutation group P 0 ; hence, two pairs of vectors x 1 , x 2 and z 1 , z 2 all belong to this group.The product of any two vectors in G x and G z is equal to 0 or 1 simultaneously; i.e., it means that two pairs of vectors of the inner product with module 2 is equal to 0.
Similarly, we consider the generator matrix based on the block matrix J 16 and 16-order permutation matrix stemmed from P 4 .
Generally, the N-order generator matrix for N = 4 m , m ≥ 2 can be generated by On the other hand, from the factorization of a large number, 4 m , in Equation ( 22), we can obtain a group A with Formula (23) whose element 4 m in the group can be taken as the generator G x .Similarly, we can form another group B based on P 4 , and take the element B i 4 m (same structure to A i 4 m ) as the generator G z .As a result, a large-length generator matrix with size 4 m by 4 m is gained.

Quasi-Jacket Block Matrices with Size 4p m by 4p m
Furthermore, to obtain the more general length of code, we define a matrix It is easily checked that Furthermore, in terms of the basic properties, we consider matrix and its inverse matrix It is obvious that the constructed matrix Q is also a quasi-Jacket matrix which can be seen as a 4p × 4p binary cocyclic matrix.It means that a quantum code with a length of 4p may be obtained.Similarly, the recursive relationships may be used to obtain the large-length block matrices with a larger size 4p m by using the Kronecker product, i.e., such that the more general case of the generator matrix with a long size is achieved.In a similar manner, the obtained large-size matrix also satisfies the conditions of self-orthogonal and commuting like the process before proving the case of size 4 m .We take two matrices based on parameters h 1 , h 2 in P e = {ω, ω 3 , . . ., ω p−2 } or in P o = {ω 2 , ω 4 , . . ., ω p−1 } as G x and G z , respectively.It is evident that the obtained matrices commute each other.

Quantum LDPC Codes Based on Cocyclic Block Matrix
To encode k information qubits into a QC quantum code C with parameter [[N, k]], one should firstly gain two N × N cocyclic circulant permutation matrices G based on the previously constructed matrices G x and G z of satisfying the orthogonal condition.According to the recursive relationship, we take any N − k rows of the two generated matrices as the generators G x and G z .It means that the concatenated matrix G N in Equation (3) meets the construction condition of quantum code.
In classic coding theory, a class of very important linear codes called LDPC code is widely explored and used in practice.Its theoretical importance mainly embodies that its good performance of linear-time decoding can achieve Shannon capacity [14,16].Its quantum version, known as quantum LDPC (QLDPC) codes, plays a very important role in quantum information.However, compared to its classical counterparts, the achieved fruits of quantum LDPC codes are still far less.The main problem is the lack of efficient methods for obtaining iterative coding.Generally, the proportion of 1s in a binary matrix is called the matrix's destiny.Given that ι 0 and ι 1 are the numbers of 0s and 1s in a matrix, respectively, the destiny is According to the construction method described in Equation ( 41), it obvious that the total number of 0s and 1s in a 4p × 4p matrix is 16p 2 .Furthermore, each of the 12 matrices I, P −h and P h in G contain p 1s; hence, there are 12p 1s in this matrix: i.e, the 1s density in G is From the ratio, it means that the constructed cocylic matrix has a good performance of lowdensity property.On the other hand, there are N(log 2N − 1) additions and 2N(log 2N − 1) multiplications for the computation complexity in processing the obtained matrix with size N = 4 m .The check matrix H of LDPC codes is generally represented by a simpler and more intuitive Tanner graph with a one-to-one map.If there are equal degrees of all variable nodes in the Tanner graph, its corresponding LDPC code is called a regular code.Otherwise, it is called an irregular code.Denote by η and ρ the variable and check the nodes' degrees corresponding respectively to the minimum weights of the row and column of the check matrix.Assume that the girth g in the graph refers to the length of the shortest cycle [36,37].Then, the following bounds should be satisfied, i.e., for girth g = 4x + 2, and for girth g = 4x, where x is an integer, it can be seen from the above bounds that the code length will increase exponentially, while the cycle length and weights also are rising.
To obtain the good performance, having no 4-cycles in the Tanner graph of the quantum LDPC codes is required.A cycle is composed of a group of interconnected vertices in a Tanner graph, with one of these vertices serving as both the starting and ending points, and passing through each vertex only once.If the length of a cycle is called the number of lines it contains, the girth of a graph is defined as the minimum cycle length in the graph.The traditional cycle detection methods are generally based on block matrix sequences [15].When using this algebraic description to detect cycles in a given check matrix, there is redundancy due to the possibility that different block matrix sequences may form the same girth.To avoid redundancy caused by repeated traversal, we propose the following method for detection.Theorem 3. Based on the cocyclic block matrices in Equations ( 38) and (43), the constructed quantum codes have no 4−cycles in the Tanner graph.
Proof.According to the proposed construction method, matrix J 4 derived from matrices α and β obviously cannot form a cycle.If it considers a quasi-Jacket block matrix sequence with a cycle of girth 2i, (j 0 , l 0 ); (j 1 , l 1 ); . . .; (j i−1 , l i−1 ); (j 0 , l 0 ) as the loop is cyclic.Here, elements (j v , l v ) represent the block matrix of row j v and column l v for 0 ≤ v ≤ i.The formation of a cycle means that any block matrix sequence (j 0 , l i−1) ; (j i−1 , l i−2 ); (j i−2 , l i−3 ); (j v+1 , l v )(j 1 , l 0 ) is different from the original sequence.
Assume new coordinates (J v , L v ) to describe any block matrix in a block matrix sequence.For a cycle with a length of 2i, the block matrices in the sequence can be located at most in the i-th row.Without losing generality, the sequence of block matrices that may generate a 2i girth is classified according to the number of possible rows.Assume that i block matrices are located in the u(2 ≤ u ≤ i) row, and u rows are marked in order from top to bottom (1, 2, . . ., u).The new row coordinate ranges J v ∈ {1, 2, . . ., u} for each block matrix, as adjacent block matrices cannot be in the same row.Therefore, we can obtain J v ∈ {1, 2, . . ., m}\J v−1 , where 2 ≤ v < i.When v = i, the i-th block matrix is adjacent not only to the i-th block matrix but also to the first block matrix, so J i ∈ {1, 2, . . ., u}\{J 1 , J i−1 } can be obtained.For column coordinates, making L 1 = 1, "1" does not mean that the block matrix is in the first column, but rather that the column in which the block matrix is located is the first occurrence in the sequence.The subsequent block matrix is marked with new column coordinates L v in the order.We implement this by setting the parameter M v = max i=1:v L i , where M v represents the number of columns distributed in the first v block matrices.Then, for the v-th block matrix, L v ∈ {1, 2, . . ., M v−1 , M v−1 + 1}\L v−1 will indicate that the column can be in a different column from the first v − 1 block matrix.For the i−th block matrix located in row u(2 ≤ u ≤ i), the first block matrix will be J 1 = L 1 = 1, i.e., the new coordinate is (1,1) and M 1 = 1.For the v-th block matrix, we may obtain that L v ∈ {1, 2, . . ., M v−1 , M v−1 + 1}\L v−1 ; hence, the new coordinates of the current block matrix sequence have (u − 1) × (v − 1) possibilities.After traversing each possibility and determining all the branches of the i-th block matrix , we can obtain that the cycle of the constructed matrix is bigger than 4 when the parameter in construction method meets m ≥ 2.
According to the above description, after determining all the branches of the i-th block matrix, a tree diagram of the i−th layer will be obtained, thereby avoiding redundancy caused by repeated traversal.
In the classical coding field, we firstly compare the performance of the classical codes obtained from the cocyclic block matrices with the classical quasi-cyclic (QC) codes in Figure 1.In light of the character of the quasi-Jacket block matrix, the gained generator matrix of classical code has sparse property; hence, the parity check matrix of low weight can be obtained.In this figure, we take the different cocyclic matrices presented in this paper to analyze their bit error ratio (BER) and consider the coding performance of a single one among the generated pair classical codes.Here, given the different prime numbers p are taken n = 5, 7 with the construction method of the quasi-Jacket block matrix in this paper, their corresponding quantum code lengths are N = 4p m = 500, 1372 for m = 3, respectively.Assuming the binary phase shift keying (BPSK) modulation is over an additive white Gaussion noise (AWGN) channel, the bit error performance of the selected quantum codes is shown in this figure.As a result, the gained codes have the good properties of a sparse matrix.In decoding the performance of LDPC, short cycles will degrade the result under an iterative algorithm.The fewer short cycles of code makes its performance better; hence, the prime numbers p = 5, 7 more than 4 are taken to avoid 4−cycles in the Tanner graph.
It can be seen from the result that the constructed quantum LDPC codes have no apparent performance in the low signal-to-noise ratio (SNR) region, but the cocyclic codes are better than quasi-cyclic codes at high SNR.Furthermore, we also can see that the difference becomes more pronounced when the code length is larger.Therefore, the designed cocyclic codes may be better applied to construct large-length quantum codes.
Generally, three types of independent Pauli errors occur in a realistic depolarizing communication channel, i.e., bit-flip errors X, phase-flip errors Z and Y, and errors that are the combination of X and Z.Assuming that the total probability denotes as f , the probabilities of the three classes of errors are equally f /3 owing to their independence.Namely, the distributed marginal flip probability can be depicted as f m = 2 f /3.As a result, the errors of types X and Z in the obtained generator matrices G x and G z are separately corrected by reason of their isomorphism.
Furthermore, in the construction of quantum LDPC codes, the binary generator matrix of a quantum code is constructed by matrices H x and H z which can be seen as the parity-check matrices of two two classical codes C x and C z , respectively.The spaces of degenerate X-codewords C ⊥ x and Z-codewords C ⊥ z are derived from the rows of matrices H x and H z .The minimum distance d = min{d x , d z } of the quantum code is defined as the smaller of two-types distances d x and d z .According to the structure of the generator matrices, at least the minimum distances are bounded above by O(N 1/2 log N) as the code length N → ∞.To achieve the bigger length of cycle, the prime numbers p = 5, 7 more than 4 (for avoiding 4-cycles) are taken in the construction method.If the input is taken in the model of the Gaussian white noise channel (GWNC), we apply the quasi-Jacket block matrix to construct generator matrices.On the other hand, the quasi-cyclic codes with similar length are also constructed.
In the following, we compare the proposed quantum cocyclic codes to the quantum Reed-Solomon coding method with the previous work in ref. [38] under the 0.5 coding rate of constructed codes in Figure 2. The designed quantum codes with lengths N = 500, 1372, respectively, and similar length of CSS-type quantum Reed-Solomon (RS) codes are taken for comparison.Here, we take the generator in terms of the similar method from Equation (31) with size N = 1024, 4096 for parameter m = 5, 6 respectively.Furthermore, similar lengths of quantum Reed-Solomon codes are constructed according to the circulant permutation of two classical codes in its construction method.Here, we take the prime number p = 31, 61, so that the codes of size N = 930, 3660, respectively, are gained for comparison.Here, the proposed quantum codes with lengths N = 930, 3660 are generated by m = 5, 6 on a binary symmetric channel, and the conventional quantum LDPC codes are obtained with similar lengths with rates R q = 0.50, respectively.For simplicity, assuming the equal error rate of two kinds of errors X and Z, only the bit error rate of flip probability f m is shown here.
In the SNR range, it is shown that the taken two shorter lengths of quantum codes have better performance from 0.2 to 0.4 dB than the longer one.However, owing to the close cycles and their bit-flip error correction capability, the higher length of code shows its superiority with increasing length.

Conclusions
The coding method of massive data has important applications in the field of information transmission and big-data processing.With the rapid development of quantum computing, this class of coding construction becomes one of the research fields of quantum information.In this paper, based on a class of circulant permutation matrices, we presented quantum cocyclic stabilizer codes in terms of two long-length block matrices.Motivated by the advantage of convenient and low-complexity implementation, a special family of Jacket block matrices is applied to the generator construction of a stabilizer.The presented method can construct long-length quantum cocyclic codes based on the iterative coding process.If the involved parameters are designed appropriately, the gained quantum codes have no 4-cycles in their generators, and the LDPC codes have good performance.In terms of the analysis for marginal flip probability f m about our proposed codes, the good performance of the sparse matrices that appeared with a bit error rate is shown in Section IV.Therefore, the obtained quantum cocyclic stabilizer codes can be potentially applied to the future quantum information field of processing massive data.

Figure 1 .
Figure1.BER performance comparison of classical cocyclic code with N = 500, 1372 and classical quasi-cyclic codes with code rate R = 0.5.To achieve the bigger length of cycle, the prime numbers p = 5, 7 more than 4 (for avoiding 4-cycles) are taken in the construction method.If the input is taken in the model of the Gaussian white noise channel (GWNC), we apply the quasi-Jacket block matrix to construct generator matrices.On the other hand, the quasi-cyclic codes with similar length are also constructed.

Figure 2 .
Figure 2.Comparison of decoding error rate of the proposed quantum cocyclic codes (red) and conventional quantum Reed-Solomon LDPC codes (blue).Here, the proposed quantum codes with lengths N = 930, 3660 are generated by m = 5, 6 on a binary symmetric channel, and the conventional quantum LDPC codes are obtained with similar lengths with rates R q = 0.50, respectively.For simplicity, assuming the equal error rate of two kinds of errors X and Z, only the bit error rate of flip probability f m is shown here.