---
title: Kraus Decomposition for Quantum Channels
url: https://www.emergentmind.com/topics/kraus-decomposition
type: topic
---

# Kraus Decomposition for Quantum Channels

The Kraus decomposition, also known as the operator-sum representation, is a fundamental structural result in quantum information theory and the theory of open quantum systems. It provides a canonical form for all linear maps that are completely positive and trace-preserving (CPTP), describing the most general evolution—both noisy and unitary—that a quantum state may undergo due to its interaction with an environment. This decomposition reveals not only the structure of physical quantum channels but also underpins a wide array of techniques across quantum physics, information processing, statistical learning, and emerging applications in geometry and representation theory.

## 1. Mathematical Definition and Derivation

Every CPTP map $\mathcal{E}$ between operator spaces $\Omega_n$ on a finite-dimensional Hilbert space can be written in the Kraus (operator-sum) form:
\[
\mathcal{E}(\rho) = \sum_{i=1}^r K_i \, \rho \, K_i^\dagger , \qquad \sum_{i=1}^r K_i^\dagger K_i = I .
\]
Here $\rho$ is a density operator, $K_i$ are "Kraus operators" ($K_i \in \mathbb{C}^{n \times n}$), and $r \leq n^2$ is the Kraus rank, minimized to match the rank of the corresponding Choi matrix [1011.1585, 1510.09081]. 

Starting from a joint system-environment unitary evolution $U$ acting on $\mathcal{H}_S \otimes \mathcal{H}_E$ and an initial uncorrelated state $\rho_S(0) \otimes \rho_E$, the reduced system state at time $t$ is:
\[
\rho_S(t) = \mathrm{Tr}_E\{ U(t) [\rho_S(0) \otimes \rho_E] U^\dagger(t) \} = \sum_k K_k \rho_S(0) K_k^\dagger .
\]
Kraus operators $K_k$ are constructed from matrix elements of $U$—for instance, $ \langle s_k|K_l|s_m\rangle = \langle s_k, e_l| U | s_m, e_0 \rangle $ for suitable system and environment bases [1510.09081, 1011.1585].

## 2. Fundamental Properties and Non-Uniqueness

The Kraus decomposition guarantees three essential properties for quantum channels:

- **Linearity:** Each $K_i$ is linear; the map $\mathcal{E}$ is linear as a result.
- **Complete Positivity:** For any extended Hilbert space and any positive semidefinite input, $(I \otimes \mathcal{E})$ maps positive semidefinite operators to positive semidefinite operators.
- **Trace Preservation:** $\sum_i K_i^\dagger K_i = I$ ensures $\operatorname{Tr} \mathcal{E}(\rho) = \operatorname{Tr} \rho$ for any $\rho$.

The Kraus representation is non-unique: alternative sets $\{K'_l\}$ connected by unitary rotations ($K'_l = \sum_n V_{ln} K_n$ for unitary $V$) describe the same channel. The minimal Kraus rank equals the rank of the Choi matrix of the map [1011.1585, 1510.09081].

For any CPTP map on dimension $d_S$, the minimal number of required Kraus operators is bounded above by $d_S^2$ [1510.09081].

## 3. Explicit Construction and Examples

Canonical construction of the Kraus decomposition utilizes the Choi-Jamiolkowski isomorphism. Given a channel $\mathcal{E}$, the Choi matrix $D_{\mathcal{E}}$ is obtained by reshaping the superoperator or by applying $\mathcal{E} \otimes \mathrm{id}$ to a maximally entangled state. Its spectral decomposition $D_{\mathcal{E}} = \sum_i \sigma_i |u_i\rangle\langle u_i|$ yields Kraus operators via unvectorization:
\[
K_i = \sqrt{\sigma_i} \, \operatorname{unvec}(|u_i\rangle) .
\]
This construction allows for systematic calculation and composition of channels [1011.1585].

### Prototype: Amplitude-Damping Channel

A two-level atom (qubit) coupled to the vacuum electromagnetic field realizes an amplitude-damping channel with Kraus operators
\[
K_0 = \begin{pmatrix} 1 & 0 \\ 0 & \sqrt{1-p} \end{pmatrix}, \quad K_1 = \begin{pmatrix} 0 & \sqrt{p} \\ 0 & 0 \end{pmatrix}, \quad p=1-e^{-\gamma t}.
\]
These satisfy $\sum_{l=0}^1 K_l^\dagger K_l = I$ [1510.09081].

### Microscopic Derivation: Generalized Depolarizing Channel

From a full open-system Hamiltonian and microscopic Lindblad equation, physically motivated generalized Kraus operators emerge, capturing both anisotropies and Lamb-shift-induced coherent rotations, refining the standard depolarizing-channel model [1512.07843].

### Non-Markovian Kraus Expansions

In non-Markovian settings, especially with reservoir-induced memory effects, the Kraus map becomes infinite, indexed by interaction order, and can be constructed via time-ordered expansions and continued-fraction solutions in Laplace space [1807.07325].

