GKSL Dilation for Quantum Fragmentation Dynamics

• GKSL Dilation is an operator-theoretic embedding that maps classical pure-jump master equations for fragmentation dynamics into the quantum Lindblad framework.
• It constructs an explicit Hilbert space representation with nonlocal Lindblad jump operators to exactly recover mass-conserving classical dynamics.
• The approach leverages a mass-weighted interpretation of fragmentation kernels, providing rigorous spectral analysis and inversion applications in non-equilibrium systems.

The Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) dilation provides a rigorous operator-theoretic construction to embed classical pure-jump master equations, arising in self-similar fragmentation dynamics, into the quantum GKSL (Lindblad) framework. This dilation is exact under a mass-weighted interpretation of the fragmentation daughter kernel and is formulated on a Hilbert space whose diagonal sector recovers the classical conservative master equation without approximation, elucidating both the structure and spectral properties of evolving particle size distributions (PSDs) (Segura, 10 Jan 2026).

1. Foundations: The Pure-Jump Master Equation in Log-Size Space

The dynamical evolution of mass distributions $m(x,t)=x\,n(x,t)$ for fragmenting particles with size $x>0$ is often encoded via population balance equations (PBEs). Introducing the logarithmic size variable $\xi = \ln x$, the normalized mass-fraction density is

$$p(\xi,t) = \frac{m(e^\xi,t)\,e^\xi}{M}, \quad M = \int_0^\infty m(x,t)\,dx,$$

satisfying $\int p(\xi,t)\,d\xi = 1$. For a self-similar fragmentation kernel

$$S(x) = k\,x^\alpha, \qquad b_m(x, y) = \frac{1}{y} \kappa\bigl(\frac{x}{y}\bigr), \qquad \int_0^1 \kappa(z)\,dz = 1,$$

the corresponding pure-jump master equation reads

$$\partial_t p(\xi,t) = (Gp)(\xi) = -\lambda(\xi)\,p(\xi,t) + \int_0^\infty \lambda(\xi+u) K(u) p(\xi+u, t)\,du,$$

where

$$\lambda(\xi) = S(e^\xi), \qquad K(u) = e^{-u} \kappa(e^{-u}), \qquad \int_0^\infty K(u)\,du = \int_0^1 z\,\kappa(z)\,dz = 1.$$

$G$ acts as the generator of a conservative pure-jump semigroup. The representation in log-size space is pivotal, as it enables a direct correspondence with the Hilbert space framework required for GKSL dilation.

2. Hilbert Space Structure and Density Operators

The embedding is formulated on the Hilbert space $\mathcal{H} = L^2(\mathbb{R}, d\xi)$ with basis $\{|\xi\rangle\}$, enforcing $\langle \xi | \xi' \rangle = \delta(\xi - \xi')$. Dynamics are described via a density operator $\rho(t)$ on $\mathcal{H}$, whose diagonal sector $p(\xi, t) = \langle \xi | \rho(t) | \xi \rangle$ retains the interpretation as mass-fraction density, mapping the classical stochastic evolution into the quantum framework in a measure-theoretically faithful manner.

3. Construction of the Exact Lindblad (GKSL) Superoperator

The quantum generator is of purely dissipative form ($H = 0$), although a commuting Hermitian component may be added without affecting the diagonal. The central object is the nonlocal Lindblad jump operator $L : \mathcal{H} \rightarrow \mathcal{H}$: $$L(\xi, \eta) = \sqrt{\lambda(\eta)}\,\sqrt{K(\eta-\xi)}\,\Theta(\eta-\xi),$$ where $\Theta$ is the Heaviside step function. For any $\psi \in \mathcal H$,

$$(L\psi)(\xi) = \int_\xi^\infty \sqrt{\lambda(\eta)K(\eta-\xi)}\,\psi(\eta)\,d\eta.$$

The master equation for the density operator is

$$\frac{d\rho}{dt} = \mathcal{L}[\rho] = L \rho L^\dagger - \frac{1}{2}\{L^\dagger L, \rho\}.$$

On diagonal density matrices, this recovers the original pure-jump generator due to the kernel properties and conservation relations, contingent on the mass-weighted normalization $\int_0^1 z\,\kappa(z)\,dz = 1$.

4. Exactness: Mass-Weighted Kernels and Faithfulness of Embedding

The exactness of the GKSL dilation crucially depends on interpreting the daughter kernel in its mass-weighted form,

$$b_m(x, y) = \frac{1}{y}\,\kappa\left(\frac{x}{y}\right),$$

which translates to $K(u) = e^{-u}\,\kappa(e^{-u})$ in log-size space. The condition $\int_0^\infty K(u)\,du = 1$ is equivalent to demanding $z\,\kappa(z)$ be a probability measure, ensuring the operator-level correspondence between the quantum and classical dynamics remains strictly valid for arbitrary jump sizes. No diffusion or small-jump approximation is invoked; the mapping handles all nonlocal jumps explicitly (Segura, 10 Jan 2026).

5. Diagonal Projection and Recovery of the Classical Generator

Assuming $\rho$ remains diagonal, $\rho(\eta, \eta') = p(\eta)\,\delta(\eta-\eta')$, the Lindblad superoperator's diagonal recovers the classical master equation: $$\langle \xi | L\rho L^\dagger |\xi\rangle = \int_0^\infty \lambda(\xi+u) K(u) p(\xi+u)\,du, \qquad \langle \xi | \tfrac{1}{2}\{L^\dagger L, \rho\} |\xi\rangle = \lambda(\xi)p(\xi).$$ Thus, $\langle\xi| \mathcal{L}[\rho] | \xi \rangle$ coincides with the gain–loss terms of $G$, establishing that the GKSL dilation, when projected onto the diagonal, generates exactly the intended classical evolution.

6. Non-Uniqueness, Generality, and Quantum Trajectory Connections

This operator-level GKSL embedding is not unique. Alternative dilations can be constructed by unitary rotations or by adding Hermitian $H_0$ terms commuting with the diagonal. Off-diagonal coherences will generally arise unless the dynamics are subject to rapid dephasing in the $\xi$ basis, wherein the diagonal sector forms a closed, effective description. The methodology parallels approaches in quantum trajectory and unraveling theories, but provides a fully explicit construction for a class of physically motivated population-balance equations, repurposing quantum stochastic tools for non-equilibrium statistical applications (Segura, 10 Jan 2026).

7. Applicability, Implications, and Spectral Framework

The GKSL dilation enables new spectral and inverse approaches to particle size distributions evolving under fragmentation, including routes to Schrodinger-type operator analysis via symmetry and log-size transformations. A plausible implication is that spectral dictionaries and inversion techniques, such as those employing parametric time series or direct steady-state inversion, gain a rigorous footing through the quantum dilation structure. Stationary distributions in the genuinely non-Hermitian regime are represented as biorthogonal products of left and right ground states, reflecting the underlying quantum-like structure. The normalization $\int_0^1 z\,\kappa(z)\,dz=1$ is the key criterion ensuring the faithfulness and applicability of this exact operator embedding in modeling complex fragmentation dynamics (Segura, 10 Jan 2026).