Two-Parameter Quantum Fisher Matrix

• The two-parameter quantum Fisher information matrix is a 2×2 metric that quantifies estimation precision by capturing both diagonal sensitivities and off-diagonal parameter couplings.
• It establishes a dual role as both a precision matrix and a Riemannian metric, linking fidelity, reparameterization, and quantum geometric structures.
• It addresses challenges like singularity and nonregularity in quantum metrology, guiding optimal measurement strategies and probe state designs.

The two-parameter quantum Fisher information matrix (QFIM) is the $2\times 2$ local information matrix associated with a quantum statistical model $\rho(\theta_1,\theta_2)$. In contemporary quantum metrology it is usually the symmetric-logarithmic-derivative (SLD) QFIM, although in the broader information-geometric formulation quantum Fisher information is a family of monotone metrics indexed by a standard operator monotone function (Petz et al., 2010). In two parameters, the QFIM already exhibits the main structural issues of multiparameter estimation: diagonal sensitivities, off-diagonal parameter coupling, invertibility versus rank deficiency, measurement compatibility, and discontinuities induced by parameter-dependent rank changes (Goldberg et al., 2021).

1. Definition and metric status

For a smooth two-parameter family $\rho(\theta_1,\theta_2)$, the SLDs $L_i$ are defined by

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$

and the SLD-QFIM is

$$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$

The corresponding matrix is

$$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$

and enters the multiparameter quantum Cramér–Rao bound through the inverse matrix when $F$ is nonsingular (Yue et al., 2016).

The broader mathematical setting is not unique. If monotonicity under quantum channels is taken as the defining requirement, then there is a whole family of quantum Fisher informations, each determined by a standard operator monotone function $f$. For a multiparameter model, the corresponding matrix is

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$

with generalized logarithmic derivatives

$\rho(\theta_1,\theta_2)$0

Within this family, the minimal monotone metric, corresponding to $\rho(\theta_1,\theta_2)$1, is the most popular in physics and reproduces the SLD form (Petz et al., 2010).

This dual role—simultaneously a statistical precision matrix and a Riemannian metric tensor—is central. In two parameters, $\rho(\theta_1,\theta_2)$2 and $\rho(\theta_1,\theta_2)$3 are the squared metric lengths of the tangent directions $\rho(\theta_1,\theta_2)$4 and $\rho(\theta_1,\theta_2)$5, while $\rho(\theta_1,\theta_2)$6 is their inner product under the chosen quantum Fisher metric (Petz et al., 2010).

2. Geometry, fidelity, and reparameterization

For mixed states, fidelity can be used as a generating function for the two-parameter QFIM. Writing

$\rho(\theta_1,\theta_2)$7

the mixed-state quantum metric $\rho(\theta_1,\theta_2)$8 satisfies

$\rho(\theta_1,\theta_2)$9

so the nearby-state expansion becomes

$\rho(\theta_1,\theta_2)$0

For a two-parameter family $\rho(\theta_1,\theta_2)$1, this reads

$\rho(\theta_1,\theta_2)$2

(Chen, 7 Nov 2025).

For pure states, the same paper identifies the modulus of the overlap as the generator of the quantum metric and the phase as the generator of the Berry curvature. In the standard normalization used there, $\rho(\theta_1,\theta_2)$3, so the two-parameter QFIM is the symmetric part of the pure-state quantum geometric tensor (Chen, 7 Nov 2025).

Reparameterization is controlled by a Jacobian congruence transformation. If two parameterizations describe the same unitary family, then for a two-parameter change of variables $\rho(\theta_1,\theta_2)$4,

$\rho(\theta_1,\theta_2)$5

The same linear rule holds for the infinitesimal generators in the $\rho(\theta_1,\theta_2)$6 unitary setting, which makes parameter changes geometrically transparent (Shemshadi et al., 2018).

The geometric interpretation becomes subtler in nonregular models. In that regime the QFIM and the Bures metric need not coincide pointwise; the Bures metric can remain continuous even when the QFIM is discontinuous because the density matrix rank changes with the parameters (Goldberg et al., 2021).

3. Invertibility, singular directions, and nonregularity

For a two-parameter Fisher matrix

$\rho(\theta_1,\theta_2)$7

invertibility is decided by

$\rho(\theta_1,\theta_2)$8

If $\rho(\theta_1,\theta_2)$9, the matrix is invertible. If $L_i$0, it is singular, and only one or zero independent parameter combinations are estimable with finite variance (Namkung et al., 2024).

