# Quantum Computing

**text to come**

<!-- Equation labels as ordinary links -->
<div id="_auto1"></div>

$$
\begin{equation}
   \Psi(x,t)=\Psi(x,t)_a+\Psi(x,t)_b,
\label{_auto1} \tag{1}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto2"></div>

$$
\begin{equation}
    |q\rangle = \alpha \vert 0 \rangle + \beta |1\rangle,
\label{_auto2} \tag{2}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto3"></div>

$$
\begin{equation}
    |q\rangle = \frac{1}{\sqrt{2}}(\vert 0 \rangle + |1\rangle),
\label{_auto3} \tag{3}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:dataord"></div>

$$
\begin{equation}
    |\psi\rangle = a_0|\psi\rangle_0+a_1|\psi\rangle_1+a_2|\psi\rangle_2+a_3|\psi\rangle_3+a_4|\psi\rangle_4+a_5|\psi\rangle_5+a_6|\psi\rangle_6+a_7|\psi\rangle_7,
\label{eq:dataord} \tag{4}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto4"></div>

$$
\begin{equation}
    \hat{H}\vert 0 \rangle=
    \frac{1}{\sqrt{2}}(\vert 0 \rangle +|1\rangle),
\label{_auto4} \tag{5}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto5"></div>

$$
\begin{equation}
    \hat{H}|1\rangle=
    \frac{1}{\sqrt{2}}(\vert 0 \rangle -|1\rangle).
\label{_auto5} \tag{6}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto6"></div>

$$
\begin{equation}
    \hat{H}\hat{H}\vert 0 \rangle=
    \hat{H}\frac{1}{\sqrt{2}}(\vert 0 \rangle +|1\rangle)
    =\vert 0 \rangle,
\label{_auto6} \tag{7}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto7"></div>

$$
\begin{equation}
    \hat{H}\hat{H}|1\rangle=
    \hat{H}\frac{1}{\sqrt{2}}(\vert 0 \rangle -|1\rangle)
    =|1\rangle.
\label{_auto7} \tag{8}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto8"></div>

$$
\begin{equation}
   \vert 0 \rangle=\left(\begin{array}{c} 1 \\ 0\end{array}   \right),
\label{_auto8} \tag{9}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto9"></div>

$$
\begin{equation}
   |1\rangle=\left(\begin{array}{c} 0 \\ 1\end{array}   \right),
\label{_auto9} \tag{10}
\end{equation}
$$

og

<!-- Equation labels as ordinary links -->
<div id="_auto10"></div>

$$
\begin{equation}
   \hat{H}=\frac{1}{\sqrt{2}}
   \left(\begin{array}{cc} 1 & 1 \\ 1& -1\end{array}   \right).
\label{_auto10} \tag{11}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto11"></div>

$$
\begin{equation}
   \hat{H}_{\mathrm{NOT}}=   
   \left(\begin{array}{cc} 0 & 1 \\ 1& 0\end{array}   \right).
\label{_auto11} \tag{12}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto12"></div>

$$
\begin{equation}
   \hat{H}_{\mathrm{NOT}}\vert 0 \rangle=   
   \left(\begin{array}{cc} 0 & 1 \\ 1& 0\end{array}\right)
   \left(\begin{array}{c} 1 \\ 0\end{array}   \right)=
   \left(\begin{array}{c} 0 \\ 1\end{array}   \right)=|1\rangle,
\label{_auto12} \tag{13}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto13"></div>

$$
\begin{equation}
   \hat{H}_{\mathrm{NOT}}|1\rangle=   
   \left(\begin{array}{cc} 0 & 1 \\ 1& 0\end{array}\right)
   \left(\begin{array}{c} 0 \\ 1\end{array}   \right)=
   \left(\begin{array}{c} 1 \\ 0\end{array}   \right)=\vert 0 \rangle.
\label{_auto13} \tag{14}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto14"></div>

$$
\begin{equation}
   \hat{H}_{\Phi}\vert 0 \rangle=e^{i\phi}\vert 0 \rangle,
\label{_auto14} \tag{15}
\end{equation}
$$

og

