Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2401.10330 · Quantum

Optimal sampling of tensor networks targeting wave function's fast decaying tails

Marco Ballarin, Pietro Silvi, Simone Montangero, Daniel Jaschke

Published Apr 18, 2025• 1 citation
arXivHTMLPDF

Abstract

We introduce an optimal strategy to sample quantum outcomes of local measurement strings for isometric tensor network states. Our method generates samples based on an exact cumulative bounding function, without prior knowledge, in the minimal amount of tensor network contractions. The algorithm avoids sample repetition and, thus, is efficient at sampling distribution with exponentially decaying tails. We illustrate the computational advantage provided by our optimal sampling method through various numerical examples, involving condensed matter, optimization problems, and quantum circuit scenarios. Theory predicts up to an exponential speedup reducing the scaling for sampling the space up to an accumulated unknown probability ϵ from O(ϵ−1) to O(log⁡(ϵ−1)) for a decaying probability distribution. We confirm this in practice with over one order of magnitude speedup or multiple orders improvement in the error depending on the application. Our sampling strategy extends beyond local observables, e.g., to quantum magic.

DOI
10.22331/q-2025-04-18-1714

Related papers

Sampling two-dimensional isometric tensor network states

Alec Dektor, Eugene Dumitrescu, Chao Yang

Sampling a quantum system's underlying probability distributions is an important computational task, e.g., for quantum advantage experiments and quantum Monte Carlo algorithms. Tensor networks are an invaluable tool for efficiently representing states of large quantum systems with limited entanglement. Algorithms for sampling one-dimensional (1D) tensor networks are well-established and utilized in several 1D tensor network methods. In this paper we introduce two novel sampling algorithms for two-dimensional (2D) isometric tensor network states (isoTNS) that generalize existing 1D tensor networks sampling algorithms. Our first proposed algorithm performs independent sampling and yields a single configuration together with its associated probability. The second algorithm employs a greedy search strategy to identify $K$ high-probability configurations and their corresponding probabilities. Numerical results demonstrate the effectiveness of these algorithms across quantum states with varying entanglement and system size.

quant-pharXiv:2602.02245•Journal of Computational Science•Feb 3, 2026•0 citations

Perfect sampling with unitary tensor networks

Andrew J. Ferris, Guifré Vidal

Tensor network states are powerful variational Ans\"atze for many-body ground states of quantum lattice models. The use of Monte Carlo sampling techniques in tensor network approaches significantly reduces the cost of tensor contractions, potentially leading to a substantial increase in computational efficiency. Previous proposals are based on a Markov chain Monte Carlo scheme generated by locally updating configurations and, as such, must deal with equilibration and autocorrelation times, which result in a reduction of efficiency. Here we propose perfect sampling schemes, with vanishing equilibration and autocorrelation times, for unitary tensor networks, namely, tensor networks based on efficiently contractible, unitary quantum circuits, such as unitary versions of the matrix product state (MPS) and tree tensor network (TTN), and the multiscale entanglement renormalization Ansatz (MERA). Configurations are directly sampled according to their probabilities in the wave function, without resorting to a Markov chain process. We consider both complete sampling, involving all the relevant sites of the system, and incomplete sampling, which only involves a subset of those sites and which can result in a dramatic (basis-dependent) reduction of sampling error.

arXiv:1201.3974•Physical Review B•Apr 30, 2012•138 citations

Quantum-Inspired Optimization of Deep Learning Hyperparameters via Tensor Networks

Jincheng Zhang

Optimizing hyperparameters in deep learning models remains a significant bottleneck in the development and deployment of effective neural networks. Traditional methods, such as grid search and random search, often struggle with the exponentially growing search space, leading to substantial computational costs and limited exploration of potentially optimal configurations. This paper proposes a novel approach to deep learning hyperparameter optimization that leverages the principles of quantum-inspired optimization, specifically employing tensor network methods, notably Matrix Product States (MPS), to efficiently approximate the complex optimization landscape. Our method utilizes quantum-inspired sampling techniques to guide the exploration, effectively mimicking the probabilistic nature of quantum measurements. We demonstrate the feasibility and potential advantages of this approach through theoretical analysis and a conceptual framework. The core of our method lies in representing the hyperparameter space as a tensor network, allowing for efficient calculations of the objective function (e.g., validation loss) without exhaustively evaluating every possible combination. This significantly reduces the computational burden compared to conventional methods. The paper outlines a detailed algorithm and discusses key considerations for implementation, including tensor network construction, sampling strategies, and convergence criteria. We propose a framework for iteratively refining the tensor network representation based on observed validation performance, ultimately converging towards a high-performing hyperparameter configuration. The theoretical underpinning draws upon concepts from quantum mechanics, particularly the idea of measurement and probability distribution, to drive the optimization process. We believe this quantum-inspired approach offers a promising avenue for accelerating deep learning hyperparameter tuning and achieving state-of-the-art performance.

•Sep 2, 2026•0 citations

Accelerating Quantum Tensor Network Simulations with Unified Path Variations and Non-Degenerate Batched Sampling

Taylor Lee Patti, Paavai Pari, Yang Gao +5 more

