$C^{12}-C^{13}$ Quantum computer with longitudinal rf-magnetic field

TL;DR

This paper demonstrates how to implement fundamental quantum operations and teleportation in a diamond structure with a linear chain of $C^{13}$ atoms using a specific magnetic field configuration.

Contribution

It introduces a novel method for quantum computing in diamond structures utilizing longitudinal rf-magnetic fields for key quantum gates and teleportation.

Findings

Successful realization of single-qubit rotations.

Implementation of CNOT and CCNOT gates.

Feasibility of quantum teleportation in the described system.

Abstract

It is shown that in the diamond $C^{12}$ structure with a linear chain of $C^{13}$ atoms which is inside a transverse static field, having a gradient along the linear chain and an rf-magnetic field in a plane with a component in the direction of the static field, one qubit rotation, the Controlled-Not (CNOT), the Controlled-Controlled-NoT (CCNOT) quantum gates, and teleportation algorithm.

Equations

$$U=\frac{{\mu }_0}{4\pi }\frac{3({\mathbf{m}}_1\cdot \mathbf{x})({\mathbf{m}}_2\cdot \mathbf{x})-{\mathbf{m}}_1\cdot {\mathbf{m}}_2}{|\mathbf{x}{|}^3},$$

$${\mathbf{m}}_i=\gamma {\mathbf{S}}_i,i=1,2.$$

$$U=\frac{2J}{\hslash }{S}_1^z{S}_2^z,J=\frac{{\mu }_0{\gamma }^2\hslash }{4\pi a^3},$$

$${\bf B}(x,t)=\bigl{(}0,b\cos\theta,-b\cos\theta+B_{0}(x)\bigr{)},\quad\quad\theta=\omega t+\varphi,$$

$$U=-\sum _{i=1}^n{\mathbf{m}}_i\cdot {\mathbf{B}}_i(t),{\mathbf{B}}_i(t)=\mathbf{B}(x_i,t),$$

$$\hat{H}=U+\frac{2J}{\hslash }\sum _{i=1}^{n-1}{S}_i^z{S}_{i+1}^z+\frac{2J^{\prime }}{\hslash }\sum _{i=1}^{n-2}{S}_i^z{S}_{i+2}^z.$$

$$\hat{H}={\hat{H}}_0+\hat{W},$$

$${\hat{H}}_0=-\sum _{i=1}^n{\omega }_i{S}_i^z+\frac{2J}{\hslash }\sum _{i=1}^{n-1}{S}_i^z{S}_{i+1}^z+\frac{2j^{\prime }}{\hslash }\sum _{i=1}^{n-2}{S}_i^z{S}_{i+2}^z$$

$$\widehat{W}=\Omega\sum_{k}^{n}\biggl{\{}\frac{1}{2i}(S_{k}^{+}-S_{k}^{-})\cos\theta+S_{k}^{z}\sin\theta\biggr{\}},$$

$$S_i^z|{\xi }_n,\dots ,{\xi }_i,\dots ,{\xi }_1⟩=\frac{(-1{)}^{{\xi }_i}\hslash }{2}|{\xi }_n,\dots ,{\xi }_i,\dots ,{\xi }_1⟩,$$

$$S_i^{+}|{\xi }_n,\dots ,{\xi }_i,\dots ,{\xi }_1⟩=\hslash {\delta }_{{\xi }_i,0}|{\xi }_n,\dots ,{\xi }_i+1,\dots ,{\xi }_1⟩,$$

$$S_i^{-}|{\xi }_n,\dots ,{\xi }_i,\dots ,{\xi }_1⟩=\hslash {\delta }_{{\xi }_i,1}|{\xi }_n,\dots ,{\xi }_i-1,\dots ,{\xi }_1⟩,$$

$${\hat{H}}_0|\xi ⟩=E_{\xi }|\xi ⟩$$

$$E_{\xi }=-\frac{\hslash }{2}\sum _{i=1}^n{\omega }_i+\frac{\hslash }{2}J\sum _{i=1}^{n-1}(-1{)}^{{\xi }_i+{\xi }_{i+1}}+\frac{\hslash }{2}{J}^{\prime }\sum _{i=1}^{n+2}(-1{)}^{{\xi }_i+{\xi }_{i+2}},$$

$$i\hslash \frac{\partial |\Psi ⟩}{\partial t}=\hat{H}|\Psi ⟩,$$

$$|\Psi (t)⟩=\sum _{\eta =1}^{2^n}{e}^{-i{E}_{\eta }t/\hslash }{D}_{\eta }(t)|\eta ⟩,$$

