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:2405.05163 · The European Physical Journal Plus

Fast Fourier transforms and fast Wigner and Weyl functions in large quantum systems

Chao Lei, Apostol Vourdas

Published May 8, 2024• 0 citations
arXivHTMLPDF

Abstract

Abstract Two methods for fast Fourier transforms are used in a quantum context. The first method is for systems with dimension of the Hilbert space $$D=d^n$$ D = d n with d an odd integer, and is inspired by the Cooley-Tukey formalism. The ‘large Fourier transform’ is expressed as a sequence of n ‘small Fourier transforms’ (together with some other transforms) in quantum systems with d-dimensional Hilbert space. Limitations of the method are discussed. In some special cases, the n Fourier transforms can be performed in parallel. The second method is for systems with dimension of the Hilbert space $$D=d_0...d_{n-1}$$ D = d 0 . . . d n - 1 with $$d_0,...,d_{n-1}$$ d 0 , . . . , d n - 1 odd integers coprime to each other. It is inspired by the Good formalism, which in turn is based on the Chinese reminder theorem. In this case also the ‘large Fourier transform’ is expressed as a sequence of n ‘small Fourier transforms’ (that involve some constants related to the number theory that describes the formalism). The ‘small Fourier transforms’ can be performed in a classical computer or in a quantum computer (in which case we have the additional well known advantages of quantum Fourier transform circuits). In the case that the small Fourier transforms are performed with a classical computer, complexity arguments for both methods show the reduction in computational time from $$\mathcal{O}(D^2)$$ O ( D 2 ) to $$\mathcal{O}(D\log D)$$ O ( D log D ) . The second method is also used for the fast calculation of Wigner and Weyl functions, in quantum systems with large finite dimension of the Hilbert space.

DOI
10.1140/epjp/s13360-024-05184-1

Related papers

Fast Quantum Fourier Transforms for a Class of Non-abelian Groups

Markus Püschel, Martin Roetteler, Thomas Beth

An algorithm is presented allowing the construction of fast Fourier transforms for any solvable group on a classical computer. The special structure of the recursion formula being the core of this algorithm makes it a good starting point to obtain systematically fast Fourier transforms for solvable groups on a quantum computer. The inherent structure of the Hilbert space imposed by the qubit architecture suggests to consider groups of order 2^n first (where n is the number of qubits). As an example, fast quantum Fourier transforms for all 4 classes of non-abelian 2-groups with cyclic normal subgroup of index 2 are explicitly constructed in terms of quantum circuits. The (quantum) complexity of the Fourier transform for these groups of size 2^n is O(n^2) in all cases.

arXiv:quant-ph/9807064•Jul 22, 1998•0 citations

Fast quantum nD Fourier and radon transforms

Valeri G. Labunets, Ekaterina V. Labunets-Rundblad, Jaakko T. Astola

Fast Classical and quantum algorithms are introduced for a wide class of non-separable nD discrete unitary K- transforms(DKT)KNn. They require a number of 1D DKT Kn smaller than in the Cooley-Tukey radix-p FFT-type approach. The method utilizes a decomposition of the nDK- transform into a product of original nD discrete Radon Transform and of a family parallel/independ 1DK-transforms. If the nDK-transform has a separable kernel, that again in this case our approach leads to decrease of multiplicative complexity by factor of n compared to the tow/column separable Cooley-Tukey p-radix approach.

Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIE•Jul 26, 2001•1 citation

Quantum systems with finite Hilbert space

Apostol Vourdas

Quantum systems with finite Hilbert space are considered, and phase-space methods like the Heisenberg–Weyl group, symplectic transformations and Wigner and Weyl functions are discussed. A factorization of such systems in terms of smaller subsystems, based on the Chinese remainder theorem, is studied. The general formalism is applied to the case of angular momentum. In this context, SU (2) coherent states are used for analytic representations. Links between the theory of finite quantum systems and other fields of research are discussed.

Reports on Progress in Physics•Feb 2, 2004•321 citations

On the Fast Fourier Transform on SU(2)

Julio Delgado, Alejandro Umaña