Quantum trajectory methods reduce the computational overhead of simulating noisy quantum systems, approximating them with $m$ stochastically sampled $2^n$-entry quantum statevectors rather than exact $2^{2n}$-entry density matrices. Recently, Pre-Trajectory Sampling with Batched Execution (PTSBE) has dramatically increased the data collection rate of these methods. While statevector PTSBE has demonstrated data collection speedups of over $10^6 \times$, tensor network implementations only achieved $\sim 15 \times$ speedup. This comparatively modest tensor network advantage stemmed from 1) contraction path recalculations, 2) sequential tensor network sampling, and 3) inflexible/unoptimized contraction hyperparameters. In this manuscript, we increase PTSBE's tensor network data collection rate to more than $10^8\times$ that of traditional trajectories methods by developing 1) error-independent unified path variation, 2) non-degenerate tensor network sampling, and 3) a flexible/optimized contraction framework. While our methods are particularly powerful for accelerating non-proportional sampling, we also demonstrate a more than $1000\times$ speedup for more general quantum simulations.

quant-pharXiv:2604.08467•Apr 10, 2026•0 citations

Solving the Sampling Problem of the Sycamore Quantum Circuits

Feng Ping Pan, Keyang Chen, Pan Zhang

We study the problem of generating independent samples from the output distribution of Google's Sycamore quantum circuits with a target fidelity, which is believed to be beyond the reach of classical supercomputers and has been used to demonstrate quantum supremacy. We propose a method to classically solve this problem by contracting the corresponding tensor network just once, and is massively more efficient than existing methods in generating a large number of uncorrelated samples with a target fidelity. For the Sycamore quantum supremacy circuit with 53 qubits and 20 cycles, we have generated 1×10^{6} uncorrelated bitstrings s which are sampled from a distribution P[over ^](s)=|ψ[over ^](s)|^{2}, where the approximate state ψ[over ^] has fidelity F≈0.0037. The whole computation has cost about 15 h on a computational cluster with 512 GPUs. The obtained 1×10^{6} samples, the contraction code and contraction order are made public. If our algorithm could be implemented with high efficiency on a modern supercomputer with ExaFLOPS performance, we estimate that ideally, the simulation would cost a few dozens of seconds, which is faster than Google's quantum hardware.

arXiv:2111.03011•Physical Review Letters•Aug 22, 2022•142 citations

Relational Quantum Circuit Synthesis via Tensor Network Optimization

Jincheng Zhang

The synthesis of efficient quantum circuits for complex computations remains a significant challenge in quantum computing. Traditional methods often struggle with the exponential growth of the solution space, particularly as the number of qubits increases. This work introduces a novel approach to quantum circuit synthesis leveraging tensor network optimization. We represent quantum circuits as tensor networks and employ techniques such as Tensor Train and Product Tensor Decomposition to efficiently approximate the optimal circuit topology. A quantum-inspired probabilistic sampling strategy is then integrated to explore the solution space effectively. The core claim is that optimizing quantum circuit synthesis for complex systems is intrinsically intractable using conventional methods. This approach offers a potentially scalable alternative, providing a pathway toward designing practical quantum circuits for a wider range of applications. The method is presented with a focus on the mathematical formalism and optimization techniques involved.

