McLachlan-projected reduced dynamics for ill-posed Schrödingerized backward diffusion

Resumo Técnico: Dinâmica Reduzida Projetada por McLachlan para a Schrödingerização Mal-Posta da Difusão Reversa

Enunciado do Problema
O artigo aborda o desafio numérico de resolver equações de difusão reversa, que são problemas de evolução mal-postos prototípicos. Nestes sistemas, as frequências espaciais altas crescem exponencialmente ao longo do tempo, fazendo com que esquemas de marcha temporal baseados em malha (como Crank–Nicolson), sem regularização explícita, sejam rapidamente sobrecarregados por ruído. Embora o framework de "Schrödingerização" tenha sido estabelecido para incorporar dinâmicas lineares não unitárias (incluindo problemas mal-postos) em dinâmicas hermitianas em um espaço estendido, uma lacuna crítica permanece: não está claro se o vetor de amplitude completo elevado deve ser propagado a cada passo de tempo, ou se uma representação reduzida de dimensão fixa é suficiente. Além disso, numéricas comparativas existentes frequentemente falham em separar os efeitos do levantamento, da truncagem do subespaço e da estimativa estatística dos elementos da matriz.

Metodologia
Os autores propõem uma abordagem de regularização estruturada aplicando o princípio variacional de McLachlan ao sistema elevado Schrödingerizado. A metodologia procede em três camadas:

1. Realização Elevada: A equação de difusão reversa semidiscreta $\dot{u}_h = A_h u_h$ (onde $A_h$ possui autovalores não negativos) é mapeada para um sistema hermitiano em um espaço estendido de dimensão $M$. Isso é formalizado via Hipótese 1, que postula a existência de um gerador hermitiano $H$, uma injeção $J$ e um mapa de recuperação $R$ tal que o estado físico é aproximado por $R e^{-itH} J u_0$ até um erro de levantamento $\varepsilon_{\text{lift}}$.
2. Projeção de McLachlan: Em vez de propagar o estado completo de dimensão $M$, os autores constroem uma dinâmica reduzida em um quadro fixo $V$ de $m$ instantâneos elevados. Ao minimizar a norma do resíduo da equação de Schrödinger dentro do span de $V$, eles derivam uma equação de evolução reduzida para o vetor de coeficientes $c(t)$:
   $$iS \dot{c} = H_K c$$
   onde $S = V^\dagger V$ é a matriz de Gram (sobreposição) e $H_K = V^\dagger H V$ é o Hamiltoniano projetado.
3. Decomposição de Erro e Estimativa Estatística: O artigo decompõe rigorosamente o erro total em componentes de discretização espacial, levantamento/recuperação, projeção e amostragem. Analisa especificamente o impacto de estimar o lápis $(S, H_K)$ via medições quânticas de tiro finito (por exemplo, testes de Hadamard), derivando limites de perturbação que dependem do número de condição $\kappa(S)$.

Principais Contribuições
O artigo faz as seguintes contribuições teóricas e metodológicas específicas:

• Unicidade e Conservação: Prova a unicidade da solução do fluxo reduzido e estabelece que o fluxo hermitiano projetado conserva a norma de Gram ($d/dt(c^\dagger S c) = 0$).
• Limites de Erro: Os autores derivam um limite de lacuna elevado–reduzido (Teorema 2) expresso em termos do defeito de projeção $\eta_V(t) = \|(I-P_V)H\Psi_m(t)\|_2$. Este defeito separa a propagação elevada completa da trajetória reduzida, fornecendo um orçamento de erro quantificável.
• Análise de Perturbação: A Proposição 1 fornece limites sobre a perturbação do gerador reduzido e dos coeficientes propagados quando as matrizes de sobreposição e Hamiltoniano são estimadas com ruído estatístico. Isso destaca o papel crítico do condicionamento da matriz de sobreposição $\kappa(S)$ na amplificação do ruído de tiro.
• Baselines Clássicas Justas: O artigo constrói explicitamente baselines clássicas justas (filtragem passa-baixa espectral ou regularização de Tikhonov) onde a frequência de corte ou a força da crista são ajustadas ao conteúdo de informação (dimensão $m$) do quadro reduzido escolhido. Isso isola o desempenho do pipeline reduzido Schrödingerizado da instabilidade de solucionadores não regularizados.

Resultados
Experimentos numéricos em um problema de difusão reversa 1D ($N_x=7$, $T=0.3$) usando Qiskit Aer com estimativa de tiro finito ($10^4$ tiros) demonstram:

• Estabilidade vs. Métodos Não Regularizados: Um esquema Crank–Nicolson não regularizado exibe um erro relativo discreto $L_2$ de $\approx 41,99\%$ devido à explosão de ruído. Em contraste, o levantamento Schrödingerizado (dimensão completa) limita o erro a $\approx 19,37\%$.
• Desempenho da Dinâmica Reduzida: A projeção de McLachlan com $m=4$ com matrizes clássicas (exatas) produz um erro de $\approx 20,02\%$. Ao usar matrizes amostradas pelo Aer (simulando ruído de hardware quântico), o erro é de $\approx 19,82\%$, com um pico transitório de $\approx 22,03\%$.
• Mascaramento de Ruído: Apesar de uma grande discrepância na norma de Frobenius entre o Hamiltoniano amostrado e o exato ($\| \Delta H \|_F / \| H_K \|_F \approx 7,08$) e uma matriz de sobreposição altamente mal-condicionada ($\kappa(S) \approx 1,7 \times 10^4$), a trajetória física recuperada não exibe explosão catastrófica. Os resultados confirmam que a camada de projeção atua como um regularizador estruturado, e o erro permanece na faixa de dezenas de por cento, consistente com os limites teóricos derivados do defeito de projeção e das perturbações de amostragem.

Significado e Afirmações
Os autores posicionam este trabalho como um avanço metodológico, em vez de uma afirmação de resolver completamente a má-postura. Eles afirmam explicitamente:

• Sem Cura para a Má-Postura: A Schrödingerização não remove a natureza mal-posta da variável física recuperada, nem a projeção supera o design ótimo de Tikhonov.
• Regularização Interpretável: A principal contribuição é enquadrar o erro de projeção como um "botão interpretável" análogo a uma largura de banda passa-baixa. Isso permite comparações "de igual para igual" entre dinâmicas reduzidas inspiradas em quântica e métodos regularizados clássicos.
• Separação de Erros: O framework isola com sucesso as contribuições do levantamento, da truncagem do subespaço e da estimativa estatística, prevenindo a confusão dessas fontes de erro distintas que frequentemente ocorre na literatura atual.
• Validação da Estrutura: Os resultados ilustram que, mesmo com ruído estatístico significativo na estimativa do gerador, a natureza estruturada da projeção de McLachlan previne a amplificação exponencial de ruído que caracteriza a difusão reversa não regularizada.