**A New Mechanism of Open System Evolution and Its Entropy Using Unitary Transformations in Noncomposite Qudit Systems**

#### **Julio A. López-Saldívar 1,2,\*, Octavio Castaños 1, Margarita A. Man'ko <sup>3</sup> and Vladimir I. Man'ko 2,3,4**


Received: 17 May 2019; Accepted: 24 July 2019; Published: 27 July 2019

**Abstract:** The evolution of an open system is usually associated with the interaction of the system with an environment. A new method to study the open-type system evolution of a qubit (two-level atom) state is established. This evolution is determined by a unitary transformation applied to the qutrit (three-level atom) state, which defines the qubit subsystems. This procedure can be used to obtain different qubit quantum channels employing unitary transformations into the qutrit system. In particular, we study the phase damping and spontaneous-emission quantum channels. In addition, we mention a proposal for quasiunitary transforms of qubits, in view of the unitary transform of the total qutrit system. The experimental realization is also addressed. The probability representation of the evolution and its information-entropic characteristics are considered.

**Keywords:** entropy; open systems; unitary evolution; qubit; qutrit

#### **1. Introduction**

The open system evolution of a qudit state is known to be the result of interactions with an environment. Usually, the states of the complete system are thought to evolve by a unitary transformation in the Hilbert space *<sup>H</sup>*<sup>ˆ</sup> <sup>=</sup> *<sup>H</sup>*<sup>ˆ</sup> *<sup>q</sup>* <sup>⊗</sup> *<sup>H</sup>*<sup>ˆ</sup> *env*, then the density operator of the composite system leads us, using the partial tracing procedure, to the density operator of the subsystem *ρ*ˆ*<sup>q</sup>* (qudit), and its evolution is induced by the unitary evolution of the complete system. In this picture, the qubit state dynamics needs the structure of the Hilbert space *H*ˆ corresponding to the presence of two subsystems, qudit and environment [1]. In this work, we suggest a new mechanism to study the open system evolution, which does not demand the complete system to have a subsystem.

We show that for any system without subsystems, there exist a unitary evolution, which due to hidden correlations in the system, evolves according to the Gorini–Kossakowski–Sudarshan–Lindblad equation [2–5]. We demonstrate this picture using the example of a qutrit (complete system without subsystems), where the open-like evolution is available for their associated qubits.

In previous works [6–10], a new method to define different qubit density matrices from a qudit system was established. This procedure uses the occupation probabilities and transition probability amplitudes for different levels of a qudit system and groups them as if there exists two levels only. This is done by mapping the qudit density matrix to the closest higher even-dimensional density matrix. The partial trace operation then is enacted on the resulting matrix in order to obtain well-defined qubit density matrices.

The obtained qubits have been used to define a new geometric representation of the *d*-dimensional qudit states through *d* Bloch vectors [10] associated with the generated qubits. Furthermore, it has been possible to describe quantum phenomena as the entanglement on a two-qubit system in terms of standard probabilities [9].

The evolution of a qutrit density matrix can provide the quantum channel, which maps the initial state *ρ*ˆ*<sup>a</sup>* onto the density matrix *ρ*ˆ *<sup>a</sup>*. The proposed open-type evolution establishes a new mechanism, which will need a special state preparation and a specific unitary operation for the qutrit system, as we will show later on. The experimental possibilities by which one can realize this new mechanism are related to superconducting circuit devices [11,12].

Most quantum computing processes consider a set of pure qubit states, which are transformed by unitary operators, also called gates, that are used to implement different computing algorithms. In this article, instead, we have density matrices (which might be describing a mixed state) of larger qudit systems. The definition of a set of qubit states from a qudit system is similar to the ideas established in [13], where the emulation of a spin system was obtained from qudit states, and in [14], where the quantum logic of qubits was simplified by the use of a higher dimensional Hilbert space; and in general, with all the procedures that make use of larger Hilbert spaces. In this work, we demonstrate that subsystems of qubits defined by larger systems can be used in quantum information. A principal foundation of quantum computation is the study of quantum channels. These channels are linked to unitary transformations of the qubit density matrix. There exist several channels that can describe the interaction between a quantum system and its environment such as the bit-flip, depolarization, spontaneous emission, phase, and amplitude damping channels. For this, the study of quantum channels has been of relevance in the error correction theory of quantum computation [15,16].

Here, we present different examples of quantum channels, which act on the associated qubits to qudit states. These quantum channels have the advantage of being represented as unitary transformations acting in the qudit system, providing the possibility to study the qubits as if they were interacting with an environment.

On the other hand, the study of the interaction of three-level systems with electromagnetic fields has led to the discovery of important phenomena, such as the presence of dark states [17] together with black resonances [18] and electromagnetically-induced transparency [19–21]. This is important to our objectives as in some cases, the herein proposed qubit quantum channels can be obtained by a unitary transformation of dark states, suggesting the possibility of checking our results experimentally.

The work is organized as follows: In Section 2, a review of the qubit density matrices that are associated with a qutrit state is given. Furthermore, the association of a unitary transform of the qutrit to the nonunitary transformations of the qubits is studied. In Section 3, the definitions of the qubit phase damping and spontaneous-emission quantum channels are reviewed. Later, the unitary transformations of a qutrit system are explicitly given, which yields the phase damping and spontaneous-emission channels on the associated qubits. A way to obtain a quasi-unitary transformation on the qubits is also explored. The change of entropy associated with the nonunitary evolution of the qubits is discussed in Section 4. Finally, some concluding remarks are given.

#### **2. Nonunitary Evolution for the Qubit Decomposition of Qutrit States**

In a previous work [10], we showed the existence of six different qubit states associated with a general qutrit density matrix:

$$
\boldsymbol{\beta} = \begin{pmatrix}
\rho\_{11} & \rho\_{12} & \rho\_{13} \\
\rho\_{21} & \rho\_{22} & \rho\_{23} \\
\rho\_{31} & \rho\_{32} & \rho\_{33}
\end{pmatrix}.
$$

To define these states, different maps of *ρ*ˆ to a 4 × 4 density matrix, with one row and one column equal to zero (in such a way that ensures an eigenvalue equal to zero), were used. Then, the partial trace of the resulting 4 × 4 matrix was performed as if it was describing a two-qubit system. The obtained qubit partial density operators can be explicitly written as:

