---
title: Dynamical Quantum Tomography
url: https://www.emergentmind.com/topics/dynamical-quantum-tomography
type: topic
---

# Dynamical Quantum Tomography

Dynamical quantum tomography denotes a family of quantum-characterization protocols in which temporal structure is part of the inverse problem. In the most common formulation, an unknown state is subjected to a known evolution and probed by a fixed or restricted measurement setup at multiple times, so that dynamics enlarge the effective measurement span and can replace a large static tomographically complete measurement family. In other formulations, the unknown state itself changes during data acquisition and is tracked online, or the reconstructed object is not a state but a process, a Liouvillian, a Hamiltonian, a reduced subsystem state, or a family of dynamical correlation functions [1605.06786][2205.06389][2504.10393].

## 1. Formal definitions and informational completeness

A standard finite-dimensional formulation uses a POVM \(P=(Q_1,\dots,Q_m)\) together with a known Heisenberg-picture evolution \(\mathcal T\). The associated dynamical measurement scheme is
\[
\mathcal T^l(P)=\Big((Q_1,\dots,Q_m),(\mathcal T(Q_1),\dots,\mathcal T(Q_m)),\dots,(\mathcal T^{l-1}(Q_1),\dots,\mathcal T^{l-1}(Q_m))\Big),
\]
with measurement data
\[
h_{\mathcal T^l(P)}(x)=\big(\operatorname{tr}(Q_i(\mathcal T^\dagger)^j(x))\big)_{\substack{i=1,\dots,m\\ j=0,\dots,l-1}}.
\]
Within this framework, a measurement scheme is \(\mathcal R\)-complete when it is injective on a restricted state set \(\mathcal R\), and informationally complete when it is injective on the full state space. For a feasible unitary \(U\in U(n)\), almost all POVMs with \(n+1\) outcomes yield informationally complete schemes \(\mathcal T_U^n(P)\); with prior information \(\mathcal R\) represented by a semi-algebraic difference set \(\mathcal D\), the condition \(l(m-1)>\dim\mathcal D\) guarantees stable \(\mathcal R\)-completeness for almost all \(m\)-outcome POVMs, while feasible invertible CPTP dynamics remove the unitary lower bound \(m\ge n\) and permit two-outcome schemes [1605.06786].

A complementary control-theoretic formulation defines the observable subspace
\[
O=\operatorname{span}\{\Phi_t(X_i)\mid X_i\in X,\ t\in T\},
\]
for a known Heisenberg evolution \(\Phi_t\) and an available observable family \(X=\{X_i\}\). Full dynamical quantum state tomography is feasible exactly when
\[
O=B(H).
\]
For Markovian dynamics, this reduces to finite Krylov-type tests:
\[
O=\operatorname{span}\{L^t(X_i)\mid X_i\in X,\ t=0,\dots,k^*\}=B(H)
\]
in continuous time, or
\[
O=\operatorname{span}\{\Phi_t(X_i)\mid X_i\in X,\ t=0,\dots,k^*\}=B(H)
\]
in discrete time, with deterministic rank tests built from observability matrices \(O_c\) and \(O_d\) [2509.24636]. This same idea appears in open-system stroboscopic tomography for phase-damping channels, where the evolution
\[
\rho(t)=D(t)\circ \rho(0)
\]
is decomposed as
\[
D(t)=\sum_{k=1}^{\mu}\lambda_k(t)A_k,
\]
so that repeated measurements at times \(t_j\) yield a linear system for the projections \(\operatorname{Tr}\{(Q_i\circ A_k^T)\rho(0)\}\); reconstruction requires \(\operatorname{rank}[\lambda_k(t_j)]=\mu\) and
\[
\bigoplus_{i=1}^r \mathcal S(Q_i;A_1,\dots,A_\mu)=B_*(\mathcal H)
\]
with \(\mathcal S(Q_i;A_1,\dots,A_\mu)=\operatorname{Span}\{Q_i\circ A_1^T,\dots,Q_i\circ A_\mu^T\}\) [1509.09318].

