---
title: Proper Quantum Scoring Rules Overview
url: https://www.emergentmind.com/topics/proper-quantum-scoring-rules
type: topic
---

# Proper Quantum Scoring Rules Overview

Proper quantum scoring rules are the noncommutative analogue of classical proper scoring rules: they are scoring mechanisms for quantum-state reports under which truthful reporting of the underlying density operator is optimal in expectation. In the finite-dimensional setting, the literature develops two closely related formulations. One makes the measurement stage explicit, defining a quantum score as a pair consisting of a payoff rule and a report-dependent measurement, with expected score \(S(\hat\rho;\rho)\) computed from the outcome distribution induced by measuring the true mixed state \(\rho\) after a report \(\hat\rho\) [2203.07469]. The other emphasizes an operator-valued state-level rule \(\mathbf S:D(H)\to \mathcal L(H)_{\rm Herm}\), evaluated directly by trace pairing as \(\operatorname{Tr}(\rho\,\mathbf S(\sigma))\), and derives properness from a convex unitarily invariant Quantum Value Functional \(V_Q(\rho)=\operatorname{Tr}(f(\rho))\) [2605.05268]. Across both formulations, the central structural result is a quantum version of the classical convex-analytic correspondence between proper scores, supporting hyperplanes, and Bregman-type divergences.

## 1. Concept and formal definitions

In the measurement-explicit formulation, a quantum mixed state is a density matrix \(\rho\in \mathrm{Dens}(X)\), where \(X=\mathbb C^n\), and a quantum score is a pair \(S=(s,\mu)\). Here \(\mu(\hat\rho)\) is a POVM chosen as a function of the reported state \(\hat\rho\), and the expected score is
\[
S(\hat\rho;\rho):=\mathbb E_{Y\sim \mu(\hat\rho)\rho}\, s(\hat\rho,Y).
\]
The score is called truthful if
\[
S(\rho;\rho)\ge S(\hat\rho;\rho)\quad\text{for all }\rho,\hat\rho\in \mathrm{Dens}(X),
\]
and strictly truthful if the inequality is strict whenever \(\hat\rho\neq \rho\) [2203.07469]. In the terminology used elsewhere, truthful and strictly truthful are the quantum analogues of proper and strictly proper scoring rules.

The operator-valued formulation replaces the classical scalar score \(S(x,q)\) with an operator-valued rule
\[
\mathbf S:D(H)\to \mathcal L(H)_{\rm Herm},
\]
defined on density operators \(\rho,\sigma\in D(H)\) over a finite-dimensional complex Hilbert space \(H\). The realized scalar score against the true state \(\rho\) is
\[
\operatorname{Tr}\big(\rho\,\mathbf S(\sigma)\big),
\]
and properness is the condition
\[
\operatorname{Tr}\big(\rho\,\mathbf S(\rho)\big)\ge \operatorname{Tr}\big(\rho\,\mathbf S(\sigma)\big)\quad \forall \rho,\sigma\in D(H).
\]
The 2026 formulation identifies the primitive convex object as a Quantum Value Functional \(V_Q:D(H)\to \mathbb R\) that is convex, unitarily invariant, and lower semicontinuous; by spectral calculus,
\[
V_Q(\rho)=\sum_i f(\lambda_i)=\operatorname{Tr}(f(\rho))
\]
for a convex \(f:[0,1]\to \mathbb R\) [2605.05268].

A notable conceptual difference from the classical case is that a density matrix does not directly generate an observable classical outcome. The 2022 framework therefore includes the measurement rule as part of the contract, while the 2026 framework packages the dependence on the report into the operator-valued score and evaluates it by trace pairing against the true state. This suggests two complementary viewpoints: one operational and measurement-level, the other convex-analytic and state-level.

## 2. Convex-analytic structure and duality

