A Novel Approach to State-to-State Transformation in Quantum Computing

1. Introduction

Quantum state transformation - changing a quantum state through quantum operations or gates - is fundamental to many quantum algorithms. It enables the manipulation and evolution of quantum information, facilitating complex computations. Mastery of state transformations is crucial for multi-qubit systems, which utilize specific gates to create entangled states and perform calculations. It is also vital for quantum error correction to ensure the coherence and fidelity of quantum information. Effective and low-noise state preparation procedures are essential for scalable distribution loading, which underpins a wide range of algorithms that provide a quantum advantage. Therefore, the important step in quantum computation in many applications is the multi-qubit superposition transformation. This is the problem of transformation of one quantum state into another one. Many papers have been published related to this problem [1,2,3,4,5,6,7,8]. Analyzing these works, it is important to note the following. The view has been formed that operations on three or more qubits are very complex for quantum computers. Therefore, when creating quantum circuits, much attention is paid to one- and two-qubit operations, or gates. As a result, each multi-qubit gate, as a unitary transformation, is represented by a chain of single-qubit and dual-qubit gates together with a chain of CNOTs, which may include the long-range CNOTs. The schemes of many simple (for a classical computer) operations have become very complex and are full of such switches of operations from one qubit to another. Therefore, much attention is paid to reducing the number of such gates in quantum circuits [9,10,11]. As is known, the number of elementary rotation gates is around the number ${2(2}^r-1)$ and CNOTs in number ${2(2}^r-r-1)$ [12,13,14].

The solution to the problem of multi-qubit state-to-state transformation $|\mathit{x}⟩\to |\mathit{y}⟩$ traditionally lies in the idea of ​​preparing the desired state from the conventional basis state, $|{0⟩}^{\oplus r}$. Here, $r\ge 1$ is the dimension of the superposition, or the number of qubits it contains. We note the simple method with Givens rotations, $\left(x,y\right)\to \left(\pm \sqrt{x^2+y^2},0\right)$. They rotate points inside the unit circle into the interval $[-\mathrm{1,1}]$ on the $X$-axis. In matrix theory, such rotations are used in each step of the QR decomposition of a square matrix, when transforming first columns $\mathit{z}$ of the submutrix into the vector of form of $(\pm ‖\mathit{z}‖,0, 0,\dots ,0)$. We consider only real vectors with norm 1, as for quantum superpositions with real amplitudes. Since it is possible to perform transformations of two vectors $\mathit{x}$ and $\mathit{y}$ into the unit vector ${\mathit{e}}_1=(\pm \mathrm{1,0}, 0,\dots ,0)$ of same dimension,

$$U_x:\mathit{x}\to {\mathit{e}}_1 \mathrm{a}\mathrm{n}\mathrm{d} U_y:\mathit{y}\to {\mathit{e}}_1,$$

(1)

the transformation $\mathit{x}\to \mathit{y}$, can be fulfilled by the unitary transform ${T=U}_y^{-1}{U}_x$. Both transforms $U_x$ and $U_y$ are unitary, and $U_y^{\prime }$ is the inverse of $U_y$. Moreover, each of them requires ${(2}^r-1)$ rotations with almost the same number of CNOTs, when implementing calculations with 1- and 2-qubit gates [8]. If CNOTs operate not on the adjacent, or nearest-neighbor, bit planes (BP), they can be implemented by a cascade of CNOTs operating on the nearest neighbor bit planes. The permutations with Gray code can be used for this purpose [14,15]. The number of all CNOTs is estimated as $(10/3)2^r$ [7].

In this work, we analyze the method of rotations and describe the best and fastest, in our opinion, way of elementary rotations for transforming the state. The quantum states, or superpositions of qubits, are considered with real amplitudes. The state-to-state transformation can be performed in one step, which is equivalent to implementing only one transformation which is similar to $U_x$, instead of two transforms $U_x$ and $U_y$ in the Eq. 1. This is the main goal of this work, not counting and reducing the number of CNOT gates; a lot of work has already been done in this direction. We introduce and describe the new concept of the quantum signal-induced heap transform (QsiHT) which is the analogue of the discrete signal-induced heap transform (DsiHT) [16,17]. With this transform, we show

• How to implement $r$-qubit state-to-state transformation with real amplitudes, by using only ${(2}^r-1)$ elementary rotations, each with only one angle. This number is half as much as the best-known estimation of ${2(2}^r-1)$ rotations.
• Visual numerical examples of preparing states $\mathit{x}\to \mathit{y}$ for two, three, and four qubits.
• How to initiate any multi-qubit superposition of qubits (without operation of tensor product).
• The importance of the path in $r$-qubit QsiHT and existence of the fast paths for effective computing the QsiHT. For large multi-qubit superpositions, there are various fast paths (with their number increasing with the number of qubits), and we are confident that among them it is possible to choose the most convenient path for implementing state transformations in the topology (architecture) of quantum systems.
• How to build the simple quantum circuits for the $r$-qubit QsiHT.

2. One- and Two-Qubit Operation

In this section, we briefly describe the concept of local gates on single qubits in the quantum superposition [18,19]. The 1- and 2-qubit operators are described by 2$\times$2 and 4$\times$4 unitary real or complex matrices, respectively. Let us find out what it means to apply a gate to one qubit in a 4-qubit superposition. Consider the conventional basis of states. In the general case, a 2-qubit superposition $|\mathit{z}⟩={\mathrm{z}}_0|00⟩+{\mathrm{z}}_1|01⟩+{\mathrm{z}}_2|10⟩+z_3|11⟩$ may be the Kronecker product of two 1-qubits $|\mathit{x}⟩={\mathrm{x}}_0|00⟩+{\mathrm{x}}_1|01⟩$ and $|\mathit{y}⟩={\mathrm{y}}_0|00⟩+{\mathrm{y}}_1|01⟩$, or maybe not. In the first case, the 2-qubit $|\mathit{z}⟩=|\mathit{x}⟩\otimes |\mathit{y}⟩$ in non-entangled, and in the second case the 2-qubit $|\mathit{z}⟩$ is entangled. Let $U_2$, $V_2$, and $U_4$ denote some unitary operations, that is, 2$\times$2, 2$\times$2, and 4$\times$4 matrices, respectively. If the operation on the non-entangled 2-qubit, that is, 2-qubit superposition, is described by

$$U_4|\mathit{z}⟩=\left(U_2\left|\mathit{x}⟩,V_2\right|\mathit{y}⟩\right)=\left(U_2|\mathit{x}⟩\right)\otimes \left(V_2|\mathit{y}⟩\right)={(U}_2\otimes V_2)\left(|\mathit{x}⟩\otimes |\mathit{y}⟩\right),$$

(2)

then, it is the qubit-wise operation, or the operation on two qubits $|\mathit{x}⟩$ and $|\mathit{y}⟩$ of $|\mathit{z}⟩$. This operation includes the operations on single qubits, too,

$$U_4|\mathit{z}⟩=\left(U_2|\mathit{x}⟩,|\mathit{y}⟩\right)=U_2|\mathit{x}⟩\otimes |\mathit{y}⟩ \mathrm{a}\mathrm{n}\mathrm{d} U_4|\mathit{z}⟩=\left(\left|\mathit{x}⟩,V_2\right|\mathit{y}⟩\right)=\left|\mathit{x}⟩\otimes V_2\right|\mathit{y}⟩.$$

(3)

As an example, consider the unitary matrix $U_2=\left[a,b;c,d\right],$ with the real coefficients and $a^2+b^2+c^2+d^2=1.$ The first operation is described in the matrix form as

$$U_4|\mathit{z}⟩=U_2|\mathit{x}⟩\otimes |\mathit{y}⟩=\left(U_2\otimes I_2\right)\left(|\mathit{x}⟩\otimes |\mathit{y}⟩\right)=\left[\begin{array}{cccc}a& & b& \\ & a& & b\\ c& & d& \\ & c& & d\end{array}\right]\left[\begin{array}{c}{x}_0{y}_0\\ x_0{y}_1\\ x_1{y}_0\\ x_1{y}_1\end{array}\right].$$

(4)

The circuit element of this gate ${U_4=U}_2\otimes I_2$ is shown in Figure 1 in part (a). It is the local gate on the first qubit. Here, $I_2$ is the 2$\times$2 identity matrix. In the diagram this gate (operator) is connected to the first abstract quantum wire (quantum bit). The circuit element shown in part (b) as the operation on the second qubit. This operation can be written as $I_2\otimes U_2.$ Note that in general the input, or the 2-qubit superposition $|\mathit{z}⟩$, of the circuit does not necessarily have to be entangled. Any 2-qubit superposition is assumed as an input.

The 1-qubit gate controlled by one qubit (or bit) are described by the operation of the Kronecker sum of matrices. The circuit elements in Figure 1(c) and (d) are described in the matrix form as follows:

$$C_1\left(U_2\right)=I_2\oplus U_2=\left[\begin{array}{cccc}1& & & \\ & 1& & \\ & & a& b\\ & & c& d\end{array}\right] \mathrm{a}\mathrm{n}\mathrm{d} C_0\left(U_2\right)=U_2\oplus I_2=\left[\begin{array}{cccc}a& b& & \\ c& d& & \\ & & 1& \\ & & & 1\end{array}\right].$$

Note that if a 2-qubit superposition $|\mathit{z}⟩$ is nonentangled, then $C_{\mathrm{0,1}}\left(U_2\right)|\mathit{z}⟩$ may be entangled. These operations are performed on the first or last two amplitudes of $|\mathit{z}⟩=|\mathit{x}⟩\otimes |\mathit{y}⟩$, but not on the qubit $|\mathit{x}⟩$ or $|\mathit{y}⟩$. The matrices of the local gates $U_2\otimes I_2$ and $I_2\otimes U_2$ are more complex than the matrices of the controlled gates $C_0\left(U_2\right)$ and $C_1\left(U_2\right)$. For instance,

$$I_2\otimes U_2=\left[\begin{array}{cccc}a& b& & \\ c& d& & \\ & & a& b\\ & & c& d\end{array}\right]=C_1\left(U_2\right)C_2\left(U_2\right).$$

3. Definition of the DsiHT

