Solving 2D Black Scholes Equation via Hermitian Block Embedding and Generalised Quantum Signal Processing

Resumo Técnico: Resolução de Equações de Black–Scholes 2D via Embutimento de Bloco Hermítico e Processamento de Sinal Quântico Generalizado

1. Definição do Problema

A precificação de derivativos financeiros europeus, como opções, é fundamentalmente modelada pela equação diferencial parcial (PDE) de Black–Scholes. Para opções multi-ativos (ex: opções de cesta que dependem de dois ativos subjacentes), o problema torna-se uma PDE bidimensional. Sob o framework padrão de ausência de arbitragem, a precificação requer a evolução do payoff terminal para trás no tempo.

Ao realizar a discretização usando diferenças finitas e um esquema de passo de tempo implícito (especificamente o Euler para trás), a PDE contínua é transformada em uma sequência de problemas de álgebra linear. A tarefa computacional central é a aplicação do inverso de uma matriz de passo de tempo, $\tilde{M}^{-1}$, a um vetor que representa o valor da opção.

• Desafio: No caso bidimensional, a matriz discretizada $\tilde{M}$ é real, mas geralmente não é Hermítica nem unitária devido aos termos de drift, derivadas mistas e condições de contorno.
• Contexto Quântico: Embora algoritmos de álgebra linear quântica como o Processamento de Sinal Quântico (QSP) e o Processamento de Sinal Quântico Generalizado (GQSP) ofereçam métodos poderosos para implementar funções de matrizes (como a inversão) via transformações polinomiais, eles tipicamente exigem que o operador de entrada seja unitário ou Hermítico. Abordagens padrão para matrizes não-Hermíticas frequentemente dependem de "embutimento de bloco" (block-encoding), o que introduz um overhead significativo em termos de registradores ancila, profundidade de circuito e preparação de estado. Trabalhos anteriores [13] aplicaram com sucesso uma abordagem de Hermitian-GQSP para a equação de Black–Scholes unidimensional ao utilizar uma transformação de similaridade diagonal para converter a matriz de passo de tempo tridiagonal em uma forma Hermítica. No entanto, este método falha em duas dimensões porque a estrutura de bloco da matriz 2D (derivada da discretização por produto tensorial e derivadas mistas) impede a existência de uma transformação de similaridade diagonal que preserve a Hermiticidade.

2. Metodologia

Os autores propõem uma estratégia de embutimento de bloco Hermítico para estender o framework GQSP para a equação de Black–Scholes bidimensional sem exigir uma transformação de similaridade diagonal ou os overheads convencionais de embutimento de bloco.

A. Embutimento de Bloco Hermítico
Em vez de tentar simetrizar a matriz não-Hermítica $\tilde{M}$ diretamente, os autores a embutem em uma matriz Hermítica maior $H$ de dimensão $2N \times 2N$ (onde $N$ é a dimensão da grade discretizada):
$$H = \begin{pmatrix} 0 & \tilde{M} \\ \tilde{M}^T & 0 \end{pmatrix}$$
Como $\tilde{M}$ é real, $H$ é real simétrica e, portanto, Hermítica. O inverso desta matriz embutida é:
$$H^{-1} = \begin{pmatrix} 0 & \tilde{M}^{-T} \\ \tilde{M}^{-1} & 0 \end{pmatrix}$$
Crucialmente, o operador de passo de tempo inverso desejado $\tilde{M}^{-1}$ aparece no bloco off-diagonal inferior esquerdo de $H^{-1}$.

B. Aproximação por Polinômio Ímpar
A estrutura de $H$ impõe uma propriedade de paridade às suas potências: potências pares de $H$ são bloco-diagonais, enquanto potências ímpares são bloco-off-diagonais.

• Os autores constroem um polinômio ímpar $p(x) = \sum c_{2m+1} x^{2m+1}$ para aproximar a função inversa $f(x) = 1/(\gamma x)$ no espectro de $H$.
• Como $p(x)$ é ímpar, $p(H)$ permanece bloco-off-diagonal. Consequentemente, o bloco inferior esquerdo de $p(H)$ aproxima $\tilde{M}^{-1}$.
• Isso permite a recuperação do passo de evolução temporal para trás aplicando a transformação polinomial a um estado de entrada embutido $\begin{pmatrix} V^n \\ 0 \end{pmatrix}$ e extraindo a metade inferior do vetor resultante.

C. Escalonamento Espectral e Implementação GQSP