<!-- Equation labels as ordinary links -->
<div id="_auto15"></div>

$$
\begin{equation}
   \hat{H}_{\Phi}|1\rangle=|1\rangle,
\label{_auto15} \tag{16}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto16"></div>

$$
\begin{equation}
   \hat{H}_{\Phi}=   
   \left(\begin{array}{cc} e^{i\phi} & 0 \\ 0& 1\end{array}   \right).
\label{_auto16} \tag{17}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto17"></div>

$$
\begin{equation}
    |q\rangle_{in} = \alpha \vert 0 \rangle_{in} + \beta |1\rangle_{in},
\label{_auto17} \tag{18}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto18"></div>

$$
\begin{equation}
    |q\rangle_{ut}=\hat{H}|q\rangle_{in} = 
    \frac{1}{\sqrt{2}}((\alpha+\beta)\vert 0 \rangle_{ut} + 
    (\alpha-\beta)|1\rangle_{ut}),
\label{_auto18} \tag{19}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto19"></div>

$$
\begin{equation} 
   |q\rangle_{ut}=\hat{H}\hat{H}|q\rangle_{in}=|q\rangle_{in}.
\label{_auto19} \tag{20}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto20"></div>

$$
\begin{equation} 
   \hat{H}\hat{\Phi}\hat{H}\vert 0 \rangle,
\label{_auto20} \tag{21}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto21"></div>

$$
\begin{equation}
    \frac{1}{\sqrt{2}}((e^{i\phi}+1)\vert 0 \rangle + 
    (e^{i\phi}-1)|1\rangle).
\label{_auto21} \tag{22}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:superent"></div>

$$
\begin{equation}
   \frac{\vert 0 \rangle_1 |1\rangle_2 + |1\rangle_1\vert 0 \rangle_2}{\sqrt{2}},
\label{eq:superent} \tag{23}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:superent1"></div>

$$
\begin{equation}
   \frac{|+\rangle_1 |-\rangle_2 - |-\rangle_1|+\rangle_2}{\sqrt{2}},
\label{eq:superent1} \tag{24}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto22"></div>

$$
\begin{equation}
    |q\rangle_{1} = \alpha_1 \vert 0 \rangle_{1} + \beta_1 |1\rangle_{1},
\label{_auto22} \tag{25}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto23"></div>

$$
\begin{equation}
     |q\rangle_{12}= |q\rangle_{1} |q\rangle_{2}
\label{_auto23} \tag{26}
\end{equation}
$$

eller

<!-- Equation labels as ordinary links -->
<div id="_auto24"></div>

$$
\begin{equation}
   |q\rangle_{12}=\alpha_1 \alpha_2\vert 0 \rangle_1\vert 0 \rangle_2 +
                  \alpha_1 \beta_2\vert 0 \rangle_1|1\rangle_2 +
                  \beta_1 \alpha_2|1\rangle_1\vert 0 \rangle_2 +
                  \beta_1 \beta_2|1\rangle_1|1\rangle_2,
\label{_auto24} \tag{27}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto25"></div>

$$
\begin{equation} 
   |q\rangle_{ut}=\hat{H}|q\rangle_{in},
\label{_auto25} \tag{28}
\end{equation}
$$

$\hat{H}^{\dagger}=\hat{H}$ og $\hat{H}^{\dagger}\hat{H}=1$,

<!-- Equation labels as ordinary links -->
<div id="_auto26"></div>

$$
\begin{equation} 
   \hat{H}^{\dagger}|q\rangle_{ut}=\hat{H}^{\dagger}\hat{H}|q\rangle_{in}=
   |q\rangle_{in}.
\label{_auto26} \tag{29}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto27"></div>

$$
\begin{equation}
   \vert 0 \rangle_{c}  \vert 0 \rangle_{t}\rightarrow
   \vert 0 \rangle_{c}  \vert 0 \rangle_{t},
\label{_auto27} \tag{30}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto28"></div>

$$
\begin{equation}
   \vert 0 \rangle_{c}  |1\rangle_{t}\rightarrow
   \vert 0 \rangle_{c}  |1\rangle_{t},
