Universal Information Hydrodynamics
A common generator for Markov, Fokker-Planck and GKLS flows
Abstract
We work in the universal information hydrodynamics (UIH) framework developed in two companion papers, in which reversible quantum dynamics and irreversible response are organised in a local metriplectic geometry driven by a convex free energy, a Fisher information metric, and an antisymmetric no-work channel. Within this setting the irreversible drift is always a Fisher gradient flow in a weighted geometry, equipped with a cost entropy inequality and a curvature coercivity bound that are invariant under the addition of compatible reversible dynamics. In the present paper we show that this structure has a canonical realisation for Fokker Planck equations, finite Markov chains and GKLS semigroups, and that it can be reconstructed operationally from process and state tomography on real quantum hardware.
On the classical side we identify, for any reversible Markov generator with stationary law , a Fisher Dirichlet operator that saturates the cost entropy inequality mode by mode and generates a Fisher gradient flow of the relative entropy in the hydrodynamic Fokker Planck limit. On the quantum side we show that, in a wide class of thermal, nonreversible and coherent GKLS models, the symmetric part of the real generator in the stationary BKM metric has a density block that coincides exactly with this classical Fisher operator, so that population dynamics are always classical Fisher information hydrodynamics while circulation and coherence effects live in the skew part . We prove a finite dimensional UIH hypocoercivity theorem in which the large time decay rate of the semigroup generated by is bounded below by a positive multiple of the Fisher spectral gap of , and we construct a simple renormalisation scheme that preserves Fisher dissipation on chosen observables and exhibits a robust diffusive Fisher universality basin. In one dimensional Fisher diffusions this scheme supports a Fisher analogue of the Jarzynski relation with renormalisation stable free energy differences.
We implement a suite of IBM Quantum experiments that realise this picture on a superconducting qubit. Process tomography of idle and driven circuits yields an effective real generator in a Pauli basis whose metric adjoint split in the BKM metric at the device stationary state produces a clean metriplectic decomposition ; independent idle depths are consistent with a single time homogeneous generator; the dissipative spectrum of sets an information theoretic speed limit for BKM Fisher decay; and a curvature test shows that the same BKM metric is the local Hessian of quantum relative entropy at the stationary state. On both synthetic models and hardware reconstructions the BKM metric, the Hamiltonian symplectic form and a compatible complex structure assemble into a Kähler triple that is approximately preserved by UIH coarse graining, while cross coherence diagnostics turn the reconstructed generators into a practical UIH universality spectrometer. All claims are backed by a reproducible Python archive that includes constructive Fokker Planck to Markov to GKLS lifts, hypocoercive decay diagnostics, fluctuation and renormalisation experiments, and hardware tomography. Taken together, these results provide evidence that the UIH programme is both mathematically coherent and experimentally visible, with emerging signatures of universality in gaps, fluctuation free energies and Kähler RG structure.
Contents
1 Introduction
Open quantum systems are most often described either at the microscopic level by GKLS master equations for density matrices or at the macroscopic level by irreversible transport equations and stochastic differential equations. Between these two descriptions lies a growing body of work in which dissipative evolution is organised as a gradient flow in an information geometry, with Fisher metrics, entropy production and hypocoercive decay replacing bare operator norms as the natural objects. In practice, however, the reversible and irreversible channels of a given open-system model are rarely assembled into a single geometric structure, and it is not clear a priori which parts of that structure are universal and which depend on microscopic details.
In two companion papers [1, 2] we introduced a universal information hydrodynamics (UIH) framework that begins to address this question. The reversible channel is written as a Hamiltonian flow on a pair of density and phase variables, while the irreversible channel is generated by a symmetric mobility operator that determines both a Fisher information metric and a weighted geometry. A single convex free energy drives the irreversible drift, and the combination of with an antisymmetric no-work operator produces a local metriplectic structure. Within this setting we proved a cost entropy inequality that links control cost to entropy production, and a curvature coercivity bound that links Fisher curvature to the irreversible drift. Both statements are invariant under the addition of any reversible flow generated by an antisymmetric operator satisfying a weighted Liouville condition.
The first aim of the present paper is to show that this abstract structure has a canonical realisation in both the Fokker–Planck and reversible Markov settings, and a numerically rigid realisation across a wide class of GKLS semigroups. On the classical side, once a reversible Markov generator and its stationary distribution are fixed, there is a unique Fisher–Dirichlet operator
that saturates the cost entropy inequality mode by mode and realises the irreversible channel as a Fisher gradient flow. In the hydrodynamic limit this operator generates an overdamped Fokker–Planck equation with the same free energy, Fisher metric and Dirichlet structure, placing discrete reversible chains and continuum Fokker–Planck equations on a single UIH ladder.
On the quantum side we show that for a broad class of GKLS models, including coherent and nonreversible examples, the density block of the metric-symmetric part of the real generator in the stationary BKM metric is exactly this same classical Fisher operator. Population dynamics are therefore always classical Fisher information hydrodynamics, with all coherence and circulation effects living in the antisymmetric part . We construct explicit lifts from continuum Fokker–Planck flows to reversible Markov chains and diagonal GKLS generators, and show that coherent Hamiltonian dressing leaves the density-sector Fisher geometry unchanged. Many distinct GKLS models therefore share the same irreversible Fisher hydrodynamics, differing only in their skew channels.
The second aim is to study the large-time behaviour of the full UIH generator . We prove a finite dimensional UIH hypocoercivity theorem in which the asymptotic decay rate of the semigroup generated by is bounded below by a positive multiple of the Fisher spectral gap of , with the proportionality constant depending only on dimensionless coupling parameters that measure the strength of relative to . Once the stationary geometry and Fisher gap are fixed, coherent circulation can reshape trajectories and induce non-normal transients but cannot arbitrarily slow the irreversible drift along density modes. We then introduce a UIH renormalisation group that coarse-grains generators while preserving Fisher dissipation on chosen observables, and we show numerically that a wide class of Fokker–Planck, Markov and GKLS models flow to a diffusive Fisher universality basin. In one-dimensional Fisher diffusions this provides a Fisher analogue of the Jarzynski relation with renormalisation-stable free energy differences.
The third aim is to demonstrate that this picture can be reconstructed operationally. Using IBM Quantum hardware we perform process and state tomography of idle and driven circuits on a superconducting qubit, reconstruct effective short-time generators via matrix logarithms, and obtain clean metriplectic decompositions in the stationary BKM metric. Independent idle depths agree with a single time-homogeneous generator. The dissipative spectrum of sets an information-theoretic speed limit for the decay of the BKM quadratic functional, and a curvature diagnostic confirms that the same BKM metric is the local Hessian of quantum relative entropy. A complex-structure test shows that the BKM metric, Hamiltonian symplectic form and a compatible complex structure organise into an approximate Kähler triple that is preserved by Fisher-preserving coarse-graining. Cross coherence diagnostics turn the reconstructed generators into a practical UIH universality spectrometer, distinguishing full UIH closures from diagonal baselines and identifying Fisher-active modes.
All results are supported by a reproducible Python archive that includes constructive Fokker–Planck to Markov to GKLS lifts, hypocoercive decay diagnostics, renormalisation experiments, Fisher–Jarzynski and Kähler RG tests, IBM -tomography, curvature and BKM speed-limit diagnostics, and permanent links to the code in Appendix A.
1.1 Reader roadmap
Section 2 recalls reversible Fisher hydrodynamics and the metriplectic framework, fixing notation for the Fisher metric, reversible current and entropy production. It emphasises that the same Fisher geometry underlies both the reversible Madelung picture and the irreversible metriplectic structure.
Section 3 introduces finite dimensional information manifolds and shows how a UIH generator on perturbations decomposes as into a Fisher-gradient part and a Hamiltonian part . It defines the metric, the Fisher gap on the traceless subspace, and the basic hypocoercive setting used in the rest of the paper.
Section 4 constructs diagonal GKLS generators from detailed-balance Markov chains and identifies the Fisher–Dirichlet density block as the unique UIH realisation of the density sector.
Section 5 takes the hydrodynamic limit to overdamped Fokker–Planck equations, showing that the same Fisher–Dirichlet operator governs diffusion in the continuum.
Section 6 treats coherent GKLS models with genuine off-diagonal operators, showing that the density sector again closes to a classical Fisher generator while circulation and coherent effects live in the skew channel .
Section 7 develops the UIH hypocoercivity theorem for finite dimensional generators, bounding the decay rate of in terms of the Fisher gap and dimensionless UIH couplings built from and .
Section 8 introduces the UIH renormalisation group, preserves Fisher dissipation on chosen observables, and identifies a diffusive Fisher universality basin in the plane.
Section 9 presents numerical tests of decay clocks, Fisher floors and coupling-parameter behaviour across representative Markov, Fokker–Planck and GKLS examples.
Section 10 implements UIH channel tomography and universality spectroscopy on IBM Quantum hardware: process tomography of idle and driven channels is used to reconstruct in the BKM geometry, split it into and , and estimate together with an emergent Kähler structure and holomorphicity defects. Section 11 summarises the UIH picture, discusses falsifiers, and outlines extensions to larger probe sectors and many-body systems.
2 Background from the reversible and dissipative Fisher geometry
This section gathers the parts of the companion work that will be used in the sequel. The aim is to make the present paper logically self contained for a reader who is willing to accept the main theorems of the previous studies, without repeating proofs. We therefore state the reversible Fisher Schrödinger structure on , the dissipative Fisher metriplectic structure on , and the way in which both are organised by a single underlying Fisher information object.
Commentary. This section only fixes the geometric stage. There is one Fisher geometry on densities, and it supports two kinds of motion. The antisymmetric operator generates the reversible Schrödinger flow at fixed entropy; the symmetric operator generates the dissipative gradient flow that relaxes entropy. Later, when we look at Markov chains, Fokker-Planck equations and GKLS generators, we are simply realising these same two Fisher quadratures in concrete models.
2.1 Reversible Fisher Schrödinger sector
The reversible companion paper [1] works with hydrodynamic variables and on configuration space, treated as probability density and phase potential. Dynamics are generated by a Poisson bracket on functionals of Dubrovin Novikov type, restricted by locality, Euclidean invariance, global phase symmetry, and reversibility. Within that class the classification reduces the bracket to the canonical form
up to equivalence inside the axioms. In particular, .
The Hamiltonian functionals of interest are of the form
where the first term is the usual kinetic energy in hydrodynamic variables, the second is a potential term, and is a local curvature functional built from and its gradient.
The equations of motion are the continuity equation and a Hamilton Jacobi equation with an extra potential coming from . Reversibility and the projective linearity requirement at the complex level restrict to be the Fisher functional on . Concretely, one finds that the only admissible local functional of the form
whose Euler Lagrange contribution can be absorbed into a linear complex Schrödinger equation after a local complexifier is the Fisher functional
with playing the role of . The associated potential
is then the usual quantum potential. With this choice, the equations of motion for are equivalent to the linear Schrödinger equation for , and the scale is fixed by many body and symmetry arguments.
At this reversible level the Fisher object is therefore a curvature term in the Hamilton Jacobi equation. It is determined by axioms that know nothing about Lindblad operators or dissipation. It is precisely this object that will reappear in the dissipative and hydrodynamic sectors.
2.2 Dissipative Fisher metriplectic sector
The dissipative companion paper moves to a density only description. The state space is a set of strictly positive densities on a domain , with fixed total mass and admissible boundary classes (periodic or no flux).
A free energy functional induces a chemical potential up to an additive constant. The irreversible part of the dynamics is assumed to be local and quadratic in at fixed , with a symmetric positive definite mobility tensor that is uniformly elliptic on the domain.
Under a short list of axioms, including a steepest descent principle at fixed state and locality of probes, one shows that the only possible irreversible drift is the continuity form
and that the instantaneous irreversible power
induces a weighted norm on conservative tangents at fixed . The associated entropy production is
and a cost entropy inequality relates , the minimal quadratic cost to realise a given tangent , and the instantaneous dissipation . A curvature coercivity bound further shows that the Hessian of at is bounded below by a curvature constant times the squared norm of .
For later use it is convenient to record a mild interpretation of the cost–entropy relation as an instantaneous dissipative speed constraint. At a fixed state the sharp inequality
| (2.1) |
shows that, for a prescribed entropy–production scalar , the free–energy decay rate cannot be increased independently of the minimal quadratic cost of the flux realising . The UIH irreversible drift is the unique direction that saturates this bound in the weighted geometry. We will only appeal to this as a local geometric way to read the inequality, rather than as a separate postulate about thermodynamic optimisation.
The reversible drift at the density level is represented by an antisymmetric operator that generates fluxes of the form
subject to a weighted Liouville constraint and the no work condition . Under these conditions the reversible sector is orthogonal to the irreversible cone in the inner product, and the instantaneous scalars , and the curvature constant depend only on and .
The Fisher information reappears in this dissipative sector as the curvature of and as the metric underlying the geometry. In particular, when is chosen to be a relative entropy functional with respect to a reference density and is taken to be a scalar mobility, one recovers the standard Fokker Planck and Wasserstein gradient flow structures.
It is convenient to record at this point a simple entropy monotonicity identity that will be used implicitly later. When the physical evolution is chosen to be the –gradient flow of the free energy , so that
the free energy is automatically a Lyapunov functional. Differentiating along the trajectory and using the continuity equation gives
where is the entropy production defined above. In particular, for the canonical choice of a relative entropy functional with respect to a fixed smooth reference density and a scalar mobility , this recovers the familiar statement that relative Shannon entropy is nonincreasing along Fokker-Planck flows, with entropy production . Within the present axioms this “Second Law” is therefore nothing more than a geometric monotonicity property of the Fisher metric: irreversibility corresponds to a nondegenerate positive operator on the tangent space, and the associated relative entropy functional decreases along its gradient flow.
The companion paper goes further and records detailed diagnostics and falsifiers, but for the present work it is enough to know that within the axioms the symmetric and antisymmetric are fixed up to simple gauges, and that the irreversible and reversible parts can be cleanly separated.
2.3 Common Fisher structure
Both companion papers [1, 2] are built around a single Fisher object. In the reversible Schrödinger sector Fisher curvature appears as a regulariser of the Hamilton Jacobi equation and is responsible for the quantum potential.
In the dissipative metriplectic sector Fisher information appears as the curvature of the free energy and as the metric that underlies the geometry and its cost entropy relations. The present paper does not introduce a new Fisher functional. It instead explores how these existing structures interact with Lindblad generators and hydrodynamic limits.
At a practical level this means that all constructions in the sequel will assume a fixed Fisher geometry on the density manifold. The symmetric mobility is the same one that appears in the metriplectic study.
The antisymmetric operator is the same one that appears in the reversible study, at least at the level of coarse grained densities. The new object that will be introduced in the next section is therefore not an extra degree of freedom. It is a complex packing of existing operators that simplifies the bookkeeping when one wants to talk about reversible and irreversible quadratures in the same breath.
In addition, the same Fisher structure reduces to the usual parametric Fisher information for one dimensional translation families. Let be a smooth, strictly positive density on and consider the translations . The associated configuration space Fisher functional in the density sector is
| (2.2) |
which induces a one dimensional Fisher metric on the parameter . The companion Fisher Kähler paper [3] shows that this metric coincides exactly with the parametric Fisher information
| (2.3) |
so that for translation channels. In the EPI language this identifies the UIH Fisher functional with the Fisher information used in Extreme Physical Information; in the present operator language it means that the single density sector Fisher geometry supporting the decomposition simultaneously controls configuration space curvature, entropy production, and parametric information for translation invariant channels.111Throughout this paper, the complex form is used when acting on wavefunctions or density matrices, while the real metric-adjoint form (with absorbing the factor of ) is used for perturbations in the metrised Hilbert space. The two conventions are reconciled explicitly in Section 9.
3 Information manifolds and UIH flows
We now introduce the language in which the rest of the paper will be written.
We regard the set of states, classical or quantum, as a manifold equipped with an information metric and a complex mobility operator. The symmetric part of the mobility encodes irreversible gradient flow with respect to a free energy, the antisymmetric part encodes reversible Hamiltonian flow, and Lindblad generators appear as particular ways of coupling these parts when the states are density matrices rather than scalar densities.
3.1 Information manifolds
An information manifold is a smooth manifold of states equipped with three pieces of structure. First, there is a free energy functional , usually a relative entropy or energy functional that plays the role of a Lyapunov function for the irreversible dynamics.
Second, there is an information metric on the tangent bundle , typically a Fisher type Riemannian metric that arises as the Hessian of or of a related divergence. Third, there is a compatibility between and that allows one to define a gradient flow of on with respect to .
In the classical density setting is the space of strictly positive densities with fixed mass on a domain . The free energy is often taken to be a relative entropy plus potential terms, and the Fisher metric can be realised through the weighted geometry defined by the elliptic operators and . In the quantum setting can be taken as the convex set of density matrices on a Hilbert space. The role of is played by the Umegaki relative entropy or a related quantity, and the Fisher metric becomes a quantum information metric on density matrices. For the explicit examples in this paper we will mostly work in the classical setting, but the quantum examples inherit the same logic.
Commentary. Here we forget about PDEs and work in a finite-dimensional laboratory. We take a vector space with a metric and split a linear operator into a symmetric part and a skew part . The symmetric is the “friction” or downhill part in the metric, the skew is the “rotation” or sideways part, and just packages them together. All of the -split statements in this section are linear algebra, but the same patterns will reappear for Markov chains, Fokker-Planck flows and GKLS semigroups.
On such a manifold the gradient is defined implicitly by the relation
for all tangent vectors . In the density setting the metric is represented in terms of the elliptic operators by
where and are potentials solving Poisson type equations , and similarly for . The metriplectic paper shows that, under the stated axioms, the irreversible drift is exactly and is represented as a continuity equation with .
3.2 Complex mobility operator
To prepare for the Lindblad discussion it is convenient to lift this gradient flow picture to the cotangent space. Instead of working directly with tangent vectors , we consider covectors that play the role of chemical potentials or Hamiltonian potentials. The metric and its inverse identify tangent and cotangent spaces, and the irreversible drift can be written schematically as where is the gradient of and is a positive selfadjoint operator with respect to .
The reversible sector meanwhile is encoded by an antisymmetric operator that maps covectors to tangent vectors in such a way that the power vanishes for all . In the density setting is represented via the weighted Liouville identity and generates divergence free fluxes. In the quantum setting corresponds to the commutator with a Hamiltonian and generates unitary evolution that preserves the von Neumann entropy.
The complex mobility operator is then defined by
Acting on a complex potential
with the chemical potential and a Hamiltonian potential, produces a complex current of the form
which can then be inserted into a continuity equation. In the density case this gives
with decomposed into reversible and irreversible parts. In the Schrödinger case, after the complexifier, the same structure can be read at the level of the wave function. In the Lindblad case, acts on matrix valued potentials and the continuity equation is replaced by the GKLS master equation.
The usefulness of is twofold. Conceptually, it packages the symmetric and antisymmetric mobility into a single object that carries both quadratures. Practically, it allows one to write diagnostics and inequalities in a compact way. For example, the cost entropy inequality and curvature coercivity bounds of the metriplectic paper can be viewed as statements about the real part of on the real axis of potentials, while reversible invariance corresponds to the imaginary part dropping out of certain scalar functionals.
In linear response one can form, for suitable observables, a susceptibility
where the inner product is taken with respect to the metric and encode the chosen probe and readout. Whenever the associated semigroup is causal and extends analytically to the upper half plane, standard dispersion theory implies a Kramers-Kronig relation between the real and imaginary parts of . In the present split this means that the dissipative and reversible quadratures encoded by and appear as the two Hilbert transform linked components of a single causal response function rather than as freely tunable additions. The concrete density sector tests in the metriplectic companion paper provide one instance of this general picture.
3.3 Schrödinger, Fokker Planck, and Lindblad as flows
Within this abstract language, the three families of dynamics that concern us can be described informally as follows.
First, Schrödinger evolution on pure states corresponds to a purely imaginary with and given by the Hamiltonian commutator in the appropriate representation. The free energy is constant, and the state moves along level sets of . At the hydrodynamic level this structure is encoded by the reversible companion paper in the canonical bracket on and the Fisher curvature.
Second, Fokker Planck evolution in the overdamped regime corresponds to a purely real with and given by the Fisher metriplectic mobility. The state moves down the gradient of with respect to the metric. The dissipative companion paper characterises this class in detail and shows that, under the stated axioms, the irreversible flow is uniquely determined by and .
Third, Lindblad evolution on density matrices mixes both parts. In the GKLS representation the generator of a contractive semigroup is a sum of a Hamiltonian commutator and a dissipator built from jump operators. When the jump operators and stationary state obey detailed balance, the dissipator can be interpreted as a gradient flow of the quantum relative entropy with respect to a quantum Fisher metric, while the Hamiltonian part is reversible.
In the special class of models studied in this paper the dissipator reduces on the diagonal to a classical Fisher gradient flow of the Kullback Leibler divergence, and the Hamiltonian part becomes a familiar reversible operator. The resulting dynamics on the diagonal sector can then be seen as a projection of a flow compatible with the classical Fisher geometry.
The rest of the paper will make these schematic statements precise for a concrete family of GKLS generators with diagonal jumps and for their hydrodynamic limits.
Section 4 will introduce the finite dimensional GKLS model, show explicitly how its diagonal sector reduces to a reversible Markov chain, and identify the equality of quantum and classical entropy decay curves.
Section 5 will construct the hydrodynamic limit and show that the limiting Fokker Planck equation is exactly a Fisher metriplectic gradient flow. Section 6 will use a simple coherent two level GKLS model to illustrate how feeds into an effective on coarse grained populations in a strong dephasing regime. The present section is therefore the conceptual backbone that allows these examples to be read as instances of a single universal information hydrodynamics. This finite dimensional picture serves as a laboratory for the UIH K split.
In the next section, we embed it into diagonal GKLS semigroups whose density sector is exactly Fisher-Dirichlet, while in Section 6 we show that the same split arises in coherent GKLS models once the BKM metric at the stationary state is used.
The continuum Fokker-Planck flows of Section 5 and the numerical and hardware tests in the appendices then demonstrate that the same structures control irreversible response from classical overdamped Langevin dynamics all the way to real, physical noisy, superconducting qubits.
4 Finite GKLS with diagonal jumps and classical reversible chains
In this section we give a fully explicit family of finite dimensional GKLS generators that act as Lindblad semigroups on density matrices and reduce exactly to classical reversible Markov chains on their diagonal sector.
For diagonal states, the quantum relative entropy with respect to the stationary state coincides identically with the classical Kullback Leibler divergence, so the entropy decay curve of the GKLS semigroup on that sector is exactly the same as that of the corresponding Markov chain. This realises the discrete Fisher gradient flow on the probability simplex as a special case of a quantum information hydrodynamics, in a way that is completely transparent at the level of matrix elements and supported numerically by the script 01_gkls_diagonal_to_markov_checks.py.
Commentary. This section answers a concrete question: what does the UIH picture look like for an honest Lindblad jump model? We start from a thermal GKLS generator, restrict the dissipator to the energy eigenbasis to obtain a classical Markov generator , and then symmetrise with respect to the stationary law to obtain the Dirichlet operator . The Fisher metric and entropy decay in the diagonal sector are not fitted by hand; they are determined entirely by the GKLS data.
4.1 Set up and notation
Let be a finite dimensional Hilbert space with orthonormal basis . A density matrix is a positive semidefinite operator on with unit trace. We will work in the matrix representation of in the basis :
Fix a strictly positive probability vector with . Think of as a discrete Gibbs distribution or stationary measure.
We introduce non negative rates for all ordered pairs . These rates will act as the jump intensities of a continuous time Markov chain in the classical picture, and as the amplitudes of Lindblad jump operators in the quantum picture. A convenient way to enforce detailed balance with respect to is to start from a symmetric matrix with zero diagonal and set
| (4.1) |
Then
| (4.2) |
so detailed balance holds automatically.
On the quantum side we define Lindblad jump operators
| (4.3) |
and we consider a GKLS generator with these jumps and no Hamiltonian term. That is, the evolution of a density matrix is governed by
| (4.4) |
We will show that this semigroup has a unique stationary state
and that its action on diagonal density matrices coincides exactly with that of a reversible Markov chain with generator built from .
4.2 Diagonal sector and reversible Markov chain
We first show that the diagonal and off diagonal matrix elements of evolve in a simple and decoupled way under (4.4). Using the definitions (4.1) and (4.3), one checks
Writing , a short computation shows that the diagonal entries satisfy
| (4.5) |
while the off diagonal entries satisfy
| (4.6) |
The derivation is standard: the term only contributes to diagonal entries, and the anticommutator term involves projectors that act diagonally in the chosen basis.
Equation (4.6) shows that the off diagonal entries decay exponentially to zero with a rate bounded below by half the total outgoing rate from the corresponding indices. They do not feed back into the diagonal entries. Equation (4.5) shows that the diagonal entries are closed under the evolution.
If we define the population vector by
then (4.5) can be written as
| (4.7) |
This is exactly the master equation of a continuous time Markov chain on , with transition rates from to . If we define a generator by
| (4.8) |
then (4.7) can be written as
| (4.9) |
where the column vector has entries . The detailed balance identity (4.2) is equivalent to
so is a reversible generator with invariant distribution .
We also note that the stationary quantum state satisfies . Indeed, in the basis we have
and using detailed balance each term in (4.4) cancels when applied to , in a manner analogous to the usual classical detailed balance identity. Thus is the unique stationary state of in the full quantum system, and its diagonal entries are the unique stationary distribution of the classical chain.
4.3 Quantum relative entropy and classical Kullback Leibler divergence
For a general density matrix the Umegaki quantum relative entropy with respect to is defined by
When both and are diagonal in the same basis, the logarithms act elementwise, so if and , then
| (4.10) |
the classical Kullback Leibler divergence from to .
Suppose now that the initial state is diagonal, so that remains diagonal under the GKLS evolution (4.4) and the corresponding population vector evolves by (4.9). Then for all times we have
| (4.11) |
The decay of quantum relative entropy along the diagonal sector of the Lindblad semigroup is therefore identical, pointwise in time, to the decay of classical relative entropy along the corresponding Markov chain.
The evolution (4.9) is generated by a reversible , so it admits a well known Dirichlet form representation. One can show that
| (4.12) |
Every term in the sum is non negative by the usual monotonicity of the logarithm, so the derivative is non positive, and is zero if and only if is constant in . Thus decays strictly along the Markov chain unless , and the same holds for the quantum relative entropy of diagonal states under the GKLS evolution.
In the language of the metriplectic companion paper, (4.12) expresses the entropy production as a discrete Fisher quadratic form in the logarithmic gradient of . Indeed, the difference is a discrete gradient in the tilted variable, and the factor is a discrete mobility. The entire finite dimensional picture is thus a discrete counterpart of the Fisher metric and geometry that appear in the continuum.
4.4 Numerical verification
The script 01_gkls_diagonal_to_markov_checks.py implements this construction for concrete choices of , and . It constructs the Lindblad superoperator in vectorised form, builds the Markov generator , and evolves both the GKLS equation and the Markov master equation from the same diagonal initial state. At a fixed set of times it compares:
-
•
the population vectors extracted from the diagonal of and from the Markov chain, and
-
•
the relative entropies and .
The script reports the maximum difference in populations and in relative entropy over the time grid. In all tested cases these differences are at the level of numerical integration error, confirming (4.9) and (4.11) at the level of floating point computations. For completeness, the script also checks positivity and normalisation of the density matrix along the evolution.
5 Hydrodynamic limit and Fisher metriplectic Fokker Planck flows
We now move from finite dimension to continuum. We keep the classical reversible chain constructed above but embed it in a family indexed by a lattice spacing and endowed with a nearest neighbour structure. In the limit the master equation of the chain converges to an overdamped Fokker Planck equation with drift determined by a potential and diffusion coefficient . We then show that this Fokker Planck equation is exactly a Fisher metriplectic gradient flow with mobility tensor and free energy , as described in the dissipative companion paper. Scripts 02_markov_to_fp_limit_checks.py and 03_fp_fisher_metriplectic_checks.py provide numerical support.
Commentary. The finite reversible chain from Section 4 is now sent to a continuum limit. Jumps turn into drift and diffusion, and the discrete Fisher-Dirichlet form becomes the quadratic form of an overdamped Langevin / Fokker-Planck operator. The same Fisher geometry and -split survive the limit. In the weighted Fisher geometry the continuum Fokker-Planck operator again takes the form , with giving entropy production and giving incompressible transport.
5.1 Reversible chain on a lattice
For simplicity we work in one spatial dimension and with a periodic or reflecting domain. Fix a macroscopic interval and a lattice spacing , and let
where . We consider indices modulo in the periodic case.
Let be a smooth potential, and fix a diffusion coefficient . Define the discrete Gibbs weights
We will write when the dependence on is clear from context.
We now define nearest neighbour rates that satisfy detailed balance with respect to and have the natural diffusive scaling in . For each set
| (5.1) |
This is the one dimensional specialisation of (4.1). It satisfies detailed balance in the sense that
The generator of the chain is then
with all other entries zero. The master equation for the population vector is
| (5.2) |
The stationary distribution of this chain is , and detailed balance ensures reversibility.
5.2 Flux form and continuum limit
To extract a continuum equation, it is convenient to rewrite (5.2) in flux form. Define the flux across the bond between sites and by
Then (5.2) can be written as
| (5.3) |
Introduce the tilted variables
so that . Using (5.1) and detailed balance, we can write the flux as
This is a discrete gradient in with a mobility prefactor .
We now define a piecewise constant density by
so that approximates . Similarly define , so that approaches the continuum Gibbs density as .
In the formal limit , one shows that
where , and
where . Thus
Using (5.3) and dividing by we have
where in the limit we can identify
The continuum equation is therefore
Since , we have , so
| (5.4) |
This is the standard overdamped Fokker Planck equation with drift potential and diffusion coefficient .
The same argument extends to higher dimensions and more general neighbour graphs. The key point is that the diffusive scaling and the detailed balance structure ensure that the macroscopic limit is a second order operator of the form , which is the continuum Fisher metriplectic drift with and free energy built from the Gibbs density.
5.3 Fokker Planck as a Fisher metriplectic K flow
We now identify the Fokker Planck equation (5.4) with a Fisher metriplectic K flow in the sense of the dissipative companion paper. The generalisation to higher dimensions is straightforward.
Let be a smooth strictly positive density on and define the continuum Gibbs density
Define the free energy functional
| (5.5) |
A variation with yields
so the chemical potential is
| (5.6) |
The constant does not affect gradients and can be ignored.
We choose the mobility tensor acting as a scalar multiple of the identity. In one dimension this gives
The irreversible Fisher metriplectic drift is then
| (5.7) |
Using (5.6), we have
so
Inserting this into (5.7) gives exactly (5.4). Thus the Fokker Planck equation obtained from the hydrodynamic limit of the reversible chains is exactly the Fisher metriplectic gradient flow of with mobility , in the class singled out by the dissipative companion paper.
The entropy production along the flow is
assuming suitable boundary conditions. This is minus the square of the norm of , and coincides with the continuum limit of the discrete Dirichlet form (4.12). In particular, the decay of along the Fokker Planck flow is the continuum counterpart of the decay of Kullback Leibler divergence along the reversible chain, and the discrete Fisher form in (4.12) converges to the continuum Fisher form .
5.4 Numerical verification
The script 02_markov_to_fp_limit_checks.py implements the reversible chains described above for a sequence of lattice spacings and compares them to a numerical solution of the Fokker Planck equation (5.4). For each spacing it:
-
•
constructs and evolves the master equation (5.2) for a fixed initial density profile,
-
•
interpolates the resulting discrete density to a common grid, and
-
•
compares it to the Fokker Planck solution at the same times.
The code reports error norms between and and estimates a convergence rate in using a log log fit. In all tested regimes the errors decay at the expected order dictated by the spatial discretisation, and the discrete relative entropy decays along curves that converge to the continuum decay of .
The script 03_fp_fisher_metriplectic_checks.py works directly at the continuum level. It discretises the Fokker Planck equation (5.4) on a fine grid and, at each time step, computes:
-
•
a finite difference approximation to ,
-
•
the Fisher metric right hand side , and
-
•
the free energy and the quadratic form .
It then checks numerically that matches to within discretisation error, and that along the evolution. These diagnostics provide a concrete numerical confirmation of the Fisher metriplectic structure at the PDE level.
Together, the reversible GKLS to Markov reduction of Section 4 and the hydrodynamic Fisher Fokker Planck limit of this section show that the entropy decay and Fisher geometry of a simple class of Lindblad semigroups can be understood entirely in terms of classical Fisher gradient flows on densities, once one restricts to the diagonal sector and passes to the macroscopic scale.
6 Coherent GKLS models, coherences, and effective Fisher chains
The examples so far have been deliberately classical in flavour. Lindblad generators with diagonal jump operators and no Hamiltonian term act as quantum lifts of reversible Markov chains, and their entropy decay on the diagonal sector is exactly classical. In this section we consider a genuinely coherent GKLS model in which the Hamiltonian and the dissipator are not simultaneously diagonal, so that populations and coherences interact. At the level of the full density matrix the evolution is first order in time and of GKLS type, but at the level of populations alone the dynamics becomes second order and non Markovian. In a strong dephasing regime one can nevertheless eliminate the coherences and obtain an effective classical Markov chain for the populations. This effective chain falls inside the Fisher metriplectic class, and its rates depend on the Hamiltonian part in a way that makes the interplay between and explicit.
Commentary. Up to now the dissipative sector came from GKLS generators that were already diagonal in the energy basis. Here we let a genuinely coherent GKLS model run and ask what the diagonal density alone “sees”. The off-diagonal coherences feed into an effective Markov generator on the populations, and its symmetric block is still a Fisher-Dirichlet operator. In UIH language, the coherent sector renormalises the effective seen by the coarse-grained density, but it does not destroy the underlying Fisher-metriplectic structure.
6.1 A qubit GKLS model with coherences
We work with a single qubit and choose a basis in which the Lindblad operator is diagonal but the Hamiltonian is not. Let be the eigenbasis of , and set
with real parameters and . The GKLS equation reads
| (6.1) |
This generator has a unique stationary state and describes coherent rotations about the axis combined with dephasing in the basis.
It is convenient to write the state in Bloch form,
where the real vector lies in the unit ball. Inserting this parameterisation into (6.1) and using the usual Pauli algebra, one finds that the Bloch components satisfy the linear system
| (6.2) | ||||
| (6.3) | ||||
| (6.4) |
The dissipator damps the and components at rate and leaves unchanged, while the Hamiltonian couples and through precession at frequency . The stationary point is at , corresponding to the maximally mixed state .
The population of the excited state is
and the ground state population is . The diagonal sector is therefore determined by the single function .
Combining (6.3) and (6.4) we can eliminate . Differentiating (6.4) gives
Using from (6.4), we obtain
| (6.5) |
This is the equation of a damped harmonic oscillator with natural frequency and damping coefficient . In terms of the population we have
| (6.6) |
Thus, although the full GKLS equation (6.1) is first order in time, the induced dynamics on the coarse grained population variable is second order and involves an inertial term . There is no exact closed first order Markovian equation for alone at finite ; any attempt to write as a function of only would require keeping track of hidden variables that encode the coherence .
6.2 Overdamped limit and effective two state chain
In the strong dephasing regime , the variables and are fast and strongly damped, while evolves slowly. On timescales large compared to but not so large that has fully relaxed, the system is effectively overdamped and one can eliminate adiabatically from the dynamics.
Formally, on the slow manifold one can set in (6.3) and solve for in terms of :
Substituting this into (6.4) gives an effective first order equation for :
Equivalently,
| (6.7) |
where is a small residual term that vanishes in the asymptotic regime at fixed . In terms of the population this reads
This is exactly the master equation of a symmetric two state Markov chain with states and , stationary distribution , and transition rates
Indeed, the master equation for such a chain is
so identifying coefficients gives and hence .
The effective generator
is reversible with respect to . The free energy functional on the population simplex is the relative entropy
Writing and suppressing the time argument, we can calculate
Thus the discrete chemical potential is
The effective Markov dynamics for can then be written as a one dimensional gradient flow
| (6.8) |
with mobility
| (6.9) |
For and the numerator and denominator have the same sign, so . At the symmetric point one can use l’Hospital’s rule to define the limit
so extends smoothly to a strictly positive function on the open interval. The dynamics (6.8) is therefore a Fisher type one dimensional gradient flow of with a state dependent scalar mobility , entirely analogous to the continuum Fisher metriplectic flows considered in the dissipative companion paper.
This example shows that even a tiny quantum system with genuine coherences produces, in a suitable regime, an effective classical chain whose entropy decay and metric structure sit inside the Fisher metriplectic class. The reversible Hamiltonian part does not appear directly in the irreversible entropy production, but it does feed into the effective mobility through the factor in (6.9).
6.3 Numerical illustration
The script 04_gkls_coherence_elimination_checks.py implements the qubit GKLS model (6.1) in the Bloch representation and compares the exact GKLS evolution to the effective Markov chain dynamics in the overdamped regime.
Concretely, the script:
- •
-
•
extracts the population and fits its long time decay to a single exponential, obtaining an effective rate for the relaxation of to ;
-
•
compares to the theoretical value predicted by (6.7);
-
•
constructs the effective two state Markov chain with rate and integrates its master equation from the same initial population, comparing the resulting to the population curve from the full GKLS dynamics in the overdamped regime.
In regimes where is large, the numerical results show that the slow tail of under the full GKLS evolution is well described by the effective exponential with rate , and that the effective two state chain reproduces the population dynamics on timescales larger than to within the numerical accuracy of the integrator. The script reports the mismatch between the fitted and theoretical rates and the maximum deviation between and over a chosen time window. These diagnostics quantify how the coherent reversible part feeds into the effective irreversible dynamics of populations when coherences are fast variables.
From the perspective of the information manifold, this example realises a projection of a quantum flow onto a one dimensional classical Fisher metriplectic flow, with an effective mobility that depends quadratically on the Hamiltonian amplitude and inversely on the dephasing rate. It provides a concrete instance of the principle that, although reversible operators do not contribute directly to entropy production, they can modify the effective irreversible channel seen by coarse grained observables.
6.4 Signatures on IBM Quantum hardware
The abstract K split and Fisher-Lindblad picture are not restricted to classic computing. To test them on actual quantum hardware we performed a series of IBM Quantum experiments.
We first carry out a K tomography experiment on a noisy idle channel of a superconducting qubit, detailed in Appendix D.1. A depth four identity circuit is reconstructed by one qubit process tomography in the Hermitian Pauli basis, yielding a channel superoperator . From its unique stationary state we build the BKM metric and restrict to the traceless Pauli block . Assuming for some effective generator , we compute by a matrix logarithm and perform the metric adjoint split
On representative runs both the symmetry residual and the skewness residual are at the level of machine precision, while the dissipative spectrum of is strictly positive. This is a direct experimental realisation of the UIH K split on hardware.
A second experiment tests semigroup scaling across two idle depths, showing that the same governs the idle dynamics over a range of time scales to within small channel level deviations compatible with sampling noise and drift, detailed in Appendix D.2.
A third experiment uses the BKM metric and the dissipative spectrum to formulate and verify an information theoretic speed limit on the hardware: the smallest eigenvalue of sets the natural decay clock for the quadratic functional , and the full flow decays faster than the pure gradient flow, as predicted by the one current two quadratures picture, detailed in Appendix D.3.
Finally, a curvature test confirms that the BKM metric extracted from is the local second order curvature of quantum relative entropy for small unitary perturbations, with measured entropy changes tracking the quadratic BKM prediction across three independent directions, detailed in Appendix D.4.
CPTP-repaired two qubit K tomography and semigroup decay.
A final suite of hardware tests targets the genuinely noncommuting two qubit idle channel of the ibm_fez backend. For a depth- identity circuit we perform full two qubit process tomography, reconstruct a completely positive and trace preserving map by a spectral CPTP repair, and extract a traceless block . The unique stationary state of determines the BKM metric and its restriction to the traceless operator space. A matrix logarithm defines an effective generator , and performing the metric-adjoint split
again yields symmetry and skew-symmetry residuals at machine precision, demonstrating that the metric adjoint picture persists in the full two qubit setting. The dissipative Fisher spectrum is then directly accessible as the eigenvalues of .
Three independent datasets (4k, 6k and 8k tomography shots) display the same structure. The Fisher dissipative spectrum obeys with slowest rates between and depending on shot count, while the spectral condition number of remains modest. For each dataset nine distinct initial conditions are constructed by perturbing the stationary state along orthonormal directions in the BKM metric. The evolution is propagated for and its quantum relative entropy is monitored as a function of time.
In all runs the relative entropy exhibits clean exponential decay. Fitting over short windows () and long windows () yields decay rates between and depending on the initial condition. Crucially, every observed decay rate is strictly more negative than the Fisher gap extracted from , with ratios ranging from to . This is the expected signature of the cost-entropy inequality: the Fisher gap provides a geometric lower bound on irreversible decay under the full K flow, and the hardware dynamics lie safely above this floor.
A particularly sharp test uses the “gap mode” initial state obtained by taking an eigenvector of corresponding to and constructing a traceless BKM-orthonormal perturbation of . If the UIH picture were not geometrically faithful for the reconstructed channel, this state would yield the most stringent challenge to the Fisher gap prediction. Instead the behaviour aligns cleanly with theory: the gap mode decays exponentially with rates again significantly faster than the Fisher gap, with no violation in any dataset.
The semigroup extrapolations to inevitably leave the positive cone for some initial conditions, as the matrix logarithm does not enforce GKLS form at intermediate times. However positivity is maintained for early times (), the decay rates are already stable by , and all relevant quantities in the Fisher-Lindblad tests depend only on and on with eigenvalues clipped inside the logarithm. The positivity excursions therefore do not affect the validity of the speed-limit comparison.
Outcome.
These two qubit tests provide an independent and far more demanding validation of the geometric content of the UIH framework. The BKM metric extracted from hardware determines a Fisher dissipative operator whose smallest eigenvalue acts as an experimentally measurable irreversible clock; the full K flow constructed from the same data respects this bound across a wide family of initial conditions; and even the eigenvector-engineered gap mode obeys the same inequality. Together these results show that the Fisher-Lindblad structure, originally derived abstractly from reversible-dissipative metriplectic theory, is realised quantitatively on a real superconducting device, further reinforcing the geometric smoking guns established above.
Together these hardware results provide geometric smoking guns for UIH on superconducting qubits: the K split, the Fisher dissipation spectrum as an irreversible clock, and the identification of the BKM metric as the local information curvature are all realised experimentally on a real device.
6.5 Hydrogen tunneling in palladium as an external Fisher-Lindblad testbed
The two state tunneling experiment of Ozawa et al. [15] on hydrogen in palladium provides a particularly clean external testbed for the Fisher-Lindblad picture developed here. In their setup a 10 nm Pd film is loaded with hydrogen at low temperature, creating a metastable population of H atoms in tetrahedral (T) interstitial sites which relax into energetically favorable octahedral (O) sites upon annealing. The dynamics is read out through the in plane resistance, and for each anneal temperature the resistance trace is accurately fitted by a single exponential,
| (6.10) |
with interpreted as the T to O hopping rate.
By combining channeling nuclear reaction analysis with density functional and path integral calculations, Ozawa et al. show that the relevant hydrogen dynamics can be modeled as a two site system with an energy splitting between an active T level and a higher O vibrational level, coupled to both phonon and conduction electron baths. The measured hopping rate exhibits three regimes: a high temperature Arrhenius law with activation energy attributed to over barrier diffusion, an intermediate Arrhenius regime with a small barrier associated with phonon assisted resonant tunneling, and a low temperature regime where
| (6.11) |
in agreement with Kondo style theory for nonadiabatic electron mediated tunneling in an asymmetric double well.
At the level of occupations the hydrogen subsystem can be represented as a reversible two state Markov chain
| (6.12) |
with detailed balance and nonzero eigenvalue . For such a chain the relative entropy
| (6.13) |
decays exponentially at late times with a universal rate
| (6.14) |
independent of the asymmetry and the initial nonequilibrium state . In particular, in the low temperature electron dominated regime the UIH entropy clock for hydrogen site information inherits the measured scaling
| (6.15) |
This two state Markov representation is exactly the classical density sector of a two level GKLS model in which the coherent tunneling matrix element and the electron induced broadening obey the Kondo scalings , . In the overdamped limit the Bloch equations reduce to an effective rate
| (6.16) |
so that the same exponent simultaneously governs the antisymmetric Hamiltonian part and the symmetric mobility in the Fisher-Lindblad decomposition .
In the present work we treat the published curve as defining the reversible two state generator and hence the Fisher Dirichlet operator for the hydrogen occupancy sector, and we use the resulting as an external example of a Fisher clock with a nontrivial bath exponent. A more stringent information geometric test would require access to the underlying time resolved relaxation traces , from which one could reconstruct directly and verify both the entropy clock relation and the saturation of the cost entropy inequality along the physical trajectory. We therefore propose PdH tunneling as a natural condensed matter target for future Fisher-Lindblad analysis once such data are available.
A self contained numerical demonstration of the Markov gap entropy clock for the Ozawa style two state PdH model is provided in the code archive as 52_pdh_entropy_clock_demo.py; see Appendix A. This script constructs the reversible two state generator from the published fit parameters, computes the relative entropy decay rate and verifies across the full temperature range. It can be used as a simple, referee ready check that the PdH tunneling testbed realises the UIH entropy clock in the classical density sector.
7 A finite dimensional UIH hypocoercivity theorem
In this subsection we formulate a finite dimensional hypocoercivity theorem in the UIH language. The setting covers the reversible Markov chains, the diagonal GKLS models and the BKM Pauli blocks that appear throughout the paper, and it is the natural place where the Fisher gap and the one current two quadratures picture combine to give a quantitative exponential decay rate for the full non normal generator.
We work on a real vector space equipped with a symmetric positive definite matrix , regarded as a Riemannian metric via . Let be a real matrix and denote its metric adjoint by
| (7.1) |
We assume that admits a metric symmetric plus metric skew split
| (7.2) |
so that is the dissipative part and the reversible part in the UIH sense. The evolution equation
| (7.3) |
is then a linear UIH flow on , with associated quadratic Fisher functional
| (7.4) |
In all examples of interest there is a distinguished stationary direction, given either by the constant vector in Markov models or by the identity operator in GKLS sectors. Let denote such a stationary direction and let be its orthogonal complement,
| (7.5) |
which is invariant under . We write for the restrictions to , and we work entirely on from now on. On the Fisher functional is strictly positive definite.
Fisher gap
On the traceless subspace the dissipative part is symmetric and negative definite, so that the spectrum of is contained in . The Fisher gap is the smallest eigenvalue
| (7.6) |
Equivalently, for all one has the Fisher Dirichlet inequality
| (7.7) |
The pure gradient flow is therefore strictly contracting in the Fisher metric, with decay rate at least . The hypocoercivity problem is to show that the full UIH flow inherits exponential decay, with a rate that is controlled from below in terms of and the reversible part .
To state the theorem in a form that is uniform over UIH models, it is convenient to use norms induced by the metric . For a linear map we define the operator norm
| (7.8) |
and the associated quadratic form norm
| (7.9) |
All norms on the finite dimensional space of operators are equivalent, so any choice here would do, but these metric norms make the estimates invariant under UIH coordinate changes.
Theorem 7.1 (Finite dimensional UIH hypocoercivity — constructive form).
Let , , be as above, with symmetric positive definite, strictly negative definite on and . Let be the Fisher gap of , and set
| (7.10) |
Then there exist constants and , depending only on , and , with the following properties.
-
1.
There exists a real matrix on solving the Sylvester equation
(7.11) For every the modified metric
(7.12) is symmetric positive definite on and equivalent to in the sense that there exist constants (independent of ) such that
(7.13) -
2.
Defining the modified quadratic functional
(7.14) one has, along solutions of the full UIH evolution , the differential inequality
(7.15) In particular,
(7.16) for all and some constant depending only on the metric equivalence constants in (7.13).
-
3.
The hypocoercive rate can be bounded from below in terms of the Fisher gap and the UIH operator norms. More precisely, there is an explicit continuous function of the dimensionless ratios and such that
(7.17) In particular, as long as the UIH data stay in a compact set of such triples with fixed Fisher gap and bounded , , the ratio admits a strictly positive uniform lower bound.
Proof.
We sketch the main steps, since the argument is a finite dimensional variant of the standard hypocoercivity construction, written in the UIH metric language.
Existence and uniqueness of a solution of the Sylvester equation (7.11) follow from the strict negativity of with respect to . Indeed, in the orthonormal basis that diagonalises , the operator has real eigenvalues in , so and the Sylvester map is invertible. Standard bounds on Sylvester equations then give an estimate of in terms of , and .
For small enough, the matrix remains positive definite and equivalent to , with constants as in (7.13). This uses only the fact that is strictly positive definite and that is finite; the bound on is quantitative in these norms.
The time derivative of along the full UIH flow is
| (7.18) |
where we have introduced the symmetric matrix
| (7.19) |
Inserting and and using , , the Sylvester equation (7.11) is designed so that all terms linear in cancel in the symmetric part. Concretely one finds
| (7.20) |
where is symmetric negative definite, with spectrum bounded above by , while and involve only commutators such as and higher order combinations. Bounds on and in terms of , and follow from the definitions.
By choosing small enough, depending only on the dimensionless ratios , , one ensures that the perturbative terms in do not destroy the strict negativity inherited from . In other words, there exists such that
| (7.21) |
which is exactly the differential inequality (7.15). The lower bound (7.17) arises from tracking the dependence of the perturbative estimates on the UIH operator norms and the Fisher gap. Finally, applying Grönwall’s lemma gives exponential decay of , and the equivalence of and yields the corresponding decay in the original Fisher norm. This completes the proof. ∎
Remark (Relation to the dimensionless formulation).
Remark (Universality across UIH sectors).
The theorem applies verbatim to the reversible Markov chains, their Fokker-Planck limits and the GKLS models treated in this paper. In each case the metric is the Fisher or BKM metric at the stationary state, is the Fisher Dirichlet operator, and is the metric skew part of the generator. The Fisher gap is the smallest positive eigenvalue of the Dirichlet operator, and the operator norms , are directly computable from the finite dimensional matrices that appear in the code archive. The lower bound (7.17) therefore shows that once these UIH quantities are controlled, the exponential decay rate of the full non normal semigroup is uniformly bounded below by a fixed multiple of the Fisher gap, independently of microscopic details. In the IBM experiments of Section 6.4 the measured short time and long time BKM entropy decay slopes sit several times above the Fisher gap, consistently with this UIH hypocoercive picture.
8 UIH renormalisation group for dissipation
In this section we introduce a renormalisation group (RG) framework for the dissipative sector of universal information hydrodynamics. The aim is to understand how Fisher geometry, cost-entropy inequalities and hypocoercive decay rates behave under coarse-graining, and to identify universality classes of irreversible behaviour.
The key observation is that the entire dissipative structure of UIH is encoded in a symmetric, positive operator on a Fisher or BKM Hilbert space, together with an antisymmetric reversible piece. Coarse-graining can therefore be expressed purely at the level of these operators, via Galerkin projection. Once this is done carefully, the usual intuition from diffusion and Markov chains becomes precise: the Fisher gap scales as a second derivative under block RG, while reversible couplings become RG irrelevant at large scales in the class of local models we consider.
The UIH hypocoercivity theorem from Section 7 then propagates to large scales and yields a universal relation between hypocoercive rates and Fisher gaps.
8.1 Discrete UIH models and Fisher gaps
We begin with a finite dimensional formulation that unifies lattice diffusion, reversible Markov chains and linearised GKLS semigroups.
Let be a finite dimensional real Hilbert space. A discrete UIH system on is specified by a pair of linear operators with the following properties:
| (8.1) |
where the transpose is taken with respect to the inner product. The full generator is
| (8.2) |
The symmetric part encodes irreversible Fisher-Dirichlet dissipation and the antisymmetric part encodes reversible Hamiltonian or transport effects. In the density sector of lattice diffusion or Markov chains the data arise from the Fisher-Dirichlet operator, while in the linearised GKLS setting they arise from the BKM metric and the symmetric part of the GKLS generator.
We assume that there is a distinguished one dimensional subspace of constants corresponding to mass conservation and that . The dynamics is then considered on the orthogonal complement
| (8.3) |
where is strictly negative definite in the models of interest.
On we define the Fisher gap by
| (8.4) |
This is the smallest positive eigenvalue of the positive operator on . In the reversible case , the Fisher gap controls the asymptotic decay rate of quadratic entropies and relative entropy in the usual way. In the non normal case , the finite dimensional UIH hypocoercivity theorem proved in Section 7 shows that the true decay rate of the semigroup generated by satisfies a bound of the form
| (8.5) |
where is a positive function of appropriate dimensionless UIH couplings and built from and commutators . We will make these couplings explicit below.
8.2 Block renormalisation on the density sector
We now define a concrete RG step on the density sector for local lattice models. This will serve both as an illustrative example and as a building block for Markov and GKLS settings.
Consider a one dimensional periodic lattice with sites labelled by . Let be a positive mobility on the bond between sites and (indices taken modulo ). The discrete analogue of the continuum operator is the symmetric matrix defined by
| (8.6) |
with periodic indices. Constants lie in the kernel of and on the orthogonal complement of constants we have a positive operator with
| (8.7) |
The generator on density deviations is and the associated quadratic Fisher-Dirichlet form is
| (8.8) |
Fix a block size dividing and define . We index coarse blocks by , with block consisting of sites . We define two linear maps between fine and coarse density deviations:
| (8.9) |
The map averages over blocks and replicates coarse values within each block.
The coarse Fisher-Dirichlet form is defined by requiring that, for any coarse deviation , the irreversible power evaluated at the coarse scale equals the fine irreversible power of its uplift :
| (8.10) |
This uniquely defines the coarse operator as
| (8.11) |
The constant coarse vector lies in the kernel of , since maps it to the constant fine vector, and on the orthogonal complement we again have a strictly positive operator with a diffusive Fisher gap .
This construction is nothing more than Galerkin projection of the Fisher-Dirichlet operator onto the coarse subspace . It is the natural UIH RG step on the dissipative sector: standard inner products and Fisher-Dirichlet forms are preserved on the coarse modes by construction.
Numerically, for heterogeneous mobilities generated by smooth random fields on the lattice and block size , we observe that the diffusive Fisher gap scales as
| (8.12) |
with the ratio typically lying within a few percent of unity across ensembles of random heterogeneity fields. This matches the continuum intuition that a second derivative operator scales like under the map .
To compare irreversible behaviour across scales, we introduce a time rescaling that keeps the Fisher gap fixed. The rescaled coarse generator is
| (8.13) |
The linear dynamics on deviations is then at the fine scale and at the coarse scale, with coarse initial data .
For quadratic entropies of the form and , numerics show that the late time decay slopes satisfy
| (8.14) |
with quantitative agreement between the two rates at the level expected from finite dimensional effects. Thus, once time is measured in Fisher units, the asymptotic dissipation clock becomes invariant under the RG map on the dissipative sector.
8.3 Markov RG as a specialisation of UIH RG
We now show that the block RG defined above specialises naturally to reversible Markov chains and coincides with the coarse-graining implicit in the GKLS to Markov ladder.
Let be the generator of a finite reversible Markov chain on a state space with stationary distribution . On the space of functions we equip the inner product
| (8.15) |
The Dirichlet form is
| (8.16) |
Reversibility means that is self adjoint with respect to this inner product after conjugation by the diagonal matrix . In terms of densities around , the associated Fisher-Dirichlet operator is
| (8.17) |
exactly as in the GKLS to Markov unification in Section 6. The positive operator has a Fisher gap equal to the spectral gap of on the orthogonal complement of constants in .
Let be a surjection that partitions the fine state space into blocks labelled by coarse states . The coarse stationary distribution is
| (8.18) |
We define the coarse projection and inclusion by
| (8.19) |
The operator is the conditional expectation of onto the sigma algebra generated by the partition, and is an orthogonal projection with respect to . The map is its adjoint between and .
We then define the coarse generator by
| (8.20) |
A short computation shows that the coarse Dirichlet form satisfies
| (8.21) |
Moreover, row sums of vanish and is stationary and reversible for the coarse chain. Thus the Markov RG is exactly the same Galerkin projection of the Fisher-Dirichlet operator that we used for diffusion, now implemented in the geometry.
From the UIH viewpoint, the data form a reversible UIH model on the density sector. The RG step with
| (8.22) |
is precisely the UIH RG on that sector. The Fisher gap coincides with the Markov gap and scales under block RG in the same way as in the diffusion examples. After rescaling time to keep fixed, the asymptotic decay of Fisher information and relative entropy in the Markov chain is invariant under RG, which matches the asymptotic decay clocks observed in the finite dimensional Markov and Fokker-Planck tests in Section 3 and Appendix E.
8.4 Fisher-Lindblad RG in BKM geometry
We now lift the same construction to the linearised GKLS setting with BKM geometry. This yields a Fisher-Lindblad RG for dissipation that acts directly on the UIH data of the GKLS semigroup.
Let be a faithful stationary state of a GKLS semigroup on a finite dimensional Hilbert space . The tangent space of traceless Hermitian perturbations at can be identified with a real Hilbert space equipped with the BKM inner product
| (8.23) |
where is the BKM metric operator. The linearised GKLS generator acts on and splits into symmetric and antisymmetric parts
| (8.24) |
with
| (8.25) |
with respect to . The symmetric part defines a BKM Dirichlet form
| (8.26) |
and the Fisher gap is the smallest positive eigenvalue of on the orthogonal complement of the identity.
Suppose we are interested only in the dynamics of a finite set of coarse observables that span a subspace . Let be the BKM orthogonal projection onto and let be the inclusion. We define the coarse UIH data by Galerkin projection:
| (8.27) |
The BKM Dirichlet form on coarse observables satisfies
| (8.28) |
so quadratic irreversible power is preserved on the coarse sector. The Fisher gap is the smallest positive eigenvalue of on the orthogonal complement of the identity within . The coarse reversible generator is simply the restriction of the reversible transport to the coarse observables.
When we restrict to diagonal perturbations of in a basis where is diagonal, with inner product inherited from the BKM metric, the above construction reduces exactly to the Markov RG in Subsection 8.3. In this sense the Fisher-Lindblad RG unifies the Markov and GKLS sectors under a single UIH RG operation at the level of the linearised Fisher or BKM geometry.
The full microscopic GKLS semigroup need not close on a smaller Hilbert space after coarse-graining. For the purposes of dissipation, however, that closure is not necessary: the UIH hypocoercivity theorems and cost-entropy inequalities are formulated directly on the linearised Fisher geometry, and it is at that level that RG acts.
8.5 UIH couplings and hypocoercive RG flow
To quantify the effect of the reversible sector on dissipation we introduce dimensionless UIH couplings built from and commutators . On a given finite dimensional UIH model with Fisher gap , we define
| (8.29) |
where is the operator norm induced by the metric on (equivalently, the BKM Dirichlet norm in the GKLS examples). These couplings measure, in Fisher units, the strength of non-normality and the magnitude of the mixed commutator that controls how the reversible and irreversible sectors interact. With this convention, and coincide with the metric operator-norm couplings used in the UIH hypocoercivity Theorem of section 9, and in the IBM spectrometer invariants below.
The finite dimensional UIH hypocoercivity theorem from the previous section can be written schematically as
| (8.30) |
where is the spectral abscissa of the generator on , and is a positive function defined on a region of coupling space that depends on model assumptions. In particular, for bounded and within a compact set, has a positive lower bound.
Given a block RG step with block size , the Galerkin projections and define a coarse UIH model with Fisher gap and couplings
| (8.31) |
We then rescale time on the coarse model by a factor so that the reference Fisher gap is held fixed along the RG flow. On the rescaled coarse system, the relevant couplings to compare with the fine system are
| (8.32) |
A single RG step with time rescaling thus induces a map
| (8.33) |
in coupling space.
For local lattice models where is a discrete second derivative and is a discrete first derivative, dimension counting suggests that under the map , with time rescaling chosen to keep the diffusive Fisher gap invariant, the couplings should scale as
| (8.34) |
for some order one constant that depends on details of the discretisation. The diffusion and Markov numerics reported above confirm this scaling in detail for one dimensional heterogeneous models with both simple and random local antisymmetric .
Iterating the RG step one obtains a flow
| (8.35) |
where is the coupling map associated with a single block RG. In the local models tested, the flow drives rapidly towards the origin as increases. After finitely many steps the couplings fall into a small neighbourhood of that depends only on the block size and dimension, not on microscopic details of or .
Combining this with the hypocoercivity bound, we obtain an RG universality statement for UIH dissipation.
Proposition 8.1 (UIH hypocoercive RG universality, informal).
Consider a class of discrete UIH models on density or BKM sectors satisfying the following conditions:
-
1.
The symmetric part is uniformly elliptic and local, of diffusive type (second order on a lattice or graph), with Fisher gap .
-
2.
The antisymmetric part is local and of strictly lower differential order than (for example, a discrete first derivative or a finite range transport term).
-
3.
Under the block RG map with block size , followed by time rescaling to keep fixed, there exist constants depending only on the model class such that
(8.36) for all models in the class.
Then the RG flow drives the couplings into a compact set containing the origin, and the hypocoercive rates satisfy a uniform lower bound
| (8.37) |
with a constant depending only on the model class, not on microscopic details of or . In particular, after rescaling time so that is fixed, the asymptotic dissipation clocks of all models in the class are universally controlled by up to an order one prefactor.
In the one dimensional heterogeneous diffusion and reversible Markov models studied numerically, the assumptions of the proposition are borne out by explicit RG iterations. The Galerkin RG on the Fisher-Dirichlet operator realises the expected scaling of , while the couplings and contract strongly under RG, flowing towards the purely diffusive fixed point .
These explicit RG iterations are implemented by the script 50_uih_rg_coupling_flow_suite.py in Appendix A. For ensembles of random finite dimensional UIH models it constructs the Galerkin coarse grained generators, rescales time to keep fixed, and tracks the flow of , and as the RG step is iterated, providing a direct numerical realisation of the coupling-space map in Section 8.
8.6 Numerical evidence and universality classes
We briefly summarise the numerical evidence that supports the UIH RG picture.
For one dimensional heterogeneous diffusion models with random smooth mobility fields on the lattice, the block RG map with yields a coarse Fisher gap satisfying
| (8.38) |
across ensembles of random realisations. After rescaling time by , the late time decay slopes of quadratic entropies computed along fine and coarse trajectories agree to within a few percent, and are both close to , as expected for reversible diffusion.
For the same diffusion matrices, supplemented by various local antisymmetric , the UIH couplings and show strong contraction under RG. For simple advection type built from nearest neighbour discrete derivatives, as well as for random band antisymmetric matrices with finite range, a single RG step followed by time rescaling reduces by a factor of order and by a factor of order with an order one constant that depends only weakly on the realisation. Two or three iterations are sufficient to drive and close to zero in the models tested.
These results indicate that, in a broad class of local one dimensional models, the reversible sector is RG irrelevant in UIH variables once time is measured in Fisher units. The large scale hypocoercive decay rates thus become universal and are essentially determined by the Fisher gap of the dissipative sector alone.
From the UIH perspective, this identifies a diffusive Fisher universality class: models whose dissipative part is diffusive and local, and whose reversible part is of lower differential order and local, flow under RG to a fixed point with , at which the hypocoercive rate is equal to the Fisher gap up to a constant factor. All such models have the same large scale dissipation clocks when time is expressed in Fisher units.
More exotic universality classes arise when one or more of the assumptions above fail. If the reversible sector carries additional conservation laws or slow modes, if is of comparable differential order to , if there are topological obstructions in higher dimensions, or if extra dynamical fields such as fluxes in Cattaneo or telegraph type models are included, then the RG flow of may approach a non zero fixed point rather than the origin. In such cases the asymptotic relationship between and may carry nontrivial dependence on a small number of fixed point couplings. The Fisher-Lindblad RG introduced above provides a natural framework in which to analyse these possibilities.
Under this RG, the Fisher gap scales like a second derivative under spatial coarse-graining, while reversible UIH couplings built from and contract strongly in the diffusive class of local models. After rescaling time so that the Fisher gap is fixed, the asymptotic dissipation clocks become RG invariants and the hypocoercive decay rates are universally controlled by the Fisher gap. This identifies a diffusive Fisher universality class for UIH, and sets the stage for exploring more exotic universality classes in models where the reversible sector remains RG relevant.
8.7 Cattaneo Fisher regularisation as a UIH stress test
The Fisher-regularised Cattaneo model of Appendix C.4 provides a useful stress test for the UIH renormalisation group picture. At the PDE level we take a one dimensional Fisher diffusion with mass , Planck constant and an effective signal speed ,
| (8.39) |
and replace it by the hyperbolic regularisation
| (8.40) |
This is the standard Cattaneo or telegraph equation, with characteristic front speed
| (8.41) |
so that Fisher diffusion is endowed with a finite propagation scale set by .
Script 30_fisher_cattaneo_relativistic_speed_checks.py
implements (8.40) on a periodic grid and tracks the position of a sharp
front, verifying that the measured propagation speed matches at the few per cent
level. See Appendix C.4 for numerical details and parameter choices.
For UIH purposes the key point is that Eq 8.40 is not an arbitrary hyperbolic correction. It is the hyperbolic regularisation of an underlying Fisher metriplectic structure built from the same density , the same Fisher functional and the same Dirichlet operator as in the purely diffusive model. One convenient formulation introduces an auxiliary flux and writes the system in first order form as
| (8.42) |
with Fisher potential . Eliminating from (8.42) yields (8.40). The irreversible channel is still driven by the same Fisher gradient ; the new ingredient is the inertial term , which delays the response of the flux and enforces finite signal speed.
From the viewpoint of the density manifold the Cattaneo model therefore leaves the local Fisher geometry unchanged. The instantaneous entropy production rate and its curvature bound are still governed by the Fisher Dirichlet form
| (8.43) |
exactly as in the diffusive model. What changes is the short time relation between and : is no longer a pure Fisher gradient flow, but is mediated by a flux variable with its own inertial dynamics. The UIH question is whether this extra degree of freedom changes the large scale irreversible clock when we coarse grain.
To answer this, consider the diffusive block renormalisation of Section 8.2 applied to the Cattaneo equation. Introduce coarse grained variables
| (8.44) |
with diffusive dynamic exponent . Under this rescaling the derivatives transform as
| (8.45) |
Substituting into (8.40) and multiplying by gives
| (8.46) |
The inertial term carries an explicit factor . At fixed microscopic and , coarse graining to larger and larger spatial blocks corresponds to , so the rescaled equation flows towards
| (8.47) |
In other words, the hyperbolic correction is RG irrelevant in the diffusive UIH sense. The Cattaneo model and the Fisher diffusion live in the same UIH universality class for the coarse grained density: they share the same Fisher Dirichlet operator, the same late time entropy decay clock and the same diffusive scaling limit. The only difference is the finite signal cone at microscopic scales.
This analysis is fully consistent with the numerical role of the Cattaneo test.
Script 30_fisher_cattaneo_relativistic_speed_checks.py
is designed to probe the existence and scale of the Fisher light cone, not to
redefine the irreversible clock. The RG calculation above shows that once
coarse grained on scales with , the Cattaneo dynamics are
indistinguishable from Fisher diffusion at the level of the density sector.
The UIH renormalisation group therefore treats Fisher Cattaneo as a physically
motivated hyperbolic regularisation that leaves the diffusive Fisher universality class
intact.
8.8 Scope of the diffusive UIH class and anomalous directions
The renormalisation group construction in this section has been deliberately conservative. We have worked within a single density sector, with local, symmetric, second order Dirichlet forms and short range interactions, and we have imposed diffusive scaling , from the outset. Within this admissible class the combination of analytical symmetry arguments and numerical tests points to a single diffusive UIH universality class on densities.
It is natural to ask what kinds of “anomalous” behaviour lie outside this class, and how they might fit into a broader UIH programme. Two directions are worth highlighting as motivation for future work.
First, one can relax the locality assumptions on the Dirichlet form while keeping a Fisher-type metric on densities. In a translationally invariant setting an admissible nonlocal Fisher Dirichlet operator can be specified by a symbol in Fourier space, with the local Laplacian case corresponding to . Jump processes with heavy tails and Lévy flights naturally lead to fractional generators with symbols of the form
| (8.48) |
and associated fractional diffusion equations with dynamic exponent . From a UIH standpoint these models still admit a Fisher metric and a metriplectic split, but the admissible scaling laws and coarse graining maps are different. The diffusive blocking transformation used here suppresses any contribution at small ; a fractional universality class would require a distinct RG scheme tuned to and to nonlocal Dirichlet forms. We do not pursue such constructions in the present paper. The fractional direction is mentioned here only to emphasise that the uniqueness statements in the main text apply within the local, second order Fisher-Dirichlet class, and that genuinely new universality classes are expected once one admits long range jumps and anomalous scaling.
Second, one can enlarge the state space from a single density to a multi current hydrodynamic manifold carrying, for example, mass, momentum and energy densities. The reversible side then supports ballistic or sound like modes with dynamic exponent , and the dissipative sector is built from a collection of coupled Fisher metrics and Dirichlet forms. Even in such settings the density block of the dissipative operator is expected to obey the same UIH constraints as in the single component case, and the Fisher Dirichlet form on densities should still provide a coercive entropy decay floor. What changes is that the dominant relaxation of observables can be mediated by ballistic channels on intermediate scales, with the diffusive Fisher sector taking over only at very long times. The qutrit and Fokker-Planck asymptotic decay clock experiments already show that once one isolates a single conserved mode with a fixed Markov gap, the Fisher decay rate is universal across very different microscopic realisations. Extending this logic to full hydrodynamic closures is a natural next step, but lies beyond the single density scope of this first UIH paper.
Both of these directions suggest that the diffusive UIH universality class studied here should be seen as the base point of a larger landscape. Local Fisher diffusion with second order generators and short range interactions gives a unique, robust fixed point under diffusive coarse graining, and provides the natural density sector for Fisher-Lindblad unification. Fractional and multi current extensions require modified RG schemes and additional structure, and are left as programme items for future work.
8.9 Summary
The UIH renormalisation group for dissipation developed in this section ties together three strands.
First, at the level of definitions, the Fisher metric on densities and the canonical Dirichlet operator fix a natural class of irreversible generators. A block coarse graining and rescaling with diffusive exponent map admissible microscopic models back into the same class, with renormalised parameters. Within this single density, local, second order Fisher-Dirichlet sector the only stable fixed point is the diffusive Fisher class: short range Markov chains, their continuum Fokker-Planck limits and Fisher-Lindblad GKLS density sectors all flow towards the same coarse grained description.
Second, the code archive demonstrates that this picture is not an abstraction. Markov chains and Fokker-Planck models with matched gaps share the same entropy decay clock; finite GKLS generators reduce to canonical Fisher Dirichlet operators on densities; IBM hardware experiments realise the same Fisher geometry and irreversible clocks in genuine quantum devices; and the asymptotic decay clock tests confirm that the Markov spectral gap acts as a universal irreversible timescale across discrete and continuum realisations. The RG language provides a unifying interpretation of these results: it identifies the Fisher-Dirichlet diffusion as the universal fixed point of the density sector under admissible coarse graining.
Third, the Fisher Cattaneo test shows that hyperbolic regularisation with a finite propagation speed is compatible with this UIH picture. The Cattaneo system shares the same density level Fisher metric and Dirichlet operator as its diffusive counterpart. Under diffusive coarse graining the inertial term is suppressed by and flows to zero, so that Fisher Cattaneo and Fisher diffusion lie in the same UIH universality class on densities. The hyperbolic correction enforces a light cone at microscopic scales without altering the coarse grained irreversible clock.
Taken together, these ingredients support a simple but strong claim. Within the admissible class considered here, the combination of Fisher metric, Dirichlet form and diffusive RG defines a unique universal density sector for information hydrodynamics. The remaining sections and appendices embed this density sector into the Fisher-Lindblad GKLS framework, document the numerical and experimental tests in detail, and outline how more exotic sectors such as fractional and multi current hydrodynamics can be approached within the same UIH mindset.
9 UIH hypocoercivity and Fisher decay floors
The renormalisation group analysis in the previous section identifies a universal diffusive class for the density sector of admissible UIH models. Hypocoercivity provides the complementary structural statement at fixed scale: given a Fisher metric and a symmetric Dirichlet operator with spectral gap , how much can a reversible sector slow down irreversible relaxation? In this section we reformulate hypocoercivity in purely UIH terms and show that for finite dimensional metriplectic generators the large time decay rate of perturbations is uniformly bounded below by a positive multiple of , with a prefactor that depends only on two dimensionless UIH couplings. We then explain how the same structure extends to abstract Hilbert spaces, and how the finite dimensional tests in the code archive should be read as Galerkin approximations to that general statement.
9.1 Metric split and Fisher gap
Let be a finite dimensional real vector space equipped with a positive definite symmetric metric , represented in coordinates by a symmetric positive definite matrix (also denoted ). The associated inner product and norm are
| (9.1) |
The metric adjoint of a linear operator is the unique operator satisfying
| (9.2) |
so that in coordinates .
We consider a real generator split into metric symmetric and skew parts
| (9.3) |
By construction and . The case of interest is when is negative semidefinite in the metric inner product and encodes the dissipative Fisher-Dirichlet operator associated to some entropy functional, while is the metric skew part of the generator that preserves entropy to first order.
In the earlier complex-mobility discussion (3) we wrote for the complex Schrödinger/Lindblad generator acting on wavefunctions and density matrices. Here we work with the associated real metrised generator on perturbations, so the factor of is absorbed into the definition of and the split is understood in the metric-adjoint sense.
We assume that acts on a codimension one subspace of metric mean zero perturbations, orthogonal to a distinguished stationary direction (the constant density mode or stationary state). On the symmetric operator is strictly positive definite. The Fisher gap associated to is defined as
| (9.4) |
Equivalently, is the smallest positive eigenvalue of on . In typical applications is the Fisher-Dirichlet operator of a reversible Markov generator, a Fokker-Planck diffusion or a GKLS dissipator in BKM metric, and is the Markov gap driving Fisher entropy decay in the purely reversible case.
To quantify the size of the reversible sector we equip with the induced operator norm
| (9.5) |
and define
| (9.6) |
The two adimensional UIH couplings associated to are
| (9.7) |
They compare the strength of the reversible sector and its commutator with the Fisher-Dirichlet part to the fundamental Fisher gap . In reversible models and ; in nonreversible flows and are typically of order one.
The key hypocoercive question is: given , and , how small can the actual decay rate of perturbations under the full semigroup be? The finite dimensional UIH hypocoercivity theorem below shows that, under these metric assumptions, the large time decay rate is always bounded below by a positive multiple of that depends only on and .
This theoretical statement is complemented by a large-scale Monte Carlo scan over finite dimensional UIH models, implemented in the script 49_uih_hypocoercivity_coupling_scan.py in Appendix A. There we generate random triples with a single stationary mode, compute , and the dimensionless couplings , on the mean-zero subspace, and map out the distribution of across coupling space. The numerical results show that remains an order-one multiple of throughout the sampled regime, in line with the finite dimensional UIH hypocoercivity theorem.
Quantumness of the commutator coupling.
In the abstract hypocoercive setting above the coupling is defined purely in terms of the Fisher metric and the pair , and measures the size of the commutator in Fisher units. In concrete GKLS realisations, where is the BKM tangent space at a faithful stationary state , it is convenient to separate a canonical “classical” contribution from a genuinely quantum part.
Let denote the density sector spanned by diagonal perturbations in the eigenbasis of , and let be its BKM orthogonal complement, spanned by coherent off diagonal modes. Appendix E.6 and the GKLS tests in Appendix E show that the restriction of to coincides with the classical Fisher Dirichlet operator associated to an effective Markov generator on the stationary eigenvalues of . Denote by and the orthogonal projections onto and respectively, and write
We define the density sector commutator and its remainder by
The operator measures the non commutativity of and seen purely inside the classical Fisher shadow, while collects all contributions that involve coherent directions or density coherence mixing.
Using the same metric operator norm as in the definition of , we introduce dimensionless couplings
By construction one has for purely classical UIH models, such as finite Markov chains in Fisher coordinates, and for diagonal GKLS lifts where populations evolve autonomously. In genuinely coherent GKLS models the numerics of Appendix E show that is typically non zero, and reflects the strength of the coupling between dissipative and coherent degrees of freedom in the BKM geometry. In this sense can be viewed as an optional diagnostic of quantumness for UIH systems, while records the remaining classical contribution to the commutator coupling.
9.2 A finite dimensional UIH hypocoercivity theorem
We work on the codimension one subspace on which is strictly positive. The finite dimensional UIH hypocoercivity theorem proceeds by two steps: solving a Sylvester equation to construct a metric correction , and using it to define a deformed energy functional that decays at an improved rate.
Theorem 9.1 (Finite dimensional UIH hypocoercivity).
Let be as above, with , , , and suppose that on the mean zero subspace the Fisher gap . Let and be the metric operator norms of and , and let be the dimensionless couplings defined in (9.7).
Then there exists an explicit positive function and a constant such that for all mean zero initial data the solution satisfies
| (9.8) |
for all . Moreover:
-
1.
In the reversible case one has and , so that (9.8) reduces to .
-
2.
For in any bounded set there exist universal positive constants and such that and . In particular, for fixed and the hypocoercive decay rate is uniformly comparable to .
The proof uses a standard Villani-type hypocoercivity argument rewritten in UIH language [8]. On the reduced space the spectrum of lies in . Consider the Sylvester equation
| (9.9) |
for an unknown operator . The right hand side is metric symmetric and bounded, with norm controlled by . Because the spectra of and are separated by at least , the Sylvester map is invertible on . One can therefore solve (9.9) and obtain a unique symmetric with a bound of the form
| (9.10) |
where the implied constant is dimensionless.
For sufficiently small the deformed metric
| (9.11) |
is still positive definite and equivalent to , so that and are mutually bounded by constants depending only on and . Define the associated energy functional on mean zero perturbations
| (9.12) |
For any solution of one can differentiate with respect to time and, using the identities , and the defining equation (9.9) for , derive an estimate of the form
| (9.13) |
for a concrete function that can be written down explicitly in terms of the bounds on , and . Grönwall’s inequality then gives a uniform exponential decay bound in the deformed norm , and the equivalence of metrics translates this back into (9.8).
Theorem 9.1 thus expresses hypocoercivity entirely in UIH variables. The only inputs are the Fisher metric and gap of , the metric operator norm of the reversible sector and of its commutator with , and the associated adimensional couplings . All information about the coordinate representation, the underlying lattice or continuum, and the microscopic details of the model has been absorbed into these three quantities.
9.3 Markov, Fokker-Planck and GKLS examples
The abstract setting above is designed to encompass the finite dimensional models appearing in the code archive and appendices.
For reversible and nonreversible Markov chains on a finite state space with strictly positive stationary distribution , the natural metric is the classical Fisher metric on mean zero perturbations,
| (9.14) |
and the symmetric operator is the classical Fisher-Dirichlet operator , where is the Markov generator and its adjoint in the weighted inner product. The Fisher gap is then the usual Markov gap. Nonreversible Markov chains with the same and the same symmetrised generator share the same and but differ in the skew part . The finite dimensional theorem above shows that, provided and are controlled, the large time decay rate of Fisher information in nonreversible chains cannot be arbitrarily slower than that of their reversible counterparts. This is exactly the behaviour observed in the Markov decay clock experiments.
For finite dimensional Galerkin truncations of Fokker-Planck equations with smooth confining potentials, one obtains the same structure. Expanding the density in a finite orthonormal basis in , the Fisher metric is the quadratic form induced by the continuum Dirichlet form, is the truncated symmetric diffusion operator, and encodes the antisymmetric drift. Theorem 9.1 provides a rigorous finite dimensional bound on the gap between the continuum Fisher gap and the realised decay rate of entropy along the truncated flow. The numerical tests in the archive probe exactly this relation in concrete one dimensional and multi dimensional models.
For GKLS generators in finite dimension the relevant metric is the BKM metric at a faithful stationary state , acting on traceless Hermitian perturbations . The inner product reads
| (9.15) |
and coincides with the Fisher metric associated to quantum relative entropy. Expressing the GKLS generator in a real basis adapted to this metric produces a split with metric symmetric negative and metric skew. The finite dimensional UIH hypocoercivity theorem then gives a lower bound on the decay rate of BKM relative entropy in terms of the Fisher gap of and the couplings .
The IBM hardware experiments provide a direct illustration of this mechanism. Idle channel tomography on a noisy qubit, combined with BKM reconstruction, produces an effective K generator and its split into and . The eigenvalues of the metric symmetrised operator give a Fisher spectrum and, in particular, a Fisher gap . The measured decay rates of BKM relative entropy for a selection of initial states lie uniformly above a floor set by this gap, and the gap mode saturates the bound in the reversible limit. In the UIH language of Theorem 9.1, these hardware channels have moderate and , and hence a hypocoercive decay rate that is a controlled multiple of .
9.4 Abstract Hilbert space formulation outlook
The finite dimensional result above is structurally identical to the general hypocoercivity theorems of Villani and coauthors [8]. The UIH contribution is to express all constants and assumptions in terms of the Fisher metric, the Fisher-Dirichlet gap and the two dimensionless couplings built from and . This packaging admits a direct generalisation to infinite dimensional Hilbert spaces.
Let be a real Hilbert space with inner product induced by a strictly positive bounded metric operator , and let be a self adjoint negative operator with a spectral gap on the orthogonal complement of its kernel. Let be a skew adjoint operator that is relatively bounded with respect to , and suppose the commutator extends to a bounded operator on . Under these assumptions one can solve the operator Sylvester equation
| (9.16) |
on the reduced space by the convergent integral representation
| (9.17) |
and obtain a bounded symmetric correction whose norm is controlled by . For sufficiently small the deformed metric is equivalent to , and the same energy estimate as in the finite dimensional case leads to an abstract hypocoercivity bound of the form
| (9.18) |
for all in the mean zero subspace of . The constants and depend only on the dimensionless ratios and .
The continuum diffusion and GKLS models considered in this paper fit naturally into this abstract framework: Fokker-Planck generators with confining potentials and regular coefficients have self adjoint diffusion parts with spectral gaps in , and their antisymmetric drifts and commutators are bounded in the Dirichlet norm; finite dimensional GKLS semigroups with faithful stationary states yield bounded operators on the BKM Hilbert space, where the abstract theorem reduces to the finite dimensional one.
A full Hilbert space treatment with detailed functional analytic assumptions and proofs would take us too far afield here, and we refer to the hypocoercivity literature for such developments. For the purposes of the present UIH paper the finite dimensional Theorem 9.1, together with the Markov, Fokker-Planck and GKLS examples above, is sufficient to justify the interpretation of as a universal Fisher decay floor and to connect the numerical and experimental results to a general structural inequality.
9.5 Fisher gap as intrinsic time beyond finite dimension
The finite dimensional hypocoercivity theorem can be reformulated in a way that is independent of coordinates and dimension, and which makes the rôle of the Fisher gap as an intrinsic time unit explicit.
Let be a real separable Hilbert space equipped with a strictly positive definite metric operator . Let be a densely defined, selfadjoint, non positive operator on with a one dimensional kernel spanned by a stationary mode , and assume that has a spectral gap on the orthogonal complement
We interpret as the Fisher-Dirichlet operator of the irreversible sector and as its Fisher gap.
Let be a skew adjoint operator on with and such that and the commutator are bounded on . Writing
we obtain a linear UIH flow
with Fisher quadratic
Under the assumptions above one can adapt the abstract hypocoercivity estimates of Villani and coauthors to the present Fisher packaging. There exist positive constants and , depending only on the dimensionless couplings
such that for all and all one has the exponential decay bound
| (9.19) |
In particular the Fisher gap provides an intrinsic irreversible timescale: after the reparametrisation
the slow Fisher modes decay at least as fast as independently of the microscopic realisation of and or the choice of basis.
This abstract formulation covers the finite dimensional matrix models, Galerkin truncations of confining Fokker–Planck equations, and quantum GKLS generators in the BKM metric at a faithful stationary state. In each case the Fisher gap of the symmetric part of the generator defines a natural “entropy clock” which is stable under Fisher preserving renormalisation group steps: after coarse graining the generator and rescaling time by the ratio of Fisher gaps, the large time decay rate of remains of order one. In this sense the Fisher gap plays the rôle of an intrinsic time unit for UIH flows, beyond any specific choice of coordinates or finite dimensional representation.
10 UIH channel tomography and universality spectroscopy
The previous sections identified three structural layers of universal information hydrodynamics in finite dimensions: an algebraic layer where a UIH model is a real vector space of perturbations with metric and generator ; a spectral layer where the Fisher gap of on the traceless subspace and the hypocoercive decay scale of organise irreversible relaxation; and an RG layer where and the couplings built from and flow under coarse graining, with a diffusive Fisher universality class characterised by contraction of and a hypocoercive rate controlled by . In this section we show that, once a noisy quantum device admits process tomography for its idle and simple driven channels, these same UIH objects can be reconstructed experimentally. A single reconstructed channel with a full rank stationary state determines a finite dimensional UIH model in the BKM geometry at , together with a Fisher gap, hypocoercive decay scale and UIH couplings. Varying the device, idle depth, drive and coarse graining then turns -tomography into a concrete universality spectrometer for dissipation: different noise processes appear as points and flows in the plane.
Throughout this section we work at the level of linearised Fisher geometry and real Bloch coordinates, as in Section 6 and Appendix E. All algebraic identities are direct consequences of the metriplectic framework and the BKM geometry. All numerical statements refer to the finite dimensional GKLS tests of Appendix E. The hardware statements are experimental and are supported by the IBM K tomography and BKM curvature suite of Appendix E, together with the additional K tomography and semigroup scaling scripts listed in Appendix A.
10.1 From process tomography to UIH generators
Let be a trace preserving completely positive map on with . Suppose that admits a strictly positive stationary state and that we are given a tomographic reconstruction of in an operator basis adapted to . Concretely, for qubits and two qubit systems we use Pauli type bases and work with Bloch coordinates as in the Bloch GKLS examples of Appendix E. The tomographically reconstructed channel is then a real matrix such that
with the first component corresponding to the identity fixed, , and the remaining components encoding deviations from .
To interpret as a time slice of a continuous semigroup we perform a semigroup scaling test, following the finite dimensional analysis of Appendix B and the IBM idle depth experiments in Appendix E. One prepares families of circuits that implement idles of length , reconstructs the corresponding channels for several integers , and diagonalises each in a common eigenbasis when possible. Semigroup compatibility is tested by checking that the logarithms of the nontrivial eigenvalues scale approximately linearly in ,
with slopes that are independent of within experimental tolerances. The IBM idle experiments show that for modest depths the dominant eigenvalues lie close to a straight line in and that the deviations are consistent with sampling noise and small non Markovian corrections.
Once such a scaling regime is identified we fix a reference time step and define the effective real generator
on the traceless subspace of Bloch coordinates, using a principal branch of the logarithm near the identity. In practice we work directly at the level of the real Bloch representation and verify that the resulting generator has negative real parts in its spectrum, as in the GKLS and IBM K tomography scripts.
The BKM metric at is computed once and for all in the same operator basis. In the eigenbasis of the BKM metric is diagonal on the matrix units , with weights determined by the eigenvalues of ; this defines a positive diagonal matrix on the vectorised operator space. Transforming back to the chosen operator basis gives a positive definite metric matrix on Bloch coordinates. This procedure is implemented in detail in Appendix E for coherent qubit models and in the IBM BKM curvature experiments.
With and in hand we form the metric adjoint
and split the generator into its symmetric and antisymmetric parts in the BKM metric,
The algebraic identities of Appendix E and the GKLS examples there guarantee that, for GKLS dynamics linearised at a full rank stationary state, is strictly negative definite and symmetric in the BKM metric and is skew. The IBM K tomography experiments show that the same structure appears for a noisy idle channel on hardware: the reconstructed , , and satisfy the metric symmetry and skew symmetry conditions to within numerical residuals, and the stationary direction is isolated and orthogonal to the traceless subspace.
Thus a single channel with a full rank fixed point and a mild semigroup property on a suitable idle family determines a finite dimensional UIH quadruple in the sense of Section 7, with all objects reconstructed from experimental data by linear algebra.
10.2 Fisher gap, hypocoercive scale and UIH couplings
On the traceless BKM Hilbert space the dissipative block is symmetric and negative definite. The Fisher gap of the channel is defined exactly as in the abstract theory,
where the spectrum is taken on . The corresponding Fisher Dirichlet inequality
is an algebraic identity once and are computed, and provides an information theoretic curvature floor on the dissipative contraction in the BKM metric.
The full generator has spectrum contained in the left half plane, with a spectral abscissa
In the linearised regime a perturbation evolves according to , and the BKM relative entropy between and decays approximately as at late times. The finite dimensional UIH hypocoercivity theorem of Section 7 shows that there exists a positive function of a small number of operator norms such that
where the UIH couplings are
and the norms are induced by the BKM metric. In particular, when and lie inside a compact region of parameter space, has a strictly positive lower bound and the hypocoercive decay rate remains uniformly comparable to the Fisher gap.
For a reconstructed channel these quantities become experimentally accessible scalars. The IBM K tomography experiments of Appendix A.3 follow exactly this route. First, the Fisher gap of the dissipative block is extracted from the smallest positive eigenvalue of on the traceless Bloch space. Second, a catalogue of initial perturbations is evolved on hardware, the BKM relative entropy is reconstructed from tomographic snapshots and the late time slope of is fitted, giving an empirical that is independent of the initial condition within error bars. In all runs the measured decay rates lie above the Fisher gap, with the gap mode initial condition producing the closest approach to and more generic perturbations decaying faster. Third, the operator norms of and are evaluated numerically in the BKM metric, giving explicit values of and .
Taken together, these measurements implement a complete UIH diagnostic for a single channel: they identify a Fisher curvature floor , a hypocoercive decay scale , and a pair of dimensionless couplings that quantify the strength of the reversible sector in Fisher units. The analytical hypocoercivity theorem then provides a structural inequality between these experimentally derived scalars.
10.3 Semigroup scaling and RG flow from channels
The Fisher Lindblad RG of Section 8 acts directly on the Fisher or BKM geometry and on the pair by Galerkin projection onto coarse subspaces of observables, followed by a time rescaling that keeps the Fisher gap fixed. For a reconstructed channel this RG step can be implemented numerically at the level of the Bloch representation, without reference to a particular microscopic Lindblad dilation.
Let be a linear map that selects a coarse family of observables and let be its BKM adjoint, so that for an induced coarse metric . The RG step defines the coarse operators
as in Section 8.4. The Fisher Dirichlet form is preserved on coarse observables and the Fisher gap of is computed on . Time is then rescaled by the factor and the rescaled coarse couplings
define an RG map
in coupling space. In diffusive Fisher universality classes the RG flow contracts to the origin, while more exotic reversible sectors can generate nontrivial fixed points with .
On a finite device the most natural coarse observables are low weight operators and block averages of local observables. For a single qubit idle channel the only nontrivial coarse subspace is the full traceless Bloch space, so the RG step is trivial. For two qubit channels the coarse observables include the single qubit Pauli operators and a subset of two qubit correlators; the corresponding projectors and adjoints are constructed explicitly in the two qubit GKLS and IBM tests of Appendix E. For larger systems, where process tomography is restricted to a fixed operator basis, the same construction applies on the subspace reachable by the tomographic probes.
In parallel with spatial coarse graining, the semigroup scaling test provides a time coarse graining. Eigenvalues of at idle length define effective generators ; comparing the associated across probes how the UIH data transform under time blocking. In a Markovian regime the triples coincide up to statistical errors. Deviations indicate non Markovian memory or slow drift in the hardware noise model. In either case, the channel family produces a discrete RG trajectory in the UIH coupling space.
10.4 Hardware as a universality spectrometer
The IBM experiments reported in Appendix E realise all three layers of this construction for concrete noisy channels. For a single qubit idle channel, process tomography in the Pauli basis yields a real Bloch matrix with a unique full rank stationary state. The BKM metric at this state is reconstructed from the eigenvalues of and agrees with finite difference estimates of the Hessian of quantum relative entropy. The resulting metric matrix is well conditioned and positive definite. Using the semigroup scaling protocol on increasing idle depths one finds a clear linear regime for the logarithms of the dominant eigenvalues, allowing an effective generator to be defined.
Splitting in the BKM metric gives and with numerical residuals at the level of the GKLS tests in Appendix B. The Fisher gap of is strictly positive and the associated Fisher Dirichlet inequality holds mode by mode. A catalogue of initial perturbations is then implemented by preparing different input states and evolving under the idle channel. The BKM relative entropy to decays exponentially in all runs, with late time slopes that cluster above and are independent of the initial state within error bars. The gap mode initial condition, constructed by perturbing along an approximate Fisher eigenvector of , gives a decay rate that nearly saturates the Fisher floor.
For two qubit idles and simple driven channels, the same pipeline can be carried out in a larger Bloch space. The stationary state is again full rank, the BKM metric is well conditioned, and the reconstructed generator exhibits a clean split into a strictly dissipative symmetric part and a reversible part. Projecting to the density sector in the eigenbasis of reproduces an effective classical Markov generator with a Fisher Dirichlet operator that coincides with the density block of , as in the GKLS ensembles of Appendix B. The measured Fisher gaps and hypocoercive scales line up with the Markov spectral gaps to the same tolerances seen in the numerical suite.
In UIH language, each such channel provides a point in a low dimensional parameter space of universality data. Comparing different devices, drives, idle lengths and coarse observables maps out a cloud of points and RG trajectories in this space. The diffusive Fisher universality class identified in Section 8 corresponds to channels whose measured couplings contract under coarse graining and whose hypocoercive decay rates are essentially set by the Fisher gap. More exotic classes would manifest as channels whose remain of order unity under Information RG or flow to nontrivial fixed points, with hardware decoherence spectra that cannot be reduced to purely diffusive Fisher curvature.
We do not claim such exotic universality has been observed. The present hardware experiments simply show that, on contemporary noisy devices, the idle and simple driven channels that are accessible to full tomography fall firmly inside the diffusive Fisher regime and obey the UIH inequalities to the stated tolerances. What the construction provides is a systematic way to look for deviations: once a channel can be tomographically reconstructed and its BKM geometry computed, the associated UIH data can be measured and compared across platforms. In this sense a quantum processor equipped with process tomography is not only a computer but also a small universality spectrometer for dissipation, reading out the Fisher curvature, reversible couplings and hypocoercive scales of its own noise.
10.5 UIH spectrometer on IBM Quantum hardware
Track 3 tests the UIH picture on real quantum hardware by treating an IBM two qubit noisy channel as an unknown element of a dissipative universality class and then extracting its UIH invariants directly from tomography data. The aim is not to model the microscopic hardware in detail, but to ask whether the effective generator sits inside a UIH hypocoercive basin and whether a Fisher preserving renormalisation group flow drives it towards a stable slow sector with universal parameters.
10.5.1 Two qubit idle channel and BKM generator reconstruction
We consider a fixed two qubit device and prepare a family of tomographic circuits that implement an idle channel over a prescribed time interval followed by an informationally complete measurement. The circuits and calibration settings are recorded in JSON files, while the processed tomography results are stored as .npz archives. For each run we reconstruct a completely positive trace preserving map acting on density matrices and its real superoperator representation
| (10.1) |
in a Pauli like basis, normalised so that the identity is an eigenvector with eigenvalue one. Restricting to the dimensional traceless subspace produces a real matrix , while logarithmic regularisation with a small spectral shift yields an effective generator
| (10.2) |
with the physical idle time and the regularised channel. This acts on the traceless space in an orthonormal BKM basis at the reconstructed stationary state , so that the BKM metric becomes the identity and Fisher geometry is Euclidean to first order.
The resulting archives store, for each tomography run, the regularised generator on the dimensional BKM orthonormal traceless space, together with the underlying stationary state and metric data used to construct it. For the purposes of the UIH spectrometer, we treat as the unique effective GKLS generator for that run.
The numerical UIH spectrometer built on these archives is implemented in the IBM
toolchain script 48_ibmq_uih_spectrometer_suite.py in the code archive
(Appendix A), which reads the *uih_split_bkm.npz files,
extracts in the BKM orthonormal traceless basis and computes the
associated UIH invariants and slow-mode RG data.
10.5.2 UIH invariants for the IBM channel
Given on the BKM orthonormal traceless space, the symmetric and skew parts
| (10.3) |
play the role of the Fisher Dirichlet operator and reversible generator respectively. The Fisher gap of the IBM channel is defined as the smallest strictly positive eigenvalue of on the traceless space,
| (10.4) |
while the hypocoercive rate of the full generator is
| (10.5) |
The strength of the reversible sector and its noncommutativity with the dissipative sector are captured by the operator norms
| (10.6) |
and the associated dimensionless UIH couplings
| (10.7) |
These quantities are invariant under similarity transformations that preserve the BKM orthonormal structure and are therefore intrinsic to the IBM channel at the reconstructed stationary state.
We apply this spectrometer to three tomography runs of the same physical idle channel, with effective shot counts of , and . The resulting invariants are summarised. The Fisher gaps sit consistently in the few milli range, while remains of the same order but fluctuates between approximately half and twice across runs. The UIH couplings indicate a strongly hypocoercive regime: the reversible part has operator norm of order but lies in the high tens, and the commutator norm produces in the hundreds to low thousands.
| run | |||||
|---|---|---|---|---|---|
These measurements place the IBM channel firmly in a weakly dissipative but strongly hypocoercive regime. The Fisher gap sets a diffusive entropy clock of order a few hundred idle time units, while the reversible part is tens of times stronger than the dissipative drift in Fisher units and substantially noncommuting with it. The variation of across runs is consistent with the expectation that small reconstruction differences in can significantly change how and interfere on the slowest modes when the couplings and are large.
These scalar invariants are extracted by the script
48_ibmq_uih_spectrometer_suite.py in Appendix A,
applied to the reconstructed *uih_split_bkm.npz files for each run.
10.5.3 Slow mode renormalisation and fixed point behaviour
To connect these hardware measurements to the UIH renormalisation group picture, we define a simple Fisher preserving coarse graining scheme based on the slow modes of . For each run we diagonalise on the traceless BKM space and select the eigenvectors associated to the four smallest strictly positive eigenvalues. Writing these orthonormal eigenvectors as the columns of a matrix , the coarse generator in the slow Fisher sector is
| (10.8) |
We compute the Fisher gap and UIH invariants of using the same formulas, then rescale by a scalar factor so that its Fisher gap matches the microscopic . This defines a single RG step that preserves the entropy clock but integrates out fast Fisher modes. The same procedure can be iterated.
In practice this slow-mode renormalisation is implemented by the spectrometer, rg and rg2 commands of the IBM UIH spectrometer script 48_ibmq_uih_spectrometer_suite.py in Appendix A, which apply the same invariant extraction and Fisher-gap matching procedure directly to the reconstructed matrices.
For the IBM channel, a single step with slow modes has two striking effects. First, the Fisher gap remains fixed by construction but moves into a narrow band, with ratios
| (10.9) |
for the , and runs respectively. Despite the raw spread in at the microscopic level, the slow sector decays at approximately two to three times the Fisher gap in all three experiments. Second, the UIH couplings undergo a strong contraction under RG. For the same runs one finds
| (10.10) |
and
| (10.11) |
The slow Fisher manifold therefore flows from a strongly hypocoercive microscopic regime, with large couplings, into a moderate coupling basin where both and are of order one to ten and is stable across runs.
Applying a second RG step to the coarse slow sector, using the same four mode scheme and Fisher matching, leaves the invariants unchanged to numerical precision: the Fisher gaps, hypocoercive rates and couplings at step two coincide with those at step one. In other words, the four dimensional slow Fisher sector of this IBM channel is a fixed point of the chosen renormalisation group transformation. The flow
| (10.12) |
contracts the UIH couplings from large values into a small basin and then stabilises, while the ratio converges to an experiment defined constant of order two to three.
The pattern of strong coupling contraction and stabilisation of observed here mirrors the synthetic UIH RG flows generated by the ensemble script 50_uih_rg_coupling_flow_suite.py in Appendix A, in which large families of random finite dimensional UIH models are coarse grained along their Fisher-slow modes and exhibit the same approach to a moderate-coupling basin with an order-one hypocoercive constant.
From the UIH perspective this provides an operational demonstration of a dissipative universality class on real hardware. The IBM idle channel is weakly diffusive at the Fisher level but strongly hypocoercive microscopically. Under Fisher preserving coarse graining along the slow modes of , it flows into a four dimensional slow sector with a stable hypocoercive constant and moderate couplings, and this sector is invariant under further RG steps. The hardware thus realises, within experimental uncertainty, the type of UIH renormalisation group structure predicted by the theoretical analysis in Sections 4 and 6, with a measurable universality constant
| (10.13) |
for the slow dissipative dynamics of the reconstructed two qubit channel.
10.6 Summary
The UIH hypocoercivity framework developed in this section can be summarised as follows. Given a Fisher metric on perturbations around a stationary state, the symmetric Fisher-Dirichlet operator and its gap determine the irreversible clock in the reversible limit. Introducing a reversible sector deforms the dynamics without changing the metric. Theorem 9.1 shows that, on finite dimensional shells, the large time decay rate of perturbations under the full generator is bounded below by a positive multiple of , with a prefactor depending only on the dimensionless UIH couplings and . In particular, in all examples studied in this paper the Fisher gap provides a robust lower bound on entropy decay rates, and the reversible sector cannot slow relaxation by more than an order one factor controlled by and .
This hypocoercive structure is realised concretely in reversible and nonreversible Markov chains, in finite dimensional Galerkin truncations of Fokker-Planck equations, and in GKLS generators in BKM geometry. The IBM hardware experiments show that the same mechanism survives in genuine quantum devices: the smallest positive eigenvalue of the BKM Dirichlet operator inferred from idle channel tomography acts as a Fisher decay floor for BKM relative entropy on hardware, and the gap mode saturates this bound in the reversible limit. The abstract Hilbert space formulation indicates that these finite dimensional observations are not artefacts of truncation, but manifestations of a general UIH hypocoercivity inequality for Fisher-Dirichlet generators with bounded reversible couplings.
In macroscopic scalar reductions such as the Fisher halo models of Ref. [dunkley2025gravity], the same operator data enter through an effective Fisher temperature that weights a bounded entropy channel, so that, in simple two-scale toy constructions, can be treated as a dimensionless parameter built from the hypocoercive couplings and the Fisher gap, providing a natural way to organise cusp–core behaviour and susceptibility trends within the same UIH hypocoercivity framework. A detailed microscopic derivation of this mapping is left for future work.
11 Fisher-Jarzynski Fluctuation Structure and RG Delta-Free-Energy Tests
The reversible-dissipative decomposition of a UIH flow, endows any one-dimensional Fisher-diffusion model with a natural family of free energies parameterised by a control variable entering the potential or mobility. This setting is sufficiently rich to realise a full Fisher analogue of the Jarzynski relation and to test its compatibility with the RG coarse-graining map developed earlier in the paper.
11.1 Fisher free energy and the work functional
Let be the stationary density of the Fokker-Planck operator and define the Fisher free energy
normalised so that the Fisher entropy and the potential energy are treated on equal footing. If is a smooth ramp from to , a trajectory naturally carries a Fisher-thermodynamic work functional
Following the usual argument for diffusion processes, but now with Fisher mobility, one obtains a Jarzynski identity of the form
| (11.1) |
valid for arbitrary finite-time ramps. No assumptions beyond the basic UIH regularity hypotheses (symmetry and positivity of , normalisability of , and standard ellipticity) are required.
11.2 RG invariance of the free energy difference
A key question for UIH is whether the free energy difference, which is a strictly microscopic object, is preserved under the Gaussian UIH coarse-graining operator . This maps densities to densities by convolution with a Gaussian of width , inducing a family of coarse-grained stationary states . Define the corresponding coarse-grained free energies in the same Fisher form.
For small the RG fixed point calculation gives
that is, the free energy difference is invariant up to fourth order in the coarse-graining scale. The absence of quadratic corrections follows from the symmetry of the Gaussian kernel and the fact that all terms integrate to total derivatives in the Fisher metric.
11.3 Numerical validation
The script 53_uih_rg_fluctuation_suite.py in the code archive implements this test for a canonical one-dimensional Fisher diffusion with a quartic potential. For a ramp , the microscopic free energy change is . After coarse graining the stationary densities with kernels of width , the numerical free energy differences satisfy
with a log-log slope of approximately , in excellent agreement with the theoretical prediction of quartic suppression of RG error. The same script performs an independent UIH-Jarzynski check by sampling trajectories under a time-dependent potential. The estimate obtained from (11.1) differs from by less than , consistent with finite sample variance of the exponential work average.
11.4 Interpretation
These tests support a central structural claim of the UIH framework: for Fisher diffusions, the free energy difference behaves as an RG invariant and encodes the leading thermodynamic information of the flow. The Jarzynski identity provides a stochastic realisation of the same quantity, and both characterisations agree with high numerical accuracy. This validates the use of as a Fisher appropriate observable for UIH renormalisation and prepares the ground for the full metriplectic RG in the quantum and Markov settings.
12 Kähler structure and holomorphic RG for finite dimensional K flows
The operator has so far been treated in its information geometric form, where induces a Fisher metric and encodes the Hamiltonian part of the flow. In finite dimensions this structure can be sharpened into an explicit Kähler triple, and the UIH coarse graining map can be tested for holomorphicity in a clean toy model.
12.1 A canonical Kähler triple on a lattice K flow
Consider a one dimensional periodic lattice with sites and real tangent space coordinatised as The complex structure is taken to act on each site by a ninety degree rotation
so that in matrix form with A site dependent positive metric is chosen diagonal in this basis,
so that commutes with and defines a Kähler metric. The corresponding symplectic form is then
which is automatically antisymmetric and nondegenerate. The triple satisfies and for all tangent vectors so it realises a standard Kähler structure on the finite dimensional K flow state space.
12.2 Gaussian coarse graining and induced Kähler data
The UIH coarse graining in this setting is built from a Gaussian convolution operator on the lattice. Let denote the matrix that smooths a lattice field by convolution with a Gaussian of width on the periodic chain. The full coarse graining operator on is then
which mixes neighbouring sites but acts identically on the and components. This commuting action of with preserves the complex structure at the linear level.
The coarse grained Kähler data are defined by pushforward of the metric and symplectic form,
and the effective complex structure is then obtained as the unique solution of In other words, we reconstruct from the requirement that the Kähler compatibility relation holds at the coarse grained level.
12.3 Holomorphicity diagnostics
The script 54_uih_kahler_rg_suite.py implements this construction on a modest lattice for a range of scales For each it computes and evaluates three diagnostic norms:
where denotes the Frobenius norm. These errors measure, respectively, the defect of from a genuine complex structure, the defect from Kähler compatibility, and the failure of the RG map to act holomorphically.
In the translation invariant model defined above the numerics show that all three quantities remain at the level of floating point roundoff for every tested :
within numerical noise. This is consistent with the analytic observation that commutes exactly with and that the site dependent weights are copied identically into the and slots, so the Kähler triple is preserved by construction.
12.4 Role in the UIH RG picture
This finite dimensional test serves two purposes. First, it demonstrates that the UIH K flow admits a natural Kähler geometric formulation even in simple lattice models, with the Fisher metric and Hamiltonian structure packaged into a single compatible triple Second, it shows that in a controlled setting the Gaussian UIH coarse graining can be exactly holomorphic, preserving the complex structure and the associated symplectic form at all scales.
In more general UIH applications, such as the quantum and Markov examples treated later in this section, one expects to be small but nonzero, decaying at least quadratically in for small coarse graining scales. The lattice calculation implemented in 54_uih_kahler_rg_suite.py provides a clean baseline in which the holomorphic RG limit is realised exactly, against which these more physical deviations can be measured.
13 Tomography–level UIH spectrometry on IBM Quantum hardware
The UIH formalism provides a natural set of operator diagnostics that can be applied directly to reconstructed generators from experiment. In this subsection we apply the cross–coherence machinery of Section 14 to dynamical generators obtained from two–qubit quantum process tomography carried out on the IBM Fez backend. The aim is not to test any particular physical model, but to assess whether the experimentally reconstructed behaves as a coherent UIH operator under the same coarse–graining probes that stabilise the classical, Markov, and GKLS examples elsewhere in the paper.
13.1 Experimental generators and the UIH split
Each tomographic run provides an real Liouville–space generator with after removing the trace component. Depending on the analysis pipeline, several versions of the same dataset are available: a raw reconstruction (4kshots.npz, 6kshots.npz, 8kshots.npz), a CPTP–projected version, and a UIH–motivated – splitting that produces arrays in the _uih_split_bkm.npz files. For the present test only the real part of is required; the BKM metric and symplectic form become relevant in later sections.
To compare datasets and ensure reproducibility, all numerical results here use the script 56_uih_tomography_spectrometer.py from the code archive. Given a chosen tomography file, the script:
-
1.
loads a generator key (‘‘K’’, ‘‘K_reg’’, ‘‘K_cptp’’, ‘‘K_real’’, etc.);
-
2.
simulates a synthetic state trajectory using the eigenbasis of ;
-
3.
builds beta maps by regressing against the full closure and the diagonal baseline ;
-
4.
evaluates cross–coherence under Gaussian index–space smoothing with
13.2 Beta maps: evidence for a stable UIH closure
Across all datasets with sufficient shot count, the beta maps exhibit a remarkably stable pattern. For the 6kshots_uih_split_bkm and 8kshots_uih_split_bkm files we obtain:
with the small scatter in consistent with numerical floating point noise. This means along the hardware generator is almost exactly described by the full closure , with negligible contribution from any diagonal–only approximation.
At lower shot counts (e.g. 4kshots_uih_split_bkm) the same qualitative structure persists, but the diagonal baseline becomes less suppressed, at the level of in norm. This increases the volatility of and produces certain pathological sign–flip effects under strong smoothing. Importantly, this does not affect the behaviour of the full closure: remains pinned at unity with vanishing variance.
13.3 Cross–coherence: scale stability of the hardware generator
The cross–coherence curves reinforce the same picture. For high–shot datasets we find:
which shows that the full closure retains a coherent activation profile under moderate smoothing in the 15–dimensional index space. The modest dip at intermediate scales is familiar from the synthetic tests and reflects the competition between initial activation and the subsequent mixing of components.
By contrast, carries no dynamical information. For the stable datasets it remains close to zero in absolute value for , while for the noisiest dataset it alternates between as a trivial consequence of the diagonal beta map being dominated by only a few components. This provides a clean separation between genuine operator structure (full closure) and an intentionally weak baseline.
13.4 Interpretation and UIH significance
These tomographic tests illustrate that a reconstructed superconducting qubit generator passes the same UIH closure diagnostics that validate the model in classical, Markov, and GKLS settings. The full operator structure is picked up with unit weight across all Bloch components, and its coherence is stable under moderate RG–style coarse graining. The diagonal baseline serves as a useful null model whose behaviour differs by orders of magnitude.
Two conclusions follow. First, UIH operator analysis extends naturally to experimental quantum data, with no modification to the core machinery. Second, the holomorphic and Fisher–geometric structures treated earlier in the paper have practical signatures: in a real hardware reconstruction the dominant coherent closure is the full off–diagonal K channel, precisely as predicted by the UIH metriplectic picture.
The script 56_uih_tomography_spectrometer.py and the accompanying figures and data (ibmq_beta_maps.npz, ibmq_cross_coherence.npz) are provided in the code archive to support full reproducibility.
14 Cross–coherence diagnostics for UIH operator closures
The UIH framework treats a dynamical generator not only as an operator but as a closure for the time derivative of a state vector or probability density. In finite dimensions this closure can be tested directly by treating the index of the state as a synthetic spatial coordinate and asking whether the measured time derivative is stably approximated by the operator action under controlled smoothing. The resulting “cross–coherence” is a highly efficient diagnostic for UIH compatibility, and is particularly useful when comparing alternative approximations or reconstructions of a generator.
14.1 Closure maps and regression
Let be a time–resolved state (for example the time series of Bloch components obtained from a simulated or experimental K flow). Given a fixed generator we define two operator closures acting indexwise:
The diagonal closure serves as a null baseline: it retains the correct scaling dimensions but discards all off–diagonal structure. A direct least–squares regression for each index gives coefficients
yielding two “beta maps” and that quantify how much each closure participates in the true derivative structure. For a well–behaved UIH generator we expect with small scatter and to be negligible.
14.2 Index–space coarse graining
To probe the robustness of these maps we apply Gaussian coarse graining in the index direction:
where denotes the discrete Fourier transform over the index variable. This smoothing mixes neighbouring components in a controlled way and introduces a scale analogous to spatial coarse graining in field theories. Repeating the beta–map regression using gives coarse grained beta maps and
14.3 Cross–coherence
The stability of the closure across scales is quantified by the Pearson correlation across index space:
A generator whose structure is genuinely encoded in the off–diagonal couplings should satisfy remaining high for modest , while the diagonal baseline typically loses coherence or exhibits trivial sign–flip behaviour due to its lack of spatial structure.
14.4 Numerical example
The script 55_uih_cross_coherence_spectrometer.py implements this diagnostic on a synthetic UIH Fokker–Planck evolution, with a time series of snapshots on a -point grid. The full closure is recovered with and a sharply peaked activation in the operator spectrum. The cross–coherence behaves as expected:
showing that the full operator structure survives smoothing while the diagonal approximation does not.
14.5 Role in UIH
This diagnostic complements the Fisher–RG and Kähler–RG results by testing a different aspect of operator universality: rather than focusing on free energy or geometric compatibility, the cross–coherence measures the stability of the dynamical closure itself. In contexts where the generator is reconstructed from data (for example quantum process tomography, Markovian inference, or nonparametric hydrodynamic fits), the cross–coherence spectrometer provides a one–page, model–agnostic test that the operator extracted from experiment has the correct UIH–compatible structure and does not collapse to a trivial diagonal form.
The same methodology will be applied in the following subsection to the IBM Quantum tomography data, where the reconstructed K flow passes this stability test with high precision.
14.6 Holomorphic Kähler RG from IBM tomography
The IBM two qubit tomography data also allow a direct test of the Kähler structure and holomorphic RG picture for a real hardware generator. Starting from the BKM orthonormal Bloch representation described above, the npz files ibmq_two_qubit_k_tomography_ibm_fez_20251119_225707_4kshots_uih_split_bkm.npz, ibmq_two_qubit_k_tomography_ibm_fez_20251120_090502_6kshots_uih_split_bkm.npz and ibmq_two_qubit_k_tomography_ibm_fez_20251120_094042_8kshots_uih_split_bkm.npz provide a real generator on the traceless Bloch sector of dimension . We split
in the BKM orthonormal metric. The symmetric part is expected to be negative semidefinite, with a single zero mode corresponding to the stationary state.
We first remove this stationary Fisher mode. Diagonalising and identifying the eigenvector with eigenvalue closest to zero, we construct an orthonormal basis for and restrict
to obtain a reduced generator on with . The Fisher metric on this reduced space is defined as
with a small eigenvalue floor applied if needed to ensure positive definiteness. For all three datasets the spectrum of exhibits three very small eigenvalues of order created by the regularisation and an active band of eigenvalues between approximately and .
We therefore define a Fisher active subspace by diagonalising and selecting those eigenvectors whose eigenvalues satisfy with . Writing for the diagonal matrix of active eigenvalues and for the restriction of to the corresponding block, we obtain
The active block is symmetrised and diagonalised as , and we form the inverse square root on the positive spectrum,
with a small cutoff to suppress numerical null directions. The metric polar factor
defines an approximate complex structure and symplectic form on . Embedding and as blocks in the full metric eigenbasis and rotating back by yields microscopic tensors and on the reduced dimensional space. On the active block we diagnose the Kähler quality via
with . For the three IBM Fez idle channels with k, k and k shots we find , and respectively, while is of order . Thus the microscopic generator supports an emergent Kähler triple on its Fisher active modes, and the Kähler defect decreases as the tomography statistics improve.
To probe holomorphicity under information RG we reuse the Gaussian index coarse graining introduced for the synthetic Kähler RG suite. On the reduced dimensional space this is a real matrix of the form
where is the discrete Fourier transform on the index space and runs over the discrete frequencies. For each scale we form
transform these to the metric eigenbasis, and restrict to the active block to obtain and . The coarse complex structure is reconstructed by enforcing Kähler compatibility on the active block,
implemented as a linear solve. We then evaluate three diagnostics,
| (14.1) | ||||
| (14.2) | ||||
| (14.3) |
where and denote the active blocks of and .
On all three datasets remains at numerical noise level by construction, while stays in the range to across the interval , with the k tomography run exhibiting the smallest microscopic defect. The holomorphicity defect is of order to at the smallest scales and grows to order by . This behaviour matches the UIH picture of a genuine microscopic Kähler complex structure on the Fisher active modes, which the Gaussian information coarse graining respects almost holomorphically at small before gradually breaking holomorphicity as more structure is integrated out. In particular, the combination of an emergent Kähler triple and a small band at short scales constitutes direct hardware evidence for the hidden complex structure underlying the one current two quadratures decomposition of the IBM idle channel.
15 Discussion and outlook
The three papers in this series describe what we have called an information hydrodynamics viewpoint on quantum theory and irreversible dynamics. The reversible companion paper shows that a Fisher information metric on densities, together with a canonical Poisson bracket on , singles out Schrödinger dynamics as the unique reversible hydrodynamics compatible with a small set of continuity and covariance axioms.
The entropy geometry paper shows that the same metric data support a metriplectic structure with cost-entropy inequalities and curvature coercivity, and that in simple settings a scalar “Fisher gravity” sector can be attached to density fields. The present work adds a concrete density sector: we exhibit finite dimensional GKLS generators, reversible Markov chains and Fokker-Planck limits whose dissipative dynamics are all realised by the same Fisher-Dirichlet operator on densities, prove a finite dimensional hypocoercivity theorem for the resulting generators , construct an information preserving renormalisation group picture for coarse graining, and show how these structures appear in quantum hardware experiments.
Read together, the three UIH papers can be viewed as different projections of a single complex -flow on pairs : the reversible Madelung work fixes the antisymmetric channel , the entropy geometry paper fixes the symmetric Fisher channel and its scalar gravity sector, and the present operator paper realises the full structure in finite dimensional, Markov, Fokker–Planck and GKLS settings.
It is therefore natural to speak of a single universal information hydrodynamics -flow, with the reversible, metriplectic and operator papers fixing its structure on and the Fisher halo gravity work providing a scalar gravitational channel built from the same Fisher part .
It is also useful to note that once a Fisher information metric has been fixed as the geometry on densities, the usual probabilistic weighting of observables is no longer an independent postulate but the natural dual pairing. In the classical setting, a tangent vector is paired with a potential by , so expectation values arise as the canonical duality between states and observables in the Fisher geometry. In the wavefunction representation this becomes , so the Born rule for position measurements can be viewed as a consequence of the chosen information geometry rather than an extra axiom.
UIH view of decoherence and measurement.
The operator picture developed above also suggests a simple language for the quantum to classical crossover.
In all of our examples the complex generator has a metric skew part that transports information without dissipation and a Fisher symmetric part that generates entropy production and drives the state toward its stationary distribution. Regimes in which dominates over in an appropriate operator norm are therefore nearly reversible: information functionals built from the BKM metric are approximately conserved and the dynamics is well approximated by unitary evolution. When dominates, the dynamics reduces to gradient flow of relative entropy and contracts states rapidly onto a small manifold of stationary or metastable configurations.
Coupling a microscopic system to a large environment can then be viewed, in UIH language, as a change in the effective generator in which the Fisher sector becomes large, so that off diagonal coherences in a suitable basis are quickly suppressed while the stationary state is preserved. In this sense decoherence and idealised measurement processes correspond to high regimes of the same -flow that governs coherent evolution, rather than requiring a separate dynamical postulate.
15.1 UIH as a natural geometry for dissipative quantum dynamics
The first broad conclusion is that the Fisher and BKM metrics and associated Dirichlet forms appear as very natural, and perhaps canonical, geometric data for dissipative quantum dynamics.
On the algebraic side, any finite GKLS semigroup with a faithful stationary state admits a BKM metric on observables and a Dirichlet form constructed from the dissipator and . In the classes of models treated here, this Dirichlet form coincides with the Fisher–Dirichlet operator obtained from the stationary distribution in the diagonal basis, both for diagonal GKLS generators and for coherent GKLS models whose density sector reduces to a reversible Markov chain. The continuum Fokker–Planck limits of these chains realise the same Fisher metric and Dirichlet operator at the level of densities, and the associated free energy functionals exhibit consistent decay behaviour.
A further structural feature is that, on suitable subspaces, the BKM metric and the Hamiltonian sector often admit a compatible complex structure. In the finite dimensional setting this gives a Kähler triple in which the symmetric part and skew part of are metric adjoints, the Hamiltonian dynamics are symplectic with respect to , and the complex structure relates the two. In the synthetic models studied here, and in the IBM experiments, this Kähler structure persists under Fisher preserving coarse graining, with only modest defects on the Fisher active modes. This suggests that Kähler geometry is not a special property of exactly solvable models, but an emergent organising structure for a broad class of dissipative quantum dynamics.
On the numerical side, ensembles of finite dimensional K flows and GKLS models exhibit the same splitting into a symmetric Dirichlet part and a skew part that is antiadjoint with respect to the BKM metric. The experiments on IBM hardware show that idle and driven channels reconstructed by process and state tomography admit a regularised generator whose symmetric part is well modelled by a Fisher–Dirichlet operator in the BKM geometry, while the skew part captures coherent dressing and non-normality. The Kähler diagnostics indicate that, on the slow subspace selected by Fisher geometry, the reconstructed generators are close to Kähler in the same sense as the synthetic models. These observations do not prove that every Markovian open quantum system admits such a realisation, but they strongly support the following working hypothesis.
For a broad class of GKLS generators with faithful stationary states, the density sector can be represented on a Fisher or BKM manifold in which the dissipative and reversible contributions are the symmetric and skew parts of a single operator , with the symmetric part realising a canonical Dirichlet form and, on suitable subspaces, admitting a Kähler structure together with the Hamiltonian sector.
15.2 Fisher gaps, hypocoercivity and limits on control
A central technical result of this paper is the hypocoercivity theorem for finite dimensional UIH generators. For a K flow with , where is symmetric negative semidefinite with spectral gap on the orthogonal complement of the stationary state and is BKM antiadjoint, one can bound the large time decay rate of the semigroup from below by
where depends only on dimensionless coupling parameters that measure the strength of relative to . In the ensembles studied here, and in the concrete GKLS and Markov models, the ratio remains of order one across a wide range of couplings.
In the GKLS–Markov–Fokker–Planck ladder, we further observed that the Fisher gap of the reversible Markov generator matches the late time decay rate of the corresponding Fokker–Planck generator, even when the curvature spectra of the two differ markedly. In the hardware experiments, effective Fisher gaps extracted from K tomography and BKM geometry control the observed relaxation rates of slow modes on IBM devices. In one dimensional diffusions, the same Fisher gap also governs the decay of free energy differences appearing in the Fisher–Jarzynski relation, so that the fluctuation structure and the linear relaxation rates are tied to a common geometric scale.
Taken together, these results support the view that the Fisher gap provides a natural decay clock for irreversible dynamics. Once the stationary geometry and Dirichlet form are fixed, the skew part can reshape trajectories, generate oscillatory transients and induce non-normal amplification, but it cannot indefinitely delay relaxation along density modes. Coherent driving can reduce relative to , but the hypocoercivity theorem and numerics indicate that for realistic couplings this reduction remains within a modest factor.
From a control perspective, this suggests that for fixed stationary state and microscopic couplings there are intrinsic limits to how much purely Hamiltonian control can suppress entropy production rates along given modes without leaving the UIH framework or introducing strong memory effects. Making such statements precise in concrete hardware models would require additional work, but the present results already indicate that UIH quantities such as the Fisher gap and coupling parameters, together with Kähler diagnostics on slow subspaces, function as natural invariants constraining control and decoherence engineering.
From a geometric point of view it is convenient to reparametrise time using the Fisher gap. The abstract hypocoercivity estimates of Section 9 show that for a UIH flow with one can bound the Fisher functional along the trajectory by
with positive constants and that depend only on the dimensionless couplings built from , and . Writing
this can be restated as
so that simply fixes the unit in which irreversible relaxation is measured. In this minimal sense we refer to the Fisher gap as defining an “entropy clock” or “irreversible clock”: different microscopic generators with the same stationary geometry and Fisher gap have comparable late time decay when expressed in the rescaled time variable , even if their reversible sectors and transient behaviour differ substantially. In the present paper we do not attempt to relate this Fisher clock to thermodynamic or relativistic notions of time; it is used only as an intrinsic parameter for comparing decay rates within the information geometric framework.
15.3 Universality, renormalisation and device characterisation
The renormalisation constructions in this paper are deliberately simple, but they already show several robust features. In the finite dimensional setting, Galerkin projections onto a chosen set of observables preserve the Fisher dissipation quadratic form exactly on those observables, and induce effective generators with the same UIH structure. In reversible Markov models and their Fokker–Planck limits, block decimation and rescaling steps produce flows in a small space of effective parameters that leave the Dirichlet form on slow modes invariant and tend to contract the couplings towards a diffusive regime.
The fluctuation layer shows a similar rigidity. In one dimensional Fisher diffusions, coarse graining procedures that preserve Fisher dissipation on slow observables also preserve free energy differences in the Fisher–Jarzynski relation, so that the same renormalisation group that organises gaps and couplings also leaves certain non-linear fluctuation quantities invariant. On the geometric side, the Kähler diagnostics indicate that the effective complex structure on Fisher active subspaces is stable under the same coarse graining operations, up to small and controlled defects.
The numerical GKLS examples and the IBM experiments are consistent with a simple picture in which a variety of microscopic generators flow under such coarse graining into a basin where the large time decay rate is controlled by the Fisher gap, Fisher–Jarzynski free energy differences are stable, and the Kähler triple remains approximately intact. We refer to this informally as a diffusive Fisher universality class. At this stage we do not claim a classification theorem, but the evidence suggests that Fisher regularised universality classes, characterised jointly by gaps, couplings, fluctuation free energies and Kähler data, may provide a useful organising language for non-equilibrium dynamics, in an analogous role to universality classes in equilibrium statistical mechanics.
Once this viewpoint is adopted, K tomography and UIH reconstruction become a form of universality spectroscopy for quantum devices. In the IBM experiments reported here, we reconstruct K on suitable subspaces, split it into symmetric and skew parts in the BKM metric, estimate Fisher gaps and hypocoercive couplings, apply Kähler diagnostics, and track cross coherence under coarse graining. The tomography data indicate that the idle and driven channels studied sit inside a diffusive Fisher basin with of order a few, with fluctuation and Kähler properties aligned with the synthetic ensembles. Cross coherence diagnostics distinguish full UIH closures from diagonal baselines and make the approach to the Fisher basin visible as a flow of correlation functions. It appears likely that similar UIH based characterisations can be extended to other platforms and to higher dimensional subspaces, providing a complementary figure of merit to conventional gate fidelities and noise parameters. The location of a device in Fisher coupling and Kähler diagnostic space, together with its effective Fisher gaps and fluctuation data, may offer a concise summary of its irreversible structure.
15.4 Outlook
We close by outlining a few directions that appear particularly natural in light of the present results.
Extended systems and open quantum field theory.
The hypocoercivity framework and the UIH splitting have been formulated in a way that depends only on metric adjoints and Dirichlet forms, not on finite dimensionality. Many generators used in open quantum field theory and kinetic theory have a similar structure: a selfadjoint dissipator with a spectral gap on suitable subspaces, together with a relatively bounded skew part implementing coherent dynamics. It would therefore be natural to investigate UIH realisations of open field theoretic models and to ask whether existing results on thermalisation and transport can be reformulated in Fisher geometric terms. In particular, one could seek analogues of the Fisher gap bounds, fluctuation relations and Kähler diagnostics in infinite dimensional settings, and study how they behave under spatial coarse graining.
Non-equilibrium universality beyond diffusion.
The examples studied here mostly fall into a diffusive Fisher basin. There exist, however, many systems with anomalous transport, glassy dynamics or strongly non-normal generators where relaxation is not well described by simple diffusion. The UIH renormalisation scheme provides a way to define effective generators and coupling flows for such systems. It would be interesting to identify models, either classical or quantum, whose coarse grained behaviour does not tend towards the diffusive Fisher regime and to examine whether they realise distinct UIH universality classes, potentially with different relationships between Fisher gaps, fluctuation free energies and Kähler structure. This could provide a common language for a range of non-equilibrium phenomena that are currently treated separately.
UIH structured model reduction and control design.
The Galerkin constructions in this paper show that, for a chosen set of observables, one can construct effective generators that preserve the Fisher dissipation quadratic form exactly on those observables. This suggests a systematic UIH theory of model reduction, in which reduced models are selected by their ability to preserve Fisher geometry, Dirichlet structure and, where relevant, Kähler data on the modes of interest. It also suggests a route to UIH guided control design, in which Fisher gaps, coupling parameters and Kähler diagnostics of are used as high level design targets for noise tailoring and dissipative state preparation. Developing these ideas would require combining the present linear theory with more detailed hardware specific modelling, but the underlying geometric structure is already in place.
UIH anomalies and the search for new behaviour.
Finally, it is worth noting that the UIH framework is constrained enough that its failure modes should be informative. In our view a useful long term goal is to formulate and search for UIH anomalies in experimental or numerical data: channels and dynamical regimes for which no reasonable choice of stationary state and metric makes the symmetric part of the reconstructed K into a Dirichlet form, for which hypocoercive estimates fail persistently even after refinements, for which Fisher–Jarzynski free energy differences are not stable under coarse graining, or for which no approximate Kähler structure can be found on any Fisher active subspace. One may also encounter cases where cross coherence never stabilises under renormalisation for any plausible closure, or where coarse graining destroys GKLS structure in ways not attributable to known non-Markovian effects. Such cases would indicate that some combination of the UIH assumptions breaks down at the scales being probed, and would help delineate the domain of validity of information hydrodynamics. Whether any such anomalies will be found remains to be seen, but the present work provides a concrete set of tools with which to look for them.
Summary
Information geometry, Kähler structure and metriplectic dynamics provide a natural backbone for describing irreversible quantum evolution, from finite dimensional models to real hardware. The specific realisations we have exhibited represent only the first steps in what appears to be a much wider programme. We expect that extending these constructions to more complex systems, and deliberately searching for their breakdown, will clarify both the scope and the limitations of the UIH picture.
Appendix
Appendix A Code Archive & Data
All numerical and experimental claims in this paper are backed by a reproducible Python archive, distributed as a single compressed file uih_archive.zip. The archive can be downloaded from
https://github.com/feuras/uhi_archive/ - DOI 10.5281/zenodo.17651171
Raw and processed data from IBM hardware experiments available from
https://zenodo.org/records/17672347 - DOI 10.5281/zenodo.17672346
and contains self-contained scripts with all parameters exposed near the top of each file, together with lightweight helper modules for common linear algebra, plotting and data management. The scripts are intended to be read as companions to the main text: each file implements the concrete model, geometry and diagnostics described in the corresponding section or appendix.
The archive is organised into three main families.
A.1 Core GKLS-Markov-Fokker-Planck chain
These scripts construct the basic UIH ladder from diagonal GKLS generators to reversible Markov chains and their Fokker-Planck limits.
- 01_gkls_diagonal_to_markov_checks.py
-
Implements the finite dimensional thermal GKLS jump model. Starting from diagonal Lindblad jump operators it constructs the Lindblad superoperator and the associated classical generator on the energy basis. It evolves both the density matrix and the Markov chain, comparing populations and the entropy curves and , and reports population and entropy discrepancies over the time grid. This numerically confirms that the diagonal sector of the GKLS semigroup coincides with the reversible chain and that the entropy decay curves match to numerical precision.
- 02_markov_to_fp_limit_checks.py
-
Implements the reversible lattice chains and Fokker-Planck limit of Section 5. For a sequence of lattice spacings it constructs reversible generators with Gibbs stationary law, integrates the master equation, and compares interpolated densities with a direct numerical solution of the limiting Fokker-Planck equation. The script estimates convergence rates in the spacing and compares discrete and continuum free energy decay , providing a concrete Markov-to-Fokker-Planck UIH ladder.
- 03_fp_fisher_metriplectic_checks.py
-
Works directly at the continuum PDE level for the overdamped Fokker-Planck equation of Section 5. It discretises the PDE on a fine spatial grid, computes numerical , the Fisher right-hand side , and the free energy functional along the evolution, checking that and within discretisation error. This confirms the Fisher-metriplectic structure at the PDE level.
- 04_gkls_coherence_elimination_checks.py
-
Numerical coherence-elimination test for a driven, damped qubit GKLS model. The script implements a coherently driven, strongly dephased qubit, integrates the full GKLS evolution in Bloch coordinates, and extracts the slow exponential tail of the excited-state population in the overdamped regime. It fits the late-time decay rate and compares it to the analytical coherence-elimination prediction , constructs the corresponding effective two-state Markov chain with rate , and overlays the Markov population curve with the GKLS dynamics. The diagnostics quantify the accuracy of the coherence-elimination approximation and identify the parameter range where the reduced Markov description is reliable.
- 05_markov_fp_free_energy_gap_checks.py
-
Free-energy gap and irreversible clock test for a reversible Markov chain and its Fokker-Planck limit. The script builds a simple reversible lattice Markov generator and its continuum Fokker-Planck approximation, computes the symmetrised generators and their spectral gaps, and evolves a large ensemble of random initial densities. For each run it tracks the free energy decay, fits late-time exponential rates, and compares them to the theoretical irreversible clock set by the smallest non-zero eigenvalue of the symmetrised Markov generator. The results illustrate that the Fisher curvature gap provides a coercive decay floor without necessarily fixing the dominant relaxation rate, and that the Markov and Fokker-Planck models share the same irreversible clock in the appropriate scaling regime.
- 30_fisher_cattaneo_relativistic_speed_checks.py
-
Implements the Fisher-regularised Cattaneo model. It solves the hyperbolic diffusion equation with Fisher regularisation on a periodic domain, tracks the propagation of a sharp front, and fits the front position to extract a measured signal speed . The script compares to the effective light speed , reporting relative errors at the few-percent level and confirming that the UIH Fisher regularisation supports a finite propagation speed.
- 31_uih_asymptotic_decay_clock_qutrit_fp_checks.py
-
Realises the asymptotic decay clock experiment. It constructs a reversible qutrit generator and a high-resolution reversible chain approximating a Fokker-Planck limit, computes the symmetrised generators and their spectral gaps, and evolves many random initial data. For each run it forms the instantaneous decay rate for the classical relative entropy and fits a late-time window, showing that in both models and that the Markov gap acts as a universal irreversible clock.
A.2 Finite dimensional flows and Fisher-Lindblad suite
These scripts implement the finite dimensional Fisher-Lindblad unification tests, built around real -generators, Fisher metrics and their symmetric and skew parts. They are designed to be run as a numerical suite.
- 06_gkls_fp_G_unification_checks.py
-
Constructs a discretised one-dimensional Fokker-Planck model, extracts the reversible Markov generator and its stationary distribution , and forms the canonical Fisher operator . It verifies that the irreversible drift can be written both as and as and checks the cost-entropy identity mode by mode, providing the canonical irreversible slice used throughout.
- 07_gkls_to_markov_G_unification_checks.py
-
Implements the diagonal GKLS jump model of Section 4 in matrix form and restricts the dissipator to the diagonal sector to obtain a classical generator with stationary law . It compares the resulting Fisher operator with the density block of the symmetric part of the real GKLS generator, confirming that the latter coincides with the canonical Fisher Dirichlet operator.
- 08_gkls_fp_nonrev_qutrit_checks.py
-
Builds nonreversible qutrit GKLS models whose density sector induces effective nonreversible Markov chains. It tests that the symmetric density block is still of Fisher-Dirichlet form and that the entropy and Fisher decays are controlled by the Markov spectral gap even away from detailed balance.
- 09_gkls_random_qubit_density_tests.py
- 10_gkls_random_qutrit_density_tests.py
- 11_gkls_random_qutrit_density_fp_limit_checks.py
-
Generate random ensembles of thermal and driven qubit and qutrit GKLS generators, compute their stationary states, and project to the density sector. They test, across many random samples, that the symmetric density block always matches a canonical Fisher operator of the form and that the global entropy decay follows the Markov gap while the Fisher gap supplies a curvature floor.
- 12_bloch_bkm_metric_checks.py
-
Computes the exact BKM metric on Bloch space for qubits, compares it to the discretised Fisher metric inferred from finite differences of relative entropy, and verifies numerically that the BKM metric is the local Hessian of quantum relative entropy at a chosen stationary state.
- 13_bloch_bkm_k_split_checks.py
-
Works with real Bloch generators in the BKM metric and performs the metric adjoint split , checking that is metric symmetric and metric skew at the level of numerical residuals, and that the Fisher production identity holds along the flow.
- 14_bloch_bkm_entropy_decay_checks.py
-
Evolves traceless Bloch vectors under both the full flow and the pure flow for random initial conditions, fits late-time decay rates of the BKM quadratic functional, and confirms that both rates scale with the smallest dissipative eigenvalue of , as predicted by the UIH irreversible clock picture.
- 15_qutrit_markov_vs_FP_universal_gap_checks.py
-
Constructs a reversible qutrit Markov chain and a high-resolution reversible chain approximating a Fokker-Planck limit, computes their gaps and tracks entropy decay, showing that both models share the same asymptotic decay rate once rescaled, reinforcing the universality of the reversible clock.
- 16_gkls_diagonal_to_markov_checks.py
-
Implements diagonal GKLS models with detailed balance and their mapping to reversible Markov chains. Starting from diagonal jump operators, it extracts the population generator from the GKLS superoperator and compares it to the constructed Markov generator . It evolves populations under both the GKLS semigroup and for random initial states and checks that off diagonal entries of remain negligible, with optional parallel ensembles.
- 17_gkls_fisher_dirichlet_checks.py
-
Tests Fisher-Dirichlet matching for diagonal GKLS models with detailed balance. For random reversible Markov models it builds a purely dissipative GKLS lift with jump operators , constructs the BKM metric at , forms the metric adjoint and symmetric part , and restricts to diagonal perturbations. It compares the classical Fisher-Dirichlet form with the density sector quadratic form for many mass conserving perturbations and reports maximal discrepancies.
- 18_gkls_coherent_density_sector_checks.py
-
Studies coherent GKLS models that share a fixed reversible Markov sector. The populations follow a reversible generator with stationary distribution , while diagonal Hamiltonians and dephasing operators generate nontrivial coherence dynamics that leave the populations and unchanged. The script checks that the population generator extracted from the full GKLS equals and that the Fisher-Dirichlet energy on the density sector, computed from the symmetric part in the BKM metric, matches the classical Fisher Dirichlet form. This probes density sector universality in the presence of coherent and dephasing dynamics.
- 19_gkls_nonrev_density_sector_checks.py
-
Extends the density sector Fisher geometry tests to non reversible Markov chains and diagonal GKLS lifts. The script constructs random non reversible Markov generators with column sums zero, computes the stationary distribution , and builds a diagonal GKLS generator with jump rates . It then forms the BKM metric, the metric adjoint and symmetric part , and compares the classical Fisher Dirichlet energy to the GKLS density sector Dirichlet for random mass conserving perturbations, testing that the Fisher-Dirichlet structure persists without detailed balance.
- 20_gkls_nonrev_decay_clock.py
-
Probes the Fisher decay clock for non reversible Markov chains and their GKLS lifts. For each random non reversible generator with stationary , it builds a diagonal GKLS model, constructs the symmetric Fisher Laplacian on the tilt variables with edge weights , and computes its spectral gap . It evolves both the Markov chain and the GKLS model from random diagonal initial states, tracks classical and quantum relative entropy and Dirichlet energies, and fits late time decay rates, comparing them to to test whether the Fisher symmetric gap still clocks entropy decay when detailed balance is broken.
- 21_gkls_nonrev_rate_vs_spectrum.py
-
Compares effective decay rates of entropy and Fisher energy in non reversible Markov chains and their GKLS lifts with different spectral gaps. Using the same class of non reversible generators as in scripts 19 and 20, it computes the gap of , the gap of the symmetric part in the metric, and the Fisher Laplacian gap . It then measures decay rates of entropy and Dirichlet energies from time series and compares them to these spectral scales, clarifying which gap is the relevant UIH irreversible clock in the non reversible regime.
- 22_gkls_nondiagonal_coherent_density_checks.py
-
Analyses a single explicit non diagonal GKLS model with genuine coherences for a two level system. The model has a Hamiltonian and dissipators given by amplitude damping and dephasing, leading to a unique stationary state with non zero coherences in the computational basis. The script diagonalises , transforms the generator into this eigenbasis, extracts the induced density generator , verifies reversibility and stationarity, and compares the GKLS density sector Dirichlet form in the BKM metric with the classical Fisher Dirichlet built from . This provides an explicit coherent, non diagonal test of reversible UIH density hydrodynamics.
- 23_gkls_random_qubit_density_ensemble.py
-
Builds an ensemble of genuinely coherent qubit GKLS models with random Hamiltonians and Lindblad operators for amplitude damping, excitation and dephasing. For each accepted model it finds a full rank, nondegenerate stationary state , diagonalises it, and transforms the generator into the eigenbasis. It then extracts the effective two state density generator , checks stationarity and approximate detailed balance, constructs the BKM metric and symmetric part , and compares the GKLS density sector Dirichlet form to the classical Fisher Dirichlet energy. The script summarises residuals and mismatches across the ensemble.
- 24_gkls_random_qutrit_density_ensemble.py
-
Extends the random coherent density sector tests to qutrits. For three level systems with random Hermitian Hamiltonians and Lindblad operators for nearest neighbour jumps and dephasing, the script finds full rank, nondegenerate stationary states, diagonalises , and extracts effective three state density generators . It verifies stationarity and small detailed balance residuals, constructs the BKM metric and symmetric density sector operator , and compares the GKLS density sector Dirichlet form with the classical Fisher Dirichlet built from over many mass conserving perturbations.
- 25_fp_to_markov_to_gkls_realisation.py
-
Provides a constructive realisation of a Fisher metriplectic Fokker-Planck flow as a Markov chain and then as a GKLS semigroup. It starts from an overdamped Fokker-Planck model on a periodic one dimensional domain with stationary density , discretises the domain to build a reversible nearest neighbour generator , and lifts to a diagonal GKLS generator. It then constructs the BKM metric and symmetric part , extracts an effective density generator from the GKLS, and compares and the density sector Dirichlet form to their classical counterparts, demonstrating an explicit FP to Markov to GKLS UIH ladder.
- 26_gkls_coherent_dressing_fp_chain.py
-
Explores coherent dressing of the FP to Markov to GKLS realisation. Starting from the reversible chain and diagonal GKLS generator of script 25, it builds a tight binding Hamiltonian on the same lattice and forms the Hamiltonian superoperator . The total generator leaves the same diagonal stationary state invariant. The script verifies that the effective density generator extracted from matches and that the density sector Dirichlet form in the BKM metric agrees with both the classical Fisher Dirichlet and the dissipative only GKLS Dirichlet, showing explicitly that many coherent GKLS dressings share the same irreversible Fisher hydrodynamics on densities.
- 27_gkls_density_sector_tomography.py
-
Implements density sector tomography for a diagonal jump GKLS generator built from a reversible three state Markov chain. After lifting a reversible generator to a diagonal GKLS model and confirming that is stationary, it extracts the exact density generator and then reconstructs a tomographic estimate using only diagonal states and their instantaneous responses. It compares to and , builds the BKM metric and symmetric part , and checks that the Fisher Dirichlet built from matches both the classical Dirichlet and the GKLS density sector Dirichlet. This shows that UIH density hydrodynamics can be reconstructed from density response data alone.
- 28_gkls_bloch_k_split_checks.py, 29_gkls_bloch_metriplectic_split.py
-
Work directly with real Bloch-space representations of GKLS generators for dissipative qubits, reconstruct the BKM metric, perform the metriplectic split , and verify that the symmetric and skew blocks satisfy the metric adjoint conditions with small residuals. They implement the one current two quadratures picture in full Bloch space and document the associated entropy and Fisher decay diagnostics.
- 32_uih_two_quadratures_k_flow_checks.py
-
Realises the abstract one current two quadratures experiment. It constructs a random symmetric positive definite metric , builds and satisfying the -adjoint conditions, and compares the decay of under the pure gradient flow and the full flow . It checks production identities and fits late-time decay rates, confirming the UIH picture that redistributes modes without changing the dissipative spectrum.
- 33_gkls_bloch_two_quadratures_k_flow_checks.py
-
Two-quadratures -flow test in a concrete driven, damped qubit GKLS model expressed in Bloch coordinates and equipped with the BKM metric at the stationary state. The GKLS generator induces a real matrix on Bloch space and a positive matrix ; restricting to the traceless sector, the script constructs the metric adjoint , forms the split and , and verifies symmetric and skew. It then studies the quadratic functional , computes the dissipative spectrum from , identifies , evolves both and via matrix exponentials, and fits late-time decay rates, illustrating the UIH one-current two-quadratures structure in a physical qubit example.
- 33a_uih_two_quadratures_visual_explorer.py
-
Provides the interactive visual explorer. It implements the flow for a user-controlled mixing parameter , plotting trajectories in the metric space and the decay of the quadratic functional . This serves as an educational tool to visualise how and combine in a single current with two quadratures.
A.3 IBM Quantum hardware suite
The following scripts implement the IBM Quantum experiments.
- 34_uih_k_tomography_ibmq_qubit_test.py
-
Carries out process tomography on idle circuits for a one-qubit IBM backend, reconstructs the quantum channels in the Pauli basis, performs a matrix logarithm to infer an effective real generator , and checks consistency across idle depths. It computes the BKM metric at the stationary state and performs the metriplectic split in the hardware BKM geometry.
- 35_qapi_semigroup_scaling.py
-
Tests the semigroup property and time-homogeneous scaling of the hardware channel family by comparing powers of the reconstructed channel with channels obtained at different idle depths, quantifying deviations from perfect semigroup behaviour and feeding into the error budget.
- 37_qapi_bkm_speed_limit_test.py
-
Implements the BKM speed limit diagnostics. Using the reconstructed and BKM metric , it computes the dissipative spectrum of , extracts the smallest eigenvalue , and compares the late-time decay of the BKM quadratic functional under the hardware flow and the pure flow, confirming that both are controlled by the UIH clock .
- 38_qapi_bkm_curvature_test.py
-
Performs the curvature test. It prepares small unitary perturbations around the hardware stationary state, reconstructs the perturbed density matrices, computes the exact quantum relative entropy and its quadratic BKM prediction , and compares the two across several directions. The script reports ratios and confirms that the BKM metric extracted from the device is the local curvature of quantum relative entropy.
- 40_qapi_two_qubit_k_tomography_ibmq_test.py
-
Performs a full two-qubit K tomography experiment on an IBM superconducting backend, extending the one-qubit diagnostics of the earlier scripts to a 16 dimensional operator space. It constructs idle identity circuits on a chosen two-qubit pair, calls the Qiskit QAPI / Runtime process tomography tools to reconstruct the noisy idle channel as a completely positive trace preserving map in the Hermitian two-qubit Pauli basis, and extracts the corresponding Pauli transfer matrix . From it infers an effective generator via a matrix logarithm on the traceless sector, computes the stationary state , and saves all objects (tomography data, stationary state, generator blocks and basic condition numbers) into
.npzfiles for the downstream metric split and decay-clock diagnostics used in Section 6.4. - 41_qapi_two_qubit_k_uih_metric_split_checks.py
-
Loads the two-qubit tomography data from
40_qapi_two_qubit_k_tomography_ibmq_test.pyand implements the UIH metric adjoint split on the 15 dimensional traceless sector. From the stationary state it builds the BKM metric block , forms the effective generator , and computes the metric adjoint . It then defines and , reporting symmetry and skewness residuals for and , the spectrum of the generalised dissipative operator , and basic conditioning data. This provides the primary algebraic smoking gun that the two-qubit idle dynamics realises a metriplectic K split in the hardware BKM geometry. - 42_qapi_two_qubit_k_bkm_fisher_split_checks.py
-
Refines the two-qubit diagnostics by constructing the Fisher dissipative operator associated with the BKM metric at . Starting from and it forms the Fisher operator in natural BKM coordinates, computes its generalised eigenvalues with respect to , and identifies the Fisher gap and higher modes that act as irreversible clocks. The script checks that the dissipative spectrum is strictly positive, quantifies condition numbers and eigenvalue spreads, and compares the Fisher spectrum for several tomography runs to assess robustness against sampling noise and backend drift. The resulting Fisher gaps and eigenmodes are saved for use by the CPTP repair and decay-clock analysis.
- 43_qapi_two_qubit_k_cptp_repair_and_bkm_checks.py
-
Implements a CPTP repair pipeline for the noisy two-qubit channel and repeats the BKM and Fisher diagnostics on the repaired data. It takes a raw Pauli transfer matrix from tomography, projects it onto the nearest completely positive trace preserving map using a convex / spectral repair routine, and recomputes the effective generator on the traceless sector together with the stationary state . From it builds the BKM metric , re-evaluates the K split and Fisher operator, and reports channel-level consistency checks such as , metric adjoint residuals and dissipative spectra. These diagnostics show that the UIH K split and Fisher-clock structure survive explicit CPTP enforcement on the reconstructed two-qubit channel.
- 44_qapi_two_qubit_k_cptp_semigroup_decay_checks.py
-
Uses the CPTP-repaired generator and BKM metric to perform a full semigroup decay-clock test on the two-qubit idle dynamics. For several tomography runs (different shot budgets) it loads , and the Fisher operator , computes the Fisher dissipative spectrum and gap, and evolves a set of initial density matrices under the semigroup out to long times. At each time it evaluates the BKM relative entropy , fits short- and long-time decay rates from , and compares these fitted slopes to the Fisher gap extracted from . The script reports slope-to-gap ratios across multiple initial states and datasets, together with trace and eigenvalue bounds along the flow, providing the dynamic two-qubit confirmation that the Fisher gap acts as an irreversible decay clock on real IBM hardware.
- 48_ibmq_uih_spectrometer_suite.py
-
UIH spectrometer and RG toolkit for IBM two-qubit tomography data. The script reads BKM-orthonormal split files
*uih_split_bkm.npzproduced by the IBM K-tomography pipeline, extracts the regularised generator on the traceless BKM basis, and defines and . For each file it computes the Fisher gap from the smallest positive eigenvalue of , the hypocoercive decay rate from the smallest positive decay rate of , and the UIH couplings , , together with the raw norms and . Thespectrometercommand produces a table of these invariants across all files, while thergandrg2commands implement one and two slow-mode RG steps respectively: they project onto the Fisher-slow sector of , rescale the coarse generator to match the microscopic Fisher gap, and report how , and flow under coarse-graining. This turns the IBM channel tomography data into a practical UIH “universality spectrometer” for hardware noise. - 56_uih_tomography_spectrometer.py
-
Applies the UIH spectrometer to real hardware generators extracted from two–qubit quantum process tomography on IBM Fez. The script loads a reconstructed Bloch–space generator from any of the tomography npz files (raw, CPTP projected, or UIH–split), simulates the synthetic evolution , regresses beta maps, and evaluates cross coherence under index smoothing. It reports the stability of the full closure ( with variance ) and the strong suppression of the diagonal baseline for high–shot datasets, together with the expected pathologies of the diagonal channel at lower shot counts. The resulting beta maps and coherence curves reproduce the IBM analysis in Section 13.
- 57_ibmq_kahler_rg_from_tomography.py
-
Holomorphic Kähler RG diagnostics for IBM two qubit tomography. The script reads a *uih_split_bkm.npz file produced by the IBM K tomography pipeline, extracts the regularised generator on the BKM orthonormal Bloch sector, and removes the stationary Fisher mode to work on the reduced traceless space. It builds the Fisher metric , diagonalises it to identify an active Fisher subspace, and constructs an approximate Kähler triple on this subspace via a metric polar decomposition of the Hamiltonian block. Using the same Gaussian index coarse grainer as in the synthetic Kähler RG suite, the script then pushes and across scales, reconstructs coarse complex structures , and computes the diagnostics , and . The output npz file records these quantities as functions of together with the microscopic Kähler defect, demonstrating that the IBM idle generator supports an emergent Kähler structure on its Fisher active modes and that the Gaussian information RG acts approximately holomorphically at small scales.
A.4 Hypocoercivity, Fisher decay floors and flow
- 49_uih_hypocoercivity_coupling_scan.py
-
Fully synthetic finite-dimensional UIH hypocoercivity coupling scan. The script generates random UIH models with a single stationary mode in dimension , takes , constructs as a symmetric negative semidefinite matrix with a one-dimensional kernel, and as a skew matrix with . It then restricts to the mean-zero subspace and, for each model, computes the Fisher gap of , the hypocoercive rate from the spectral abscissa of , the metric-induced norms and , and the corresponding dimensionless couplings , . Sweeping over dimensions and reversible-sector strength, the script records the distribution of as a function of and writes all diagnostics to a compressed
.npzfile. The results demonstrate on a large random ensemble that consistently provides a coercive decay floor and that remains an order-one quantity controlled by the UIH couplings. - 50_uih_rg_coupling_flow_suite.py
-
UIH renormalisation group coupling-flow suite on random finite-dimensional models. Starting from a synthetic UIH model of dimension with a single stationary direction, the script restricts to the mean-zero sector , computes , , and as in 49_uih_hypocoercivity_coupling_scan.py, and then applies a Fisher-preserving mode-space RG map: it diagonalises on , selects a fixed number of Fisher-slow modes, projects onto their span, and iterates this coarse-graining. An ensemble of models is evolved in parallel across many worker processes, and the script reports ensemble-averaged flows of , , and as functions of the RG step. The output shows that strongly hypocoercive microscopic models with large flow under RG into a four-dimensional Fisher-slow sector with order-one couplings and a stable band of , providing a numerical realisation of the UIH renormalisation group picture.
A.5 Further
- 51_brownian_trap_uhi_clock_checks.py
-
UIH irreversibility clock test for an optically trapped Brownian particle. The script reads the vanMameren-raw.txt time series, which records the position of a colloidal particle in a harmonic optical trap at a sampling rate of . It centres the data, fits an Ornstein-Uhlenbeck drift and diffusion coefficient from the increments, and constructs a coarse-grained Markov generator on a position grid.
The spectral gap of this generator is compared to an empirical entropy decay rate extracted from a non-equilibrium ensemble built by conditioning on large displacements and following the relaxation of this ensemble back to equilibrium. For each time lag the script estimates a coarse-grained density , computes the relative entropy with respect to the stationary histogram, and fits a late-time exponential decay. The observed decay rate is found to be of the same order as and , where is the fitted OU drift and is the first non-zero eigenvalue of the discrete generator, providing an experimental classical confirmation of the UIH irreversible clock.
- 52_pdh_entropy_clock_demo.py
-
Markov gap entropy clock test for an Ozawa style hydrogen tunneling model. The script builds a reversible two state Markov generator for the hydrogen sites in palladium, using the three term temperature dependent fit for the hopping rate reported by Ozawa et al. (Sci. Adv. 10, eady8495, 2024) and a Boltzmann factor with level splitting for the stationary weights and .
For a representative grid of temperatures spanning the phonon dominated and electron mediated regimes, the script generates synthetic trajectories , computes the two state relative entropy and fits a late time exponential decay rate from . The resulting ratios are reported and plotted as a function of , and are found to be unity to numerical precision. This provides a referee ready demonstration, in an Ozawa calibrated two state model, of the UIH Markov gap entropy clock relation , with the low temperature scaling inherited from the published exponent .
- 53_uih_rg_fluctuation_suite.py
-
Implements the Fisher–Jarzynski fluctuation tests and the RG invariance checks for a one dimensional UIH Fisher diffusion. The script constructs stationary densities , , evaluates the microscopic free energy difference , and computes its coarse–grained counterparts for a family of Gaussian RG kernels . It verifies the predicted quartic suppression via a log–log slope close to , and independently runs a UIH Jarzynski simulation over thousands of diffusion trajectories to confirm that reproduces the same to within sampling error. This file provides the reproducible numerical backbone for Section 11.
- 54_uih_kahler_rg_suite.py
-
Builds an explicit Kähler triple on a dimensional lattice K flow and tests the holomorphicity of the UIH coarse graining map. The script constructs the complex structure , a site dependent Kähler metric , and the symplectic form , then pushes these objects through the Gaussian index–space coarse graining to obtain . It evaluates the compatibility defects , , and the holomorphicity commutator , all of which remain at machine precision in this exactly holomorphic toy model. This script supports the Kähler–RG analysis in Section 12.
- 55_uih_cross_coherence_spectrometer.py
-
Implements the UIH cross–coherence spectrometer on a synthetic Fokker–Planck trajectory with snapshots on a point grid. Treating the index of as a spatial coordinate, it regresses onto the full UIH closure and the diagonal baseline , producing beta maps and . Under Gaussian index coarse graining it computes the cross coherence and , demonstrating that the full closure survives smoothing with high coherence while the diagonal approximation quickly decorrelates. This serves as the reference implementation for Section 14.
Each script is written to run in a single command, printing summary diagnostics to standard output and, where appropriate, saving small data files or plots in the local directory. Together, the archive provides a complete numerical and experimental realisation of the universal information hydrodynamics picture described in the main text.
Appendix B Lindblad superoperator and Markov generator in matrix form
For completeness we collect the explicit matrix representations of the GKLS superoperator and the classical generator used in the numerical scripts. This is convenient for readers who want to inspect the discrete operators directly or reproduce the computations in other environments.
B.1 Vectorisation and basis
Let with orthonormal basis . Every density matrix on can be represented as a complex matrix with entries . We choose the column stacking vectorisation
where the component corresponding to the pair is
with a fixed ordering convention such as .
For any matrices of compatible size we recall the standard identity
where denotes the Kronecker product and is complex conjugation.
B.2 Lindblad superoperator
The GKLS generator with jump operators and no Hamiltonian term is
Vectorising, we obtain a linear operator on such that
Using the Kronecker identity we have
and
Therefore
| (B.1) |
In the diagonal jump case of Section 4 we have , so is diagonal. The matrix then has a particularly sparse structure, which is exploited in the script 01_gkls_diagonal_to_markov_checks.py.
B.3 Classical generator
On the classical side the Markov generator is an matrix with elements
and the population vector evolves according to
The detailed balance condition can be written as
where is the diagonal matrix with entries and the star denotes adjoint. This shows that is selfadjoint on with respect to the inner product . In this representation the Dirichlet form controlling the entropy decay is
with .
The scripts treat and as concrete sparse matrices and use standard ODE solvers to integrate the associated linear systems.
Appendix C Numerical parameters and tolerances
The numerical scripts are designed to be simple to read and modify. This appendix records the parameter choices and tolerances used in the default runs, so that published figures can be reproduced to within floating point variation. The scripts expose these parameters near the top of each file as global constants or easily editable variables.
Subsections C.5 to E.19 treat finite dimensional flows, diagonal GKLS models, nonreversible qutrit ensembles and constructive Fokker-Planck to Markov to GKLS realisations.
Subsection C.9 adds a simple soft-matter benchmark, applying the same entropy-clock diagnostics to an independent optical-trap Brownian trajectory.
Subsections D.1 to D.4 document IBM Quantum hardware experiments, while Subsection E presents a Fisher-Lindblad unification suite that links Fisher curvature, Markov gaps and entropy decay rates across continuum, Markov and GKLS sectors.
C.1 Time stepping and grids
For the finite dimensional GKLS and Markov chain checks the natural time scale is set by the spectral gap of the generator. In the diagonal jump example the nonzero rates are of order one, so the relaxation time is also of order one. The script 01_gkls_diagonal_to_markov_checks.py therefore uses a final time of order a few relaxation times and a time grid of several hundred points. A standard adaptive ODE solver such as solve_ivp with relative and absolute tolerances of order suffices to resolve the evolution.
For the lattice hydrodynamic limit in 02_markov_to_fp_limit_checks.py, the spatial grid spacing is chosen as a sequence . The macroscopic time horizon is fixed, for instance , and the ODE solver is run with step size control that ensures a temporal resolution sufficient for all grid spacings. The PDE reference solution is computed on the finest spatial grid, and either interpolated to coarser grids or recomputed on each grid for consistency.
For the continuum Fisher metriplectic checks in 03_fp_fisher_metriplectic_checks.py, the Fokker Planck equation is discretised in space using second order finite differences and integrated in time using an unconditionally stable scheme such as backward Euler or Crank Nicolson. The spatial grid typically has a few hundred points, and the time step is chosen so that the discrete solution remains stable and resolves the decay of the free energy accurately. The script reports the discrepancy between the numerical time derivative and the Fisher right hand side as a function of time.
For the coherent qubit GKLS example in 04_gkls_coherence_elimination_checks.py, the Bloch equations are integrated over a time interval that covers several dephasing times and several periods , with chosen large to realise the overdamped regime. An adaptive ODE solver with tight tolerances suffices.
C.2 Tolerance choices
Each script declares tolerances against which numerical identities are tested. The choices can be adjusted, but a representative set is as follows.
In the GKLS to Markov reduction script:
-
•
population consistency tolerance , defined as the maximum over time of the infinity norm of the difference between the GKLS and Markov population vectors;
-
•
entropy consistency tolerance , defined as the maximum over time of the absolute difference between the quantum relative entropy and the classical Kullback Leibler divergence;
-
•
mass conservation tolerance , defined as the maximum deviation of the total population from unity.
In the Markov to Fokker Planck limit script:
-
•
error norms between the discrete density and the PDE solution are computed in an or sense, and a convergence rate is estimated from a log log fit;
-
•
a convergence rate threshold is used as a loose lower bound for first order behaviour;
-
•
the script prints the values of and the estimated rate, and declares a passing run when the rate exceeds and the finest grid error falls below a specified small value.
In the Fisher metriplectic Fokker Planck script:
-
•
the discrepancy is monitored and required to fall below a tolerance such as at all times;
-
•
the discrepancy between the numerical free energy derivative and the Fisher quadratic form, , is also monitored and required to stay below a similar tolerance.
In the coherent GKLS elimination script:
-
•
the fitted relaxation rate for the population decay is compared to the theoretical value , and the relative error is required to be small, e.g. less than a few percent in the regimes where is large;
-
•
the maximum deviation between the populations under the full GKLS evolution and the effective two state chain over a chosen time window is reported, and runs where this deviation is small compared to the absolute change in population are identified as clean overdamped examples.
These values are indicative and can be tuned. The important point is that each script makes its own error thresholds explicit and reports discrepancies quantitative enough for a referee to assess the strength of the numerical support.
C.3 K tomography and operational equivalence
The script 05_K_tomography_on_lattice.py plays a slightly different role from the other numerical checks. Rather than starting from a known metric operator and reversible operator and verifying their consequences, it inverts the process and shows that and can be reconstructed, up to operational equivalence, from a finite set of Fisher style probes.
The set up is deliberately simple. On a one dimensional periodic lattice of size we construct a symmetric positive definite matrix as an identity plus a small multiple of the discrete Laplacian, and a skew symmetric matrix as the centred finite difference convection operator. For a family of probe potentials in a low dimensional real Fourier subspace we form the synthetic response data
Stacking the data in matrix form gives a linear system with , where contains the probes and the responses. Since the number of probes is small compared to , the system is highly underdetermined: only the action of on the probe subspace is constrained and all components orthogonal to that subspace are left in a large nullspace.
The script first computes the minimal norm solution using a pseudoinverse of , then projects the reconstructed blocks onto the metriplectic class by setting
where and denote the left and right blocks of . By construction is symmetric, is skew symmetric, and the pair defines a candidate complex mobility .
The diagnostics in the script make two points explicit. First, the structural constraints are satisfied to machine precision: and lie at the level of floating point roundoff. Second, the forward map generated by reproduces all probe responses to high accuracy, both on the training probes used to build and on an independent test set of probes drawn from the same Fourier subspace:
stays at the level of numerical noise when measured across all . In contrast, entrywise differences and need not be small, which is exactly what one expects in an underdetermined inverse problem where the data constrain only the probe subspace.
This behaviour makes a conceptual point that is central to the universal information hydrodynamics viewpoint. The complex mobility is only defined operationally up to its action on the directions that are actually probed. Many different microscopic pairs can produce exactly the same responses in a given family of Fisher probes, and are therefore indistinguishable at the level of density dynamics. The K tomography script illustrates this explicitly: it recovers a canonical metriplectic representative that lies in the same operational equivalence class as and matches all observable data, but it does not and need not coincide with the original matrices off the probe subspace. From the point of view of the density manifold, it is exactly this equivalence class that matters, not any particular microscopic representative.
C.4 Hyperbolic Fisher Cattaneo speed test
To illustrate that Fisher regularisation admits a bona fide relativistic propagation scale, we implemented a one dimensional Fisher Cattaneo test in the script 30_fisher_cattaneo_relativistic_speed_checks.py. The model is the hyperbolic regularisation of a Fisher diffusion with mass , Planck constant and an effective signal speed ,
| (C.1) |
for which the characteristic front speed is
| (C.2) |
We work on the periodic domain with , using grid points and grid spacing . The parameters are fixed at
| (C.3) |
so that , and . Time integration is performed with a second order explicit scheme with
| (C.4) |
The initial condition is a sharply localised bump centred at the origin, normalised so that . During the evolution we track the position of the right moving front using a fixed relative amplitude threshold, namely the smallest for which
| (C.5) |
A linear fit of the front position on the interval yields a measured speed
| (C.6) |
The relative error stays at the few percent level despite the rather modest choice of grid and time step. This hyperbolic Fisher Cattaneo test therefore confirms that the Fisher regularised dynamics supports a finite propagation speed consistent with the effective light speed , providing a simple numerical realisation of relativistic signal bounds within the UIH framework.
C.5 Asymptotic decay clocks in reversible Markov chains
The script 31_uih_asymptotic_decay_clock_qutrit_fp_checks.py probes the UIH prediction that for reversible Markov chains the spectral gap of the symmetrised generator fixes an asymptotic information decay clock that is robust under discretisation and microscopic realisation.
Given a finite state reversible generator with stationary distribution , we consider the classical relative entropy
| (C.7) |
and define the instantaneous decay rate
| (C.8) |
Writing and , reversibility implies that is symmetric with non positive spectrum and a simple zero eigenvalue. The Markov spectral gap is
| (C.9) |
UIH predicts that for generic initial data the late time behaviour is governed by
| (C.10) |
so that acts as a universal asymptotic clock for information decay in the reversible sector.
The script constructs two examples.
Qutrit reversible chain.
A three state reversible generator is built from random symmetric conductances and a random stationary distribution with full support. The gap is computed from the symmetrised operator as above, giving
| (C.11) |
We sample many random initial distributions , evolve them using diagonalisation of , and evaluate on a uniform time grid up to . Finite difference estimates of are formed on the interior of the grid and fitted in a late time window. Across a broad set of initial data the fitted rates cluster around
| (C.12) |
with transient overshoots at early times when multiple modes contribute significantly.
High resolution FP like chain.
We then construct a nearest neighbour reversible generator on a periodic lattice of size , with a smooth multi well potential and Gibbs stationary distribution at inverse temperature . Symmetric conductances to neighbours enforce detailed balance. The raw gap is computed, and the generator is rescaled so that
| (C.13) |
Repeating the same entropy based diagnostics, we again find that for a wide range of initial conditions the late time decay rates approach
| (C.14) |
with early time rates that can be larger when higher modes are excited.
The key observation is that the Markov spectral gap sets a common asymptotic information decay scale across two very different realisations: a low dimensional qutrit density sector of a GKLS model, and a high resolution Fokker Planck limit. This supports the UIH identification of as a universal irreversible clock in the reversible Markov sector.
C.6 Finite dimensional K flow with one current and two quadratures
The script 32_uih_two_quadratures_k_flow_checks.py tests the UIH picture of a single current with two quadratures in a finite dimensional setting, using an abstract real vector , a fixed symmetric positive definite metric , and an -adjoint mobility operator
| (C.15) |
with symmetric and skew in the geometry.
In this subsection we adopt the convention that positive eigenvalues of correspond to decay and evolve with ; in the IBM Bloch and GKLS sections below we instead treat as the real GKLS generator with negative dissipative spectrum and identify the metriplectic mobility with . The two sign conventions are equivalent up to this overall minus sign.
We construct:
-
•
A random symmetric positive definite metric with condition number of order .
-
•
A symmetric matrix and a skew symmetric matrix in the Euclidean sense, and then enforce the metric adjoint conditions numerically:
(C.16)
The key quadratic functional is
| (C.17) |
and we consider two linear flows on :
| (C.18) |
The script computes
| (C.19) |
along both flows using high resolution matrix exponentials, and compares the numerical derivative of with the analytic expression at each time. We then form log linear fits of in a late time window to extract decay rates and for the pure flow and full flow respectively.
The diagnostics reported are:
-
•
Dimension .
-
•
Condition number of of order .
-
•
Norm of the antisymmetric part of of order , corresponding to a relative symmetry residual .
-
•
Norm of the symmetric part of of order , with relative skewness residual .
-
•
The generalised eigenvalues of are all positive, with
(C.20)
For five random initial conditions satisfying , we find:
-
•
The production residual is at most of order , with relative errors around to .
-
•
The fitted decay rate for the pure flow matches the spectral prediction
(C.21) to numerical precision.
-
•
The full flow has a larger decay rate
(C.22) essentially independent of the chosen initial condition.
These tests confirm in a simple finite dimensional setting that:
-
1.
The metric adjoint split yields a genuine no work direction , in the sense that does not appear in the production of the quadratic information functional .
-
2.
The smallest positive eigenvalue of the dissipative operator in the geometry sets the natural decay clock for the pure gradient flow.
-
3.
Adding the reversible component can accelerate decay by mixing eigenmodes, without changing the spectrum of the dissipative channel.
C.7 GKLS qubit K split in Bloch BKM geometry
The script 33_gkls_bloch_two_quadratures_k_flow_checks.py realises the same K split picture in a concrete GKLS model for a driven, dissipative qubit, expressed in Bloch coordinates and equipped with the BKM metric.
We start from a four dimensional Bloch vector with enforcing unit trace. The full linearised GKLS generator in this basis has the form
| (C.23) |
where the block on the traceless coordinates encodes both Hamiltonian precession and dissipative damping. The stationary Bloch vector is obtained by solving , giving a mixed thermal state with nontrivial coherences.
At this stationary state we compute the four dimensional BKM metric and its traceless block using the exact matrix formula for qubits in Bloch form. The diagnostics reported are:
-
•
A well conditioned metric with .
-
•
Metric adjoint residuals and of order to .
-
•
A strictly positive dissipative spectrum for with
(C.24)
We then consider the quadratic functional
| (C.25) |
for traceless Bloch displacements , and compare the decay for two flows
| (C.26) |
where the latter isolates the gradient component.
Using exact matrix exponentials on a time grid up to with 400 samples, and five random initial conditions normalised to , we find:
-
•
The production identity holds to relative errors below .
-
•
The decay rate for the pure flow satisfies
(C.27) within a few tenths of a percent across all initial conditions.
-
•
The full flow decays faster, with fitted rates
(C.28)
This GKLS example confirms that the metriplectic K split is not restricted to abstract finite dimensional models. It arises naturally in a realistic qubit Lindbladian once the BKM metric at the stationary state is used to define the adjoint. The dissipative spectrum again controls the irreversible decay clock, while the Hamiltonian component accelerates decay only by rotating within the metric geometry.
C.8 Visual two quadratures explorer
The script 33_uih_two_quadratures_visual_explorer.py provides an interactive visual companion to the finite dimensional K flow tests. It implements a two dimensional parametrised family of flows
| (C.29) |
for a fixed metric , fixed symmetric part and fixed skew part , with a user controlled mixing parameter . The overall minus sign matches the convention in which the real GKLS generator is in the abstract metriplectic notation, so that positive eigenvalues of correspond to decay.
The tool plots on the same axes:
-
•
The trajectories in the metric space.
-
•
The corresponding quadratic functional .
-
•
The decomposition of the current into gradient and Hamiltonian quadratures, as measured by and .
By varying the initial condition and the mixing angle , one can visually confirm the algebraic statements of the previous subsections: the gradient component fixes entropy production and determines the asymptotic decay scale, while the Hamiltonian component alters the paths and can change how quickly the asymptotic regime is reached, but does not directly contribute to the production of . This explorer is not used as a quantitative diagnostic, but as a geometric illustration of the one current two quadratures picture.
C.9 Brownian trap entropy clock benchmark
As a classical benchmark for the UIH irreversibility picture we analyse an optical trap experiment for a single colloidal particle. The dataset vanMameren-raw.txt consists of measurements of the particle position along one axis, sampled at a fixed rate , so that successive samples are separated by . Only positions along one axis, in nanometres, are recorded. The data, provided by van Mameren and Schmidt, correspond to a particle in a harmonic optical potential and are well modelled by an Ornstein-Uhlenbeck (OU) process
| (C.30) |
with Gaussian stationary state .
We use the raw time series as supplied, without trimming individual segments or discarding samples, and apply only centring and coarse-graining onto a fixed position grid. We first fit an OU drift and diffusion coefficient from the increments. Writing and regressing gives an estimate for the drift rate and
| (C.31) |
for the diffusion coefficient. For the van Mameren data we find per time step, corresponding to a physical relaxation rate and a correlation time of order .
To place this system in the UIH framework we construct a coarse-grained Markov generator for the position coordinate. We discretise the range of observed positions into intervals with edges , assign each sample to a bin, and estimate the one-step transition matrix
| (C.32) |
from the observed counts. The corresponding discrete generator is . The long-time occupancy histogram defines a stationary distribution that we interpret as a discrete approximation to . Diagonalising we find a leading eigenvalue and a first non-zero eigenvalue , whose negative sets the spectral gap of the coarse-grained Fokker-Planck operator. For the present data the gap is per time step, corresponding to a physical rate .
A single long OU trajectory is already stationary, so global relaxation is not directly visible in . To probe the irreversible clock we therefore construct a non-equilibrium ensemble from the same time series. We select all indices for which , where is the empirical standard deviation of the stationary distribution and is a fixed threshold (we use ), and for each such index follow the subsequent positions for a fixed maximal lag . At each lag this yields an ensemble of positions , which defines an empirical density by coarse-graining these values into the same position bins. By construction is biased towards the tails of the stationary distribution.
Using the stationary histogram as a reference, we compute the relative entropy
| (C.33) |
with a small regularisation to avoid empty bins. As increases the ensemble relaxes back towards equilibrium and decays. Plotting against reveals a clear linear regime at intermediate times, from which we extract an empirical decay rate by least-squares fitting. For the van Mameren data we obtain per time step, corresponding to .
In the exact continuum OU model with Gaussian initial data one has at late times, so that the entropic clock is fixed by the spectral gap . Here we work with a thresholded, non-Gaussian ensemble and a coarse-grained generator , so we do not expect exact agreement. Nevertheless, for the present benchmark we find
| (C.34) |
so that lies between and up to factors of order unity. We view this as a simple sanity check: a soft-matter Brownian system, analysed with no tuning beyond fixed thresholds and coarse-graining, exhibits an irreversible information clock controlled by the same Fisher-Dirichlet gap that organises the Markov and GKLS examples elsewhere in the paper. This classical example does not add new physics, but it shows that the UIH irreversibility picture extends in a straightforward way beyond the quantum test cases.
Appendix D IBM Quantum Computer Experiments
D.1 IBM Quantum K tomography diagnostic
The script 34_uih_k_tomography_ibmq_qubit_test.py implements a full K tomography experiment on an IBM superconducting qubit, demonstrating that the noisy idle dynamics on hardware realises a metriplectic K split in the BKM metric at the device stationary state.
We use the Qiskit [17] Runtime API to select an available one qubit backend, in our runs ibm_fez, and construct a simple idle circuit consisting of a depth four sequence of identity gates on a single qubit. This idle circuit is passed to the ProcessTomography routine in qiskit-experiments, which reconstructs the corresponding quantum channel as a completely positive trace preserving map in Pauli transfer matrix form , expressed in the Hermitian Pauli basis
.
From the reconstructed we compute:
-
•
The stationary state as the fixed point of , either by iterating the channel or by extracting the eigenvector of with eigenvalue one and reshaping.
-
•
The BKM metric at , using the closed form qubit expression in the Pauli basis. Its traceless block is well conditioned, with condition number of order unity.
-
•
The traceless sector superoperator obtained by restricting to the Pauli components.
Assuming the idle channel approximates a short time propagator for some effective generator , we compute
| (D.1) |
using the matrix logarithm and discarding small imaginary parts, which have norm below in the reported run. We then define the metric adjoint
| (D.2) |
and the symmetric and antisymmetric parts
| (D.3) |
The hardware diagnostics are:
-
•
The symmetry residual is of order .
-
•
The skewness residual is also of order .
-
•
The eigenvalues of are strictly positive, with
(D.4)
This constitutes a direct experimental realisation of the UIH K split on quantum hardware: the effective generator extracted from IBM process tomography decomposes into a dissipative Fisher BKM gradient part and a Hamiltonian part that is metric skew and hence does no work in the BKM geometry. All ingredients are obtained from tomography and basic linear algebra, with no model assumptions about the microscopic origin of the noise.
D.2 Semigroup scaling across idle depths
The script 35_qapi_semigroup_scaling.py complements the K tomography experiment by testing the semigroup property of the reconstructed idle channels at two different depths, thus probing the time homogeneity of the effective K on hardware.
We consider two idle circuits,
| (D.5) |
on the same IBM backend, again using ibm_fez in our runs. For each depth we run a one qubit process tomography experiment using ProcessTomography with 2048 shots per setting, yielding two superoperators
| (D.6) |
for the depth two and depth eight idling channels respectively. Restricting to the traceless Pauli block we obtain matrices and .
We interpret these as short time propagators
| (D.7) |
for some effective generator . Using the principal matrix logarithm we define
| (D.8) |
after verifying that the imaginary parts of are small. The script reports:
-
•
Imaginary parts of and with norms of order or smaller.
-
•
Frobenius norms , , and mismatch corresponding to a relative deviation of order .
At the channel level, however, the semigroup prediction can be tested more directly by comparing
| (D.9) |
in Frobenius and operator norms. The reported run finds
| (D.10) |
with a similar small mismatch in operator norm.
These results support the picture that a single effective generator , of the type reconstructed in the K tomography test of Subsection D.1, controls the dissipative idle dynamics across a range of time scales, with deviations that are small at the channel level and compatible with finite sampling noise and drift. This is the semigroup counterpart of the geometric K split, showing that the same governs both the metric decomposition and the time scaling of the irreversible flow.
While the inferred generators and differ at the tens of percent level, the induced channels agree to within a few percent in Frobenius norm; this is the operationally relevant diagnostic, and the discrepancy at the level is compatible with sampling noise and the logarithm’s sensitivity to small spectral perturbations.
D.3 BKM speed limit test on IBM Quantum
The script 37_qapi_bkm_speed_limit_test.py combines IBM process tomography with the BKM metric to probe an information theoretic speed limit on hardware expressed in terms of the dissipative spectrum in the BKM geometry.
We reuse the idle depth four channel reconstructed as in subsection D.1, and its traceless block . The stationary state is computed from , and the BKM metric and traceless block are built at . In the reported run for the speed limit test has Bloch vector length around , so it is neither pure nor maximally mixed, and is well conditioned with condition number .
The effective traceless generator is defined via
| (D.11) |
whose imaginary part is negligible. The metric adjoint split
| (D.12) |
again yields metric symmetry and skewness residuals at the level of .
The dissipative spectrum is extracted from ,
| (D.13) |
which has strictly positive eigenvalues
| (D.14) |
The UIH prediction is that the smallest sets a natural decay scale for the quadratic functional
| (D.15) |
under the pure gradient flow , with an effective clock of order .
To test this, we compare two flows on traceless Bloch displacements :
| (D.16) |
for five random initial conditions normalised to . The evolution is computed by exact matrix exponentials on a time grid up to
| (D.17) |
with 400 samples. For each trajectory we check:
-
•
The production identity , whose numerical residuals stay at the level of in absolute units and in relative terms.
-
•
Late time decay rates and obtained by linear fits of over the final third of the time window.
The fitted rates are of order
| (D.18) |
The precise prefactor depends on the definition of the time unit associated with the idle depth and on how close the chosen time window is to the strictly asymptotic regime.
The important point is that both rates scale with the UIH clock extracted from the BKM curvature, and that the full flow decays faster than the pure flow, as expected when the reversible channel mixes eigenmodes without changing the dissipative spectrum. Together with the K split and semigroup tests of Sections D.1 and D.2 this provides a dynamic demonstration for the UIH speed limit mechanism on hardware.
D.4 BKM curvature test on IBM Quantum
Finally, the script 38_qapi_bkm_curvature_test.py performs a direct experimental test of the statement that the BKM metric at the stationary state is the local curvature of quantum relative entropy, by comparing the exact relative entropy to its quadratic approximation for small unitary perturbations around the hardware stationary state.
We again use a one qubit IBM backend (ibm_fez) and prepare a stationary state via an idle depth four circuit and tomography. In the reported curvature run is close to pure, with Bloch vector of length very close to one, pointing near the north pole, and eigenvalues
| (D.19) |
The BKM metric and its traceless block are computed from using the exact qubit formula. The eigenvalues of are approximately
| (D.20) |
with condition number . This reflects the strong curvature anisotropy near a nearly pure state.
We then generate three small perturbations of by applying single qubit rotations of a fixed small angle around each Pauli axis:
| (D.21) |
For each of the four states we perform simple Pauli tomography by measuring the three expectation values of with 8192 shots per axis. This yields Bloch vectors , and thus displacements
| (D.22) |
For each displacement we evaluate:
-
1.
The true quantum relative entropy using the eigenvalues and logarithms of the reconstructed density matrices.
-
2.
The quadratic BKM prediction
(D.23)
In the reported run the results are:
-
•
For the X rotation, and , with ratio .
-
•
For the Y rotation, , , ratio .
-
•
For the Z rotation, , , ratio .
The mean ratio over the three directions is
| (D.24) |
with standard deviation . Given the relatively large rotation angle and the strong curvature anisotropy near a nearly pure state, this level of agreement is consistent with the expected truncation error in the second order expansion of the relative entropy.
Crucially, the dependence of on the displacement direction and magnitude is entirely fixed by the BKM metric at , and the hardware data follow this direction dependence. This experiment therefore provides a strong geometric demonstration: the BKM metric extracted from the IBM device via its stationary state curvature yields a quantitatively accurate local approximation of quantum relative entropy, confirming the UIH identification of the BKM metric as the local information curvature of the device.
Appendix E Fisher-Lindblad numerical unification
This appendix assembles a numerical Fisher-Lindblad unification suite based on finite GKLS models, reversible and nonreversible Markov chains, and constructive Fokker-Planck discretisations. The scripts 06_gkls_fp_G_unification_checks.py through 29_gkls_bloch_metriplectic_split.py test, in increasing generality, the following claims:
-
•
diagonal and coherent GKLS generators induce a canonical Fisher Dirichlet operator on the density sector that coincides with the classical Dirichlet form, even in nonreversible chains;
-
•
the cost-entropy inequality of the metriplectic theory is saturated mode by mode in the natural Fisher coordinates;
-
•
entropy and Fisher energy decays are governed by the Markov spectral gap, with Fisher curvature gaps providing a coercive floor but not the dominant rate;
-
•
continuum Fokker-Planck Fisher flows can be realised as the density sector of explicit GKLS semigroups, and the one current two quadratures split extends to full Bloch space generators in coherent qubit models.
The remainder of this appendix documents these tests script by script.
E.1 Numerical suite and scope
Throughout this subsection we use the same sign convention as in the Bloch K split: the real generator is defined by , so its dissipative spectrum has negative real parts and the associated metriplectic mobility is .
The IBM K tomography and BKM curvature experiments of Sections D.1-D.4 showed that a noisy idle channel on hardware realises a metriplectic split in the BKM metric at the device stationary state, and that the smallest dissipative eigenvalue of sets a natural information theoretic decay clock. These are genuinely experimental statements: both the metric and the generator are reconstructed from process tomography and simple linear algebra, with no microscopic model for the noise.
In this subsection we move to controlled numerical GKLS and Fokker-Planck models to pin down the canonical irreversible slice and to show how it arises from Lindblad dynamics, reversible Markov chains and their continuum limits.
E.2 Canonical irreversible slice and the Fokker-Planck limit
We start from a one dimensional overdamped Langevin model and its Fokker-Planck generator on a periodic domain, discretised on a fine grid. From the discretisation we extract a reversible Markov generator with stationary density and define the canonical Fisher operator
| (E.1) |
By construction the irreversible drift for a test potential in the density representation can be written in two equivalent ways,
| (E.2) |
where denotes pointwise multiplication. Script 06_gkls_fp_G_unification_checks.py verifies three facts at high resolution.
First, is symmetric to numerical precision in the weighted inner product and its skew part is purely trace like, so is a genuine Fisher Dirichlet operator. Second, the action of on a catalogue of Fisher probes reproduces the discretised Fokker-Planck drift to relative accuracy of order , modulo the expected overall diffusion scale.
Third, if one attempts a naive tomographic reconstruction from the same probes without imposing structure, then is symmetric but fails to match and only fits the drift on the probe subspace to order . The conclusion is that the canonical irreversible slice is (E.1). Tomography identifies an equivalence class of symmetric operators that agree on the probed directions, and the Fisher metriplectic theory singles out the Markovian representative as the one that descends from a local Fokker-Planck structure.
E.3 GKLS jump models and classical Fisher geometry
We then move to finite GKLS models. Script 07_gkls_to_markov_G_unification_checks.py takes a four level thermal GKLS jump model with Hamiltonian , Lindblad jump operators implementing upward and downward transitions between energy levels, and a Gibbs stationary state . Restricting the GKLS dissipator to the diagonal in the energy basis produces a classical Markov generator with stationary probabilities .
The script checks that
| (E.3) |
that detailed balance holds, and that the classical Dirichlet operator
| (E.4) |
is symmetric and negative definite on the subspace orthogonal to constants. Two representations of the quadratic form agree to machine precision:
| (E.5) |
for the natural choice of conjugate variables and . Thus the irreversible piece of the GKLS dynamics has a unique classical Fisher realisation on densities and the canonical Fisher operator is again .
E.4 Cost entropy inequality in finite Fisher geometry
Script 08_cost_entropy_inequality_markov_checks.py uses the same four state thermal chain to probe the cost entropy inequality in the finite dimensional Fisher geometry. For the canonical Fisher metric defined by the quadratic forms
| (E.6) | ||||
| (E.7) |
are evaluated in the eigenbasis of the metric operator . For each positive eigenmode and for random linear combinations in the positive eigenspace the ratio
| (E.8) |
is equal to to numerical precision.
Remark (finite dimensional thermodynamic metric).
For later reference it is convenient to rewrite the minimal cost as a quadratic form on the tangent space of the probability simplex. In the canonical Fisher geometry defined by we set
| (E.9) |
which is strictly positive on the subspace orthogonal to constants. In this notation the cost entropy inequality takes the form
| (E.10) |
with equality when is collinear with the Fisher gradient direction . Thus can be read as a thermodynamic length type metric on the finite state Fisher manifold picked out by , in the limited sense that it controls pathwise cost at fixed entropy production. In the present paper this plays only a bookkeeping role; no general thermodynamic geometry on the simplex is developed beyond the identities already established.
This shows that in the natural normal form for the Fisher geometry the inequality derived in the metriplectic framework is in fact an equality mode by mode. The cost entropy inequality is therefore not an external constraint but an encoded property of the Fisher metric associated with .
E.5 Metric K splitting for coherent qubits and qutrits
Scripts 09_gkls_K_splitting_qubit_checks.py and 11_gkls_qutrit_K_and_G_unification_checks.py lift the analysis to genuinely quantum systems with coherences. For a thermal qubit and a thermal qutrit with fixed Hamiltonians and jump operators, we construct the real generator on a linearised state space consisting of populations and real and imaginary parts of coherences. The stationary state induces a Fisher information metric on , and we compute the adjoint .
The metriplectic splitting
| (E.11) |
is then obtained by taking the symmetric and antisymmetric parts in the inner product,
| (E.12) |
The scripts confirm that:
| (E.13) | ||||
| (E.14) |
where is the dissipative block arising from the GKLS dissipator and is the Hamiltonian commutator block. The residuals and are of order . Moreover is metric self adjoint and is metric skew adjoint:
| (E.15) |
to the same numerical tolerance. This provides an explicit GKLS realisation of the abstract metriplectic structure with the metric induced by .
E.6 Classical density sector as quantum Fisher block
Script 10_gkls_density_block_G_unification_checks.py takes the coherent qubit example and projects the symmetric part onto the density sector, giving a block acting on the populations. Independently the population dynamics of the GKLS model define a two state Markov generator with stationary distribution . The canonical classical Fisher operator is again .
The script shows that
| (E.16) |
to numerical precision, and that both the drift representation and the Dirichlet form coincide with their Fisher counterparts . Script 11 extends this to the qutrit, where the population block of the symmetric part extracted from the full 9 dimensional GKLS generator reproduces the three level Markov generator, its Dirichlet form and its cost entropy structure. In other words, the classical information hydrodynamics on densities is literally the density block of the quantum metriplectic generator obtained from the GKLS .
E.7 Entropy decay, Fisher curvature and spectral gaps
The next group of scripts analyses how the various gaps control entropy decay. Given a reversible Markov generator with stationary distribution and symmetrised generator with , we define the Markov spectral gap as the smallest positive value of over the spectrum of . The Fisher Dirichlet operator has associated curvature operator and Fisher curvature gap given by the smallest positive eigenvalue of .
Scripts 12_qutrit_markov_entropy_decay_vs_G_gap_checks.py and 13_qutrit_full_quantum_entropy_decay_checks.py investigate the entropy decay of the three state Markov chain induced by the qutrit GKLS model and of the full nine dimensional GKLS generator.
In the Markov case the relative entropy and Fisher quadratic both decay at a rate very close to , while is significantly smaller. In the full GKLS case the decay rates extracted from the Hilbert Schmidt distance and the quantum relative entropy cluster within a few percent of , even though the GKLS generator has additional coherent modes.
This supports a picture in which the density sector Fisher geometry controls entropy production and the dominant decay rate is set by the Markov spectral gap, not by the smallest positive Fisher curvature eigenvalue. The Fisher gap measures local reversible curvature while the Markov gap measures the global relaxation timescale.
E.8 Universality of the density sector across GKLS families
Script 14_qutrit_GKLS_family_universal_density_sector_checks.py considers three different GKLS families on the same three level system. All models share the same energy levels and jump rates but differ in their dephasing structure: uniform projector dephasing, non uniform projector dephasing and non local diagonal dephasing. For each family the script extracts the full real generator , the symmetric part in the induced Fisher metric and the population block of , then compares these with the classical three state Markov generator and its Fisher operator.
Across all families the following quantities are identical up to numerical tolerance: the classical generator , the canonical Fisher operator , the Fisher curvature gap and the Markov spectral gap . In contrast the GKLS spectral gap and the entropy decay rates of the full quantum dynamics depend strongly on the dephasing structure. The density sector geometry is therefore universal for a given collection of jump rates. It is insensitive to how coherences are generated and destroyed, and is determined entirely by the irreversible jump channel.
E.9 Universal Fisher decay clock across discrete and continuum chains
Finally, script 15_qutrit_markov_vs_FP_universal_gap_checks.py connects the density sector of a finite GKLS model to a continuum like Fokker-Planck chain. On one side we have the three state Markov generator inherited from the qutrit GKLS density block with gap . On the other side we construct a high resolution reversible nearest neighbour Markov chain on a periodic lattice representing a discrete Laplacian with uniform stationary distribution. After rescaling so that the Markov gap matches that of the qutrit chain, , we obtain a second chain with very different microscopic structure but the same spectral gap.
For each chain we then evaluate the Fisher Dirichlet quadratic
| (E.17) |
with and , for a catalogue of random initial conditions, and fit an exponential envelope at late times. The empirical rates for the qutrit chain and the FP like chain both cluster extremely tightly around
| (E.18) |
with relative deviations of order for the qutrit chain and of order for the FP like chain. At the same time the Fisher curvature gaps and differ by more than an order of magnitude.
This establishes a universal Fisher decay clock: once the Markov gaps are matched, the Fisher Dirichlet energy decays at essentially the same rate in a three state chain derived from a quantum GKLS model and in a high resolution discrete Fokker-Planck chain. The Fisher curvature gap controls the shape of the local metric, but the global relaxation timescale is set by the spectral gap of the reversible Markov generator.
Taken together, scripts 06 to 15 show that the Fisher metriplectic operator that appears in our axiomatic construction is not an abstract choice. In thermal GKLS models it is canonically realised as the symmetric part of the real generator in the stationary Fisher metric, and its density block is exactly the classical Fisher Dirichlet operator . The cost entropy inequality is saturated mode by mode in the natural Fisher coordinates, the density sector geometry is universal across GKLS families that share the same jump rates, and the actual entropy decay rates are governed by the Markov spectral gap and can be matched between discrete and continuum models. This provides an operational Fisher-Lindblad unification that ties the abstract metriplectic structure directly to GKLS dynamics, reversible Markov chains and Fokker-Planck limits.
Scripts 16 through 29 extend this Fisher-Lindblad picture to diagonal GKLS lifts, coherent dressing, and tomographic reconstruction of both the density sector and the full Bloch generator.
E.10 Diagonal GKLS models and exact Markov reduction on the density sector
Script 16_gkls_diagonal_to_markov_checks.py tests the most basic link between GKLS dynamics and classical Markov chains in the density sector. We work with three level models where the Lindblad operators are diagonal in a preferred basis, so that the populations obey a closed master equation. For each random instance the script constructs both
-
•
the full GKLS generator acting on density matrices, and
-
•
the classical generator acting on populations,
and evolves an ensemble of random initial states under both descriptions.
Two key diagnostics are monitored across the ensemble: the maximal difference between the GKLS induced population drift and , and the mismatch between quantum and classical relative entropies evaluated on diagonal states. Both errors stay at machine precision, with maximal generator and trajectory discrepancies of order and relative entropy mismatches below . The norm of the off diagonal block of the GKLS induced population generator is also at numerical floor, confirming that no hidden coherences leak into the density sector.
This establishes that for diagonal jump GKLS models the restriction of to diagonal density matrices is exactly the classical Markov generator with stationary distribution , and that on this sector the quantum relative entropy reduces identically to the classical Kullback-Leibler divergence. In the UIH framework this means that on the density sector the abstract mobility is fixed as the classical Markov generator , so the irreversible slice is completely determined by the GKLS semigroup.
E.11 Density sector Fisher Dirichlet equality for diagonal GKLS
Script 17_gkls_fisher_dirichlet_checks.py turns from trajectories to the instantaneous Fisher geometry. For the same class of three level detailed balance GKLS models as in script 16, we construct
-
•
the classical Fisher Dirichlet form
(E.19) with , and
-
•
the density sector quantum Dirichlet form
(E.20) where is the symmetric part of the real generator in the BKM metric at .
For each random reversible model an ensemble of mass conserving perturbations is sampled and both quadratic forms are evaluated. Across the ensemble the maximal absolute difference is of order , while the maximal relative error stays at the level. The ranges of the classical and quantum Dirichlet values match to machine precision.
This shows that for diagonal detailed balance GKLS models the density block of the symmetric operator coincides exactly with the classical Fisher Dirichlet operator . The Fisher metric structure on populations that enters our axioms is therefore canonically realised by thermal GKLS dynamics, and in UIH the density block of is literally the classical Fisher Dirichlet operator, not a free modelling choice.
E.12 Coherent diagonal GKLS models and robustness of the density sector
Script 18_gkls_coherent_density_sector_checks.py adds coherent Hamiltonian dynamics and dephasing to the diagonal jump models of scripts 16 and 17. The GKLS generator now contains a nontrivial Hamiltonian part and off diagonal density matrix elements are genuinely excited along the trajectories, but the stationary state remains diagonal in the jump basis.
The script checks three properties across an ensemble of such models:
-
•
the induced population drift from GKLS matches the classical generator to numerical precision,
-
•
the stationary distribution is the same for both descriptions, and
-
•
the quantum and classical Fisher Dirichlet forms on densities still agree.
As in script 17, the maximal discrepancy between and remains at the level, with relative errors below .
This shows that the density sector Fisher Dirichlet equality is robust under the addition of coherent Hamiltonian evolution and dephasing, provided the Lindblad operators are diagonal in the stationary basis. Coherences are present in the full GKLS dynamics, but the density sector geometry and entropy production are still governed exactly by the classical Fisher structure.
E.13 Nonreversible diagonal GKLS models and universal density sector geometry
Script 19_gkls_nonrev_density_sector_checks.py drops detailed balance and considers fully nonreversible three level classical generators with unique stationary distributions . These are lifted to diagonal GKLS generators via the same class of jump operators as before. The aim is to test whether the density sector Fisher Dirichlet equality survives in the absence of reversibility.
For each random nonreversible generator the script constructs and as in scripts 17 and 18, and samples mass conserving perturbations . The maximal stationarity residuals for both and the GKLS lift sit at , confirming that is stationary in both descriptions. The discrepancy remains at or below, with maximal relative errors of order .
The density block of the symmetric operator therefore coincides with the classical Fisher Dirichlet operator even when the underlying Markov generator is nonreversible. Nonreversibility affects the antisymmetric part of the real generator, but the instantaneous Fisher metric on populations is universal and independent of the presence or absence of detailed balance.
E.14 Entropy and Fisher energy decay in nonreversible GKLS chains
Script 20_gkls_nonrev_decay_clock.py turns back to dynamics in the nonreversible setting. For each random nonreversible three state generator and its diagonal GKLS lift, the script evolves an ensemble of trajectories from random initial populations and records four decay diagnostics:
| (E.21) |
where is the classical Kullback-Leibler divergence , is the quantum relative entropy , and are the classical and GKLS Fisher Dirichlet energies of the density perturbation.
For each trajectory an exponential envelope is fitted at late times. Across the ensemble, the mean decay rates of all four diagnostics coincide to high precision,
| (E.22) |
while the Fisher curvature gap of the symmetric density sector operator lies in the much smaller range to . This shows that in nonreversible models the Fisher gap provides only a coercive lower bound. The actual decay of entropy and Fisher energy is governed by the full nonnormal generator and can be several times faster than , while the quantum and classical descriptions agree perfectly on the density sector.
E.15 Nonreversible GKLS: decay rates versus generator spectra
Script 21_gkls_nonrev_rate_vs_spectrum.py refines the analysis of script 20 by comparing decay rates with three spectral objects. For each nonreversible generator the script computes
-
•
the Markov spectral gap ,
-
•
the Fisher Laplacian gap of the symmetric density sector operator built from the Fisher Dirichlet form, and
-
•
the fitted decay rates of .
Over an ensemble of models we find on average and , while the mean decay rates cluster around . The ratios satisfy
| (E.23) |
with very small variance across the ensemble, and all four diagnostics share the same rate within numerical error.
This makes the hypocoercive structure explicit in the Fisher language: the symmetric Fisher operator supplies the universal Dirichlet form and a spectral gap , but the presence of an antisymmetric nonreversible part produces an effective decay scale set by the full Markov generator, comparable to and largely independent of . The quantum GKLS lift inherits this entire picture in its density sector.
E.16 Driven qubit GKLS with coherences: explicit coherent example
Script 22_gkls_nondiagonal_coherent_density_checks.py leaves the diagonal jump class and studies a physically standard driven qubit with genuinely coherent dynamics. The model has Hamiltonian
| (E.24) |
with , , amplitude damping with , and dephasing with . The stationary state has significant off diagonal coherence in the computational basis, with coherence norm of order .
The script diagonalises as and transforms the GKLS generator into this eigenbasis. Restricting to the subspace of diagonal matrices in the basis yields an effective two state generator with stationary distribution . The column sums of vanish to machine precision, the stationarity residual is of order , and the detailed balance condition holds to within .
As in the diagonal case, the BKM metric at and the symmetric part of the real generator are used to form a density sector Dirichlet form. For a catalogue of mass conserving perturbations one finds , with relative errors at the level. Thus even for a fully coherent driven qubit with non diagonal Lindbladians, the eigenbasis density sector is exactly a reversible two state Fisher Markov system, and the density block of is again the classical Fisher Dirichlet operator.
E.17 Random qubit GKLS ensemble and universal Fisher density sectors
Script 23_gkls_random_qubit_density_ensemble.py upgrades the single example of script 22 to an ensemble of random qubit GKLS models. Each model has
-
•
a random Hamiltonian ,
-
•
three Lindblad operators , , , with rates sampled in a moderate range.
For each random model the stationary state is computed and models are accepted only if is full rank, nondegenerate and has nontrivial coherence in the computational basis. In a sample of 20 attempts all 20 models are accepted, with coherence norms ranging from to .
In the eigenbasis of the generator is transformed to and restricted to the diagonal subspace to give . For all accepted models the stationarity residual , the column sum residuals of , and the stationarity residual remain at the level or below. Detailed balance holds exactly for two state chains, as confirmed numerically.
The BKM metric at and the symmetric part are used to compute on the density sector, and compared to the classical Fisher Dirichlet for across random . Over the ensemble the maximal absolute discrepancy is of order , and the maximal relative error is below . This shows that for a broad, physically natural class of random qubit GKLS generators, the eigenbasis density sector is always exactly a two state Fisher Markov model, and the density block of is universally the classical Fisher Dirichlet operator. In the UIH picture this gives a robust, ensemble level confirmation that the irreversible mobility is canonically fixed by GKLS dynamics once the stationary state is known.
E.18 Random qutrit GKLS ensemble and nonreversible density sectors
Script 24_gkls_random_qutrit_density_ensemble.py extends the random ensemble test from qubits to three level systems. Each model has a random Hermitian Hamiltonian on , together with jump operators corresponding to nearest neighbour ladders and with random up and down rates, and three dephasing projectors with random strengths. This family includes genuinely coherent and generically nonreversible open qutrit dynamics.
For each random model the stationary state is computed, and only models with full rank, nondegenerate eigenvalues and nontrivial coherence in the computational basis are retained. In a sample of 10 models all 10 are accepted, with coherence norms between and . In the eigenbasis the GKLS generator yields a three state effective generator on the diagonal subspace. Column sum residuals are of order and stays below , so is a valid Markov generator with stationary .
Unlike the qubit case, three state Markov chains can be nonreversible. The script therefore monitors the maximal detailed balance residual , which reaches values of order in the ensemble. The density sector Markov chains induced by the GKLS models are thus genuinely nonreversible in general.
Despite this, the BKM plus density sector Dirichlet form coincides with the classical Fisher Dirichlet of to numerical precision. Across all accepted models and random perturbations , the maximal absolute discrepancy is of order , with maximal relative error around . This confirms that in coherent, generically nonreversible qutrit GKLS models the eigenbasis density sector is always governed by a classical Fisher Dirichlet structure, and that the antisymmetric circulation in the Markov generator resides entirely in the non gradient part of the real GKLS operator. For UIH this means that even in nonreversible quantum chains the density sector mobility is still the canonical Fisher Dirichlet operator, while all nonreversible circulation is pushed into the channel.
E.19 Constructive Fokker-Planck to Markov to GKLS realisation
Script 25_fp_to_markov_to_gkls_realisation.py closes the loop from continuum Fisher flows to discrete Markov chains and back to a GKLS semigroup. We start from a one dimensional periodic domain with potential
| (E.25) |
which defines a continuum stationary density . The associated overdamped Fokker-Planck equation
| (E.26) |
is a standard example of a Fisher-gradient flow of the relative entropy .
The script discretises on a periodic lattice with points and constructs a reversible nearest neighbour Markov generator with column sums zero and stationary vector via a detailed balance conductance scheme. The stationary distribution is fixed directly from the potential as , and the symmetric nearest neighbour rates are chosen such that . The resulting generator satisfies exactly and , confirming that is stationary on the lattice.
From this we build a diagonal jump GKLS generator on an -dimensional Hilbert space with computational basis . For each off diagonal entry with we introduce a Lindblad operator
| (E.27) |
and form the Lindblad superoperator on via the standard GKLS prescription, with no Hamiltonian part. In this construction, any diagonal density matrix evolves as , and the stationary state satisfies . Extracting the effective density generator by acting on the diagonal basis matrices shows that the entries of match those of to high accuracy, with .
To compare Dirichlet structures we take in the computational basis and build the BKM metric at , which reduces to scalar weights
| (E.28) |
assembled into a diagonal metric on the vectorised space. The metric adjoint is then defined by , and the symmetric part furnishes the irreversible Fisher operator.
For a catalogue of random mass conserving density perturbations the script evaluates both the classical Fisher Dirichlet form
| (E.29) |
with , and the GKLS density sector Dirichlet form
| (E.30) |
Over random perturbations the Dirichlet values lie in the range , with a maximal absolute discrepancy and maximal relative error of order .
This example provides a fully constructive realisation of a Fisher-metriplectic Fokker-Planck flow as the density sector of a GKLS semigroup. Starting from a continuum gradient flow defined by a potential , we discretise to a reversible Markov generator , lift to a diagonal jump GKLS generator , and recover both the Markov dynamics and the Fisher Dirichlet quadratic form on densities from the BKM symmetric part . The irreversible geometry of the continuum Fokker-Planck equation is thus concretely embedded as the density block of the GKLS real generator, giving a direct continuum realisation of the UIH density sector mobility .
Proposition E.1 (Density sector universality of the Fisher operator).
Let be a GKLS generator on a finite dimensional matrix algebra with a full rank stationary state . Let denote the eigenvalue vector of in a basis where , and place the BKM metric at as in Section 6.
Define the metric adjoint by
and the symmetric Fisher operator
For any mass conserving density perturbation with , set
Then the GKLS density sector Dirichlet form
coincides exactly with the classical Fisher Dirichlet form of the reversible Markov generator associated to via its restriction to the diagonal sector,
with and reversible in the geometry.
In particular,
for all mass conserving perturbations . The symmetric Fisher operator on densities is therefore uniquely fixed by the underlying reversible Markov chain and is insensitive to coherent Hamiltonian dressing: many distinct GKLS generators with the same density sector share the same Fisher metriplectic hydrodynamics on densities.
E.20 Coherent Hamiltonian dressing of the Fokker-Planck GKLS chain
Script 26_gkls_coherent_dressing_fp_chain.py starts from the same reversible Markov generator and diagonal jump GKLS generator constructed in script 25, and then adds a genuinely coherent Hamiltonian dressing without changing the density sector hydrodynamics.
On the lattice of sites we define a tight binding Hamiltonian
| (E.31) |
with periodic wrap. The associated Hamiltonian superoperator on is
| (E.32) |
so that corresponds to . The total GKLS generator is then
| (E.33) |
which describes a coherent tight binding evolution on top of the dissipative Fokker-Planck-like chain.
The dissipative generator is built exactly as in script 25 from the reversible Markov generator via diagonal jump operators . Its stationary state is with the Gibbs-like stationary vector induced by the potential. Numerically, , confirming stationarity for the purely dissipative chain. After adding the Hamiltonian part one finds : the full GKLS dynamics now generates off diagonal coherences at , so is no longer the fixed point of , even though it remains the natural centre for the dissipative geometry.
Despite this, when we restrict to the density sector the dynamics is entirely unchanged. Extracting the effective density generators and from and by acting on the diagonal basis matrices shows that both coincide with to numerical precision:
| (E.34) |
and the two effective generators agree exactly at machine precision on the density sector.
To probe the Fisher geometry we again place the BKM metric at the dissipative reference state . The BKM weights are as in script 25, assembled into a diagonal metric on the vectorised space. For each of and we build the metric adjoint and the symmetric part , obtaining and . The density sector GKLS Dirichlet forms are then
| (E.35) |
with , while the classical Fisher Dirichlet is computed from and as in script 25.
Over a catalogue of random mass conserving perturbations , the three Dirichlet values lie in the same range and satisfy
| (E.36) |
with maximal relative errors of order , and to within numerical precision for all tested perturbations. In other words, the tight binding Hamiltonian dressing leaves both the density sector Markov generator and the Fisher Dirichlet quadratic form completely unchanged.
This makes the separation between reversible and irreversible channels explicit in the FP-Markov-GKLS correspondence. The symmetric Fisher operator on densities is fixed by the dissipative chain and is insensitive to coherent Hamiltonian dressing, while the antisymmetric part and the full GKLS stationary state can vary substantially. Many distinct coherent GKLS models thus share the same Fisher-metriplectic hydrodynamics on densities, differing only in the reversible quadrature.
E.21 Density sector tomography for a GKLS generator
Script 27_gkls_density_sector_tomography.py demonstrates that the density sector of a GKLS generator can be reconstructed purely from density responses, and that the resulting Fisher Dirichlet form coincides with both the classical and GKLS constructions.
We first build a random reversible three state Markov generator with a strictly positive stationary distribution using a symmetric conductance construction. Random positive weights are drawn and normalised, and a symmetric non negative conductance matrix is sampled. The off diagonal entries of are then defined by detailed balance,
| (E.37) |
with diagonal entries chosen so that the column sums vanish. In a typical run this yields, for example,
| (E.38) |
with and .
From this generator we build a diagonal jump GKLS superoperator on a three dimensional Hilbert space in the computational basis , introducing Lindblad operators
| (E.39) |
and forming the standard GKLS combination on . The stationary state is confirmed to be fixed by the GKLS evolution with .
The exact density sector generator is then extracted directly from the GKLS dynamics by acting on the diagonal basis matrices and reading off the induced drift of the diagonal entries. The resulting matrix matches the original generator to numerical precision,
| (E.40) |
with and .
To mimic an operational setting where only density responses are accessible, the script performs a simple tomography of the density generator. A catalogue of random full support probability vectors is drawn, each defining a diagonal state . For each probe state the GKLS response is computed and the diagonal drift is extracted. Stacking these into matrices and , with probes, we solve the least squares system
| (E.41) |
for a reconstructed generator , and softly enforce column sum zero by subtracting the mean column sum from each column. The reconstructed generator agrees with both the original and GKLS effective generators at machine precision,
| (E.42) |
with column sums and stationarity residuals again at the level.
To connect with the Fisher-metriplectic structure we place the BKM metric at . The BKM weights are
| (E.43) |
assembled into a diagonal metric on the vectorised space. The metric adjoint is defined by and the symmetric part defines the GKLS Fisher operator.
For a collection of random mass conserving perturbations the classical Fisher Dirichlet form associated with and with the reconstructed ,
| (E.44) |
and the GKLS density sector Dirichlet form
| (E.45) |
are compared. Over random perturbations the maximal absolute and relative discrepancies satisfy
| (E.46) | ||||||
| (E.47) |
with and agreeing to the same level.
This script shows that for a finite GKLS model in the diagonal jump class, the density sector generator can be recovered purely from density responses, and that the Fisher Dirichlet geometry reconstructed from this tomographic coincides with the GKLS Fisher geometry defined by the BKM symmetric part . The irreversible information hydrodynamics on densities is therefore an operational object: it can be inferred from macroscopic density probes alone, and it matches the metriplectic structure seen by the underlying Lindblad generator. In UIH language, the density sector mobility is fully reconstructible from macroscopic probes and agrees with the canonical Fisher Dirichlet operator extracted from GKLS.
E.22 Full Bloch space tomography of a driven damped qubit
Script 28_gkls_full_bloch_tomography.py upgrades the density sector tomography of script 27 to the full operator space. We consider a driven, damped qubit with Hamiltonian and Lindblad operators
| (E.48) |
and GKLS generator
| (E.49) |
We work in the Hermitian operator basis and represent any Hermitian operator by coordinates defined through
| (E.50) |
In this basis the GKLS generator is a real matrix such that .
For a fixed choice , , and the script first constructs the exact generator by acting with on the basis elements and reading off the coordinates of . The resulting matrix is
| (E.51) |
with the first row vanishing exactly, as required by trace preservation.
The same generator is then reconstructed tomographically from state responses. A collection of random Bloch vectors with is sampled, each defining a strictly positive density matrix . For each probe state we compute coordinates and their GKLS derivatives . Stacking these into matrices and with probes, we solve the least squares system
| (E.52) |
giving . The reconstructed generator is
| (E.53) |
with , a trace preserving row norm , and purely real entries within numerical precision. Testing on additional random states shows that the predicted Bloch derivatives match the true GKLS derivatives with maximal residual .
This script demonstrates that for a coherent driven qubit with dissipation the full real GKLS generator on operator space can be reconstructed from state responses in the Pauli basis, not just its density sector. It extends the density level Markov tomography of script 27 to a complete Bloch space tomography that resolves both the irreversible and coherent channels encoded in .
E.23 Metriplectic split of a driven qubit in Bloch coordinates
Script 29_gkls_bloch_metriplectic_split.py takes the coherent, damped qubit model of script 28_gkls_full_bloch_tomography.py and performs a full metriplectic decomposition of the GKLS generator in the Pauli basis with respect to the BKM metric at the stationary state.
In particular the traceless block and the corresponding blocks and realise exactly the real of our metriplectic axioms specialised to this driven, damped qubit.
We work with the same Hamiltonian and jumps
| (E.54) |
with parameters , , , . In the Hermitian basis any density operator is represented as , with coordinates . As in script 28 the GKLS evolution induces a real matrix on the Bloch coordinates via , and for the chosen parameters one recovers
| (E.55) |
with the first row vanishing exactly, confirming trace preservation.
The stationary Bloch vector is obtained by solving the linear fixed point equation . Writing with and using , the script solves with the traceless block of and the driving column from . Numerically one finds
| (E.56) |
corresponding to the stationary state
| (E.57) |
which agrees with the coherently dressed stationary state already seen in script 22. Diagonalising gives eigenvalues and , together with a unitary whose columns are the eigenvectors.
In the eigenbasis of the BKM metric is diagonal on the matrix units , with weights
| (E.58) |
For the present model this yields
| (E.59) |
which is positive definite. Vectorising matrices in the eigenbasis produces a diagonal metric on the four dimensional operator space. The script then transforms each Pauli basis element into the eigenbasis, , vectorises to obtain , and assembles the BKM metric matrix in Bloch coordinates as
| (E.60) |
The resulting matrix is symmetric and strictly positive definite, with
| (E.61) |
and eigenvalues approximately , confirming a well conditioned BKM geometry on the Bloch space.
With this metric in hand the script defines the metric adjoint of by
| (E.62) |
and splits the generator into symmetric and skew parts
| (E.63) |
By construction one has . The diagnostics confirm this to numerical precision, with . The key metriplectic identities also hold:
| (E.64) |
with symmetry residuals and skew residuals . This shows that is the symmetric dissipative channel and the antisymmetric reversible channel with respect to the BKM metric at .
Restricting to the traceless subspace spanned by gives blocks
| (E.65) |
| (E.66) |
The symmetric block is negative definite and captures the irreversible contraction of Bloch vectors toward in the BKM geometry, while the antisymmetric block encodes the effective rotation generated jointly by the Hamiltonian and the coherent part of the dissipator. In particular the off diagonal entries of reflect the precession of the Bloch vector around the driven axis and the shifts induced by the jump operators.
This script therefore provides a fully explicit metriplectic decomposition of a nontrivial qubit GKLS generator: the complex GKLS flow on operators is represented as a real Bloch space generator which splits as with symmetric and antisymmetric in the BKM metric at the stationary state, extending the density sector Fisher-Dirichlet identifications to the full space of observables including coherences.
Appendix F Notation index
For ease of reference we list some of the main symbols used throughout the paper.
| classical probability density on configuration space | |
| phase field in the reversible Schrödinger sector | |
| quantum density matrix evolving under a GKLS equation | |
| stationary Gibbs density | |
| stationary distribution of a finite Markov chain | |
| population of state at time | |
| free energy functional, typically a relative entropy | |
| chemical potential | |
| Fisher information functional or metric | |
| symmetric mobility operator in the metriplectic sector | |
| antisymmetric reversible operator satisfying no work conditions | |
| UIH one-current operator (G gradient, J reversible; i labels the reversible quadrature) | |
| weighted elliptic operators defining geometry | |
| weighted space with metric induced by | |
| classical Markov generator matrix | |
| GKLS generator on density matrices | |
| matrix representation of under vectorisation | |
| transition rate from to in a finite chain | |
| diffusion coefficient in the Fokker Planck equation | |
| potential in the overdamped Langevin setting | |
| Hamiltonian drive and relaxation/dephasing rates in the qubit examples. |
The same symbol may appear in slightly different guises when passing between classical, quantum, and hydrodynamic descriptions. In particular, is always a convex functional whose gradient drives the irreversible dynamics, is always the symmetric mobility that sets the Fisher metric, and is always the antisymmetric no work operator that generates reversible flow.
Acknowledgements
All work was produced independently without affiliation or funding.
Author’s email is contact@nomogenetics.com.
References
- Dunkley [2025-1] J. R. Dunkley. The Converse Madelung Question: Schrödinger Equation from Minimal Axioms. arXiv:2511.03552 (2025).
- Dunkley [2025-2] J. R. Dunkley. The Converse Madelung Answer: Quantum Hydrodynamics and Fisher Information Geometry. 10.5281/zenodo.17643885 (2025).
- Dunkley [2025-5] J. R. Dunkley. Intrinsic Fisher-Kähler Information Geometry. 2025.
- Gorini et al. [1976] V. Gorini, A. Kossakowski and E. C. G. Sudarshan. Completely positive dynamical semigroups of -level systems. Journal of Mathematical Physics 17(5):821-825 (1976).
- Lindblad [1976] G. Lindblad. On the generators of quantum dynamical semigroups. Communications in Mathematical Physics 48(2):119-130 (1976).
- Jordan et al. [1998] R. Jordan, D. Kinderlehrer and F. Otto. The variational formulation of the Fokker-Planck equation. SIAM Journal on Mathematical Analysis 29(1):1-17 (1998).
- Ambrosio et al. [2008] L. Ambrosio, N. Gigli and G. Savare. Gradient Flows in Metric Spaces and in the Space of Probability Measures. Birkhauser, Basel (2008).
- Villani [2009-1] C. Villani. Optimal Transport: Old and New. Springer, Berlin (2009).
- Maas [2011] J. Maas. Gradient flows of the entropy for finite Markov chains. Journal of Functional Analysis 261(8):2250-2292 (2011).
- Mielke [2014] A. Mielke. On the relation between gradient flows and the large-deviation principle, with applications to Markov chains and diffusion. Potential Analysis 41(4):1293-1327 (2014).
- Carlen and Maas [2017] E. A. Carlen and J. Maas. Gradient flow and entropy inequalities for quantum Markov semigroups with detailed balance. Journal of Functional Analysis 273(5):1810-1869 (2017).
- Petz [1996] D. Petz. Monotone metrics on matrix spaces. Linear Algebra and its Applications 244:81-96 (1996).
- Petz and Toth [1993] D. Petz and G. Toth. The Bogoliubov inner product in quantum statistics. Letters in Mathematical Physics 27(3):205-216 (1993).
- Villani [2009-2] C. Villani. Hypocoercivity. Memoirs of the American Mathematical Society 202(950) (2009).
- Ozawa et al. [2024] T. Ozawa, Y. Kumagai, T. Kume, T. Oda and H. Kadowaki. Observation of resonant tunneling of proton from octahedral to tetrahedral sites in palladium. Journal of Physics and Chemistry of Solids 185:111741 (2024).
- Volpe and Volpe [2013] G. Volpe and G. Volpe. Simulation of a Brownian particle in an optical trap. American Journal of Physics 81(3):224-230 (2013).
- Abraham et al. [2019] H. Abraham et al. Qiskit: An open-source framework for quantum computing. Zenodo, doi:10.5281/zenodo.2562111 (2019).
- Frieden [1998] B. R. Frieden. Physics from Fisher Information: A Unification. Cambridge University Press (1998).
- Wolfram [2002] S. Wolfram. A New Kind of Science. Wolfram Media (2002).