A unified Cramér–Rao formulation replaces the inverse by the Moore–Penrose pseudoinverse,

$L_i$1

and for simultaneous estimation uses

$L_i$2

In the rank-one case, if $L_i$3, then

$L_i$4

so only the combination $L_i$5 is estimable (Namkung et al., 2024).

The physical meaning of singularity depends on its origin. One possibility is genuine probe insensitivity: the state encodes only one effective combination of $L_i$6. Another is a coordinate singularity, where the physical state manifold is regular but the chosen chart becomes degenerate. A third is nonregularity from parameter-dependent rank changes, where the QFIM can become discontinuous and the usual QCRB need not hold in general (Goldberg et al., 2021).

This distinction matters because the remedies differ. If singularity is coordinate-induced, the correct response is reparameterization. If it is physical, the probe state must be changed. The paper on nonregular models explicitly warns that pseudoinverse manipulations may be overly optimistic and that changing measurement protocols can rectify a singular classical FIM but never a singular QFIM, since the latter is already measurement-optimized (Goldberg et al., 2021).

A paradigmatic rank-one two-parameter QFIM is

$L_i$7

which arises when only $L_i$8 is encoded. Another is

$L_i$9

for which only the sum $\partial_i\rho=\frac12(L_i\rho+\rho L_i),$0 is estimable (Namkung et al., 2024).

4. Closed-form qubit constructions

For mixed qubits, an explicit invertibility criterion is available. In the mixed-qubit model

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$1

with eigenvectors parameterized by a mixing angle $\partial_i\rho=\frac12(L_i\rho+\rho L_i),$2 and a constant phase $\partial_i\rho=\frac12(L_i\rho+\rho L_i),$3, the two-parameter QFIM takes the form

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$4

with

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$5

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$6

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$7

Its determinant is

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$8

Hence the QFIM is singular exactly when

$\partial_i\rho=\frac12(L_i\rho+\rho L_i),$9

equivalently when the Jacobian of the map $$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$0 vanishes. In Bloch-vector form the same condition becomes

$$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$1

The paper treats invertibility as a necessary condition for simultaneous estimability, not a sufficient one (Yue et al., 2016).

Its dissipative-qubit example makes the criterion explicit. For a decaying two-level system with parameters $$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$2 and $$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$3, the QFIM satisfies

$$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$4

so the original parameter pair is not jointly identifiable; after orthogonal diagonalization, only one linear combination has nonzero Fisher information (Yue et al., 2016).

For unitary single-qubit $$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$5 processes, the QFIM acquires a geometric form. Writing the unitary generators as

$$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$6

and characterizing the probe by eigenvalues $$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$7 and Bloch direction $$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$8, the paper gives

$$F_{ij}=\frac12\operatorname{Tr}\!\big[\rho(L_iL_j+L_jL_i)\big],\qquad i,j\in\{1,2\}.$$9

In the two-parameter case this is a projected Gram matrix of the two generator Bloch vectors in the plane orthogonal to $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$0. The off-diagonal element is the projected inner product, so $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$1 when the effective generator directions are orthogonal in that transverse plane. The paper further emphasizes that linear independence of the bare generators is not sufficient for invertibility, because projection onto the probe-dependent plane can still produce a singular $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$2 matrix (Shemshadi et al., 2018).

5. Computation and measurement-induced Fisher geometry

Analytical evaluation of a two-parameter QFIM is often difficult because standard formulas require diagonalization of $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$3. One major alternative is vectorization. If

$$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$4

then

$$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$5

and

$$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$6

This representation is used directly in two-parameter Heisenberg $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$7 models and in qubit Bloch-vector calculations, because it bypasses spectral decomposition (Bakmou et al., 2019).

A more general framework dispenses with both diagonalization and orthonormal bases. For a linearly independent, possibly non-orthogonal basis $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$8 spanning the support of $$F=\begin{pmatrix} F_{11} & F_{12}\ F_{21} & F_{22} \end{pmatrix}, \qquad F_{12}=F_{21},$$9, with Gram matrix $F$0, operator multiplication and trace are represented as

$F$1

The SLD Lyapunov equations are then solved by two matrix inverses, and the QFIM is obtained from

$F$2

This method is designed for arbitrary-rank density matrices and is particularly useful when the state is naturally expressed in non-orthogonal states, as in discrete quantum imaging or coherent-state models (Fiderer et al., 2020).

