---
title: 'Mutually Unbiased Bases: Theory & Applications'
url: https://www.emergentmind.com/topics/mutually-unbiased-bases-mubs
type: topic
---

# Mutually Unbiased Bases: Theory & Applications

Mutually unbiased bases (MUBs) are families of orthonormal bases in a finite-dimensional Hilbert space that realize maximal complementarity: if a system is prepared in a vector from one basis, then measurement in any different basis from the family yields a uniform outcome distribution. Formally, two orthonormal bases \(B=\{\lvert e_i\rangle\}_{i=1}^d\) and \(C=\{\lvert f_j\rangle\}_{j=1}^d\) are mutually unbiased when
\[
|\langle e_i,f_j\rangle|^2=\frac1d
\quad\text{for all }i,j.
\]
A set of MUBs is complete if it contains \(d+1\) bases, which is the maximal possible number in dimension \(d\) [1309.6557][2508.15412]. MUBs sit at the intersection of finite-field methods, generalized Pauli groups, complex Hadamard matrices, projective and Hjelmslev geometries, operator algebras, and several operational tasks in quantum information.

## 1. Definition, maximality, and the existence problem

The defining condition of mutual unbiasedness is equivalent to uniform transition probabilities between distinct bases, or, in projective language, to the rank-one projector relation \(\operatorname{Tr}(P_iQ_j)=1/d\) for all \(i,j\) [2508.15412]. A classical bound shows that no set of MUBs in \(\mathbb{C}^d\) can contain more than \(d+1\) bases, and a set achieving that bound is called complete [2508.15412].

Complete sets are known to exist whenever \(d\) is a prime power, \(d=p^n\). Standard constructions proceed either through finite fields, in the style of Wootters–Fields, or through commuting classes of the generalized Pauli group, as in Bandyopadhyay, Boykin, Roychowdhury, and Vatan [1309.6557]. In prime dimension \(p\), the complete set of \(p+1\) stabilizer bases is essentially unique up to unitary equivalence, whereas in prime-power dimension \(q=p^n\) many inequivalent complete stabilizer MUBs exist [1503.00003].

Outside prime-power dimensions, the existence problem remains central and unresolved in general. The smallest non-prime-power case, \(d=6\), is the canonical example: only three MUBs are known, and strong evidence suggests that no more than three simultaneous MUBs exist [1309.6557][1304.4081]. A stronger nonexistence statement for completeness was established by relating MUBs to mutually orthogonal Latin squares: if a complete set of \(N+1\) MUBs exists in \(\mathbb{C}^N\), then a complete set of \(N-1\) MOLS of order \(N\) must exist; since no complete set of MOLS exists for \(N=6\), a complete set of seven MUBs cannot exist in dimension six [2511.03537]. The same implication rules out complete MUB sets in dimension ten as well [2511.03537].

## 2. Algebraic constructions in prime-power dimension

A standard algebraic route begins from generalized Pauli operators. On \(\mathbb{C}^p\), one defines
\[
X\lvert k\rangle=\lvert k+1 \bmod p\rangle,\qquad
Z\lvert k\rangle=\omega^k\lvert k\rangle,\qquad
\omega=e^{2\pi i/p}.
\]
On \(n\) qupits, tensor products of \(X^lZ^m\) generate the generalized Pauli group, and commuting classes of size \(p^n-1\) yield MUBs through their common eigenbases [1309.6557]. This stabilizer perspective underlies most known explicit complete constructions and also the symmetry analysis of stabilizer MUBs [1503.00003].

A compact reformulation is given by the graph-state formalism. For \(n\) \(p\)-level systems, with \(p\) prime, a generalized graph is specified by a symmetric adjacency matrix \(A\in\mathbb{Z}_p^{n\times n}\), with self-loops allowed. The associated graph state is
\[
\lvert G\rangle=\prod_{i\le j}U_{i,j}^{A_{i,j}}\lvert+\rangle^{\otimes n},
\]
where \(U_{i,i}\) is a one-body phase gate and \(U_{i,j}\) for \(i\neq j\) is a controlled-phase gate [1309.6557]. The corresponding graph-state basis is generated by local \(Z\)-shifts. The decisive criterion is that if \(A_r\) and \(A_s\) are symmetric adjacency matrices, then the associated bases are mutually unbiased whenever
\[
\det(A_r-A_s)\neq 0\quad \text{mod }p.
\]
This converts mutual unbiasedness into invertibility of matrix differences over \(\mathbb{Z}_p\) [1309.6557].