A more specialized channel-based construction uses a time-dependent average channel
\[
\Psi_t(\rho)=\sum_{i=0}^{d^2-1}\mu_i(t)M_i\rho_0 M_i^\dagger
\]
built from Weyl–Heisenberg operators \(M_\alpha\). Measuring one fixed projector \(|\phi\rangle\langle\phi|\) at \(d^2\) time instants gives a linear system whose unknowns are
\[
\operatorname{tr}[\rho_0\,M_i^\dagger|\phi\rangle\langle\phi|M_i]
=\operatorname{tr}[\rho_0|\phi_i\rangle\langle\phi_i|].
\]
If the Gram matrix
\[
\mathcal A=\big[|\langle\phi_j|\phi_k\rangle|^2\big]_{0\le j,k\le d^2-1}
\]
is nonsingular, the projectors \(\{|\phi_\alpha\rangle\langle\phi_\alpha|\}\) are informationally complete, and the corresponding normalized family is an IC-POVM [2402.05431].

## 2. Online state tracking and evolving-state tomography

A distinct branch of dynamical tomography addresses the case in which the unknown state changes while data are being acquired. In this setting, batch estimators such as maximum-likelihood estimation and least-squares tomography are poorly matched to the experiment because they require a tomographically complete data set before returning an estimate. Matrix-exponentiated gradient tomography instead updates the estimate after each measurement record:
\[
\hat{\rho}_{t+1}
=
\frac{\exp(\log(\hat{\rho_t})-\eta_t\nabla L_t)}
{\operatorname{tr}\!\left[\exp(\log(\hat{\rho_t})-\eta_t\nabla L_t)\right]},
\]
with single-measurement loss
\[
L_t=(\operatorname{tr}(\hat\rho_t X_t)-y_t)^2,
\qquad
\nabla L_t=2(\operatorname{tr}(\hat\rho_t X_t)-y_t)X_t,
\]
and, in the photonic implementation, a multi-outcome version
\[
L_t=\sum_{i=1}^d\left(\operatorname{tr}(\hat\rho_tX_t^{(i)})-y_t^{(i)}\right)^2.
\]
Because the update acts through \(\log\hat\rho_t\) and matrix exponentiation, positivity and unit trace are preserved at every step, so the estimate remains a valid density matrix without ad hoc projection [2205.06389].

The experimental realization was carried out on a photonic qutrit encoded in transverse spatial modes. Measurements were chosen sequentially from informationally complete families based on mutually unbiased bases or generalized Pauli operator measurements, and the evolving target state was taken as
\[
\ket{\psi_t}=\exp(-i\sigma\omega t)\ket{\psi_0},
\]
with stationary, structured, and random-Hermitian trajectories. With a constant learning rate \(\eta=5\), the protocol reached infidelity below \(10\%\) in about 4 iterations in high-count qutrit experiments and in roughly 13–16 iterations for low-count MUB tomography, while mean infidelities remained around \(3.4\%-5.6\%\), corresponding to fidelities around \(94\%-96.6\%\). Under added background noise up to approximately \(N_{\text{back}}\approx 2\) kHz, performance stayed comparable to the no-added-noise case within uncertainty bars, and the authors interpret this as corresponding to \(\mathrm{SNR}\approx 0.05\) [2205.06389].

A related pure-state line assumes a known time-independent Hamiltonian and projective intensity measurements, treating reconstruction as phase retrieval under Schrödinger evolution and aiming to decrease the number of distinct projectors by exploiting the known unitary dynamics [1507.00715]. This suggests that “dynamical” can refer either to online tracking of a drifting state or to using known evolution as an informational resource even when the target state itself is static.

## 3. Process, channel, and Liouvillian tomography