In this section, we describe the concept of the DsiHT with Givens rotations. We note the following feature of many unitary transformations used in engineering and science. The quantum Fourier transform (QFT) is well-known in quantum computation [20,21,22]. This transform is used in the solution of such problems, as the Shor’s algorithm for integer factorization, the quantum phase estimation algorithm, systems of linear equations and others [23,24,25,26]. The QFT is the quantum analogue of the $2^r$-point discrete Fourier transform, where integer $r>1$. The basis functions of the transform are $2^r$complex exponential functions on the unit circle. Many other discrete unitary transformations also use sets of defined in advance basis functions. We mention for example, the cosine, Hadamard, hartley, cosine and slant transforms [27,28,29]. The $N$-point discrete signal-induced heap transform (DsiHT) generates the basis by a given signal of length $N$. This signal is called the generator of the transform [30,31]. The DsiHT with its simple and fast algorithm (for any order $N$ of transformation) is successfully used in many applications in grayscale and color image processing. In addition, the DsiHT can be used for QR-decomposition of square matrices [32]. Together with the generator, the transform uses a parameter which defines the order of processing the input signals. This parameter, or path, is very important and allows us to obtain different QR-decompositions, with diagonal matrices $R$ for unitary matrices.

The key of the DsiHT is the generator together with the path. The number of generators can be more than one, but in this work, we focus on the transforms using one pre-selected generator. The selection of the generator depends on the application. For example, in image enhancement, the generators can be defined by the mean or median values ​​along the rows or columns of an image result in very good quality images. The generator-signal is denoted by $\mathit{x}=(x_0,x_1,x_2,\dots , x_{N-1})$. The length of the signal is considered to be a power of 2, $N=2^r, r>1.$ The basic building blocks that make up the transform are 2$\times$2 transforms, which can be linear or nonlinear [30]. In this work, we consider 2$\times$2 rotation matrices, or Givens rotations. The transform is calculated by using $(N-1)$ Givens rotations by angles ${\vartheta }_{\mathrm{k}}$. In this case, the signal-generator is presented in its angular representation

$$\mathit{x}\to A_{\mathit{x}}=\left\{{\vartheta }_1,{\vartheta }_2,{\vartheta }_3,\dots ,{\vartheta }_{N-1}\right\}$$

(5)

plus, the norm, $||\mathit{x}||$. The signal is considered real. The DsiHT is a unitary transform, $H$, and when applying to the generator, it results in the vector $(||\mathit{x}||,\mathrm{0,0},\dots , 0 )$. For a generator with the unit norm or after normalization, this $N$-dimensional vector $(\mathrm{1,0},0,\dots , 0)$ presents the computational basis state $|00\dots 0⟩=|{0⟩}^r.$ Thus, the following takes place: $H\mathit{x}=|{0⟩}^{\oplus r}$.

To define the DsiHT, $(N-1)$ rotations are performed on the generator and at each stage of the transformation one of the two rotation outputs can be reset to zero. For instance, the generator data can be processed in order (it is the path) $x_0$ and $x_1$, then the updated $x_0$ with $x_2$, the next updated $x_0$ with $x_3$, and so on. The illustration of the calculation with the 8-point generator is given in Figure 2. Here, $N=8$ and 2$\times$2 transforms $T_k,$ $k=1:7,$ are rotations by angles ${\vartheta }_1,{\vartheta }_2,{\vartheta }_3, \dots , {\vartheta }_7$. The DsiHT with this regular path is called the DsiHT with weak wheel-carriage [16,31]. The last value is the updated 7 times value of $x_0$ with the splash of energy of the generator, as if we were collecting energy, $E\left[\mathit{x}\right]=‖x‖,$ in one heap.

The DsiHT can be considered as a two-level unitary transformation. First, the angles ${\vartheta }_k$ are calculated from the generator. Then, at each stage $k$ of calculation, the transform $T_k$ is applied to the input signal $\mathit{z}=(z_0,z_1,z_2,\dots ,z_{N-1})$ in the way (path) as for the generator. For the DsiHT with the week carriage-wheel, this two-level transform is illustrated in Figure 3, for the $N=4$ case. The results of calculations are

$$x\to (x_0^{\left(3\right)},\mathrm{0,0},0)=(‖x‖,\mathrm{0,0},0), z\to (z_0^{\left(3\right)},z_1^{\prime },z_2^{\prime },z_3^{\prime }).$$

4. The DsiHT with Strong 2-Wheel Carriage

We will consider the concept of the strong DsiHT, or the DsiHT with strong 2-wheel carriage (for more detail, see [16,17]). The path of this DsiHT differs from the path used in the above DsiHT with the week wheel-carriage. As an example, the block-diagram of the 8-point strong DsiHT on its generator is shown in Figure 4. The case of real signals is considered. In this example, on the first step, the angle of the rotation is calculated from the conditions

$$T_1\left[\begin{array}{c}{x}_6\\ x_7\end{array}\right]=\left[\begin{array}{c}{x}_0^{\left(1\right)}\\ 0\end{array}\right] \mathrm{a}\mathrm{n}\mathrm{d} {|x}_0^{\left(1\right)}|=\sqrt{x_6^2+x_7^2}.$$

(6)

This transform is the Givens rotation calculated by

$$T_1=T_{{\phi }_1}=\left[\begin{array}{rr}\cos {\vartheta }_1& -{\sin \vartheta }_1\\ {\sin \vartheta }_1& \cos {\vartheta }_1\end{array}\right],$$

(7)

with the angle ${\vartheta }_1=-\arctan \left(x_7{/x}_6\right)$. If $x_6=0,$ the angle ${\vartheta }_1=-\pi /2,$ or $\pi /2.$ Thus, the angle of rotation is defined from the angular equation $x_6\sin {\vartheta }_1+x_7\cos {\vartheta }_1=0.$Then, the first output, $x_0^{\left(1\right)}$, of the transform is calculated in Eq. 6. On the next step, the second angle of rotation, $T_2=T_{{\vartheta }_2}$ is calculated from the conditions

$$T_2\left[\begin{array}{c}{x}_5\\ x_0^{\left(1\right)}\end{array}\right]=\left[\begin{array}{c}{x}_0^{\left(2\right)}\\ 0\end{array}\right] \mathrm{a}\mathrm{n}\mathrm{d} {|x}_0^{\left(2\right)}|=\sqrt{x_5^2+{\left(x_0^{\left(1\right)}\right)}^2}.$$

Thus, ${\vartheta }_2=-\arctan \left(x_5/x_0^{\left(1\right)} \right).$ Continuing similar calculations, the last angle ${\vartheta }_7=-\arctan \left(x_0^{\left(6\right)}{/x}_0\right)$ of rotation $T_7=T_{{\phi }_7}$ is calculated and then the last value $x_0^{\left(7\right)},$

$$T_7\left[\begin{array}{c}{x}_0\\ x_0^{\left(6\right)}\end{array}\right]=\left[\begin{array}{c}{x}_0^{\left(7\right)}\\ 0\end{array}\right].$$

The value $x_0^{\left(7\right)}$ contains information of the energy of the generator.

The input 8-point signal $\mathit{z}=\left(z_0,z_1,z_2,\dots , z_7\right)$ is processed by the same seven rotations in the same order (path). This 2-level strong carriage-wheal DsiHT with processing the input signal $\mathit{z}$ with the generator is shown in Figure 5.

In the general case, the $N$-point DsiHT uses $(N-1)$ basic transformations $T_{{\vartheta }_k}, k=1:(N-1)$, regardless of the path of the transform. The transform of the signal-generator equals

$$H:\mathit{x}=\left(x_0,x_1,x_2,\dots ,x_{N-1}\right)\to \left(y_0,\mathrm{0,0},\dots ,0\right),$$

(8)

where $y_0$ is equal to plus/minus energy $E\left[\mathit{x}\right]=‖\mathit{x}‖$ of the signal. We consider $y_0=\sqrt{x_0^2+x_1^2+\cdots {+x}_{N-1}^2}$; if it is negative, the last angle in rotation can be changed as ${\vartheta }_{N-1}={\vartheta }_{N-1}+\pi .$ Thus, the transformation, as well as the signal-generator, are uniquely determined by the data of $E\left[\mathit{x}\right]$ and the angles of rotations. In other words, the transformation is decoded into the vector $\left(‖\mathit{x}‖,{\vartheta }_1,{\vartheta }_2, \dots ,{\vartheta }_{N-1}\right).$ For the signal-generator with norm 1, the angular representation holds, $\mathit{x}\to A_{\mathit{x}}=\left\{{\vartheta }_1,{\vartheta }_2, \dots ,{\vartheta }_{N-1}\right\}.$

The output is the transform of the input

$$H\left(\mathit{z}\right)=\left(z_0^{\prime },z_1^{\prime },z_2^{\prime },\dots ,z_{N-2}^{\prime },z_{N-1}^{\prime }\right), z_0^{\prime }=z_0^{\left(N-1\right)}.$$

(9)

Note that the above algorithm can be split by two parts. In the first part, all angles ${\vartheta }_k$ of rotation can be calculated from the signal-generator. This is about the angular representation of the generator, $A_{\mathit{x}}$. In the second part, the DsiHT of the input signal $\mathit{z}$ is calculated, $H\left(\mathit{z}\right).$

We now describe the DsiHT in the matrix form. For $i=1:\left(N-1\right),$we denote by $R_{\vartheta ;(i,i+1)}$ the block-diagonal matrix with all 1 in the diagonal and the matrix of rotation $R_{\vartheta }$ in the cell $\left(i,i+1\right)\times (i,i+1)$, that is,

$$R_{\vartheta ;(i,i+1)}=I_i\oplus R_{\vartheta }\oplus I_{N-i-2} \mathrm{a}\mathrm{n}\mathrm{d} R_{\vartheta ;\left(\mathrm{0,1}\right)}=R_{\vartheta }\oplus I_{N-2}.$$

(10)

The matrix $R_{\vartheta ;(i,i+1)}$describes the rotation on neighbor bit planes $i$ and $(i+1)$. For instance, in the $N=4$ case, consider the rotations