• Escalonamento: A matriz $H$ é reescalonada por um fator $\gamma = \max |\lambda_i(H)|$ para definir $A = H/\gamma$, garantindo que o espectro esteja dentro de $[-1, 1]$. O espectro é simétrico em relação a zero, situando-se em $[-1, -\delta] \cup [\delta, 1]$, onde $\delta$ é determinado pelo número de condição de $\tilde{M}$.
• GQSP: A matriz Hermítica reescalonada $A$ é embutida em um operador unitário $U = A + i\sqrt{I-A^2}$. O algoritmo GQSP é então usado para implementar a transformação polinomial ímpar $p(A)$ sobre $U$.
• Saída: O circuito quântico gera um estado onde o registro inferior contém a solução aproximada $V^{n+1} = \tilde{M}^{-1}V^n$.

3. Principais Contribuições

• Extensão para 2D: O artigo estende o framework Hermitian-GQSP de equações de Black–Scholes unidimensionais para bidimensionais, superando a obstrução estrutural que impede o uso de transformações de similaridade diagonal em dimensões superiores.
• Abordagem Livre de Embutimento de Bloco: Os autores demonstram um método para resolver o problema inverso para matrizes não-Hermíticas usando GQSP sem o pesado overhead das técnicas padrão de embutimento de bloco, baseando-se em um escalonamento de dilatação Hermítica específica e na estrutura de polinômio ímpar.
• Prova de Princípio: O trabalho fornece uma prova de princípio de que técnicas modernas de álgebra linear quântica podem ser aplicadas a problemas de precificação de opções multidimensionais.

4. Resultados Numéricos

Os autores realizaram simulações numéricas para opções de compra (call) de cesta europeia de dois ativos usando Python e PennyLane.

• Configuração: Duas simulações foram conduzidas com parâmetros financeiros variados (volatilidade, correlação, preço de exercício, etc.) e tamanhos de grade (grade de $32 \times 32$ pontos, embutida em $2048$ dimensões).
• Comparação: As soluções baseadas em GQSP foram comparadas com:
  1. Um método clássico de diferenças finitas de Euler para trás.
  2. Uma aproximação polinomial clássica aplicada diretamente à matriz hermitizada.
• Descobertas:
  • As soluções GQSP mostraram concordância próxima tanto com a aproximação polinomial clássica quanto com o benchmark de Euler para trás.
  • Os resultados confirmaram que o embutimento de bloco Hermítico captura com precisão a dinâmica do operador não-Hermítico original.
  • Fontes de Erro: A análise identificou duas fontes primárias de erro:
    1. Erro de Aproximação Polinomial: Dependente do gap espectral $\delta$. Na segunda simulação, um gap espectral maior resultou em menores erros de aproximação.
    2. Erro de Discretização: Devido à necessidade de usar um único passo de tempo grande (conforme implementado como um solver de passo único), o erro de discretização foi significativo (até \approx \0,09$ na primeira simulação e $\approx 60\%$ maior na segunda). Este erro é inerente ao benchmark clássico também.
  • Ruído: A simulação quântica exibiu um piso de ruído (observado em torno de $10^{-2}$ na segunda simulação), que não estava presente na aproximação polinomial clássica.

5. Significância e Alegações

O artigo alega que o método proposto demonstra a viabilidade de combinar embutimentos de bloco Hermíticos com GQSP para resolver problemas de Black–Scholes multidimensionais.

• Escopo Modesto: Os autores declaram explicitamente que a implementação atual é um "solver de inverso de passo único". Eles não reivindicam uma vantagem quântica total ainda, observando que o método é limitado pela necessidade de polinômios de grau mais alto quando o gap espectral é pequeno e pelos erros de discretização de passos de tempo grandes.
• Potencial Futuro: A significância reside em fornecer um caminho para algoritmos quânticos para derivativos multi-ativos onde computadores clássicos enfrentam a "maldição da dimensionalidade". Os autores sugerem que a verdadeira vantagem quântica pode emergir em 3+ dimensões, onde os métodos clássicos de diferenças finitas tornam-se computacionalmente intratáveis.
• Limitações: O artigo reconhece que trabalhos futuros são necessários para desenvolver esquemas de evolução temporal multi-passo eficientes, reduzir o piso de ruído e minimizar os erros de discretização através de melhores estratégias de grade ou implementações de GQSP multi-passo.

Em resumo, o artigo estabelece uma base teórica e numérica para aplicar GQSP à precificação de opções 2D ao introduzir um embutimento de bloco Hermítico que evita a necessidade de transformações de similaridade diagonal, oferecendo uma rota viável para algoritmos quânticos para derivativos financeiros de alta dimensão.