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

Simulation of quantum algorithms using classical probabilistic bits and circuits

Deniz D. Yavuz, Anirudh Yadav

Published Jul 26, 2023• 0 citations
arXivHTMLPDF

Abstract

We discuss a new approach to simulate quantum algorithms using classical probabilistic bits and circuits. Each qubit (a two-level quantum system) is initially mapped to a vector in an eight dimensional probability space (equivalently, to a classical random variable with eight probabilistic outcomes). The key idea in this mapping is to store both the amplitude and phase information of the complex coefficients that describe the qubit state in the probabilities. Due to the identical tensor product structure of combining multiple quantum systems as well as multiple probability spaces, $n$ qubits are then mapped to a tensor product of $n$ 8-dimensional probabilistic vectors (i.e., the Hilbert space of dimension $2^n$ is mapped to a probability space of dimension $8^n$). After this initial mapping, we show how to implement the analogs of single-qubit and two-qubit gates in the probability space using correlation-inducing operations on these classical random variables. The key defining feature of both the mapping to the probability space and the transformations in this space (i.e., operations on the random variables) is that they are not linear, but instead affine. Using this architecture, the evolution of the $2^n$ complex coefficients of the quantum system can be tracked in the joint fully-correlated probabilities of the polynomial number of random variables. We then give specific procedures for implementing (1) the Deutsch-Jozsa algorithm, and (2) the Quantum Fourier Transform in the probability space. Identical to the Quantum case, simulating the Quantum Fourier Transform in the probability space requires $O(n)$ probabilistic bits and $O(n^2)$ (i.e., quadratic in the number of quantum bits) operations.

DOI
10.48550/arxiv.2307.14452

Related papers

Mapping of Quantum Systems to the Probability Simplex

Deniz D. Yavuz, Anirudh Yadav

We start with the simplest quantum system (a two-level system, i.e., a qubit) and discuss a one-to-one mapping of the quantum state in a two-dimensional Hilbert space to a vector in an eight dimensional probability space (probability simplex). We then show how the usual transformations of the quantum state, specifically the Hadamard gate and the single-qubit phase gate, can be accomplished with appropriate transformations of the mapped vector in the probability simplex. One key defining feature of both the mapping to the simplex and the transformations in the simplex is that they are not linear. These results show that both the initial state and the time evolution of a qubit can be fully captured in an eight dimensional probability simplex (or equivalently using three classical probabilistic bits). We then discuss multi-partite quantum systems and their mapping to the probability simplex. Here, the key tool is the identical tensor product structure of combining multiple quantum systems as well as multiple probability spaces. Specifically, we explicitly show how to implement an analog of the two-qubit controlled-not (CNOT) gate in the simplex. We leave it an open problem how much the quantum dynamics of $N$ qubits can be captured in a probability simplex with $3N$ classical probabilistic bits. Finally, we also discuss the equivalent of the Schrodinger's equation for the wavefunction (in a Hilbert space of arbitrary dimension), which dictates the time evolution of the vectors in the simplex.

arXiv:2301.06572•Jan 16, 2023•0 citations

Lecture Notes on Quantum Algorithms

Muhammad Faryad

These notes begin in Chapter 1 with a review of linear algebra and the postulates of quantum mechanics, leading to an explanation of single- and multi-qubit gates. Chapter 2 explores the challenge of constructing arbitrary quantum states from given initial states, and introduces circuits for building oracles. Chapter 3 presents foundational algorithms such as entanglement creation, quantum teleportation, Deutsch-Jozsa, Bernstein-Vazirani, and Simon's algorithm. Chapters 4 and 5 cover algorithms based on the quantum Fourier transform, including phase estimation, period finding, factoring, and logarithm computation. These chapters also include complexity analysis and detailed quantum circuits suitable for implementation in code. Chapter 6 introduces Grover's algorithm for quantum search and amplitude amplification, including its realization via Hamiltonian simulation and a method for derandomization. Chapter 7 discusses basic techniques for Hamiltonian simulation, such as Lie-Trotter decomposition, sparse Hamiltonians, and the linear combination of unitaries. It also provides example circuits for simulating Hamiltonians expressed as linear combinations of Pauli operators. Chapter 8 introduces variational quantum algorithms, and Chapter 9 presents an algorithm for simulating fermionic many-particle systems, with an emphasis on molecular Hamiltonians. It also outlines the key transformations needed to map a molecular Hamiltonian to a form suitable for simulation on a quantum computer.