$$R_{\vartheta ;\left(\mathrm{0,1}\right)}=\left[\begin{array}{cccc}\cos \vartheta & -\sin \vartheta & 0& 0\\ \sin \vartheta & \cos \vartheta & 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right]={\left[\begin{array}{cc}{R}_{\vartheta }& 0\\ 0& I_2\end{array}\right]=R}_{\vartheta }\oplus I_2 \mathrm{a}\mathrm{n}\mathrm{d} R_{R_{\vartheta ;\left(\mathrm{1,2}\right)}}=1\oplus R_{\vartheta }\oplus 1.$$

The bit planes 0 and 1, or 00 and 01, are adjacent planes, that is, the plane numbers differ by only one bit. The bit planes 1 and 2, or 01 and 10, are not adjacent planes.

In the matrix form, the DsiHT, $H_N$, is written as the multiplication of $(N-1)$ block-diagonal matrices with rotations on adjacent planes,

$$H_N=R_{{\vartheta }_{N-1;\left(\mathrm{0,1}\right)}}{R}_{{\phi }_{\vartheta -2;\left(\mathrm{1,2}\right)}}{R}_{{\phi }_{\vartheta -3;\left(\mathrm{2,3}\right)}}\dots R_{{\vartheta }_{2;\left(N-3,N-2\right)}}{R}_{{\vartheta }_{1;\left(N-2,N-1\right)}}.$$

(11)

In the $N=4$ case, the 4-point strong DsiHTs with three rotations is written as

$$H_4=R_{{\vartheta }_{3;\left(\mathrm{0,1}\right)}}{R}_{{\vartheta }_{2;\left(\mathrm{1,2}\right)}}{R}_{{\vartheta }_{1;\left(\mathrm{2,3}\right)}}=\left(R_{{\vartheta }_3}\oplus I_2\right)\left({1\oplus R}_{{\vartheta }_2}\oplus 1\right)\left(I_2\oplus R_{{\vartheta }_1}\right).$$

(12)

Now, we consider an example with the 4-point DsiHT generated by a signal $\mathit{x}$.

Example  1:

Let $\mathit{x}$ be the vector-signal $\mathit{x}=(\mathrm{1,2},\mathrm{2,1})$ with energy $E\left[\mathit{x}\right]=\sqrt{10}$ and let the input signal be $\mathit{z}=(\mathrm{1,2},\mathrm{4,2})$ with energy $E\left[\mathit{z}\right]=5$. The DsiHT with the generator $\mathit{x}$ has the following matrix:

$$H_4=\left[\begin{array}{rrrr}0.3162& 0.6325& 0.6325& 0.3162\\ -0.9487& 0.2108& 0.2108& 0.1054\\ 0& -0.7454& 0.5963& 0.2981\\ 0& 0& -0.4472& 0.8944\end{array}\right]=D\left[\begin{array}{rrrr}1& 2& 2& 1\\ -9& 2& 2& 1\\ 0& -2.5& 2& 1\\ 0& 0& -1& 2\end{array}\right].$$

(13)

Here, the diagonal matrix $D=\{$0.3162,0.1054,0.2981,0.4472$\}$. The angles for this transformation are $\left\{{\vartheta }_1,{\vartheta }_2,{\vartheta }_3\right\}=\{-26.5651°, -48.1897°, -71.5651°\}$. One can note that the normalized vector-generator $(\mathrm{1,2},\mathrm{2,1})/\sqrt{10}$ lies in the first row of this matrix. The transformation is described by the following decomposition by three rotations:

$$\begin{array}{l}{H}_4=R_{{\vartheta }_{3;\left(\mathrm{0,1}\right)}}{R}_{{\vartheta }_{2;\left(\mathrm{1,2}\right)}}{R}_{{\vartheta }_{1;\left(\mathrm{2,3}\right)}}\\ =\underset{{\vartheta }_3=-71.56°, { \vartheta }_2=-48.19°, { \vartheta }_1=-26.56° }{\underbrace{\left[\begin{array}{cccc}\cos {\vartheta }_3& -\sin {\vartheta }_3& 0& 0\\ \sin {\vartheta }_3& \cos {\vartheta }_3& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right]\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& \cos {\vartheta }_2& -\sin {\vartheta }_2& 0\\ 0& \sin {\vartheta }_2& \cos {\vartheta }_2& 0\\ 0& 0& 0& 1\end{array}\right]\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& \cos {\vartheta }_1& -\sin {\vartheta }_1\\ 0& 0& \sin {\vartheta }_1& \cos {\vartheta }_1\end{array}\right].}}\end{array}$$

Consider the permutations $P=(\mathrm{0,3},\mathrm{2,1})$ and $P\prime =(\mathrm{0,1},\mathrm{2,3})$ with the matrices

$$P=\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\\ 1& 0& 0& 0\end{array}\right] \mathrm{a}\mathrm{n}\mathrm{d} P\prime =\left[\begin{array}{cccc}0& 0& 0& 1\\ 1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\end{array}\right],$$

respectively. These two permutations describe the cyclic shifts of a 4-point vector. The following holds for the rotation on the first- and second-bit planes: $R_{{\vartheta }_{2;\left(\mathrm{1,2}\right)}}=P{R}_{{\vartheta }_{2;\left(\mathrm{2,3}\right)}}P\prime$, that is,

$$R_{{\vartheta }_{2;\left(\mathrm{1,2}\right)}}=\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& 1\\ 1& 0& 0& 0\end{array}\right]\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& \cos {\vartheta }_2& -\sin {\vartheta }_2\\ 0& 0& \sin {\vartheta }_2& \cos {\vartheta }_2\end{array}\right]\left[\begin{array}{cccc}0& 0& 0& 1\\ 1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\end{array}\right].$$

Therefore, the matrix of the 4-point DsiHT is described by three controlled gates as

$$H_4=\left(R_{{\vartheta }_3}\oplus I_2\right)P\prime \left(I_2\oplus R_{{\vartheta }_2}\right)P\prime {(I}_2\oplus R_{{\vartheta }_1}).$$

(14)

The transform of the generators is $H_4{\mathit{x}}^{\prime }={\left(\mathrm{1,0},\mathrm{0,0}\right)}^{\prime }$, and for the input signal, $H_4\mathit{z}\prime =\left(\mathrm{4.7434,0.5270,1.4907,0}\right)\prime$. The normalized vectors, $\mathit{x}=\mathit{x}/\sqrt{10}$ and $\mathit{z}=\mathit{z}/5$ can be considered as 2-qubit superpositions

$$|\mathit{x}⟩=(|00⟩+2|01⟩+2|10⟩+|11⟩)/\sqrt{10}, \mathrm{a}\mathrm{n}\mathrm{d}$$

$$|\mathit{z}⟩=(|00⟩+2|01⟩+4|10⟩+2|11⟩)/5.$$

The 2-qubit QsiHT of these superpositions are $H_4|\mathit{x}⟩=|00⟩$ and

$$H_4|\mathit{z}⟩=(4.7434|00⟩+0.5270|01⟩+1.4907|10⟩+0|11⟩)/5.$$

The quantum scheme of the 2-qubit quantum QsiHT, is given in Figure 6. The flowchart sorts denote the permutations.

IBM’s Qiskit framework is an open-source library in Python for quantum computing developed by IBM [33,34]. It provides tools for building quantum circuits, simulating their behavior, and running them on actual quantum hardware through IBM’s backend. IBM’s Qiskit framework was used to construct and simulate the 2-qubit QsiHT quantum circuit designed to prepare the normalized quantum state corresponding to the vector $x=(\mathrm{1,2},\mathrm{2,1})\prime$. The circuit was built using the decomposition of rotation angles through the QsiHT to encode the target amplitudes and its inverse was applied to map the target state onto a computational basis. The simulation involved extracting both the ideal state vector amplitudes and the corresponding measurement probabilities. These were compared against shot-based measurements from simulated runs to assess the accuracy of the circuit. The mean-square root errors (MSRE) between the expected and observed results were computed to quantify the performance of the state preparation. The results of simulation of this circuit in Qiskit are shown in Table 1.

5. Qubits Initiation by the DsiHT

In this section, we consider the inverse DsiHT and describe a simple quantum scheme for initiation of a multi-qubit superposition. Then in Section 6 we show how important it is to choose such a path to make this scheme more efficient.

Let $x$ be the $2^r$-D unit vector $\mathit{x}=(x_0,x_1, \dots ,x_{2^r-1})$ which represents the $r$-qubit superposition

$$|\mathit{x}⟩=x_0|0⟩+x_1|1⟩+\dots +x_{2^r-1}|2^r-1⟩, {\left|x_0\right|}^2+{\left|x_1\right|}^2+\dots +{\left|x_{2^r-1}\right|}^2=1,$$

Here, the components of the vector are real and complex, and $\{|k⟩, k=0:2^r-1\}$ is the conventional basis of states. We consider the case when $\mathit{x}$ is real. For the DsiHT generated by this vector, we have