$$i{\dot{D}}_{\xi }(t)=-i\frac{\Omega \cos \theta }{2}\sum _{\eta =1}^{2^n}\sum _{k=1}^n{e}^{i{\omega }_{\xi \eta }t}{D}_{\eta }(t)⟨\xi |S_k^{+}-S_k^{-}|\eta ⟩+\frac{\Omega (-1{)}^{{\xi }_k}}{2}{D}_{\xi }(t)\sin \theta ,$$

$${\omega }_{\xi \eta }=\frac{E_{\xi }-E_{\eta }}{\hslash }.$$

$$\tau ={\omega }_{0}t,{\omega }_0=2\pi kHz.$$

$$1-qubit:{\omega }_1=100,\Omega =0.1$$

$$2-qubits:{\omega }_1=100,{\omega }_2=200,\Omega =0.1,J=2$$

$$3-qubits:{\omega }_1=100,{\omega }_2=200,{\Omega }_3=300,\Omega =0.1,J=2,J^{\prime }=1/4,$$

$$1-qubit:D_1(0)=1,D_2(0)=0$$

$$2-qubits:D_1(0)=\frac{1}{\sqrt{16}},D_2(0)=\sqrt{\frac{2}{16}},D_3(0)=\sqrt{\frac{5}{16}},D_4(0)=\sqrt{\frac{8}{16}}$$

$$3-qubits:D_1(0)=D_2(0)=\sqrt{\frac{1}{8}},D_3(0)=D_4(0)=0,D_5(0)=\sqrt{\frac{3}{8}},D_7(0)=\sqrt{\frac{2}{8}},D_8(0)=\frac{1}{\sqrt{8}}$$

$$\displaystyle\frac{1}{\sqrt{2}}\bigl{(}|0\rangle+|1\rangle\bigr{)}\quad\quad\hbox{(1-qubit)}$$

$$|\psi ⟩=c_1|0⟩+c_2|1⟩,c_1=\sqrt{\frac{3}{8}},c_2=\sqrt{\frac{5}{8}},$$

$$|{\Psi }_0⟩=|\psi ⟩\otimes |0⟩\otimes |0⟩=c_1|000⟩+c_2|100⟩.$$

$$T=\bigl{(}H_{3}\cdot CN_{32}\bigr{)}\cdot\bigl{(}CN_{21}\cdot H_{2}\bigr{)},$$

$$T|\Psi_{0}\rangle=\frac{1}{2}\biggl{\{}|00\rangle\otimes|\psi\rangle+|01\rangle\otimes N_{3}|\psi\rangle+|10\rangle\otimes Z_{3}|\psi\rangle+|11\rangle\otimes N_{3}Z_{3}|\psi\rangle\biggr{\}},$$

$$N_3|\psi ⟩=c_1|⟩+c_2|0⟩,Z_3|\psi ⟩=c_1|0⟩-c_1|1⟩,N_3{Z}_3|\psi ⟩=c_2|0⟩-c_1|1⟩.$$

$$R_{ij}(\tau ),$$

Full text

Keywords: Quantum computer, diamond, universal quantum gates

PACS: 03.67.-a, 03.67.Ac

1 Introduction