The same paper shows how a single symmetric matrix \(Q\) with irreducible characteristic polynomial over \(\mathbb{Z}_p\) generates a matrix representation of the finite field \(\mathbb{F}_{p^n}\):
\[
S=\left\{\sum_{i=0}^{n-1}a_iQ^i:\vec a\in\mathbb{Z}_p^n\right\}.
\]
Because \(S\) is a field, any two distinct matrices in \(S\) differ by an invertible matrix, so the corresponding graph-state bases are pairwise mutually unbiased. Together with the computational basis, they form a complete set of \(p^n+1\) MUBs [1309.6557]. In many cases \(Q\) may be chosen tridiagonal with unit super- and subdiagonals and diagonal entries specified by a single vector \(\vec d\in\mathbb{Z}_p^n\), so one graph, or even one vector, suffices to encode the full complete set [1309.6557].

This algebraic economy is closely related to an older unifying viewpoint in terms of modules and finite geometries. Several standard complete constructions—planar-function, Alltop, symplectic-spread, and Galois-ring constructions—admit a description in which the non-standard MUB vectors form a module under componentwise multiplication, while the corresponding exponent patterns realize subspaces of projective geometries \(PG(\cdot)\) or projective Hjelmslev geometries \(PHG(\cdot)\) [1206.0531]. This suggests that the known prime-power constructions share a common higher-dimensional geometric substrate rather than merely a counting analogy.

## 3. Geometric, operator-theoretic, and classification frameworks

Beyond explicit construction, MUBs admit several complementary structural descriptions. One is geometric: the space of unordered projective orthonormal bases of \(\mathbb{C}^n\) can be modeled as
\[
\mathcal{M}_n = U(n)/\mathrm{C}_n,
\]
where \(\mathrm{C}_n=S_n\ltimes T^n\) is the monomial unitary group [2508.15412]. A \(U(n)\)-equivariant embedding into the Grassmannian of \((n-1)\)-planes in \(\mathfrak{su}(n)\) induces a metric
\[
D(p,q)=\sqrt{\,n-1-\sum_{i,j}\Bigl(\operatorname{Tr}(P_iQ_j)-\frac1n\Bigr)^2\,},
\]
with \(D(p,q)=\sqrt{n-1}\) if and only if the two bases are mutually unbiased [2508.15412]. In this framework, the familiar relation between a basis unbiased to the standard basis and a complex Hadamard matrix becomes a statement about the orbit decomposition of the subspace \(\mathcal V(e)\) under the left action of \(\mathrm C_n\) [2508.15412].

The same geometric language organizes extension problems. Given a list of \(k\) MUBs, the candidate extensions form the intersection
\[
\mathcal N\bigl(q^{(1)},\dots,q^{(k)}\bigr)
 = \bigcap_{m=1}^k \mathcal V\bigl(q^{(m)}\bigr),
\]
and inequivalent extensions are precisely the distinct orbits of this space under the simultaneous stabilizer of the given list [2508.15412]. In dimension \(4\), this orbit analysis yields new symmetries that reduce the parameter space of MUB triples by a factor of \(4\) [2508.15412].

A second structural viewpoint replaces bases by commuting classes of matrices. There is a one-to-one correspondence between MUBs in \(\mathbb{C}^n\) and maximal commuting classes of orthogonal normal matrices in \(M_n(\mathbb{C})\): from \(m\) MUBs one obtains \(m\) commuting classes, each consisting of \(n\) commuting orthogonal normal matrices, and conversely a maximal commuting basis of \(M_n(\mathbb{C})\) ensures a complete set of MUBs [2507.02492]. The extendability problem then becomes equivalent to finding real points on a quadratic affine variety encoding orthogonality and unbiasedness constraints, with the unbiasedness-to-standard-basis sector forming a complete intersection prime ideal [2507.02492].

