Terletsky–Margenau–Hill Quasiprobability Approach

• The Terletsky–Margenau–Hill approach is a quasiprobability framework that extends classical joint probability to noncommuting quantum observables through symmetrized operator ordering.
• It reproduces exact quantum marginals and moments, providing a rigorous tool to detect quantum coherence, contextuality, and measurement incompatibility.
• Applications include quantum thermodynamics and work distribution analyses, where the approach offers experimental access to coherent effects via its inherent negativity.

The Terletsky–Margenau–Hill (TMH or MH) quasiprobability approach provides a mathematically rigorous and physically motivated framework for generalizing classical joint probability distributions to the quantum regime, particularly for noncommuting observables and sequential measurements. Developed originally by Terletsky, Margenau, and Hill, the MH construction gives a real-valued but sign-indefinite (“quasi”) distribution whose marginals and moments coincide with symmetrized operator orderings, and whose negativity directly signals quantum coherence, incompatibility, and contextuality. This approach has found foundational and practical use in quantum thermodynamics, quantum measurement theory, descriptions of energy distributions, process tomography, and studies of nonclassicality emerging in quantum dynamics.

1. Formal Structure of the Margenau–Hill Quasiprobability

The TMH quasiprobability associates to any pair (or finite tuple) of Hermitian operators (observables) a signed measure that reproduces the exact quantum marginals and (for symmetrized polynomials) quantum expectation values. For observables $A$, $B$ with spectral projectors $\Pi_a^A$, $\Pi_b^B$, and system density operator $\rho$, the two-point Margenau–Hill quasiprobability is defined as

$$q_{MH}(a, b) = \frac{1}{2} \mathrm{Tr}\left[(\Pi_b^B \Pi_a^A + \Pi_a^A \Pi_b^B)\, \rho\right] = \Re\, \mathrm{Tr}(\Pi_b^B \Pi_a^A \rho)$$

This function is real-valued, linear in $\rho$, normalized ($\sum_{a,b} q_{MH}(a,b) = 1$), and recovers the Born rule probabilities as marginals: $$\sum_a q_{MH}(a, b) = \mathrm{Tr}(\Pi_b^B \rho), \quad \sum_b q_{MH}(a, b) = \mathrm{Tr}(\Pi_a^A \rho)$$ Negativity of $q_{MH}$ is possible if $B$0, $B$1, or $B$2 do not mutually commute, and witnesses quantum nonclassicality in the statistics of sequential or joint measurements (Bizzarri et al., 2024, Shukla et al., 1 Jan 2026).

The construction is readily generalized to $B$3 observables $B$4 via the fully symmetrized Margenau–Hill correspondence, associating classical monomials with their fully symmetrized operator products. The corresponding multi-point characteristic function is

$B$5

whose inverse Fourier transform yields $B$6 (Vasudevrao et al., 2021, Sabbagh et al., 2024).

2. Physical Requirements and Uniqueness

Among various quantum quasiprobability proposals (e.g., Wigner, Kirkwood–Dirac, Full Counting statistics), the MH form is uniquely characterized by five stringent physical conditions:

1. Linearity and normalization: $B$7 for some operator measure $B$8, and $B$9.
2. Support: Distribution is supported only on genuine transitions, i.e., differences of eigenvalues of the relevant observables.
3. First-law identity: The mean value reproduces quantum expectation; for work, $\Pi_a^A$0.
4. Time-reversal symmetry: The distribution is invariant under reversal of drive/unitary protocols.
5. Positivity of variance: The second moment is always nonnegative, $\Pi_a^A$1 for all $\Pi_a^A$2.

Only the symmetrized MH kernel—i.e., $\Pi_a^A$3—survives imposition of all of these (Pei et al., 2023).

Alternative quasiprobabilities like the Kirkwood–Dirac (KD) distribution generally violate at least positivity of the second moment, and Full Counting statistics may fail time-reversal symmetry or yield unphysical support (Pei et al., 2023, Díaz et al., 2020).

3. Applications in Quantum Work and Thermodynamics

The MH quasiprobability has been most extensively adopted for the definition of work distributions in quantum thermodynamics, where standard two-projective-measurement (TPM) protocols fail in the presence of initial coherence. The MH distribution delivers a physically meaningful and experimentally reconstructable quasiprobability for quantum work,

$\Pi_a^A$4

ensuring correct average work and variance, and reducing to TPM when $\Pi_a^A$5 (Pei et al., 2023, Díaz et al., 2020).

The characteristic function formalism ensures that moments of $\Pi_a^A$6 can be computed via symmetrized polynomials: $\Pi_a^A$7 The influence of coherence is explicit: all differences between MH and TPM results (for mean and variance) are controlled by the $\Pi_a^A$8-coherence $\Pi_a^A$9 in the initial energy basis. Notably, the MH approach admits negative average entropy production in certain regimes, signaling reversible protocols enabled by quantum coherence (Díaz et al., 2020).

The resulting distribution exhibits convergence to the classical, positive, Liouville work distribution in the $\Pi_b^B$0 limit. This is manifest in exactly solvable models such as the breathing harmonic oscillator, where the discrete MH distribution increasingly approximates the continuous classical outcome (Pei et al., 2023).

4. Contextuality, Negativity, and Temporal Extensions

