Table of Contents
Fetching ...

Quantum state preparation without coherent arithmetic

Sam McArdle, András Gilyén, Mario Berta

TL;DR

This work presents a QSVT-based framework for preparing quantum states whose amplitudes follow a known function, without relying on amplitude oracles or handcrafted coherent arithmetic. By block-encoding a sine-based template and applying a degree-d polynomial via QSVT, the method constructs |Psi_f> from a low-cost encoding with only 3 ancilla qubits and resource usage governed by the discretized L2-filling fraction. It delivers concrete, favorable complexity for smooth functions that admit polynomial or Fourier approximations, and demonstrates Gaussian and Kaiser window state preparations with favorable ancilla counts compared to black-box approaches. The approach is versatile, extensible to priors, non-smooth functions, and multivariate cases, and has practical implications for phase estimation, quantum chemistry, and financial modeling.

Abstract

We introduce a versatile method for preparing a quantum state whose amplitudes are given by some known function. Unlike existing approaches, our method does not require handcrafted reversible arithmetic circuits, or quantum table reads, to encode the function values. Instead, we use a template quantum eigenvalue transformation circuit to convert a low cost block encoding of the sine function into the desired function. Our method uses only 4 ancilla qubits (3 if the approximating polynomial has definite parity), providing order-of-magnitude qubit count reductions compared to state-of-the-art approaches, while using a similar number of gates if the function can be well represented by a polynomial or Fourier approximation. Like black-box methods, the complexity of our approach depends on the 'L2-norm filling-fraction' of the function. We demonstrate the algorithmic utility of our method, including preparing Gaussian and Kaiser window states.

Quantum state preparation without coherent arithmetic

TL;DR

This work presents a QSVT-based framework for preparing quantum states whose amplitudes follow a known function, without relying on amplitude oracles or handcrafted coherent arithmetic. By block-encoding a sine-based template and applying a degree-d polynomial via QSVT, the method constructs |Psi_f> from a low-cost encoding with only 3 ancilla qubits and resource usage governed by the discretized L2-filling fraction. It delivers concrete, favorable complexity for smooth functions that admit polynomial or Fourier approximations, and demonstrates Gaussian and Kaiser window state preparations with favorable ancilla counts compared to black-box approaches. The approach is versatile, extensible to priors, non-smooth functions, and multivariate cases, and has practical implications for phase estimation, quantum chemistry, and financial modeling.

Abstract

We introduce a versatile method for preparing a quantum state whose amplitudes are given by some known function. Unlike existing approaches, our method does not require handcrafted reversible arithmetic circuits, or quantum table reads, to encode the function values. Instead, we use a template quantum eigenvalue transformation circuit to convert a low cost block encoding of the sine function into the desired function. Our method uses only 4 ancilla qubits (3 if the approximating polynomial has definite parity), providing order-of-magnitude qubit count reductions compared to state-of-the-art approaches, while using a similar number of gates if the function can be well represented by a polynomial or Fourier approximation. Like black-box methods, the complexity of our approach depends on the 'L2-norm filling-fraction' of the function. We demonstrate the algorithmic utility of our method, including preparing Gaussian and Kaiser window states.
Paper Structure (19 sections, 16 theorems, 76 equations, 1 figure, 2 tables)

This paper contains 19 sections, 16 theorems, 76 equations, 1 figure, 2 tables.

Key Result

Theorem 1

Given a degree $d$ definite-parity function $h(y)$ such that $|h(y)|_{\mathrm{max}}^{y \in [-1,1]} \leq 1$, which approximates $f(\cdot)$ as where $\tilde{f}(y) := h(\sin(y/a))$, then we can prepare a quantum state $|\Psi_{\tilde{f}}\rangle$ that is no more than $\epsilon$-far from $|\Psi_f\rangle$ in trace distance using a quantum circuit requiring $\mathcal{O}\left( \frac{n d}{\mathcal{F}_{\ti

Figures (1)

  • Figure 1: The quantum circuit implementing QSVT-based state preparation. We define $R_y(\theta) := e^{-i\theta Y}$, $R_z(\theta) = \mathrm{Diag}(1, e^{i\theta})$. a) The circuit $U_{\mathrm{sin}}$ that block-encodes $\sum_{x} \mathrm{sin}(2x/N) \ketbra{x}{x}$ by applying a Hadamard test circuit to a directionally controlled phase gradient GidneyBlog2017 (see Lemma \ref{['Lemma:BlockEncodeSin']}). This circuit requires (n+1) $Z$ rotations, and CNOT chains that can be implemented in $\mathcal{O}(\log(n))$ depth low2018trading, and can be further optimized for fault-tolerant implementation in e.g. the surface code . b) The circuit $U_{\tilde{f}}$ that block-encodes $\sum_x \tilde{f}(\bar{x}) \ketbra{x}{x}$ by applying QSVT to $U_{\mathrm{sin}}$. The angles $\theta_i$ correspond to the pre-computed QSVT-angles for the desired polynomial. c) The (exact) amplitude-amplification circuit which block encodes $\ketbra{\Psi_{\tilde{f}}}{\bar{0}}$, including an additional qubit to adjust the amplitude (see Appendix \ref{['App:AmpAmp']}).

Theorems & Definitions (31)

  • Theorem 1
  • proof
  • Theorem 2
  • Lemma 1: Amplitude amplification
  • proof
  • Theorem 3: Exact amplitude amplification
  • proof
  • Lemma 2: see grinshpan2009AnalysisNotes
  • Lemma 3: see grinshpan2009AnalysisNotes
  • Definition 1
  • ...and 21 more