$$H{\mathit{x}}^{\mathbf{\text{'}}}=H{\left(x_0,x_1,\dots ,x_{2^r-1}\right)}^{\prime }={\left(\mathrm{1,0},\dots ,0\right)}^{\prime }={|0⟩}^{\otimes r}.$$

(15)

In Example 1, such DsiHT is described for the $N=4$ case, or the 2-qubit superposition $|\mathit{x}⟩=(|0⟩+2|1⟩+2|2⟩+|3⟩)/\sqrt{10}$. This unitary transformation $H_4=H$ is real. Therefore, its inverse transformation is defined by the matrix being transpose to $H_4$, that is, $H_4^{-1}=H_4^{\prime }.$ We obtain $H_4^{\prime }|00⟩=|\mathit{x}⟩$. Using the property of the inverse operation (transposition), the matrix $H_4^{\prime }$ can be composed by the matrices of rotation by angles ${\{-\vartheta }_3,-{\vartheta }_2,-{\vartheta }_1\}$ as follows:

$$\begin{gathered} H_4^{\prime }={\left[R_{{\vartheta }_{3;\left(\mathrm{0,1}\right)}}{R}_{{\vartheta }_{2;\left(\mathrm{1,2}\right)}}{R}_{{\vartheta }_{1;\left(\mathrm{2,3}\right)}}\right]}^{\prime }={\left[R_{{\vartheta }_{1;\left(\mathrm{2,3}\right)}}\right]}^{\prime }{\left[R_{{\vartheta }_{2;\left(\mathrm{1,2}\right)}}\right]}^{\prime }{\left[R_{{\vartheta }_{3;\left(\mathrm{0,1}\right)}}\right]}^{\prime }=R_{{-\vartheta }_{1;\left(\mathrm{2,3}\right)}}{R}_{{-\vartheta }_{2;\left(\mathrm{1,2}\right)}}{R}_{-{\vartheta }_{3;\left(\mathrm{0,1}\right)}} \\ =\left(I_2\oplus R_{{-\vartheta }_1}\right)\left(P\left(I_2\oplus R_{{-\vartheta }_2}\right)P\prime \right)\left(R_{{-\vartheta }_3}\oplus I_2\right). \end{gathered}$$

To initiate the 2-qubit superposition $|\mathit{x}⟩$ from the first basis state $|00⟩$, the quantum circuit of the inverse 2-qubit DsiHT shown in Figure 7 can be used. The angular representation of the signal $\mathit{x}$ (or the superposition $|\mathit{x}⟩$) thus allows this superposition to be initiated by three rotations.

In the general case of $r>1$, the unitary matrix of the DsiHT generated by the real signal is real. Therefore, the inverse DsiHT is described by the matrix transpose to the matrix of the DsiHT and

$$H^{-1}\left({|0⟩}^{\otimes r}\right)=H^{\prime }{\left(\mathrm{1,0},\dots ,0\right)}^{\prime }={\left(x_0,x_1,\dots ,x_{2^r-1}\right)}^{\prime }=\mathit{x}\prime =|\mathit{x}⟩.$$

Example  2:

Consider the 8-D unit vector $\mathit{x}=\left(\mathrm{1,1},2,-2,-\mathrm{1,3},\mathrm{1,2}\right)/5$ and the corresponding 3-qubit superposition

$$|\mathit{x}⟩={\displaystyle \frac{1}{5}}\left(1\left|000⟩+1\right|001⟩+2\left|010⟩-2\right|011⟩-1\left|100⟩+3\right|101⟩+1\left|110⟩+2\right|111⟩\right).$$

For the 8-point strong DsiHT with this generator, the transform matrix can be written as

$$H_8=\left(\prod _{i=0}^6{R}_{{\vartheta }_{7-i;\left(i,i+1\right)}}\right)=\left(R_{{\vartheta }_{7;\left(\mathrm{0,1}\right)}}{R}_{{\vartheta }_{6;\left(\mathrm{1,2}\right)}}{R}_{{\vartheta }_{5;\left(\mathrm{2,3}\right)}}{R}_{{\vartheta }_{4;\left(\mathrm{3,4}\right)}}{R}_{{\vartheta }_{3;\left(\mathrm{4,5}\right)}}{R}_{{\vartheta }_{2;\left(\mathrm{5,6}\right)}}{R}_{{\vartheta }_{1;\left(\mathrm{6,7}\right)}}\right).$$

(16)

Seven rotations by the following angles (in degrees) are used:

$$\left\{{\vartheta }_k;k=1:7\right\}=\left\{-78.46°,-78.22°,65.35°,-62.68°,75.03°,-36.69°,-63.43°\right\}.$$

(17)

The transform in the matrix form can be written as

$$H_8{\mathit{x}}^{\prime }=D\left[\begin{array}{rrrrrrrr}1& 1& 2& -2& -1& 3& 1& 2\\ -24& 1& 2& -2& -1& 3& 1& 2\\ 0& -23& 2& -2& -1& 3& 1& 2\\ 0& 0& -9.5& -2& -1& 3& 1& 2\\ 0& 0& 0& 7.5& -1& 3& 1& 2\\ 0& 0& 0& 0& 14& 3& 1& 2\\ 0& 0& 0& 0& 0& -5/3& 1& 2\\ 0& 0& 0& 0& 0& 0& -2& 1\end{array}\right]\times {\displaystyle \frac{1}{5}}\left[\begin{array}{r}1\\ 1\\ 2\\ -2\\ -1\\ 3\\ 1\\ 2\end{array}\right]=\left[\begin{array}{c}1\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]=|000⟩.$$

(18)

Here, the diagonal matrix

$$D=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left\{0.2,0.0408,0.0426,-0.0957,-0.1185,0.0690,0.3586,0.4472\right\}.$$

The matrix of the inverse 8-point DsiHT is calculated by seven rotations in reverse order, as

$$H_8^{\prime }=R_{{-\vartheta }_{1;\left(\mathrm{6,7}\right)}}{R}_{{-\vartheta }_{2;\left(\mathrm{5,6}\right)}}{R}_{{-\vartheta }_{3;\left(\mathrm{4,5}\right)}}{R}_{{-\vartheta }_{4;\left(\mathrm{3,4}\right)}}{R}_{{-\vartheta }_{5;\left(\mathrm{2,3}\right)}}{R}_{{-\vartheta }_{6;\left(\mathrm{1,2}\right)}}{R}_{-{\vartheta }_{7;\left(\mathrm{0,1}\right)}}$$

The inverse transform of the first basis state $|000⟩$ can be written with the diagonal matrix $D$ as

$$H_8^{\prime }\left[\begin{array}{c}1\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]=\left[\begin{array}{rrrrrrrr}1& -24& 0& 0& 0& 0& 0& 0\\ 1& 1& -23& 0& 0& 0& 0& 0\\ 2& 2& 2& -9.5& 0& 0& 0& 0\\ -2& -2& -2& -2& 7.5& 0& 0& 0\\ -1& -1& -1& -1& -1& 14& 0& 0\\ 3& 3& 3& 3& 3& 3& -5/3& 0\\ 1& 1& 1& 1& 1& 1& 1& -2\\ 2& 2& 2& 2& 2& 2& 2& 1\end{array}\right]D\times \left[\begin{array}{c}1\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\\ 0\end{array}\right]={\displaystyle \frac{1}{5}}\left[\begin{array}{r}1\\ 1\\ 2\\ -2\\ -1\\ 3\\ 1\\ 2\end{array}\right]=|\mathit{x}⟩.$$

(19)

This diagonal matrix $D$ was introduced to better illustrate the process of constructing this matrix $H_8$. The first basis function of this transformation is the generator $\mathit{x}$ itself (the first column of the matrix). Each next basis function is result of motion of the previous basis functions (for mode detail, see [31]).

Any $r$-qubit superposition $|\mathit{x}⟩$ can be obtained from the zero state ${|0⟩}^{\otimes r}$ by the unitary inverse heap transform $H_{2^r}^{\prime }$ generated by $\mathit{x}.$ The DsiHT requires only $(2^r-1)$ rotations and it is defined by the angular representation of $\mathit{x}$. Not all of these Givens rotations act on adjacent bit planes, what is defined by the path of the transformation. However, there are paths that allow us to perform all calculations only on adjacent bit planes.

6. Path in the DsiHT

In this section, we discuss the path which is an important characteristic of the DsiHT. The path, or the order in which components of both generator and input signal are processed, can be chosen in different ways (for more detail, see [17]). Two more example paths are shown in Figure 8. They are paths with partitioning into pairs. The path shown in part (a) resembles the path that is used in the 8-point fast Fourier transform (FFT) with decimation-in-frequency [35]. Only the 8-point DsiHT uses 7 rotations (as butterflies) rather than $4\cdot 3=12$ butterflies in the FFT. These rotations $T_k$ are denoted by the circles with numbers $k$ in the figure. The angles ${\vartheta }_k$ of these rotations’ changes when the path of the DsiHT changes. We will call the path in part (a) fast path #1 (or simply fast path) and fast path #2 is the path in part (b). The small open circles in these diagrams are used to indicate that these transform outputs are zero when the input is the generator $\mathit{x}.$

The path can be defined in such a way that the computation complexity of the transformation and its inverse will be minimized. For instance, it is possible to define the path that allows for simplification of quantum circuits of the $r$-qubit QsiHT with only controlled gates of rotation. No permutations are required. The representation of the QsiHT matrix is still in the form shown in Eq. 11, only the order of processing the input data and therefore the angles of rotations must be changed accordingly. The $r$-qubit QsiHT always requires maximum ${(2}^r-1)$ rotations on different bit-planes (some of 2$\times$2-rotations may be trivial). In the case of the strong QsiHT, for many rotations, these bit planes do not differ by only one bit and require additional permutations, as shown in the $R_{{\vartheta }_2;\mathrm{1,2}}$ rotation example above. Other paths require rotations in other planes. Therefore, it is necessary to find better ways to simplify the calculations [36]. To illustrate this fact, we first consider the fast path which is shown in Figure 8(a), for the 8-point DsiHT, or 3-qubit QsiHT. Then, the example for 2-qubit QsiHT is presented and compared with Example 1.

The notation $R_{\vartheta ;(i,j)}$ will be used for the rotation $R_{\vartheta }$ by the angle $\vartheta$ in planes number $i$and $j$, where integer $j\ne i\in \{\mathrm{0,1},\dots ,r-1\}$. The 3-qubit DsiHT with the fast path is calculated by seven rotations in the following way. On stage 1, four rotations are used, as shown below with the corresponding bit-planes

$$R_{{\vartheta }_1;\left(\mathrm{0,1}\right)}\to \left(\begin{array}{c}000\\ 001\end{array}\right), R_{{\vartheta }_2;\left(\mathrm{2,3}\right)}\to \left(\begin{array}{c}010\\ 011\end{array}\right), R_{{\vartheta }_3;\left(\mathrm{4,5}\right)}\to \left(\begin{array}{c}100\\ 101\end{array}\right), R_{{\vartheta }_4;\left(\mathrm{6,7}\right)}\to \left(\begin{array}{c}110\\ 111\end{array}\right).$$

(20)

The circuit elements for these controlled rotation gates are shown on Figure 9.

The second stage uses two rotations, and the last stage uses one rotation, as shown below

$$\underbrace{R_{{\vartheta }_5;\left(\mathrm{0,2}\right)}\to \left(\begin{array}{c}000\\ 010\end{array}\right),R_{{\vartheta }_6;\left(\mathrm{4,6}\right)}\to \left(\begin{array}{c}100\\ 110\end{array}\right)}, \mathrm{a}\mathrm{n}\mathrm{d} \underbrace{R_{{\vartheta }_7;\left(\mathrm{0,4}\right)}\to \left(\begin{array}{c}000\\ 100\end{array}\right)}.$$

(21)

The circuit elements of these rotation gates are shown in Figure 10.

Summarizing the above reasonings, we obtain the quantum circuit for the 3-qubit QsiHT with the fast path which is shown in Figure 11. Seven controlled gates of rotations with 2 control qubits are used: no additional permutations. All rotations operate on adjacent bit planes, that is, on bit planes that differ by only one bit. Such bit planes are also called adjacent qubits.

Based on this scheme, the circuit for the inverse 3-qubit QsiHT can be easily composed and it is shown in Figure 12. The signs of all angles change, and seven controlled gates of rotation are executed in reverse order. Here, the input is a 3-qubit superposition $|\mathit{y}⟩$ and the output is $H_8^{\prime }|\mathit{y}⟩$. In the case when the input is the first basis state, $|\mathit{y}⟩=|000⟩$, the output of the circuit is the 3-qubit superposition, $|\mathit{x}⟩$. Thus, to initiate a 3-qubit superposition maximum of seven rotation gates are only required.

Now, let us analyze this circuit and compare it with the circuit of the strong 3-qubit DsiHT.

Example  3:

Consider the quantum circuit of the 3-qubit QsiHT with the fast path, which is given in Figure 13. The generator $\mathit{x}=\left(\mathrm{1,1},2,-2,-\mathrm{1,3},\mathrm{1,2}\right)$ and its 3-qubit $|\mathit{x}⟩$ is the same as in Example 2. The angular representation of the generator with the fast path equals to

$$A_{\mathit{x}}=\left\{{\vartheta }_1,{\vartheta }_2,{\dots ,\vartheta }_7\right\}=\left\{\underbrace{-45°,45°,71.5651°,-63.4349°},\underbrace{-63.4349°,35.2644°},\underbrace{50.7685°}\right\}.$$

The matrix of this transformation can be written as

$$H_{8;\mathrm{f}\mathrm{a}\mathrm{s}\mathrm{t}}=D\left[\begin{array}{rrrrrrrr}1& 1& 2& -2& -1& 3& 1& 2\\ -1& 1& 0& 0& 0& 0& 0& 0\\ -2& -2& 1& -1& 0& 0& 0& 0\\ 0& 0& 2& 2& 0& 0& 0& 0\\ 2& 2& 3& -3& 1& -3& -1& -2\\ 0& 0& 0& 0& 3& 1& 0& 0\\ 0& 0& 0& 0& 1& -3& 2& 4\\ 0& 0& 0& 0& 0& 0& -2& 1\end{array}\right]$$

(22)

where $D=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left\{0.2,0.7071,0.3162,0.3536,0.1633,0.3162,0.1826,0.4472\right\}.$This matrix has 32 zero coefficients compared to 21 zeros in the matrix $H_8$ for the strong DsiHT, which is given in Eq. 18.

Now we describe the quantum circuit for the 3-qubit strong DsiHT. This transform is composed of seven rotations

$$H_8=R_{{\vartheta }_{7;\left(\mathrm{0,1}\right)}}{R}_{{\vartheta }_{6;\left(\mathrm{1,2}\right)}}{R}_{{\vartheta }_{5;\left(\mathrm{2,3}\right)}}{R}_{{\vartheta }_{4;\left(\mathrm{3,4}\right)}}{R}_{{\vartheta }_{3;\left(\mathrm{4,5}\right)}}{R}_{{\vartheta }_{2;\left(\mathrm{5,6}\right)}}{R}_{{\vartheta }_{1;\left(\mathrm{6,7}\right)}}.$$

The angles ${\vartheta }_k;k=1:7$ are from the set in Eq. 17. The block-diagram of the transform is shown in Figure 14.

Four matrices from the above list, namely $R_{{\vartheta }_{7;\left(\mathrm{0,1}\right)}}$, $R_{{\vartheta }_{5;\left(\mathrm{2,3}\right)}}$, $R_{{\vartheta }_{3;\left(\mathrm{4,5}\right)}}$, and $R_{{\vartheta }_{1;\left(\mathrm{6,7}\right)}}$ are described by the controlled rotation gates similar to those shown in Figure 9, only with different angles ${\vartheta }_7, {\vartheta }_5$, ${\vartheta }_3$, and ${\vartheta }_1$ in parts (a)-(d), respectively. Three rotations by the angles ${\vartheta }_2,$ ${\vartheta }_4$, and ${\vartheta }_6$ operate on the nonadjacent bit-planes 5 and 6, 3 and 4, and 1 and 2, respectively. For example, in binary representation, the numbers $5=\left(101\right)$ and $6=\left(110\right)$ differ by two bits. Therefore, these three rotations should be fulfilled by the rotation gates on nearest-neighbor bit-planes. The circuit elements for other three rotation gates $R_{{\vartheta }_{2;\left(\mathrm{5,6}\right)}}$, $R_{{\vartheta }_{4;\left(\mathrm{3,4}\right)}}$, and $R_{{\vartheta }_{6;\left(\mathrm{1,2}\right)}}$ can be described as follows.

1)

By the permutation $\left(\mathrm{4,5}\right)$, the matrix of the rotation on bit-planes 5 and 6 is presented as

$$R_{{\vartheta }_{2;\left(\mathrm{5,6}\right)}}=\left[\begin{array}{ccc}{I}_4& & \\ & X& \\ & & I_2\end{array}\right]R_{{\vartheta }_{2;\left(\mathrm{4,6}\right)}}\left[\begin{array}{ccc}{I}_4& & \\ & X& \\ & & I_2\end{array}\right].$$

(23)

The permutation $\left(\mathrm{4,5}\right)$ is the controlled NOT, that is, $I_4\oplus X\oplus I_2.$ The controlled gate for the rotation on bit-planes 4 and 6 is shown in Figure 10 in part (b). Therefore, the operation $R_{{\vartheta }_{2;\left(\mathrm{5,6}\right)}}$ can be fulfilled by using the rotation gate $R_{{\vartheta }_{2;\left(\mathrm{4,6}\right)}}$. The circuit element of the operation in Eq. 19 is shown in Figure 15.

2)

By two permutations $\left(\mathrm{2,3}\right)$ and $\left(\mathrm{4,6}\right)$, the matrix of the rotation on bit-planes 3 and 4 is presented as

$$R_{{\vartheta }_{4;\left(\mathrm{3,4}\right)}}=\left[\begin{array}{ccc}{I}_2& & \\ & X& \\ & & P_{\mathrm{0,2}}\end{array}\right]R_{{\vartheta }_{4;\left(\mathrm{2,6}\right)}}\left[\begin{array}{ccc}{I}_2& & \\ & X& \\ & & P_{\mathrm{0,2}}\end{array}\right], P_{\mathrm{0,2}}=\left[\begin{array}{cccc}0& 0& 1& 0\\ 0& 1& 0& 0\\ 1& 0& 0& 0\\ 0& 0& 0& 1\end{array}\right].$$

(24)

Here, $P_{\mathrm{0,2}}$ is the permutation (0,2). The circuit element of this rotation is shown in Figure 16. The operation is reduced to the rotation gate on the adjacent bit-planes $2=\left(010\right)$ and $6=\left(110\right).$

3)