The special unitary group SU(2) plays a fundamental role in the description of symmetries in quantum mechanics, theoretical physics, and spherical signal processing. In this paper, we address the computational challenges of performing spectral analysis on this non-abelian compact Lie group. We present the Fourier Transform (FT) on SU(2) and develop a Fast Fourier Transform (FFT) algorithm inspired by the classical Cooley-Tukey divide-and-conquer scheme. Our approach efficiently discretizes the group using Euler angles, applying a two-dimensional FFT on the angular variables and exploiting the recursive properties of Jacobi polynomials. We provide an analysis of the computational complexity, demonstrating that our FFT-based method significantly outperforms the direct computation of the FT. This algorithm serves as a foundational tool for understanding the implementation of the FFT on SU(2), a key component in numerical simulations and advanced data analysis for high-performance computing applications on curved manifolds and quantum systems.

physics.comp-phCross-listedarXiv:2605.23923•Apr 21, 2026

Unitarily inequivalent local and global Fourier transforms in multipartite quantum systems

Chao Lei, Apostol Vourdas

Abstract A multipartite system comprised of n subsystems, each of which is described with ‘local variables’ in $${{\mathbb {Z}}}(d)$$ Z ( d ) and with a d-dimensional Hilbert space H(d), is considered. Local Fourier transforms in each subsystem are defined and related phase space methods are discussed (displacement operators, Wigner and Weyl functions, etc). A holistic view of the same system might be more appropriate in the case of strong interactions, which uses ‘global variables’ in $${{\mathbb {Z}}}(d^n)$$ Z ( d n ) and a $$d^n$$ d n -dimensional Hilbert space $$H(d^n)$$ H ( d n ) . A global Fourier transform is then defined and related phase space methods are discussed. The local formalism is compared and contrasted with the global formalism. Depending on the values of d, n the local Fourier transform is unitarily inequivalent or unitarily equivalent to the global Fourier transform. Time evolution of the system in terms of both local and global variables, is discussed. The formalism can be useful in the general area of Fast Fourier transforms.

arXiv:2301.12137•Quantum Information Processing•Jan 19, 2023•2 citations

Interesting Open Problem Related to Complexity of Computing the Fourier Transform and Group Theory

Nir Ailon

The Fourier Transform is one of the most important linear transformations used in science and engineering. Cooley and Tukey's Fast Fourier Transform (FFT) from 1964 is a method for computing this transformation in time $O(n\log n)$. From a lower bound perspective, relatively little is known. Ailon shows in 2013 an $Ω(n\log n)$ bound for computing the normalized Fourier Transform assuming only unitary operations on pairs of coordinates is allowed. The goal of this document is to describe a natural open problem that arises from this work, which is related to group theory, and in particular to representation theory.

