^{1}

^{*}

^{1}

^{*}

We make a numerical study of decoherence on the teleportation algorithm implemented in a linear chain of three nuclear spins system. We study different types of environments, and we determine the associated decoherence time as a function of the dissipative parameter. We found that the dissipation parameter to get a well defined quantum gates (without significant decoherence) must be within the range of γ≤4×10^{-4} for not thermalized case, which was determined by using the purity parameter calculated at the end of the algorithm. For the thermalized case the decoherence is stablished for very small dissipation parameter, making almost not possible to implement this algorithm for not zero temperature.

It is well known now that in the real world the interaction of the system with the environment is almost unavoidable. The study of this type of systems implies a many bodies problem which is unsolvable within any picture of the quantum mechanics. At this moment, there are two approaches to attack this type of problems. The first one consists on to look for the phenomenological classical dissipative system and to get its associated Hamiltonian, then to proceed to do the usual quantization of the system [

In the first part of our work, we describe the model and the Hamiltonian of our quantum system, and we must point out that, although this Hamiltonian will be time explicitly dependent, if we consider weak interaction between our system and the environment (the characteristic times of the quantum system are much longer than those of the environment) as a first approximation, the above mentioned Markovian-Lindblad master type equation can be still used for our study [11,17,18]. One needs to mention that even this linear chain of three nuclear spins model for solid state quantum computer has not been built yet, it has been very useful for theoretical studies about implementation of quantum gates and quantum algorithms [19-21] which can be extrapolated to other solid state quantum computers [22,23]. After doing this, we establish the five cases to be considered with the quantum-environment system: independent environment interaction (A), pure dephasing interaction (B), correlated dissipation interaction (C), dephasing correlated interaction (D), and thermalization case (E). The analyticcal dynamical systems of the reduced density matrix elements are obtained for these cases, and the results of the numerical simulations are presented. Finally, we study the behavior of the purity parameter for the teleportation algorithm.

The Hamiltonian that describes the ideal insulated system of a linear chain of N paramagnetic atoms with nuclear spin one half inside the magnetic field

where and are the amplitude, the angular frequency and the phase of the RF-field, and represents the amplitude of the z-component of the magnetic field, is given by [

where represent the magnetic moment of the kthnucleus, which it is given in terms of the nuclear spin as, with being the proton gyromagnetic ratio and being the jth-component of the spin operator, represents the magnetic field, Equation (1) valuated at the location of the kth-nuclear spin. The parameters and represent the coupling constant at first and second neighbor interaction. The angle between the linear chain and the z-component of the magnetic field is chosen as to eliminate the dipole-dipole interaction between the spins.

This Hamiltonian can be written as

where and are defined as

and

where we have used the relation, with the operator written in terms of Pauli matrixes as. Here we have that: is the Larmor frequency of the kth-spin, is the Rabi frequency, and represents the ascend operator or the descend operator. The Hamiltonian is diagonal in the basis with (zero for the exited state and zero for the exited state ). In this work, we consider that the action of the spin operators on its respective qubit is given by

, and. The eigenvalues of in this basis are given by

The elements of this basis forms a register of - qubits with a total number of registers which is the dimensionality of our Hilbert space. The allowed transition of one state to another one is gotten by choosing the angular frequency of the RF-field, , as the associated angular frequency due to the energy difference of these two levels, and by choosing the normalized evolution time with the proper time duration (so called RFfield pulse). The set of selected pulses defines the quantum gates and quantum algorithm. The energy difference between two eigenstates of is

with depending on the state of the first and second neighbor of the kth-quit. Any gate is realized applying RF-pulses of rectangular shape and choosing the radio frequency in resonance with the desired transition,. The unitary evolution of the system is denoted by, where denotes the type of pulse of duration [

Denoting by our three-qubit register, where

represents the state to be teleported by Alice, , to Bob, , the initial state of our system is given by

where one has used also the decimal notation,

, meaning that the initial reduced density matrix has the following expression

Then, a resonant -pulse between the states and, , with frequency, and a resonant -pulse between the states and, , with frequency, are applied to get the state

To put Alice and Bob in an entangled state, a CNOT gate between these elements is done, that is, it is applied a resonant π-pulse between the states and with frequency, , and a resonant π- pulse between the states and with frequency, , getting the state

Now, to get the teleportation of the state to Bob, we make an Hadamar gate to the state, that is, we use a resonant -pulse between the states and

with frequency, , a resonant -pulse between the states and with frequency, , a resonant -pulse between the states and with frequency and phase, , and a resonant -pulse between the sates and with frequency and phase, , to get finally the state (binary notation)

where and N represent the identity operator,

, z-Pauli matrix, , and NOT quantum gate, (i = 0,1), acting on the teleported state at Bob location, meaning that all the final density matrix elements would be (without interaction with environment) different from zero. For our numerical studies, we have selected the following coefficients for the state.

In summary, one has the following eight pulses to get the teleportation algorithm,

Now, to deal with the non ideal situation where the effect of the environment is taken into account, we make use of the Lindblad type master equation for the evolution of the reduced density matrix

where the first part on the right side denotes the usual von Neuman unitary evolution of the reduced density matrix, and the second term represents the Lindblad part (non unitary) evolution. This second term has different expression for different consideration of the system-environment interaction. For the qubits interacting independently with the environment (case (A)), this term has the following form [

where is the dissipative parameter associated to the jth-qubit.

For the pure dephasing interaction case, where the qubits independently dephase to their respective bath with a dephasing rate, the Lindblad term is given by

For the independent-qubit-correlated case (qubits interact with the environment collectively), the Lindblad operator is written as

where one has that is the decay rate of case (A). In this case, the decay of the state of a qubit has an effect on the other qubits.

For the qubit-correlated and dephasing case, with as the decay rate of the correlated dephasing, the Lindblad operator is given by

In this case, the decay of one qubit affects too the other qubits.

Finally, we define another environmental description for the Lindblad term related to a thermalization process of the system. This system-environment interaction induces an energy absorption process leading the system into a mixed thermalized state. The environment is now at a certain finite temperature, and it can be thought as field radiation modes contained in a cavity where the central system lies. The Lindblad term has the following form [11,24]

where the damping factors are now functions of the temperature and the characteristic frequencies related to the eigenenergies of the closed system. They have the form

and the function

represents the Planck’s distribution function, and are phenomenological damping factors which depend on the cavity, the eigenfrequencies of the system and the strength of the coupling between the system and the environment. We can manipulate the dissipation parameters by considering a low or high strength of the coupling between the system and the environment, and also some other phenomenological parameters like the volume of the cavity. In this way, we have some freedom to modulate the damping factors. We want to point out that if we go to temperature equals to zero, (the thermal vacuum), the case (A) is recovered since for,

and,.

The dynamical system for each case for the reduced density matrix elements is deduced from Equation (15) as

where and represent the elements of the basis of the Hilbert space,

and.

In our case, one has that, the dimensionality of our Hilbert space is eight, and the explicit equations for the dynamical system of each case be see in the appendix of reference [24,25]. We have considered that it is not necessary to repeat those equations in this manuscript.

Our registers are made up of three qubits with, (also denoted as, do not confuse with the type of environment) or written them with decimal notation, , and so on. The parameters used for our simulation are taken from [

The selected teleported state is defined by the coefficients (13) of the state. Assuming that the environment acts homogeneously on the qubits, the damping parameter can be the same for each qubit, and the damping parameter for correlated cases at second neighbors can be one order of magnitude weaker that at first neighbors. Thus, the dissipative coefficients appearing for the cases (A), (B), (C), and (D) are taken of the following way

where is the free common parameter which takes into account the interaction with the environment. The reduced density matrix is then made up of complex elements, and if the initial reduced density matrix is given by (9).

behavior as (B) case has. Therefore they are not presented on the figures. The thermalized case is presented for T = 2˚K and the destruction of the algorithm is evident for very low dissipation.

high dissipation ranges, and for the case B (independent).

In

We have made a numerical study of decoherence on the teleportation algorithm implemented in a linear chain of three nuclear spins system. We have studied different types of environments, and we have determined the associated decoherence time as a function of the dissipative parameter. We have used the purity value at the end of the algorithm as a quality factor to determine the behavior of the teleportation algorithm. With this parameter and with the selection of the other parameters as (24), we have found that the dissipation parameter to get a well defined quantum gates (without significant decoherence) must be within the range of for the non thermalization case. For high dissipation parameter we observed the expected recovery of the purity since must be conserved. With the selected dissipative coefficients, the cases (C) and (D), corresponding to correlation between spins, have little contribution to the cases (A) and (B) (independent and dephasing cases), and the most danger situation corresponds to the A-independent case when thermalization is not taken into account. However, when thermalization is taken into account, we have shown that even for very low dissipation parameter and very low temperature, the destruction of the algorithm occurs.