Geometric Quantum Discord (GQD)

• Geometric Quantum Discord (GQD) is a metric-based measure that quantifies non-classical correlations in bipartite and multipartite quantum states by minimizing the distance to zero-discord states.
• GQD employs diverse metrics—including Hilbert–Schmidt, trace, and Bures distances—to capture distinct operational, contractivity, and computational features, enhancing analytical tractability.
• Its applications span quantum state discrimination, metrology, and multipartite systems, offering insights into quantum correlations and dynamics under various noise channels.

Geometric quantum discord (GQD) is a distance-based quantifier of non-classical correlations in a bipartite or multipartite quantum state, defined by measuring how far the state lies from an appropriate set of zero-discord, or classical, states. In its original finite-dimensional form, GQD replaces the optimization over local measurements appearing in entropic quantum discord by an optimization over distances in state space, typically to the set of classical-quantum states (Hassan et al., 2010). Subsequent work showed that this geometric idea is not tied to a unique metric: Hilbert–Schmidt, trace, and Bures constructions lead to distinct discord functionals with different monotonicity, operational, and computational properties (Paula et al., 2013, Spehner et al., 2013).

1. Classical-quantum geometry and the basic definition

For a bipartite state $\rho$, the zero-discord states with respect to measurements on subsystem $A$ are classical-quantum states of the form

$$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$

or, equivalently,

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$

with $\{\Pi_k^a\}$ orthogonal projectors on the measured subsystem and $\rho_k^b$ arbitrary states of the unmeasured subsystem (Paula et al., 2013, Hassan et al., 2010). In the Bures formulation these are denoted $A$-classical states and written as

$$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$

emphasizing the asymmetry inherited from the choice of measured subsystem (Spehner et al., 2013).

The Hilbert–Schmidt version introduced in the early literature is

$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$

so that GQD is the minimum squared Hilbert–Schmidt distance to the zero-discord set (Paula et al., 2013). Equivalent measurement-based forms were also derived: $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$ which makes explicit that GQD can be understood as the minimal disturbance induced by an optimal local projective measurement (Doustimotlagh, 2014).

This geometric reformulation was motivated by the computational difficulty of entropic discord. In the formulation $A$0, the classical correlation term requires an optimization over local measurements, and exact formulas are available only for restricted state families. GQD replaces that information-theoretic optimization by a distance minimization to the nearest zero-discord state, which is often more tractable analytically and numerically (Doustimotlagh, 2014).

2. Hilbert–Schmidt GQD and analytically solvable families

For two-qubit states, the Hilbert–Schmidt GQD admits the well-known closed form

$A$1

where $A$2, $A$3 with $A$4, and $A$5 is the largest eigenvalue of $A$6 (Qiang et al., 2014, Hassan et al., 2012). In Bloch-operator language for a general $A$7 bipartite system, the problem can be written as an optimization over the coefficient matrix $A$8,

$A$9

which underlies several general bounds and exact results (Hassan et al., 2010).

For Bell-diagonal two-qubit states,

$$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$0

the Hilbert–Schmidt discord simplifies to

$$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$1

and the physically allowed region is the tetrahedron in $$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$2-space, with the zero-discord Bell-diagonal states lying on the coordinate axes (Yao et al., 2013). This yields a geometric picture in which level sets of constant discord are three intersecting cylinders truncated by the tetrahedron (Yao et al., 2013).

Several later works pushed exact Hilbert–Schmidt formulas beyond the original two-qubit setting. A rigorous lower bound stronger than the earlier Luo–Fu bound was derived for arbitrary finite-dimensional bipartite states in terms of the coherence vector $$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$3, correlation matrix $$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$4, and the eigenvalues of

$$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$5

with the bound becoming exact for $$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$6 systems measured on the qubit side (Hassan et al., 2010). More recently, an explicit exact analytical value for the GQD of an arbitrary two-qudit state was obtained via the Bloch vector of the measured subsystem, the correlation matrix, and the spectrum of

$$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$7

with the important conclusion that a previously known lower bound is in fact the exact value for every two-qudit state (Loubenets et al., 2024).

Further exact or closed formulas were found for structured state families. For two-qubit $$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$8 states, both geometric discord and geometric global quantum discord were derived analytically, with the latter satisfying $$\chi=\sum_k p_k\,|k\rangle\langle k|\otimes \rho_k,$$9 for all $\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$0 states (Qiang et al., 2014). For the multiqubit family

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$1

the multipartite geometric discord takes the exact form

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$2

showing an explicit $\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$3 scaling for fixed correlation coefficients (Zhu et al., 2021).

3. Metric choices, consistency, and the criticism of the Hilbert–Schmidt geometry

A central development in the subject was the realization that the Hilbert–Schmidt version is not a good measure of quantum correlations, because it may increase or change under local reversible operations on the unmeasured subsystem (Paula et al., 2013). For the local ancilla map

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$4

one has

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$5