Despite of the great experimental and theoretical efforts, quantum computer realization for 1000 qubits remains being the goal of computer science and physical science. The two main problem faced on this goal are technological limitations [(for example number of controlled ions trapped, maximum number of signal that can be register in NMR system), and the decoherence due to interaction of any quantum system with the environment [1, 2, 3, 4]. Several quantum systems have been proposed for quantum computation [5, 6, 7, 8, 9, 10, 11], an some of the attractive (due to their high decoherence time [12]) are nuclear spin systems. One resent proposal of this type [13] is the $C^{12}-C^{13}$ diamond structure for a possible quantum computer, where a chain of $C^{13}$ stable isotopes of spin one half are embedded periodically along a given line of the diamond lattice. A transverse static field with some gradient along this line, and transversally to it, and rf-magnetic field are applied. The universal CNOT and CCNOT quantum gates are obtained with a single $\pi$-pulses (evolution time is $\tau=\pi/\Omega$, being $\Omega$ the Rabi’s frequency), and Hadamar’s gate is obtained with a $\pi/2$-pulse (evolution time of the quantum system is $\tau=\pi/2\Omega$). In this paper, we show that even with an rf-magnetic field having a component along the direction of the static magnetic field, it is possible to obtain the CNOT, CCNOT universal quantum gates with a $2\pi$-pulse, the Hadamar’s gate with a $\pi$-pulse. In addition, the implementation of teleportation algorithm using three qubits is presented with its value of goodness (different of Fidelity [14]) of the reproduction of the algorithm and its associated Boltzmann-Shanon’s entropy [15, 16].

2 Analytical Approach

The diamond structure $C^{12}$ (spin zero) with the line of isotopes $C^{13}$ (spin one half) along the x-axis is shown in Figure 1. The spin-spin interaction is due to their dipole magnetic moments interaction [17]

[TABLE]

where $|{\bf x}|$ is the distance separation between the magnetic moments ${\bf m}_{1}$ and ${\bf m}_{2}$ which are related with the nuclear spin as

[TABLE]

being $\gamma$ the gyromagnetic ratio ($\gamma\approx 2.675\times 10^{8}rad/T\cdot s$). Assuming Ising interaction, this energy can be written as

[TABLE]

where $\hbar$ is the Planck’s constant divided by $2\pi$ ($\hbar\approx 1,054571818\times 10^{-34}J\cdot s$) and $a=|{\bf x}|\sim 10^{-10}m$ is the $C^{13}-C^{13}$ nuclear separation in the diamond crystal. The propose magnetic field is

[TABLE]

where $\omega$, $\varphi$ and $b$ are the frequency, the phase and the amplitude of the rf-magnetic field, $B_{0}(x)$ is the static magnetic field transverse to the line of $C^{13}$’s. The interaction of the $C^{13}$’s with the magnetic field is given by[17]

[TABLE]

where $x_{i}$ is the position of the ith-atoms $C^{13}$. The Hamiltonian of the system considering Ising interaction at first and second neighbor is

[TABLE]

After rearranging terms, this Hamiltonian can be written as

[TABLE]

where the operators $\widehat{H}_{0}$ and $\widehat{W}$ are defined as

[TABLE]

and

[TABLE]

being $\omega_{i}=\gamma B_{0}(x_{i})$ the Larmore’s frequency associated to the ith-atoms $C^{13}$, $\Omega=\gamma b$ is the so called Rabi’s frequency, and $S_{k}^{+}$ and $S_{k}^{-}$ are the ascend and descend operators, $S_{k}^{\pm}=S_{k}^{x}\pm S_{k}^{y}$. The n-qubits registers $\{|\xi\rangle=|\xi_{n},\dots,\xi_{1}\rangle\}$ ($\xi_{i}=0,1$ for digital notation, and $\xi=1,\dots,2^{n}$ for decimal notation) form the basis for the Hilbert space of dimensionality $2^{n}$, our qubit is just the spin one half of the nucleus of the atom $C^{13}$. Now, having the known operations

[TABLE]

[TABLE]

and

[TABLE]

being $\delta_{ij}$ the Kronecker’s delta. The solution of the eigenvalue problems

[TABLE]

is given by the eigenvalues

[TABLE]

and the eigenfunctions are just $\{|\xi\rangle\}$, for $\xi=1,\dots,2^{n}$. Therefore, proposing a solution of the Schrödinger’s equation

[TABLE]

of the form

[TABLE]

we get the equation for the coefficients $D_{\xi}(t)$ as

[TABLE]

where the frequency $\omega_{\xi\eta}$ has been defined as

[TABLE]

We can make all system without units by defining the new time evolution as

[TABLE]

Then the constant $\Omega/\omega_{0}$, $J/\omega_{0}$ and $\omega_{i}/\omega_{0}$ become without units, but the equation (17) would be exactly of the same form. So, from now on, we will talk about the time evolution in terms of the variable (19), and all the constant must be though as its value times $2\pi kHz$.

3 Quantum Universal Gates

Equations (17) are solved numerically using Runge-Kutta method at fourth order. We made the simulation of the NOT quantum gate (1-qubit), Controlled-Not gate (2-qubits) and Controlled-Controlled-Not gate (3-qubits) to show the realization of a quantum computer with this magnetic field configuration in the $C^{12}-C^{13}$ diamond structure. The parameters used in our simulations are given by

[TABLE]

and the initial conditions are

[TABLE]

Figure 2 shows the probabilities $|D_{i}(\tau)|^{2}$ for the system to be in the state $|i\rangle$ as a function of $\tau$. This curves show the realization of the quantum gates NOT (transition: $|0\rangle\longleftrightarrow|1\rangle$), CNOT (transition: $|10\rangle\longleftrightarrow|11\rangle$), and CCNOT (transition: $|110\rangle\longleftrightarrow|111\rangle$) using a $2\pi$-pulse ($\tau=2\pi/\Omega$). The probabilities of the other states (CNOT and CCNOT) remain constant within and error of the order $10^{-5}$. In addition, We must observe that a superposition of states (Hadamar quantum gate) is gotten for a time evolution of $\tau=\pi/\Omega$ ($\pi$-pulse), that is the states of the form

[TABLE]

where some phases could be as factors on each state. We want to mention that if the initial conditions for CNOT and CCNOT are $D_{3}(0)=1$ and $D_{7}(0)=1$ and all other coefficients are zero, the superposition of states $|10\rangle+|11\rangle$ and $|110\rangle+|111\rangle$ are obtained too.

4 Realization of teleportation algorithm

Quantum teleportation [18] is one of the most fundamental element in quantum computation and quantum information which deserves its realization for any quantum computer proposal. Therefore, this algorithm will be studied in this section with three qubits configuration, where the first state will represent Bob, the second state will be Alice and the third state will be the teleported state. From the Alice location, we want to transport this state (in general, an unknown state even by Alice) to Bob place. If the state we want to transport is

[TABLE]

the initial state of our 3-qubits system is then

[TABLE]

The ”ideal” (mathematical) teleportation algorithm is represented by quantum circuit shown on Figure 3 (going from left to right),

or written in term of operators we have

[TABLE]

where $CN_{ij}$ is the CNOT operation between the qubit ”i” and the cubit ”j” ($CN_{ij}|\xi_{n},\dots,\xi_{i},\dots,\xi_{j},\dots,\xi_{1}\rangle=|\xi_{n},\dots,\xi_{i},\dots,\xi_{j}\oplus\xi_{i},\dots,\xi_{1}\rangle$), and $H$ is the Hadamar operator (1-qubit). The first term, \bigl{(}CN_{21}\cdot H_{2}\bigr{)}, makes the entanglement between Alice(second qubit) and Bob (first qubit), and the term \bigl{(}H_{3}\cdot CN_{32}\bigr{)} represents the Alice operations to produce the teleportation. The result of the application of the operator (25) to the initial state (24) can be written as

[TABLE]

where $N_{3}$ (NOT gate), $Z_{3}$ and $N_{3}Z_{3}$ act in the following way

[TABLE]

Thus, the state $|\psi\rangle$ has been teleported to Bob’s place, and this has occurred instantly (independently how separated are Alice and Bob). Of course, Alice needs to notify (at the speed of light or less) to Bob about the operation he needs to do in order to get the state $|\psi\rangle$ ($N_{3}^{2}=Z_{3}^{2}=I$, the identity operator). Notice then that at the end of the algorithm, we have four states with a probability $|c_{1}|^{2}/4$ and four states with probability $|c_{2}|^{2}/4$. However, in our non-ideal system, we need to apply the correct type of pulse with the desired frequency transition in order to get the state what we want. Let us denote by

[TABLE]

the pulse of length $\tau$ done in the system with a resonant frequency $\omega=(E_{i}-E_{j})/\hbar$. Our non-ideal teleportation operator is given by

[TABLE]

which represents $\pi$ and $2\pi$ pulses where the Rabi’s frequency has been omitted. The last square parenthesis represents the entanglement, and the other one represents the teleportation. The result of this application to the initial state (24) can be seen Figure 4, where we have gotten four states with the probability $3/32$ and four states with probability $5/32$.

The goodness of the reproducibility of the teleportation algorithm can be measured by comparing the probabilities of the states between the real case with the ideal case at the end of the algorithm, and this can be done through the quantity

[TABLE]

which goes from zero (excellent reproducibility) to one (non reproducibility at all). In our case, the value obtained is $G=7.4\times 10^{-4}$, meaning that we have a very good reproduction of the teleportation algorithm.

On the other hand, it is of interest to see how the information of the system is lost during the teleportation algorithm since starting with 2 states we finish with 8 states. Therefore, there must be an increasing on Boltzmann-Shanon entropy and its average value,

[TABLE]

Of course, the maximum value this entropy can have is went all the probabilities are the same (1/8), $S_{max}=\ln 8\approx 2.08$. Figure 5 shows the evolution of the entropy and its average value, and as we expected, there is an increasing on the entropy and its average value due to the increasing of the number of states involved in the quantum dynamics of teleportation algorithm. The oscillations shown on this figure are due to CNOT and CCNOT operations which change the states but does increases the number of states. The rapidly increasing on the entropy is due to Hadamar’s operators which increases the number of states involved in the quantum dynamics.

5 Conclusions and Comments

We have shown that a solid state quantum computer can be realized using $C^{12}-C^{13}$ structure in a diamond with even a longitudinal rf-magnetic field (one of its component is along the static field). The universal CNOT and CCNOT quantum gates and the one qubit rotation were well stablished on the simulations with multiples of a $\pi$-pulse. The teleportation algorithm was also implemented using these quantum gates in a 3-qubits configuration and using 10 pulses, and its associated Boltzmann-Shanon’s entropy was determined. We do not think decoherence [19, 20, 21, 22, 23, 24, 25, 26] could be of much consideration on this system. However, to make the escalation to many qubits system and the read out [27] the final computation result are things to need addressed in the near future.