A third viewpoint concerns symmetry. Complete stabilizer MUBs are often group covariant, but sharp covariance is exceptionally rare. Among stabilizer MUBs, only the cases \(d=2\) and \(d=4\) are sharply covariant, and this remains true even if antiunitary symmetries are allowed [1503.00003]. This sharply distinguishes MUBs from SICs, for which all known examples are sharply covariant [1503.00003].

## 4. Operational roles in quantum information

MUBs are central to tomography, uncertainty, communication, and certification. In state tomography, a complete set of \(d+1\) MUBs is informationally complete and minimal among orthonormal-basis measurement schemes. In the star-product formulation, if \(p_{a\alpha}=\operatorname{Tr}(\rho\,\Pi_{a\alpha})\) are the probabilities of the MUB projectors \(\Pi_{a\alpha}\), then
\[
\rho=\sum_{a=0}^{d}\sum_{\alpha=0}^{d-1} p_{a\alpha}\,\Pi_{a\alpha}-I,
\]
so knowledge of all MUB probabilities suffices for exact reconstruction [1008.2675]. For \(n\)-qubit systems, a complete set consists of \(2^n+1\) bases, compared with \(4^n\) Pauli observables, and an explicit Clifford-circuit synthesis with the three-stage form
\[
U(j)=U_{CZ}(j)\,U_S(j)\,H^{\otimes n}
\]
constructs every MUB circuit in \(O(n^3)\) time [2311.11698].

The same \(n\)-qubit MUB circuits can be used as the ensemble in classical shadow tomography. Uniform MUB sampling yields an informationally complete channel
\[
\mathcal M_u(\rho)=\frac{\rho+\operatorname{tr}(\rho)I}{2^n+1},
\qquad
\mathcal M_u^{-1}(X)=(2^n+1)X-\operatorname{tr}(X)I,
\]
but for general observables the variance can be exponential in \(n\) [2310.09644]. For approximately MUBs-average observables or states, and for MUBs-sparse instances under biased sampling, the variance becomes polynomial in \(n\) [2310.09644].

In prepare-and-measure self-testing, the \(2^d\to1\) QRAC has optimal average success probability
\[
\bar p_Q(d)=\frac12\Bigl(1+\frac1{\sqrt d}\Bigr),
\]
and this optimum is achieved if and only if Bob’s two measurements are rank-1 projective and mutually unbiased [1803.00363]. This gives a device-independent-style certification, within a dimension-bounded setting, of a pair of MUBs from operational performance alone [1803.00363].

MUBs also support entanglement detection. For bipartite \(d\times d\) systems, if \(I_m\) denotes the sum of mutual predictabilities across \(m\) local MUB settings, then all separable states satisfy
\[
I_m\le 1+\frac{m-1}{d},
\]
and in particular \(I_{d+1}\le 2\) for a complete set [1309.2271]. In the magic simplex of Bell-diagonal states, this criterion detects PPT-bound entanglement when a complete set of MUBs is used; for \(d=3,\dots,9\) the paper reports explicit PPT states with \(I_{d+1}>2\) [1309.2271].

At the same time, operational behavior can depend on which MUB set and which outcome labeling are used. Under coherent measurement choice, certain relabelings of Wootters–Fields MUBs can make the measurement outcome perfectly guessable, with \(P_g=1\) and \(H_{\min}(X|E)=0\), even though the measurements remain mutually unbiased in the usual overlap sense [2006.08226]. This does not alter the mathematical definition of MUBs, but it shows that some security statements depend on implementation details beyond pairwise overlaps.

## 5. Approximate and generalized variants

Because exact complete sets are unavailable in many dimensions, several approximate notions have been developed. A \(\beta\)-AMUB is a family of orthonormal bases for which all cross-basis overlaps satisfy \(|\langle\psi_i^l,\psi_j^m\rangle|\le \beta/\sqrt d\), while an Almost Perfect MUB (APMUB) imposes the sharper biangular spectrum
\[
\Delta=\left\{0,\frac{\beta}{\sqrt d}\right\},
\qquad
\beta=1+O(d^{-\lambda}),\ \lambda>0,
\]
with sparse basis vectors whose nonzero components all have equal magnitude [2402.03964][2512.12828]. For composite dimensions of the form \(d=(q-e)(q+f)\), with \(q\) a prime power and \(0\le f\le e\), constructions based on resolvable block designs yield \(O(\sqrt d)\) APMUBs, often far more than the known exact MUBs in those dimensions [2402.03964].