Negativity in the TMH quasiprobability is directly tied to quantum contextuality and is operationally significant: any observed $\Pi_b^B$1 signals the breakdown of any noncontextual hidden-variable model compatible with the measurement sequence (Shukla et al., 1 Jan 2026, Bizzarri et al., 2024). This diagnosis is sharpened via metrics such as the total negativity $\Pi_b^B$2 and the minimal time of appearance of negativity in dynamical scenarios (“first-time negativity”, FTN), which distinguishes coherent quantum dynamics from classical or decohered ones.

The temporal extension of MH quasiprobabilities is formalized by considering multi-time measurement protocols and associating the MH distribution with the real part of temporal Kirkwood–Dirac (KD) quasiprobabilities: $\Pi_b^B$3 Here, $\Pi_b^B$4 denotes the ordered sequence of projective outcomes under possibly non-commuting evolution channels, and $\Pi_b^B$5 is experimentally accessible via interferometric or Bloch tomography protocols (Jia et al., 8 Jan 2026).

These temporal (and spatiotemporal) MH distributions satisfy key consistency properties (reduction to marginals, symmetry), and provide a unifying operational construct for a variety of “temporal state” and pseudo-density operator formalisms.

5. Joint Measurability, Measurement Incompatibility, and Symmetrization

The MH approach naturally encodes the boundaries of joint measurability for sets of noncommuting observables. By introducing an unsharpness parameter $\Pi_b^B$6—which modulates the degree of “fuzzification”—the transition from quasiprobability to proper probability is marked by the positivity of all elements of the associated operator-valued measure. This yields precise and sometimes tight bounds on joint measurability in finite-dimensional systems:

• For orthogonal qubit observables, $\Pi_b^B$7 is necessary and sufficient for joint measurability.
• For qutrits or two-qubit composite systems, stricter bounds such as $\Pi_b^B$8 are sufficient (Vasudevrao et al., 2021).

Increased Hilbert space dimension and operator noncommutativity generically move systems away from joint measurability, reflected in persistent negativity of $\Pi_b^B$9. Thus, MH negativity offers an explicit quantification of measurement incompatibility, with direct operational meaning.

6. Relation to Other Quasiprobability Constructs

While the Wigner and Kirkwood–Dirac quasiprobabilities are more commonly employed in continuous-variable and foundational scenarios, the MH approach is distinguished by three key features:

• It is always real-valued; its negativity is a sufficient witness of nonclassicality, without the ambiguities associated with complex-valued KD distributions or the marginalization peculiarities of Wigner functions (Jia et al., 8 Jan 2026, Sabbagh et al., 2024).
• Its support is explicitly tied to the spectra of the underlying operators.
• It builds a rigorous bridge between discrete, spectrally supported quasimeasures and the continuous Wigner formalism, with the MH particle approximations converging to the Wigner distribution (as described via the Lie–Trotter formula and classical Mehler–Heine limits for polynomials) (Sabbagh et al., 2024).

The MH approach also coincides with the real part of weak value formulas and is naturally associated with the symmetric Hadamard product in operator orderings. As a result, measurement protocols based on weak measurements or temporal tomography can access the MH distribution directly (Bizzarri et al., 2024, Jia et al., 8 Jan 2026).

7. Energy Distribution, Transport, and Thermodynamic Uncertainty

The TMH quasiprobability delivers a unique resolution for defining local quantum energy densities, especially in settings where coordinate and energy operators do not commute. Derivations from the Dirac energy–momentum tensor establish that the only physically consistent definition of local energy density (in the nonrelativistic limit) is via the TMH marginal of energy and position operators: $\rho$0 where $\rho$1 is the TMH phase-space quasiprobability (Stepanyan et al., 2023).

This formulation coincides with the real part of the weak value of the energy observable and recovers the classical hydrodynamic (Madelung) energy including quantum potential corrections. It is locally conserved and, in systems with spin, reveals additional “holographic” energy contributions associated with boundary terms.

In open quantum systems, the TMH quasiprobability quantifies dynamical fluctuations and forms the basis for quantum thermodynamic uncertainty relations (TURs). The variance of observable changes, captured through the TMH moments, bounds the entropy production rate. Only negativity in the TMH distribution enables certain regimes (e.g., dissipationless heat currents) unattainable in classical systems, highlighting its role as a stronger criterion than mere presence of coherence (Yoshimura et al., 20 Aug 2025).

---

References:

• "Exploring quasiprobability approach to quantum work in the presence of initial coherence: Advantages of the Margenau-Hill distribution" (Pei et al., 2023)
• "Quantum work statistics with initial coherence" (Díaz et al., 2020)
• "Margenau-Hill operator valued measures and joint measurability" (Vasudevrao et al., 2021)
• "Particle approximations of Wigner distributions for n arbitrary observables" (Sabbagh et al., 2024)
• "First appearance of quasiprobability negativity in quantum many-body dynamics" (Shukla et al., 1 Jan 2026)
• "Quasiprobability distributions with weak measurements" (Bizzarri et al., 2024)
• "Energy densities in quantum mechanics" (Stepanyan et al., 2023)
• "Quasiprobability Thermodynamic Uncertainty Relation" (Yoshimura et al., 20 Aug 2025)
• "Temporal Kirkwood-Dirac Quasiprobability Distribution and Unification of Temporal State Formalisms through Temporal Bloch Tomography" (Jia et al., 8 Jan 2026)