\label{_auto28} \tag{31}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto29"></div>

$$
\begin{equation}
   |1\rangle_{c}  \vert 0 \rangle_{t}\rightarrow
   |1\rangle_{c}  |1\rangle_{t},
\label{_auto29} \tag{32}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto30"></div>

$$
\begin{equation}
   |1\rangle_{c}  |1\rangle_{t}\rightarrow
   |1\rangle_{c}  \vert 0 \rangle_{t}.
\label{_auto30} \tag{33}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto31"></div>

$$
\begin{equation}
   \hat{H}_{\mathrm{CNOT}}=
    \left(\begin{array}{cccc} 1 & 0 & 0 &0  \\ 
                              0 & 1 & 0 &0  \\
                              0& 0 & 0 &1  \\
                              0 & 0 & 1 &0\end{array}\right),
\label{_auto31} \tag{34}
\end{equation}
$$

$\vert 0 \rangle_c\vert 0 \rangle_t$, $\vert 0 \rangle_c|1\rangle_t$, $|1\rangle_c\vert 0 \rangle_t$ 
$|1\rangle_c|1\rangle_t$.

<!-- Equation labels as ordinary links -->
<div id="_auto32"></div>

$$
\begin{equation}
   \vert 0 \rangle_c\vert 0 \rangle_t=\left(\begin{array}{c} 1  \\ 
                                                 0  \\
                                                 0  \\
                                                 0  \end{array}\right),
\label{_auto32} \tag{35}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto33"></div>

$$
\begin{equation}
   \vert 0 \rangle_c|1\rangle_t=\left(\begin{array}{c} 0  \\ 
                                                 1  \\
                                                 0  \\
                                                 0  \end{array}\right),
\label{_auto33} \tag{36}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto34"></div>

$$
\begin{equation}
   |1\rangle_c\vert 0 \rangle_t=\left(\begin{array}{c} 0  \\ 
                                                 0  \\
                                                 1  \\
                                                 0  \end{array}\right),
\label{_auto34} \tag{37}
\end{equation}
$$

og

<!-- Equation labels as ordinary links -->
<div id="_auto35"></div>

$$
\begin{equation}
   |1\rangle_c|1\rangle_t=\left(\begin{array}{c} 0  \\ 
                                                 0  \\
                                                 0  \\
                                                 1  \end{array}\right).
\label{_auto35} \tag{38}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto36"></div>

$$
\begin{equation}
   \hat{H}_{\mathrm{CNOT}}|1\rangle_c|1\rangle_t=
    \left(\begin{array}{cccc} 1 & 0 & 0 &0  \\ 
                              0 & 1 & 0 &0  \\
                              0& 0 & 0 &1  \\
                              0 & 0 & 1 &0\end{array}\right)
                              \left(\begin{array}{c} 0  \\ 
                                                 0  \\
                                                 0  \\
                                                 1  \end{array}\right)=
\left(\begin{array}{c} 0  \\ 
                                                 0  \\
                                                 1  \\
                                                 0  \end{array}\right),
\label{_auto36} \tag{39}
\end{equation}
$$

$|1\rangle_c\vert 0 \rangle_t$, dvs. target-biten forandrer verdi fra bit '1'
til bit '0' når kontroll-biten har verdi bit '1'.




## Hamiltonians

A general two-body Hamiltonian for fermionic system 
can be written as

<!-- Equation labels as ordinary links -->
<div id="eq:twobodyH"></div>

$$
\begin{equation}
\label{eq:twobodyH} \tag{40}
H = E_0 + \sum_{ij=1} E_{ij} a^\dag_i a_j
+\sum_{ijkl = 1} V_{ijkl} a^\dag_i a^\dag_j a_l a_k,
\end{equation}
$$

