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:2403.16729 · Physical Review Applied

Quantum state preparation for probability distributions with reflection symmetry using matrix-product states

Yuichi Sano, Ikko Hamamura

Published Aug 25, 2025• 3 citations
arXivHTMLPDF

Abstract

Quantum circuits for loading probability distributions into quantum states are essential subroutines in quantum algorithms used in physics, finance engineering, and machine learning. The ability to implement these with high accuracy in low-depth quantum circuits is a critical issue. We propose a novel quantum state preparation method for probability distribution with reflection symmetry using matrix product states. By considering reflection symmetry, our method reduces the entanglement of probability distributions and improves the accuracy of approximations by matrix product states. As a result, we improved the accuracy by two orders of magnitude over existing methods using matrix product states. Our approach, characterized by linear scalability with qubit count, is highly advantageous for noisy quantum devices. Also, our demonstration results reveal that the approximation accuracy in tensor networks depends heavily on the bond dimension, with minimal reliance on the number of qubits. Our method is demonstrated for a normal distribution encoded into 10 and 20 qubits on a real quantum processor.

DOI
10.1103/yqj3-hyxv

Related papers

Entanglement scaling in matrix product state representation of smooth functions and their shallow quantum circuit approximations

Vladyslav Bohun, I. V. Lukin, Mykola Luhanko +4 more

Encoding classical data in a quantum state is a key prerequisite of many quantum algorithms. Recently, matrix product state (MPS) methods emerged as the most promising approach for constructing shallow quantum circuits approximating input functions, including probability distributions, with only linear number of gates. We derive rigorous asymptotic expansions for the decay of entanglement across bonds in the MPS representation depending on the smoothness of the input function, real or complex. We also consider the dependence of the entanglement on localization properties and function support. Based on these analytical results, we construct an improved MPS-based algorithm yielding shallow and accurate encoding quantum circuits. By using tensor cross interpolation, we are able to construct utility-scale quantum circuits in a compute- and memory-efficient way. We validate our methods by loading heavy-tailed distributions, including Lévy, important in finance, but they apply to any smooth function inputs. We test the performance of the resulting quantum circuits by executing and sampling from them on IBM quantum devices, for up to 156 qubits.

arXiv:2412.05202•Physical Review Research•Mar 2, 2026•2 citations

Minimizing entanglement entropy for enhanced quantum state preparation

Anonymous, William Steadman, Ossi Niemimäki +1 more

Quantum state preparation is an important subroutine in many quantum algorithms. The goal is to encode classical information directly to the quantum state so that it is possible to leverage quantum algorithms for data processing. However, quantum state preparation of arbitrary states scales exponentially in the number of two-qubit gates, and this makes quantum state preparation a very difficult task on quantum computers, especially on near-term noisy devices. This represents a major challenge in achieving quantum advantage. We present and analyze a novel two-step state preparation method where we first minimize the entanglement entropy of the target quantum state, thus transforming the state to one that is easier to prepare. The state with reduced entanglement entropy is then represented as a matrix product state, resulting in a high accuracy preparation of the target state. Our method is suitable for NISQ devices and we give rigorous lower bounds on the accuracy of the prepared state in terms of the entanglement entropy. We benchmark our method with 2D normal distribution and Ricker wavelet states with 6--20 qubits.

arXiv:2507.22562•Physical Review A•Aug 10, 2026•0 citations

Quantum state preparation of normal distributions using matrix product states

Jason Iaconis, Sonika Johri, Elton Yechao Zhu

Abstract State preparation is a necessary component of many quantum algorithms. In this work, we combine a method for efficiently representing smooth differentiable probability distributions using matrix product states with recently discovered techniques for initializing quantum states to approximate matrix product states. Using this, we generate quantum states encoding a class of normal probability distributions in a trapped ion quantum computer for up to 20 qubits. We provide an in depth analysis of the different sources of error which contribute to the overall fidelity of this state preparation procedure. Our work provides a study in quantum hardware for scalable distribution loading, which is the basis of a wide range of algorithms that provide quantum advantage.

arXiv:2303.01562•npj Quantum Information•Jan 25, 2024•37 citations

Exact log-depth preparation of highly entangled matrix product states

Keisuke Murota, Frédéric Sauvage, Marco Ballarin +2 more

