A key task in many quantum-computing applications, e.g., quantum simulation and quantum state tomography (QST), is to partition an arbitrary set of operators into mutually commuting subsets for efficient measurements. However, brute-force approaches to this task quickly become intractable as the number and dimensionality of operators grow. Here, we reformulate operator partitioning as a graph-coloring (GC) problem and develop an efficient computational framework to solve it, balancing accuracy and efficiency. Our framework enables leveraging a range of GC algorithms, which we benchmark for operator partitioning. Then, we demonstrate their utility in optimizing QST experiments, where determining non-overlapping data acquisition settings for QST is a major challenge, and prioritizing among these settings, i.e., selecting the experiments that provide the most information. We further show how to perform these experiments by synthesizing Clifford circuits for joint measurement of commuting Pauli operators in multi-qubit systems. We validate our framework across multi-qubit (up to five qubits), multi-qutrit (up to three qutrits), and hybrid qubit-qutrit systems. Our results show that heuristic GC methods substantially reduce the number of required measurement settings for QST and enable priority-based scheduling that maximizes the information gain per experiment. The optimization converges within minutes on a student-grade laptop, providing speedups of several orders of magnitude over brute-force methods already for these relatively small quantum systems. This demonstrates the potential of GC heuristics as a scalable and practical tool for characterization of noisy intermediate-scale quantum devices. We have made the Python implementation of our GC framework to optimize and schedule QST experiments publicly available at https://github.com/ssm8015/QST_GT.git.
A key task in many quantum-computing applications, e.g., quantum simulation and quantum state tomography (QST), is to partition an arbitrary set of operators into mutually commuting subsets for efficient measurements. However, brute-force approaches to this task quickly become intractable as the number and dimensionality of operators grow. Here, we reformulate operator partitioning as a graph-coloring (GC) problem and develop an efficient computational framework to solve it, balancing accuracy and efficiency. Our framework enables leveraging a range of GC algorithms, which we benchmark for operator partitioning. Then, we demonstrate their utility in optimizing QST experiments, where determining non-overlapping data acquisition settings for QST is a major challenge, and prioritizing among these settings, i.e., selecting the experiments that provide the most information. We further show how to perform these experiments by synthesizing Clifford circuits for joint measurement of commuting Pauli operators in multi-qubit systems. We validate our framework across multi-qubit (up to five qubits), multi-qutrit (up to three qutrits), and hybrid qubit-qutrit systems. Our results show that heuristic GC methods substantially reduce the number of required measurement settings for QST and enable priority-based scheduling that maximizes the information gain per experiment. The optimization converges within minutes on a student-grade laptop, providing speedups of several orders of magnitude over brute-force methods already for these relatively small quantum systems. This demonstrates the potential of GC heuristics as a scalable and practical tool for characterization of noisy intermediate-scale quantum devices. We have made the Python implementation of our GC framework to optimize and schedule QST experiments publicly available at https://github.com/ssm8015/QST_GT.git.
Abstract We present a framework that formulates the quest for the most efficient quantum state tomography (QST) measurement set as an optimization problem which can be solved numerically, where the optimization goal is the maximization of the information gain. This approach can be applied to a broad spectrum of relevant setups including measurements restricted to a subsystem. To illustrate the power of this method we present results for the six-dimensional Hilbert space constituted by a qubit–qutrit system, which could be realized e.g. by the 14N nuclear spin-1 and two electronic spin states of a nitrogen-vacancy center in diamond. Measurements of the qubit subsystem are expressed by projectors of rank three, i.e. projectors on half-dimensional subspaces. For systems consisting only of qubits, it was shown analytically that a set of projectors on half-dimensional subspaces can be arranged in an informationally optimal fashion for QST, thus forming so-called mutually unbiased subspaces. Our method goes beyond qubits-only systems and we find that in dimension six such a set of mutually-unbiased subspaces can be approximated with a deviation irrelevant for practical applications.
Quantum state tomography is essential for characterizing quantum systems, but it becomes infeasible for large systems due to exponential resource scaling. Overlapping tomography addresses this challenge by reconstructing all k-body marginals using few measurement settings, enabling the efficient extraction of key information for many quantum tasks. While optimal schemes are known for qubits, the extension to higher-dimensional qudit systems remains largely unexplored. Here, we investigate optimal qudit overlapping tomography, constructing local measurement settings from generalized Gell-Mann matrices. By establishing a correspondence with combinatorial covering arrays, we present two explicit constructions of optimal measurement schemes. For n-qutrit systems, we prove that pairwise tomography requires at most 8+56⌈log8n⌉ measurement settings, and provide an explicit scheme achieving this bound. Furthermore, we develop an efficient algorithm to determine the optimal order of these measurement settings, minimizing the experimental overhead associated with switching configurations. Compared to the worst-case ordering, our optimized schedule reduces switching costs by approximately 50\%. These results provide a practical pathway for efficient characterization of qudit systems, facilitating their application in quantum communication and computation.
Quantum state tomography (QST) is a crucial tool for characterizing quantum states. However, performing full QST becomes impractical for reconstructing multiqubit density matrices since datasets and computational costs grow exponentially with qubit number. To bypass full QST, we introduce selective and efficient QST (SEEQST)—a method that enables efficient estimation of multiple selected elements of an arbitrary multiqubit density matrix. We show that any N -qubit density matrix can be partitioned into 2 N disjoint subsets, each containing 2 N elements. With SEEQST, any such subset can be accurately estimated from just two experiments with only single-qubit measurements. The complexity for estimating any subset remains constant regardless of Hilbert-space dimension, so, if desired, SEEQST can reconstruct the full density matrix, using 2 N + 1 − 1 experiments, where standard methods would use 3 N experiments. We provide a circuit decomposition for the SEEQST experiments, demonstrating that their maximum circuit depth scales logarithmically with N assuming all-to-all connectivity.
Abstract Quantum state tomography (QST) is a widely employed technique for characterizing the state of a quantum system. However, it is plagued by two fundamental challenges: computational and experimental complexity grows exponentially with the number of qubits, rendering experimental implementation and data post-processing arduous even for moderately sized systems. Here, we introduce gradient-descent (GD) algorithms for the post-processing step of QST in discrete- and continuous-variable systems. To ensure physically valid state reconstruction at each iteration step of the algorithm, we use various density-matrix parameterizations: Cholesky decomposition, Stiefel manifold, and projective normalization. These parameterizations have the added benefit of enabling a rank-controlled ansatz, which simplifies reconstruction when there is prior information about the system. We benchmark the performance of our GD-QST techniques against state-of-the-art methods, including constrained convex optimization, conditional generative adversarial networks, and iterative maximum likelihood estimation. Our comparison focuses on time complexity, iteration counts, data requirements, state rank, and robustness against noise. We find that rank-controlled ansatzes in our stochastic mini-batch GD-QST algorithms effectively handle noisy and incomplete data sets, yielding significantly higher reconstruction fidelity than other methods. Simulations achieving full-rank seven-qubit QST in under three minutes on a standard laptop, with 18 GB of RAM and no dedicated GPU, highlight that GD-QST is computationally more efficient and outperforms other techniques in most scenarios, offering a promising avenue for characterizing noisy intermediate-scale quantum devices. Our Python code for GD-QST algorithms is publicly available at github.com/mstorresh/GD-QST .
Measuring the expectation value of Pauli operators on prepared quantum states is a fundamental task in a multitude of quantum algorithms. Simultaneously measuring sets of operators allows for fewer measurements and an overall speedup of the measurement process. We investigate the task of partitioning a random subset of Pauli operators into simultaneously-measurable parts. Using heuristics from coloring random graphs, we give an upper bound for the expected number of parts in our partition. We go on to conjecture that allowing arbitrary Clifford operators before measurement, rather than single-qubit operations, leads to a decrease in the number of parts which is linear with respect to the lengths of the operators. We give evidence to confirm this conjecture and comment on the importance of this result for a specific near-term application: speeding up the measurement process of the variational quantum eigensolver.
Quantum computers solve ever more complex tasks using steadily growing system sizes.Characterizing these quantum systems is vital, yet becoming increasingly challenging.The gold-standard method for this task is quantum state tomography (QST), capable of fully reconstructing a quantum state without prior knowledge.The measurement and classical computing costs, however, increase exponentially with the number of constituents (e.g., qubits)-a daunting bottleneck given the scale of existing and nearterm quantum devices.Here, we demonstrate a scalable and practical QST approach that only uses a single measurement setting, namely symmetric informationally complete (SIC) positive operator-valued measures (POVMs).We implement these nonorthogonal measurements on an ion trap quantum processor by utilizing additional energy levels within each ion-without requiring ancillary ions to assist in measurements.More precisely, we locally map the SIC POVM to orthogonal states embedded in a higher-dimensional system, which we read out using repeated in-sequence detections, thereby providing full tomographic information in every shot.Combining this SIC tomography with the recently developed randomized measurement toolbox ("classical shadows") proves to be a powerful combination.SIC tomography alleviates the need for choosing measurement settings at random ("derandomization"), while classical shadows enable the estimation of arbitrary polynomial functions of the density matrix orders of magnitudes faster than standard methods.The latter enables in-depth entanglement characterization, which we experimentally showcase on a five-qubit absolutely maximally entangled state.Moreover, the fact that the full tomography information is available in every shot enables online QST in real time (i.e., while the experiment is running).We demonstrate this on an eight-qubit entangled state (which has 2 8 • 2 8 -1 = 65 535 degrees of freedom), as well as for fast state identification.All in all, these features single out SIC-based classical shadow estimation as a highly scalable and convenient tool for quantum state characterization.