The quality of approximate families can be quantified by several measures. These include the \(t\)-coherence
\[
\Omega_t^{(l,m)}=\sum_{i,j}|\langle\psi_i^l,\psi_j^m\rangle|^{2t},
\]
the max and RMS deviations \(\tau\) and \(\sigma\) from \(1/\sqrt d\), Bengtsson’s distance \(D^2\), and the projective 2-design defect through the normalized frame potential \(\Phi_2\) [2512.12828]. For exact MUB pairs, \(\Omega_2=1\), \(D^2=1\), and \(\tau=\sigma=0\); for APMUBs with \(\Delta=\{0,\beta/\sqrt d\}\), these quantities can be computed solely from \(d\) and \(\beta\), without using the internal construction details [2512.12828].

A different generalization relaxes bases to measurements. Mutually unbiased measurements (MUMs) are projective measurements satisfying
\[
P_a=d\,P_aQ_bP_a,\qquad
Q_b=d\,Q_bP_aQ_b
\]
for all outcomes \(a,b\) [2204.11886]. MUMs are strictly more general than MUBs because the number of outcomes need not equal the Hilbert-space dimension, and there exist MUMs that are not direct sums of MUBs [2204.11886]. This broader complementarity notion supports superdense coding protocols and admits constructions from quaternionic Hadamard matrices [2204.11886].

## 6. Special dimensions, implementations, and open directions

Dimension \(6\) remains the most studied exceptional case. Experimentally, three MUBs have been implemented and tested for single-photon six-dimensional states, encoded either in a hybrid polarization–orbital-angular-momentum Hilbert space or in pure orbital angular momentum. The measured \(18\times18\) overlap matrices showed close agreement with the ideal MUB block structure, with reported similarities \(S=(99.19\pm0.04)\%\) and \(S=(99.06\pm0.04)\%\) in the two encodings [1304.4081]. These experiments confirm the practical availability of MUB triples in \(d=6\), but do not alter the theoretical evidence that no fourth basis exists [1304.4081].

The nonexistence question in \(d=6\) has also been reformulated numerically through Bell inequalities. For every pair \((n,d)\), one can construct a Bell inequality whose maximal quantum violation in dimension \(d\) is attainable if and only if \(n\) MUBs exist in that dimension [2203.09429]. Three numerical approaches—see-saw optimization, nonlinear semidefinite programming, and Monte Carlo—correctly recover the known prime-power cases and all suggest that four MUBs do not exist in dimension six; the numerical optimizers coincide with the “four most distant bases” found earlier by optimizing a distance measure [2203.09429]. The same Monte Carlo evidence suggests that at most three MUBs exist in dimension ten [2203.09429].

There are also infinite-dimensional analogues, but their behavior depends strongly on the underlying field. For states over the real line \(L^2(\mathbb R)\), only three MUBs are known, whereas over the \(p\)-adic numbers \(L^2(\mathbb Q_p)\), for every prime \(p>2\) there are at least \(p+1\) mutually unbiased generalized bases in the rigged Hilbert-space sense [1109.0060]. The paper attributes this contrast to the different harmonic analysis of \(\mathbb R\) and \(\mathbb Q_p\), especially the ultrametric structure and finite residue field in the \(p\)-adic case [1109.0060].

Several open directions remain active. Non-prime-power dimensions continue to motivate both nonexistence criteria and alternative constructions [2511.03537][2203.09429]. The classification of inequivalent complete sets, particularly under local Clifford or graph operations, remains open even in prime-power dimensions [1309.6557]. This suggests that MUBs are best understood not as a single construction problem, but as a network of equivalent extremal problems in finite geometry, operator theory, algebraic geometry, and quantum information.

Source: https://www.emergentmind.com/topics/mutually-unbiased-bases-mubs