A translator who never misses a whisper ✛

Original: "Computational framework for quantum state tomography of spin ensembles"
· Jiwoo Seo, Vesna F. Mitrović
arXiv:2609.11868v1 · 2026-09-10 · CC BY 4.0 · 2 min
The open-source Tomography-NMR package turns raw NMR spectra into density matrices with accuracy above 99%.

A courtroom interpreter catches every word, every breath, every whisper in the corner. They have no right to miss a detail, because the meaning of the whole conversation depends on it. Quantum state tomography works in much the same way: it reconstructs the full quantum state of a system from a set of measurements, leaving no detail out.

For a long time this "interpreter" was locked inside the black box of commercial software. Researchers saw only the result, never the process itself: how exactly raw data from the instrument turn into a density matrix, the mathematical description of a quantum state. The Tomography-NMR package opens that box. It takes the time-domain signal, applies a Fourier transform, corrects the phase, and then extracts peak intensities from the spectrum.

Three ways to extract peaks: height at a single point, numerical integration over a fixed window, and systematic parameter sweeping for maximum accuracy. The first two work blind, while the third checks itself against a known reference.

Every doublet, that is, a pair of peaks from coupled spins, encodes two numbers: the sum and the difference of the left and right peaks. These numbers are directly tied to the coefficients in the expansion of the state over the product operator basis. The resulting matrix is then subjected to physical constraints: it must be Hermitian, have a trace equal to one, and be positive semidefinite. That is how a correct quantum state is born from a spectrum.

A test on 20 experimental states prepared on two phosphorus-31 nuclei in an ADP molecule showed an average accuracy from 0.975 to 0.995. Basis states were reconstructed with an accuracy of 0.995, outputs of the Hadamard gate with 0.992, and Bell states with 0.980. Optimized integration for individual states reaches 0.998, whereas methods that do not know the target stay around 0.98.

Errors on the density matrix elements are estimated at the level of 10⁻³–10⁻⁴. It is like an interpreter who gets one word in ten thousand wrong.

In the past, checking someone else's result meant trusting closed-source code or reproducing it from scratch. Now anyone can download Tomography-NMR, run the same data through it, and compare the density matrices. For experimental quantum information science this means a shift in practice: transparency and reproducibility become the norm rather than the exception. Students get a hands-on teaching tool, and researchers get a common language for comparing results.

The interpreter no longer sits behind soundproof glass. They have stepped into the room, and anyone can hear exactly how they turn whispers into meaning.

The package is called Tomography-NMR and is distributed under the MIT license, which allows free use in academic and commercial projects. Validation was performed on the ADP molecule, where two phosphorus-31 nuclei are coupled through a J interaction of about 22 Hz.

\hat{\rho} = \sum_{i,j \in \{o,x,y,z\}} c_{ij} \hat{B}_{ij}
The density matrix is represented as a linear combination of 16 basis operators with coefficients c_ij.
\mathrm{Re}(L+R) = \mathrm{Tr}(\hat{\rho}\tilde{\hat{I}}^1_x), \quad \mathrm{Re}(L-R) = \mathrm{Tr}(\hat{\rho}2\tilde{\hat{I}}^1_x \tilde{\hat{I}}^2_z)
The spectral peak intensities directly encode the expectation values of operators, which makes it possible to extract the expansion coefficients.
F_{\mathrm{proj}}(\hat{\rho}_m, \hat{\rho}_{th}) = \frac{\mathrm{Tr}(\hat{\rho}_m \hat{\rho}_{th})}{\sqrt{\mathrm{Tr}(\hat{\rho}_m^2) \mathrm{Tr}(\hat{\rho}_{th}^2)}}
The normalized inner product between the measured and theoretical density matrices, used to assess reconstruction quality.
Concepts
tomography NMR spectroscopy quantum snapshot quantum algorithm density matrix mri numerical simulation image reconstruction decoherence free subspace renormalization decoder rydberg atom simulator computed tomography coherent state two level system neural network polarimetry kernel density estimation nuclear spin spectroscopy mixed state error suppression