Fisher Metric: Definition and Applications

• The Fisher metric is a Riemannian metric on statistical models quantifying distinguishability between neighboring probability distributions, defined by specific integrals of score functions.
• It induces geodesics, distances, and curvature, essential for statistical inference and optimal transportation, and can be used in various fields such as quantum theory and Riemannian geometry.
• Applications of the Fisher metric include functional data registration, quantum state analysis, and Ricci flow geometry.

The Fisher metric, or Fisher–Rao metric, is a Riemannian metric on a statistical model that quantifies the local distinguishability of neighboring probability distributions. For a smooth family of normalized densities $p(x;\theta)$, it is defined by

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$

or equivalently,

$$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$

Under the usual regularity conditions it is also the negative expected Hessian of the log-likelihood,

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$

The metric induces geodesics, distances, curvature, and natural notions of statistical speed and information loss. Its square-root representation identifies probability densities with a submanifold of an $L^2$ sphere, while its monotonicity under probabilistic maps underlies its characterization by Čencov-type theorems. In quantum theory, the corresponding Hermitian tensor combines a Riemannian Fisher or Fubini–Study metric with a symplectic component; in geometric analysis, analogous constructions measure the distinguishability of heat-kernel measures.

1. Definition and fundamental structures

Let $\mathcal P$ be a statistical model consisting of probability distributions $p_\theta$ indexed by $\theta=(\theta^1,\ldots,\theta^d)$. A tangent vector corresponds to an infinitesimal perturbation of the density. In a finite sample space, if $X$ and $Y$ have coordinate representations $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$0 and $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$1 with zero total sum, the Fisher metric is

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$2

Writing the score, or $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$3-representation, as

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$4

gives

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$5

Because score functions have zero expectation, this is also their covariance.

The metric is positive semidefinite:

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$6

It is positive definite when the model is identifiable in the relevant tangent directions. In non-identifiable models, distinct tangent vectors can produce the same score and the metric can be degenerate.

For a parameter curve $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$7, the Fisher length is

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$8

where $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$9. The Fisher–Rao distance is the infimum of such lengths over curves joining two distributions. Unlike the Kullback–Leibler divergence, it is symmetric, satisfies the triangle inequality, and separates points under the usual regularity assumptions. Locally,

$$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$0

so relative entropy is locally a squared Fisher distance rather than a distance itself (Miyamoto et al., 2023).

The Fisher metric also has a square-root form:

$$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$1

Consequently, the map

$$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$2

pulls back the $$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$3 metric to the Fisher metric. For categorical distributions, this map embeds the simplex into the positive orthant of a sphere of radius $$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$4. The resulting geometry makes the Fisher–Rao distance a spherical arc length rather than the Euclidean chord associated with the Hellinger distance (Miyamoto et al., 2023).

2. Distances, geodesics, and geometric models

The Fisher distance generally requires solving a geodesic boundary-value problem. In one dimension it is obtained by direct integration:

$$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$5

For categorical distributions $$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$6 and $$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$7,

$$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$8

For a multinomial model with $$g_{ij}(\theta) = \int \frac{\partial_i p(x;\theta)\,\partial_j p(x;\theta)} {p(x;\theta)}\,d\mu(x).$$9 trials, the distance is multiplied by $$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$0:

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$1

Many two-parameter location-scale models have Fisher metrics of the form

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$2

After rescaling coordinates, this becomes a constant multiple of the Poincaré half-plane metric. The corresponding distance is hyperbolic:

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$3

Gaussian, Laplace, Cauchy, Student, logistic, generalized Gaussian, inverse Gaussian, Gumbel, Fréchet, Weibull, Pareto, and power-function models provide examples in which parameter transformations reduce the metric to spherical, hyperbolic, logarithmic, or matrix geometries. Wishart models yield the affine-invariant metric on the cone of positive-definite matrices:

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$4

For product distributions,

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$5

the Fisher metric is block diagonal and the squared distance is additive:

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$6

Closed-form distances are exceptional. Gamma and Beta families generally have nonconstant negative curvature and do not usually admit elementary distance formulas (Miyamoto et al., 2023).

3. Invariance, information loss, and uniqueness

A measurable statistic or Markov kernel maps one statistical model to another. If $$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$7 is a statistic, the processed distribution is $$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$8. The Fisher metric is monotone:

$$g_{ij}(\theta) = -\mathbb E_\theta \left[ \frac{\partial^2\log p(X;\theta)} {\partial\theta^i\partial\theta^j} \right].$$9

Thus data processing cannot increase Fisher information. The proof is based on the $L^2$0 contraction of conditional expectation:

$L^2$1

If the statistic is sufficient, no information is lost and equality holds. In this sense, sufficient statistics act isometrically on the statistical model. The same principle applies to probabilistic mappings between more general, possibly singular models. Diffeological statistical models replace an ordinary parameter manifold with a family of plots and can accommodate boundaries, corners, self-intersections, varying local dimensions, and other singularities. For an almost $L^2$2-integrable diffeological model, tangent scores lie in $L^2$3, and the Fisher quadratic form is

$L^2$4

Probabilistic mappings preserve almost $L^2$5-integrability, and Fisher monotonicity continues to hold (Lê, 2019).

Čencov’s theorem characterizes the Fisher metric, up to a positive multiplicative constant, by invariance under suitable Markov embeddings on finite probability simplices. For general separable metrizable sample spaces, the uniqueness result requires monotonicity together with strong continuity in a mixed topology on the space of measures and their $L^2$6 tangent fields. Under these assumptions,

$L^2$7

for a universal constant $L^2$8. Monotonicity and strong continuity also imply invariance under sufficient statistics (Lê, 2013).

On conditional probability polytopes, the precise normalization depends on the selected structural principle. For a row-stochastic matrix $L^2$9, the three principal metrics are

$\mathcal P$0

the unscaled product metric

$\mathcal P$1

and the marginal-weighted metric

$\mathcal P$2

The first arises from invariance under homogeneous conditional embeddings, the second from exponential-family embeddings of Cartesian products, and the third from embedding conditional distributions into joint distributions with fixed marginal $\mathcal P$3 (Montufar et al., 2014).

4. Dual, generalized, and continuous formulations

The inverse Fisher metric is naturally a metric on the cotangent bundle. For a finite sample space, every cotangent vector is the differential of an expectation:

$\mathcal P$4

Random variables represent cotangent vectors modulo constants because constants have zero differential. Covariance therefore defines the Fisher co-metric:

$\mathcal P$5

In particular,

$\mathcal P$6

The inverse Fisher information matrix is the coordinate form of this cotangent metric. The Cramér–Rao inequality becomes a minimum-norm extension property and can be written

$\mathcal P$7

Under Markov maps, the co-metric is monotone in the dual direction:

$\mathcal P$8

For random variables this is conditional-expectation contraction of variance (Nagaoka, 2023).

The Fisher metric is also the Hessian of relative entropy. For the Kullback–Leibler divergence,

$\mathcal P$9

More generally, if a relative entropy group is generated by a function $p_\theta$0, its Hessian at coincidence is

$p_\theta$1

The Boltzmann–Gibbs and Kaniadakis classes reproduce the ordinary Fisher metric, while the Tsallis class gives

$p_\theta$2

The Levi–Civita connection and Ricci tensor are unchanged by this constant rescaling, whereas scalar curvature is divided by the same factor (Gomez et al., 2018).

For smooth positive densities on a compact manifold, the Fisher–Rao metric is the unique diffeomorphism-invariant metric on probability densities up to a positive constant. On the larger space of all positive densities, diffeomorphism-invariant metrics have the form

$p_\theta$3

The first term governs fixed-volume variations and the second governs total-mass variations. Under the square-root map, these metrics become radial warped metrics on an $p_\theta$4 cone. On probability densities, the radial direction disappears and the geometry is a constant-curvature sphere (Bruveris et al., 2016).

5. Quantum, geometric, and holographic interpretations

For a normalized amplitude

$p_\theta$5

the probability density contains amplitude information while $p_\theta$6 contains phase information. Pulling back the Hermitian tensor of projective Hilbert space gives

$p_\theta$7

Its real part is

$p_\theta$8

where

$p_\theta$9

is the classical Fisher metric. If $\theta=(\theta^1,\ldots,\theta^d)$0, then

$\theta=(\theta^1,\ldots,\theta^d)$1

Thus the classical Fisher metric is the phase-independent sector of the quantum geometric metric, while parameter-dependent phase contributes additional metric and symplectic information (Facchi et al., 2010).

For density matrices, the symmetric logarithmic derivative is defined by

$\theta=(\theta^1,\ldots,\theta^d)$2

and the SLD quantum Fisher tensor is

$\theta=(\theta^1,\ldots,\theta^d)$3

For pure states, $\theta=(\theta^1,\ldots,\theta^d)$4, giving

$\theta=(\theta^1,\ldots,\theta^d)$5

