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:2601.17724

Quantum-Inspired Algorithms beyond Unitary Circuits: the Laplace Transform

Noufal Jaseem, Sergi Ramos-Calderer, Gauthameshwar S., Dingzu Wang, José Ignacio Latorre, Dario Poletti

quant-phmath-phphysics.data-anSubmitted Jan 25, 2026• Announced Jan 27, 2026• 0 citations• Revised Mar 18, 2026 (v2)
v1 · Jan 25, 2026v2 · Mar 18, 2026v3 · Sep 4, 2026
arXivHTMLPDF

Abstract

Quantum-inspired algorithms can deliver substantial speedups over classical state-of-the-art methods by executing quantum algorithms with tensor networks on conventional hardware. Unlike circuit models restricted to unitary gates, tensor networks naturally accommodate non-unitary maps. This flexibility lets us design quantum-inspired methods that start from a quantum algorithmic structure, yet go beyond unitarity to achieve speedups. Here we introduce a tensor-network approach to compute the discrete Laplace transform, a non-unitary, aperiodic transform (in contrast to the Fourier transform). We encode a length-$N$ signal on two paired $n$-qubit registers and decompose the overall map into a non-unitary exponential Damping Transform followed by a Quantum Fourier Transform, both compressed in a single matrix-product operator. This decomposition admits strong MPO compression to low bond dimension resulting in significant acceleration. We demonstrate simulations up to $N=2^{30}$ input data points, with up to $2^{60}$ output data points, and quantify how bond dimension controls runtime and accuracy, including precise and efficient pole identification.

Comments
9 pages

Related papers

A Polylogarithmic-Time Quantum Algorithm for the Laplace Transform

Akash Kumar Singh, Ashish Kumar Patra, Anurag K. S. V. +3 more

We introduce a quantum algorithm to perform the Laplace transform on quantum computers. Already, the quantum Fourier transform (QFT) is the cornerstone of many quantum algorithms, but the Laplace transform or its discrete version has not seen any efficient implementation on quantum computers due to its dissipative nature and hence non-unitary dynamics. However, a recent work has shown an efficient implementation for certain cases on quantum computers using the Taylor series. Unlike previous work, our work provides a completely different algorithm for doing Laplace Transform using Quantum Eigenvalue Transformation and Lap-LCHS, very efficiently at points which form an arithmetic progression. Our algorithm can implement $N \times N$ discrete Laplace transform in gate complexity that grows as $O((log\,N)^3)$, ignoring the state preparation cost, where $N=2^n$ and $n$ is the number of qubits, which is a superpolynomial speedup in number of gates over the best classical counterpart that has complexity $O(N\cdot log\,N)$ for the same cases. Also, the circuit width grows as $O(log\,N)$. Quantum Laplace Transform (QLT) may enable new Quantum algorithms for cases like solving differential equations in the Laplace domain, developing an inverse Laplace transform algorithm on quantum computers, imaginary time evolution in the resolvent domain for calculating ground state energy, and spectral estimation of non-Hermitian matrices.

quant-pharXiv:2512.17980•Dec 22, 2025•0 citations

Who can compete with quantum computers? Lecture notes on quantum inspired tensor networks computational techniques

Xavier Waintal, Chen-How Huang, Christoph W. Groth

This is a set of lectures on tensor networks with a strong emphasis on the core algorithms involving Matrix Product States (MPS) and Matrix Product Operators (MPO). Compared to other presentations, particular care has been given to disentangle aspects of tensor networks from the quantum many-body problem: MPO/MPS algorithms are presented as a way to deal with linear algebra on extremely (exponentially) large matrices and vectors, regardless of any particular application. The lectures include well-known algorithms to find eigenvectors of MPOs (the celebrated DMRG), solve linear problems, and recent learning algorithms that allow one to map a known function into an MPS (the Tensor Cross Interpolation, or TCI, algorithm). The lectures end with a discussion of how to represent functions and perform calculus with tensor networks using the "quantics" representation. They include the detailed analytical construction of important MPOs such as those for differentiation, indefinite integration, convolution, and the quantum Fourier transform. Three concrete applications are discussed in detail: the simulation of a quantum computer (either exactly or with compression), the simulation of a quantum annealer, and techniques to solve partial differential equations (e.g. Poisson, diffusion, or Gross-Pitaevskii) within the "quantics" representation. The lectures have been designed to be accessible to a first-year PhD student and include detailed proofs of all statements.

