{"title":"Generalized Quantum Signal Processing","authors":"Danial Motlagh, Nathan Wiebe","doi":"10.1103/prxquantum.5.020368","DOIUrl":null,"url":null,"abstract":"Quantum signal processing (QSP) and quantum singular value transformation (QSVT) currently stand as the most efficient techniques for implementing functions of block-encoded matrices, a central task that lies at the heart of most prominent quantum algorithms. However, current QSP approaches face several challenges, such as the restrictions imposed on the family of achievable polynomials and the difficulty of calculating the required phase angles for specific transformations. In this paper, we present a generalized quantum signal processing (GQSP) approach, employing general SU(2) rotations as our signal-processing operators, rather than relying solely on rotations in a single basis. Our approach lifts all practical restrictions on the family of achievable transformations, with the sole remaining condition being that <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mo stretchy=\"false\">|</mo></mrow><mi>P</mi><mrow><mo stretchy=\"false\">|</mo></mrow><mo>≤</mo><mn>1</mn></math>, a restriction necessary due to the unitary nature of quantum computation. Furthermore, GQSP provides a straightforward recursive formula for determining the rotation angles needed to construct the polynomials in cases where <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>P</mi></math> and <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>Q</mi></math> are known. In cases where only <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>P</mi></math> is known, we provide an efficient optimization algorithm capable of identifying in under a minute of GPU time, a corresponding <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>Q</mi></math> for polynomials of degree on the order of <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><msup><mn>10</mn><mn>7</mn></msup></math>. We further illustrate GQSP simplifies QSP-based strategies for Hamiltonian simulation, offer an optimal solution to the <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>ϵ</mi></math>-approximate fractional query problem that requires <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi mathvariant=\"script\">O</mi></mrow><mrow><mo>(</mo><mo stretchy=\"false\">(</mo><mn>1</mn><mo>/</mo><mi>δ</mi><mo stretchy=\"false\">)</mo><mo>+</mo><mi>log</mi><mo></mo><mo stretchy=\"false\">(</mo><mn>1</mn><mo>/</mo><mi>ϵ</mi><mo stretchy=\"false\">)</mo><mo>)</mo></mrow></math> queries to perform where <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi mathvariant=\"script\">O</mi></mrow><mo stretchy=\"false\">(</mo><mn>1</mn><mo>/</mo><mi>δ</mi><mo stretchy=\"false\">)</mo></math> is a proved lower bound, and introduces novel approaches for implementing bosonic operators. Moreover, we propose a novel framework for the implementation of normal matrices, demonstrating its applicability through synthesis of diagonal matrices, as well as the development of a new algorithm for convolution through synthesis of circulant matrices using only <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi mathvariant=\"script\">O</mi></mrow><mo stretchy=\"false\">(</mo><mi>d</mi><mi>log</mi><mo></mo><mi>N</mi><mo>+</mo><msup><mi>log</mi><mn>2</mn></msup><mo></mo><mi>N</mi><mo stretchy=\"false\">)</mo></math> 1 and 2-qubit gates for a filter of lengths <math display=\"inline\" overflow=\"scroll\" xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>d</mi></math>.","PeriodicalId":501296,"journal":{"name":"PRX Quantum","volume":"364 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"PRX Quantum","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1103/prxquantum.5.020368","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Quantum signal processing (QSP) and quantum singular value transformation (QSVT) currently stand as the most efficient techniques for implementing functions of block-encoded matrices, a central task that lies at the heart of most prominent quantum algorithms. However, current QSP approaches face several challenges, such as the restrictions imposed on the family of achievable polynomials and the difficulty of calculating the required phase angles for specific transformations. In this paper, we present a generalized quantum signal processing (GQSP) approach, employing general SU(2) rotations as our signal-processing operators, rather than relying solely on rotations in a single basis. Our approach lifts all practical restrictions on the family of achievable transformations, with the sole remaining condition being that , a restriction necessary due to the unitary nature of quantum computation. Furthermore, GQSP provides a straightforward recursive formula for determining the rotation angles needed to construct the polynomials in cases where and are known. In cases where only is known, we provide an efficient optimization algorithm capable of identifying in under a minute of GPU time, a corresponding for polynomials of degree on the order of . We further illustrate GQSP simplifies QSP-based strategies for Hamiltonian simulation, offer an optimal solution to the -approximate fractional query problem that requires queries to perform where is a proved lower bound, and introduces novel approaches for implementing bosonic operators. Moreover, we propose a novel framework for the implementation of normal matrices, demonstrating its applicability through synthesis of diagonal matrices, as well as the development of a new algorithm for convolution through synthesis of circulant matrices using only 1 and 2-qubit gates for a filter of lengths .