Quantum computation, in a more strict fundamental understanding, is the movement of a point in a 2n-dimensional complex space, where the initial and final positions of the point correspond to the input and output of the computation.

Therefore, the process of quantum computation is nothing more than a linear transformation of a $2^n$ - dimensional vector from a state $\mid \psi >$ to $\mid \varphi >$.

Thus, the coordinates of vector $\mid \psi >$ are the initial conditions, and the coordinates of vector $\mid \varphi >$ are the result of the computation - the output.

An ideal quantum computation is then one that performs a direct linear transformation from $\mid \psi >$ to $\mid \varphi >$ in a single step - a singular transform, a single computation execution.

Ideal here means that there is no more efficient, in terms of computational power, computational expense of the quantum computer’s work.

To make ideal quantum computation technically possible, it is necessary to find a group of unitary matrices, the elements of which form a universal set of quantum gates for the direct transformation from $\mid \psi >$ to $\mid \varphi >$. As is known, quantum computations use quantum bits or qubits instead of classical bits, which have two basis states $\mid 0>$ и $\mid 1>$. All other states of a qubit are defined as a linear combination of basis states with complex coefficients, that is where ${\lambda }_1\in C,{\lambda }_2\in C$ and $\mid {\lambda }_1{\mid }^2+\mid {\lambda }_2{\mid }^2=1,C$ - being the space of complex numbers.

The basis states of a qubit $\mid 0>$ and $\mid 1>$ are also denoted using vectors ${(1,0)}^T$ and ${(0,1)}^T$ respectively, where the index T denotes the transposition sign.

A two-qubit quantum system consists of two qubits and has four basis states $\mid 00>={(1,0,0,0)}^T,\mid 01>={(0,1,0,0)}^T,\mid 10>={(0,0,1,0)}^T$ and $\mid 11>={(0,0,0,1)}^T$. Then, an arbitrary state of a two-qubit quantum system can be written as [7,8] where ${\lambda }_i\in C,i=1,2,3,4,$ and $\mid {\lambda }_1{\mid }^2+\mid {\lambda }_2{\mid }^2+\mid {\lambda }_3{\mid }^2+\mid {\lambda }_4{\mid }^2=1.$ A similar record for an arbitrary state of a two-qubit system looks like [7,8] i.e., any quantum state of a two-qubit quantum system is uniquely determined by the complex amplitudes ${\lambda }_1,{\lambda }_2,{\lambda }_3,{\lambda }_4$.

Similarly, n-qubit quantum systems, consisting of n qubits, are defined. Such a system would have $2^n$ basis states and any of its states are determined as a linear combination of basis states with complex amplitudes. In this work, we only consider two-qubit quantum systems, although all results can in principle be generalized for an n-qubit quantum system. However, this would require the development of the basics of $2^n$-dimensional mathematics with commutative multiplication. Here, we rely on four-dimensional mathematics, the foundations of which were laid by the great Kazakh mathematician Abenov M.M. [1].

Quantum computations consist in the sequential application of unitary U operators to the quantum state of the quantum system, which are called gates. Unitary operators or gates applied to an n-qubit quantum system are represented in the form of a matrix of size $2^n\times 2^n$. For example, unitary operators for a two-qubit system have the form of a $4\times 4$ matrix U:

Note that if a certain matrix (4) is a two-qubit gate, then the Hermitian conjugate matrix to it is also a gate.

The main two-qubit gates, or binary operators, are SWAP, CNOT, CZ, represented by matrices [4,8]

Note that two-qubit gates can be applied to any n-qubit quantum system ($n\ge 2$), with any two qubits from the n-qubit system selected to which the two-qubit gate is applied, and the identity operator, which is defined by the unit matrix of the required dimension, is applied to the other qubits.

In addition to the gates indicated in (6), one can construct many different two-qubit gates, but a general description of all possible two-qubit gates is absent. In the work [3], with the participation of one of the authors, an Abelian group of commutating two-qubit gates was constructed. However, the constructed group of gates is not complete, as none of the gates (6) belongs to the mentioned group. This work is a continuation of [3], and here we have constructed all other commutative groups of two-qubit gates, as well as defined the correspondence between quantum states and gates of two-qubit quantum systems, introduced the concept of a unitary state and shown the possibility of transitioning in one step from any unitary state to any basis state, and also from any unitary state to any other unitary state. The obtained results open up new possibilities for constructing quantum algorithms not only for two-qubit quantum systems.

