REVIEW 4 major objections 6 minor 2 cited by
Harnessing CUDA-Q's MPS for Tensor Network Simulations of Large-Scale Quantum Circuits
T0 review · 4 major / 6 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read CUDA-Q's matrix-product-state simulator runs 60-qubit circuits on a single GPU, where state-vector simulation cannot fit any existing system.
desk verdict A solid single-GPU benchmark of CUDA-Q MPS with useful runtime data, but the 60-qubit result is runtime-only; accuracy is validated for one 10-qubit circuit. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the matrix product state (MPS), a tensor network shaped as a chain of tensors in which each qubit corresponds to one tensor and adjacent tensors share a bond of dimension $\chi$. Starting from a full state tensor, repeated singular value decompositions produce the MPS, and truncating singular values—here capped at max bond 64 with absolute and relative cutoffs of $10^{-5}$—yields a compact approximation whose memory cost is $d n \chi^2$, only linear in the number of qubits $n$ for fixed $\chi$. The same SVD truncation is what both enables large-qubit simulation and limits accuracy, since $\chi$ must be large enough to capture the circuit's entanglement. The paper uses the entanglement ratio (two-qubit gates over total gates) to characterize circuits and finds constant asymptotic ratios between 0.25 and 1.0 across the five benchmarks, which lets it connect circuit structure to MPS scalability.
What would settle it
Run the same CUDA-Q MPS implementation on a 20- to 30-qubit version of each test circuit with the default max bond 64 and compare the sampled outcome distribution to an exact state-vector simulation; if the four most likely outcomes do not match the exact reference for any of the circuits, the default approximation is already untrustworthy at small scale, so the 60-qubit results would require a larger bond dimension that may not fit on the GPU. A cleaner test is to increase the max bond from 64 to, say, 256 at 60 qubits for one circuit: if the top outcomes change materially, the reported 60-qubit feasibility is not accompanied by stable results.
Extended reading notes
Core claim
The central discovery is that MPS simulation on a single GPU is able to simulate 60-qubit circuits in all evaluated scenarios, and 90 qubits for the GHZ circuit, whereas state-vector simulation would not be able to execute any of those simulations on any currently existing system. The paper establishes this empirically on a test system with 96 GB of GPU memory: state-vector simulation supports up to 33 qubits, while MPS reaches 60 qubits for QAOA, quantum volume, QFT, and GHZ, and 90 for GHZ. Runtime scaling in the high-qubit regime is either linear in the number of qubits (QFT, GHZ) or a power law $t = \alpha n^\beta$ (QAOA, quantum volume), in contrast to the exponential cost of state-vector simulation. The paper also shows that for qubit counts below about 12, state-vector simulation is faster, and that the SVD phase dominates MPS execution (about 70% of runtime) while leaving the GPU only about a third utilized. Correctness is checked on a 10-qubit QAOA circuit, where the four most probable sampled states match the state-vector reference down to bond dimension 16.
Load-bearing premise
The load-bearing premise is that the five benchmark circuits represent real quantum workloads well enough in entanglement structure; if a real target algorithm is more entangled, the fixed bond dimension of 64 will either produce wrong answers or require a bond dimension too large for a single GPU, and the 60-qubit feasibility claim will not carry over.
Editorial extensions
If this is right
- A 96 GB single GPU can simulate 60-qubit versions of low-entanglement circuits, and 90-qubit GHZ, where no currently existing system could execute those as state-vector simulations.
- For these circuits the per-shot MPS runtime grows linearly or polynomially with qubit count, so the practical limit is set by circuit structure rather than by the exponential state size.
- State-vector simulation remains the right choice below roughly 12 to 13 qubits; tensor-network and MPS backends only win when the state vector no longer fits in memory.
- Approximation is controllable through the bond dimension: the default max bond 64 preserved the top four QAOA outcomes, the top outcomes were still intact at bond 16, and they dropped to at most two matches at bond 8 or below.
- GPU profiling shows the SVD phase dominates MPS time but under-utilizes the GPU, with about 33% average activity and less than 1% Tensor Core utilization, meaning current MPS implementations leave headroom for optimization.
Reading between the lines
- The paper's five benchmarks all have bounded entanglement ratios, so the 60-qubit result should be read as "low-entanglement circuits," not as a general license: any real workload with volume-law entanglement will require a larger bond dimension and will eventually hit a memory wall, only at a different height.
- The correctness check at 10 qubits does not by itself certify the 60-qubit outputs; a practical user should verify stability of top outcomes with respect to bond dimension on the target circuit before trusting a large MPS run.
- Because each shot in MPS costs a fresh contraction, the method is most attractive when the number of samples is modest or when approximate samples can drive an outer optimization loop, as in variational algorithms; the paper notes but does not test this use.
- If the GPU under-utilization seen in profiling is fixed, the constants in the runtime fits would improve substantially without changing the scaling laws, with the largest benefit going to the SVD-dominated MPS path.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper evaluates CUDA-Q's matrix product state (MPS) and exact tensor network simulators on an NVIDIA Grace Hopper system, comparing them with the state vector backend for five circuits (GHZ, QFT, Quantum Volume, QAOA, and counterfeit coin). The authors measure runtime scaling from low qubit counts up to 60 qubits for four circuits and 90 qubits for GHZ, report that MPS runtime scales linearly or polynomially while state vector runtime scales exponentially, and profile MPS and tensor network executions to show that SVD iterations dominate and under-utilize the GPU. They also propose a correctness validation approach based on comparing the most likely sampled outcomes against state vector results for a single 10-qubit QAOA circuit with varying bond dimension limits.
Significance. The runtime measurements are credible and the public benchmark repository (Ref. [13]) is a valuable resource. The paper's main contribution, demonstrating that a single GPU can execute MPS simulations of 60-qubit circuits for four common circuit families, is directly supported by the wall-clock data. The honest discussion of limitations (the unverified 12-to-13 qubit gap, the application-dependence of the number of relevant outcomes, and the inability to identify the role of CPU-side tasks) is commendable. However, the significance of the headline claim is weakened by the absence of any accuracy validation at large n: the only correctness check is for one 10-qubit circuit, and the entanglement-ratio metric is never quantitatively linked to MPS truncation error. As a result, the paper establishes runtime feasibility but not that the large-scale simulations are accurate.
major comments (4)
- [Section IV-B, Fig. 5] The claim that "MPS simulation is able to simulate 60-qubit circuits in all evaluated scenarios on a single GPU" is supported only by wall-clock runtime measurements; no accuracy check is reported for n ≥ 35. The only correctness validation (Section V) is a single 10-qubit QAOA instance with random parameters. Because MPS is approximate, a simulation that completes quickly may return an arbitrarily inaccurate distribution. The paper should either provide fidelity or truncation-error estimates at large n, or explicitly restrict the claim to runtime feasibility.
- [Section III-A and Section V] The default MPS parameters (CUDAQ_MPS_MAX_BOND=64, absolute and relative cutoffs 1e-5) are used throughout, and the correctness experiment shows that for one 10-qubit QAOA instance the top-4 outcomes are preserved down to χmax=16. This does not establish that χmax=64 is sufficient for QFT, QV, or QAOA at 60 qubits, where the required bond dimension depends on the entanglement structure and can grow with depth. The paper should report convergence checks (e.g., varying χmax at intermediate sizes or computing singular value spectra) to bound the approximation error.
- [Section V-B, Table III] The correctness assessment is restricted to the four most probable states from one instance, with no quantitative measure such as fidelity, total variation distance, or statistical confidence. The paper itself notes that the number of relevant outcomes is application-dependent. This is too narrow a basis for generalizing correctness to the large-scale simulations. A more robust validation (e.g., multiple circuits, varying sizes, and a proper distance metric) is needed.
- [Section IV-B, Fig. 5] The scaling "models" are empirical fits (t = αn^β, t = an+b) with parameters estimated from the measured data; calling this "analytically model[ing] the scalability" in the contribution list overstates the contribution. The fits have no uncertainty quantification, and the extrapolation of state vector runtime beyond the measured range relies on the exponential fit. The memory argument alone suffices to show state vector infeasibility at 60 qubits; the fitting language should be corrected.
minor comments (6)
- [Section III-B, Fig. 3] The legend spells "Couterfeit Coin" instead of "Counterfeit Coin"; the typo should be corrected.
- [Section IV-A] The performance gap from 12 to 13 qubits is speculatively attributed to "cache" or "a change of behavior of the simulator," and the authors state they could not verify it in documentation; this should be presented as an unresolved observation rather than a conclusion.
- [Section III-B] The statement that "This property would allow us to correlate this metric with tensor network simulation performance" is not followed by any actual correlation analysis; either add such an analysis or remove the implication.
- [Section IV-C2] The suggestion that CPU-side tasks are "ill-suited for GPU parallelization" is not directly evidenced by the profiling; it should be clearly presented as a hypothesis.
- [Section V-B, Table III] The table caption mentions "Green color" but the table may render in grayscale; ensure that the color coding is visible in print.
- [Section II-B] The expression "the number of parameters to represent a tensor network is dnχ²" is stated without specifying that it applies to an MPS; clarify that n is the number of qubits and the expression is for a fixed bond dimension.
Circularity Check
No significant circularity: the runtime and correctness results are measured against independent baselines, not derived from assumptions containing the conclusions.
full rationale
The paper's central claims are empirical. The 60-qubit MPS feasibility statement is a direct runtime observation (Section IV-B, Fig. 5), and the scaling laws t = α n^β and t = a n + b are fits to measured data, not predictions derived from the inputs; the state-vector 'projected runtime' is an extrapolation of observed exponential cost and is not used to define the MPS result. Correctness of MPS is validated against an independent exact state-vector reference on a 10-qubit QAOA instance (Section V-A), and the truncation-parameter study in Section V-B directly measures the effect of χ_max against that reference. Self-citations (refs [3], [8], [11], [15]) appear only in background, related-work, and a hardware-capacity statement; none carries the load of the MPS feasibility or accuracy claims. The paper itself limits correctness validation to a single small circuit, stating 'In this section, we focus on a single 10-qubit QAOA circuit' (Section V), which is an acknowledged validation gap for the high-qubit runs, not a circularity, since those runs are reported as runtime and approximate simulations rather than as results derived from the 10-qubit check. No equation or fitted parameter is defined in terms of the claim it supports, so no circular step can be exhibited.
Assumptions & free parameters
free parameters (5)
- QAOA power-law exponent and coefficient =
alpha = 0.0017, beta = 3.1644
- Quantum Volume power-law exponent and coefficient =
alpha = 0.0008, beta = 3.5786
- QFT linear-fit slope and intercept =
a = 5.53, b = -155.16
- GHZ linear-fit slope and intercept =
a = 2.11, b = -50.82
- State-vector exponential fit constants =
2E-10 exp(0.742 n) for QAOA/QV, 2E-10 exp(0.6867 n) for QFT
assumptions (4)
- standard math MPS memory requirement for n qubits with bond dimension chi is 2 n chi^2 complex values (Section II-B).
- standard math State vector simulation requires 2^n complex values (Section II-B).
- domain assumption Entanglement ratio (fraction of two-qubit gates) is a proxy for computational cost of tensor network simulation (Section III-B).
- ad hoc to paper Default MPS parameters (max bond 64, cutoff 1e-5) are sufficient to produce correct results for the evaluated circuits (Section III-A, V).
Cite this review
Pith. "Pith review of Harnessing CUDA-Q's MPS for Tensor Network Simulations of Large-Scale Quantum Circuits." pith.science (2026). https://pith.science/paper/B62V7YIR
@misc{pith2026250115939,
author = {Pith},
title = {Pith review of: Harnessing CUDA-Q's MPS for Tensor Network Simulations of Large-Scale Quantum Circuits},
year = {2026},
howpublished = {\url{https://pith.science/paper/B62V7YIR}},
note = {Machine review of arXiv:2501.15939}
}
read the original abstract
Quantum computer simulators are an indispensable tool for prototyping quantum algorithms and verifying the functioning of existing quantum computer hardware. The current largest quantum computers feature more than one thousand qubits, challenging their classical simulators. State-vector quantum simulators are challenged by the exponential increase of representable quantum states with respect to the number of qubits, making more than fifty qubits practically unfeasible. A more appealing approach for simulating quantum computers is adopting the tensor network approach, whose memory requirements fundamentally depend on the level of entanglement in the quantum circuit, and allows simulating the current largest quantum computers. This work investigates and evaluates the CUDA-Q tensor network simulators on an Nvidia Grace Hopper system, particularly the Matrix Product State (MPS) formulation. We compare the performance of the CUDA-Q state vector implementation and validate the correctness of MPS simulations. Our results highlight that tensor network-based methods provide a significant opportunity to simulate large-qubit circuits, albeit approximately. We also show that current GPU-accelerated computation cannot fully utilize GPU efficiently in the case of MPS simulations.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 2 Pith papers
-
Benchmarking Zero-Setup Quantum Circuit Simulators
GPU-accelerated zero-setup simulators show sub-quadratic MPS bond-dimension scaling and up to 1,400× PPS speedups, uniquely reaching fine truncation accuracy on the 127-qubit kicked Ising circuit.
-
Intelligence-Guided Adaptive Purification for DDoS-Resilient Quantum Networks: A CUDA-Q based Study
IDS-driven adaptive purification in a simulated 8-node quantum repeater chain restores fidelity-qualified entanglement delivery under SSDP-induced degradation (0.098 to 0.344 above-target; oracle 0.335).
Reference graph
Works this paper leans on
-
[13]
https://github.com/KTH-ScaLab/cudaq-benchmarks/
-
[1]
Simulating the Sycamore quantum supremacy circuits,
F. Pan and P. Zhang, “Simulating the Sycamore quantum supremacy circuits,” Mar. 2021, arXiv:2103.03074
arXiv 2021
-
[2]
mpiQulacs: A Scalable Distributed Quantum Computer Simulator for ARM-based Clusters,
A. Tabuchi, S. Imamura, M. Yamazaki, T. Honda, A. Kasagi, H. Nakao, N. Fukumoto, and K. Nakashima, “mpiQulacs: A Scalable Distributed Quantum Computer Simulator for ARM-based Clusters,” in 2023 IEEE In- ternational Conference on Quantum Computing and En- gineering (QCE), vol. 01, Sep. 2023, pp. 959–969
work page 2023
-
[3]
What is quantum parallelism, anyhow?
S. Markidis, “What is quantum parallelism, anyhow?” in ISC High Performance 2024 Research Paper Proceedings (39th International Conference) . Prometeus GmbH, 2024, pp. 1–12
work page 2024
-
[4]
cuTN-QSVM: cuTensorNet-accelerated Quantum Sup- port Vector Machine with cuQuantum SDK,
K.-C. Chen, T.-Y . Li, Y .-Y . Wang, S. See, C.-C. Wang, R. Wille, N.-Y . Chen, A.-C. Yang, and C.-Y . Lin, “cuTN-QSVM: cuTensorNet-accelerated Quantum Sup- port Vector Machine with cuQuantum SDK,” May 2024, arXiv:2405.02630 [quant-ph]
arXiv 2024
-
[5]
Tensor Network Quantum Virtual Machine for Simulating Quantum Circuits at Exascale,
T. Nguyen, D. Lyakh, E. Dumitrescu, D. Clark, J. Larkin, and A. McCaskey, “Tensor Network Quantum Virtual Machine for Simulating Quantum Circuits at Exascale,” ACM Transactions on Quantum Computing, vol. 4, no. 1, pp. 6:1–6:21, Oct. 2022
work page 2022
- [6]
-
[7]
Google, “Cirq,” https://quantumai.google/cirq, 2024
work page 2024
Show all 32 references
-
[8]
Quan- tum computer simulations at warp speed: Assessing the impact of gpu acceleration: A case study with ibm qiskit aer, nvidia thrust & cuquantum,
J. Faj, I. Peng, J. Wahlgren, and S. Markidis, “Quan- tum computer simulations at warp speed: Assessing the impact of gpu acceleration: A case study with ibm qiskit aer, nvidia thrust & cuquantum,” in 2023 IEEE 19th International Conference on e-Science (e-Science) . IEEE, 202...
2023
-
[9]
Azure quantum,
Microsoft, “Azure quantum,” https://quantum.microsoft. com/, 2024
2024
-
[10]
Q. A. team and collaborators, “qsim,” https://doi.org/10. 5281/zenodo.4023103, 2020
2020
-
[11]
Enabling quantum computer simulations on amd gpus: a hip backend for google’s qsim,
S. Markidis, “Enabling quantum computer simulations on amd gpus: a hip backend for google’s qsim,” in Proceedings of the SC’23 Workshops of The International Conference on High Performance Computing, Network, Storage, and Analysis , 2023, pp. 1478–1486
2023
-
[12]
Cuda-q documentation,
Nvidia, “Cuda-q documentation,” https://nvidia.github.io/ cuda-quantum/latest/index.html
-
[14]
Tensor Networks in a Nutshell,
J. Biamonte and V . Bergholm, “Tensor Networks in a Nutshell,” Jul. 2017, arXiv:1708.00006
2017 arXiv
-
[15]
Harnessing integrated cpu-gpu system memory for hpc: a first look into grace hopper,
G. Schieffer, J. Wahlgren, J. Ren, J. Faj, and I. Peng, “Harnessing integrated cpu-gpu system memory for hpc: a first look into grace hopper,” in Proceedings of the 53rd International Conference on Parallel Processing , 2024, pp. 199–209
2024
-
[16]
SV-sim: scalable PGAS-based state vector simulation of quantum circuits,
A. Li and et al., “SV-sim: scalable PGAS-based state vector simulation of quantum circuits,” in Proceedings of the SC21, ser. SC ’21. New York, NY , USA: Association for Computing Machinery, Nov. 2021, pp. 1–14
2021
-
[17]
cuquantum sdk documentation,
Nvidia, “cuquantum sdk documentation,” https://docs. nvidia.com/cuda/cuquantum/
-
[18]
Quantum++: A modern c++ quantum computing library,
V . Gheorghiu, “Quantum++: A modern c++ quantum computing library,” PLOS ONE , vol. 13, no. 12, pp. 1– 27, 12 2018
2018
-
[19]
Exatn: Scalable gpu-accelerated high- performance processing of general tensor networks at ex- ascale,
D. I. Lyakh, T. Nguyen, D. Claudino, E. Dumitrescu, and A. J. McCaskey, “Exatn: Scalable gpu-accelerated high- performance processing of general tensor networks at ex- ascale,” Frontiers in Applied Mathematics and Statistics , vol. 8, p. 838601, 2022
2022
-
[20]
SupermarQ: A Scalable Quantum Benchmark Suite,
T. Tomesh, P. Gokhale, V . Omole, G. S. Ravi, K. N. Smith, J. Viszlai, X.-C. Wu, N. Hardavellas, M. R. Martonosi, and F. T. Chong, “SupermarQ: A Scalable Quantum Benchmark Suite,” in Proceedings of the 2022 IEEE International Symposium on High-Performance Computer Architecture...
2022
-
[21]
QASM- Bench: A Low-Level Quantum Benchmark Suite for NISQ Evaluation and Simulation,
A. Li, S. Stein, S. Krishnamoorthy, and J. Ang, “QASM- Bench: A Low-Level Quantum Benchmark Suite for NISQ Evaluation and Simulation,” ACM Transactions on Quantum Computing, vol. 4, no. 2, pp. 10:1–10:26, 2023
2023
-
[22]
Cache Blocking Technique to Large Scale Quantum Computing Simulation on Supercomput- ers,
J. Doi and H. Horii, “Cache Blocking Technique to Large Scale Quantum Computing Simulation on Supercomput- ers,” in 2020 IEEE International Conference on Quantum Computing and Engineering (QCE), Oct. 2020, pp. 212– 222
2020
-
[23]
MEMQSim: Highly Memory-Efficient and Modularized Quantum State-Vector Simulation,
B. Zhang, B. Fang, Q. Guan, A. Li, and D. Tao, “MEMQSim: Highly Memory-Efficient and Modularized Quantum State-Vector Simulation,” in Proceedings of the SC ’23 Workshops of The International Conference on High Performance Computing, Network, Storage, and Analysis. ACM, Nov. 20...
2023
-
[24]
Quantum computing simulator on a heteroge- nous HPC system,
J. Doi, H. Takahashi, R. Raymond, T. Imamichi, and H. Horii, “Quantum computing simulator on a heteroge- nous HPC system,” in Proceedings of the 16th ACM In- ternational Conference on Computing Frontiers , ser. CF ’19. New York, NY , USA: Association for Computing Machinery, A...
2019
-
[25]
Prototype of a Batched Quantum Circuit Simulator for the Vec- tor Engine,
K. Takahashi, T. Mori, and H. Takizawa, “Prototype of a Batched Quantum Circuit Simulator for the Vec- tor Engine,” in Proceedings of the SC ’23 Workshops of The International Conference on High Performance Computing, Network, Storage, and Analysis , ser. SC-W ’23. New York, N...
2023
-
[26]
DiaQ: Efficient State-Vector Quantum Simulation,
S. Chundury, J. Li, I.-S. Suh, and F. Mueller, “DiaQ: Efficient State-Vector Quantum Simulation,” Apr. 2024, arXiv:2405.01250 [quant-ph]
2024 arXiv
-
[27]
Quantum Computer Simulations at Warp Speed: Assessing the Impact of GPU Acceleration: A Case Study with IBM Qiskit Aer, Nvidia Thrust & cuQuantum,
J. Faj, I. Peng, J. Wahlgren, and S. Markidis, “Quantum Computer Simulations at Warp Speed: Assessing the Impact of GPU Acceleration: A Case Study with IBM Qiskit Aer, Nvidia Thrust & cuQuantum,” in 2023 IEEE 19th International Conference on e-Science (e-Science) , Oct. 2023, pp. 1–10
2023
-
[28]
Validating quantum-classical programming models with tensor network simulations,
A. McCaskey, E. Dumitrescu, M. Chen, D. Lyakh, and T. Humble, “Validating quantum-classical programming models with tensor network simulations,” PLOS ONE , vol. 13, no. 12, Dec. 2018, publisher: Public Library of Science
2018
-
[29]
Tensor-network-assisted variational quantum algorithm,
J. Huang, W. He, Y . Zhang, Y . Wu, B. Wu, and X. Yuan, “Tensor-network-assisted variational quantum algorithm,” Physical Review A, vol. 108, no. 5, p. 052407, Nov. 2023, publisher: American Physical Society
2023
-
[30]
Efficient Classical Simulation of Slightly En- tangled Quantum Computations,
G. Vidal, “Efficient Classical Simulation of Slightly En- tangled Quantum Computations,” Physical Review Let- ters, vol. 91, no. 14, p. 147902, Oct. 2003, publisher: American Physical Society
2003
-
[31]
Tensor networks for complex quantum sys- tems,
R. Or ´us, “Tensor networks for complex quantum sys- tems,” Nature Reviews Physics , vol. 1, no. 9, pp. 538– 550, Sep. 2019, publisher: Nature Publishing Group
2019
-
[32]
State of practice: evaluating GPU performance of state vector and tensor network methods,
M. Vallero, F. Vella, and P. Rech, “State of practice: evaluating GPU performance of state vector and tensor network methods,” Jan. 2024, arXiv:2401.06188
2024 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.