---
title: Quantum Hellinger Distance
url: https://www.emergentmind.com/topics/quantum-hellinger-distance
type: topic
---

# Quantum Hellinger Distance

Searching arXiv for recent and foundational papers on quantum Hellinger distance to ground the article in the cited literature.
arXiv search query: "quantum Hellinger distance"
Quantum Hellinger distance is a square-root-based distance on quantum states, typically defined for density matrices \(\rho\) and \(\sigma\) by
\[
d_{\mathrm{Hel}}(\rho,\sigma)=\|\sqrt{\rho}-\sqrt{\sigma}\|_2
=\left(2-2\,\mathrm{tr}\sqrt{\rho}\sqrt{\sigma}\right)^{1/2},
\]
with the overlap term \(\mathrm{tr}\sqrt{\rho}\sqrt{\sigma}\) often called the quantum affinity [1611.03449]. In the commuting case it reduces to the classical Hellinger distance between probability distributions, while in the noncommutative case it defines a contractive, Riemannian metric that has become a basis for geometric quantifiers of coherence, discord, non-classical correlation, and measurement-induced nonlocality [1611.03449, 1806.10814]. A parallel operator-theoretic line of work treats the usual geometric-mean-based expression as one member of a broader family of generalized quantum Hellinger divergences built from Kubo–Ando means, thereby linking the subject to maximal quantum \(f\)-divergences, Bregman geometry, and barycenter problems [1903.10455].

## 1. Definitions, affinity, and normalization conventions

For density matrices on a finite-dimensional Hilbert space, the survey literature defines the quantum Hellinger distance by
\[
d_{\mathrm{Hel}}(\rho,\sigma)=\left(2-2\,\mathrm{tr}\sqrt{\rho}\sqrt{\sigma}\right)^{1/2}
\]
and equivalently as the Hilbert–Schmidt norm of the difference of square roots [1611.03449]. Closely related papers instead work directly with the squared form
\[
D_H(\rho,\sigma)=\mathrm{Tr}\,(\sqrt{\rho}-\sqrt{\sigma})^2,
\]
or with a rescaled version such as
\[
D_H(\rho_1,\rho_2)=\sqrt{1-\mathrm{Tr}\big(\sqrt{\rho_2}\sqrt{\rho_1}\big)}.
\]
These formulas differ by constants and square roots, so the literature is not uniform in normalization [1806.10814, 1306.3248].

The common structural quantity is the affinity
\[
A(\rho,\sigma)=\mathrm{Tr}\big(\sqrt{\rho}\sqrt{\sigma}\big),
\]
so that \(d_H^2\) is expressed as \(2-2A(\rho,\sigma)\) in one widespread convention [2409.14560]. One source states the bound
\[
0\le D_H[\rho_1,\rho_2]\le \sqrt{2},
\]
while another notes that some authors include an extra factor \(1/\sqrt{2}\), rescaling the range to \([0,1]\) [1107.1732, 2409.14560]. The same square-root structure underlies the classical reduction: for commuting states \(\rho=\sum_k p_k|k\rangle\langle k|\) and \(\sigma=\sum_k q_k|k\rangle\langle k|\), the quantum formula collapses to
\[
\left(\sum_k(\sqrt{p_k}-\sqrt{q_k})^2\right)^{1/2},
\]
the ordinary Hellinger distance between probability vectors [1611.03449].

This dependence on \(\sqrt{\rho}\) rather than \(\rho\) is the defining technical feature of the subject. It gives a noncommutative overlap simpler than Uhlmann fidelity, yet still close enough to fidelity-based geometry to support contractivity and several resource-theoretic constructions [1408.4477, 1611.03449].

## 2. Metric structure, convexity, and information geometry

The quantum Hellinger distance is presented as a genuine metric: it is nonnegative, symmetric, vanishes exactly when \(\rho=\sigma\), and satisfies the triangle inequality [1611.03449]. It is also contractive under quantum operations,
\[
d_{\mathrm{Hel}}(\mathcal M(\rho),\mathcal M(\sigma))\le d_{\mathrm{Hel}}(\rho,\sigma),
\]
and its square is jointly convex [1611.03449]. These properties place it among the standard contractive distances used in quantum information geometry.