where $E_0$ is a constant energy term, $E_{ij}$ represent all the
one-particle terms, allowing for non-diagonal terms as well. The
one-body term can represent 
a chosen single-particle potential, the kinetic energy or other more
specialized terms such as  
those discussed in connection with the Hubbard model or the pairing Hamiltonian  
discussed below.
The two-body interaction part is given by $V_{ijkl}$ and can be any
two-body interaction, from 
Coulomb interaction to the interaction between nucleons.  
The sums run over all possible single-particle levels $N$. 
Note that
this model includes particle numbers from zero to the number of
available quantum levels, $n$. To simulate states with fixed numbers
of fermions one would have to either rewrite the Hamiltonian or
generate specialized input states in the simulation.

The algorithm which we will develop in this section and in
However, 
in our demonstrations of the quantum computing algorithm, we will limit ourselves to
two simple models, which however capture much of the important physics
in quantum mechanical 
many-body systems. We will also limit ourselves to spin $j=1/2$
systems, although our algorithm 
can also simulate  higher $j$-values, such as those which occur in nuclear, atomic and
molecular physics, it simply uses one qubit for every available
quantum state.   
These simple models are the Hubbard model and a pairing
Hamiltonian.
We start with the spin $1/2$ Hubbard model, described by the following Hamiltonian

<!-- Equation labels as ordinary links -->
<div id="_auto37"></div>

$$
\begin{equation}
H_H = \epsilon \sum_{i, \sigma} a_{i\sigma}^\dag a_{i\sigma} 
-t \sum_{i, \sigma} \left(a^\dag_{i+1, \sigma}a_{i, \sigma}
+a^\dag_{i, \sigma}a_{i+1, \sigma} \right) \notag 
\label{_auto37} \tag{41}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:hubbard"></div>

$$
\begin{equation} 
 + U \sum_{i=1} a_{i+}^\dag a_{i-}^\dag a_{i-}a_{i+},
\label{eq:hubbard} \tag{42}
\end{equation}
$$

where $a^{\dagger}$ and $a$ are fermion creation and annihilation operators, respectively.
This is a chain of sites where each site has room for one spin up
fermion and one spin down fermion. 
The number of sites is $N$, and the sums over $\sigma$ are sums over
spin up and down only.
Each site has a single-particle
energy $\epsilon$. There is a repulsive term $U$ if there is a pair
of particles at the same site. It is energetically favourable to tunnel to
neighbouring sites, described by 
the hopping terms with coupling constant $-t$.

The second model-Hamiltonian is the simple pairing Hamiltonian

<!-- Equation labels as ordinary links -->
<div id="eq:pairing"></div>

$$
\begin{equation}
   H_P=\sum_i \varepsilon_i a^{\dagger}_i a_i -\frac{1}{2} g\sum_{ij>0}
           a^{\dagger}_{i}
     a^{\dagger}_{\bar{\imath}}a_{\bar{\jmath}}a_{j},
\label{eq:pairing} \tag{43}
\end{equation}
$$