so simply attaching or removing a mixed local ancilla on the unmeasured side changes the discord value by the ancilla purity factor (Paula et al., 2013). This pathology was traced to the noncontractivity of the Hilbert–Schmidt norm under trace-preserving completely positive maps (Paula et al., 2013).

This motivated alternative geometric discords based on better-behaved distances. In the Schatten-$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$6 family,

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$7

the multiplicativity

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$8

implies

$\rho_c=\sum_k p_k\,\Pi_k^a\otimes \rho_k^b,$9

Because $\{\Pi_k^a\}$0, the trace norm is the only Schatten norm for which the value is invariant under this class of local reversible operations, making the $\{\Pi_k^a\}$1-norm geometric discord the only consistent Schatten-$\{\Pi_k^a\}$2 version in that framework (Paula et al., 2013). The trace-distance discord is therefore defined as

$\{\Pi_k^a\}$3

and satisfies the contractivity property

$\{\Pi_k^a\}$4

for any trace-preserving local operation $\{\Pi_k^a\}$5 on the unmeasured subsystem (Paula et al., 2013).

For Bell-diagonal states the trace-distance discord has the simple formula

$\{\Pi_k^a\}$6

the intermediate value among the absolute correlation coefficients, and coincides on that family with the negativity of quantumness (Paula et al., 2013). A hierarchy established in this setting is

$\{\Pi_k^a\}$7

where $\{\Pi_k^a\}$8 is entropic discord and $\{\Pi_k^a\}$9 is entanglement negativity (Paula et al., 2013).

The Bures-distance construction pursues the same consistency goal from a different direction. The Bures distance

$\rho_k^b$0

is contractive under completely positive trace-preserving maps, Riemannian, and linked to the quantum Fisher information (Spehner et al., 2013). The corresponding Bures geometric discord is

$\rho_k^b$1

and vanishes on exactly the same $\rho_k^b$2-classical states as standard discord (Spehner et al., 2013). This suggests that “geometric quantum discord” is not a single quantity but a family of metric-dependent correlation measures sharing the same zero set while differing in contractivity, ordering, and operational meaning.

4. Pure states, higher dimensions, multipartite systems, and Gaussian extensions

For pure states, several GQD variants collapse to especially transparent forms. In the Bures case, if

$\rho_k^b$3

is in Schmidt form and $\rho_k^b$4 is the largest Schmidt coefficient, then

$\rho_k^b$5

so the Bures-GQD coincides with the geometric measure of entanglement for pure states (Spehner et al., 2013). In the Hilbert–Schmidt framework used in a recent holographic analysis, the pure-state formula becomes

$\rho_k^b$6

with $\rho_k^b$7 the Schmidt coefficients and $\rho_k^b$8 the second Rényi entropy (Banerjee et al., 2023).

Multipartite generalizations proceed in more than one way. One line defines discord with respect to a specific subsystem in an $\rho_k^b$9-partite system and yields a generic tensor formula

$A$0

where $A$1 is the correlation tensor in local operator bases (Hassan et al., 2012). Another line defines total quantum correlations by sequential optimal local measurements on all parties, leading to

$A$2

with explicit computable formulas for $A$3-qubit states (Hassan et al., 2012). A different multipartite extension, geometric global quantum discord, reduces for two-qubit $A$4 states to a closed formula and coincides with total quantum correlations in the sense of Hassan and Joag (Qiang et al., 2014).

Continuous-variable Gaussian systems require further restriction because Gaussian local measurements preserve Gaussianity. The Gaussian geometric discord is defined as

$A$5

with the optimization restricted to Gaussian POVMs on one subsystem (Adesso et al., 2011). For two-mode squeezed thermal states, characterized by $A$6 in standard-form covariance matrices, the optimal Gaussian POVM is a noisy heterodyne measurement and the discord becomes

$A$7

(Adesso et al., 2011). A crucial result in this setting is that, without an energy bound, there is no universal positive lower bound on Gaussian geometric discord at fixed entropic discord; one can have arbitrarily small geometric discord and arbitrarily large entropic discord (Adesso et al., 2011). This sharpens the distinction between geometric and entropic quantifiers in infinite-dimensional systems.

5. Operational interpretations and applications

The Bures formulation has a particularly direct operational meaning. For mixed states, the maximal fidelity $A$8 to the set of $A$9-classical states is equal to the optimal success probability of an ambiguous quantum state discrimination task associated with an ensemble $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$0 extracted from $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$1 (Spehner et al., 2013). In this picture, the closer $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$2 is to the classical set, the easier it is to discriminate the associated conditional states. The closest zero-discord states are obtained from the corresponding optimal measurements, and for qubit-measured systems the fidelity can be computed from the eigenvalues of a $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$3 Hermitian matrix $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$4 (Spehner et al., 2013).

This machinery was applied to the DQC1 output state. If the eigenvalues of the unitary $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$5 are $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$6, then

