Real-time Krylov Diagonalisation for Open Quantum Systems

Resumo Técnico: Diagonalização de Krylov em Tempo Real para Sistemas Quânticos Abertos

Enunciado do Problema
O artigo aborda o desafio de simular sistemas quânticos abertos descritos por equações mestras de Lindblad, focando especificamente na estimativa do gap de Liouvilliano. Embora a Diagonalização de Subespaço Quântico (QSD) tenha surgido como um algoritmo líder não tolerante a falhas para sistemas fechados (Hamiltonianos), sua aplicação a sistemas dissipativos não é trivial. Sistemas abertos são governados por superoperadores de Liouvilliano não hermitianos e não normais, que introduzem diferenças estruturais como autovetores bi-ortogonais, autovalores complexos e dinâmicas não unitárias. Essas características complicam os métodos padrão de subespaço de Krylov, particularmente no que diz respeito à estabilidade numérica (condicionamento ruim) e à interpretação de filtros espectrais. Os autores visam adaptar a QSD em tempo real para estimar o gap de Liouvilliano, um parâmetro crítico que define o tempo de relaxação e a presença de modos lentos em transições de fase dissipativas.

Metodologia
Os autores propõem uma modificação da Diagonalização de Subespaço Quântico (QSD) em tempo real para lidar com dinâmicas de Liouvilliano. A metodologia central envolve:

1. Construção do Subespaço de Krylov: Em vez de evoluir um estado sob um Hamiltoniano unitário $H$, o sistema evolui sob o superoperador de Liouvilliano $\mathcal{L}$. Um subespaço de Krylov é gerado aplicando potências do operador de evolução temporal $e^{\mathcal{L}\tau}$ a um estado inicial $|\rho_0)$:
   $$\text{span}\{|\rho_0), e^{\mathcal{L}\tau}|\rho_0), \dots, e^{\mathcal{L}(D-1)\tau}|\rho_0)\}$$
2. Problema de Autovalores Generalizado (GEVP): Como circuitos quânticos não podem realizar reortogonalização em tempo real, o método resolve um GEVP $\bar{H}v = \lambda Sv$. Para Liouvillianos, a matriz de sobreposição $S$ e a matriz Hamiltoniana $\bar{H}$ são construídas usando autovetores esquerdo e direito bi-ortogonais.
3. Interpretação do Filtro Espectral: O método é enquadrado como a otimização de um filtro espectral. Diferentemente de sistemas Hamiltonianos, onde os autoestados são reais e a evolução temporal é oscilatória (exigindo amostragem de Nyquist estrita), as dinâmicas de Liouvilliano são dissipativas. Modos de decaimento rápido espiralam para dentro no plano complexo, permitindo tempos de evolução $\tau$ mais longos sem artefatos de aliasing. Isso permite o isolamento de modos lentos (autovalores com partes reais pequenas) mesmo com uma dimensão de Krylov reduzida $D$.
4. Tratamento da Não Normalidade: Os autores reconhecem que operadores não normais levam a alta sensibilidade a perturbações (teorema de Bauer-Ficke). A matriz de sobreposição $S$ torna-se exponencialmente mal condicionada à medida que $D$ aumenta. Para mitigar isso, o método confia na natureza dissipativa do sistema para suprimir modos rápidos, permitindo $\tau$ maiores e $D$ menores, reduzindo assim a instabilidade numérica associada a grandes dimensões de Krylov.
5. Seleção do Estado Inicial: Crucialmente, o estado inicial deve ter uma projeção não nula no setor de traço zero (coerências ou desequilíbrios de população) para excitar os modos lentos. O setor de traço 1 (estado estacionário) não contribui para o sinal de estimativa do gap neste quadro.

Contribuições Principais

• Extensão da QSD para Liouvillianos: O artigo fornece um quadro teórico para aplicar QSD em tempo real a superoperadores de Liouvilliano não hermitianos e não normais na forma de Lindblad.
• Estratégia de Estimativa de Gap: Demonstra como estimar o gap de Liouvilliano ($\Delta = -\max \Re(\lambda_k)$) explorando a natureza dissipativa do sistema para filtrar modos de decaimento rápido, contornando efetivamente algumas restrições de Nyquist que se aplicam à evolução unitária.
• Análise de Estabilidade Numérica: O trabalho destaca o compromisso entre o tempo de evolução $\tau$ e a dimensão de Krylov $D$. Argumenta que aumentar $\tau$ (em vez de $D$) é uma estratégia viável para sistemas abertos para melhorar a resolução do gap enquanto mantém a estabilidade numérica, já que os modos rápidos decaem naturalmente.
• Aplicação a Qubits Cat de Kerr: O método é aplicado a um ressonador de Kerr supercondutor acionado por dois fótons (Qubit Cat de Kerr), um sistema relevante para arquiteturas de qubits protegidos.

Resultados
Os autores aplicaram seu método a um modelo simplificado de um Qubit Cat de Kerr com perda de fóton único.

• Parâmetros de Simulação: As simulações foram realizadas com um máximo de 30 fótons e uma dimensão de Krylov $D=20$, resultando em uma matriz de Krylov de $900 \times 20$.
• Reconstrução do Gap: O método reconstruiu com sucesso o gap de Liouvilliano através de diferentes forças de acionamento. Os resultados mostraram que aumentar o tempo de evolução $\tau$ melhorou a precisão da estimativa do gap, particularmente em regimes onde o comutador de Liouvilliano (uma medida de não normalidade) era alto.
• Identificação de Modo Lento: As simulações confirmaram o surgimento de um modo lento correspondente ao decaimento de coerências e desequilíbrios de população. No regime de "ruído enviesado", o gap encolhe exponencialmente com a amplitude do estado coerente ($\propto e^{-\alpha^2}$), indicando a formação de um subespaço computacional protegido.
• Representação de Wigner: Os autoestados reconstruídos da matriz de Krylov foram usados para visualizar a representação de Wigner do modo lento, mostrando a transição de um regime amortecido para um regime com dois lóbulos distintos (quebra de simetria) e de volta a um regime comprimido em alta acionamento.

Significado e Alegações
O artigo afirma demonstrar a viabilidade de usar QSD em tempo real para estudar sistemas quânticos abertos, especificamente para estimar o gap de Liouvilliano no contexto de transições de fase dissipativas. Os autores posicionam este trabalho como uma adaptação necessária de algoritmos da era NISQ para o projeto e simulação de hardware quântico, que inerentemente envolve dissipação.

O significado é enquadrado em torno do potencial de usar computadores quânticos para simular modos lentos de sistemas quânticos que são computacionalmente caros de simular classicamente, particularmente no contexto de qubits protegidos como o Cat de Kerr. Os autores observam modestamente que, embora Qubits Cat de Kerr individuais sejam simuláveis classicamente, redes de tais qubits podem não o ser, sugerindo a utilidade do método para sistemas maiores e mais complexos. O trabalho é apresentado como um passo em direção à "descoberta automatizada de hardware quântico", fornecendo ferramentas para caracterizar as propriedades de ruído e a estabilidade de arquiteturas de qubits propostas. O artigo não afirma ter resolvido o problema geral da simulação de operadores não normais, mas oferece uma abordagem específica e modificada para estimativa de gap em sistemas dissipativos.