$$\begin{aligned} \label{eq:SD1} \dot{\rho}\_1 &= \begin{pmatrix} 1 - \rho\_{33} & \rho\_{13} \\ \rho\_{31} & \rho\_{33} \end{pmatrix}, \quad \dot{\rho}\_2 = \begin{pmatrix} 1 - \rho\_{22} & \rho\_{12} \\ \rho\_{21} & \rho\_{22} \end{pmatrix}, \quad \dot{\rho}\_3 = \begin{pmatrix} \rho\_{11} & \rho\_{13} \\ \rho\_{31} & 1 - \rho\_{11} \end{pmatrix}, \\\ \dot{\rho}\_4 &= \begin{pmatrix} \rho\_{22} & \rho\_{23} \\ \rho\_{32} & 1 - \rho\_{22} \end{pmatrix}, \quad \dot{\rho}\_5 = \begin{pmatrix} \rho\_{11} & \rho\_{12} \\ \rho\_{21} & 1 - \rho\_{11} \end{pmatrix}, \quad \dot{\rho}\_6 = \begin{pmatrix} 1 - \rho\_{33} & \rho\_{23} \\ \rho\_{32} & \rho\_{33} \end{pmatrix}. \end{aligned} \tag{1}$$

The qubit states can be characterized in different sets by their corresponding von Neumann entropy *Sk* = −Tr *ρ<sup>k</sup>* ln *ρk*, with *k* = 1, 2, ... , 6. These qubits correspond to the reduction of the three-level system to different two-level systems by the summation of the population probabilities of two levels into one.

When the qutrit state is transformed using a general three-dimensional unitary matrix *U*ˆ , i.e., *ρ*ˆ = *U*ˆ † *ρ*ˆ *U*ˆ , the qubits in Equation (1) are transformed in a nonunitary way. The transformed qubit density matrices can be written by the following expressions:

*ρ*ˆ <sup>1</sup> <sup>=</sup> <sup>1</sup> *D D* − *M*3,1*N*1,3 + *M*2,1*N*2,3 − *M*1,1*N*3,3 *M*3,3*N*1,3 − *M*2,3*N*2,3 + *M*1,3*N*3,3 *<sup>M</sup>*3,1*N*1,1 <sup>−</sup> *<sup>M</sup>*2,1*N*2,1 <sup>+</sup> *<sup>M</sup>*1,1*N*3,1 *<sup>M</sup>*3,1*N*1,3 <sup>−</sup> *<sup>M</sup>*2,1*N*2,3 <sup>+</sup> *<sup>M</sup>*1,1*N*3,3 , *ρ*ˆ <sup>2</sup> <sup>=</sup> <sup>1</sup> *D D* + *M*3,2*N*1,2 − *M*2,2*N*2,2 + *M*1,2*N*3,2 *M*3,3*N*1,2 − *M*2,3*N*2,2 + *M*1,3*N*3,2 <sup>−</sup>*M*3,2*N*1,1 <sup>+</sup> *<sup>M</sup>*2,2*N*2,1 <sup>−</sup> *<sup>M</sup>*1,2*N*3,1 <sup>−</sup>*M*3,2*N*1,2 <sup>+</sup> *<sup>M</sup>*2,2*N*2,2 <sup>−</sup> *<sup>M</sup>*1,2*N*3,2 , *ρ*ˆ <sup>3</sup> <sup>=</sup> <sup>1</sup> *D M*3,3*N*1,1 − *M*2,3*N*2,1 + *M*1,3*N*3,1 *M*3,3*N*1,3 − *M*2,3*N*2,3 + *M*1,3*N*3,3 *<sup>M</sup>*3,1*N*1,1 <sup>−</sup> *<sup>M</sup>*2,1*N*2,1 <sup>+</sup> *<sup>M</sup>*1,1*N*3,1 *<sup>D</sup>* <sup>−</sup> *<sup>M</sup>*3,3*N*1,1 <sup>+</sup> *<sup>M</sup>*2,3*N*2,1 <sup>−</sup> *<sup>M</sup>*1,3*N*3,1 , *ρ*ˆ <sup>4</sup> <sup>=</sup> <sup>1</sup> *D* −*M*3,2*N*1,2 + *M*2,2*N*2,2 − *M*1,2*N*3,2 −*M*3,2*N*1,3 + *M*2,2*N*2,3 − *M*1,2*N*3,3 *<sup>M</sup>*3,1*N*1,2 <sup>−</sup> *<sup>M</sup>*2,1*N*2,2 <sup>+</sup> *<sup>M</sup>*1,1*N*3,2 *<sup>D</sup>* <sup>+</sup> *<sup>M</sup>*3,2*N*1,2 <sup>−</sup> *<sup>M</sup>*2,2*N*2,2 <sup>+</sup> *<sup>M</sup>*1,2*N*3,2 , *ρ*ˆ <sup>5</sup> <sup>=</sup> <sup>1</sup> *D M*3,3*N*1,1 − *M*2,3*N*2,1 + *M*1,3*N*3,1 *M*3,3*N*1,2 − *M*2,3*N*2,2 + *M*1,3*N*3,2 <sup>−</sup>*M*3,2*N*1,1 <sup>+</sup> *<sup>M</sup>*2,2*N*2,1 <sup>−</sup> *<sup>M</sup>*1,2*N*3,1 *<sup>D</sup>* <sup>−</sup> *<sup>M</sup>*3,3*N*1,1 <sup>+</sup> *<sup>M</sup>*2,3*N*2,1 <sup>−</sup> *<sup>M</sup>*1,3*N*3,1 , *ρ*ˆ <sup>6</sup> <sup>=</sup> <sup>1</sup> *D D* − *M*3,1*N*1,3 + *M*2,1*N*2,3 − *M*1,1*N*3,3 −*M*3,2*N*1,3 + *M*2,2*N*2,3 − *M*1,2*N*3,3 *<sup>M</sup>*3,1*N*1,2 <sup>−</sup> *<sup>M</sup>*2,1*N*2,2 <sup>+</sup> *<sup>M</sup>*1,1*N*3,2 *<sup>M</sup>*3,1*N*1,3 <sup>−</sup> *<sup>M</sup>*2,1*N*2,3 <sup>+</sup> *<sup>M</sup>*1,1*N*3,3 , (2)

where *Njk* = (*ρ*ˆ*U*ˆ )*jk*, *D* is the determinant of *U*ˆ , and *Mjk* are the components of the minors of matrix *<sup>U</sup>*<sup>ˆ</sup> , i.e., its elements are the determinants after eliminating the (<sup>4</sup> <sup>−</sup> *<sup>j</sup>*)th row and (<sup>4</sup> <sup>−</sup> *<sup>k</sup>*)th column of *U*ˆ . The transformed states are characterized into different sets by their corresponding transformed entropies *S <sup>k</sup>* = −Tr *ρ <sup>k</sup>* ln *ρ <sup>k</sup>*. We emphasize that the resulting qubit density matrices are associated, in general, with a nonunitary evolution of the original qubits. This fact establishes a new mechanism to obtain the open-like system evolution in a noncomposite qutrit system. Additionally, this procedure can be extended to any qudit system, in view of the general definition of the qubit density matrices obtained from a qudit system [10].

In [9], we discussed that a two-qubit density matrix with one of its rows and columns equal to zero describes separable states, if one of the off-diagonal terms is equal to zero, for example, the state:

$$
\hat{\rho} = \begin{pmatrix}
\rho\_{11} & \rho\_{12} & \rho\_{13} & 0 \\
\rho\_{21} & \rho\_{22} & \rho\_{23} & 0 \\
\rho\_{31} & \rho\_{32} & \rho\_{33} & 0 \\
0 & 0 & 0 & 0
\end{pmatrix},
$$

is separable iff *ρ*<sup>23</sup> = 0. To show this, one can consider the previous density matrix to be in the standard two-qubit representation |00, |01, |10, and |11. It can be seen that the partial transpose operation [22] implies the change *ρ*<sup>12</sup> ↔ *ρ*21, and for this reason, the eigenvalues of *ρ*ˆ with *ρ*<sup>23</sup> = 0 are equal to the eigenvalues of its partial transpose. As the partial transpose is a nonnegative operator, then the system is separable. The separability implies the invariance of the partial density matrices under local unitary transformations. As this two-qubit density matrix has a pair of row-column with a diagonal term equal to zero, the correspondence with a qutrit density matrix can be made. On the other hand, the correspondence between two-qubit local unitary transformations and qutrit unitary transformations can be made in the same way, e.g., by eliminating one row and one column of the two-qubit local transformation. This procedure allows us to define different unitary transformations that almost leave the qubits in Expression (1) invariant.

#### **3. Phase Damping and Spontaneous-Emission Channels**

It is known that the interaction of a qubit system with an environment leads to several physical phenomena such as dissipation and decoherence in the qubit subsystem; an example of these interactions is the phase damping channel. In this channel, the evolution of the qubit plus environment (|···*q*|···*e*) is given by a unitary transformation *<sup>T</sup>*ˆ, which acts differently if the qubit is in the ground or excited state, according to the following rules: *<sup>T</sup>*ˆ(|0*q*|0*e*) = '<sup>1</sup> <sup>−</sup> *<sup>p</sup>*|0*q*|0*<sup>e</sup>* <sup>+</sup> <sup>√</sup>*p*|0*q*|1*<sup>e</sup>* and *<sup>T</sup>*ˆ(|1*q*|0*e*) = '<sup>1</sup> <sup>−</sup> *<sup>p</sup>*|0*q*|0*<sup>e</sup>* <sup>+</sup> <sup>√</sup>*p*|0*q*|2*<sup>e</sup>* with *<sup>p</sup>* being a probability, i.e., the environment subsystem goes to a superposition of the states (|0*e*, |1*e*), or to (|0*e*, |2*e*), if the environment is in |0*e*, or |1*e*, respectively [15,23]. This two-qubit unitary transformations result in a nonunitary change when the partial trace over the environment subsystem is taken:

$$
\begin{pmatrix} 1 - \rho\_{22} & \rho\_{12} \\ \rho\_{12}^\* & \rho\_{22} \end{pmatrix} \rightarrow \begin{pmatrix} 1 - \rho\_{22} & \rho\_{12}(1 - p) \\ \rho\_{12}^\*(1 - p) & \rho\_{22} \end{pmatrix} \cdot \mathbf{1}
$$

When the map is applied a very large number of times (→ ∞), it is straightforward that the initial state tends to the completely decoherent state:

$$
\begin{pmatrix} 1 - \rho\_{22} & \rho\_{12} \\ \rho\_{12}^\* & \rho\_{22} \end{pmatrix} \rightarrow \begin{pmatrix} 1 - \rho\_{22} & 0 \\ 0 & \rho\_{22} \end{pmatrix} \wedge$$

with an exponential convergence.

The other example is the spontaneous-emission (also called the amplitude-damping) quantum channel. In this channel, the dynamics of the qubit system plus the environment is determined by a unitary transform *<sup>T</sup>*ˆ, which only acts if the qubit system is in the excited state <sup>|</sup>1*q*, according to the following rules: *<sup>T</sup>*ˆ(|0*q*|0*e*) = <sup>|</sup>0*q*|0*<sup>e</sup>* and *<sup>T</sup>*ˆ(|1*q*|0*e*) = '<sup>1</sup> <sup>−</sup> *<sup>p</sup>*|1*q*|0*<sup>e</sup>* <sup>+</sup> <sup>√</sup>*p*|0*q*|1*e*, where *<sup>p</sup>* is the probability [15,23]. This channel then defines a nonunitary evolution over the qubit subsystem, which transforms the qubit density matrix as follows:

$$
\begin{pmatrix} 1 - \rho\_{22} & \rho\_{12} \\ \rho\_{12}^\* & \rho\_{22} \end{pmatrix} \rightarrow \begin{pmatrix} 1 - (1 - p)\rho\_{22} & \rho\_{12}\sqrt{1 - p} \\ \rho\_{12}^\*\sqrt{1 - p} & (1 - p)\rho\_{22} \end{pmatrix} \cdot \vec{k}
$$

If this channel is applied a very large number of times (→ ∞), the density matrix converges to a ground state, i.e.,

$$
\begin{pmatrix} 1 - \rho\_{22} & \rho\_{12} \\ \rho\_{12}^\* & \rho\_{22} \end{pmatrix} \to \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}
$$

.

In addition to these examples, there exists another type of quantum channel defined in the theory of interaction between a quantum system and an environment, which can be considered [15,23].

It is possible to demonstrate that phase damping and spontaneous-emission quantum channels for qubits *ρ*ˆ1, ... , *ρ*ˆ6 in Equation (1) can be obtained by the use of particular unitary transformations of a qutrit state *ρ*ˆ. To justify this, we assumed a two-qubit quantum system where one of the levels cannot be populated, i.e., the 4 × 4 density matrix has an eigenvalue equal to zero, e.g.,

$$
\boldsymbol{\dot{\rho}} = \begin{pmatrix}
\rho\_{11} & \rho\_{12} & \rho\_{13} & 0 \\
\rho\_{21} & \rho\_{22} & 0 & 0 \\
\rho\_{31} & 0 & \rho\_{33} & 0 \\
0 & 0 & 0 & 0
\end{pmatrix};
\tag{3}
$$

it is clear that this density matrix is separable since *ρ*<sup>23</sup> = *ρ*<sup>∗</sup> <sup>32</sup> = 0. The partial density matrices can be operated locally by unitary transformations of the form *u*ˆ1 ⊗ *u*ˆ2. When only one of the qubits is operated, i.e., when the unitary matrix corresponds to a controlled operation [15]: *u*ˆ1 = ˆ*I* or *u*ˆ2 = ˆ*I*. If *u*ˆ2 = ˆ*I*, then the unitary transformation only operates over the second qubit,

$$
\hat{u} = \begin{pmatrix} u\_{11} & u\_{12} & 0 & 0 \\ u\_{21} & u\_{22} & 0 & 0 \\ 0 & 0 & u\_{11} & u\_{12} \\ 0 & 0 & u\_{21} & u\_{22} \end{pmatrix} . \tag{4}
$$