By the permutation $\left(\mathrm{0,1}\right)$, the matrix of the rotation on bit-planes 1 and 2 is presented as

$$R_{{\vartheta }_{6;\left(\mathrm{1,2}\right)}}=\left[\begin{array}{ccc}X& & \\ & I_2& \\ & & I_4\end{array}\right]R_{{\vartheta }_{6;\left(\mathrm{0,2}\right)}}\left[\begin{array}{ccc}X& & \\ & I_2& \\ & & I_4\end{array}\right].$$

(25)

This operation is fulfilled by using the rotation gate $R_{{\vartheta }_{6;\left(\mathrm{0,2}\right)}}$. This controlled gate of the rotation on adjacent bit-planes 0 and 2 is shown in Figure 10 in part (a). The permutation $\left(\mathrm{0,1}\right)$ is the controlled NOT, $X\oplus I_2\oplus I_4.$ The circuit element of the operation in Eq. 25 is shown in Figure 17.

The entire quantum circuit for the 3-qubit QsiHT is shown on Figure 18. Together with seven controlled rotation gates, eight controlled NOTs are used. All these gates are controlled by two bits (qubits).

The circuit for the inverse 3-qubit strong QsiHT is shown in Figure 19. Here, the rotation matrix $R_{{\vartheta }_k}^{\prime }$ denotes the inverse matrix to $R_{{\vartheta }_k}$, that is, $R_{{\vartheta }_k}^{\prime }=R_{{-\vartheta }_k}$, for $k=1:7$. There is some symmetry in this scheme (as well as in Figure 18). The change of each controlled NOT from one rotation operation to the next one is accomplished by adding $\left(\mathrm{2,2}\right)$ to the bit-plane numbers. Indeed, these changes are $\left(\mathrm{0,1}\right)\to \left(\mathrm{2,3}\right)\to \left(\mathrm{4,5}\right)\to \left(\mathrm{6,7}\right)$. Quantum circuits for $r$-qubit strong QsiHT can be described similarly for the number of qubits $r\ge 4.$

The above comparison of 3-qubit QsiHTs shows how important it is to choose a path that leads to a simplified circuit. The fast path for the 3-qubit QsiHT allows the effective quantum circuit to build. This circuit is without permutations. The results of simulation of this circuit in Qiskit for the 3-qubi superposition $|\mathit{x}⟩$ in Example 3 are shown in Table 2. For comparison, the results of using the circuit in Figure 13 with the QsiHT by the fast path are also given in this table. One can note from this table that the DsiHT with the fast path improves (reduces) the computational errors.

Example  4:

Now we consider the 4-point DsiHT and its inverse transform with the fast path shown in Figure 20. This is the two-stage path with partitioning into 3 pairs.

On these two stages of calculation, the rotations are described as

$$\underbrace{R_{{\vartheta }_1;\left(\mathrm{0,1}\right)}\to \left(\begin{array}{c}00\\ 01\end{array}\right), R_{{\vartheta }_2;\left(\mathrm{2,3}\right)}\to \left(\begin{array}{c}10\\ 11\end{array}\right)}, \mathrm{a}\mathrm{n}\mathrm{d} \underbrace{R_{{\vartheta }_3;\left(\mathrm{0,2}\right)}\to \left(\begin{array}{c}00\\ 10\end{array}\right)}.$$

