Papers
Topics
Authors
Recent
Search
2000 character limit reached

Quantum Fourier Transform

Updated 10 January 2026
  • Quantum Fourier Transform is a unitary operation that converts computational basis states into periodic frequency representations, forming the backbone of quantum algorithms.
  • The standard QFT circuit uses Hadamard and controlled-phase gates with O(n²) complexity, while approximations and dynamic circuits reduce resource costs for fault tolerance.
  • QFT underpins landmark protocols such as Shor’s algorithm and phase estimation, driving advances in quantum arithmetic, simulation, and overall quantum computing performance.

The Quantum Fourier Transform (QFT) is a fundamental unitary transformation that underpins a broad array of quantum algorithms, providing the mechanism for extracting periodicity and frequency-domain information from quantum registers. Its mathematical structure is tightly connected to the discrete Fourier transform (DFT), yet it is efficiently implementable on quantum hardware using a sequence of elementary gates. The QFT is a critical component of landmark quantum algorithms such as Shor’s factoring algorithm and phase estimation, and it serves as the core for a variety of quantum protocols, arithmetic subroutines, and Hamiltonian simulation primitives.

1. Definition and Mathematical Structure

For an nn-qubit register, the QFT is defined as the unitary operator QFTn\mathrm{QFT}_n acting on the computational basis {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\} according to

QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.

In matrix terms, [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}. For an arbitrary input state ∑xαx∣x⟩\sum_x\alpha_x\ket{x}, the QFT maps this to ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}, performing an exact discrete Fourier transform on the amplitudes (Bäumer et al., 2024).

2. Circuit Decompositions and Resource Analysis

Standard Circuit Construction: The canonical unitary QFT circuit consists of:

  • For each qubit jj (from most to least significant), apply a Hadamard gate HH, then for each k>jk > j, apply a controlled-phase gate QFTn\mathrm{QFT}_n0 with angle QFTn\mathrm{QFT}_n1.
  • Finally, perform a reversal of the qubit register via swap gates.

The controlled-phase operation QFTn\mathrm{QFT}_n2 is implemented as

QFTn\mathrm{QFT}_n3

For an QFTn\mathrm{QFT}_n4-qubit system, the standard decomposition requires QFTn\mathrm{QFT}_n5 Hadamards and QFTn\mathrm{QFT}_n6 two-qubit controlled-phase gates. Exact gate count is QFTn\mathrm{QFT}_n7; limited connectivity adds SWAP overhead (Bäumer et al., 2024, Camps et al., 2020). The circuit depth is also QFTn\mathrm{QFT}_n8 in the naive sequential implementation.

Approximate QFT and Fault-Tolerance: In fault-tolerant settings, the QFT is typically approximated by dropping controlled-phase gates with angles below a certain threshold QFTn\mathrm{QFT}_n9, yielding diamond-norm error {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}0. The resource-dominant cost is quantified in T gates for Clifford+T synthesis:

  • Standard (banded) AQFT: {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}1-count {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}2, T-depth {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}3.
  • Optimized constructions remove Toffoli gates, halving leading T-count to {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}4, and using parallel adders to reduce T-depth to {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}5 (Park et al., 2022).

3. Dynamic Circuits and Semi-Classical QFT

The dynamic-circuit (Griffiths–Niu) implementation exploits mid-circuit measurements and classical feed-forward, drastically reducing hardware resource requirements for QFTs followed immediately by measurement:

  • The algorithm sequentially applies Hadamards, measures each qubit, and uses the measurement result to classically control subsequent single-qubit phase rotations, bypassing all two-qubit entangling gates.
  • The dynamic QFT needs only {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}6 mid-circuit measurements, {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}7 Hadamards, and {∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}8 classical control operations—no two-qubit gates or connectivity constraints (Bäumer et al., 2024).
  • Empirically, this approach on superconducting hardware (ibm_kyiv) enables QFT{∣x⟩:x∈{0,…,2n−1}}\{\ket{x}:x\in\{0,\dots,2^n-1\}\}9Measurement with certified process fidelity QFTn:∣x⟩⟼12n∑k=02n−1e2Ï€i xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.0 up to QFTn:∣x⟩⟼12n∑k=02n−1e2Ï€i xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.1 and QFTn:∣x⟩⟼12n∑k=02n−1e2Ï€i xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.2 up to QFTn:∣x⟩⟼12n∑k=02n−1e2Ï€i xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.3 qubits; these results exceed all previous multi-qubit QFT reports on any hardware platform. Key to this advance is a tailored dynamical decoupling protocol—feed-forward-compensated dynamical decoupling (FC-DD)—which protects idle qubits during measurement and feed-forward (Bäumer et al., 2024).

4. Alternative Physical and Algorithmic Realizations

Continuous-Variable QFT

The continuous-variable QFT (cvQFT) generalizes the transformation to bosonic mode systems. Mathematically, the cvQFT is implemented as a rotation operator QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.4, where QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.5 is the DFT matrix QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.6: QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.7 The decomposition leverages phase shifters and beam splitters, with dense implementations requiring QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.8 elements, but FFT-inspired designs reduce the cost to QFTn:∣x⟩⟼12n∑k=02n−1e2πi xk/2n∣k⟩.\mathrm{QFT}_n:\quad \ket{x}\longmapsto \frac{1}{\sqrt{2^n}}\sum_{k=0}^{2^n-1} e^{2\pi i\,xk/2^n}\ket{k}.9 (Cariolaro et al., 14 Dec 2025).

