How to Hear the Quantum 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 Tomography-NMR program turns signals from atoms into a precise description of a quantum state.

When you listen to an echo in the mountains, the time it takes for the sound to return tells you how far away the cliff is. Quantum physicists do something similar with atoms: they send out a radio pulse, and from the returning "echo" they reconstruct what state the atoms were in. This process is called quantum state tomography.

For a long time, the software for this kind of reconstruction was closed. Researchers saw only the final result, never the path itself: how a messy signal becomes a density matrix, that is, an exact mathematical description of the quantum state. The new Tomography-NMR package opens up that path. It takes the raw signal, turns it into a spectrum using a Fourier transform, and then measures the heights of the peaks.

A spectrum looks like a graph with two humps for every pair of coupled atoms. The sum and the difference of the heights of these humps give you all the information you need.

After the measurements, the program imposes physical constraints: the state has to be real, not a mathematical fantasy. Testing on 20 different states showed accuracy above 99%. It is like a translator who makes a mistake in one word out of ten thousand.

In the past, to check someone else's result you had to blindly trust closed code. Now anyone can download Tomography-NMR, run the same data, and compare the answers. This makes science more honest and teaching students more visual. The quantum whisper can finally be heard without intermediaries.

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