The indices $i$ and $j$ run over the number of levels $N$, and the label $\bar{\imath}$ 
stands for a time-reversed state. The parameter $g$ is the strength of the pairing force 
while $\varepsilon_i$ is the single-particle energy of level $i$. 
In our case
we assume that the single-particle levels are equidistant (or
degenerate) with a fixed spacing $d$. 
Moreover, in our simple model, the degeneracy of the single-particle
levels is set to $2j+1=2$, with $j=1/2$  
being the spin of the particle. This gives a set of single-particle
states with the same spin projections as 
for the Hubbard model.  Whereas in the Hubbard model we operate with
different sites with  
spin up or spin down particles, our pairing models deals thus with
levels  with double degeneracy. 
Introducing the pair-creation operator 
$S^+_i=a^{\dagger}_{im}a^{\dagger}_{i-m}$,
one can rewrite the Hamiltonian in 
Eq. ([43](#eq:pairing)) as

$$
H_P=d\sum_iiN_i+
     \frac{1}{2} G\sum_{ij>0}S^+_iS^-_j,
$$

where  $N_i=a^{\dagger}_i a_i$
is the number operator, and 
$\varepsilon_i = id$ so that the single-particle orbitals 
are equally spaced at intervals $d$. The latter commutes with the 
Hamiltonian $H$. In this model, quantum numbers like seniority 
$\cal{S}$ are good quantum numbers, and the eigenvalue problem 
can be rewritten in terms of blocks with good seniority. 
Loosely 
speaking, the seniority quantum number $\cal{S}$ is equal to 
the number of unpaired particles.
Furthermore, in  a series of papers, Richardson
obtained the exact solution of the pairing Hamiltonian, with 
semi-analytic (since there is still the need for a numerical solution) 
expressions for the eigenvalues and eigenvectors. The exact solutions
have had important consequences for several fields, from Bose condensates to
nuclear superconductivity and is currently a very active field of studies, see for example
Finally, for particle numbers up to $P \sim 20$, the above model can be 
solved exactly through numerical diagonalization and one can obtain all eigenvalues.
It serves therefore also as an excellent ground for comparison with our algorithm based
on models from quantum computing.


## Basic quantum  gates

Benioff showed that one could make a quantum mechanical Turing machine
by using various  unitary operations on a quantum system.
Benioff  demonstrated 
that a quantum computer can calculate anything a
classical computer can. To do this one needs a quantum system and
basic operations that can approximate all unitary operations
on the chosen many-body system. We describe in this subsection the basic ingredients entering 
our algorithms.

### Qubits, gates and circuits

In this article we will use the standard model of quantum information,
where
the basic unit of information is the qubit, the quantum bit. 
As mentioned in the introduction, any
suitable 
two-level quantum system can be a qubit, 
it is the smallest system there is with the
least complex dynamics.
Qubits are both abstract measures of information and physical objects.
Actual physical qubits can be ions trapped in magnetic fields where
lasers can access only two energy levels or  the nuclear spins of some of
the atoms in molecules accessed and manipulated by an NMR machine.
Several other ideas have been proposed and some tested.

The computational basis for one qubit is ${\ensuremath{\vert 0 \rangle}}$ (representing for example bit $0$) 
for the first state
and ${\ensuremath{|1\rangle}}$ (representing bit $1$) for the second, and for a set of qubits  
the tensor products of
these basis states for each qubit form a product basis. Below we write out the different
basis states for a system of $n$ qubits.

<!-- Equation labels as ordinary links -->
<div id="eq:compBasis"></div>

$$
\begin{equation}
\label{eq:compBasis} \tag{44}
{\ensuremath{\vert 0 \rangle}} \equiv {\ensuremath{\vert 00\cdots 0 \rangle}} =
      {\ensuremath{\vert 0 \rangle}} \otimes {\ensuremath{| 0 \rangle}} \otimes
          \cdots 
\otimes {\ensuremath{\vert 0 \rangle}} 
\notag 
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto38"></div>

$$
\begin{equation} 
{\ensuremath{|1\rangle}} \equiv {\ensuremath{\vert 00\cdots 1\rangle}} =
    {\ensuremath{\vert 0 \rangle}} \otimes {\ensuremath{| 0 \rangle}} \otimes
        \cdots 
\otimes {\ensuremath{|1\rangle}} 
\notag 
\label{_auto38} \tag{45}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto39"></div>

$$
\begin{equation} 
\vdots \notag 
\label{_auto39} \tag{46}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto40"></div>

$$
\begin{equation} 
{\ensuremath{|2^n-1\rangle}} \equiv {\ensuremath{|11\cdots 1\rangle}} =
    {\ensuremath{|1\rangle}} \otimes {\ensuremath{| 1\rangle}} \otimes 
\cdots \otimes {\ensuremath{|1 \rangle}}.
\notag 
\label{_auto40} \tag{47}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto41"></div>

$$
\begin{equation} 
\label{_auto41} \tag{48}
\end{equation}
$$

This is a $2^n$-dimensional system and we number the different basis
states using binary numbers corresponding to the order in which they appear in the
tensor product.

Quantum computing means to  manipulate and measure  qubits in such a
way that the results from a measurement yield the solutions to  a given problem. 
The quantum operations we need to be able to perform our simulations are 
a small set of elementary single-qubit
operations, or single-qubit gates, and one universal two-qubit gate,
in our case the so-called CNOT gate defined below.

To represent quantum computer algorithms graphically we use circuit
diagrams. In a circuit diagram each qubit is represented by a line,
and operations on the different qubits are represented by boxes.

## Number of work qubits versus number of simulation qubits

The largest possible amount of different eigenvalues is $2^s$, where
$s$ is the number of simulation qubits. The resolution in the energy
spectrum we get from measuring upon the work qubits is $2^w$, with $w$ the number of
work qubits.
Therefore the resolution per eigenvalue in a non-degenerate system is
$2^{w-s}$. The higher the degeneracy the less work qubits are needed. 




## Number of operations

Counting the number of single-qubit and $\sigma_z\sigma_z$ operations
for different sizes of systems simulated gives us an indication of the
decoherence time needed for different physical realizations of a
quantum simulator or computer. The decoherence time is an average time
in which the state of the qubits will be destroyed by noise, also called 
decoherence, while the operation time is the average time an operation takes
to perform on the given system. Their fraction is the number of
operations possible to perform before decoherence destroys the
computation. In table  we have listed the number of
gates used for the pairing model, $H_P$, and the Hubbard model, $H_H$,
for different number of simulation qubits. 

<table border="1">
<thead>
<tr><th align="center">     </th> <th align="center">$s=2$</th> <th align="center">$s=4$</th> <th align="center">$s=6$</th> <th align="center">$s=8$</th> <th align="center">$s=10$</th> <th align="center">$s=12$</th> </tr>
</thead>
<tbody>
<tr><td align="left">   $H_P$    </td> <td align="center">   9        </td> <td align="center">   119      </td> <td align="center">   333      </td> <td align="center">   651      </td> <td align="center">   1073      </td> <td align="center">   1598      </td> </tr>
<tr><td align="left">   $H_H$    </td> <td align="center">   9        </td> <td align="center">   51       </td> <td align="center">   93       </td> <td align="center">   135      </td> <td align="center">   177       </td> <td align="center">   219       </td> </tr>
</tbody>
</table>
Number of two-qubit gates used in simulating the time
  evolution operator of the pairing model, $H_P$, and the Hubbard
  model, $H_H$, for different number of simulation qubits $s$.



We list here some useful relations involving different $\sigma$ matrices,

<!-- Equation labels as ordinary links -->
<div id="_auto42"></div>

$$
\begin{equation}
\sigma_x \sigma_z = -i\sigma_y, \quad
\sigma_z \sigma_x = i\sigma_y, \quad [\sigma_x, \sigma_z]=-2i\sigma_y,
\label{_auto42} \tag{49}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto43"></div>

$$
\begin{equation}
\sigma_x \sigma_y = i\sigma_z, \quad
\sigma_y \sigma_x = -i\sigma_z, \quad [\sigma_x, \sigma_y]=2i\sigma_z,
\label{_auto43} \tag{50}
\end{equation}
$$

and

<!-- Equation labels as ordinary links -->
<div id="_auto44"></div>

$$
\begin{equation}
\sigma_y \sigma_z = i\sigma_x, \quad
\sigma_z \sigma_y = -i\sigma_x, \quad [\sigma_y, \sigma_z]=2i\sigma_x.
\label{_auto44} \tag{51}
\end{equation}
$$

For any two non-equal $\sigma$-matrices $a$ and $b$ we have

<!-- Equation labels as ordinary links -->
<div id="_auto45"></div>

$$
\begin{equation}
aba = -b.
\label{_auto45} \tag{52}
\end{equation}
$$

The Hermitian $\sigma$-matrices $\sigma_x$, $\sigma_y$ and $\sigma_z$
result in the identity matrix when squared

<!-- Equation labels as ordinary links -->
<div id="_auto46"></div>

$$
\begin{equation}
\sigma_x^2 = _1_,\quad 
\sigma_y^2 = _1_,\quad 
\sigma_z^2 = _1_,\quad 
\label{_auto46} \tag{53}
\end{equation}
$$

which can be used to obtain  simplified expressions for exponential functions involving $\sigma$-matrices

<!-- Equation labels as ordinary links -->
<div id="_auto47"></div>

$$
\begin{equation}
e^{\pm i\alpha \sigma}=\cos(\alpha) _1_ \pm i \sin(\alpha) \sigma. 
\label{_auto47} \tag{54}
\end{equation}
$$

The equations we list below are necessary for the relation between  a general unitary
transformation on a set of qubits with  a product of two-qubit unitary
transformations. We have the general equation for $a,b \in \{\sigma_x,\sigma_y, \sigma_z\}$, where $a\neq b$.

<!-- Equation labels as ordinary links -->
<div id="_auto48"></div>

$$
\begin{equation}
  e^{-i\pi/4a} b e^{i\pi/4a} = \frac{1}{2} (_1_ -ia) b ( _1_ + ia)
  \notag
\label{_auto48} \tag{55}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto49"></div>

$$
\begin{equation}  
  =  \frac{1}{2} (b + aba + i[b,a]) \notag
\label{_auto49} \tag{56}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto50"></div>

$$
\begin{equation} 
  = \frac{i}{2}[b,a].
\label{_auto50} \tag{57}
\end{equation}
$$

The more specialized equations read

<!-- Equation labels as ordinary links -->
<div id="eq:rotations1"></div>

$$
\begin{equation}
\label{eq:rotations1} \tag{58}
  e^{-i\pi/4 \sigma_x} \sigma_z e^{i\pi/4 \sigma_x} = -\sigma_y, 
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:rotations2"></div>

$$
\begin{equation} 
\label{eq:rotations2} \tag{59}
  e^{-i\pi/4 \sigma_y} \sigma_z e^{i\pi/4 \sigma_y} = \sigma_x, 
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:rotations3"></div>

$$
\begin{equation} 
\label{eq:rotations3} \tag{60}
  e^{-i\pi/4 \sigma_z} \sigma_x e^{i\pi/4 \sigma_z} = \sigma_y, 
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="eq:rotations4"></div>

$$
\begin{equation} 
\label{eq:rotations4} \tag{61}
  e^{-i\pi/4 \sigma_z} \sigma_y e^{i\pi/4 \sigma_z} = -\sigma_x. 
\end{equation}
$$

We need also different products of the operator$\sigma_z$  with the  raising and lowering operators

<!-- Equation labels as ordinary links -->
<div id="eq:pmzs"></div>

$$
\begin{equation}
\label{eq:pmzs} \tag{62}
  \sigma_+ \sigma_z = -\sigma_+ 
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto51"></div>

$$
\begin{equation} 
  \sigma_z \sigma_+ = \sigma_+, 
\label{_auto51} \tag{63}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto52"></div>

$$
\begin{equation} 
  \sigma_- \sigma_z = \sigma_-, 
\label{_auto52} \tag{64}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto53"></div>

$$
\begin{equation} 
  \sigma_z \sigma_- = -\sigma_-. 
\label{_auto53} \tag{65}
\end{equation}
$$

<!-- Equation labels as ordinary links -->
<div id="_auto54"></div>

$$
\begin{equation} 
\label{_auto54} \tag{66}
\end{equation}
$$

9d5939f749ebbcf368e9f8769de7eba8ce9e7579
ipynb
Jupyter Notebook
doc/LectureNotes/_build/jupyter_execute/basicquantumcomputing.ipynb
Schoyen/ComputationalPhysics2
9cf10ffb2557cc73c4e6bab060d53690ee39426f
CC0-1.0
2015-01-21T08:29:56.000Z
2022-03-28T07:11:53.000Z
doc/LectureNotes/_build/jupyter_execute/basicquantumcomputing.ipynb
Schoyen/ComputationalPhysics2
9cf10ffb2557cc73c4e6bab060d53690ee39426f
CC0-1.0
2020-01-18T10:43:38.000Z
2020-02-08T13:15:42.000Z
doc/LectureNotes/_build/jupyter_execute/basicquantumcomputing.ipynb
Schoyen/ComputationalPhysics2
9cf10ffb2557cc73c4e6bab060d53690ee39426f
CC0-1.0
2015-02-09T10:02:00.000Z
2022-03-07T10:44:14.000Z
Qwen/Qwen-72B
1. YES
2. YES