arXiv:2507.11565•Jul 14, 2025•0 citations

Embedding classical dynamics in a quantum computer

Dimitrios Giannakis, A. Ourmazd, Philipp Pfeffer +2 more

We develop a framework for simulating measure-preserving, ergodic dynamical systems on a quantum computer. Our approach provides an operator-theoretic representation of classical dynamics by combining ergodic theory with quantum information science. The resulting quantum embedding of classical dynamics (QECD) enables efficient simulation of spaces of classical observables with exponentially large dimension using a quadratic number of quantum gates. The QECD framework is based on a quantum feature map that we introduce for representing classical states by density operators on a reproducing kernel Hilbert space, $\mathcal{H}$. Furthermore, an embedding of classical observables into self-adjoint operators on $\mathcal{H}$ is established, such that quantum mechanical expectation values are consistent with pointwise function evaluation. In this scheme, quantum states and observables evolve unitarily under the lifted action of Koopman evolution operators of the classical system. Moreover, by virtue of the reproducing property of $\mathcal{H}$, the quantum system is pointwise-consistent with the underlying classical dynamics. To achieve a quantum computational advantage, we project the state of the quantum system onto a finite-rank density operator on a ${2}^{n}$-dimensional tensor product Hilbert space associated with $n$ qubits. By employing discrete Fourier-Walsh transforms of spectral functions, the evolution operator of the finite-dimensional quantum system is factorized into tensor product form, enabling implementation through an $n$-channel quantum circuit of size $O(n)$ and no interchannel communication. Furthermore, the circuit features a state preparation stage, also of size $O(n)$, and a quantum Fourier transform stage of size $O({n}^{2})$, which makes predictions of observables possible by measurement in the standard computational basis. We prove theoretical convergence results for these predictions in the large-qubit limit, $n\ensuremath{\rightarrow}\ensuremath{\infty}$. In light of these properties, QECD provides a consistent simulator of the evolution of classical observables, realized through projective quantum measurement, which is able to simulate spaces of classical observables of dimension ${2}^{n}$ using circuits of size $O({n}^{2})$. We demonstrate the consistency of the scheme in prototypical dynamical systems involving periodic and quasiperiodic oscillators on tori. These examples include simulated quantum circuit experiments in Qiskit Aer, as well as actual experiments on the IBM Quantum System One.

arXiv:2012.06097•Physical Review A•May 3, 2022•35 citations

Quantum compiler for classical dynamics

Dimitrios Giannakis, A. Ourmazd, Joanna Sławińska +1 more

We develop a framework for simulating measure-preserving, ergodic dynamical systems on a quantum computer. Our approach provides a new operator-theoretic representation of classical dynamics by combining ergodic theory with quantum information science. The resulting \emph{quantum compiler} enables efficient simulation of spaces of classical observables with exponentially large dimension using a quadratic number of quantum gates. The compiler is based on a quantum feature map that we introduce for representing classical states by density operators on a reproducing kernel Hilbert space, $\mathcal H$. Furthermore, an embedding of classical observables into self-adjoint operators on $\mathcal H$ is established, such that quantum mechanical expectation values are consistent with pointwise function evaluation. In this scheme, quantum states and observables evolve unitarily under the lifted action of the Koopman evolution operators of the classical system. Moreover, by virtue of the reproducing property of $\mathcal H$, the quantum system is pointwise-consistent with the underlying classical dynamics. To achieve an exponential quantum computational advantage, we project the state of the quantum system onto a finite-rank density operator on a $2^n$-dimensional tensor product Hilbert space associated with $n$ qubits. By employing discrete Fourier-Walsh transforms of spectral functions, the evolution operator of the finite-dimensional quantum system is factorized into tensor product form, enabling implementation through an $n$-channel quantum circuit with an $O(n)$ gates and no interchannel communication. Furthermore, the circuit features a quantum Fourier transform stage with $O(n^2)$ gates, which generates predictions by measurement in the standard computational basis. We implement the framework in Qiskit, and demonstrate its consistency for prototypical periodic and quasiperiodic systems.