Preparing matrix product states (MPS) on a quantum device is a key subroutine in many quantum algorithms. The most competitive methods, based on the renormalisation group, prepare translationally invariant MPS of size $L$ and bond dimension $χ$, up to an error $\varepsilon$, in circuit depth $\tilde O(χ^{4}\log(L/\varepsilon))$ or $\tilde O(χ^{6}\log\log(L/\varepsilon))$. We improve multiple aspects of these methods. First, using block-encoded correction maps, whose post-selection succeeds with constant probability, we render the preparation exact without sacrificing the scaling in $L$. Second, through a generalisation of oblivious amplitude amplification to isometries, we reduce the bond-dimension dependence, improving the depth to $\tilde O(χ^{2}\log L + χ^{4})$ or $\tilde O(χ^{2}\log\log L + χ^{4})$, and even to $\tilde O(χ^{3}\log L)$ for incoherent preparations. Finally, we extend the framework to non-translationally invariant MPS and prove logarithmic-depth exact preparation for independent and identically distributed random tensor sequences. Confirmed by numerical studies, these results constitute, to the best of our knowledge, the most efficient exact MPS preparation protocols in the relevant parameter regimes.

quant-pharXiv:2606.24475•Jun 24, 2026•0 citations

The State Preparation of Multivariate Normal Distributions using Tree Tensor Network

Hidetaka Manabe, Yuichi Sano

The quantum state preparation of probability distributions is an important subroutine for many quantum algorithms. When embedding D -dimensional multivariate probability distributions by discretizing each dimension into 2n points, we need a state preparation circuit comprising a total of nD qubits, which is often difficult to compile. In this study, we propose a scalable method to generate state preparation circuits for D -dimensional multivariate normal distributions, utilizing tree tensor networks (TTN). We establish theoretical guarantees that multivariate normal distributions with 1D correlation structures can be efficiently represented using TTN. Based on these analyses, we propose a compilation method that uses automatic structural optimization to find the most efficient network structure and compact circuit. We apply our method to state preparation circuits for various high-dimensional random multivariate normal distributions. The numerical results suggest that our method can dramatically reduce the circuit depth and CNOT count while maintaining fidelity compared to existing approaches.

arXiv:2412.12067•Quantum•May 28, 2025•6 citations

Efficient Quantum Circuits for Accurate State Preparation of Smooth, Differentiable Functions

Adam A. Holmes, A. Y. Matsuura

Effective quantum computation relies upon making good use of the exponential information capacity of a quantum machine. A large barrier to designing quantum algorithms for execution on real quantum machines is that, in general, it is intractably difficult to construct an arbitrary quantum state to high precision. Many quantum algorithms rely instead upon initializing the machine in a simple state, and evolving the state through an efficient (i.e. at most polynomial-depth) quantum algorithm. In this work, we show that there exist families of quantum states that can be prepared to high precision with circuits of linear size and depth. We focus on real-valued, smooth, differentiable functions with bounded derivatives on a domain of interest, exemplified by commonly used probability distributions. We further develop an algorithm that requires only linear classical computation time to generate accurate linear- depth circuits to prepare these states, and apply this to well-known and heavily-utilized functions including Gaussian and lognormal distributions. Our procedure rests upon the quantum state representation tool known as the matrix product state (MPS). By efficiently and scalably encoding an explicit amplitude function into an MPS, a high fidelity, linear-depth circuit can directly be generated. These results enable the execution of many quantum algorithms that, aside from initialization, are otherwise depth-efficient.