The same survey identifies it as a Riemannian distance and gives the induced metric tensor
\[
(g_{\mathrm{Hel}})_\rho(O,O)=\sum_{k,l}\frac{|\langle k|O|l\rangle|^2}{(\sqrt{p_k}+\sqrt{p_l})^2},
\]
for \(\rho=\sum_k p_k|k\rangle\langle k|\) and traceless Hermitian tangent vector \(O\) [1611.03449]. In this sense the distance is not merely algebraic; it induces a differential geometry on state space. The same source contrasts it with the Bures metric and records the inequality
\[
d_{\mathrm{Bu}}(\rho,\sigma)\le d_{\mathrm{Hel}}(\rho,\sigma)\le \sqrt{2}\,d_{\mathrm{Bu}}(\rho,\sigma),
\]
together with the fact that Bures coincides with the Fubini–Study distance on pure states whereas Hellinger does not [1611.03449].

A central operational interpretation comes through quantum uncertainty. For unitary dynamics generated by \(H\), the infinitesimal Hellinger speed is linked to the Wigner–Yanase skew information,
\[
I_{\mathrm{skew}}(\rho,H)=-\frac12\mathrm{tr}([\sqrt{\rho},H]^2)
=\frac12(g_{\mathrm{Hel}})_\rho(-i[H,\rho],-i[H,\rho]),
\]
which underlies the connection between Hellinger geometry and Local Quantum Uncertainty [1611.03449]. This distinguishes it from the Bures metric, whose infinitesimal form is tied instead to quantum Fisher information [1611.03449].

The square-root overlap also interacts naturally with relative-entropy-type quantities. The survey records
\[
d_{\mathrm{Hel}}^2(\rho,\sigma)=2-2\exp\!\left[-\tfrac12 S_{1/2,1}(\rho\|\sigma)\right],
\]
placing the metric within the \(\alpha\)-\(z\) relative entropy framework [1611.03449]. This suggests that quantum Hellinger distance occupies an intermediate position between fidelity geometry and divergence geometry, a theme made explicit in the generalized operator-mean formulation.

## 3. Operator means, generalized divergences, and barycenters