Dynamical tomography frequently targets the evolution law itself rather than an instantaneous state. One early route is weak-measurement-based process tomography, where a process
\[
\mathcal E:\mathcal L(\mathcal H_{in})\to\mathcal L(\mathcal H_{out})
\]
is expanded as
\[
\mathcal{E}(\Omega_{in})
=
\sum_{i_1,i_2=1}^{d_{in}}
\sum_{i_3,i_4=1}^{d_{out}}
\chi_{i_1 i_2 i_3 i_4}
\langle \alpha_{i_2}|\Omega_{in}|\psi_{i_1}\rangle
|\beta_{i_3}\rangle\langle \phi_{i_4}|,
\]
and each coefficient \(\chi_{i_1i_2i_3i_4}\) is linked directly to a weak-value-like quantity \(X^{i\hat A\hat Bf}\). In the proposed scheme, every process parameter is determined from only five experimental values—four pointer correlations \(r_1,r_2,r_3,r_4\) and one post-selection probability \(p_{f|i\hat A\hat B}\)—and complete tomography requires only \(d_{in}\) setups, with product input states sufficient even for multiparticle processes [1309.5780].

Another route keeps the object of reconstruction at the channel level but uses tomography to expose dynamical structure invisible to scalar coherence times. In superconducting-qubit dynamical decoupling experiments, single-qubit QPT in the \(\chi\)-matrix representation
\[
\rho_f=\sum_{nm}\chi_{nm}E_n\rho_iE_m^\dagger,
\qquad
\{E_m\}=\{I,\sigma_x,i\sigma_y,\sigma_z\},
\]
showed that dynamical decoupling suppresses dephasing but does not suppress spontaneous emission, and that pulse imperfections generate coherent residual rotations under \(XY\)-4. Robust sequences such as \(XY\)-8, \(XY\)-16\), KDD, and UR20 remove the oscillatory behavior caused by control errors, showing that process tomography can distinguish dephasing suppression, persistent relaxation, and coherent pulse-induced artifacts within a single reconstructed dynamical map [2006.10585].

More recent work moves from maps to generators. Room-temperature qutrit process tomography in a \(^{87}\mathrm{Rb}\) vapor reconstructs finite-time process matrices \(\hat P(t)\) in the Bloch–Fano basis and then infers an effective total relaxation superoperator \(\hat R_T\), separating residual Zeeman Hamiltonian terms, dephasing, and isotropic relaxation in a realistic noisy ensemble [2508.19634]. Lindblad-like quantum tomography then generalizes the semigroup paradigm to time-local non-Markovian maps by maximizing a multi-snapshot likelihood over time-local master equations with possibly negative decay rates,
\[
\dot\rho
=
-i[H(t),\rho]
+
\sum_n \gamma_n(t)
\left(
L_n(t)\rho L_n^\dagger(t)-\frac12\{L_n^\dagger(t)L_n(t),\rho\}
\right),
\]
and shows explicitly, for single-qubit dephasing, why multiple temporal snapshots are required once the dynamics is not time homogeneous [2403.19799]. Quantum Liouvillian Tomography pushes the same program to multi-qubit open dynamics by combining gradient-based QPT with regression over derivatives of Pauli-string probabilities,
\[
\frac{d}{dt}p_{l|ij}(t)
=
\operatorname{Tr}\!\left[\mathcal L_t\Lambda_t(\rho_i)M_{lj}\right],
\]
to reconstruct a time-local Liouvillian
\[
\mathcal L_t(\cdot)
=
-i[H(t),\cdot]
+
\sum_\mu \gamma_\mu(t)
\left(
J_\mu(t)\,\cdot\,J_\mu^\dagger(t)
-\frac12\{J_\mu(t)J_\mu^\dagger(t),\cdot\}
\right)
\]
and to detect non-Markovianity through negative canonical rates in idle two-qubit superconducting dynamics [2504.10393].

## 4. Many-body, subsystem, and correlation-function tomography