Bulletin of the American Physical Society•Dec 11, 2020•0 citations

Quantum computing with classical bits

CHRISTOF WETTERICH

A bit-quantum map relates probabilistic information for Ising spins or classical bits to quantum spins or qubits. Quantum systems are subsystems of classical statistical systems. The Ising spins can represent macroscopic two-level observables, and the quantum subsystem employs suitable expectation values and correlations. We discuss static memory materials based on Ising spins for which boundary information can be transported through the bulk in a generalized equilibrium state. They can realize quantum operations as the Hadamard or CNOT-gate for the quantum subsystem. Classical probabilistic systems can account for the entanglement of quantum spins. An arbitrary unitary evolution for an arbitrary number of quantum spins can be described by static memory materials for an infinite number of Ising spins which may, in turn, correspond to continuous variables or fields. We discuss discrete subsets of unitary operations realized by a finite number of Ising spins. They may be useful for new computational structures. We suggest that features of quantum computation or more general probabilistic computation may be realized by neural networks, neuromorphic computing or perhaps even the brain. This does neither require small isolated entities nor very low temperature. In these systems the processing of probabilistic information can be more general than for static memory materials. We propose a general formalism for probabilistic computing for which deterministic computing and quantum computing are special limiting cases.

arXiv:1806.05960•Nuclear Physics B•Sep 26, 2019•14 citations

Stochastic emulation of quantum algorithms

Daniel Braun, Ronny Müller

Abstract Quantum algorithms profit from the interference of quantum states in an exponentially large Hilbert space and the fact that unitary transformations on that Hilbert space can be broken down to universal gates that act only on one or two qubits at the same time. The former aspect renders the direct classical simulation of quantum algorithms difficult. Here we introduce higher-order partial derivatives of a probability distribution of particle positions as a new object that shares these basic properties of quantum mechanical states needed for a quantum algorithm. Discretization of the positions allows one to represent the quantum mechanical state of n bit qubits by 2(n bit + 1) classical stochastic bits. Based on this, we demonstrate many-particle interference and representation of pure entangled quantum states via derivatives of probability distributions and find the universal set of stochastic maps that correspond to the quantum gates in a universal gate set. We prove that the propagation via the stochastic map built from those universal stochastic maps reproduces up to a prefactor exactly the evolution of the quantum mechanical state with the corresponding quantum algorithm, leading to an automated translation of a quantum algorithm to a stochastic classical algorithm. We implement several well-known quantum algorithms, analyse the scaling of the needed number of realizations with the number of qubits, and highlight the role of destructive interference for the cost of the emulation.