(26)

The corresponding circuit elements of these controlled rotation gates are shown in Figure 21. Two gates are with 0-control bit and the gate for the second rotation is with 1-control bit. Connected together these gates compose the 2-qubit QsiHT.

It directly follows from this circuit that the inverse 2-qubit DsiHT with the fast path is described by the circuit given in Figure 22. Thus, this circuit can be used to initiate the 2-qubit superposition $|\mathit{x}⟩$ from the basis state $|00⟩$.

Compared with the quantum scheme in Figure 7, the above circuit does not require permutations, which makes this circuit more effective. The angular representations $A_{\mathit{x}}$ of the same signal $\mathit{x}$ are different, due to their different paths. For the fast path, this representation is $A_{\mathit{x}}=\left\{{\vartheta }_1,{\vartheta }_2,{\vartheta }_3\right\}=\left\{-63.4349°,-26.5651°,-45°\right\}.$ The matrix of the DsiHT is equal to

$$H_{4;\mathrm{f}\mathrm{a}\mathrm{s}\mathrm{t}}=\left[\begin{array}{rrrr}0.3162& 0.6325& 0.6325& 0.3162\\ -0.8944& 0.4472& 0& 0\\ -0.3162& -0.6325& 0.6325& 0.3162\\ 0& 0& -0.4472& 0.8944\end{array}\right]=D\left[\begin{array}{rrrr}1& 2& 2& 1\\ -2& 1& 0& 0\\ -1& -2& 2& 1\\ 0& 0& -1& 2\end{array}\right],$$

(27)

where the diagonal matrix $D=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left\{0.3162,0.4472,0.3162,0.4472\right\}.$ (from: example3_paper.m)

The new circuit for calculating this 2-qubit superposition $|\mathit{x}⟩=(|0⟩+2|1⟩+2|2⟩+|3⟩)/\sqrt{10}$ is given in Figure 23. The rotation by ${-\vartheta }_3=45°$ is $R_{45}=1/\sqrt{2}[1 -1; 1 1]$. Also, note that the 2-qubit $|\mathit{x}⟩$ is not the tensor product of qubits, that is, $|\mathit{x}⟩$ is an entangled superposition.

The results of simulation of this circuit in Qiskit are shown in Table 3. For comparison, the results of using the inverse QsiHT with strong path (with the circuit in Figure 7) are also given in this table.

We can also consider the transformation of two 2-qubits superpositions $|\mathit{y}⟩$ to $|\mathit{x}⟩$, by using two 2-qubit QsiHT.

Example  5:

Consider the same 2-qubit $|\mathit{x}⟩=(|0⟩+2|1⟩+2|2⟩+|3⟩)\prime /\sqrt{10}.$ Let the second superposition be $|\mathit{y}⟩={\left(\mathrm{1,3},-\mathrm{4,1}\right)}^{\prime }/\sqrt{27}.$ Denoting by the $H_y$ and $H_x$ the 2-qubit QsiHTs generated by $|\mathit{y}⟩$ and $|\mathit{x}⟩$, respectively, we obtain the 2-qubit $|\mathit{x}⟩$, by two steps: $1)H_y|\mathit{y}⟩=|00⟩, \mathrm{a}\mathrm{n}\mathrm{d} 2)H_x^{\prime }|00⟩=|\mathit{x}⟩.$ Calculating the angular representation of $|\mathit{y}⟩,$ $A_y=\left\{{\theta }_1,{\theta }_2,{\theta }_3\right\}=\{-71.5651°, 14.0362°, 52.5131°\}$, we obtain

$$|\mathit{x}⟩=H_x^{\prime }{H}_y|\mathit{y}⟩=\underbrace{R_{{-\vartheta }_1;\left(\mathrm{0,1}\right)}{R}_{-{\vartheta }_2;\left(\mathrm{2,3}\right)}{R}_{{-\vartheta }_3;\left(\mathrm{0,2}\right)}}\underbrace{R_{{\theta }_3;\left(\mathrm{0,2}\right)}{R}_{{\theta }_1;\left(\mathrm{2,3}\right)}{R}_{{\theta }_1;\left(\mathrm{0,1}\right)}}|\mathit{y}⟩,$$

(28)

or

$$|\mathit{x}⟩=\underbrace{R_{{-\vartheta }_1;\left(\mathrm{0,1}\right)}{R}_{-{\vartheta }_2;\left(\mathrm{2,3}\right)}}{R}_{{{\theta }_3-\vartheta }_3;\left(\mathrm{0,2}\right)}\underbrace{R_{{\theta }_1;\left(\mathrm{2,3}\right)}{R}_{{\theta }_1;\left(\mathrm{0,1}\right)}}|\mathit{y}⟩.$$

(29)

The quantum circuit for this 2-qubit transform with five controlled rotation gates is given in Figure 24.

From the circuits in Figures 11,12,18,19, and 24, it is not difficult to notice a recursive structure of the direct and inverse QsiHTs with the fast path and construct the quantum circuits, to initiate any $r$-qubit superposition by $(2^r-1)$ rotations, for $r\ge 2$. To show this, we describe the 4-qubit QsiHT.

Example  6:

Let $|\mathit{x}⟩$ be the 4-qubit superposition $x_0|0⟩+x_1|1⟩+\cdots +x_{15}|15⟩$ with real amplitudes $x_k, k=0:15$. The 4-qubit QsiHT by this generator $|\mathit{x}⟩$ can be calculated by the quantum circuit given in Figure 25. Fifteen angles are used; the angles of the angular representation of the generator, $A_{\mathit{x}}=\left\{{\vartheta }_k, k=1:15\right\}.$ To simplify the notations in this circuit, all rotations are denoted by $R_k=R_{{\vartheta }_k}.$ Four stages of calculations are shown with a total of 15 gates with 3 control bits. On the first stage, eight rotation gates and the corresponding bit-planes are

$$R_{{\vartheta }_1;\left(\mathrm{0,1}\right)}\to \left(\begin{array}{c}0000\\ 0001\end{array}\right), R_{{\vartheta }_2;\left(\mathrm{2,3}\right)}\to \left(\begin{array}{c}0010\\ 0011\end{array}\right), R_{{\vartheta }_3;\left(\mathrm{4,5}\right)}\to \left(\begin{array}{c}0100\\ 0101\end{array}\right), R_{{\vartheta }_4;\left(\mathrm{6,7}\right)}\to \left(\begin{array}{c}0110\\ 0111\end{array}\right),$$

$$R_{{\vartheta }_5;\left(\mathrm{8,9}\right)}\to \left(\begin{array}{c}1000\\ 1001\end{array}\right), R_{{\vartheta }_6;\left(\mathrm{10,11}\right)}\to \left(\begin{array}{c}1010\\ 1011\end{array}\right), R_{{\vartheta }_7;\left(\mathrm{12,13}\right)}\to \left(\begin{array}{c}1100\\ 1101\end{array}\right), R_{{\vartheta }_8;\left(\mathrm{14,15}\right)}\to \left(\begin{array}{c}1110\\ 1111\end{array}\right).$$

On the 2^nd stage, 4 rotation gates have the last control bit 0,

$$R_{{\vartheta }_9;\left(\mathrm{0,2}\right)}\to \left(\begin{array}{c}0000\\ 0010\end{array}\right), R_{{\vartheta }_{10};\left(\mathrm{4,6}\right)}\to \left(\begin{array}{c}0100\\ 0110\end{array}\right), R_{{\vartheta }_{11};\left(\mathrm{8,10}\right)}\to \left(\begin{array}{c}1000\\ 1010\end{array}\right), R_{{\vartheta }_{12};\left(\mathrm{12,14}\right)}\to \left(\begin{array}{c}1100\\ 1110\end{array}\right).$$

Two rotation gates with the last two 0-control bits are used on stage 3, and the last rotation with three control bits 0 on stage 4,

$$R_{{\vartheta }_{13};\left(\mathrm{0,4}\right)}\to \left(\begin{array}{c}0000\\ 0100\end{array}\right), R_{{\vartheta }_{14};\left(\mathrm{8,12}\right)}\to \left(\begin{array}{c}1000\\ 1100\end{array}\right), \mathrm{a}\mathrm{n}\mathrm{d} R_{{\vartheta }_{15};\left(\mathrm{0,8}\right)}\to \left(\begin{array}{c}0000\\ 1000\end{array}\right).$$

The circuit for the inverse 4-qubit DsiHT is obtained directly from the above circuit and is shown in Figure 26. All rotations are performed in reverse order and the notations $R_k^{\prime }$=$R_{{-\vartheta }_k}$ are used, for $k=1:15.$

Let us consider for example the 4-qubit superposition $|\mathit{x}⟩$ with amplitudes of the normalized vector

$$\mathit{x}=\left(1,-\mathrm{1,2},4,2,-\mathrm{1,2},6,\mathrm{1,1},4,-\mathrm{2,4},-\mathrm{1,1},2\right)/\sqrt{111}.$$

With the fast path, the angular representation of this vector is the set of angles (rounded to two decimals)

$$A_{\mathit{x}}=\left\{\begin{array}{c}\underbrace{45°,-63.43°,26.56°,-71.56°,-45°,26.56°,14.04°,-63.43°}, \\ \underbrace{-72.45°,-70.53°,-72.45°,-28.47°}, \underbrace{-55.04°,-45°}, \underbrace{-39.02°}\end{array}\right\}.$$

The results of simulation of the above circuit for the 4-qubit superposition $|\mathit{x}⟩$ are given in Table 4.

Now, we consider the circuit symbols used in [37], for the $k$-fold uniformly controlled 1-qubit gate, which is a full sequence of $k$-fold controlled gates of rotation. Figure 27 shows the equivalent circuit of Figure 25.

Here, the notation $R_I$ stands for the set of controlled rotations $R_1, R_2, \dots , R_8$. Also, $R_{II}$ stands for rotations $R_9, R_{10}, R_{11}, R_{12}$, and $R_{III}$ for two rotations $R_{13}, R_{14}$, and $R_{IV}=R_{15}$. This circuit does not use any permutations. With these notations for the $k$-fold uniformly controlled 1-qubit gates, the quantum scheme in Figure 26 can be presented, as shown in Figure 28.