In many-body settings, dynamical tomography is often feasible only because the target states occupy a structured submanifold of Hilbert space. Matrix product state tomography exploits the fact that out-of-equilibrium states of one-dimensional systems with finite-range interactions remain efficiently approximable by MPS for any fixed evolution time. The protocol reconstructs local reduced density matrices on contiguous blocks of size \(k\), requiring at most \(3^k\) local Pauli bases per \(k\)-site block for qubits, and then fits a global MPS with a certified fidelity lower bound
\[
\langle \Psi_c^k|\rho_{\mathrm{lab}}|\Psi_c^k\rangle\ge F_c^k.
\]
In a trapped-ion simulator, this enabled reconstruction of dynamical states of up to 14 spins, with \(F_c^1=0.98\pm0.01\) for the initial product state in the 8-spin experiment and \(F_c^3>0.8\) up to \(t=3\) ms for triplet-block tomography, while the eventual collapse of the certificate tracked the physical spread of correlations beyond the chosen block size [1612.08000].

A more targeted subsystem variant appears in lattice gauge theory, where quench dynamics in the Schwinger model were characterized not only through Loschmidt echoes and non-equal-time correlators but also through time-resolved entanglement tomography of reduced density matrices. The reconstructed subsystem state
\[
\rho_A(t)=\operatorname{Tr}_B[\rho(t)]
\]
was analyzed through the second Rényi entropy
\[
S_A^{(2)}(t)=-\log_2\{\operatorname{Tr}_A(\rho_A^2(t))\}
\]
and through a fitted entanglement Hamiltonian ansatz \( \rho_A(t)=e^{-H_A(t)} \), making it possible to extract time-dependent entanglement spectra and entanglement Hamiltonians during a dynamical quantum phase transition [2210.03089].

The same logic extends from states to families of observables. Fermionic-Adapted Shadow Tomography reformulates dynamical commutators and anti-commutators into expectation values compatible with shadow methods, so that many correlators can be estimated simultaneously with at most two-copy measurements and uncontrolled Hamiltonian simulation. For example, the commutator
\[
\mathrm{tr}\!\left(\rho\,[e^{iHt}c_i^\dagger c_j e^{-iHt},\,c_k^\dagger c_l]\right)
\]
is reduced to ordinary expectations on three effective evolved states, while retarded Green’s functions
\[
G^R_{ab}(t)=-i\Theta(t)\,\mathrm{tr}\!\left(\rho\{c_i(t),c_j^\dagger\}\right)
\]
are treated through probabilistically prepared branches \(\rho_\pm\). The resulting sample complexities improve from \(\mathcal O(n^4\log n/\epsilon^2)\) to \(\mathcal O(n^3\log n/\epsilon^2)\) or \(\mathcal O(n^2\log n/\epsilon^4)\) in the commutator case, and from \(\mathcal O(n^2\log n/\epsilon^2)\) to \(\mathcal O(n\log n/\epsilon^4)\) in the anti-commutator regime \(n\ge 1/\epsilon^2\) [2508.03192].

## 5. Spectroscopic and Hamiltonian-learning formulations

In ultrafast molecular spectroscopy, dynamical tomography is formulated as the reconstruction of a time-dependent molecular density operator from measured angularly resolved observables. A maximum-entropy approach represents the reconstructed state as
\[
\hat{\rho}_{ME}(\vec{\lambda},\hat f)
=
\frac{1}{Z(\lambda_1,\ldots,\lambda_k)}
\exp\!\left\{-\sum_k \lambda_k \hat f_k\right\},
\]
or, for the NH\(_3\) case study,
\[
\hat{\rho}_{ME}(\Omega,\vec{\lambda}(t))
=
\frac{
\exp\left\{
-\lambda_{00}(t)\hat\beta_{00}
-\lambda_{20}(t)\hat\beta_{20}
-\lambda_{40}(t)\hat\beta_{40}
\right\}
}{
Z(\lambda_{00}(t),\lambda_{20}(t),\lambda_{40}(t))
}.
\]
The measured anisotropy coefficients
\[
\beta_{LM}(\epsilon;t)
=
\sum_{ij}\sum_{KQS}
C^{LM}_{KQS}(i,j;\epsilon)\,A^K_{QS}(i,j;t)
\]
link experimental photoelectron observables to molecular angular distribution moments, enabling time-resolved reconstruction of the lab-frame density matrix \(\hat\rho(\Omega,t)\), visualization of charge migration, and extraction of the electronic subsystem entropy \(S(\rho)=-\mathrm{Tr}(\rho\ln\rho)\) [2407.16630].