A major extension replaces the geometric mean in the standard formula by an arbitrary Kubo–Ando mean. For positive definite operators \(A,B\in B(\mathcal H)_{++}\), the generalized quantum Hellinger divergence is defined by
\[
D_\mu(A,B)=\operatorname{Tr}\Big((1-c(\mu))A+c(\mu)B-A\,f_\mu\,B\Big),
\]
or equivalently
\[
\phi(A,B)=\operatorname{Tr}\big((1-c)A+cB-A\sigma B\big),
\]
where \(c\in(0,1)\) is the weight of the Kubo–Ando mean \(\sigma\) [1903.10455]. The usual Hellinger-type quantity studied by Bhatia, Gaubert, and Jain is recovered by taking \(\sigma=\#\), \(c=\tfrac12\), and the arcsine measure, giving
\[
\operatorname{Tr}\left(\frac12A+\frac12B-A\#B\right)
\]
as the special geometric-mean case [1903.10455].

The same paper shows that these generalized divergences are maximal quantum \(f\)-divergences. Writing
\[
g_\mu(x)=(1-c(\mu))+c(\mu)x-f_\mu(x),
\]
one has
\[
D_\mu(A,B)=\operatorname{Tr}\Big(A\,g_\mu(A^{-1/2}BA^{-1/2})\Big),
\]
and because \(g_\mu\) is operator convex, the divergence is jointly convex and satisfies the data processing inequality for completely positive trace-preserving maps [1903.10455]. It is also related to operator-valued Bregman divergences, which yields nonnegativity, definiteness, vanishing first derivative at the diagonal, and positive second derivative at the diagonal [1903.10455].

This framework leads naturally to barycenters. For positive definite operators \(A_1,\dots,A_m\) with weights \(w_1,\dots,w_m\), the barycenter is the unique minimizer of
\[
X\mapsto \sum_{j=1}^m w_j D_\mu(A_j,X).
\]
Its characterization is the matrix equation
\[
c(\mu)\,I = \sum_{j=1}^m w_j \int_{[0,1]} \big((1-\lambda)A_j^{-1}X+\lambda I\big)^{-1} \,d\mu(\lambda),
\]
which has a unique positive definite solution [1903.10455].

An important correction follows. The claim that the barycenter for the geometric-mean-based Hellinger divergence is the weighted multivariate \(1/2\)-power mean is true in the commuting case but false in general [1903.10455]. A numerical \(2\times2\) counterexample shows that the minimizer of the Hellinger divergence need not coincide with the weighted \(1/2\)-power mean when the matrices do not commute [1903.10455]. This is a precise noncommutative obstruction rather than a matter of notation.

## 4. Quantum coherence, discord, and measurement-induced nonlocality

Quantum Hellinger distance has been used extensively to define distance-to-free-set resource measures. For coherence in a fixed computational basis, one paper defines
\[
C_H(\rho)=\min_{\delta\in\mathcal I}D_H(\rho,\delta),
\qquad
D_H(\rho,\delta)=\mathrm{Tr}\,(\sqrt{\rho}-\sqrt{\delta})^2,
\]
and derives the closed form
\[
C_H(\rho)=2\left(1-\sqrt{\sum_{k=0}^{d-1}\langle k|\sqrt{\rho}|k\rangle^2}\right),
\]
with the optimal incoherent state
\[
\delta_0=\sum_k \frac{\langle k|\sqrt{\rho}|k\rangle^2}{\sum_{k'}\langle k'|\sqrt{\rho}|k'\rangle^2}|k\rangle\langle k|
\]
[1806.10814]. The same paper proves faithfulness, convexity, and strong monotonicity under incoherent selective operations, so \(C_H\) satisfies the full set of standard coherence axioms [1806.10814].

For non-classical correlation, the same work defines
\[
D(\rho_{AB})=\min_{\{\Pi_A^k\}}\sum_k D_H\!\left(\rho_{AB},(\Pi_A^k\otimes I_B)\rho_{AB}(\Pi_A^k\otimes I_B)\right),
\]
proves local-unitary invariance and contractivity under CPTP maps on the unmeasured subsystem, and shows
\[
D(\rho_{AB})=0 \iff \rho_{AB}=\sum_k \lambda_k |k\rangle\langle k|\otimes \rho_k,
\]
so the measure vanishes exactly on classical-quantum states [1806.10814]. For qubit–qudit states an analytic formula is obtained,
\[
D(\rho_{AB})=1-\lambda_{\max},
\]
with \(\lambda_{\max}\) the largest eigenvalue of a matrix built from \(\sqrt{\rho_{AB}}\) and Pauli operators [1806.10814].

A broader geometric program uses Hellinger distance to define geometric discord, measurement-induced geometric discord, and discord of response [1510.06995]. In that framework,
\[
D_H^G(\rho)=2-2\max_{\sigma_A\in\mathcal C_A}\mathrm{tr}\sqrt{\rho}\sqrt{\sigma_A},
\]
and for pure states with Schmidt coefficients \(\mu_i\),
\[
D_H^G(|\Psi\rangle)=2-2\,K(|\Psi\rangle)^{-1/2},
\qquad
K(|\Psi\rangle)=\left(\sum_i\mu_i^2\right)^{-1}
\]
[1611.03449]. For qubit–qudit systems, both the Hellinger geometric discord and the Hellinger discord of response are fully computable, and the latter is connected to Local Quantum Uncertainty [1510.06995, 1611.03449].

A related line quantifies measurement-induced nonlocality by
\[
N_H(\rho)=1-\min_{\Pi^a}\mathrm{tr}\big[\sqrt{\rho}\,\Pi^a(\sqrt{\rho})\big],
\]
where the optimization is over locally invariant von Neumann measurements on subsystem \(a\) [2007.08126]. This Hellinger-distance MIN is presented as a bona fide measure of nonlocal correlation, is resistant to the local ancilla problem, and coincides with skew-MIN in the formulation used there [2007.08126]. For a bipartite pure state with Schmidt coefficients \(s_i\),
\[
N_H(|\psi\rangle\langle\psi|)=1-\sum_i s_i^2,
\]
and closed formulas are also obtained for general \(2\otimes n\) mixed states [2007.08126].

The multipartite extension appears in a Hellinger-based generalized geometric discord defined by
\[
D^H=\frac12\min_\sigma \|\sqrt{\rho}-\sqrt{\sigma}\|^2
=1-\max_\sigma \mathrm{Tr}\!\left[\sqrt{\rho}\sqrt{\sigma}\right],
\]
with exact evaluation for bipartite pure states via Schmidt decomposition and direct extensions to symmetric multipartite settings [1412.8098]. For permutation-invariant and translation-invariant states, the nearest classical state is proposed to inherit the same symmetry [1412.8098].

## 5. Gaussian states and open-system distinguishability

In continuous-variable Gaussian settings, Hellinger distance has been used to define a geometric Gaussian discord by minimizing over product Gaussian states [1408.4477]. A central caveat is built into that construction: because the zero-discord Gaussian states are precisely the product Gaussian states, the resulting Hellinger quantity measures all intermode correlations, not only the genuinely quantum part, and is therefore an upper bound for the geometric discord in the unrestricted sense [1408.4477]. The corresponding affinity for \(n\)-mode Gaussian states is written explicitly in terms of the covariance matrices of the square-root states, and in the two-mode case the closest Gaussian product state can be determined exactly [1408.4477].

For symmetric two-mode squeezed thermal states,
\[
D_H(\hat\rho_{ST})=[\tanh(r)]^2,
\]
while for mode-mixed thermal states analogous closed forms are obtained [1408.4477]. The same work argues that, in these Gaussian families, the Hellinger-based quantity behaves more like a geometric analogue of quantum mutual information than a strict measure of quantum discord [1408.4477]. This is a recurrent source of confusion in the Gaussian literature, and the distinction is explicit in the source.

In open-system dynamics, the quantum Hellinger distance has also been compared with trace distance, Bures distance, and Jensen–Shannon divergence as a witness of initial system–environment correlations [1107.1732, 1306.3248]. The main conclusion is metric dependence. In an infinite-environment dephasing model, the Hellinger distance can decrease at early times, reach a minimum, and then rise again, but it remains below its initial value; unlike the trace distance, it does not show the same correlation-induced growth above the initial distinguishability [1107.1732]. In the bosonic dephasing model of a later comparative study,
\[
\Delta D_H(\lambda,t)\le 0 \quad \forall \lambda,t,
\]
so the Hellinger distance does not witness initial correlations there [1306.3248]. In a spin-bath model it can increase, but much less often than the trace distance [1306.3248].

These results rule out a common overgeneralization: contractivity and geometric regularity do not imply that a distance is an effective dynamical witness of initial correlations in reduced open-system dynamics [1107.1732, 1306.3248].

## 6. Random states and matrix-analytic extensions

The quantum Hellinger distance has also been studied statistically for random density matrices. For two states \(\rho_1,\rho_2\), one paper defines
\[
d_H(\rho_1,\rho_2)=\sqrt{\operatorname{tr}\big(\sqrt{\rho_1}-\sqrt{\rho_2}\big)^2}
=\sqrt{2-2A(\rho_1,\rho_2)},
\]
with \(A(\rho_1,\rho_2)=\operatorname{tr}(\sqrt{\rho_1}\sqrt{\rho_2})\), and derives exact mean and variance of the squared distance \(D_H=d_H^2\) when one or both states are random [2409.14560]. The random ensembles considered are the Hilbert–Schmidt and Bures–Hall ensembles, and the analysis reduces the first two cumulants of \(D_H\) to the first two moments of the affinity [2409.14560]. Matching those cumulants yields a gamma-distribution approximation for the law of \(D_H\), which is reported to agree well with Monte Carlo simulations [2409.14560]. The same work states that, for matched parameters, the average distance is smallest for two Hilbert–Schmidt states, largest for two Bures–Hall states, and intermediate for mixed Hilbert–Schmidt/Bures–Hall pairs [2409.14560].

A related matrix-analysis literature studies “matrix versions of the Hellinger distance” on the cone of positive definite matrices [1901.01378]. The literal square-root analogue is
\[
d_1(A,B)=\|A^{1/2}-B^{1/2}\|_2,
\]
while the Bures–Wasserstein metric is
\[
d_2(A,B)=\Big[\mathrm{tr}(A+B)-2\mathrm{tr}(A^{1/2}BA^{1/2})^{1/2}\Big]^{1/2}
\]
[1901.01378]. Two further constructions replace the overlap term by the Pusz–Woronowicz geometric mean \(A\#B\) and the log-Euclidean mean \(\exp((\log A+\log B)/2)\), giving quantities whose squares are divergences rather than metrics [1901.01378]. The squared forms
\[
\Phi_3(A,B)=\mathrm{tr}(A+B)-2\mathrm{tr}(A\#B),
\qquad
\Phi_4(A,B)=\mathrm{tr}(A+B)-2\mathrm{tr}\exp\!\big(\tfrac{\log A+\log B}{2}\big)
\]
are shown to be jointly convex and strictly convex in each variable separately, and their barycenters are characterized by nonlinear matrix equations [1901.01378].

This broader matrix setting does not identify a single noncommutative Hellinger geometry. Instead, it exhibits a family of square-root- and mean-based noncommutative distances and divergences, with the standard quantum Hellinger distance corresponding to the Hilbert–Schmidt norm of square-root differences and the generalized Kubo–Ando framework giving a systematic extension beyond that basic case [1901.01378, 1903.10455].

Source: https://www.emergentmind.com/topics/quantum-hellinger-distance