By means of this type of unitary matrix, one can define an operation in the qutrit system that approximately only affects *ρ*ˆ2. This is done by ignoring the fourth row and the fourth column of (3); the resulting qutrit state is then operated by the unitary matrix resulting from the elimination of the fourth row and the fourth column of Equation (4). For the operator to be still unitary, the (3, 3) entry must be replaced by one. Following these and other analogous arguments, we study the application of the unitary transforms:

$$\begin{aligned} \bullet \bullet = \left( \begin{array}{cccc} u\_{11} & u\_{12} & 0 \\ u\_{21} & u\_{22} & 0 \\ 0 & 0 & 1 \end{array} \right) , \quad \bullet \bullet = \left( \begin{array}{cccc} u\_{11} & 0 & u\_{12} \\ 0 & 1 & 0 \\ u\_{21} & 0 & u\_{22} \end{array} \right) , \quad \bullet \bullet = \left( \begin{array}{cccc} 1 & 0 & 0 \\ 0 & u\_{11} & u\_{12} \\ 0 & u\_{21} & u\_{22} \end{array} \right) \end{aligned} \tag{5}$$

on the qutrit density matrices:

$$\begin{aligned} \label{eq:SD1} \mathcal{O}\_1 = \begin{pmatrix} \rho\_{11} & \rho\_{12} & \rho\_{13} \\ \rho\_{21} & \rho\_{22} & 0 \\ \rho\_{31} & 0 & \rho\_{33} \end{pmatrix}, \quad \begin{array}{c} \mathcal{O}\_2 = \begin{pmatrix} \rho\_{11} & 0 & \rho\_{13} \\ 0 & \rho\_{22} & \rho\_{23} \\ \rho\_{31} & \rho\_{32} & \rho\_{33} \end{pmatrix}, \quad \mathcal{O}\_3 = \begin{pmatrix} \rho\_{11} & \rho\_{12} & 0 \\ \rho\_{21} & \rho\_{22} & \rho\_{23} \\ 0 & \rho\_{23} & \rho\_{33} \end{pmatrix}. \end{aligned} \tag{6}$$

The unitary transformations in Equation (5) can be enacted on any of the density matrices in Equation (6), which define a nonunitary transformation of the qubits defined in Equation (1). These qubit transformations are found by the substitution of Equations (5) and (6) into Equation (2), e.g., the unitary transformation *U*ˆ † <sup>1</sup>*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup> results in the following transformations of the qubits:

*ρ*ˆ <sup>1</sup> = 1 − *ρ*<sup>33</sup> *ρ*<sup>13</sup> *u*<sup>∗</sup> 11 *<sup>ρ</sup>*<sup>31</sup> *<sup>u</sup>*<sup>11</sup> *<sup>ρ</sup>*<sup>33</sup> , *ρ*ˆ <sup>2</sup> = 1 − *u*<sup>∗</sup> <sup>12</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>12</sup> <sup>−</sup> *<sup>u</sup>*<sup>∗</sup> <sup>22</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>22</sup> *<sup>u</sup>*<sup>∗</sup> <sup>11</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>12</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>21</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>22</sup> *u*∗ <sup>12</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>11</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>22</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>21</sup> *<sup>u</sup>*<sup>∗</sup> <sup>12</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>12</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>22</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>22</sup> , *ρ*ˆ <sup>3</sup> = *u*∗ <sup>11</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>11</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>21</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>21</sup> *<sup>ρ</sup>*<sup>13</sup> *<sup>u</sup>*<sup>∗</sup> 11 *ρ*<sup>31</sup> *u*<sup>11</sup> 1 − *u*<sup>∗</sup> <sup>11</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>11</sup> <sup>−</sup> *<sup>u</sup>*<sup>∗</sup> <sup>21</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>21</sup> , *ρ*ˆ <sup>4</sup> = *u*∗ <sup>12</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>12</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>22</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>22</sup> *<sup>ρ</sup>*13*u*<sup>∗</sup> 12 *ρ*31*u*<sup>12</sup> 1 − *u*<sup>∗</sup> <sup>12</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>12</sup> <sup>−</sup> *<sup>u</sup>*<sup>∗</sup> <sup>22</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>22</sup> , *ρ*ˆ <sup>5</sup> = *u*∗ <sup>11</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>11</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>21</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>21</sup> *<sup>u</sup>*<sup>∗</sup> <sup>11</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>12</sup> + *<sup>u</sup>*<sup>∗</sup> <sup>21</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>22</sup> *u*∗ <sup>12</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)11(*r*11*u*<sup>11</sup> + *<sup>r</sup>*12*u*21) + *<sup>u</sup>*<sup>∗</sup> <sup>22</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>21</sup> <sup>1</sup> <sup>−</sup> *<sup>u</sup>*<sup>∗</sup> <sup>11</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>11</sup> <sup>−</sup> *<sup>u</sup>*<sup>∗</sup> <sup>21</sup>(*σ*ˆ1*U*<sup>ˆ</sup> <sup>1</sup>)<sup>21</sup> , *ρ*ˆ <sup>6</sup> = 1 − *ρ*<sup>33</sup> *ρ*<sup>13</sup> *u*<sup>∗</sup> 12 *<sup>ρ</sup>*<sup>31</sup> *<sup>u</sup>*<sup>12</sup> *<sup>ρ</sup>*<sup>33</sup> .

*Entropy* **2019**, *21*, 736

From these results, one can notice that the transformed qubits *ρ*ˆ <sup>1</sup> and *ρ*ˆ <sup>6</sup> correspond to the phase damping channel of *ρ*ˆ1 with different damping parameters. Furthermore, the qubit states *ρ*ˆ <sup>2</sup>, *ρ*ˆ 5 can be seen as quasi-unitary transformations of the initial states *ρ*ˆ2, *ρ*ˆ5, respectively. In a similar way, one can obtain all the possible unitary transformations of the density matrices in Equation (6). These transformations lead to the identification of two types of quantum channels: the phase damping and a quasi-unitary operation described below.

The unitary transformation over the density matrices *σ*ˆ1, *σ*ˆ2, and *σ*ˆ3 results in a change over their associated qubits *ρ*ˆ1, ... , *ρ*ˆ6, to *ρ*ˆ <sup>1</sup>, ... , *ρ*ˆ <sup>6</sup>, which denote the qubits after the transformation. We have found the following interesting expressions:

*U*ˆ † 1*σ*ˆ <sup>1</sup>*U*<sup>ˆ</sup> <sup>1</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>1</sup> = 1 − *ρ*<sup>33</sup> *u*<sup>∗</sup> <sup>11</sup>*ρ*<sup>13</sup> *<sup>u</sup>*11*ρ*<sup>31</sup> *<sup>ρ</sup>*<sup>33</sup> , *ρ*ˆ <sup>6</sup> = 1 − *ρ*<sup>33</sup> *u*<sup>∗</sup> <sup>12</sup>*ρ*<sup>13</sup> *<sup>u</sup>*12*ρ*<sup>31</sup> *<sup>ρ</sup>*<sup>33</sup> ; *U*ˆ † 2*σ*ˆ <sup>1</sup>*U*<sup>ˆ</sup> <sup>2</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>2</sup> = 1 − *ρ*<sup>22</sup> *u*<sup>∗</sup> <sup>11</sup>*ρ*<sup>12</sup> *<sup>u</sup>*11*ρ*<sup>21</sup> *<sup>ρ</sup>*<sup>22</sup> , *ρ*ˆ <sup>4</sup> = 1 − *ρ*<sup>33</sup> *u*<sup>∗</sup> <sup>12</sup>*ρ*<sup>12</sup> *<sup>u</sup>*12*ρ*<sup>21</sup> *<sup>ρ</sup>*<sup>33</sup> ; *U*ˆ † <sup>2</sup>*σ*ˆ2*U*<sup>ˆ</sup> <sup>2</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>2</sup> = 1 − *ρ*<sup>22</sup> *u*<sup>∗</sup> <sup>21</sup>*ρ*<sup>32</sup> *<sup>u</sup>*21*ρ*<sup>23</sup> *<sup>ρ</sup>*<sup>22</sup> , *ρ*ˆ <sup>4</sup> = *ρ*<sup>22</sup> *u*22*ρ*<sup>23</sup> *u*∗ <sup>22</sup>*ρ*<sup>32</sup> <sup>1</sup> <sup>−</sup> *<sup>ρ</sup>*<sup>22</sup> ; (7) *U*ˆ † <sup>3</sup>*σ*ˆ2*U*<sup>ˆ</sup> <sup>3</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>3</sup> = *ρ*<sup>11</sup> *u*22*ρ*<sup>13</sup> *u*∗ <sup>22</sup>*ρ*<sup>31</sup> <sup>1</sup> <sup>−</sup> *<sup>ρ</sup>*<sup>11</sup> , *ρ*ˆ <sup>5</sup> = *ρ*<sup>11</sup> *u*21*ρ*<sup>13</sup> *u*∗ <sup>21</sup>*ρ*<sup>31</sup> <sup>1</sup> <sup>−</sup> *<sup>ρ</sup>*<sup>11</sup> ; *U*ˆ † <sup>1</sup>*σ*ˆ3*U*<sup>ˆ</sup> <sup>1</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>1</sup> = 1 − *ρ*<sup>33</sup> *u*<sup>∗</sup> <sup>21</sup>*ρ*<sup>23</sup> *<sup>u</sup>*21*ρ*<sup>32</sup> *<sup>ρ</sup>*<sup>33</sup> , *ρ*ˆ <sup>6</sup> = 1 − *ρ*<sup>33</sup> *u*<sup>∗</sup> <sup>22</sup>*ρ*<sup>23</sup> *<sup>u</sup>*22*ρ*<sup>32</sup> *<sup>ρ</sup>*<sup>33</sup> ; *U*ˆ † <sup>3</sup>*σ*ˆ3*U*<sup>ˆ</sup> <sup>3</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>3</sup> = *ρ*<sup>11</sup> *u*12*ρ*<sup>12</sup> *u*∗ <sup>12</sup>*ρ*<sup>21</sup> <sup>1</sup> <sup>−</sup> *<sup>ρ</sup>*<sup>11</sup> , *ρ*ˆ <sup>5</sup> = *ρ*<sup>11</sup> *u*11*ρ*<sup>12</sup> *u*∗ <sup>11</sup>*ρ*<sup>21</sup> <sup>1</sup> <sup>−</sup> *<sup>ρ</sup>*<sup>11</sup> .