A different spectroscopic variant targets the generator of the dynamics rather than the state: Hamiltonian learning in engineered quantum magnets from spatially resolved dynamical response. There the basic data are local spectral functions
\[
S^{aa}_n(\omega)
=
\langle \mathrm{GS}|S_n^a\,\delta(\omega-H+E_{\mathrm{GS}})\,S_n^a|\mathrm{GS}\rangle
=
\sum_\alpha
|\langle \alpha|S_n^a|\mathrm{GS}\rangle|^2
\delta(\omega-E_\alpha+E_{\mathrm{GS}}),
\]
measured across several impurity configurations and related to tunneling spectroscopy through
\[
S^{aa}_n(\omega)\sim \frac{d^2I}{dV^2},
\qquad
\int_0^V S(\omega)\,d\omega \sim \frac{dI}{dV}.
\]
A supervised network then infers couplings \((J_2,J_{\mathrm Z},J_3,J_{\mathrm{DMI}})\) in an effective spin Hamiltonian with Heisenberg exchange, anisotropy, and Dzyaloshinskii–Moriya interaction. In the noiseless case, the reported fidelities are \(0.99\) for \(J_2\), \(J_{\mathrm Z}\), and \(J_3\), and \(0.96\) for \(J_{\mathrm{DMI}}\); with noise strength \(\chi=1.0\), the multi-impurity protocol improves the corresponding fidelities to \(0.98\), \(0.98\), \(0.98\), and \(0.89\), respectively [2510.18613].

## 6. Conceptual boundaries, limitations, and open problems

The term “dynamical quantum tomography” is not uniform across the literature. Some work uses “dynamical” in the modern sense of exploiting physical time evolution to gain informational completeness or to track a changing state, while other work casts the reconstruction algorithm itself as a dynamical system. In the latter category, pure-state tomography can be reformulated through a physical imposition operator whose fixed points are the states compatible with the measured distributions; in that framework, multiple compatible reconstructions appear as bifurcations, but the underlying tomography problem remains static rather than time-resolved [1401.1481]. A related distinction concerns the inferred object: state-tracking methods infer \(\rho_t\), process and Liouvillian methods infer \(\mathcal E_t\) or \(\mathcal L_t\), and correlation-function or Hamiltonian-learning methods infer selected observables or model parameters rather than a full density matrix [2205.06389][2504.10393][2510.18613].

The main technical limitations are also formulation-dependent. Online MEG tomography demonstrates robust empirical tracking but does not derive new nonstationary tracking bounds, and its responsiveness depends on the learning-rate tradeoff between noise averaging and adaptation [2205.06389]. Time-local generator methods such as Lindblad-like tomography and QLT assume that the reduced dynamics admit a time-local master equation; they are therefore broader than semigroup tomography but narrower than fully general process-tensor reconstructions, and they are especially sensitive to derivative estimation and discretization choices [2403.19799][2504.10393]. Structured many-body protocols gain efficiency only while entanglement and correlation length remain controlled; in MPS tomography, the block size \(k\) must eventually grow with time, so the method is efficient in system size at fixed time but not generically efficient in evolution time itself [1612.08000].

Taken together, these results suggest three persistent design principles. First, temporal diversity can substitute for measurement diversity, but only when the known dynamics genuinely expands the observable span. Second, structural priors—low entanglement, restricted rank, operator sparsity, time-locality, or Hamiltonian parametrization—are not incidental conveniences but the main reason scalable dynamical tomography is presently possible. Third, open-system dynamics can enlarge tomographic reach in ways that unitary dynamics cannot: for time-homogeneous Markovian unitary evolution, a single nontrivial observable is insufficient when \(d>2\), whereas generic open dynamics can make even single-observable tomography feasible [2509.24636]. This suggests that future progress will likely come less from a single universal protocol than from a continued stratification of tomography problems by the dynamical structure they exploit [1605.06786][2509.24636].

Source: https://www.emergentmind.com/topics/dynamical-quantum-tomography