Kirkwood-Dirac Pseudo-Distribution

• The Kirkwood–Dirac pseudo-distribution is a quasiprobability representation that extends joint probability to incompatible quantum measurements while allowing negative or complex values.
• It maintains proper Born rule marginals and quantifies quantum resources like nonclassicality, contextuality, and coherence through its distinctive negativity.
• Operational protocols such as weak measurement schemes and moment-based witnesses enable its experimental verification and application in quantum metrology and thermodynamics.

The Kirkwood–Dirac (KD) pseudo-distribution is a canonical quasiprobability representation of quantum states relative to two orthonormal bases. Originally introduced to extend the notion of joint probability to incompatible measurements in quantum mechanics, the KD distribution is now central to quantum foundations, metrology, quantum information theory, and the study of nonclassicality. A distinctive property of the KD distribution is that its entries can be negative or nonreal, in sharp contrast to genuine probability distributions, and such negativity or complex character serves as a resource indicator for quantum advantage, contextuality, and coherence.

1. Formal Definition and Fundamental Properties

Given a $d$-dimensional Hilbert space and two orthonormal bases $A = \{|a_j\rangle\}_{j=1}^d$ and $B = \{|b_k\rangle\}_{k=1}^d$, the KD pseudo-distribution of a density operator $\rho$ is defined as

$$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$

For a pure state $\rho=|\psi\rangle\langle\psi|$, this specializes to

$$Q_{jk}(|\psi\rangle) = \langle a_j|\psi\rangle\,\langle\psi|b_k\rangle\,\langle b_k|a_j\rangle.$$

The KD matrix obeys the following properties:

• Normalization: $\sum_{j,k} Q_{jk} = \operatorname{Tr}\rho = 1$;
• Correct marginals (Born probabilities): $\sum_k Q_{jk} = |\langle a_j|\psi\rangle|^2$ and $\sum_j Q_{jk} = |\langle\psi|b_k\rangle|^2$;
• Quasiprobability: KD entries can be negative or complex; only in special cases does the KD array correspond to a classical probability distribution (every entry real and nonnegative) (Xu, 2022).

For general observables $A = \{|a_j\rangle\}_{j=1}^d$0 and $A = \{|a_j\rangle\}_{j=1}^d$1 with spectral decompositions into nondegenerate projectors, the definition generalizes to

$A = \{|a_j\rangle\}_{j=1}^d$2

(Arvidsson-Shukur et al., 2024).

This family further admits generalization to sequences of noncommuting observables or even POVMs, and is operationally equivalent to weak-value statistics for appropriately defined measurement sequences.

2. Classicality, Nonclassicality, and Support Uncertainty

A pure state $A = \{|a_j\rangle\}_{j=1}^d$3 is termed KD-classical (w.r.t. bases $A = \{|a_j\rangle\}_{j=1}^d$4) if all KD matrix elements are real and nonnegative, $A = \{|a_j\rangle\}_{j=1}^d$5 for all $A = \{|a_j\rangle\}_{j=1}^d$6; otherwise it is KD-nonclassical (Xu, 2022).

The structure and characterization of KD-classical pure states is governed by the following theorems:

• Theorem 1 (Algebraic criterion): $A = \{|a_j\rangle\}_{j=1}^d$7 is KD-classical (for $A = \{|a_j\rangle\}_{j=1}^d$8) if and only if there exist index sets $A = \{|a_j\rangle\}_{j=1}^d$9, real phases $B = \{|b_k\rangle\}_{k=1}^d$0, $B = \{|b_k\rangle\}_{k=1}^d$1, and positive amplitudes such that the phase of each nonzero transition matrix entry $B = \{|b_k\rangle\}_{k=1}^d$2 satisfies $B = \{|b_k\rangle\}_{k=1}^d$3 and $B = \{|b_k\rangle\}_{k=1}^d$4 admits both $B = \{|b_k\rangle\}_{k=1}^d$5- and $B = \{|b_k\rangle\}_{k=1}^d$6-basis expansions with the given supports and phases, ensuring all nonzero $B = \{|b_k\rangle\}_{k=1}^d$7 are positive real.
• Support-uncertainty bound: For any pure state, define $B = \{|b_k\rangle\}_{k=1}^d$8 and $B = \{|b_k\rangle\}_{k=1}^d$9. Then, $\rho$0 KD-classical implies $\rho$1; conversely, $\rho$2 necessitates KD nonclassicality for completely incompatible bases (Xu, 2022, Bievre, 2022, Bievre, 2021).

For the special case of mutually unbiased bases (MUBs)—i.e., $\rho$3 for all $\rho$4—KD-classicality is tightly characterized by $\rho$5 for pure states. This resolves the De Bièvre conjecture for the discrete Fourier transform (DFT) case (Xu, 2022).

3. Geometric and Convex Structure of KD-Positive States

The set of KD-positive states—states for which the KD distribution is a bona fide probability distribution—forms a convex, closed polytope. Its extremal points include the eigenstates of the respective observables and, in certain cases, additional states depending on the interplay of the bases (Langrenez et al., 2023, Langrenez et al., 2024). Explicit structural results include:

• Dimension and vertices: For $\rho$6, all bases, and for $\rho$7, generic bases, as well as for prime-dimensional DFT pairs, every KD-positive state is a convex combination of the $\rho$8 basis projectors (Langrenez et al., 2023, Xu, 2024).
• Block structure Theorem: For a pure KD-classical state, the nonzero-support submatrix of $\rho$9 can be partitioned into $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$0 nonoverlapping nonnegative blocks. The sum of corresponding supports is constrained: $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$1, and $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$2.
• Mixed KD-positive states: For certain non-prime or degenerate cases, there exist mixed KD-positive states not decomposable as convex mixtures of pure KD-classical states; these are characterized geometrically by extended faces of the KD-positive set (Langrenez et al., 2023, Langrenez et al., 2024).

Measures for KD non-positivity (negativity) include the total variation

$$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$3

and its convex-roof extension, which serves as a faithful witness for mixtures of pure KD-classical states (Langrenez et al., 2024).

4. Connections to Quantum Resources and Operational Applications

Negative or nonreal KD entries directly witness nonclassical behavior in various quantum resource paradigms:

• Quantum metrology: Negative KD-values enable quantum phase estimation protocols to exceed classical precision bounds via mechanisms such as weak-value amplification (Xu, 2022, Arvidsson-Shukur et al., 2024).
• Quantum chaos and scrambling: Out-of-time-ordered correlators (OTOCs) are expressible in terms of KD-like quasiprobabilities; robust negativity under open-system evolution is an indicator of information scrambling (Arvidsson-Shukur et al., 2024).
• Quantum thermodynamics: KD distributions for noncommuting energy measurements describe "work" distributions with correct marginals, even capturing effects due to initial coherence. Negativity signals the presence of contextuality and nonclassical heat/work flows (Arvidsson-Shukur et al., 2024, Pezzutto et al., 10 Mar 2025).
• Contextuality and foundations: KD negativity or nonreality underpins measurement-disturbance, weak-value anomalies, and violations of classical bounds in Leggett–Garg or consistent-histories frameworks. Nonnegativity of KD arrays is equivalent to the existence of noncontextual ontological models for the process; thus, KD negativity is a necessary signature of contextuality but not sufficient (Schmid et al., 2024, Spriet et al., 3 Nov 2025).

The KD pseudo-distribution's ability to encode quantum coherence and contextuality forms the basis for resource-theoretic monotones and quantifiers, such as KD mana for computational power (Thio et al., 9 Jun 2025) or entanglement monotones for bipartite systems (Budiyono, 7 Jan 2025). In the case of MUBs, KD-based coherence monotones reproduce and generalize established resource-theoretic quantities (Liu et al., 2024).

5. Measurement Protocols and Practical Witnesses

Direct measurement of the KD distribution is possible via weak-value protocols or systematic moment reconstruction:

• Weak measurement schemes: Weak measurement of one observable, followed by strong postselection on a complementary observable, samples conditional KD distributions; repeated over all settings reconstructs the full KD table (Jordan et al., 5 Feb 2026, Lostaglio et al., 2022). In discrete settings, Vandermonde-matrix inversion relates observable moments to KD entries.
• Moment-based witnesses: Simpler certification of nonclassicality can be achieved by checking statistical moment inequalities (e.g., $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$4 for the second and third KD moments, violation indicating negativity), or more generally, via negativity of Hankel determinants constructed from KD moments. These witnesses are experimentally friendly and require only shadow tomography or randomized measurements for efficient evaluation (Chakrabarty et al., 9 Jun 2025).
• Operational identification: KD negativity and nonreality map to accessible quantities such as average and variance in quantum thermodynamic protocols, direct state reconstruction, and resource monotones, thus providing both theoretical and experimental routes to quantifying nonclassicality.

6. Hierarchy and Limitations: Classical, KD, and Postquantum Sets

KD quasiprobabilities occupy a strict intermediary inside a broader hierarchy:

• Classical joint distributions $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$5 KD quasiprobabilities $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$6 postquantum quasiprobabilities: Every classical joint distribution is a KD distribution, but not vice versa; KD arrays obey strict sum, marginal, and $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$7 pointwise bounds, while postquantum distributions satisfy only normalization and a $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$8 bound (Liu et al., 12 Apr 2025).
• Universal norm bounds: For any KD distribution $$\mathrm{KD}_\rho(j,k) = \langle a_j|\rho|b_k\rangle\,\langle b_k|a_j\rangle.$$9, $\rho=|\psi\rangle\langle\psi|$0 and $\rho=|\psi\rangle\langle\psi|$1, providing global nonclassicality witnesses (Liu et al., 12 Apr 2025).
• Boundary cases: For two observables with complete incompatibility (e.g., MUBs in prime dimensions), KD-classical pure states are restricted solely to the basis eigenstates, and the set of KD-positive states forms their convex hull (Xu, 2024, Langrenez et al., 2023).

7. Conceptual Significance and Uniqueness

The KD pseudo-distribution is uniquely characterized among all Born-compatible quasiprobability representations: it is the only representation where the conditional expectation of an observable (conditioned on the outcome of a second, incompatible observable) coincides with the $\rho=|\psi\rangle\langle\psi|$2-optimal predictor by functions of the latter (Spriet et al., 3 Nov 2025). Only the KD representation sustains the "pull-through" property essential for statistical calculus of joint and conditional probabilities in the quantum regime.

The operational and mathematical completeness of the KD pseudo-distribution underpins its central role in modern quantum theory, quantifying the quantum-classical boundary, resource-theoretic utility, and foundational structure of quantum processes.