## 4. Constraints, Reduction, and Resource Theories

The structure and number of Kraus operators are crucial in specialized resource theories. For incoherent operations (IO) and strictly incoherent operations (SIO), further criteria are imposed: each $K_n$ maps incoherent states to incoherent states (in IO), and both $K_n$ and $K_n^\dagger$ preserve incoherence (in SIO). Recent results refine the maximal number of required incoherent Kraus operators: four suffices for qubits (improving upon the prior five), and 32 for single qutrit IO (from 39); SIO for qutrits can be realized with 13 Kraus operators (from 15). These reductions simplify state-transformation characterizations and have practical implications for quantum coherence manipulation [2005.01083].

| Scenario                | Previous Bound | Improved Bound | Reference      |
|-------------------------|---------------|---------------|---------------|
| Qubit IO                | 5             | 4             | [2005.01083]  |
| Qutrit IO               | 39            | 32            | [2005.01083]  |
| Qutrit SIO              | 15            | 13            | [2005.01083]  |

## 5. Computational, Learning, and Tomographic Aspects

The Kraus decomposition underpins both theoretical and operational approaches to quantum process identification:

- **Quantum Process Tomography (QPT):** Rather than reconstructing a full $N^2 \times N^2$ Choi matrix, low-rank Kraus ansatzes directly optimize a physically admissible set of $k \ll N^2$ operators, employing Riemannian gradient descent on the Stiefel manifold. This approach guarantees CPTP structure at every iteration, scales to larger systems than standard semidefinite programming, and handles compressed (informationally incomplete) measurement sets accurately [2208.00812].

- **Low-Rank Learning:** In regression or matrix completion contexts, models that respect Kraus structure—expressing a map $\Phi(X) = \sum_{j=1}^r A_j X A_j^\dagger$—can be learned efficiently via SGD or block coordinate descent. Generalization is characterized via a pseudo-dimension bound depending on the Kraus rank $r$, and stacking layers maintains the complete-positivity and structure of the learned map [2007.00935].

| Application            | Kraus Role                  | Algorithmic Features                 | Reference      |
|------------------------|-----------------------------|--------------------------------------|---------------|
| Quantum process ID     | Parameterization, exact CP  | Gradient descent on Stiefel manifold | [2208.00812]  |
| Matrix regression      | Low-rank CP map             | SGD, backprop, generalization bound  | [2007.00935]  |

## 6. Generalizations and Emerging Applications

Recent trends extend the formalism beyond standard quantum information:

- **Geometric and Representation-Theoretic Extensions:** The w-Kraus decomposition generalizes the Kraus form to arbitrary quadratic geometries (elliptic, Lorentzian, etc.), with completeness constraints reflecting the appropriate inner product and positivity notion [2605.10317].

- **Knowledge Embedding Models:** Relations in (e.g.) knowledge-graph embedding are modeled as CPTP maps whose structure is dictated by three axioms: linearity, trace preservation, and complete positivity. The Kraus decomposition is then inevitable, and Kraus rank provides a natural per-relation complexity measure. Closure under composition (multi-hop reasoning) and handling of high-fan-out relations are guaranteed by the CPTP structure [2605.10317].

| Generalization         | Key Principle               | Structure/Constraint                             | Reference     |
|-----------------------|-----------------------------|--------------------------------------------------|--------------|
| w-Kraus (geometry)    | CPTP in arbitrary geometry  | $\sum_i K_i^\dagger W K_i = W$                   | [2605.10317] |
| Knowledge graph       | CPTP channel for relations  | Completeness $\sum_i K_i^\dagger K_i = I$        | [2605.10317] |

## 7. Computational Methods and Automation

Canonical methods for obtaining Kraus decompositions from superoperators or Choi matrices are well developed. Standard practice involves vectorization (column-stacking), reshuffling (mapping between process and Choi representations), and singular value (Schmidt) decomposition. Automation of these procedures is available via platforms such as Mathematica, allowing symbolic and numerical calculations for arbitrary channels and composite systems [1011.1585]. The superoperator formalism facilitates channel composition via mapping composition, and Kraus composition is simply $K_i L_j$ for composed channels.

## References

- "The Kraus representation for the dynamics of open quantum systems" [1510.09081]
- "Singular value decomposition and matrix reorderings in quantum information theory" [1011.1585]
- "Generalized Kraus operators for the one-qubit depolarizing quantum channel" [1512.07843]
- "The reduction of the number of incoherent Kraus operations for qutrit systems" [2005.01083]
- "Partial Trace Regression and Low-Rank Kraus Decomposition" [2007.00935]
- "Gradient-descent quantum process tomography by learning Kraus operators" [2208.00812]
- "Continued-fraction representation of the Kraus map for non-Markovian reservoir damping" [1807.07325]
- "Kraus map closed-form solution for general master equation dynamics" [2603.11207]
- "Relations Are Channels: Knowledge Graph Embedding via Kraus Decompositions" [2605.10317]

Source: https://www.emergentmind.com/topics/kraus-decomposition