Under a certain physical process, the state of a given system evolves into another state. The time evolution of a closed system is described by unitary operators. For an open quantum system, which interacts with its environment, the evolution is not unitary any longer. Dynamical processes of open quantum systems are described by a special kind of supermaps called quantum operations. A quantum operation transforms density operators to other density operators while preserving the elementary properties of density operators. In particular, density operators are positive, so a quantum operation needs to preserve positivity. However, it turns out that merely preserving positivity is not sufficient and a much stronger condition is required. Imagine that a system has interacted with its surroundings and established an entanglement with them. Physically, one expects the operation to preserve the properties, positivity in particular, of the density operator of the whole containing the system and surroundings. Essentially, a quantum operation needs to preserve not only the positivity of density operators of a given system but also all density operators of any extended system including the system itself and its surrounding systems. Mathematically, such a condition is satisfied by completely positive supermaps.
package is loaded to use the demonstrations in this documentation.
In[1072]:=
Needs["QuantumMob`Q3`"]
Definition
A supermap is a linear mapping of linear operators. Recall that linear operators on a vector space themselves form a vector space
ℒ()
. A considerable amount of interest in supermaps came with the booming of quantum information theory in the 1990s when it became clear that supermaps are important in the study of entanglement. Since then, mathematical theories on supermaps have been developed at a notably fast pace and applied to a wide range of subjects in quantum computation and quantum information.
involving a single unitary operator. Naturally it describes the unitary dynamics of a closed system. Note that if there are more than one Kraus elements involved in the Kraus representation, the associated dynamics is not unitary even if the Kraus elements are all unitary. For example, consider a supermap
ℱ(ρ)=(1-p)ρ+
p
3
(XρX+YρY+ZρZ)
.
Here, all Kraus elements (
I,X,Y,Z
up to normalization factors) are unitary. However,
is the matrix of the expansion coefficients. Putting it back to the Kraus representation, we have
ℱ(ρ)=
Σ
μν
E
μ
ρ
†
E
ν
(M
†
M
)
μν
.
The square matrix
M
†
M
of size
(dim)×(dim)
is positive semidefinite, and can be decomposed into
M
†
M
=VΛ
†
V
, where
V
is a unitary matrix and Λ is a diagonal matrix with all elements non-negative. We now define new Kraus elements
′
F
ν
:=
Σ
μ
(V
Λ
)
μν
.
Then, it is clear that they are mutually orthogonal. Furthermore,
ℱ(ρ)=
Σ
μ
′
F
ν
ρ
′†
F
μ
,
where
N≤(dim)×(dim)
. Therefore, it is noted that a set of mutually orthogonal Kraus elements is optimal in the sense that it has no more elements than
(dim)×(dim)
.
Unitary Freedom of Kraus Elements
(to be completed)
Choi Isomorphism
The Choi isomorphism is a one-to-one correspondence between supermaps and (usual) operators. It allows us to inspect a supermap in terms of the corresponding operator. It can reveal additional properties of the supermap that is not immediately clear from the supermap itself.
Unitary Representation
A quantum operation can always be regarded as a unitary operator on an extended system, which involves an “environment” in addition to the original “system”. Although it is not particularly useful for practical applications, the unitary representation provides a clear physical insight into the underlying physical processes described by the quantum operation.
Measurements as Quantum Operations
Examples
So far, we discussed a general description of quantum noisy processes in terms of the corresponding quantum operations. Let us now present some examples for a single qubit. We consider some limiting cases that allow us to easily grasp the physical meaning of the Kraus elements.
Phase Damping
The phase damping process is a decoherence process without involving any relaxation of energy or change in the population over the states. In this sense, it may be regarded as a pure decoherence process without involving any energy relaxation. For this reason, it is also referred to as a dephasing process.
The unitary representation of a phase damping process in a single-qubit system is given by
or, more explicitly,
Due to the entanglement, the final state of the system alone cannot be a pure state and coherence in the initial state has been lost through the process. With the prescription above, the Kraus elements are given by
and the corresponding quantum operation is given by
Note that they are not orthogonal,
It is therefore more convenient and efficient to choose mutually orthogonal Kraus elements
In short, the quantum operation for the phase damping process is written as
It is also convenient to expand the density operator into
The coherence has disappeared completely.
The phase damping process is specified by two Kraus elements.
Here, p is the probability for the phase to be flipped.
The supermap transforms the above density operator as follows.
Amplitude Damping
In the unitary representation, the process is described by the overall unitary operator such that
They are already orthogonal to each other. The quantum operation describing the amplitude damping process reads as
With the expansion of the density operator ρ,
the above operation reads as
The amplitude damping process is specified by two Kraus elements.
Here, p is the probability for the phase to be flipped.
The supermap transforms the above density operator as follows.
Depolarizing
From the above unitary representation, we can get the Kraus elements
The Kraus elements are already orthogonal to each other. One can also check that they satisfy the completeness relation
as they should. In the Kraus representation with the above Kraus elements, a density operator ρˆ is transformed under the decoherence process as
In terms of the components in the expansion
This implies that under the process, the “spin” polarization (i.e., the Bloch vector corresponding to the resulting density operator)
The depolarizing process is specified by three Kraus elements.
Here, p is the probability for the phase to be flipped.
The supermap transforms the above density operator as follows.