In the study [2], all spaces of four-dimensional numbers $R^4$ with commutative multiplication are examined. There are, in total, six such spaces, denoted as $M_2,M_3,M_4,M_5,M_6,M_7$. To each four-dimensional number $Z=(z_1,z_2,z_3,$$z_4)\in R^4$ from any of these spaces, a certain $4\times 4$ matrix M, is associated, with its elements being the components of the four-dimensional number Z, and this mapping is bijective. Moreover, this bijection is a homomorphism with respect to the multiplication operation of four-dimensional numbers, meaning the group of matrices forms a commutative group with identity. The results obtained in work [2] can be transferred to the case of four-dimensional numbers $Z\in C^4$. Retaining the same designations for the spaces $M_j(j=2,3,...,7)$ describe the necessary properties of these spaces for the case of complex-valued four-dimensional numbers. As we will see below, in each of the spaces $M_2,...,M_7$, there exist two groups of matrices corresponding to one operation of commutative multiplication of four-dimensional numbers [2]. The operations of addition $X+Y$ and subtraction $X-Y$ of four-dimensional numbers $X\in C^4$ and $Y\in C^4$ are defined as component-wise addition and subtraction. The multiplication operation of four-dimensional numbers can be defined in various ways, among which we are only interested in commutative multiplication. All ways of defining commutative multiplication are given in [2], where for each method, the corres-ponding space of four-dimensional numbers $M_j(j=2,3,...,7)$ is defined. Without going into details, we will go through these spaces and generalize the results needed for our purposes to the case of complex-valued four-dimensional numbers.

Let’s consider the space $M_2$, in which the multiplication of numbers $X=(x_1,x_2,x_3,x_4)$ and $Y=(y_1,y_2,y_3,y_4)$ is defined as follows: where $Z(z_1,z_2,z_3,z_4)=X·Y$. If we set $x_j=a_j+b_{j}i$, $y_j=c_j+d_{j}i$, $j=1,2,3,4$, where i is the imaginary unit, then this multiplication can be rewritten as where

One can readily verify that the multiplication of four-dimensional numbers defined in this way is commutative. To the four-dimensional number $X=(x_1,x_2,x_3,x_4)\in C^4$ we associate the matrix The mapping $F:X\to M_{20}\left(X\right)$ is bijective and onto. Indeed, two different numbers X and Y correspond to different matrices, and for any matrix in the form of (8), a corresponding four-dimensional number from $C^4$ can be found.

Theorem 1. The set of all matrices in the form of (8) is closed with respect to the operations of matrix addition, subtraction, multiplication, and multiplication by scalar. For the mapping $F:X\to M_{20}\left(X\right)$ the relationships $F(X\pm Y)=F\left(X\right)\pm F\left(Y\right)$, $F\left(XY\right)=F\left(X\right)F\left(Y\right)$ hold for any $X\in C^4,Y\in C^4$.

The proof is conducted by direct verification.

Thus, there is a bijection between the space of four-dimensional numbers and the space of matrices of the form (8), which preserves arithmetic operations, meaning the existing bijection is a homomorphism. From Theorem 1, it also follows that the operation of matrix multiplication of the form (8) is commutative.

It is further noted that if we multiply the j-th row and j-th column of matrix (8) by -1, we obtain another matrix with the same properties as the matrix $M_{20}$, that is, the statements of Theorem 1 remain valid. Moreover, if we multiply the j-th row and three columns of matrix $M_{2}0$, with indexes not equal to j, by -1, we also get a matrix corresponding to the multiplication of four-dimensional numbers (7) and possessing the properties of matrix $M_{2}0$. The matrix transposed to $M_{2}0$ also possesses all the properties of matrix $M_{2}0$. To describe such operations, let us denote by $M_{20}^{(j,k)}$, where j and k are one, two, or three indices with values from 1 to 4, the matrix obtained by multiplying by -1 the rows with numbers from index j and columns with numbers from index k. For example, $M_{20}^{(24,134)}$ is a matrix obtained from matrix $M_{2}0$ by multiplying the second and fourth rows by -1, and also by multiplying the first, third, and fourth columns by -1. Let’s describe all possible operations that lead to matrices for which the statements of Theorem 1 are valid. It is easy to verify that such operations are operations of the following types: $M_{20}^{(j,j)}$, $M_{20}^{(j,klm)}$, $M_{20}^{(klm,j)}$, $M_{20}^{(klm,klm)}$, $M_{20}^{(jk,jk)}$, $M_{20}^{(jk,lm)}$, where $j,k,l,m$ are pairwise distinct indices with values from 1 to 4. In addition, the operation of transposing a matrix also does not change its properties.