arXiv:1907.07471•Jul 17, 2019•0 citations
Abstract Two methods for fast Fourier transforms are used in a quantum context. The first method is for systems with dimension of the Hilbert space D=dnD=d^nD=dn D = d n with d an odd integer, and is inspired by the Cooley-Tukey formalism. The ‘large Fourier transform’ is expressed as a sequence of n ‘small Fourier transforms’ (together with some other transforms) in quantum systems with d-dimensional Hilbert space. Limitations of the method are discussed. In some special cases, the n Fourier transforms can be performed in parallel. The second method is for systems with dimension of the Hilbert space D=d0...dn−1D=d_0...d_{n-1}D=d0​...dn−1​ D = d 0 . . . d n - 1 with d0,...,dn−1d_0,...,d_{n-1}d0​,...,dn−1​ d 0 , . . . , d n - 1 odd integers coprime to each other. It is inspired by the Good formalism, which in turn is based on the Chinese reminder theorem. In this case also the ‘large Fourier transform’ is expressed as a sequence of n ‘small Fourier transforms’ (that involve some constants related to the number theory that describes the formalism). The ‘small Fourier transforms’ can be performed in a classical computer or in a quantum computer (in which case we have the additional well known advantages of quantum Fourier transform circuits). In the case that the small Fourier transforms are performed with a classical computer, complexity arguments for both methods show the reduction in computational time from O(D2)\mathcal{O}(D^2)O(D2) O ( D 2 ) to O(Dlog⁡D)\mathcal{O}(D\log D)O(DlogD) O ( D log D ) . The second method is also used for the fast calculation of Wigner and Weyl functions, in quantum systems with large finite dimension of the Hilbert space.
An algorithm is presented allowing the construction of fast Fourier transforms for any solvable group on a classical computer. The special structure of the recursion formula being the core of this algorithm makes it a good starting point to obtain systematically fast Fourier transforms for solvable groups on a quantum computer. The inherent structure of the Hilbert space imposed by the qubit architecture suggests to consider groups of order 2^n first (where n is the number of qubits). As an example, fast quantum Fourier transforms for all 4 classes of non-abelian 2-groups with cyclic normal subgroup of index 2 are explicitly constructed in terms of quantum circuits. The (quantum) complexity of the Fourier transform for these groups of size 2^n is O(n^2) in all cases.
Fast Classical and quantum algorithms are introduced for a wide class of non-separable nD discrete unitary K- transforms(DKT)KNn. They require a number of 1D DKT Kn smaller than in the Cooley-Tukey radix-p FFT-type approach. The method utilizes a decomposition of the nDK- transform into a product of original nD discrete Radon Transform and of a family parallel/independ 1DK-transforms. If the nDK-transform has a separable kernel, that again in this case our approach leads to decrease of multiplicative complexity by factor of n compared to the tow/column separable Cooley-Tukey p-radix approach.
Quantum systems with finite Hilbert space are considered, and phase-space methods like the Heisenberg–Weyl group, symplectic transformations and Wigner and Weyl functions are discussed. A factorization of such systems in terms of smaller subsystems, based on the Chinese remainder theorem, is studied. The general formalism is applied to the case of angular momentum. In this context, SU (2) coherent states are used for analytic representations. Links between the theory of finite quantum systems and other fields of research are discussed.
The special unitary group SU(2) plays a fundamental role in the description of symmetries in quantum mechanics, theoretical physics, and spherical signal processing. In this paper, we address the computational challenges of performing spectral analysis on this non-abelian compact Lie group. We present the Fourier Transform (FT) on SU(2) and develop a Fast Fourier Transform (FFT) algorithm inspired by the classical Cooley-Tukey divide-and-conquer scheme. Our approach efficiently discretizes the group using Euler angles, applying a two-dimensional FFT on the angular variables and exploiting the recursive properties of Jacobi polynomials. We provide an analysis of the computational complexity, demonstrating that our FFT-based method significantly outperforms the direct computation of the FT. This algorithm serves as a foundational tool for understanding the implementation of the FFT on SU(2), a key component in numerical simulations and advanced data analysis for high-performance computing applications on curved manifolds and quantum systems.
Abstract A multipartite system comprised of n subsystems, each of which is described with ‘local variables’ in Z(d){{\mathbb {Z}}}(d)Z(d) Z ( d ) and with a d-dimensional Hilbert space H(d), is considered. Local Fourier transforms in each subsystem are defined and related phase space methods are discussed (displacement operators, Wigner and Weyl functions, etc). A holistic view of the same system might be more appropriate in the case of strong interactions, which uses ‘global variables’ in Z(dn){{\mathbb {Z}}}(d^n)Z(dn) Z ( d n ) and a dnd^ndn d n -dimensional Hilbert space H(dn)H(d^n)H(dn) H ( d n ) . A global Fourier transform is then defined and related phase space methods are discussed. The local formalism is compared and contrasted with the global formalism. Depending on the values of d, n the local Fourier transform is unitarily inequivalent or unitarily equivalent to the global Fourier transform. Time evolution of the system in terms of both local and global variables, is discussed. The formalism can be useful in the general area of Fast Fourier transforms.
The Fourier Transform is one of the most important linear transformations used in science and engineering. Cooley and Tukey's Fast Fourier Transform (FFT) from 1964 is a method for computing this transformation in time O(nlog⁡n)O(n\log n)O(nlogn). From a lower bound perspective, relatively little is known. Ailon shows in 2013 an Ω(nlog⁡n)Ω(n\log n)Ω(nlogn) bound for computing the normalized Fourier Transform assuming only unitary operations on pairs of coordinates is allowed. The goal of this document is to describe a natural open problem that arises from this work, which is related to group theory, and in particular to representation theory.