In most of the cases, the resulting qubits *ρ*ˆ *<sup>j</sup>* correspond to the phase damping quantum channel of *ρ*ˆ*j*, as can be seen in Expression (8). In this channel, the probability amplitudes given by the original off-diagonal terms of the qubits are multiplied by a number. The damping parameters are associated with different entries of the unitary transformation *ujk*, which in general are complex numbers. When the unitary transformation correspond to a real matrix, then the expression for the standard phase damping map is obtained. As you can see in Equation (8), in some cases, the unitary transformations leads to the quantum channel of another qubit, e.g., after the application of *U*ˆ <sup>1</sup> to *σ*ˆ1, the qubit *ρ*ˆ <sup>6</sup> is the phase damping channel of *ρ*ˆ1. Furthermore, in some other cases, the obtained density matrices correspond to transformations similar to the phase damping channel of matrices outside the ones in Equation (1), e.g., *ρ*ˆ <sup>4</sup> after the application of *<sup>U</sup>*<sup>ˆ</sup> <sup>2</sup> to *<sup>σ</sup>*ˆ1. Although these matrices seem unrelated, they have the same form as the phase damping channel. In the case of *U*ˆ being a rotation matrix with a time-dependent angle *θ* = *ωt*, the original qubit states can be recovered at the time *t* = 2*πl*/*ω*, *l* = 0, 1, 2, . . ..

The unitary transformations (*U*ˆ 1, *U*ˆ 2, *U*ˆ 3) previously described can also lead to quasi-unitary transformations of the qubits. In particular, for the unitary transformation *U*ˆ † <sup>1</sup>*σ*ˆ1*U*<sup>ˆ</sup> 1, one gets the quasi-unitary transformations:

$$\begin{array}{rcl}\hat{\rho}'\_2 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_2 \hat{\mathcal{U}} + \rho\_{33} \left( \begin{array}{c} |u\_{12}|^2 & -u\_{11}^\* u\_{12} \\ -u\_{11} u\_{12}^\* & -|u\_{12}|^2 \end{array} \right), \\\ \hat{\rho}'\_5 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_5 \hat{\mathcal{U}} + \rho\_{33} \left( \begin{array}{c} -|u\_{21}|^2 & -u\_{21}^\* u\_{22} \\ -u\_{21} u\_{22}^\* & |u\_{21}|^2 \end{array} \right), \end{array} \tag{8}$$

with <sup>U</sup><sup>ˆ</sup> <sup>=</sup> *u*<sup>11</sup> *u*<sup>12</sup> *<sup>u</sup>*<sup>21</sup> *<sup>u</sup>*<sup>22</sup> being a two-dimensional unitary transformation. For the other qubits, one can also define quasi-unitary transformations as follows:

(a) From the qutrit unitary transformation *U*ˆ † <sup>1</sup>*σ*ˆ3*U*<sup>ˆ</sup> 1,

$$\begin{array}{rcl} \rho'\_2 &=& \mathcal{U}^\dagger \rho\_2 \mathcal{U} + \rho\_{33} \left( \begin{array}{cc} -|u\_{12}|^2 & u\_{11}^\* u\_{12} \\ u\_{11} u\_{12}^\* & |u\_{12}|^2 \end{array} \right) \\\ \rho'\_5 &=& \mathcal{U}^\dagger \rho\_5 \mathcal{U} + \rho\_{33} \left( \begin{array}{cc} |u\_{21}|^2 & u\_{21}^\* u\_{22} \\ u\_{21} u\_{22}^\* & -|u\_{21}|^2 \end{array} \right) . \end{array} \tag{9}$$

(b) For the transformation *U*ˆ † <sup>2</sup>*σ*ˆ1*U*<sup>ˆ</sup> 2,

$$\begin{array}{rcl} \hat{\rho}'\_1 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_1 \hat{\mathcal{U}} + \rho\_{22} \left( \begin{array}{c} |u\_{12}|^2 & -u\_{11}^\* u\_{12} \\ -u\_{11} u\_{12}^\* & -|u\_{12}|^2 \end{array} \right), \\\ \hat{\rho}'\_3 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_3 \hat{\mathcal{U}} + \rho\_{22} \left( \begin{array}{c} -|u\_{21}|^2 & -u\_{21}^\* u\_{22} \\ -u\_{21} u\_{22}^\* & |u\_{21}|^2 \end{array} \right). \end{array} \tag{10}$$

(c) For the transformation *U*ˆ † <sup>2</sup>*σ*ˆ2*U*<sup>ˆ</sup> 2,

$$\begin{array}{rcl} \hat{\rho}'\_1 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_1 \hat{\mathcal{U}} + \rho\_{22} \left( \begin{array}{c} |u\_{12}|^2 & -u\_{11}^\* u\_{12} \\ -u\_{11} u\_{12}^\* & -|u\_{12}|^2 \end{array} \right) \\ \hat{\rho}'\_3 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_3 \hat{\mathcal{U}} + \rho\_{22} \left( \begin{array}{c} -|u\_{21}|^2 & -u\_{21}^\* u\_{22} \\ -u\_{21} u\_{22}^\* & |u\_{21}|^2 \end{array} \right) \end{array} \tag{11}$$

(d) From *U*ˆ † <sup>3</sup>*σ*ˆ2*U*<sup>ˆ</sup> 3,

$$\begin{array}{rcl}\hat{\rho}'\_4 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_4 \hat{\mathcal{U}} + \rho\_{11} \left( \begin{array}{cc} -|u\_{21}|^2 & -u\_{12}^\* u\_{22} \\ -u\_{21} u\_{22}^\* & |u\_{21}|^2 \end{array} \right), \\\ \hat{\rho}'\_6 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_6 \hat{\mathcal{U}} + \rho\_{11} \left( \begin{array}{cc} -|u\_{12}|^2 & u\_{11}^\* u\_{12} \\ |u\_{11} u\_{12}^\* & |u\_{12}|^2 \end{array} \right). \end{array} \tag{12}$$

(e) Finally, for *U*ˆ † <sup>3</sup>*σ*ˆ3*U*<sup>ˆ</sup> 3,

$$\begin{array}{rcl}\hat{\rho}'\_4 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_4 \hat{\mathcal{U}} + \rho\_{11} \begin{pmatrix} -|u\_{21}|^2 & -u\_{12}^\* u\_{22} \\ -u\_{21} u\_{22}^\* & |u\_{21}|^2 \end{pmatrix},\\\hat{\rho}'\_6 &=& \hat{\mathcal{U}}^\dagger \hat{\rho}\_6 \hat{\mathcal{U}} + \rho\_{11} \begin{pmatrix} -|u\_{12}|^2 & u\_{11}^\* u\_{12} \\ |u\_{11} u\_{12}^\* & |u\_{12}|^2 \end{pmatrix},\end{array} \tag{13}$$

For all the cases, <sup>U</sup><sup>ˆ</sup> is a two-dimensional unitary transformation.

As in the phase-damping case, one can think of a rotation matrix with a time-dependent angle *θ* = *ωt* as the unitary operation, i.e.,

$$
\hat{\mathcal{U}} = \begin{pmatrix}
\cos(\omega t) & -\sin(\omega t) \\
\sin(\omega t) & \cos(\omega t)
\end{pmatrix} \hat{\mathcal{I}}
$$

which, in the case where *t* ≈ 0, results in the following transformations:

$$\boldsymbol{\delta}\_{\dot{j}}^{\prime} = \boldsymbol{\hat{\mathcal{U}}}^{\dagger} \boldsymbol{\hat{\rho}}\_{\dot{j}} \boldsymbol{\hat{\mathcal{U}}} - \rho\_{kk} \boldsymbol{\omega} \,\mathrm{t} \,\boldsymbol{\hat{\sigma}}\_{\dot{x}} + \mathcal{O}\left(\boldsymbol{t}^{2}\right), \tag{14}$$

where *σ*ˆ*<sup>x</sup>* is the Pauli matrix and *ρkk* is a diagonal component of *ρ*ˆ, which depends on *j*. Its value is *k* = 2 for *j* = 1, 3, *k* = 3 for *j* = 2, 5, and *k* = 1 for *j* = 4, 6. It is necessary to point out that, for *ρ*ˆ 5 associated with *U*ˆ † <sup>1</sup>*σ*ˆ3*U*<sup>ˆ</sup> 1, we need to replace *<sup>ρ</sup>*<sup>33</sup> with <sup>−</sup>*ρ*<sup>33</sup> in Equation (14).

In the case where the density matrices correspond to states, where one of the accessible levels is not occupied, i.e.,

$$
\begin{aligned}
\mathfrak{d}\_4 = \begin{pmatrix}
\rho\_{11} & \rho\_{12} & 0 \\
\rho\_{21} & \rho\_{22} & 0 \\
0 & 0 & 0
\end{pmatrix}, \quad \mathfrak{d}\_5 = \begin{pmatrix}
\rho\_{11} & 0 & \rho\_{13} \\
0 & 0 & 0 \\
\rho\_{31} & 0 & \rho\_{33}
\end{pmatrix}, \quad \mathfrak{d}\_6 = \begin{pmatrix}
0 & 0 & 0 \\
0 & \rho\_{22} & \rho\_{23} \\
0 & \rho\_{32} & \rho\_{33}
\end{pmatrix},
\end{aligned} \tag{15}
$$

we obtain the expressions:

*U*ˆ † <sup>2</sup>*σ*ˆ4*U*<sup>ˆ</sup> <sup>2</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>5</sup> = *ρ*11|*u*11| <sup>2</sup> *ρ*12*u*<sup>∗</sup> 11 *ρ*21*u*<sup>11</sup> 1 − *ρ*11|*u*11| 2 , *ρ*ˆ <sup>6</sup> = 1 − *ρ*11|*u*12| <sup>2</sup> *ρ*21*u*<sup>12</sup> *ρ*12*u*<sup>∗</sup> <sup>12</sup> *ρ*11|*u*12| 2 , *U*ˆ † <sup>3</sup>*σ*ˆ4*U*<sup>ˆ</sup> <sup>3</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>1</sup> = 1 − *ρ*22|*u*12| <sup>2</sup> *ρ*12*u*<sup>12</sup> *ρ*21*u*<sup>∗</sup> <sup>12</sup> *ρ*22|*u*12| 2 , *ρ*ˆ <sup>2</sup> = 1 − *ρ*22|*u*11| <sup>2</sup> *ρ*12*u*<sup>11</sup> *ρ*21*u*<sup>∗</sup> <sup>11</sup> *ρ*22|*u*11| 2 , *U*ˆ † <sup>1</sup>*σ*ˆ5*U*<sup>ˆ</sup> <sup>1</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>3</sup> = *ρ*11|*u*11| <sup>2</sup> *ρ*13*u*<sup>∗</sup> 11 *ρ*31*u*<sup>11</sup> 1 − *ρ*11|*u*11| 2 , *ρ*ˆ <sup>4</sup> = *ρ*11|*u*12| <sup>2</sup> *ρ*13*u*<sup>∗</sup> 12 *ρ*31*u*<sup>12</sup> 1 − *ρ*11|*u*12| 2 , (16) *U*ˆ † <sup>3</sup>*σ*ˆ5*U*<sup>ˆ</sup> <sup>3</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>1</sup> = 1 − *ρ*33|*u*22| <sup>2</sup> *ρ*13*u*<sup>22</sup> *ρ*31*u*<sup>∗</sup> <sup>22</sup> *ρ*33|*u*22| 2 , *ρ*ˆ <sup>2</sup> = 1 − *ρ*33|*u*21| <sup>2</sup> *ρ*13*u*<sup>21</sup> *ρ*31*u*<sup>∗</sup> <sup>21</sup> *ρ*33|*u*21| 2 , *U*ˆ † <sup>1</sup>*σ*ˆ6*U*<sup>ˆ</sup> <sup>1</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>3</sup> = *ρ*22|*u*21| <sup>2</sup> *ρ*23*u*<sup>∗</sup> 21 *ρ*32*u*<sup>21</sup> 1 − *ρ*22|*u*21| 2 , *ρ*ˆ <sup>4</sup> = *ρ*22|*u*22| <sup>2</sup> *ρ*23*u*<sup>∗</sup> 22 *ρ*32*u*<sup>22</sup> 1 − *ρ*22|*u*22| 2 , *U*ˆ † <sup>2</sup>*σ*ˆ6*U*<sup>ˆ</sup> <sup>2</sup> <sup>⇒</sup> *<sup>ρ</sup>*<sup>ˆ</sup> <sup>5</sup> = *ρ*33|*u*21| <sup>2</sup> *ρ*32*u*<sup>∗</sup> 21 *ρ*23*u*<sup>21</sup> 1 − *ρ*33|*u*21| 2 , *ρ*ˆ <sup>6</sup> = 1 − *ρ*33|*u*22| <sup>2</sup> *ρ*23*u*<sup>22</sup> *ρ*32*u*<sup>∗</sup> <sup>22</sup> *ρ*33|*u*22| 2 .

These transformations in many of the cases can represent the spontaneous-emission quantum channel. As in the other examples studied above, when the unitary matrices are rotated by angle *θ* = *ωt*, the original qubit systems can be recovered at times *t* = 2*πl*/*ω*; *l* = 0, 1, 2, .... It is important to mention that the states represented by Equation (15) correspond to three-level systems, where one of the levels is a dark state, and then only two of the levels can be populated, which have been experimentally obtained [24]. These kinds of systems have been of relevance as they can be created by two-photon processes in a three-level system [25] or by the adiabatic variation of the Rabi frequencies associated with the transitions between the three states [26]. For example, to obtain the state *σ*ˆ4, one can think of an atomic Λ-type three-level system (|1, |2, |3), which interacts with an environment [26]; see Figure 1. The Hamiltonian associated with this system can be written in the form:

$$
\hat{H} = \begin{pmatrix}
\omega\_1 & 0 & \omega\_{13} \\
0 & \omega\_2 & \omega\_{23} \\
\omega\_{13} & \omega\_{23} & 0
\end{pmatrix},
$$

where *ω*1,2 are the energies of the states |1, |2, respectively. By considering the energy of the ground state |3 equal to zero, *ω*<sup>13</sup> and *ω*<sup>23</sup> are the transition energies. Taking the zero energy in the ground state <sup>|</sup>3, we can make the replacements *<sup>ω</sup>*<sup>13</sup> <sup>→</sup> *<sup>ω</sup>*<sup>1</sup> *<sup>e</sup>*−*iω*1*<sup>t</sup>* and *<sup>ω</sup>*<sup>23</sup> <sup>→</sup> *<sup>ω</sup>*<sup>2</sup> *<sup>e</sup>*−*iω*2*<sup>t</sup>* . The time evolution of the density matrix can be obtained by the expression:

$$\frac{d}{dt}\mathfrak{\dot{\rho}} = i[\mathfrak{\dot{\rho}}, \hat{H}] + \mathfrak{\dot{\rho}}',\tag{17}$$

where the matrix *ρ*ˆ is given by the interaction of the original density matrix with the environment:

$$
\rho' = \begin{pmatrix}
\gamma\_{31}\rho\_{33} & -\gamma'\rho\_{12} & -\gamma\_1\rho\_{13} \\
\end{pmatrix}
$$

where the parameters *γ*31, *γ*32, and *γ* are the spontaneous-emission rates, which must satisfy *γ* = *γ*<sup>31</sup> + *γ*32, and the relaxation terms for the coherence components are named *γ*<sup>1</sup> and *γ*2, which also satisfy *γ* = *γ*<sup>1</sup> + *γ*2. The resulting differential Equation (17) can be reduced by considering that the variation of the parameters *ρ*13, *ρ*23, and *ρ*<sup>33</sup> over time is smaller compared to the spontaneous emission and decoherence terms *γ*<sup>31</sup> and *γ*32; this is called the adiabatic hypothesis. Under this hypothesis, it is possible to obtain a state with *ρ*<sup>13</sup> = *ρ*<sup>23</sup> = *ρ*<sup>33</sup> = 0, as the solution of the evolution of the density matrix *σ*ˆ4 discussed above.