For pure-state variational circuits, a different route uses Stein’s identity. There the paper identifies the QFIM with four times the Fubini–Study metric tensor and uses the overlap Hessian

$F$3

to build stochastic estimators of the full matrix. The stated motivation is that direct construction scales as $F$4 with the number of parameters $F$5, whereas the Stein framework reduces the computational complexity to a constant; the approximation error of the estimate scales as $F$6 with the number of perturbation samples $F$7 (Halla, 24 Feb 2025).

The relation between quantum and measurement-specific classical Fisher matrices can also be expressed geometrically. For a fixed POVM $F$8, with SLDs $F$9 defined at a full-rank reference state $f$0, the induced classical Fisher matrix is

$f$1

whereas the quantum one is

$f$2

Using the $f$3-dependent inner product $f$4 and the frame operator $f$5 of the POVM, these become

$f$6

In a two-parameter model, the best and worst estimated linear combinations are the generalized eigenvectors of the pair $f$7. The paper does not establish a full multiparameter matrix inequality or analyze Holevo-type compatibility, but it does show that a fixed informationally complete measurement acts as a contraction on tangent space and therefore degrades every nontrivial parameter direction relative to the SLD optimum (Saini et al., 17 Dec 2025).

6. Representative applications and interpretive issues

In quantum interferometry, the two-parameter QFIM appears when both arm phases are treated as unknown. Writing the sum and difference phases as

$f$8

the pure-state QFIM is

$f$9

and the effective two-parameter Fisher information relevant for the differential phase is the Schur complement

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$0

This framework is the correct one when the common phase is a nuisance parameter and no external reference is available. The same paper argues, however, that the single-parameter QFI retains physical meaning when an external phase reference is genuinely part of the setup, so the two-parameter QFIM is not a universal replacement for single-parameter analysis (Ataman, 2020).

Two-qubit noisy systems provide explicit jointly estimable parameter pairs. In the model

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$1

the parameters $$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$2 and $$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$3 are estimated simultaneously after open-system evolution. The QFIM

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$4

has a generally nonzero cross term

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$5

so the two parameters are statistically coupled. Yet the same model satisfies

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$6

which means that the multiparameter QCRB is locally saturable. The paper further identifies the $$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$7-coherence

$$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$8

and the purity $$J(\theta)_{ij} = \operatorname{Tr}\,\big(J_{\rho(\theta)}^f\big)^{-1}\!\big(\partial_i\rho(\theta)\big)\,\partial_j\rho(\theta),$$9 as key resources for optimal multiparameter estimation (Jahromi et al., 2018).

Thermal Heisenberg $\rho(\theta_1,\theta_2)$00 models furnish another class of exact two-parameter QFIMs. The paper studies $\rho(\theta_1,\theta_2)$01 in the anisotropic model and $\rho(\theta_1,\theta_2)$02 in the isotropic model with magnetic field, computes the full matrices by vectorization, and finds that the off-diagonal terms are generally nonzero. In both cases the relevant SLDs commute, so the QCRB is saturable, and the authors conclude that simultaneous estimation is always advantageous over separate estimation for the studied thermal states (Bakmou et al., 2019).

Discrete quantum imaging yields especially clear two-parameter subproblems. For two incoherent point sources with relative and centroid coordinates, the lowest-order paraxial QFIM shows that relative coordinates are decoupled from the relative intensity, while unequal source brightness couples relative and centroid coordinates. For three equidistant sources with unknown intensities $\rho(\theta_1,\theta_2)$03, the paper gives an explicit $\rho(\theta_1,\theta_2)$04 QFIM,

$\rho(\theta_1,\theta_2)$05

making the intensity-coupling term explicit without diagonalizing the density matrix (Fiderer et al., 2020).

Across these examples, several misconceptions are corrected by the literature. Invertibility is only a necessary condition for meaningful joint estimation at the QFIM level, not a sufficient condition for attainability (Yue et al., 2016). Vanishing off-diagonal entries imply statistical independence in the QFIM, but not by themselves full multiparameter compatibility (Fiderer et al., 2020). Conversely, nonzero off-diagonal terms indicate coupling, but they do not automatically preclude saturability, since commuting or weakly commuting SLDs can still exist (Jahromi et al., 2018). This suggests that the two-parameter QFIM should be read not merely as a precision matrix, but as a compact diagnostic of local quantum geometry, parameter identifiability, and measurement structure.