arXiv:2005.04351•IEEE International Conference on Quantum Computing and Engineering (QCE)•Oct 1, 2020•53 citations
Quantum circuits for loading probability distributions into quantum states are essential subroutines in quantum algorithms used in physics, finance engineering, and machine learning. The ability to implement these with high accuracy in low-depth quantum circuits is a critical issue. We propose a novel quantum state preparation method for probability distribution with reflection symmetry using matrix product states. By considering reflection symmetry, our method reduces the entanglement of probability distributions and improves the accuracy of approximations by matrix product states. As a result, we improved the accuracy by two orders of magnitude over existing methods using matrix product states. Our approach, characterized by linear scalability with qubit count, is highly advantageous for noisy quantum devices. Also, our demonstration results reveal that the approximation accuracy in tensor networks depends heavily on the bond dimension, with minimal reliance on the number of qubits. Our method is demonstrated for a normal distribution encoded into 10 and 20 qubits on a real quantum processor.
Encoding classical data in a quantum state is a key prerequisite of many quantum algorithms. Recently, matrix product state (MPS) methods emerged as the most promising approach for constructing shallow quantum circuits approximating input functions, including probability distributions, with only linear number of gates. We derive rigorous asymptotic expansions for the decay of entanglement across bonds in the MPS representation depending on the smoothness of the input function, real or complex. We also consider the dependence of the entanglement on localization properties and function support. Based on these analytical results, we construct an improved MPS-based algorithm yielding shallow and accurate encoding quantum circuits. By using tensor cross interpolation, we are able to construct utility-scale quantum circuits in a compute- and memory-efficient way. We validate our methods by loading heavy-tailed distributions, including Lévy, important in finance, but they apply to any smooth function inputs. We test the performance of the resulting quantum circuits by executing and sampling from them on IBM quantum devices, for up to 156 qubits.
Quantum state preparation is an important subroutine in many quantum algorithms. The goal is to encode classical information directly to the quantum state so that it is possible to leverage quantum algorithms for data processing. However, quantum state preparation of arbitrary states scales exponentially in the number of two-qubit gates, and this makes quantum state preparation a very difficult task on quantum computers, especially on near-term noisy devices. This represents a major challenge in achieving quantum advantage. We present and analyze a novel two-step state preparation method where we first minimize the entanglement entropy of the target quantum state, thus transforming the state to one that is easier to prepare. The state with reduced entanglement entropy is then represented as a matrix product state, resulting in a high accuracy preparation of the target state. Our method is suitable for NISQ devices and we give rigorous lower bounds on the accuracy of the prepared state in terms of the entanglement entropy. We benchmark our method with 2D normal distribution and Ricker wavelet states with 6--20 qubits.
Abstract State preparation is a necessary component of many quantum algorithms. In this work, we combine a method for efficiently representing smooth differentiable probability distributions using matrix product states with recently discovered techniques for initializing quantum states to approximate matrix product states. Using this, we generate quantum states encoding a class of normal probability distributions in a trapped ion quantum computer for up to 20 qubits. We provide an in depth analysis of the different sources of error which contribute to the overall fidelity of this state preparation procedure. Our work provides a study in quantum hardware for scalable distribution loading, which is the basis of a wide range of algorithms that provide quantum advantage.
Preparing matrix product states (MPS) on a quantum device is a key subroutine in many quantum algorithms. The most competitive methods, based on the renormalisation group, prepare translationally invariant MPS of size LLL and bond dimension χχχ, up to an error ε\varepsilonε, in circuit depth O~(χ4log⁡(L/ε))\tilde O(χ^{4}\log(L/\varepsilon))O~(χ4log(L/ε)) or O~(χ6log⁡log⁡(L/ε))\tilde O(χ^{6}\log\log(L/\varepsilon))O~(χ6loglog(L/ε)). We improve multiple aspects of these methods. First, using block-encoded correction maps, whose post-selection succeeds with constant probability, we render the preparation exact without sacrificing the scaling in LLL. Second, through a generalisation of oblivious amplitude amplification to isometries, we reduce the bond-dimension dependence, improving the depth to O~(χ2log⁡L+χ4)\tilde O(χ^{2}\log L + χ^{4})O~(χ2logL+χ4) or O~(χ2log⁡log⁡L+χ4)\tilde O(χ^{2}\log\log L + χ^{4})O~(χ2loglogL+χ4), and even to O~(χ3log⁡L)\tilde O(χ^{3}\log L)O~(χ3logL) for incoherent preparations. Finally, we extend the framework to non-translationally invariant MPS and prove logarithmic-depth exact preparation for independent and identically distributed random tensor sequences. Confirmed by numerical studies, these results constitute, to the best of our knowledge, the most efficient exact MPS preparation protocols in the relevant parameter regimes.
The quantum state preparation of probability distributions is an important subroutine for many quantum algorithms. When embedding D -dimensional multivariate probability distributions by discretizing each dimension into 2n points, we need a state preparation circuit comprising a total of nD qubits, which is often difficult to compile. In this study, we propose a scalable method to generate state preparation circuits for D -dimensional multivariate normal distributions, utilizing tree tensor networks (TTN). We establish theoretical guarantees that multivariate normal distributions with 1D correlation structures can be efficiently represented using TTN. Based on these analyses, we propose a compilation method that uses automatic structural optimization to find the most efficient network structure and compact circuit. We apply our method to state preparation circuits for various high-dimensional random multivariate normal distributions. The numerical results suggest that our method can dramatically reduce the circuit depth and CNOT count while maintaining fidelity compared to existing approaches.
Effective quantum computation relies upon making good use of the exponential information capacity of a quantum machine. A large barrier to designing quantum algorithms for execution on real quantum machines is that, in general, it is intractably difficult to construct an arbitrary quantum state to high precision. Many quantum algorithms rely instead upon initializing the machine in a simple state, and evolving the state through an efficient (i.e. at most polynomial-depth) quantum algorithm. In this work, we show that there exist families of quantum states that can be prepared to high precision with circuits of linear size and depth. We focus on real-valued, smooth, differentiable functions with bounded derivatives on a domain of interest, exemplified by commonly used probability distributions. We further develop an algorithm that requires only linear classical computation time to generate accurate linear- depth circuits to prepare these states, and apply this to well-known and heavily-utilized functions including Gaussian and lognormal distributions. Our procedure rests upon the quantum state representation tool known as the matrix product state (MPS). By efficiently and scalably encoding an explicit amplitude function into an MPS, a high fidelity, linear-depth circuit can directly be generated. These results enable the execution of many quantum algorithms that, aside from initialization, are otherwise depth-efficient.