Bloch–Gell–Mann Representation in Qudit Systems

Updated 28 December 2025

The Bloch–Gell–Mann representation is a formalism that expands density matrices using traceless Hermitian generators to provide a complete geometric and algebraic characterization of d-level quantum states.
It employs generalized Gell–Mann matrices to decompose quantum states, facilitating precise state reconstruction, quantum dynamics analysis, and applications in quantum machine learning.
This unified approach integrates algebraic structure, positivity constraints, and dynamical evolution across dimensions, extending from qubits to multi-qudit and higher spin systems.

The Bloch–Gell–Mann representation is a generalization of the well-known Bloch sphere formalism for qubits that provides a complete geometric and algebraic characterization of finite-dimensional quantum states (qudits). By expanding density matrices and observables in the orthonormal basis of traceless Hermitian generators of $\mathrm{SU}(d)$ —the generalized Gell–Mann matrices—any $d$ -level system can be parametrized by its so-called Bloch–Gell–Mann vector. This approach enables a unified description of state space geometry, operator expansions, dynamical evolution, and operational applications, and has broad utility in quantum information science, quantum foundations, and quantum machine learning (Valtinos et al., 2023, Loubenets, 2019, Loubenets et al., 2020, Chew et al., 2020, Giraud et al., 2014).

1. Definition of Generalized Gell–Mann Matrices

Let $d$ be the Hilbert space dimension. The $d^2-1$ traceless Hermitian generators $\{\lambda_i\}$ of $\mathfrak{su}(d)$ are constructed as follows (Loubenets, 2019, Loubenets et al., 2020):

Symmetric off-diagonal generators: $\lambda^{(s)}_{m k} = |m\rangle\langle k| + |k\rangle\langle m|$ for $1\leq m<k\leq d$ ;
Antisymmetric off-diagonal generators: $\lambda^{(a)}_{m k} = -i|m\rangle\langle k| + i|k\rangle\langle m|$ ;
Diagonal generators: For $l=1,\dots,d-1$ ,

$d$ 0

These matrices are orthonormal under the Hilbert–Schmidt inner product: $d$ 1 The case $d$ 2 recovers the usual Pauli matrices; $d$ 3 yields the Gell–Mann matrices, explicitly: $d$ 4 (Valtinos et al., 2023, Loubenets et al., 2020).

2. Expansion of Quantum States: Bloch–Gell–Mann Parametrization

Any Hermitian, trace-one operator (density matrix) $d$ 5 on $d$ 6 admits the expansion (Loubenets, 2019, Loubenets et al., 2020): $d$ 7 where $d$ 8 are real and form the $d$ 9-dimensional Bloch–Gell–Mann vector. The normalization ensures $d$ 0 and Hermiticity follows from that of the $d$ 1. The collection $d$ 2 forms an orthonormal basis for the Hilbert–Schmidt space of Hermitian matrices (Loubenets, 2019).

For $d$ 3 (qutrits), for instance,

$d$ 4

(Valtinos et al., 2023).

3. Structure Constants and Algebraic Properties

The Gell–Mann matrices close under commutation and anticommutation, with structure constants $d$ 5 (totally antisymmetric) and $d$ 6 (totally symmetric) (Loubenets, 2019, Loubenets et al., 2020): $d$ 7 where

$d$ 8

These constants enter both dynamical equations and the closure relations for the product of generators, central to the structure of $d$ 9 (Loubenets et al., 2020).

4. Geometry and Constraints of the Generalized Bloch Body

The set of physically valid Bloch–Gell–Mann vectors forms a convex subset of $d^2-1$ 0, defined by positivity constraints, Hermiticity, and normalization (Valtinos et al., 2023, Loubenets, 2019):

Purity bound: $d^2-1$ 1.
Positivity: All eigenvalues of $d^2-1$ 2 are $d^2-1$ 3. For $d^2-1$ 4, the positivity region is the unit ball; for $d^2-1$ 5, it is a proper convex subset of the corresponding ball. For example, for qutrits ( $d^2-1$ 6) the “Bloch body” is inside a ball of radius $d^2-1$ 7 but its boundary is a high-dimensional, algebraically complicated manifold defined by higher-order invariants such as $d^2-1$ 8 (Valtinos et al., 2023, Loubenets, 2019).

For three-qubit (su(8)) systems, the physical region is a convex compact subset of the $d^2-1$ 9 ball in $\{\lambda_i\}$ 0 (Chew et al., 2020). In general, boundary points corresponding to pure states form a complex projective manifold $\{\lambda_i\}$ 1.

5. Extension to Tensors and Higher Spins

The Bloch–Gell–Mann formalism admits generalization to arbitrary spin- $\{\lambda_i\}$ 2 systems using tensor expansions (Weinberg matrices) (Giraud et al., 2014). For spin- $\{\lambda_i\}$ 3 ( $\{\lambda_i\}$ 4), any density matrix $\{\lambda_i\}$ 5 may be expanded as

$\{\lambda_i\}$ 6

where the $\{\lambda_i\}$ 7 form a symmetric, Hermitian tensor basis and the coefficients are real. This construction recovers the standard Gell–Mann generators as a subset for $\{\lambda_i\}$ 8, while for $\{\lambda_i\}$ 9 it reduces to the Pauli matrices. This tensor formalism naturally encodes transformation properties under $\mathfrak{su}(d)$ 0 and partial trace operations, and is fundamental for the characterization of phenomena such as quantum polarization and anticoherence (Giraud et al., 2014).

6. Applications: Quantum Information and Quantum Machine Learning

The Bloch–Gell–Mann representation is essential in several areas:

Quantum tomography: Extraction of Bloch vector components allows full state reconstruction from expectation values of the generalized Gell–Mann operators (Loubenets, 2019, Chew et al., 2020).
Quantum dynamics: The von Neumann equation for a general qudit reduces to a differential equation for the Bloch vector, governed by the structure constants $\mathfrak{su}(d)$ 1:

$\mathfrak{su}(d)$ 2

for a Hamiltonian $\mathfrak{su}(d)$ 3 (Loubenets et al., 2020).

Quantum machine learning: The Gell–Mann feature map encodes $\mathfrak{su}(d)$ 4-dimensional data directly into a qutrit using the exponential map $\mathfrak{su}(d)$ 5. This yields richer high-dimensional embeddings and more expressive circuit ansätze for classification, support vector machines, and variational quantum algorithms compared to qubit-based feature maps (Valtinos et al., 2023).

7. Specific Implementations and Coordinate Mapping

For composite and higher-dimensional systems, the generalized Gell–Mann matrices provide a bridge to the Pauli-tensor basis, particularly for multi-qubit systems ( $\mathfrak{su}(d)$ 6). For example, $\mathfrak{su}(d)$ 7 generators can be explicitly expanded as linear combinations of three-qubit Pauli operators, with the mapping coefficients determined by change-of-basis matrices (Chew et al., 2020). Measuring all $\mathfrak{su}(d)$ 8 Pauli correlators enables reconstruction of the Bloch–Gell–Mann vector for states of $\mathfrak{su}(d)$ 9 qubits. This basis change is instrumental in experimental tomography, Hamiltonian learning, and quantum simulation.