Quadratic QFT and Analog Simulation

The quadratic QFT (QQFT), defined on cold-atom lattices via second-quantized bosonic operators, is realized by concatenating local unitaries (on-site potentials and nearest-neighbor tunnelings) programmable by digital micromirror devices. QQFT generalizes QFT to the many-body sector (particle-number conserving), reproducing QFT in the single-particle sector. It enables programmable Hamiltonian engineering, realizing nontrivial models such as 1D Poincaré spacetime crystals and 2D perfectly flat Chern bands in [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}0 steps, robust against both white and colored noise (Wang et al., 2022).

Photonic and Resonator Implementations

QFT protocols have been adapted to photonic systems (OAM encoding, cavity-QED mediated transformations) and to circuit QED via coupled oscillators and engineered cross-Kerr interactions:

  • Photonic OAM QFT achieves an efficient [QFTn]k,x=12ne2Ï€i xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}1 scaling in optical elements for high-dimensional Hilbert spaces.
  • Cavity QED can mediate QFT among photonic qubits using a single atom as the ancilla and adjustable time-delay feedback, with tunable phase gates via Stark shifting (Shi et al., 2021).
  • Superconducting resonator-based QFTs use state transfer and cross-Kerr evolution between two harmonic oscillators to perform the [QFTn]k,x=12ne2Ï€i xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}2-dimensional transform, yielding a fully quantum result without mid-circuit measurements or feed-forward (Chen et al., 2019, Kysela, 2021).

5. Applications and Algorithmic Impact

The QFT is central to the following classes of quantum algorithms and primitives:

  • Period-finding, order-finding, and integer factoring: Shor’s algorithm is built around the QFT to extract hidden periodicities efficiently (Kashani et al., 2022).
  • Quantum phase estimation: QFT translates eigenvalue information from the quantum register to measurable bit outcomes.
  • Quantum arithmetic: Addition, multiplication, and computation of weighted averages can be efficiently implemented in the QFT basis. For example, the Draper adder achieves modular addition with [QFTn]k,x=12ne2Ï€i xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}3 gates and requires no ancilla carries, a substantial shift from classical reversible-arithmetic strategies (Ruiz-Perez et al., 2014).
  • Simulation of circulant Hamiltonians: QFT diagonalizes circulant matrices, permitting linear-system and time-evolution simulations (Zhou et al., 2015).

6. Entanglement Properties and Classical Simulability

Although the operator QFT is maximally entangling due to the bit-reversal layer, the "core" QFT (without bit-reversal) produces operator Schmidt coefficients that decay exponentially, imposing a constant upper bound on the entanglement generated across any partition, regardless of system size (Chen et al., 2022). This property enables classical simulation of the QFT on matrix product states (MPS) of low bond-dimension in [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}4 time, permitting compressive speedups over the FFT for data that are MPS-compressible. However, this advantage does not extend to arbitrary ("maximally entangled") data.

7. Generalizations and Advanced Techniques

Generalized QFTs (GQFT) parameterized by bilinear phase matrices [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}5 have been developed for nontrivial group structures, notably extending the QFT to applications in non-Abelian hidden subgroup problems. The quantum Haar transform, required for quantum wavelet transforms and Schur duality, can be generated with explicit recursive formulas and implemented with [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}6 gates (Shao, 2017).

Alternative circuit factorizations—such as composite CZ[QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}7SWAP blocks—have yielded constant-factor improvements in two-qubit-gate resource for the QFT. These approaches preserve the [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}8 scaling but optimize practical performance in quantum pipeline settings (e.g., as subcircuits in Harrow-Hassidim-Lloyd solvers) (Romero et al., 11 Jul 2025).


Summary Table: Quantum Fourier Transform—Key Properties

Aspect Standard QFT Circuit Dynamic (Semi-Classical) QFT Fault-Tolerant AQFT (Optimized)
Two-qubit gate count [QFTn]k,x=12ne2πi xk/2n[\mathrm{QFT}_n]_{k,x} = \frac{1}{\sqrt{2^n}}e^{2\pi i\,xk/2^n}9 ∑xαx∣x⟩\sum_x\alpha_x\ket{x}0 ∑xαx∣x⟩\sum_x\alpha_x\ket{x}1 or ∑xαx∣x⟩\sum_x\alpha_x\ket{x}2
Mid-circuit measurements None ∑xαx∣x⟩\sum_x\alpha_x\ket{x}3 ∑xαx∣x⟩\sum_x\alpha_x\ket{x}4 classical; single-shot
Classical feedforward None ∑xαx∣x⟩\sum_x\alpha_x\ket{x}5 operations Optional
Max. tested qubits (2024) ∑xαx∣x⟩\sum_x\alpha_x\ket{x}6 (∑xαx∣x⟩\sum_x\alpha_x\ket{x}7 fidelity, unitary) ∑xαx∣x⟩\sum_x\alpha_x\ket{x}8 (∑xαx∣x⟩\sum_x\alpha_x\ket{x}9 fidelity, dynamic) ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}0 in asymptotic analysis
T-count (AQFT, error ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}1) ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}2 N/A ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}3
T-depth (AQFT, error ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}4) ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}5 N/A ∑k(12n∑xαxe2πi xk/2n)∣k⟩\sum_{k}(\frac{1}{\sqrt{2^n}}\sum_x\alpha_x e^{2\pi i\,xk/2^n})\ket{k}6 (minimized)

References

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Quantum Fourier Transform.