Another way to obtain these types of systems is the case where the environmental interaction is neglected, i.e., *ρ*ˆ = 0 in Equation (17). The corresponding Schrödinger equation is *i <sup>d</sup>*|*ψ dt* <sup>=</sup> *<sup>H</sup>*<sup>ˆ</sup> <sup>|</sup>*ψ*, with <sup>|</sup>*ψ* <sup>=</sup> *<sup>a</sup>*1(*t*)*e*−*iω*1*<sup>t</sup>* <sup>|</sup>1 <sup>+</sup> *<sup>a</sup>*2(*t*)*e*−*iω*2*<sup>t</sup>* |2 + *a*3(*t*)|3, which in view of the initial conditions *<sup>a</sup>*1(0) = <sup>√</sup> *<sup>ω</sup>*<sup>2</sup> *ω*<sup>2</sup> 1+*ω*<sup>2</sup> 2 , *<sup>a</sup>*2(0) = <sup>−</sup> <sup>√</sup> *<sup>ω</sup>*<sup>1</sup> *ω*<sup>2</sup> 1+*ω*<sup>2</sup> 2 , *a*3(0) = 0 leads to the solution:

$$a\_1(t) = \frac{\omega\_2}{\sqrt{\omega\_1^2 + \omega\_2^2}}, \quad a\_2(t) = -\frac{\omega\_1}{\sqrt{\omega\_1^2 + \omega\_2^2}}; \quad a\_3(t) = 0, 1$$

so the level |3 is never populated.

The density matrices *σ*ˆ5 and *σ*ˆ6 can be obtained by means of analogous procedures applied to the V and Ξ configurations of the three-level system depicted in Figure 1.

**Figure 1.** State configurations for the V- (**left**), the Λ- (**center**), and the Ξ-level (**right**) systems.

It is also important to mention that the unitary transformations defined by the matrices *U*ˆ 1, *U*ˆ 2, and *U*ˆ <sup>3</sup> in Equation (5) can be generated experimentally by different proposed mechanisms, such as sliding mode control [27], adiabatic passage [28–30], and the robust control scheme [31,32]. We want to emphasize that the resulting quasi-unitary evolutions and the different quantum channels obtained in our work can have applications in quantum computing and quantum information theories. We think so because the quasi-unitary operations discussed here could be used as approximations to the standard quantum gates, and furthermore, the obtained quantum channels could also be used in the quantum correction algorithms found in the literature.

#### **4. Probability Representation of the Qubit-State Evolution**

In the quantum tomographic approach of qubit states [33,34], the states are identified with tomographic probability distributions. In the case of the minimal number of probability parameters, the density matrix of the qubit (spin-1/2) state reads [6]:

$$\hat{\rho} = \left( \begin{array}{c} p\_3 \\ p\_1 - 1/2 + i(p\_2 - 1/2) \end{array} \right) , \quad \frac{3}{1 - p\_3} \left( \begin{array}{c} p\_1 - 1 \\ p\_3 \end{array} \right) , \quad \sum\_{j=1}^3 \left( p\_j - \frac{1}{4} \right)^2 \le \frac{1}{4} , \tag{18}$$

where 0 ≤ *pk*, ≤ 1 with *k* = 1, 2, 3 are the probabilities to obtain the value +1/2 in the *x*, *y*, *z* axis, respectively. Thus, any qubit state can be identified through the probabilities *p*1, *p*2, and *p*3, i.e., given the density operator, one can get the set *ρ*ˆ ↔ *p*1, *p*2, *p*<sup>3</sup> and vice versa. In the case of qubits (1) associated with the qutrit state, the evolution of the probabilities after the unitary operation of the qutrit is determined by Equation (2). For example, we have a probabilistic representation corresponding to *ρ*ˆ <sup>5</sup> in the first formula of Equation (17), i.e.,

$$p\_3 \to p\_3|u\_{11}|^2, \quad p\_1 - 1/2 - i(p\_2 - 1/2) \to (p\_1 - 1/2 - i(p\_2 - 1/2))u\_{11}^\*. \tag{19}$$

The change of probabilities can be characterized by the evolution of the Tsallis and Shannon entropies. For example, in (19), the unitary matrix parameter *u*<sup>11</sup> determines the evolution of the Shannon entropy related to a coin probability distribution (*p*3, 1 − *p*3) (assume that we have two nonideal classical coins I and II in such a game as coin flipping, coin tossing, or heads (up, ⊕) or tails (down, ), which is the practice of throwing a coin in the air and checking which side is showing when it lands, in order to choose between two alternatives *Pk* or (1 − *Pk*); *k* = 1, 2). This evolution is of the form:

$$S(\hat{\mathcal{U}}) = -p\_3|\boldsymbol{\mu}\_{11}|^2 \ln \left( p\_3|\boldsymbol{\mu}\_{11}|^2 \right) - \left( 1 - p\_3|\boldsymbol{\mu}\_{11}|^2 \right) \ln \left( 1 - p\_3|\boldsymbol{\mu}\_{11}|^2 \right) \dots$$

This entropy, as a function of the unitary evolution applied to the qutrit state, characterizes some aspects of the open dynamics of qubits. We point out that, as for *p*3, there exist other classical entropic characteristics associated with the evolution of *p*<sup>1</sup> and *p*<sup>2</sup> given by Equation (19).

#### **5. Concluding Remarks**

A new mechanism to study the open system evolution of a noncomposite qudit system was established. As an example of the general procedure, we considered a qutrit system. Associated with the qutrit system, one can define different qubit density matrices, which evolve in an open-like way when a unitary transformation is enacted on the qutrit.

The application of the resulting transformations for the qubits within the qutrit was also discussed. The quasi-unitary transformations obtained here might be used as an approximation to quantum gates, whereas the quantum channels could be employed in quantum correction protocols.

Different types of quantum channels can be observed using the qubit decomposition of a qutrit system. In particular, the phase damping and the spontaneous-emission channels were obtained using a unitary transformation acting on specific qutrit density matrices. The phase damping channel was obtained when a unitary transformation of the density matrix with one off-diagonal term equal to zero was performed. A spontaneous-emission channel can be observed by unitary transformations acting over a dark state, i.e., a three-level state where one of the levels cannot be populated.

In addition to these channels, quasi-unitary transformations of the qubit states can be defined. This was also done by the application of a unitary matrix to the generic qutrit state.

The entropy evolution of the tomographic-probability distributions determined by the system of qubits was discussed.

We can extend our analysis to other qudit systems without subsystems since, o an arbitrary spin-*j* density matrix and the spin unitary evolution, one can associate the smaller spin *j* < *j* evolution.

The possible experimental implementation of the procedure was also addressed, given that there exist several proposed ways to generate the unitary transformations such as by sliding mode control [27], adiabatic passage [28–30], or the robust control scheme [31,32].

**Author Contributions:** The original idea was given by J.A.L.-S. All authors contributed equally to the conception, design, and methodology of this study. All authors contributed equally to the analysis of the results and the conclusions. All authors contributed equally to the final writing of the manuscript.

**Funding:** This work was partially supported by DGAPA-UNAM (under Project IN101619).

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


[CrossRef]


c 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
