Nussbaum–Szkoła Distributions in Quantum Testing

• Nussbaum–Szkoła distributions are classical representations derived from quantum states’ spectral decompositions, pairing eigenvalues with basis-overlap intensities.
• They enable the precise translation of quantum quantities—such as Chernoff overlaps and Petz-type f-divergences—into equivalent classical forms for error analysis.
• Their application spans from finite-dimensional systems to separable Hilbert spaces, advancing quantum hypothesis testing and geometric state discrimination.

Searching arXiv for recent and foundational papers on Nussbaum–Szkoła distributions to ground the article. Nussbaum–Szkoła distributions are a canonical classical representation associated with a pair of quantum states. Given spectral decompositions of two states, they assign probability weights on pairs of eigen-indices by combining eigenvalues with basis-overlap intensities. Their defining significance is that key quantum distinguishability quantities become exactly classical for the resulting distributions. In particular, for two states $\rho$ and $\sigma$, the Chernoff overlap $\operatorname{Tr}(\rho^s\sigma^{1-s})$ is reproduced by the classical overlap $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$, and more generally a broad class of Petz-type quantum $f$-divergences equals the corresponding classical $f$-divergence of the Nussbaum–Szkoła pair (Li, 2015, Androulakis et al., 2023). This construction has become a central analytical device in quantum hypothesis testing, asymptotic error analysis, and the study of quantum divergences, with extensions from finite-dimensional matrix algebras to separable Hilbert spaces and semifinite von Neumann algebras (Cheng et al., 4 Jun 2026, Anastasiadis et al., 21 Apr 2026).

1. Definition and basic construction

For two finite-dimensional quantum states $\rho$ and $\sigma$ with spectral decompositions

$$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$

the Nussbaum–Szkoła distributions are defined on index pairs $(i,j)$ by

$\sigma$0

The same construction is described in equivalent notation in later works, for example

$\sigma$1

for spectral decompositions $\sigma$2 and $\sigma$3 (Androulakis et al., 2023, Roy et al., 8 Apr 2026).

Normalization follows from the orthonormality of the eigenbases: $\sigma$4 for normalized states. In a more general trace-class setting, one has $\sigma$5 and $\sigma$6 (Androulakis et al., 2023, Cheng et al., 4 Jun 2026).

The overlap matrix

$\sigma$7

is doubly stochastic: $\sigma$8 This matrix encodes the relative orientation of the eigenbases. Degeneracies are handled by choosing any orthonormal basis inside each degenerate eigenspace. The resulting NS distributions may depend on that basis choice, but several divergence values computed from them do not (Androulakis et al., 2023).

In separable Hilbert spaces, the same construction remains discrete because trace-class density operators admit countable spectral decompositions. Continuous spectra do not arise for trace-class density operators in that setting (Androulakis et al., 2023).

2. Chernoff identity and the original hypothesis-testing role

The defining identity behind the hypothesis-testing applications is

$\sigma$9

As a consequence, the classical Chernoff information of $\operatorname{Tr}(\rho^s\sigma^{1-s})$0,

$\operatorname{Tr}(\rho^s\sigma^{1-s})$1

exactly matches the quantum Chernoff quantity

$\operatorname{Tr}(\rho^s\sigma^{1-s})$2

This identity is the classical backbone of the Nussbaum–Szkoła lower-bound technique in quantum hypothesis testing (Li, 2015).

The same point is expressed in later formulations through the equality of quantum and classical Chernoff distances: $\operatorname{Tr}(\rho^s\sigma^{1-s})$3 and

$\operatorname{Tr}(\rho^s\sigma^{1-s})$4

with $\operatorname{Tr}(\rho^s\sigma^{1-s})$5 (Cheng et al., 4 Jun 2026).

This exact transfer from quantum to classical Chernoff overlaps explains why NS distributions entered the subject through quantum state discrimination. Classical large-deviation and Chernoff machinery can be applied to $\operatorname{Tr}(\rho^s\sigma^{1-s})$6, and the resulting statements then yield quantum lower bounds or exact asymptotic exponents. In the binary i.i.d. setting, this underlies the quantum Chernoff bound. In the multi-hypothesis setting, it supports the passage from pairwise exponents to the multiple quantum Chernoff distance (Li, 2015).

A related property emphasized in more recent work is that NS distributions preserve Petz Rényi divergences for $\operatorname{Tr}(\rho^s\sigma^{1-s})$7, and this supports the transfer of Hoeffding- and Stein-type asymptotic analyses to the classical side as well (Cheng et al., 4 Jun 2026).

3. Multiple-state discrimination and the Nussbaum–Szkoła conjecture

For an arbitrary finite ensemble $\operatorname{Tr}(\rho^s\sigma^{1-s})$8 on a finite-dimensional Hilbert space and an arbitrary prior $\operatorname{Tr}(\rho^s\sigma^{1-s})$9 independent of $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$0, one tests the i.i.d. hypotheses $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$1. The optimal average error probability satisfies

$\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$2

with

$\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$3

