^{1}

^{1}

^{2}

^{2}

Through the study of the factorization conditions of a wave function made up of two, three and four qubits, we propose an analytical expression which can characterize entangled states in terms of the coefficients of the wave function and density matrix elements.

In quantum mechanics, multiple dimensional or multiple particles systems are characterized by the tensor product of the Hilbert subspaces [

When a quantum systems is made up of several quantum subsystems, where the ith-subsystem is characterized by a Hamiltonian H i and a Hilbert space ε i , corresponding a two-states and having a basis { | ξ i 〉 } ξ i = 0 , 1 , the Hilbert space is written as the tensorial product of each subsystem, ε = ε n ⊗ ⋯ ⊗ ε 1 (for n- subsystems), The state in each subsystems is defined as a qubit in quantum computation and information theory [

| ψ i 〉 = a i | 0 〉 + b i | 1 〉 , | a i | 2 + | b i | 2 = 1 , (1)

where { | 0 〉 , | 1 〉 } is the basis of the two-states Hilbert subspace ε i . A general state | Ψ 〉 in the Hilbert space ε can be written as

| Ψ 〉 = ∑ ξ C ξ | ξ 〉 , with ∑ ξ | C ξ | 2 = 1 , (2)

where C ′ ξ s are complex numbers, and | ξ 〉 is an element of the basis of ε ,

| ξ 〉 = | ξ n ⋯ ξ 1 〉 = | ξ n 〉 ⊗ ⋯ ⊗ | ξ 1 〉 , ξ k = 0 , 1 k = 1 , ⋯ , n . (3)

A full factorized state in this space is

| Ψ 〉 = | ψ n 〉 ⊗ ⋯ ⊗ | ψ 1 〉 , (4)

where | ψ k 〉 is given by (1).

For n = 2 , one has a general state | Ψ 〉 ∈ ε ,

| Ψ 〉 = C 1 | 00 〉 + C 2 | 01 〉 + C 3 | 10 〉 + C 4 | 11 〉 , (5)

where we have chosen to use decimal notation for the coefficients. Let us assume that this state can be written as

| Ψ 〉 = | ψ 2 〉 ⊗ | ψ 1 〉 , (6)

with | ψ k 〉 , k = 1 , 2 given by (1). Then, substituting (1) in (6) and equaling coefficients with (5), it follows that

C 1 = a 1 a 2 , C 2 = a 1 b 2 , C 3 = b 1 a 2 , C 4 = b 1 b 2 . (7)

From these expressions one obtains a single condition for factorization,

C 1 C 4 = C 2 C 3 . (8)

Thus, if this condition is not satisfied, the state (5) represents an entangled state. So, one can use the following known expression [

C Ψ ( 2 ) = 2 | C 1 C 4 − C 2 C 3 | . (9)

For n = 3 , a general state in the Hilbert space ε is

| Ψ 〉 = C 1 | 00 〉 + C 2 | 001 〉 + C 3 | 010 〉 + C 4 | 011 〉 + C 5 | 100 〉 + C 6 | 101 〉 + C 7 | 110 〉 + C 8 | 111 〉 , (10)

where, we have chosen again decimal notation for the coefficients. Assuming that this wave function can be written as | Ψ 〉 = | ψ 3 〉 ⊗ | ψ 2 〉 ⊗ | ψ 1 〉 with | ψ k 〉 given by (1), and after some identifications (as before) and rearrangements, one gets

C 1 C 2 = C 3 C 4 = C 5 C 6 = C 7 C 8 = a 3 b 3 , b 3 = 0 (11a)

C 1 C 3 = C 2 C 4 = C 5 C 7 = C 6 C 8 = a 2 b 2 , b 2 = 0 (11b)

C 1 C 5 = C 2 C 6 = C 3 C 7 = C 4 C 8 = a 1 b 1 , b 1 = 0 (11c)

These expression reflex a parallelism between the complex vectors ( C 1 , C 3 , C 5 , C 7 ) and ( C 2 , C 4 , C 6 , C 8 ) , the vectors ( C 1 , C 2 , C 5 , C 6 ) and ( C 3 , C 4 , C 7 , C 8 ) , and the vectors ( C 1 , C 2 , C 3 , C 4 ) and ( C 5 , C 6 , C 7 , C 8 ) . In addition, they bring about the following eight independent relations

C 1 C 4 − C 2 C 3 = 0 , C 1 C 6 − C 2 C 5 = 0 , C 1 C 8 − C 2 C 7 = 0 , C 3 C 6 − C 4 C 5 = 0 (12)

C 3 C 8 − C 4 C 7 = 0 , C 5 C 8 − C 6 C 7 = 0 , C 1 C 7 − C 3 C 5 = 0 , C 2 C 8 − C 4 C 6 = 0. (13)

If one of these expressions is not satisfied, the wave function (10) represents an entangled state. Therefore, one can propose the following expression to characterize an entangled state

C Ψ ( 3 ) = 2 | C 1 C 4 − C 2 C 3 | + 2 | C 1 C 6 − C 2 C 5 | + 2 | C 1 C 8 − C 2 C 7 | + 2 | C 3 C 6 − C 4 C 5 | + 2 | C 3 C 8 − C 4 C 7 | + 2 | C 5 C 8 − C 6 C 7 | + 2 | C 1 C 7 − C 3 C 5 | + 2 | C 2 C 8 − C 4 C 6 | . (14)

For n = 4 , a general state in the Hilbert space ε is of the form

| Ψ 〉 = C 1 | 0000 〉 + C 2 | 0001 〉 + C 3 | 010 〉 + C 4 | 0011 〉 + C 5 | 0100 〉 + C 6 | 0101 〉 + C 7 | 0110 〉 + C 8 | 0111 〉 + C 9 | 100 〉 + C 10 | 1001 〉 + C 11 | 110 〉 + C 12 | 1011 〉 + C 13 | 1100 〉 + C 14 | 1101 〉 + C 15 | 1110 〉 + C 16 | 1111 〉 . (15)

Assuming this function can be expressed as | Ψ 〉 = | ψ 4 〉 ⊗ | ψ 3 〉 ⊗ | ψ 2 〉 ⊗ | ψ 1 〉 with | ψ k 〉 , k = 1 , 2 , 3 , 4 given by (1), and after some identifications and rearrangements, one gets

C 1 C 2 = C 3 C 4 = C 5 C 6 = C 7 C 8 = C 9 C 10 = C 11 C 12 = C 13 C 14 = C 15 C 16

C 1 C 3 = C 2 C 4 = C 5 C 7 = C 6 C 8 = C 9 C 11 = C 10 C 12 = C 13 C 15 = C 14 C 16

C 1 C 5 = C 2 C 6 = C 3 C 7 = C 4 C 8 = C 9 C 13 = C 10 C 14 = C 11 C 15 = C 12 C 16

C 1 C 9 = C 2 C 10 = C 3 C 11 = C 4 C 12 = C 5 C 13 = C 6 C 14 = C 7 C 15 = C 8 C 16 ,

expressing similar parallelism we mentioned before. Each row gives us 28 relations, having a total of 112 possible relations, and from these relations, one can get the following 36 independent conditions

C 1 C 4 − C 2 C 3 = 0 C 4 C 13 − C 7 C 10 = 0 C 1 C 6 − C 2 C 5 = 0 C 4 C 14 − C 6 C 12 = 0 (16a)

C 1 C 8 − C 3 C 6 = 0 C 4 C 15 − C 3 C 16 = 0 C 1 C 10 − C 2 C 9 = 0 C 4 C 16 − C 8 C 12 = 0 (16b)

C 1 C 11 − C 3 C 9 = 0 C 5 C 8 − C 6 C 7 = 0 C 1 C 12 − C 2 C 11 = 0 C 5 C 14 − C 6 C 13 = 0 (16c)

C 1 C 14 − C 9 C 6 = 0 C 5 C 15 − C 7 C 13 = 0 C 1 C 15 − C 5 C 11 = 0 C 5 C 16 − C 7 C 14 = 0 (16d)

C 2 C 8 − C 4 C 6 = 0 C 6 C 11 − C 5 C 12 = 0 C 2 C 12 − C 4 C 10 = 0 C 6 C 15 − C 16 C 8 = 0 (16e)

C 2 C 13 − C 5 C 10 = 0 C 6 C 16 − C 8 C 14 = 0 C 2 C 14 − C 6 C 10 = 0 C 7 C 16 − C 8 C 15 = 0 (16f)

C 2 C 16 − C 10 C 8 = 0 C 7 C 12 − C 8 C 11 = 0 C 3 C 8 − C 4 C 7 = 0 C 9 C 12 − C 10 C 11 = 0 (16g)

C 3 C 15 − C 7 C 11 = 0 C 9 C 14 − C 10 C 11 = 0 C 3 C 13 − C 11 C 5 = 0 C 9 C 15 − C 11 C 13 = 0 (16h)

C 10 C 16 − C 12 C 14 = 0 C 11 C 16 − C 12 C 15 = 0 C 10 C 15 − C 11 C 14 = 0 C 13 C 16 − C 14 C 15 = 0. (16i)

Again, if one of these expression fail to happen, (15) will represents an entangled state. Thus, one can propose the following expression to characterize an entangled state made up of 4-qubits basis

C Ψ ( 4 ) = 2 | C 1 C 4 − C 2 C 3 | + 2 | C 4 C 13 − C 7 C 10 | + 2 | C 1 C 6 − C 2 C 5 | + 2 | C 4 C 14 − C 6 C 12 | + 2 | C 1 C 8 − C 3 C 6 | + 2 | C 4 C 15 − C 3 C 16 | + 2 | C 1 C 10 − C 2 C 9 | + 2 | C 4 C 16 − C 8 C 12 | + 2 | C 1 C 11 − C 3 C 9 | + 2 | C 5 C 8 − C 6 C 7 | + 2 | C 1 C 12 − C 2 C 11 | + 2 | C 5 C 14 − C 6 C 13 | + 2 | C 1 C 14 − C 9 C 6 | + 2 | C 5 C 15 − C 7 C 13 | + 2 | C 1 C 15 − C 5 C 11 | + 2 | C 5 C 16 − C 7 C 14 | + 2 | C 2 C 8 − C 4 C 6 | + 2 | C 6 C 11 − C 5 C 11 | + 2 | C 2 C 12 − C 4 C 10 | + 2 | C 6 C 15 − C 16 C 8 | + 2 | C 2 C 13 − C 5 C 10 | + 2 | C 6 C 16 − C 8 C 14 | + 2 | C 2 C 14 − C 6 C 10 | + 2 | C 7 C 16 − C 8 C 15 | + 2 | C 2 C 16 − C 10 C 8 | + 2 | C 7 C 12 − C 8 C 11 | + 2 | C 3 C 8 − C 4 C 7 | + 2 | C 9 C 12 − C 10 C 11 | + 2 | C 3 C 15 − C 7 C 11 | + 2 | C 9 C 14 − C 10 C 11 | + 2 | C 3 C 13 − C 11 C 5 | + 2 | C 9 C 15 − C 11 C 13 | + 2 | C 10 C 16 − C 12 C 14 | + 2 | C 11 C 16 − C 12 C 15 | + 2 | C 10 C 15 − C 11 C 14 | + 2 | C 13 C 16 − C 14 C 15 | . (17)

As we can see from these examples, the number of conditions needed to characterize a factorized state (or entangled state) grows exponentially with the number of qubits. So, characterization of an entangled state for n-qubits in general becomes a very hard work. Now, in terms of the density matrix elements, one could have the characterization of entangled states made up of 2, 3 and 4 qubits as

C ρ ( 2 ) = 2 ρ 11 ρ 44 + ρ 22 ρ 33 − 2 Re ( ρ 12 ρ 43 ) . (18)

C ρ ( 3 ) = 2 ρ 11 ρ 44 + ρ 22 ρ 33 − 2 Re ( ρ 12 ρ 43 ) + 2 ρ 11 ρ 66 + ρ 22 ρ 55 − 2 Re ( ρ 12 ρ 65 ) + 2 ρ 11 ρ 88 + ρ 22 ρ 77 − 2 Re ( ρ 12 ρ 87 ) + 2 ρ 33 ρ 66 + ρ 44 ρ 55 − 2 Re ( ρ 34 ρ 65 ) + 2 ρ 33 ρ 88 + ρ 44 ρ 77 − 2 Re ( ρ 34 ρ 87 ) + 2 ρ 55 ρ 88 + ρ 66 ρ 77 − 2 Re ( ρ 56 ρ 87 ) + 2 ρ 11 ρ 77 + ρ 33 ρ 55 − 2 Re ( ρ 13 ρ 75 ) + 2 ρ 22 ρ 88 + ρ 44 ρ 66 − 2 Re ( ρ 24 ρ 86 ) . (19)

For other considerations of entanglement multiqubit entanglement see [

| W 〉 = C 2 | 001 〉 + C 3 | 010 〉 + C 5 | 100 〉 , | C 2 | 2 + | C 3 | 2 + | C 5 | 2 = 1 , (21)

the possible values of C ( 3 ) ( Ψ ) . As we can see, there are four possible maxima corresponding to the values C 2 = C 3 = C 5 = ± 1 / 3 and four maxima corres-

ponding to the values C 5 = 0 , C 2 = C 3 = ± 1 / 2 , related to semi-factorized state | 0 〉 ⊗ ( C 2 | 01 〉 + C 3 | 10 〉 ) .

| G H Z 〉 = C 1 | 000 〉 + C 8 | 111 〉 , | C 1 | 2 + | C 8 | 2 = 1. (22)

The maximum value of C ( 3 ) ( Ψ ) is gotten for C 1 = C 8 = 1 / 2 , as one would expect.

For a Hilbert space ε generated by n-qubits, C ( n ) defines a continuous function C ( n ) : ε → ℜ + with ℜ + = [ 0 , + ∞ ) , c ≥ 1 . Since the coefficients of the wave function defines a compact set on the real space ℜ 2 n due to the relation ∑ i | C i | 2 = 1 , the image of this compact set is a compact set in ℜ + [

Following Lloyd’s idea [

H s = − ∑ k = 1 N ω k S k z − 2 J ℏ ∑ k = 1 N − 1 S k z S k + 1 z − 2 J ′ ℏ ∑ k = 1 N − 2 S k z S k + 2 z , (23)

where N is the number of nuclear spins in the chain (or qubits), J and J ′ are the coupling constant of the nucleus at first and second neighbor. Using the

basis of the register of N-qubits, { | ξ N , ⋯ , ξ 1 〉 } with ξ k = 0 , 1 , one has that S k z | ξ k 〉 = ( − 1 ) ξ k ℏ | ξ k 〉 / 2 . Therefore, the Hamiltonian is diagonal on this basis,

and its eigenvalues are

E ξ = − ℏ 2 ∑ k = 1 N ( − 1 ) ξ k ω k − J ℏ 2 ∑ k = 1 N − 1 ( − 1 ) ξ k + ξ k + 1 − J ′ ℏ 2 ∑ k = 1 N − 2 ( − 1 ) ξ k + ξ k + 2 . (24)

Consider now that the environment is characterized by a Hamiltonian H e and its interacting with the quantum system with Hamiltonian H s . Thus, the total Hamiltonian would be H = H s + H e + H s e , where H s e is the part of the Hamiltonian which takes into account the interaction system-environment, and the equation one would need to solve, in terms of the density matrix, is [

i ℏ ∂ ρ t ∂ t = [ H , ρ t ] , (25)

where ρ t = ρ t ( s , e ) is the density matrix which depends on the system and environment coordinates. The evolution of the system is unitary, but it is not possible to solve this equation due to a lot of degree of freedom. Therefore, under some approximations and tracing over the environment coordinates [

i ℏ ∂ ρ ∂ t = [ H s , ρ ] + ∑ i = 1 I { V i ρ V i † − 1 2 V i † V i ρ − 1 2 ρ V i V i † } (26)

where V i are called Kraus’ operators. This equation is not unitary and Markovian (without memory of the dynamical process). This equation can be written in the interaction picture, through the transformation ρ ˜ = U ρ U † with U = e i H s t / ℏ , as

i ℏ ∂ ρ ˜ ∂ t = L ˜ ( ρ ˜ ) , (27)

where L ˜ ( ρ ˜ ) is the Lindblad operator

L ˜ ( ρ ˜ ) = ∑ i = 1 I { V ˜ i ρ ˜ V ˜ i † − 1 2 V ˜ i † V ˜ i ρ ˜ − 1 2 ρ ˜ V ˜ i V ˜ i † } (28)

with V ˜ = U V U † . The explicit form of Lindblad operator is determined by the type of environment to consider [

L ˜ ( ρ ˜ ) = 1 2 i ℏ ∑ k N γ k ( 2 S ˜ k − ρ ˜ S ˜ k + − S k + S k − ρ ˜ − ρ ˜ S ˜ k + S ˜ k − ) (29)

where S ˜ k + and S ˜ k − are the ascend and descend operators such that

S ˜ k ± = U S k ± U † = S k ± e ± i Ω ^ k t

where Ω ^ k is defined as Ω ^ k = w k + J ℏ ( S k + 1 z + S k − 1 z ) + J ′ ℏ ( S k + 2 z + S k − 2 z ) . . The solutions of the equations are

ρ 11 ( t ) = ρ 11 ( 0 ) + ρ 22 ( 0 ) + ρ 33 ( 0 ) + ρ 44 ( 0 ) + ρ 55 ( 0 ) + ρ 66 ( 0 ) + ρ 77 ( 0 ) + ρ 88 ( 0 ) − ( ρ 55 ( 0 ) + ρ 66 ( 0 ) + ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − γ 1 t − ( ρ 33 ( 0 ) + ρ 44 ( 0 ) + ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − γ 2 t − ( ρ 22 ( 0 ) + ρ 44 ( 0 ) + ρ 66 ( 0 ) + ρ 88 ( 0 ) ) e − γ 3 t + ( ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 2 ) t + ( ρ 66 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 3 ) t + ( ρ 44 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 2 + γ 3 ) t − ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 14 ( t ) = [ ρ 14 ( 0 ) + γ 1 ρ 58 ( 0 ) e − i ϕ 14 γ 1 2 + ( j + j ′ ) 2 ( 1 − e ( i ( j + j ′ ) − γ 1 ) t ) ] e − 1 2 ( γ 2 + γ 3 ) t ;

ρ 16 ( t ) = [ ρ 16 ( 0 ) + γ 2 ρ 38 ( 0 ) e − i ϕ 16 γ 2 2 + 4 j 2 ( 1 − e ( 2 i j − γ 2 ) t ) ] e − 1 2 ( γ 1 + γ 3 ) t ;

ρ 17 ( t ) = [ ρ 17 ( 0 ) + γ 3 ρ 28 ( 0 ) e − i ϕ 17 γ 3 2 + ( j + j ′ ) 2 ( 1 − e [ i ( j + j ′ ) − γ 3 ] t ) ] e − 1 2 ( γ 1 + γ 2 ) t ;

ρ 18 ( t ) = ρ 18 ( 0 ) e − 1 2 ( γ 1 + γ 2 + γ 3 ) t ;

ρ 22 ( t ) = ( ρ 22 ( 0 ) + ρ 44 ( 0 ) + ρ 66 ( 0 ) + ρ 88 ( 0 ) ) e − γ 3 t − ( ρ 66 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 3 ) t − ( ρ 44 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 2 + γ 3 ) t + ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 23 ( t ) = [ ρ 23 ( 0 ) + γ 1 ρ 67 ( 0 ) e − i ϕ 23 γ 1 2 + ( j − j ′ ) 2 ( 1 − e ( i ( j − j ′ ) − γ 1 ) t ) ] e − 1 2 ( γ 2 + γ 3 ) t ;

ρ 25 ( t ) = [ ρ 25 ( 0 ) + ρ 47 ( 0 ) ( 1 − e − γ 2 t ) ] e − 1 2 ( γ 1 + γ 3 ) t ;

ρ 27 ( t ) = ρ 27 ( 0 ) e − 1 2 ( γ 1 + γ 2 + γ 3 ) t ;

ρ 28 ( t ) = ρ 28 ( 0 ) e − 1 2 ( γ 1 + γ 2 + 2 γ 3 ) t ;

ρ 33 ( t ) = ( ρ 33 ( 0 ) + ρ 44 ( 0 ) + ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − γ 2 t − ( ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 2 ) t − ( ρ 44 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 2 + γ 3 ) t + ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 35 ( t ) = [ ρ 35 ( 0 ) + γ 3 ρ 46 ( 0 ) e − i ϕ 35 γ 3 2 + ( j − j ′ ) 2 ( 1 − e − [ i ( j − j ′ ) + γ 3 ] t ) ] e − 1 2 ( γ 1 + γ 2 ) t ;

ρ 36 ( t ) = ρ 36 ( 0 ) e − 1 2 ( γ 1 + γ 2 + γ 3 ) t ;

ρ 38 ( t ) = ρ 38 ( 0 ) e − 1 2 ( γ 1 + 2 γ 2 + γ 3 ) t ;

ρ 44 ( t ) = ( ρ 44 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 2 + γ 3 ) t − ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 45 ( t ) = ρ 45 ( 0 ) e − 1 2 ( γ 1 + γ 2 + γ 3 ) t ;

ρ 46 ( t ) = ρ 46 ( 0 ) e − 1 2 ( γ 1 + γ 2 + 2 γ 3 ) t ;

ρ 47 ( t ) = ρ 47 ( 0 ) e − 1 2 ( γ 1 + 2 γ 2 + γ 3 ) t ;

ρ 48 ( t ) = ρ 48 ( 0 ) e − 1 2 ( γ 1 + 2 γ 2 + 2 γ 3 ) t ;

ρ 55 ( t ) = ( ρ 55 ( 0 ) + ρ 66 ( 0 ) + ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − γ 1 t − ( ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 2 ) t − ( ρ 66 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 3 ) t + ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 58 ( t ) = ρ 58 ( 0 ) e − 1 2 ( 2 γ 1 + γ 2 + γ 3 ) t ;

ρ 66 ( t ) = ( ρ 66 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 3 ) t − ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 67 ( t ) = ρ 67 ( 0 ) e − 1 2 ( 2 γ 1 + γ 2 + γ 3 ) t ;

ρ 68 ( t ) = ρ 68 ( 0 ) e − 1 2 ( 2 γ 1 + γ 2 + 2 γ 3 ) t ;

ρ 77 ( t ) = ( ρ 77 ( 0 ) + ρ 88 ( 0 ) ) e − ( γ 1 + γ 2 ) t − ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t ;

ρ 78 ( t ) = ρ 78 ( 0 ) e − 1 2 ( 2 γ 1 + 2 γ 2 + γ 3 ) t ;

ρ 88 ( t ) = ρ 88 ( 0 ) e − ( γ 1 + γ 2 + γ 3 ) t .

where ϕ i j are given by

ϕ 14 = tan − 1 ( j + j ′ γ 1 ) ; ϕ 16 = tan − 1 ( 2 j γ 2 ) ; ϕ 17 = tan − 1 ( j + j ′ γ 3 ) ; ϕ 23 = tan − 1 ( j − j ′ γ 1 ) ; ϕ 35 = tan − 1 ( j − j ′ γ 3 ) .

In our case, we have three qubits space { | ξ 3 ξ 2 ξ 1 〉 } ξ i = 0 , 1 , and our parameter in units 2 πMHz are

ω 1 = 400 ; ω 2 = 200 ; ω 3 = 100 J = 10 ; J ′ = 0.4

γ 1 = 0.05 ; γ 2 = 0.05 ; γ 3 = 0.05.

the time is normalized by the same factor of 2 πMHz , and we include in this study the entangle state

| Ψ 1 〉 = 1 2 ( | 000 〉 + | 111 〉 + | 001 〉 + | 110 〉 ) . (30)

evolves in a mixed state and finishes in the pure ground state ( | 000 〉 ) , after sharing energy with the environment. The states | G H Z 〉 and | Ψ 1 〉 behave more robust than the state | W 〉 , C ρ ( 3 ) grows since other entangled stated con- tribute to this function. In contrast, starting with the entangled state | W 〉 , there are not other entangled states which make contribution to the function C ρ ( 3 ) in the dynamics, and one sees an exponential decay.

We have studied the full factorization of a state made up of up to 4-qubits basic states. We have seen that there is an indication that the number of conditions to characterize a factorized state grows exponentially with the number of qubits. For two, three and four basic qubits, we showed the conditions in order to have a factorized state, and if any of one of these conditions fails, one gets instead an entangled state. Therefore, an entangled state is also characterized by the complement of each of this conditions, and the resulting expression has been denoted by C ( n ) ( n = 2 , 3 , 4 ) . This non-negative function expressed in terms of the coefficients of the wave function or in terms of the density matrix elements represents a measurement of the entanglement of any wave function made up of basic n-qubits ( n = 2 , 3 , 4 ) . Using this function, we study the decay of entangled states | W 〉 , | G H Z 〉 , and | Ψ 1 〉 due to interaction with environment, and we noticed a great different behavior of the function C ρ ( 3 ) , indicating some type of robustness behavior of the states | G H Z 〉 and | Ψ 1 〉 . The main reason for this different behavior is that the entangled states | G H Z 〉 and | Ψ 1 〉 contain the ground state | 000 〉 , which is the final state in the dynamics.

López, G.V., Montes, G., Avila, M. and Rueda-Paz, J. (2017) About Factorization of Quantum States with Few Qubits. Journal of Applied Mathematics and Physics, 5, 469-480. https://doi.org/10.4236/jamp.2017.52042