arXiv:2109.07777•New Journal of Physics•Jan 13, 2022•1 citation
We discuss a new approach to simulate quantum algorithms using classical probabilistic bits and circuits. Each qubit (a two-level quantum system) is initially mapped to a vector in an eight dimensional probability space (equivalently, to a classical random variable with eight probabilistic outcomes). The key idea in this mapping is to store both the amplitude and phase information of the complex coefficients that describe the qubit state in the probabilities. Due to the identical tensor product structure of combining multiple quantum systems as well as multiple probability spaces, nnn qubits are then mapped to a tensor product of nnn 8-dimensional probabilistic vectors (i.e., the Hilbert space of dimension 2n2^n2n is mapped to a probability space of dimension 8n8^n8n). After this initial mapping, we show how to implement the analogs of single-qubit and two-qubit gates in the probability space using correlation-inducing operations on these classical random variables. The key defining feature of both the mapping to the probability space and the transformations in this space (i.e., operations on the random variables) is that they are not linear, but instead affine. Using this architecture, the evolution of the 2n2^n2n complex coefficients of the quantum system can be tracked in the joint fully-correlated probabilities of the polynomial number of random variables. We then give specific procedures for implementing (1) the Deutsch-Jozsa algorithm, and (2) the Quantum Fourier Transform in the probability space. Identical to the Quantum case, simulating the Quantum Fourier Transform in the probability space requires O(n)O(n)O(n) probabilistic bits and O(n2)O(n^2)O(n2) (i.e., quadratic in the number of quantum bits) operations.
We start with the simplest quantum system (a two-level system, i.e., a qubit) and discuss a one-to-one mapping of the quantum state in a two-dimensional Hilbert space to a vector in an eight dimensional probability space (probability simplex). We then show how the usual transformations of the quantum state, specifically the Hadamard gate and the single-qubit phase gate, can be accomplished with appropriate transformations of the mapped vector in the probability simplex. One key defining feature of both the mapping to the simplex and the transformations in the simplex is that they are not linear. These results show that both the initial state and the time evolution of a qubit can be fully captured in an eight dimensional probability simplex (or equivalently using three classical probabilistic bits). We then discuss multi-partite quantum systems and their mapping to the probability simplex. Here, the key tool is the identical tensor product structure of combining multiple quantum systems as well as multiple probability spaces. Specifically, we explicitly show how to implement an analog of the two-qubit controlled-not (CNOT) gate in the simplex. We leave it an open problem how much the quantum dynamics of NNN qubits can be captured in a probability simplex with 3N3N3N classical probabilistic bits. Finally, we also discuss the equivalent of the Schrodinger's equation for the wavefunction (in a Hilbert space of arbitrary dimension), which dictates the time evolution of the vectors in the simplex.
These notes begin in Chapter 1 with a review of linear algebra and the postulates of quantum mechanics, leading to an explanation of single- and multi-qubit gates. Chapter 2 explores the challenge of constructing arbitrary quantum states from given initial states, and introduces circuits for building oracles. Chapter 3 presents foundational algorithms such as entanglement creation, quantum teleportation, Deutsch-Jozsa, Bernstein-Vazirani, and Simon's algorithm. Chapters 4 and 5 cover algorithms based on the quantum Fourier transform, including phase estimation, period finding, factoring, and logarithm computation. These chapters also include complexity analysis and detailed quantum circuits suitable for implementation in code. Chapter 6 introduces Grover's algorithm for quantum search and amplitude amplification, including its realization via Hamiltonian simulation and a method for derandomization. Chapter 7 discusses basic techniques for Hamiltonian simulation, such as Lie-Trotter decomposition, sparse Hamiltonians, and the linear combination of unitaries. It also provides example circuits for simulating Hamiltonians expressed as linear combinations of Pauli operators. Chapter 8 introduces variational quantum algorithms, and Chapter 9 presents an algorithm for simulating fermionic many-particle systems, with an emphasis on molecular Hamiltonians. It also outlines the key transformations needed to map a molecular Hamiltonian to a form suitable for simulation on a quantum computer.
We develop a framework for simulating measure-preserving, ergodic dynamical systems on a quantum computer. Our approach provides an operator-theoretic representation of classical dynamics by combining ergodic theory with quantum information science. The resulting quantum embedding of classical dynamics (QECD) enables efficient simulation of spaces of classical observables with exponentially large dimension using a quadratic number of quantum gates. The QECD framework is based on a quantum feature map that we introduce for representing classical states by density operators on a reproducing kernel Hilbert space, H\mathcal{H}H. Furthermore, an embedding of classical observables into self-adjoint operators on H\mathcal{H}H is established, such that quantum mechanical expectation values are consistent with pointwise function evaluation. In this scheme, quantum states and observables evolve unitarily under the lifted action of Koopman evolution operators of the classical system. Moreover, by virtue of the reproducing property of H\mathcal{H}H, the quantum system is pointwise-consistent with the underlying classical dynamics. To achieve a quantum computational advantage, we project the state of the quantum system onto a finite-rank density operator on a 2n{2}^{n}2n-dimensional tensor product Hilbert space associated with nnn qubits. By employing discrete Fourier-Walsh transforms of spectral functions, the evolution operator of the finite-dimensional quantum system is factorized into tensor product form, enabling implementation through an nnn-channel quantum circuit of size O(n)O(n)O(n) and no interchannel communication. Furthermore, the circuit features a state preparation stage, also of size O(n)O(n)O(n), and a quantum Fourier transform stage of size O(n2)O({n}^{2})O(n2), which makes predictions of observables possible by measurement in the standard computational basis. We prove theoretical convergence results for these predictions in the large-qubit limit, $n\ensuremath{\rightarrow}\ensuremath{\infty}$. In light of these properties, QECD provides a consistent simulator of the evolution of classical observables, realized through projective quantum measurement, which is able to simulate spaces of classical observables of dimension 2n{2}^{n}2n using circuits of size O(n2)O({n}^{2})O(n2). We demonstrate the consistency of the scheme in prototypical dynamical systems involving periodic and quasiperiodic oscillators on tori. These examples include simulated quantum circuit experiments in Qiskit Aer, as well as actual experiments on the IBM Quantum System One.
We develop a framework for simulating measure-preserving, ergodic dynamical systems on a quantum computer. Our approach provides a new operator-theoretic representation of classical dynamics by combining ergodic theory with quantum information science. The resulting \emph{quantum compiler} enables efficient simulation of spaces of classical observables with exponentially large dimension using a quadratic number of quantum gates. The compiler is based on a quantum feature map that we introduce for representing classical states by density operators on a reproducing kernel Hilbert space, H\mathcal HH. Furthermore, an embedding of classical observables into self-adjoint operators on H\mathcal HH is established, such that quantum mechanical expectation values are consistent with pointwise function evaluation. In this scheme, quantum states and observables evolve unitarily under the lifted action of the Koopman evolution operators of the classical system. Moreover, by virtue of the reproducing property of H\mathcal HH, the quantum system is pointwise-consistent with the underlying classical dynamics. To achieve an exponential quantum computational advantage, we project the state of the quantum system onto a finite-rank density operator on a 2n2^n2n-dimensional tensor product Hilbert space associated with nnn qubits. By employing discrete Fourier-Walsh transforms of spectral functions, the evolution operator of the finite-dimensional quantum system is factorized into tensor product form, enabling implementation through an nnn-channel quantum circuit with an O(n)O(n)O(n) gates and no interchannel communication. Furthermore, the circuit features a quantum Fourier transform stage with O(n2)O(n^2)O(n2) gates, which generates predictions by measurement in the standard computational basis. We implement the framework in Qiskit, and demonstrate its consistency for prototypical periodic and quasiperiodic systems.
A bit-quantum map relates probabilistic information for Ising spins or classical bits to quantum spins or qubits. Quantum systems are subsystems of classical statistical systems. The Ising spins can represent macroscopic two-level observables, and the quantum subsystem employs suitable expectation values and correlations. We discuss static memory materials based on Ising spins for which boundary information can be transported through the bulk in a generalized equilibrium state. They can realize quantum operations as the Hadamard or CNOT-gate for the quantum subsystem. Classical probabilistic systems can account for the entanglement of quantum spins. An arbitrary unitary evolution for an arbitrary number of quantum spins can be described by static memory materials for an infinite number of Ising spins which may, in turn, correspond to continuous variables or fields. We discuss discrete subsets of unitary operations realized by a finite number of Ising spins. They may be useful for new computational structures. We suggest that features of quantum computation or more general probabilistic computation may be realized by neural networks, neuromorphic computing or perhaps even the brain. This does neither require small isolated entities nor very low temperature. In these systems the processing of probabilistic information can be more general than for static memory materials. We propose a general formalism for probabilistic computing for which deterministic computing and quantum computing are special limiting cases.
Abstract Quantum algorithms profit from the interference of quantum states in an exponentially large Hilbert space and the fact that unitary transformations on that Hilbert space can be broken down to universal gates that act only on one or two qubits at the same time. The former aspect renders the direct classical simulation of quantum algorithms difficult. Here we introduce higher-order partial derivatives of a probability distribution of particle positions as a new object that shares these basic properties of quantum mechanical states needed for a quantum algorithm. Discretization of the positions allows one to represent the quantum mechanical state of n bit qubits by 2(n bit + 1) classical stochastic bits. Based on this, we demonstrate many-particle interference and representation of pure entangled quantum states via derivatives of probability distributions and find the universal set of stochastic maps that correspond to the quantum gates in a universal gate set. We prove that the propagation via the stochastic map built from those universal stochastic maps reproduces up to a prefactor exactly the evolution of the quantum mechanical state with the corresponding quantum algorithm, leading to an automated translation of a quantum algorithm to a stochastic classical algorithm. We implement several well-known quantum algorithms, analyse the scaling of the needed number of realizations with the number of qubits, and highlight the role of destructive interference for the cost of the emulation.