This resolves the conjecture of Nussbaum and Szkoła that the optimal asymptotic error exponent for multiple quantum hypotheses equals the minimum over pairwise Chernoff distances (Li, 2015).

In that proof, NS distributions enter through a lower-bound mechanism: each quantum pair $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$4 is mapped to classical NS distributions $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$5, classical Chernoff bounds are applied to those pairs, and appropriate pairwise contributions are aggregated. The paper crystallizes this through a one-shot lower bound

$\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$6

for spectral decompositions $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$7. For density matrices with rank-one spectral projectors, the overlap term reduces to $\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$8, and the summand matches the classical quantity appearing in the NS construction (Li, 2015).

The same line of development was later sharpened in arbitrary separable Hilbert spaces. A dimension-free one-shot upper bound in terms of pairwise errors shows

$\sum_{i,j} P_{ij}^s Q_{ij}^{1-s}$9

where $f$0 is the NS pair for $f$1. In the i.i.d. regime this yields

$f$2

with

$f$3

This removes the dimension-dependent prefactor present in earlier finite-dimensional results and establishes achievability of the multiple Chernoff distance in arbitrary separable Hilbert spaces (Cheng et al., 4 Jun 2026).

These results jointly situate NS distributions as the pairwise classical objects governing the asymptotics of multiple quantum testing: they provide the lower-bound exponent, guide the structure of pairwise decompositions, and remain effective beyond finite dimensions (Li, 2015, Cheng et al., 4 Jun 2026).

4. Equality for quantum $f$4-divergences

A major generalization of the NS framework shows that, for a very broad class of quantum $f$5-divergences, the quantum quantity equals the classical $f$6-divergence of the corresponding NS distributions. Let $f$7 be convex or concave, and define the classical Csiszár $f$8-divergence by

$f$9

with the standard conventions for zero entries. For density operators $f$0, the quantum $f$1-divergence considered in (Androulakis et al., 2023) is the Petz-type quasi-entropy defined through the relative modular operator $f$2: $f$3 The main theorem states that

$f$4

where $f$5 are the NS distributions of $f$6 (Androulakis et al., 2023).

This equality holds in finite and infinite dimensions, without faithfulness assumptions, and support mismatches are accounted for by boundary terms on both sides. The equivalence

$f$7

determines when the divergence reduces to the pure-ratio form without boundary terms (Androulakis et al., 2023).

The mechanism behind the equality is spectral. The relative modular operator acts on rank-one operators $f$8 with eigenvalues $f$9, and the spectral measure weight tested against $\rho$0 is precisely $\rho$1. Thus the operator functional calculus reduces to the scalar summation defining the classical $\rho$2-divergence of the NS pair (Androulakis et al., 2023).

This yields direct quantum versions of many classical inequalities. Representative examples transferred in (Androulakis et al., 2023) include: $\rho$3

$\rho$4

and

$\rho$5

A Pinsker-type continuity statement is also obtained: for strictly convex $\rho$6 with $\rho$7, there exists $\rho$8 with $\rho$9 such that

$\sigma$0

and therefore $\sigma$1 implies $\sigma$2 (Androulakis et al., 2023).

A later extension proves the same reduction for normal states on a semifinite von Neumann algebra. If $\sigma$3 is semifinite with faithful normal semifinite trace $\sigma$4, and $\sigma$5 are normal states, then there exist a $\sigma$6-finite measured space $\sigma$7 and nonnegative measurable functions $\sigma$8 such that

$\sigma$9

for every convex $$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$0. In finite dimensions this recovers the familiar discrete formulas $$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$1 and $$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$2 (Anastasiadis et al., 21 Apr 2026).

5. Structural and geometric perspectives

The NS mapping has also been used beyond the standard $$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$3-divergence class. A 2026 work introduces a quantum relative-$$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$4-entropy

$$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$5

with the convention $$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$6 whenever $$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$7 (Roy et al., 8 Apr 2026).

For the NS-type distributions

$$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$8

the paper proves the exact correspondence

$$\rho=\sum_i \lambda_i |v_i\rangle\langle v_i|,\qquad \sigma=\sum_j \mu_j |w_j\rangle\langle w_j|,$$9

where

$(i,j)$0

The work emphasizes that this divergence lies outside the quantum $(i,j)$1-divergence class, yet still admits an exact classical reduction through NS-type distributions (Roy et al., 8 Apr 2026).

Within that framework, the overlap matrix $(i,j)$2 is interpreted as encoding the relative geometry of the eigenbases. The divergence depends on the eigenvalues and on $(i,j)$3, is invariant under simultaneous unitary conjugation, and is additive under tensor products because the corresponding NS distributions factorize: $(i,j)$4 This suggests that the NS formalism is not limited to the traditional Petz-type setting, but can also expose geometric content in divergences whose classical analogues are not Csiszár $(i,j)$5-divergences (Roy et al., 8 Apr 2026).

That same work derives generalized convexity statements on commuting subclasses using multiplicative mixtures

$(i,j)$6