If we combine these two quantum circuits, we obtain a simple circuit for preparing the state $|\mathit{x}⟩$ from another state $|\mathit{y}⟩$. This diagram is shown in Figure 29. Here, in order to not complicate the notations in circuit elements, two parts of the circuits are colored in distinct colors. The angles of rotations in both parts are different. The first part of this circuit is used the set of rotations with 15 angles of the angular representation $A_{\mathit{y}}$ of the vector $\mathit{y}.$ The second part of the citcuit is composed by the set of inverse rotations with the 15 angles of the angular representation $A_{\mathit{x}}$ of the vector $\mathit{x}.$ This circuits requires 30 angles of these two vectors, or only 15 angles of $A_{\mathit{x}}$, if the state $|y⟩=|0000⟩.$

Then, two controlled gates in the middle of the circuit can be united as one rotation. Therefore, the total number of rotation gates in this circuit with 4 qubit input equals to $g\left(4\right)=2\left(8+4+2+1\right)-1=2\times \left(2^4-1\right)-1=2\times 2^4-3=29.$In general case of $r$-qibits, these two numbers are equal to $g\left(r\right)=2\times 2^r-3.$

It should be noted that the path with partitioning into pairs, which is shown in part (b) in Figure 8, can also be used for the 8-point DsiHT. However, its implementation for the 3-qubit QsiHT is facing difficulties. The first four rotations $R_{{\vartheta }_k}, k=1:4,$ operate on inputs and outputs on different bit-planes. This stage is described by the matrix

$$H_{8;1}=\left[\begin{array}{rrrrrrrr}{\mathrm{c}}_1& 0& 0& 0& -{\mathrm{s}}_1& 0& 0& 0\\ {\mathrm{s}}_1& 0& 0& 0& {\mathrm{c}}_1& 0& 0& 0\\ 0& {\mathrm{c}}_2& 0& 0& 0& {-\mathrm{s}}_2& 0& 0\\ 0& {\mathrm{s}}_2& 0& 0& 0& {\mathrm{c}}_2& 0& 0\\ 0& 0& {\mathrm{c}}_3& 0& 0& 0& -{\mathrm{s}}_3& 0\\ 0& 0& {\mathrm{s}}_3& 0& 0& 0& {\mathrm{c}}_3& 0\\ 0& 0& 0& {\mathrm{c}}_4& 0& 0& 0& -{\mathrm{s}}_4\\ 0& 0& 0& {\mathrm{s}}_4& 0& 0& 0& {\mathrm{c}}_4\end{array}\right],$$

where $c_k=\cos {\vartheta }_k$ and $s_k=\sin {\vartheta }_k, k=1:4.$ Also, on stage 2, the 5^th rotation uses the inputs number 0 and 4, but outputs number 0 and 2. The 6^th rotation uses the inputs number 2 and 6, and outputs number 4 and 6. These rotations are not on the same pair of adjacent bit-planes. Thus, we assume that this path is not suitable for the 3-qubit QsiHT. However, there are other effective paths for the $r$-qubit QsiHT, which will be shown in the next section. The larger the number of qubits, the more such paths can be found.

7. The DsiHT in the General Form

In this section, we consider the concept of the DsiHT in the general definition [16] and show a way to accomplish the state-t-state transformation in a fast way. Namely, the state transformation will be calculated by using only one QsiHT. This is a one-step process and does not require the use of reverse DsiHT in the method shown above in Figure 29 for 4-qubit superpositions.

The DisHT by rotations in its original definition consists of the following 2$\times$2 elementary rotations [16,30]

$$T_{\vartheta }\left[\begin{array}{c}x\\ y\end{array}\right]=\left[\begin{array}{c}\cos \left(\vartheta \right)x-\sin \left(\vartheta \right)y\\ a\end{array}\right].$$

(30)

Here, the parameter $a$ is given and the angle is calculated by

$$\vartheta =\arccos \left({\displaystyle \frac{a}{\sqrt{x^2+y^2}}}\right)+\arctan \left({\displaystyle \frac{x}{y}}\right),\left(a^2\le x^2+y^2\right),$$

(31)

$$\left(\vartheta =\arcsin \left({\displaystyle \frac{a}{x}}\right),\mathrm{i}\mathrm{f}y=0\right).$$

Thus, the angle of rotation is defined from the angular equation $\sin \left(\vartheta \right)x+\cos \left(\vartheta \right)y=a.$ The $N$-point DsiHT generated by the vector $\mathit{x}$ is defined with a given vector of parameters $\mathit{a}=\{a_1,a_2,\dots ,a_{N-1}\}$. It is clear that in the case, when the vector parameter is zero, $\mathit{a}=\{\mathrm{0,0},\dots ,0\}$, the DsiHT is the transform described above in Section 1, Section 2, Section 3, Section 4, Section 5 and Section 6.

As an example, Figure 30 shows the calculation of the 8-point strong DsiHT with the given parameters $\{a_1,a_2,\dots ,a_7\}$. The value $x_0^{\left(7\right)}$ contains the information of the energy of the generator minus the number $‖\mathit{a}‖=\sqrt{a_1^2+a_2^2+\dots {+a}_7^2}$. It is assumed that $‖\mathit{a}‖\le ‖\mathit{x}‖$. Thus, $x_0^{\left(7\right)}=\pm \sqrt{{‖\mathit{x}‖}^2-{‖\mathit{a}‖}^2}$. When performing the DsiHT, the angles of rotations $T_k=T_{{\vartheta }_k }, k=1:7,$are calculated by Eq. 31.

The DsiHT with the vector-parameter is defined similarly for other paths. The vector-parameter changes the angular representation of the generator, that is, the set of angles of rotations. Figure 31 shows the diagrams for calculating the 4- and 8-point DsiHTs with the fast path in part (a) and (b) respectively.

Example 7 (2-qubit transformation)

Consider two 4-D vectors $\mathit{x}=(\mathrm{1,2},-\mathrm{1,3})$ and $\mathit{y}=(4,-\mathrm{1,1},2)$ after normalization by $‖\mathit{x}‖=\sqrt{15}$ and $‖\mathit{y}‖=\sqrt{22.}$ The corresponding 2-qubit superpositions are

$$|\mathit{x}⟩={\displaystyle \frac{1}{\sqrt{15}}}\left(\left|0⟩+2\right|1⟩-\left|2⟩+3\right|3\right)$$

and

$$|\mathit{y}⟩={\displaystyle \frac{1}{\sqrt{22}}}\left(4|0⟩-|1⟩+|2⟩+2|3⟩\right).$$

Both 2-qubits are entangled. The vector parameter for the DsiHT is taken to be equal to

$$\mathit{a}=\left(y\left(1\right),y\left(2\right),y\left(3\right)\right)={\displaystyle \frac{1}{\sqrt{22}}}\left(-\mathrm{1,1},2\right).$$

The norm of this vector equals to $‖\mathit{a}‖=\sqrt{1+1+4}/\sqrt{22}$ $=\sqrt{6}/\sqrt{22}$ =$\sqrt{3/11}$. The DsiHT signal-flow graph with the generator $\mathit{x}$ and this vector parameter is shown in Figure 32.

The matrix of the DsiHT is equal to

$$H_4=\left[\begin{array}{rrrr}-0.2984& -0.2664& 0.7012& -0.5901\\ 0.6661& -0.7459& 0& 0\\ 0.6836& 0.6105& 0.3061& -0.2576\\ 0& 0& 0.6439& 0.7651\end{array}\right], \det \left(H_4\right)=1.$$

(32)

The angular representation of the generator is equal to $A_{\mathit{x}}=\left\{{\vartheta }_1,{\vartheta }_2,{\vartheta }_3\right\}=\left\{138.2359°,40.0829°,-66.4201°\right\}.$ The transform of the generator equals to

$$H_4{\mathit{x}}^{\prime }={\displaystyle \frac{1}{\sqrt{22}}}\left(-4,-\mathrm{1,1},2\right)\prime .$$

The amplitude of the first state $|0⟩$ is $-4/\sqrt{22}$ not $4/\sqrt{22}$, as in the state$|\mathrm{y}⟩.$ Therefore, with the additional 0-controlled phase shift gate

$$P_{-}\left(\pi \right)=\left[\begin{array}{rr}-1& 0\\ 0& 1\end{array}\right],$$

(33)

we obtain the following matrix equation:

$${\left(P_{-}\left(\pi \right)\oplus I_2\right)H}_4{\mathit{x}}^{\prime }=\left[\begin{array}{cccc}-1& & & \\ & 1& & \\ & & 1& \\ & & & 1\end{array}\right]H_4{\mathit{x}}^{\prime }={\displaystyle \frac{1}{\sqrt{22}}}\left[\begin{array}{r}4\\ -1\\ 1\\ 2\end{array}\right]=\mathit{y}\mathbf{\text{'}}.$$

Thus, with three rotations $\{R_{{\vartheta }_1;\left(\mathrm{0,1}\right)},R_{{\vartheta }_2;\left(\mathrm{2,3}\right)},R_{{\vartheta }_3;\left(\mathrm{0,2}\right)}\}$ plus the 0-controlled phase shift gate, the following 2-qubit transformation of states holds:

$$|\mathit{x}⟩={\displaystyle \frac{1}{\sqrt{15}}}\left(\left|0⟩+2\right|1⟩-\left|2⟩+3\right|3\right)\to |\mathit{y}⟩={\displaystyle \frac{1}{\sqrt{22}}}\left(4|0⟩-|1⟩+|2⟩+2|3⟩\right).$$

The corresponding circuit for this 2-qubit transform is shown in Figure 33.

Table 5 shows the results of measurement of the 2-qubit $|\mathit{y}⟩$ by implementing this circuit in Qiskit.

The same circuit can be used for obtaining the 2-qubit $|\mathit{y}⟩$ from other 2-qubits. Only the angles of three rotation gates will be changed. For example, the transformation of 2-qubits

$$|{\mathit{x}}_1⟩={\displaystyle \frac{1}{\sqrt{46}}}\left(\left|0⟩+5\right|1⟩+2\left|2⟩+4\right|3⟩\right)\to |\mathit{y}⟩={\displaystyle \frac{1}{\sqrt{22}}}\left(4|0⟩-|1⟩+|2⟩+2|3⟩\right)$$

requires a new set of angles, $A_{{\mathit{x}}_1}=\left\{{\vartheta }_1,{\vartheta }_2,{\vartheta }_3\right\}=\left\{52.3485°,33.8996°,-21.7498°\right\}$ in the circuit in Figure 33.

We also can consider another path for such 2-qubit QsiHT, which is illustrated in Figure 34. The order of the last two rotations changes, and this results in a change in the rotation angles.

In this case, the matrix of the 2-qubit QsiHT transform is equal to

$$H_4=\left[\begin{array}{rrrr}-0.3870& -0.7662& 0.3925& -0.3303\\ -0.4614& 0.6426& 0.4680& -0.3939\\ -0.7983& 0& -0.4608& 0.3878\\ 0& 0& 0.6439& 0.7651\end{array}\right],\det \left(H_4\right)=1.$$

(34)

The angular representation of the generator for this transform is equal to $A_{\mathit{x}}=\left\{{\vartheta }_1,{\vartheta }_2,{\vartheta }_3\right\}=\left\{40.0829°,233.9712°,50.0125°\right\}.$ The corresponding circuit is shown in Figure 35.

It should be noted for comparison, that the existent estimation of rotation gates in number of ${\mu }_1\left(r\right)=2^{r+1}-2$ [14,15,37] gives us the number 6, for the $r=2$ case. However, some operations, such as global phase and normalization, were not considered, and these publications do not contain a single illustrative example of implementing the schemes $|\mathit{x}⟩\to |\mathit{y}⟩$. The circuits in Figure 33 and Figure 35 use one phase shift gate and only three controlled elementary rotations, that is, $\mu \left(r\right)=2^r-1$ and this estimation is valid for integers $r\ge 2.$

Example 8 (3-qubit preparation)

Consider two 8-D vectors $\mathit{x}=(\mathrm{1,2},-\mathrm{1,3},\mathrm{1,2},-3,-1)$ and $\mathit{y}=(8,-\mathrm{1,1},\mathrm{2,1},-\mathrm{2,3},1)$ after normalization them by $‖\mathit{x}‖=\sqrt{30}$ and $‖\mathit{y}‖=\sqrt{85,}$ respectively. These two vectors represent the 3-qubit superpositions $|\mathit{x}⟩$ and $|\mathit{y}⟩$,

$$|\mathit{x}⟩={\displaystyle \frac{1}{\sqrt{30}}}\left(\left|0⟩+2\right|1⟩-\left|2⟩+3\right|3⟩+\left|4⟩+2\right|5⟩-3\left|6⟩-\right|7⟩\right)$$

and

$$|\mathit{y}⟩={\displaystyle \frac{1}{\sqrt{85}}}\left(8|0⟩-|1⟩+\left|2⟩+2\right|3⟩+\left|4⟩-2\right|5⟩+3\left|6⟩+\right|7⟩\right).$$

Our goal is to obtain the vector $\mathit{y}$ from $\mathit{x}$ using only one 8-point DsiHT, or 3-qubit QsiHT. Therefore, the vector-parameter for the 8-point DsiHT with the generator $|\mathit{x}⟩$ is taken to be equal to

$$\mathit{a}=\left(y\left(1\right),y\left(2\right),y\left(3\right),\dots ,y\left(7\right)\right)={\displaystyle \frac{1}{\sqrt{85}}}\left(-\mathrm{1,1},\mathrm{2,1},-\mathrm{2,3},1\right).$$

The norm of this vector is equal to $‖\mathit{a}‖=\sqrt{85-64}/\sqrt{85}$ $=\sqrt{21/85}=0.4971$. The angular representation of the vector $\mathit{x}$ equals to

$$\begin{array}{l}{A}_{\mathit{x}}=\left\{{\vartheta }_k;k=1:7\right\}\\ =\left\{\underbrace{131.9727°,49.4953°,148.6630°,-29.2633°},\underbrace{-63.0605°,-87.9558°},\underbrace{-48.5902°}\right\}.\end{array}$$

The matrix of the 8-point DsiHT has 32 zero coefficients and is equal to

$$H_8=\left[\begin{array}{rrrrrrrr}-0.2004& -0.2228& 0.3830& -0.4484& -0.0229& -0.0139& 0.6539& 0.3664\\ 0.7435& -0.6688& 0& 0& 0& 0& 0& 0\\ 0.5962& 0.6628& 0.2943& -0.3445& 0& 0& 0& 0\\ 0& 0& 0.7604& 0.6495& 0& 0& 0& 0\\ 0.2272& 0.2526& -0.4343& 0.5084& -0.0202& -0.0123& 0.5767& 0.3231\\ 0& 0& 0& 0& 0.5201& -0.8541& 0& 0\\ 0& 0& 0& 0& 0.8536& 0.5197& 0.0311& 0.0174\\ 0& 0& 0& 0& 0& 0& -0.4888& 0.8724\end{array}\right].$$

The transform of the generator $\mathit{x}$ equals to

$$H_8{\mathit{x}}^{\prime }={\displaystyle \frac{1}{\sqrt{85}}}\left(-8,-\mathrm{1,1},\mathrm{2,1},-\mathrm{2,3},1\right)\prime .$$

Here, the sign of the first component should be changed, because $y\left(0\right)=8.$ As in Example 5, we can use the phase shift gate $P_{-}\left(\pi \right)$ with two control qubits, that is $P_{-}\left(\pi \right)\oplus I_6$, and reconstruct the sign of $y\left(0\right).$We obtain the given vector $\mathit{y}$,

$${\left(P_{-}\left(\pi \right)\oplus I_6\right)H}_8{\mathit{x}}^{\prime }=\left[\begin{array}{cccc}-1& & & \\ & 1& & \\ & & I_2& \\ & & & I_4\end{array}\right]H_8{\mathit{x}}^{\prime }={\displaystyle \frac{1}{\sqrt{85}}}\left[\begin{array}{r}8\\ -1\\ 1\\ 2\\ 1\\ -2\\ 3\\ 1\end{array}\right]=\mathit{y}\mathbf{\text{'}}.$$

Thus, with seven rotations plus the 0-controlled phase shift gate, the following 3-qubit transformation holds: $|\mathit{x}⟩\to |\mathit{y}⟩.$ The circuit for this transformation is shown in Figure 36. In this figure, the rotations in three stages are defined by

$$R_I=\left\{R_{{\vartheta }_1;\left(\mathrm{0,1}\right)},R_{{\vartheta }_2;\left(\mathrm{2,3}\right)},R_{{\vartheta }_3;\left(\mathrm{4,5}\right)},R_{{\vartheta }_4;\left(\mathrm{6,7}\right)}\right\}, R_{II}=\left\{R_{{\vartheta }_5;\left(\mathrm{0,2}\right)},R_{{\vartheta }_6;\left(\mathrm{4,6}\right)}\}, R_{III}=\{R_{{\vartheta }_7;\left(\mathrm{0,4}\right)}\right\}.$$

The results of simulation of this circuit in Qiskit are shown in Table 6.

When comparing with the known estimation ${\mu }_1\left(r\right)=2^{r+1}-2$ for the $r=3$ case, we obtain the number 14, respectively. The circuit of $|\mathit{x}⟩\to |\mathit{y}⟩$ transformation in Figure 34 includes seven elementary rotations, that is, $\mu \left(3\right)=7.$The quantum circuits for the 4-qubit state-to-state transformation, which are given in Figure 25, Figure 26, Figure 27 and Figure 28, only angles ${\vartheta }_k, k=1:7,$ should be changed according to the vector parameter $\mathit{a}=\left(y\left(1\right),y\left(2\right),\dots ,y\left(7\right)\right)$. This vector changes the angular representation of $A_{\mathit{x}}$ of the generator. Similarly, with the complete set of $k$-fold uniformly controlled 1-qubit gates on $n$ qubits, the quantum circuit can be described for the $r$-qubit state-to-state transformation, when $r\ge 2.$

8. Conclusion

In this work, we describe in detail a new method of state-to-state transformation with only elementary rotations, by using the concept of the DsiHT in its general form. This transformation is generated by the given superposition and the set of angular equations with the given parameters. The analogue of this transform, namely the QsiHT, is described, and the quantum circuits of its implementation were considered. Amplitudes of quantum superpositions are considered real. On the examples with the DsiHT with strong and fast paths, the importance of this characteristic for the transform is shown. The QsiHT with the fast path allows us to build the circuit for $r$-qubit state-to-state transformation with only ${(2}^r-1)$ elementary rotations, each with only one angle. The number of such possible paths increases with the number $r$ of qubits and they can be found, and the best ones can be selected for a specific implementation of the proposed method. The problem of state preparation with complex amplitudes can also be considered and solved by the complex QsiHT [32,36].

It is important to mention the following. We present a “permutation-free method” for transforming one quantum superposition into another. This is our vision for building quantum circuits for multi-qubit gates. If there is a desire (or need) to modify the proposed quantum circuits to circuits with (not controlled) rotations, then the following should be considered. As shown in [14], the $(r-1)$-fold uniformly controlled 1-qubit gates on $n$-qubits can be implemented by $2^{r-1}$ CNOTs controlled by one qubit. Therefore, no more than $2^{r-1}+2^{r-2}+\dots +2+1=2^r-1$ one-qubit CNOTs are required, or seven CNOTs, for the realization of thee 3-qubit QsiHT in the circuit in Figure 36. One controlled 1-qubit phase shift gate is needed only to change the sign of the first amplitude in the transformed superposition $|\mathit{y}⟩$. This is the trivial operation of the sign change, $y_0\to -y_0$, on the first bit. Note that the probability $p_0={|y_0|}^2$ of observing the first state $|0⟩$ in $|\mathit{y}⟩$ does not require this operation. The best-known number of CNOTs is ${2(2}^r-r-1)$ [14,37], or eight for the 3-qubit case, and 22 and 52 CNOTs for the 4- and 5-qubit cases, respectively. For 4-and 5-qubit state-to-state transformations, the modified circuits will require no more than $2^4-1=15$ and $2^5-1=31$ CNOTs, respectively.