The resulting tensor is proportional to the Fubini–Study metric. For a mixed two-level system with eigenvalues $\theta=(\theta^1,\ldots,\theta^d)$6 and $\theta=(\theta^1,\ldots,\theta^d)$7, the fixed-eigenvalue orbit has

$\theta=(\theta^1,\ldots,\theta^d)$8

The orbit collapses at the maximally mixed state $\theta=(\theta^1,\ldots,\theta^d)$9, while the pure-state limit recovers the ordinary projective geometry (Ercolessi et al., 2012).

For bipartite quantum systems, the local-unitary quantum Fisher tensor has block form

$X$0

The off-diagonal block is the covariance matrix between observables on the two subsystems. It vanishes for product states and is generally nonzero for entangled states. In spin-network Hilbert spaces, the quantity

$X$1

serves as a geometric entanglement monotone. For fixed-spin maximally entangled link states, it is tied to the spin Casimir and hence to the area carried by the link (Chirco et al., 2017).

Fisher geometry also appears in holographic constructions. If $X$2 is the probability distribution formed from squared Schmidt coefficients,

$X$3

then

$X$4

For the entropy-derived construction,

$X$5

up to normalization. In a $X$6-dimensional conformal field theory, the special logarithmic entropy case yields the constant-time hyperbolic metric

$X$7

whose geodesics are semicircles and whose regulated lengths reproduce the logarithmic form of the Ryu–Takayanagi entropy (Matsueda, 2014).

A related holographic result identifies the relative-entropy Hessian of a vacuum ball with bulk canonical energy. For a ball-shaped region $X$8,

$X$9

where $Y$0 is the canonical energy of the corresponding perturbation of the AdS Rindler wedge. Positivity of relative entropy implies

$Y$1

The first-order version yields the linearized Einstein equations, while the second-order version constrains quadratic gravitational and matter perturbations (Lashkari et al., 2015).

6. Applications beyond finite statistical models

The Fisher metric is used whenever a parameterized family of distributions, fields, or measures is available. In functional data registration, the Fisher–Rao metric on absolutely continuous functions is

$Y$2

The square-root velocity function

$Y$3

maps this variable metric to the standard $Y$4 metric:

$Y$5

Time warps act isometrically by

$Y$6

Quotienting by the warping group separates phase variability from amplitude variability, enables elastic distances, and supports Karcher means and registration through dynamic programming (Srivastava et al., 2011).

For structured optical beams, normalized measured intensity is treated as a probability density:

$Y$7

For transverse displacement,

$Y$8

The metric is sensitive to intensity gradients, nodal structure, local curvature, and spatial oscillations. Hermite–Gaussian modal order, Laguerre–Gaussian vortex charge, and the radial oscillation parameter of finite-energy Bessel–Gauss beams are all reported to increase displacement Fisher information under the specified measurement model (Sumaya-Martinez et al., 29 Dec 2025).

In Ricci-flow geometry, the pointed heat-kernel probability measure

$Y$9

defines a scale-dependent Fisher metric on the source manifold:

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$00

It satisfies

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$01

and is nonincreasing with scale:

$$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$02

The defect $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$03 measures the failure of heat-kernel measures to behave like translated Euclidean Gaussians. Its trace is related to pointed Nash entropy, while equality in a nonzero direction forces Euclidean splitting. Large Fisher eigenvalues generate heat-kernel splitting maps and provide curvature-regularity criteria (Chow et al., 14 Aug 2026).

Fisher geometry can also be placed on Lie groups. In a Souriau-type formulation, the geometric temperature $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$04 belongs to a Lie algebra, the dual heat belongs to its dual, and a cocycle modifies the orbit metric through Lie brackets and central-extension terms. For the explicit $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$05 and $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$06 constructions, the generalized metric is positive on the relevant coadjoint-orbit directions; in $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$07 the stabilizer direction produces a zero eigenvalue on the full Lie algebra but not on the two-dimensional orbit tangent space (Assandje et al., 25 Sep 2025).

The Fisher metric is therefore not restricted to ordinary finite-dimensional likelihood models. It appears as a metric on probability simplices, spaces of densities, singular and diffeological models, conditional distributions, quantum-state manifolds, functional-data quotient spaces, optical intensity distributions, heat-kernel measures, and coadjoint orbits. Across these settings, its defining principle remains the same: infinitesimal geometric length is determined by the $$g_{ij}(\theta) = \int p(x;\theta)\, \frac{\partial\log p(x;\theta)}{\partial\theta^i} \frac{\partial\log p(x;\theta)}{\partial\theta^j}\,d\mu(x),$$08 sensitivity of a normalized statistical object to changes in the parameters or base point.