Estimativa da fase quântica e algoritmo HHL - conhecimento dos valores próprios necessários?

10

O algoritmo de estimativa de fase do quantum (QPE) calcula uma aproximação do valor próprio associado a um dado vector próprio de um portão quântico U .

Formalmente, vamos |ψ ser um autovetor de U , QPE nos permite encontrar |θ~ , o melhor m bits de aproximação 2mθ tais que θ[0,1) e

U|ψ=e2πiθ|ψ.

O algoritmo HHL ( artigo original ) assume como entrada uma matriz A que satisfaz

eiAt is unitary 
e um estado quântico |b e calcula |x que codifica a solução do sistema linear Ax=b .

Observação : Cada matriz transposta conjugada statisfy a condição em A .

Para fazer isso, o algoritmo HHL usa o QPE na porta quântica representada por U=eiAt . Graças aos resultados de álgebra linear, sabemos que, se {λj}j são os valores próprios de A , em seguida, {eiλjt}j são os valores próprios de U . Esse resultado também é afirmado nos algoritmos de sistemas lineares quânticos: um primer (Dervovic, Herbster, Mountney, Severini, Usher & Wossnig, 2018) (página 29, entre as equações 68 e 69).

Com a ajuda de QPE, o primeiro passo do algoritmo HLL vai tentar estimar , tal que e i 2 ¸ q = e i λ j t . Isso nos leva à equação 2 π θ = λ j t + 2 k π ,θ[0,1)ei2πθ=eiλjt ie θ = λ j t

2πθ=λjt+2kπ,kZ, θ[0,1)
Analisando um pouco as implicações das condições k Z e θ [ 0 , 1 ) , concluí que se λ j t
θ=λjt2π+k,kZ, θ[0,1)
kZθ[0,1)(ou seja,k0), o algoritmo de estimativa de fase falha em prever o valor próprio correto.λjt2π[0,1)k0

Porém, como pode ser qualquer matriz eremita, podemos escolher livremente seus autovalores e, particularmente, poderíamos escolher autovalores arbitrariamente grandes para A, de modo que o QPE falhe ( λ j tAA).λjt2π[0,1)

No Projeto de Circuitos Quânticos para Resolução de Sistemas Lineares de Equações (Cao, Daskin, Frankel & Kais, 2012), eles resolvem esse problema simulando , sabendo que os autovalores deAsão{1,2,4,8}. Elesnormalizadosda matriz (e os seus valores próprios) para evitar o caso ondeλjteiAt16A{1,2,4,8}.λjt2π[0,1)

Por outro lado, parece que o parâmetro poderia ser usado para fazer essa normalização.t

Aλjt2π[0,1)

Nelimee
fonte
t=1k=2t=10<(λ/2π)+2<14π<λ<2πλ=3πλ/2π0
λθ[0,1)λ<0λλ2kπk=λ2πλ2π

Respostas:

6

Atλt2πANQ

min i a i i -j i | a i j | -NQ. a i j A

maxiaii+ji|aij|NQ,
miniaiiji|aij|NQ.
aijA

Dentro dos valores de , , se você estiver preocupado com uma matriz grande (por exemplo, qubits), enquanto a soma da linha pode ser fácil de calcular (porque não há muitas entradas), o máximo de todas as linhas pode demorar tempo (porque há linhas), haverá várias maneiras de obter boas aproximações (por exemplo, amostragem ou uso do conhecimento da estrutura do problema). Na pior das hipóteses, você provavelmente pode usar a pesquisa de Grover para acelerar um pouco.Q n 2 nNQn2n

DaftWullie
fonte
11
Grover não é uma melhoria: mesmo que possamos usar o algoritmo, ainda precisaremos de consultas , que destroem a melhoria exponencial do HHL em relação aos métodos clássicos e o substituem por uma aceleração quadrática. Portanto, a única esperança que resta é amostrar (introduzir outra fonte de erros) ou rezar e esperar que o problema nos permita estimar os limites superior / inferior. Parece uma grande falha do algoritmo para mim. O(N)
Nelimee 04/07
2
Claro, eu só queria dizer que Grover oferece uma aceleração de raiz quadrada em comparação com a maneira ingênua de obter o máximo. Claro que isso tem um impacto ruim no tempo de execução geral.
DaftWullie