quant-pharXiv:2601.03035•SciPost Physics Lecture Notes•Jan 7, 2026•0 citations

A Quantum-Inspired Algorithm for Wave Simulation Using Tensor Networks

Kevin Lively, Vittorio Pagni, Gonzalo Camacho

We present an efficient classical algorithm based on the construction of a unitary quantum circuit for simulating the Isotropic Wave Equation (IWE) in one, two, or three dimensions. Using an analogy with the massless Dirac equation, second order time and space derivatives in the IWE are reduced to first order, resulting in a Schrödinger equation of motion. Exact diagonalization of the unitary circuit in combination with Tensor Networks allows simulation of the wave equation with a resolution of $10^{13}$ grid points on a laptop. A method for encoding arbitrary analytical functions into diagonal Matrix Product Operators is employed to prepare and evolve a Matrix Product State (MPS) encoding the solution. Since the method relies on the Quantum Fourier Transform, which has been shown to generate small entanglement when applied to arbitrary MPSs, simulating the evolution of initial conditions with sufficiently low bond dimensions to high accuracy becomes highly efficient, up to the cost of Trotterized propagation and sampling of the wavefunction. We conclude by discussing possible extensions of the approach for carrying out Tensor Network simulations of other partial differential equations such as Maxwell's equations.

arXiv:2504.11181•Apr 15, 2025•0 citations

A Quantum-Inspired Algorithm for Wave Simulation Using Tensor Networks

Kevin Lively, Vittorio Pagni, Gonzalo Camacho

We present an efficient classical algorithm based on the construction of a unitary quantum circuit for simulating the Isotropic Wave Equation (IWE) in one, two, or three dimensions. Using an analogy with the massless Dirac equation second order time and space derivatives in the IWE are reduced to first order, resulting in a Schrödinger equation of motion. Exact diagonalization of the unitary circuit in combination with Tensor Networks allows simulation of the wave equation with a resolution of 1013grid points on a laptop. A method for encoding arbitrary analytical functions into diagonal Matrix Product Operators (MPOs) is employed to prepare and evolve a Matrix Product State (MPS) encoding the solution. Since the method relies on the Quantum Fourier Transform, which has been shown to generate small entanglement when applied to arbitrary MPSs, simulating the evolution of initial conditions with sufficiently low bond dimensions to high accuracy becomes highly efficient, up to the cost of Trotterized propagation and sampling of the wavefunction. We conclude by discussing possible extensions of the approach for carrying out Tensor Network simulations of other partial differential equations such as Maxwell's equations.

IEEE International Conference on Quantum Computing and Engineering (QCE)•Aug 30, 2025•2 citations

Low-rank tensor decompositions of quantum circuits

Patrick Gelß, Stefan Klus, Knebel, Sebastian +2 more

Quantum computing is arguably one of the most revolutionary and disruptive technologies of this century. Due to the ever-increasing number of potential applications as well as the continuing rise in complexity, the development, simulation, optimization, and physical realization of quantum circuits is of utmost importance for designing novel algorithms. We show how matrix product states (MPSs) and matrix product operators (MPOs) can be used to express certain quantum states, quantum gates, and entire quantum circuits as low-rank tensors. This enables the analysis and simulation of complex quantum circuits on classical computers and to gain insight into the underlying structure of the system. We present different examples to demonstrate the advantages of MPO formulations and show that they are more efficient than conventional techniques if the bond dimensions of the wave function representation can be kept small throughout the simulation.

arXiv:2205.09882•May 19, 2022•1 citation

Complexity of tensor network simulation for noisy quantum circuits

Yuguo Shao, Zishuo Zhao, Song Cheng +1 more