The number of different operations of the form $M_{20}^{(j,j)}$ is four when $j=1,2,3,4$, respectively, we get four new matrices:

But as we can easily notice, that is, these matrices lie in one group. Similarly, it can be shown that the matrices from the groups $M_{20}^{(j,klm)}$, $M_{20}^{(klm,j)}$, $M_{20}^{(klm,klm)}$, totaling 12 (4 in each group), lie in the same group. In addition, we include in this group all transposed matrices of this group, as it is easy to check that they also lie in this group. And the matrices from the groups $M_{20}^{(jk,jk)}$, $M_{20}^{(jk,lm)}$, totaling 12 (6 in each group) and the transposed matrices to them, also lie in one group, but different from the first group. For example, the matrices $M_{20}^{(12,12)}$ and $M_{20}^{(12,34)}$ respectively have the form: from which it follows that $M_{20}^{(12,12)}(x_1,x_2,x_3,x_4)=-M_{20}^{(12,34)}(x_1,x_2,x_3,x_4)$, but from $M_{20}^{(12,12)}$ it is impossible to obtain $M_{20}^{(1,1)}$ or another matrix from the first group. Thus, there are two groups of matrices that are closed with respect to the operations of addition, multiplication, and these operations are commutative. Any matrix from the corresponding group can be taken as a representative of these groups. As a representative of the second group, we take the matrix $M_{20}$, and as a representative of the first group, we take, for example, the matrix $M_{20}^{(1,1)}$, which we denote by $M_{21}$:

Thus, to each four-dimensional number $X\in C^4$ two matrices $M_{20}$ and $M_{21}$ can be associated, that is, to define two mappings, $F_{20}:X\to M_{20}\left(X\right)$ and $F_{21}:X\to M_{21}\left(X\right)$, which are bijective and onto. The products of matrices from one class are closed with respect to the operations of addition and multiplication, and the multiplication operation corresponds to the multiplication of four-dimensional numbers (7).

Note. Other transformations $M_{20}^{(j,k)}$ of the matrix $M_{20}$ can be considered and used to build unitary operators.

Now, let us consider the space $M_3$, where the multiplication operation of four-dimensional numbers $X=(x_1,x_2,x_3,x_4)$ and $Y=(y_1,y_2,y_3,y_4)$ is defined as follows [2]: where $Z=(z_1,z_2,z_3,z_4)=X·Y$. A detailed exposition of the algebra and analysis over the four-dimensional space of real numbers $M_3$ is presented in monograph [1]. To the four-dimensional number $X=(x_1,x_2,x_3,x_4)\in C^4$ we associate the matrix The mapping $F_{30}:X\to M_{30}\left(X\right)$ is bijective and onto. Indeed, two different numbers X and Y correspond to different matrices, and for any matrix of form (10), one can find the corresponding four-dimensional number from $C^4$. For the matrix $M_{30}$ and the mapping $F_{30}$ , the statements of Theorem 1 hold true. Similarly to the previous case, by considering transformations of the form $M_{30}^{(j,j)}$, $M_{30}^{(j,klm)}$, $M_{30}^{(klm,j)}$, $M_{30}^{(klm,klm)}$, $M_{30}^{(jk,jk)}$, $M_{30}^{(jk,lm)}$, where $j,k,l,m$ are pairwise distinct indices with values from 1 to 4, we find that there exist two groups of matrices corresponding to the commutative multipli-cation in the space $M_3$ and satisfying the conditions of Theorem 1. One group is represented by the matrix (10), and the other group by the following matrix $M_{31}\left(X\right)$: Similarly, for the space $M_4$, with commutative multiplication [2] we obtain two groups of matrices $M_{40}\left(X\right)$ and $M_{41}\left(X\right)$:

Proceeding with similar reasoning and corresponding calculations for spaces $M_5,M_6,M_7$, , the multiplication operations of which are defined respectively as [2] we determine the corresponding groups of matrices satisfying the conditions of Theorem 1 for the indicated commutative multiplication operations: for the space $M_5$, for the space $M_6$, for the space $M_7$. Note that the foundations of the four-dimensional space of real numbers M5 can be found in works [5,6].

A careful examination of the obtained matrices reveals that

This implies that there are in fact six independent groups of matrices corresponding to the six spaces of four-dimensional numbers $M_j,(j=2,3,...,7)$. As such matrices, we shall take $M_{20}$, $M_{30}$, $M_{40}$, $M_{50}$, $M_{60}$ and $M_{70}$, defined by equations (8), (10) - (14). Each of these matrices is bijective to the space of four-dimensional complex-valued numbers, closed with respect to the multiplication operation, and forms an abelian group with respect to the matrix multiplication operation. Moreover, as evident from the construction, no other abelian groups of matrices exist.

In the previous section, we constructed six abelian groups of matrices with elements formed from the components of a four-dimensional number. Based on these matrices, it is possible to construct gates for two-qubit quantum systems; in other words, under certain additional conditions, the constructed matrices transform into unitary operators. Let us formulate the corresponding conditions.

Theorem 2. Let $X=(x_1,x_2,x_3,x_4)\in C^4$. Let the components of the complex numbers $x_j=u_j+v_{j}i(j=1,2,3,4)$ satisfy the conditions: Then the matrix $M_{20}\left(X\right)$ is a two-qubit gate.

Proof. The system (15) is consistent and has an infinite number of solutions. Consider the Hermitian conjugate matrix $M_{20}^{*}$ to the matrix $M_{20}$: and multiply it by the matrix $M_{20}$: Using relations (15) we obtain that $M_{20}·M_{20}^{*}=E$, where E is the identity matrix of size $4\times 4$, hence $M_{20}$ is a unitary matrix.

Thus, when conditions (15) are met, the group of matrices $M_{20}\left(X\right)$ forms a commutative group of gates for a two-qubit quantum system.

Corollary 1. The group of commutative gates with real elements is of the form where Corollary 2. The group of commutative gates with purely imaginary elements is of the form where Theorem 3. Let $X=(x_1,x_2,x_3,x_4)\in C^4$. Let the components of the complex numbers $x_j=u_j+v_{j}i(j=1,2,3,4)$ satisfy the conditions: Then the matrix $M_{30}\left(X\right)$ is a two-qubit gate.

The proof is analogous to the proof of Theorem 2.

This theorem coincides with Theorem 1 from work [3]. Analogous to Corollary 1 and Corollary 2 to Theorem 2, we can write the forms of gates with real and imaginary elements for the matrix $M_{30}$.

Similarly, we consider matrices $M_{40},M_{50},M_{60}$ and $M_{70}$ and list the corres-ponding conditions for their unitarity.

Theorem 4. Let $X=(x_1,x_2,x_3,x_4)\in C^4$. Let the components of the complex numbers $x_j=u_j+v_{j}i(j=1,2,3,4)$ satisfy the conditions: Then the matrix $M_{40}\left(X\right)$ is a two-qubit gate.

Theorem 5. Let $X=(x_1,x_2,x_3,x_4)\in C^4$. Let the components of the complex numbers $x_j=u_j+v_{j}i(j=1,2,3,4)$ satisfy the conditions: Then the matrix $M_{50}\left(X\right)$ is a two-qubit gate.

Theorem 6. Let $X=(x_1,x_2,x_3,x_4)\in C^4$. Let the components of the complex numbers $x_j=u_j+v_{j}i(j=1,2,3,4)$ satisfy the conditions: Then the matrix $M_{60}\left(X\right)$ is a two-qubit gate.

Theorem 7. Let $X=(x_1,x_2,x_3,x_4)\in C^4$. Let the components of the complex numbers $x_j=u_j+v_{j}i(j=1,2,3,4)$ satisfy the conditions: Then the matrix $M_{70}\left(X\right)$ is a two-qubit gate.