$$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$7

and the geometric discord is largest when the eigenvalues of $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$8 are uniformly distributed on the unit circle modulo symmetry with respect to the origin (Spehner et al., 2013).

Another operational branch uses the modified $$\sigma_{A\text{-cl}}=\sum_i q_i\,|\alpha_i\rangle\langle\alpha_i|\otimes \rho_{B|i},$$9-based geometric discord,

$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$0

which behaves better under local operations on the unmeasured subsystem and, for all $$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$1 systems, is equivalent up to a constant factor to local quantum uncertainty (Cordero et al., 2021). In a lossy N00N-state phase-estimation problem, this quantity satisfies

$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$2

hence

$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$3

so the discord exactly measures the fraction of ideal quantum Fisher information surviving the loss channel (Cordero et al., 2021).

Further applications are dynamical and thermodynamic. In two-spin-$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$4 dimers in multiple quantum NMR, Hilbert–Schmidt GQD, entropic discord, and measurement-induced non-locality can all be evaluated analytically on an $$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$5-state family, and they show very similar dependence on inverse temperature and evolution time (Doustimotlagh, 2014). In holographic settings, the pure-state Hilbert–Schmidt GQD of the thermofield double state becomes

$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$6

so nonzero GQD signals non-factorization of the thermal partition function and is tied in that work to the presence of the Einstein–Rosen bridge (Banerjee et al., 2023). A more recent application to private-state-based quantum key distribution derives a lower bound on the distillable secret key rate in terms of the GQD of $$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$7 and $$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$8, showing that secret key can still be generated even when the relevant shield states are separable or PPT entangled (Jain et al., 5 Sep 2025).

6. Dynamics, ordering relations, and unresolved issues

Because GQD is metric-dependent, its dynamical behavior is likewise geometry-dependent. In structured bosonic reservoirs, both trace-distance and Bures-distance discords can be preserved well, and may even be improved and generated by a noisy common reservoir (Hu et al., 2014). For Bell-diagonal states under independent phase-flip channels, Hilbert–Schmidt GQD evolves as

$$D_G(\rho)=\min_{\chi\in\Omega_0}\|\rho-\chi\|_2^2, \qquad \|X\|_2=\sqrt{\operatorname{tr}(X^\dagger X)},$$9

with $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$0, $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$1, and $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$2, producing a nonanalytic sudden-change point when the maximizing coefficient switches branches (Yao et al., 2013). In that same Bell-diagonal setting, if one requires the Hilbert–Schmidt GQD to remain constant for a finite period under the phase-flip channel, the initial state must be separable (Yao et al., 2013).

Multipartite open-system studies reach a different conclusion about robustness. For three-qubit and six-qubit Werner–GHZ states under amplitude damping, phase damping, depolarizing, and flipping channels, multipartite geometric discord is more fragile than global entropic discord, and depolarizing noise is the most destructive channel among those considered (Ramzan, 2012). In another multiqubit family, local phase noise acting on a single qubit produces sudden changes of multipartite geometric discord at

$$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$3

when $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$4, reflecting a switch in the dominant coefficient in the optimization formula (Zhu et al., 2021).

Ordering relations between geometric and entropic measures are only partly stable. For Bell-diagonal states one has the hierarchy $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$5, but monotonic ordering between $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$6, $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$7, and entropic discord is preserved only on highly symmetric families such as $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$8-symmetric states with $$D_A^g(\rho)=\min_{\{\Pi_k^A\}} \left\| \rho-\sum_k(\Pi_k^A\otimes I)\rho(\Pi_k^A\otimes I) \right\|^2,$$9; it breaks down on less symmetric $A$00-symmetric families with $A$01 (Paula et al., 2013). In Gaussian systems, fixed entropic discord does not control Hilbert–Schmidt geometric discord without an energy constraint (Adesso et al., 2011). This suggests that discord orderings are geometry-sensitive rather than universal.

Two persistent misconceptions are therefore corrected by the literature. First, GQD is not a single canonical quantity; Hilbert–Schmidt, trace, Bures, and $A$02-based versions define distinct measures with distinct mathematical properties (Paula et al., 2013, Spehner et al., 2013, Cordero et al., 2021). Second, computational simplicity does not by itself justify a metric choice. The Hilbert–Schmidt geometry delivers many closed formulas and useful visualizations, but its noncontractivity under local operations on the unmeasured subsystem remains a decisive limitation for interpreting it as a bona fide measure of quantum correlations (Paula et al., 2013). The subsequent development of trace-distance and Bures-distance discords, together with exact formulas for arbitrary two-qudit states and operational links to discrimination, metrology, and non-factorization, has turned GQD from a convenient surrogate for entropic discord into a technically diverse framework for quantifying nonclassical correlations across finite-dimensional, Gaussian, and multipartite quantum systems (Loubenets et al., 2024, Spehner et al., 2013, Banerjee et al., 2023).