A central result of the finite-dimensional theory is that expected score is affine in the true state. In the measurement-based framework,
\[
S(\hat\rho;\rho)=\langle Z_S(\hat\rho),\rho\rangle,
\qquad
Z_S(\hat\rho)=\sum_{y\in Y(\hat\rho)} \mu(\hat\rho)_y\, s(\hat\rho,y),
\]
so every quantum score induces a linear functional of \(\rho\) for fixed report \(\hat\rho\) [2203.07469]. This is the quantum analogue of the classical fact that expected score is linear in the true distribution.

The main characterization theorem in that framework states that a quantum score \(S\) is truthful if and only if there exists a convex function \(F:\mathrm{Dens}(X)\to \mathbb R\) and a choice of subgradients \(d_\rho\in \partial F(\rho)\) such that
\[
S(\hat\rho;\rho)=F(\hat\rho)+\langle d_{\hat\rho},\,\rho-\hat\rho\rangle.
\]
Strict truthfulness corresponds to strict convexity. The truthfulness gap is then the generalized Bregman divergence
\[
S(\rho;\rho)-S(\hat\rho;\rho)
=
F(\rho)-F(\hat\rho)-\langle d_{\hat\rho},\rho-\hat\rho\rangle
\ge 0
\]
[2203.07469].

The 2026 theory states an analogous correspondence in terms of Quantum Value Functionals. For differentiable spectral functionals,
\[
\mathbf S(\sigma)=\nabla V_Q(\sigma)=f'(\sigma),
\]
and the induced quantum Bregman divergence is
\[
D_{V_Q}(\rho\|\sigma)
=
V_Q(\rho)-V_Q(\sigma)-\operatorname{Tr}\!\big(\mathbf S(\sigma)(\rho-\sigma)\big).
\]
Convexity yields \(D_{V_Q}(\rho\|\sigma)\ge 0\), and the paper uses the identity
\[
\operatorname{Tr}(\rho\mathbf S(\rho))-\operatorname{Tr}(\rho\mathbf S(\sigma))\ge D_{V_Q}(\rho\|\sigma)\ge 0
\]
to encode properness [2605.05268].

The duality theorem in the 2026 formulation centers on the correspondence
\[
\text{convex unitarily invariant }V_Q
\quad\longleftrightarrow\quad
\text{proper quantum scoring rule }\mathbf S=\partial V_Q.
\]
The explicit theorem establishes the forward direction: if \(V_Q\) is closed, convex, and unitarily invariant, then \(\mathbf S(\sigma)=\nabla V_Q(\sigma)\) is proper. The abstract describes this as a “complete duality theory,” but the detailed statement presented there formalizes the forward construction most explicitly [2605.05268].

## 3. Spectral structure, operator convexity, and divergences