•Aug 28, 2026•0 citations
We introduce an optimal strategy to sample quantum outcomes of local measurement strings for isometric tensor network states. Our method generates samples based on an exact cumulative bounding function, without prior knowledge, in the minimal amount of tensor network contractions. The algorithm avoids sample repetition and, thus, is efficient at sampling distribution with exponentially decaying tails. We illustrate the computational advantage provided by our optimal sampling method through various numerical examples, involving condensed matter, optimization problems, and quantum circuit scenarios. Theory predicts up to an exponential speedup reducing the scaling for sampling the space up to an accumulated unknown probability ϵ from O(ϵ−1) to O(log⁡(ϵ−1)) for a decaying probability distribution. We confirm this in practice with over one order of magnitude speedup or multiple orders improvement in the error depending on the application. Our sampling strategy extends beyond local observables, e.g., to quantum magic.
Sampling a quantum system's underlying probability distributions is an important computational task, e.g., for quantum advantage experiments and quantum Monte Carlo algorithms. Tensor networks are an invaluable tool for efficiently representing states of large quantum systems with limited entanglement. Algorithms for sampling one-dimensional (1D) tensor networks are well-established and utilized in several 1D tensor network methods. In this paper we introduce two novel sampling algorithms for two-dimensional (2D) isometric tensor network states (isoTNS) that generalize existing 1D tensor networks sampling algorithms. Our first proposed algorithm performs independent sampling and yields a single configuration together with its associated probability. The second algorithm employs a greedy search strategy to identify KKK high-probability configurations and their corresponding probabilities. Numerical results demonstrate the effectiveness of these algorithms across quantum states with varying entanglement and system size.
Tensor network states are powerful variational Ans\"atze for many-body ground states of quantum lattice models. The use of Monte Carlo sampling techniques in tensor network approaches significantly reduces the cost of tensor contractions, potentially leading to a substantial increase in computational efficiency. Previous proposals are based on a Markov chain Monte Carlo scheme generated by locally updating configurations and, as such, must deal with equilibration and autocorrelation times, which result in a reduction of efficiency. Here we propose perfect sampling schemes, with vanishing equilibration and autocorrelation times, for unitary tensor networks, namely, tensor networks based on efficiently contractible, unitary quantum circuits, such as unitary versions of the matrix product state (MPS) and tree tensor network (TTN), and the multiscale entanglement renormalization Ansatz (MERA). Configurations are directly sampled according to their probabilities in the wave function, without resorting to a Markov chain process. We consider both complete sampling, involving all the relevant sites of the system, and incomplete sampling, which only involves a subset of those sites and which can result in a dramatic (basis-dependent) reduction of sampling error.
Optimizing hyperparameters in deep learning models remains a significant bottleneck in the development and deployment of effective neural networks. Traditional methods, such as grid search and random search, often struggle with the exponentially growing search space, leading to substantial computational costs and limited exploration of potentially optimal configurations. This paper proposes a novel approach to deep learning hyperparameter optimization that leverages the principles of quantum-inspired optimization, specifically employing tensor network methods, notably Matrix Product States (MPS), to efficiently approximate the complex optimization landscape. Our method utilizes quantum-inspired sampling techniques to guide the exploration, effectively mimicking the probabilistic nature of quantum measurements. We demonstrate the feasibility and potential advantages of this approach through theoretical analysis and a conceptual framework. The core of our method lies in representing the hyperparameter space as a tensor network, allowing for efficient calculations of the objective function (e.g., validation loss) without exhaustively evaluating every possible combination. This significantly reduces the computational burden compared to conventional methods. The paper outlines a detailed algorithm and discusses key considerations for implementation, including tensor network construction, sampling strategies, and convergence criteria. We propose a framework for iteratively refining the tensor network representation based on observed validation performance, ultimately converging towards a high-performing hyperparameter configuration. The theoretical underpinning draws upon concepts from quantum mechanics, particularly the idea of measurement and probability distribution, to drive the optimization process. We believe this quantum-inspired approach offers a promising avenue for accelerating deep learning hyperparameter tuning and achieving state-of-the-art performance.
Quantum trajectory methods reduce the computational overhead of simulating noisy quantum systems, approximating them with mmm stochastically sampled 2n2^n2n-entry quantum statevectors rather than exact 22n2^{2n}22n-entry density matrices. Recently, Pre-Trajectory Sampling with Batched Execution (PTSBE) has dramatically increased the data collection rate of these methods. While statevector PTSBE has demonstrated data collection speedups of over 106×10^6 \times106×, tensor network implementations only achieved ∼15×\sim 15 \times∼15× speedup. This comparatively modest tensor network advantage stemmed from 1) contraction path recalculations, 2) sequential tensor network sampling, and 3) inflexible/unoptimized contraction hyperparameters. In this manuscript, we increase PTSBE's tensor network data collection rate to more than 108×10^8\times108× that of traditional trajectories methods by developing 1) error-independent unified path variation, 2) non-degenerate tensor network sampling, and 3) a flexible/optimized contraction framework. While our methods are particularly powerful for accelerating non-proportional sampling, we also demonstrate a more than 1000×1000\times1000× speedup for more general quantum simulations.
We study the problem of generating independent samples from the output distribution of Google's Sycamore quantum circuits with a target fidelity, which is believed to be beyond the reach of classical supercomputers and has been used to demonstrate quantum supremacy. We propose a method to classically solve this problem by contracting the corresponding tensor network just once, and is massively more efficient than existing methods in generating a large number of uncorrelated samples with a target fidelity. For the Sycamore quantum supremacy circuit with 53 qubits and 20 cycles, we have generated 1×10^{6} uncorrelated bitstrings s which are sampled from a distribution P[over ^](s)=|ψ[over ^](s)|^{2}, where the approximate state ψ[over ^] has fidelity F≈0.0037. The whole computation has cost about 15 h on a computational cluster with 512 GPUs. The obtained 1×10^{6} samples, the contraction code and contraction order are made public. If our algorithm could be implemented with high efficiency on a modern supercomputer with ExaFLOPS performance, we estimate that ideally, the simulation would cost a few dozens of seconds, which is faster than Google's quantum hardware.
The synthesis of efficient quantum circuits for complex computations remains a significant challenge in quantum computing. Traditional methods often struggle with the exponential growth of the solution space, particularly as the number of qubits increases. This work introduces a novel approach to quantum circuit synthesis leveraging tensor network optimization. We represent quantum circuits as tensor networks and employ techniques such as Tensor Train and Product Tensor Decomposition to efficiently approximate the optimal circuit topology. A quantum-inspired probabilistic sampling strategy is then integrated to explore the solution space effectively. The core claim is that optimizing quantum circuit synthesis for complex systems is intrinsically intractable using conventional methods. This approach offers a potentially scalable alternative, providing a pathway toward designing practical quantum circuits for a wider range of applications. The method is presented with a focus on the mathematical formalism and optimization techniques involved.