and obtains a generalized convexity result for Petz–Rényi divergence for $(i,j)$7, complementing the known convexity for $(i,j)$8 (Roy et al., 8 Apr 2026).

6. Operational refinements in binary and asymmetric testing

In binary Bayesian testing, the NS map compares the optimal quantum Bayes error with the optimal classical Bayes error of the associated NS pair. For priors $(i,j)$9, let $\sigma$00 and $\sigma$01. The optimal quantum Bayes error is

$\sigma$02

while the optimal classical Bayes error for the NS distributions of $\sigma$03 is

$\sigma$04

A dimension-free result shows

$\sigma$05

for all separable Hilbert spaces and arbitrary priors. Nussbaum–Szkoła had established the lower bound $\sigma$06; the matching upper bound is the new contribution in (Cheng et al., 4 Jun 2026).

The same paper introduces the Petz–Nussbaum–Szkoła trace harmonic mean

$\sigma$07

where $\sigma$08 denotes the Kubo–Ando weighted harmonic mean. It proves

$\sigma$09

and

$\sigma$10

At $\sigma$11, this yields a direct sandwich for the quantum Bayes error, and combining it with the scalar comparison between $\sigma$12 and $\sigma$13 gives the factor-of-two theorem (Cheng et al., 4 Jun 2026).

For asymmetric hypothesis testing, a converse bound based on the NS mapping takes the form

$\sigma$14

Here $\sigma$15 is the optimal type-II error under a type-I constraint $\sigma$16, and $\sigma$17 is the classical Neyman–Pearson trade-off for the NS pair. This single one-shot inequality yields unified converses in the small-, moderate-, and large-deviation regimes by choosing $\sigma$18 appropriately and importing classical asymptotic results for $\sigma$19 (Lizarribar-Carrillo et al., 20 Jan 2026).

The transferred asymptotic forms include: $\sigma$20 in the fixed-$\sigma$21 regime when $\sigma$22,

$\sigma$23

for $\sigma$24, and the Hoeffding-type exponent bound

$\sigma$25

when $\sigma$26 (Lizarribar-Carrillo et al., 20 Jan 2026).

These refinements strengthen the operational interpretation of NS distributions. They are not merely a vehicle for asymptotic exponents; they can bound the exact one-shot error trade-off and provide accurate finite-blocklength approximations (Lizarribar-Carrillo et al., 20 Jan 2026, Cheng et al., 4 Jun 2026).

7. Special cases, caveats, and conceptual position

Several simplifying scenarios clarify the NS construction.

If $\sigma$27 and $\sigma$28 commute, one can choose a common eigenbasis. Then

$\sigma$29

and the NS distributions reduce to the ordinary classical eigenvalue distributions. In that case the Chernoff distance is classical,

$\sigma$30

and classical finite-blocklength performance directly characterizes the quantum problem (Li, 2015, Androulakis et al., 2023).

For pure states $\sigma$31 and $\sigma$32,

$\sigma$33

for all $\sigma$34, hence

$\sigma$35

The corresponding NS distributions have a single nonzero entry, and the classical Chernoff information equals the quantum one (Li, 2015).

When supports are pairwise disjoint, overlap terms vanish and the NS construction becomes trivial in the sense that the relevant classical overlaps are zero. In such cases Gram–Schmidt-based constructions directly yield optimal measurements in multiple-state testing (Li, 2015).

At the same time, several caveats recur across the literature.

First, NS distributions are not obtained by measuring the two states with a single POVM. Measured $\sigma$36-divergences are produced by a fixed measurement and typically lower bound the corresponding quantum divergence. NS distributions are different: they depend simultaneously on both states through their eigenbases and serve as an analytical construct yielding exact equalities for Petz-type $\sigma$37-divergences (Androulakis et al., 2023).

Second, degeneracies matter at the level of the distributions themselves. Different orthonormal bases within degenerate eigenspaces can lead to different NS pairs, even though the divergence values relevant to the main equalities remain basis-independent (Androulakis et al., 2023). In asymmetric converse bounds, the validity of the inequalities does not depend on the basis choice, though numerical tightness may (Lizarribar-Carrillo et al., 20 Jan 2026).

Third, the exact NS equality does not extend unchanged to all quantum divergences. For example, the equality $\sigma$38 generally does not hold for sandwiched Rényi divergence; the exact correspondence applies to Petz-type Rényi divergences and to the Petz-type $\sigma$39-divergence class (Androulakis et al., 2023).

Finally, the domain of the theory has widened substantially. What began as a finite-dimensional tool for Chernoff-type quantum hypothesis testing now applies in separable Hilbert spaces for both hypothesis-testing exponents and $\sigma$40-divergence equalities, and in semifinite von Neumann algebras through a joint spectral representation of left and right multiplication operators (Cheng et al., 4 Jun 2026, Anastasiadis et al., 21 Apr 2026). This suggests a durable structural role for Nussbaum–Szkoła distributions: they provide a precise classical model of the joint spectral mismatch between two quantum states, and in that model a wide range of quantum distinguishability problems become classical.