We aim to rigorously address how local noise affects classical simulability of quantum dynamics benchmarked by tensor-network methods. Using operator entanglement entropy (OEE), we prove the following: (1) For single-qubit depolarizing noise on arbitrary circuits, tensor networks with $\mathrm{poly}(n)$ bond dimension suffice for fixed absolute Hilbert-Schmidt error after $\order{1}$ depth, while relative error demands $\order{\log n}$ depth; and this bound is optimal. (2) For single-qubit depolarizing noise on 1D local circuits, the existence of whole-trajectory error-bounded matrix product operator (MPO) of $\mathrm{poly}(n)$ bond dimension at all depths. (3) For general single-qubit noise in 1D brickwall circuits, random two-design gates with contraction coefficient $c<1/3$ yield an $\order{1}$ OEE plateau with probability $1-Te^{-Ω(n)}$, while arbitrary gates with $c<1/48$ give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth $\mathrm{poly}(n)$ average boundary-bond dimensions for projected entangled pair operators~(PEPO) across every cut -- under depolarizing noise at either absolute or relative accuracy, and under general noise with strong contraction at absolute accuracy. Our results establish a rigorous connection between certain noise models, circuit types, and their classical simulability.

cond-mat.str-elCross-listedarXiv:2606.00474•Jun 2, 2026•0 citations
Quantum-inspired algorithms can deliver substantial speedups over classical state-of-the-art methods by executing quantum algorithms with tensor networks on conventional hardware. Unlike circuit models restricted to unitary gates, tensor networks naturally accommodate non-unitary maps. This flexibility lets us design quantum-inspired methods that start from a quantum algorithmic structure, yet go beyond unitarity to achieve speedups. Here we introduce a tensor-network approach to compute the discrete Laplace transform, a non-unitary, aperiodic transform (in contrast to the Fourier transform). We encode a length-NNN signal on two paired nnn-qubit registers and decompose the overall map into a non-unitary exponential Damping Transform followed by a Quantum Fourier Transform, both compressed in a single matrix-product operator. This decomposition admits strong MPO compression to low bond dimension resulting in significant acceleration. We demonstrate simulations up to N=230N=2^{30}N=230 input data points, with up to 2602^{60}260 output data points, and quantify how bond dimension controls runtime and accuracy, including precise and efficient pole identification.
We introduce a quantum algorithm to perform the Laplace transform on quantum computers. Already, the quantum Fourier transform (QFT) is the cornerstone of many quantum algorithms, but the Laplace transform or its discrete version has not seen any efficient implementation on quantum computers due to its dissipative nature and hence non-unitary dynamics. However, a recent work has shown an efficient implementation for certain cases on quantum computers using the Taylor series. Unlike previous work, our work provides a completely different algorithm for doing Laplace Transform using Quantum Eigenvalue Transformation and Lap-LCHS, very efficiently at points which form an arithmetic progression. Our algorithm can implement N×NN \times NN×N discrete Laplace transform in gate complexity that grows as O((log N)3)O((log\,N)^3)O((logN)3), ignoring the state preparation cost, where N=2nN=2^nN=2n and nnn is the number of qubits, which is a superpolynomial speedup in number of gates over the best classical counterpart that has complexity O(N⋅log N)O(N\cdot log\,N)O(N⋅logN) for the same cases. Also, the circuit width grows as O(log N)O(log\,N)O(logN). Quantum Laplace Transform (QLT) may enable new Quantum algorithms for cases like solving differential equations in the Laplace domain, developing an inverse Laplace transform algorithm on quantum computers, imaginary time evolution in the resolvent domain for calculating ground state energy, and spectral estimation of non-Hermitian matrices.
This is a set of lectures on tensor networks with a strong emphasis on the core algorithms involving Matrix Product States (MPS) and Matrix Product Operators (MPO). Compared to other presentations, particular care has been given to disentangle aspects of tensor networks from the quantum many-body problem: MPO/MPS algorithms are presented as a way to deal with linear algebra on extremely (exponentially) large matrices and vectors, regardless of any particular application. The lectures include well-known algorithms to find eigenvectors of MPOs (the celebrated DMRG), solve linear problems, and recent learning algorithms that allow one to map a known function into an MPS (the Tensor Cross Interpolation, or TCI, algorithm). The lectures end with a discussion of how to represent functions and perform calculus with tensor networks using the "quantics" representation. They include the detailed analytical construction of important MPOs such as those for differentiation, indefinite integration, convolution, and the quantum Fourier transform. Three concrete applications are discussed in detail: the simulation of a quantum computer (either exactly or with compression), the simulation of a quantum annealer, and techniques to solve partial differential equations (e.g. Poisson, diffusion, or Gross-Pitaevskii) within the "quantics" representation. The lectures have been designed to be accessible to a first-year PhD student and include detailed proofs of all statements.
We present an efficient classical algorithm based on the construction of a unitary quantum circuit for simulating the Isotropic Wave Equation (IWE) in one, two, or three dimensions. Using an analogy with the massless Dirac equation, second order time and space derivatives in the IWE are reduced to first order, resulting in a Schrödinger equation of motion. Exact diagonalization of the unitary circuit in combination with Tensor Networks allows simulation of the wave equation with a resolution of 101310^{13}1013 grid points on a laptop. A method for encoding arbitrary analytical functions into diagonal Matrix Product Operators is employed to prepare and evolve a Matrix Product State (MPS) encoding the solution. Since the method relies on the Quantum Fourier Transform, which has been shown to generate small entanglement when applied to arbitrary MPSs, simulating the evolution of initial conditions with sufficiently low bond dimensions to high accuracy becomes highly efficient, up to the cost of Trotterized propagation and sampling of the wavefunction. We conclude by discussing possible extensions of the approach for carrying out Tensor Network simulations of other partial differential equations such as Maxwell's equations.
We present an efficient classical algorithm based on the construction of a unitary quantum circuit for simulating the Isotropic Wave Equation (IWE) in one, two, or three dimensions. Using an analogy with the massless Dirac equation second order time and space derivatives in the IWE are reduced to first order, resulting in a Schrödinger equation of motion. Exact diagonalization of the unitary circuit in combination with Tensor Networks allows simulation of the wave equation with a resolution of 1013grid points on a laptop. A method for encoding arbitrary analytical functions into diagonal Matrix Product Operators (MPOs) is employed to prepare and evolve a Matrix Product State (MPS) encoding the solution. Since the method relies on the Quantum Fourier Transform, which has been shown to generate small entanglement when applied to arbitrary MPSs, simulating the evolution of initial conditions with sufficiently low bond dimensions to high accuracy becomes highly efficient, up to the cost of Trotterized propagation and sampling of the wavefunction. We conclude by discussing possible extensions of the approach for carrying out Tensor Network simulations of other partial differential equations such as Maxwell's equations.
Quantum computing is arguably one of the most revolutionary and disruptive technologies of this century. Due to the ever-increasing number of potential applications as well as the continuing rise in complexity, the development, simulation, optimization, and physical realization of quantum circuits is of utmost importance for designing novel algorithms. We show how matrix product states (MPSs) and matrix product operators (MPOs) can be used to express certain quantum states, quantum gates, and entire quantum circuits as low-rank tensors. This enables the analysis and simulation of complex quantum circuits on classical computers and to gain insight into the underlying structure of the system. We present different examples to demonstrate the advantages of MPO formulations and show that they are more efficient than conventional techniques if the bond dimensions of the wave function representation can be kept small throughout the simulation.
We aim to rigorously address how local noise affects classical simulability of quantum dynamics benchmarked by tensor-network methods. Using operator entanglement entropy (OEE), we prove the following: (1) For single-qubit depolarizing noise on arbitrary circuits, tensor networks with poly(n)\mathrm{poly}(n)poly(n) bond dimension suffice for fixed absolute Hilbert-Schmidt error after $\order{1}$ depth, while relative error demands $\order{\log n}$ depth; and this bound is optimal. (2) For single-qubit depolarizing noise on 1D local circuits, the existence of whole-trajectory error-bounded matrix product operator (MPO) of poly(n)\mathrm{poly}(n)poly(n) bond dimension at all depths. (3) For general single-qubit noise in 1D brickwall circuits, random two-design gates with contraction coefficient c<1/3c<1/3c<1/3 yield an $\order{1}$ OEE plateau with probability 1−Te−Ω(n)1-Te^{-Ω(n)}1−Te−Ω(n), while arbitrary gates with c<1/48c<1/48c<1/48 give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth poly(n)\mathrm{poly}(n)poly(n) average boundary-bond dimensions for projected entangled pair operators~(PEPO) across every cut -- under depolarizing noise at either absolute or relative accuracy, and under general noise with strong contraction at absolute accuracy. Our results establish a rigorous connection between certain noise models, circuit types, and their classical simulability.