^{1}

^{*}

^{1}

^{*}

I
n
the diamond
C^{12}
(nuclear spin zero) structure with a linear chain of
C^{13}
(nuclear spin one half) atoms is applied
to
a transverse static field with respect
of
this linear chain, having a gradient along the linear chain, and it is also applied
to
an rf-magnetic field in a plane with a component in the direction of the static field. It is shown that one qubit rotation, the Controlled-Not (CNOT), the Controlled-Controlled-NoT (CCNOT) quantum gates, and teleportation algorithm can be implemented on this structure using integer multiples of electromagnetic π-pulses. Therefore, a quantum computer can be constructed in this form.

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 problems faced on this goal are technological limitations (for example, the number of controlled ions trapped, or the maximum number of signal that can be registered by the NMR system), and the decoherence due to interaction of any quantum system with the environment [^{12} - C^{13} diamond structure for a possible quantum computer, where a chain of C^{13} stable isotopes of spin one half is 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 π-pulses (evolution time is τ = π / Ω , being Ω the Rabi’s frequency), and Hadamar’s gate is obtained with a π/2-pulse (evolution time of the quantum system is τ = π / 2 Ω ). 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π-pulse, the Hadamar’s gate with a π-pulse. In addition, the implementation of teleportation algorithm using three qubits is presented with its value of goodness (different of Fidelity [

The diamond structure C^{12} (spin zero) with the line of isotopes C^{13} (spin one half) along the x-axis is shown in

U = μ 0 4 π 3 ( m 1 ⋅ x ) ( m 2 ⋅ x ) − m 1 ⋅ m 2 | x | 3 , (1)

where | x | is the distance separation between the magnetic moments m 1 and m 2 which are related with the nuclear spin as

m i = γ S i , i = 1 , 2. (2)

being γ the gyromagnetic ratio ( γ ≈ 2.675 × 10 8 rad / T ⋅ s ). Assuming Ising interaction, this energy can be written as

U = 2 J ℏ S 1 z S 2 z , J = μ 0 γ 2 ℏ 4 π a 3 , (3)

where ℏ is the Planck’s constant divided by 2π ( ℏ ≈ 1.054571818 × 10 − 34 J ⋅ s ) and a = | x | ~ 10 − 10 m is the C^{13} - C^{13} nuclear separation in the diamond crystal. The propose magnetic field is

B ( x , t ) = ( 0, b c o s θ , − b c o s θ + B 0 ( x ) ) , θ = ω t + φ , (4)

where ω , φ 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 [

U = − ∑ i = 1 n m i ⋅ B i ( t ) , B i ( t ) = B ( x i , t ) , (5)

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

H ^ = U + 2 J ℏ ∑ i = 1 n − 1 S i z S i + 1 z + 2 J ′ ℏ ∑ i = 1 n − 2 S i z S i + 2 z . (6)

After rearranging terms, this Hamiltonian can be written as

H ^ = H ^ 0 + W ^ , (7)

where the operators H ^ 0 and W ^ are defined as

H ^ 0 = − ∑ i = 1 n ω i S i z + 2 J ℏ ∑ i = 1 n − 1 S i z S i + 1 z + 2 j ′ ℏ ∑ i = 1 n − 2 S i z S i + 2 z (8)

and

W ^ = Ω ∑ k n { 1 2 i ( S k + − S k − ) cos θ + S k z sin θ } , (9)

being ω i = γ B 0 ( x i ) the Larmore’s frequency associated to the ith-atoms C^{13}, Ω = γ b is the so called Rabi’s frequency, and S k + and S k − are the ascend and descend operators, S k ± = S k x ± S k y . The n-qubits registers { | ξ 〉 = | ξ n , ⋯ , ξ 1 〉 } ( ξ i = 0 , 1 for digital notation, and ξ = 1 , ⋯ , 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

S i z | ξ n , ⋯ , ξ i , ⋯ , ξ 1 〉 = ( − 1 ) ξ i ℏ 2 | ξ n , ⋯ , ξ i , ⋯ , ξ 1 〉 , (10)

S i + | ξ n , ⋯ , ξ i , ⋯ , ξ 1 〉 = ℏ δ ξ i ,0 | ξ n , ⋯ , ξ i + 1, ⋯ , ξ 1 〉 , (11)

and

S i − | ξ n , ⋯ , ξ i , ⋯ , ξ 1 〉 = ℏ δ ξ i ,1 | ξ n , ⋯ , ξ i − 1, ⋯ , ξ 1 〉 , (12)

being δ i j the Kronecker’s delta. The solution of the eigenvalue problems

H ^ 0 | ξ 〉 = E ξ | ξ 〉 (13)

is given by the eigenvalues

E ξ = − ℏ 2 ∑ i = 1 n ω i + ℏ 2 J ∑ i = 1 n − 1 ( − 1 ) ξ i + ξ i + 1 + ℏ 2 J ′ ∑ i = 1 n + 2 ( − 1 ) ξ i + ξ i + 2 , (14)

and the eigenfunctions are just { | ξ 〉 } , for ξ = 1 , ⋯ , 2 n . Therefore, proposing a solution of the Schrödinger’s equation

i ℏ ∂ | Ψ 〉 ∂ t = H ^ | Ψ 〉 , (15)

of the form

| Ψ ( t ) 〉 = ∑ η = 1 2 n e − i E η t / ℏ D η ( t ) | η 〉 , (16)

we get the equation for the coefficients D ξ ( t ) as

i D ˙ ξ ( t ) = − i Ω cos θ 2 ∑ η = 1 2 n ∑ k = 1 n e i ω ξ η t D η ( t ) 〈 ξ | S k + − S k − | η 〉 + Ω ( − 1 ) ξ k 2 D ξ ( t ) sin θ , (17)

where the frequency ω ξ η has been defined as

ω ξ η = E ξ − E η ℏ . (18)

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

τ = ω 0 t , ω 0 = 2 π k H z . (19)

Then the constant Ω / ω 0 , J / ω 0 and ω i / ω 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 π k H z .

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

1 -qubit : ω 1 = 100 , Ω = 0.1

2 -qubits : ω 1 = 100 , ω 2 = 200 , Ω = 0.1 , J = 2

3 -qubits : ω 1 = 100 , ω 2 = 200 , ω 3 = 300 , Ω = 0.1 , J = 2 , J ′ = 1 / 4 (20)

and the initial conditions are

1 -qubit : D 1 ( 0 ) = 1 , D 2 ( 0 ) = 0

2 -qubits : D 1 ( 0 ) = 1 16 , D 2 ( 0 ) = 2 16 , D 3 ( 0 ) = 5 16 , D 4 ( 0 ) = 8 16

3 -qubits : D 1 ( 0 ) = D 2 ( 0 ) = 1 8 , D 3 ( 0 ) = D 4 ( 0 ) = 0 , D 5 ( 0 ) = 3 8 , D 7 ( 0 ) = 2 8 , D 8 ( 0 ) = 1 8 (21)

^{-5}. In addition, we must observe that a superposition of states

(Hadamar quantum gate) is gotten for a time evolution of τ = π / Ω (π-pulse), that is the states of the form

1 2 ( | 0 〉 + | 1 〉 ) ( 1 -qubit ) (22)

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 〉 + | 11 〉 and | 110 〉 + | 111 〉 are obtained too.

Quantum teleportation [

| ψ 〉 = c 1 | 0 〉 + c 2 | 1 〉 , c 1 = 3 8 , c 2 = 5 8 , (23)

the initial state of our 3-qubits system is then

| Ψ 0 〉 = | ψ 〉 ⊗ | 0 〉 ⊗ | 0 〉 = c 1 | 000 〉 + c 2 | 100 〉 . (24)

The “ideal” (mathematical) teleportation algorithm is represented by quantum circuit shown on

T = ( H 3 ⋅ C N 32 ) ⋅ ( C N 21 ⋅ H 2 ) , (25)

where C N i j is the CNOT operation between the qubit “i” and the cubit “j”

C N i j | ξ n , ⋯ , ξ i , ⋯ , ξ j , ⋯ , ξ 1 〉 = | ξ n , ⋯ , ξ i , ⋯ , ξ j ⊕ ξ i , ⋯ , ξ 1 〉 , (26)

and H is the Hadamar operator (1-qubit). The first term, ( C N 21 ⋅ H 2 ) , makes the entanglement between Alice(second qubit) and Bob (first qubit), and the term ( H 3 ⋅ C N 32 ) 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

H ^ 0 = − ∑ i = 1 n ω i S i z + 2 J ℏ ∑ i = 1 n − 1 S i z S i + 1 z + 2 j ′ ℏ ∑ i = 1 n − 2 S i z S i + 2 z (27)

where N 3 (NOT gate), Z 3 and N 3 Z 3 act in the following way

N 3 | ψ 〉 = c 1 | 1 〉 + c 2 | 0 〉 , Z 3 | ψ 〉 = c 1 | 0 〉 − c 1 | 1 〉 , N 3 Z 3 | ψ 〉 = c 2 | 0 〉 − c 1 | 1 〉 . (28)

Thus, the state | ψ 〉 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 | ψ 〉 ( 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

R i j ( τ ) , (29)

the pulse of length τ done in the system with a resonant frequency ω = ( E i − E j ) / ℏ . Our non-ideal teleportation operator is given by

T ˜ = [ R 26 ( π ) R 37 ( π ) R 78 ( π ) R 15 ( π ) R 78 ( 2 π ) R 57 ( 2 π ) ] [ R 78 ( 2 π ) R 34 ( 2 π ) R 57 ( π ) R 13 ( π ) ] , (30)

which represents π and 2 π 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

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

G = ∑ i = 1 8 | | D i ( i d e a l ) | 2 − | D i ( r e a l ) | 2 | , (31)

which goes from zero (excellent reproducibility) to one (non reproducibility at all). In our case, the value obtained is G = 7.4 × 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-Shannon entropy and its average value,

S ( τ ) = − ∑ i = 1 8 | D i ( τ ) | 2 ln | D i ( τ ) | 2 , 〈 S ( T ) 〉 = 1 T ∫ 0 T S ( τ ) d τ . (32)

Of course, the maximum value this entropy can have is went all the probabilities are the same (1/8), S max = ln 8 ≈ 2.08 .

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 components is along the static field). The universal CNOT and CCNOT quantum gates and the one qubit rotation were well established on the simulations with multiples of a π-pulse. The teleportation algorithm was also implemented using these quantum gates in a 3-qubits configuration and using 10 pulses, and its associated Boltzmann-Shannon’s entropy was determined.

The authors declare no conflicts of interest regarding the publication of this paper.

López, G.V. and Lizarraga, J. (2020) C^{12} - C^{13} Diamond Quantum Computer with Longitudinal rf-Magnetic Field. Journal of Modern Physics, 11, 1123-1131. https://doi.org/10.4236/jmp.2020.118070