For spectral functionals \(V_Q(\rho)=\operatorname{Tr}(f(\rho))\), the Fréchet derivative can be computed through the Daleckii–Krein formula. If \(\sigma=U\Lambda U^\dagger\), then
\[
DV_Q(\sigma)[H]
=
\operatorname{Tr}\!\left(
U\big(f^{[1]}(\Lambda)\circ (U^\dagger H U)\big)U^\dagger
\right),
\]
with first divided difference matrix
\[
f^{[1]}_{ij}
=
\begin{cases}
\dfrac{f(\lambda_i)-f(\lambda_j)}{\lambda_i-\lambda_j}, & \lambda_i\neq \lambda_j,\\[1ex]
f'(\lambda_i), & \lambda_i=\lambda_j.
\end{cases}
\]
Identifying the Hilbert–Schmidt gradient yields \(\mathbf S(\sigma)=f'(\sigma)\) [2605.05268].

The operator-convex case is singled out because it provides stronger order-theoretic and information-theoretic structure. If \(f\) is operator convex, then
\[
f(\rho)\succeq f(\sigma)+f'(\sigma)(\rho-\sigma)
\]
in Löwner order. Under the additional symmetry condition that
\[
\sigma^{-1/2}\rho\sigma^{-1/2}
\quad\text{commutes with }\sigma,
\]
the induced Bregman divergence coincides with Petz’s quantum \(f\)-divergence,
\[
D_f(\rho\|\sigma)
=
\operatorname{Tr}\!\left[
\sigma^{1/2} f\!\left(\sigma^{-1/2}\rho\sigma^{-1/2}\right)\sigma^{1/2}
\right]
\]
[2605.05268].

The 2022 theory isolates a particularly elegant subclass, spectral scores. Writing a reported state as
\[
\hat\rho=\sum_{i=1}^n \lambda_i(\hat\rho)\,|v_i\rangle\langle v_i|,
\]
a spectral score measures in the reported eigenbasis,
\[
\mu(\hat\rho)_i=|v_i\rangle\langle v_i|,
\]
and pays according to a classical scoring rule on the eigenvalue vector,
\[
s(\hat\rho,i)=\hat s(\lambda(\hat\rho),i).
\]
Its expected score is
\[
S(\hat\rho;\rho)=\sum_{i=1}^n \langle v_i|\rho|v_i\rangle\,\hat s(\lambda(\hat\rho),i).
\]
The paper proves that a spectral quantum score is truthful if and only if the underlying classical scoring rule \(\hat s\) is proper, and that a finite truthful quantum score is equivalent to a spectral score precisely when it is unitary-invariant [2203.07469].

These two lines of work converge on the same structural theme: unitary symmetry singles out spectral convex potentials and spectral proper scores, while noncommutativity forces one to distinguish ordinary convexity from operator convexity and to treat divergence notions that coincide classically as inequivalent in the quantum setting.

## 4. Relation to classical proper scoring rules

The quantum theories preserve the classical pattern with systematic replacements. The classical simplex of distributions becomes the convex set \(D(H)\) or \(\mathrm{Dens}(X)\); scalar probabilities become density matrices; Euclidean pairings become Hilbert–Schmidt trace pairings; and classical gradients or subgradients become Fréchet derivatives or matrix subgradients [2605.05268].

In the commuting case, the quantum theory reduces to the classical one. If \(\rho\) and \(\sigma\) commute, they are simultaneously diagonalizable:
\[
\rho=U\operatorname{diag}(p_i)U^\dagger,
\qquad
\sigma=U\operatorname{diag}(q_i)U^\dagger.
\]
Then
\[
V_Q(\rho)=\sum_i f(p_i),
\qquad
\mathbf S(\sigma)=U\operatorname{diag}(f'(q_i))U^\dagger,
\]
and
\[
\operatorname{Tr}(\rho\,\mathbf S(\sigma))=\sum_i p_i f'(q_i).
\]
The corresponding quantum Bregman divergence becomes
\[
D_{V_Q}(\rho\|\sigma)
=
\sum_i \big(f(p_i)-f(q_i)-f'(q_i)(p_i-q_i)\big),
\]
which is exactly the classical Bregman divergence [2605.05268].

The connection to classical convex-geometric results is not merely formal. The classical paper “Proper Scoring Rules and Domination” generalizes a theorem of Predd et al. to non-additive scoring rules on a finite state space, proving that under strict propriety and a closure or continuity condition, every non-probabilistic credence is pointwise dominated by some probabilistic one [2102.02260]. That work is entirely classical, but its relevance to the quantum setting lies in its emphasis that additivity is not essential: propriety can be studied through convex geometry of the global score set rather than coordinatewise decompositions. The 2022 and 2026 quantum theories likewise rely on convex structure, linearity in the true state, and non-additive state-level scoring objects rather than classical eventwise additivity.

A further distinction is that the 2022 framework makes report-dependent measurement explicit. Classically, one directly observes an outcome drawn from the true distribution; quantum mechanically, the report may also determine how the outcome distribution is generated. This report-dependent measurement is the main conceptual novelty of the operational formulation [2203.07469].

## 5. Estimation, tomography, and minimax theory

A major extension of proper quantum scoring rules is their use as loss functions for quantum state estimation. In the 2026 framework, the loss of estimating \(\rho\) by \(\hat\rho_n\) is
\[
\operatorname{Tr}(\rho\mathbf S(\rho))-\operatorname{Tr}(\rho\mathbf S(\hat\rho_n)),
\]
and the corresponding minimax risk over \(n\) copies is
\[
R_n^*(V_Q)
=
\inf_{\hat\rho_n}\sup_{\rho\in D(H)}
\mathbb E_{\rho^{\otimes n}}
\!\left[
\operatorname{Tr}(\rho\mathbf S(\rho))
-
\operatorname{Tr}(\rho\mathbf S(\hat\rho_n))
\right].
\]
This recasts tomography as minimax estimation under incentive-compatible scoring-rule regret [2605.05268].

The asymptotic lower bound derived there is the Quantum Cramér–Rao–McCarthy Bound,
\[
R_n^*(V_Q)
\ge
\frac{1}{2n}
\sup_{\rho\in D(H)}
\operatorname{Tr}\!\left(
f''(\rho)\cdot I_{\mathrm{SLD}(\rho)}^{-1}
\right)
+
o(n^{-1}).
\]
The ingredients are the spectral curvature \(f''(\rho)\) and the Symmetric Logarithmic Derivative Quantum Fisher Information. If \(\partial_i\rho_\theta=\frac12(\rho_\theta L_i+L_i\rho_\theta)\), then
\[
(I_{\mathrm{SLD}})_{ij}=\operatorname{Tr}(\rho_\theta L_iL_j).
\]
The theorem therefore links minimax risk under McCarthy-type incentives to local information geometry and the curvature of the generating functional [2605.05268].

The mechanism of the bound is a second-order expansion of the Bregman divergence. For \(V_Q(\rho)=\operatorname{Tr}(f(\rho))\),
\[
D_{V_Q}(\rho_\theta\|\hat\rho_n)
\approx
\frac12(\hat\theta_n-\theta)^T \nabla_\theta^2 V_Q(\theta)(\hat\theta_n-\theta),
\]
with second variation
\[
\delta^2 V_Q(\rho)[H,H]
=
\sum_{i,j} f^{[2]}(\lambda_i,\lambda_j)\,|\langle u_i,H u_j\rangle|^2.
\]
In the asymptotic approximation used there, this Hessian is represented by \(f''(\rho)\) in the eigenbasis of \(\rho\), so local risk behaves like
\[
\frac{1}{2n}\operatorname{Tr}\!\left(f''(\rho)I_{\mathrm{SLD}}^{-1}\right).
\]
A sharper curvature therefore corresponds to harsher local penalization of estimation error [2605.05268].

For the logarithmic generator \(f(t)=t\log t\), one has \(f''(t)=1/t\), hence \(f''(\rho)=\rho^{-1}\) on the support of \(\rho\). The paper states that asymptotic equality is achieved in this case by the maximum likelihood estimator when \(\rho\) is nearly pure or the parameterization is locally orthogonal in the QFI metric [2605.05268].

## 6. Examples, resources, and scope

The canonical example is the quantum logarithmic score. In the spectral functional approach, for
\[
f(t)=t\log t,
\]
the induced score is
\[
\mathbf S(\sigma)=\log \sigma+\mathbb I,
\]
and the truthful expected score is
\[
\operatorname{Tr}\big(\rho(\log\rho+\mathbb I)\big)
=
\operatorname{Tr}(\rho\log\rho)+1
=
-S(\rho)+1,
\]
where \(S(\rho)=-\operatorname{Tr}(\rho\log\rho)\) is von Neumann entropy [2605.05268]. In the 2022 spectral-score formulation, the corresponding expected score is
\[
S(\hat\rho;\rho)=\operatorname{Tr}(\rho\log \hat\rho),
\]
the truthful value is \(\operatorname{Tr}(\rho\log\rho)\), and the truthfulness gap is
\[
S(\rho;\rho)-S(\hat\rho;\rho)
=
\operatorname{Tr}\bigl(\rho(\log\rho-\log\hat\rho)\bigr)
=
D(\rho\|\hat\rho),
\]
giving an operational interpretation of von Neumann relative entropy as the regret from misreporting under the logarithmic score [2203.07469].

A quadratic analogue also appears naturally. With \(f(t)=t^2\),
\[
V_Q(\rho)=\operatorname{Tr}(\rho^2),
\qquad
\mathbf S(\sigma)=2\sigma,
\]
which is the Hilbert–Schmidt quadratic member of the spectral class [2605.05268]. The 2022 paper also gives a Brier-type realization with expected score
\[
S(\hat\rho;\rho)=2\langle \hat\rho,\rho\rangle-\langle \hat\rho,\hat\rho\rangle,
\]
which equals \(S(\rho;\rho)-\|\hat\rho-\rho\|_2^2\) and is therefore strictly truthful [2203.07469].

The qubit coherence example offers a concrete operational illustration. For
\[
\rho=|+\rangle\langle +|
=
\frac12
\begin{pmatrix}
1&1\\
1&1
\end{pmatrix}
=
\frac12(\mathbb I+\sigma_x),
\]
a classical forecaster restricted to the computational \(Z\)-basis effectively sees \(\Delta_Z(\rho)=\mathbb I/2\), while a quantum forecaster measuring in the \(X\)-basis identifies the state perfectly. Under the logarithmic score, the truthful quantum report yields
\[
\operatorname{Tr}(\rho\mathbf S(\rho))=1,
\]
whereas the best diagonal classical report gives
\[
\operatorname{Tr}(\rho\mathbf S(\sigma))
=
\log \frac12+1.
\]
The forecasting gap is therefore
\[
\Gamma
=
1-\left(\log\frac12+1\right)
=
\log 2,
\]
which matches the relative entropy of coherence
\[
C_Z(\rho)=S(\Delta_Z(\rho))-S(\rho)=\log 2
\]
[2605.05268].

More generally, the 2026 work assigns operational value to quantum resources in proper-scoring-rule risk units. For noncommuting output channels it proves a separation between classical fixed-basis strategies and fully quantum joint-measurement strategies:
\[
\Gamma_n = \Omega\!\left(\frac{d\log d}{n}\right)-O(n^{-3/2}),
\]
and states that quantum strategies can reduce sample complexity by a factor \(O(d)\) in the worst case. For the logarithmic score, it also gives the coherence-risk tradeoff
\[
R_n^{\text{classical}}
\ge
R_n^{\text{quantum}}
+
\frac{C(\rho)}{n}
+
O(n^{-3/2}),
\]
provided the dephased state \(\Delta(\rho)\) is the information projection of \(\rho\) onto the diagonal submanifold with respect to the relevant Bregman geometry. The paper summarizes qualitative resource scalings as coherence \(O(C\cdot d/n)\), entanglement \(O(E\cdot d^2/n)\), and adaptivity \(O(d^2\log d/n)\) [2605.05268].

The present theory is explicitly finite-dimensional and state-based. Its assumptions include density operators as the report space, unitarily invariant convex generators, and—when asymptotic estimation results are derived—i.i.d. copies, smooth local parameterization, differentiability, \(f''>0\), and LAN/Holevo/SLD asymptotic machinery [2605.05268]. The 2022 framework likewise remains finite-dimensional, uses finite POVMs or PVMs, and permits extended-real scores to accommodate logarithmic boundary behavior [2203.07469]. Reported limitations include asymptotic rather than exact nonasymptotic analysis for much of the minimax theory, commutation symmetry conditions in some divergence-identification statements, restriction to states rather than channels or processes, omission of computational constraints, and the absence of an infinite-dimensional theory [2605.05268]. Open directions highlighted in the literature include extensions to quantum channels, more systematic treatment of adaptive protocols, experimental validation on near-term devices, strategic interaction among multiple quantum agents, non-i.i.d. settings, infinite-dimensional systems, and robustness under noisy channels using optimized \(f\)-divergences [2605.05268].

Source: https://www.emergentmind.com/topics/proper-quantum-scoring-rules