Simulating Electron Transfer on Noisy Quantum Computers

Resumo Técnico: Simulando a Transferência de Elétrons em Computadores Quânticos Ruidosos

Definição do Problema
O projeto de materiais energéticos de próxima geração, como os destinados a baterias e fotovoltaicos orgânicos, requer a simulação precisa de processos de transferência de elétrons (ET) fora do equilíbrio. Teorias cinéticas clássicas frequentemente falham em capturar taxas de transferência em regimes caracterizados por altos sobrepotenciais ou camadas de solvatação complexas, particularmente onde a coerência eletrônica-vibracional (vibônica) de longa duração desempenha um papel crítico. Embora modelos simples de spin-boson tenham sido realizados em hardware quântico, simular redes eletrônicas estendidas com ambientes vibracionais locais permanece um desafio fundamental. Isso é especialmente verdadeiro para dinâmicas não-markovianas impulsionadas por tempos de vida vibracionais finitos, onde métodos clássicos como a Equação de Hierárquica de Movimento (HEOM) ou o Tensor-TEDOPA enfrentam severas limitações de escalonamento conforme o tamanho do sistema aumenta.

Metodologia
Os autores apresentam uma estrutura digital-analógica para simular sistemas quânticos abertos governados por Hamiltonianos com acoplamento linear-vibônico (LVC) e ambientes vibracionais estruturados. A abordagem utiliza o formalismo de pseudomódulos, onde a dinâmica não-markoviana de uma rede eletrônica acoplada a um banho estruturado é projetada ao acoplar cada sítio eletrônico a um oscilador harmônico local e amortecido (pseudomódulo).

Principais componentes metodológicos incluem:

• Exploração de Recursos de Hardware: Em vez de lutar contra o ruído do hardware, a estrutura explora a dissipação intrínseca dos qubits supercondutores para emular o relaxamento vibracional. A taxa de amortecimento efetiva ($\Gamma_{QC}$) é determinada pelos tempos de coerência $T_1$ e $T_2$ dos qubits que codificam os osciladores.
• Mapeamento de Qubits: Uma cadeia eletrônica de $N$ sítios é mapeada para $N$ qubits, e $N$ osciladores locais são mapeados para $N$ qubits adicionais usando uma codificação bóson-para-spin (truncada para dois níveis de energia por oscilador).
• Mitigação de Erros: Um esquema de mitigação de erro específico do modelo é empregado para filtrar fontes de ruído incompatíveis com o sistema aberto alvo. Isso envolve o descarte de disparos de medição (pós-processamento) que violam a conservação de partículas no subsistema eletrônico ou excedem um limite de excitações vibracionais. Isso filtra o ruído de despolarização enquanto retém o amortecimento de amplitude intrínseco do hardware.
• Trotterização: A evolução temporal é simulada usando circuitos trotterizados com um passo de tempo $\Delta t = 4$ fs, equilibrando a profundidade do circuito contra o erro de Trotter e a resolução espectral do Hamiltoniano.

Resultos Principais
A estratégia foi validada no processador quântico supercondutor IBM Aachen ao simular um modelo microscópico de ET de um doador para uma cadeia de aceitadores.

• Resolução de Ressonâncias: A equipe conseguiu resolver tanto os espectros de transferência eletrônica quanto os vibrônicos para uma cadeia de $N=5$ sítios. Eles demonstraram que o hardware quântico pode distinguir entre ressonância puramente eletrônica e ressonância vibrônica (mediada pela interação coerente elétron-oscilador), sendo que esta última aparece apenas após vários ciclos vibracionais.
• Escalonamento: A abordagem foi escalada para 10 sítios eletrônicos (20 qubits no total), representando uma escala sem precedentes para dinâmica química em computadores quânticos. As simulações reproduziram dinâmicas não-markovianas e o acúmulo de população nos sítios aceitadores.
• Amortecimento Efetivo: Os experimentos alcançaram um tempo de vida vibracional efetivo ($1/\Gamma_{QC}$) entre 50 e 150 fs. Embora seja mais curto que os tempos de vida alvo de modos de estiramento C=C intramolecular ($>500$ fs), aproxima-se do regime onde métodos clássicos perturbativos tornam-se imprecisos.
• Caracterização de Ruído: A análise revelou que a taxa de amortecimento efetiva é limitada principalmente pelo qubit com o menor tempo de coerência no circuito. O número de disparos válidos restantes após a mitigação de erro escala exponencialmente com o número de qubits, mas polinomialmente com os passos de tempo dentro da janela de coerência.

Significância e Alegações
O artigo posiciona este trabalho como um benchmark portátil e orientado à aplicação para simular estados emaranhados de longa duração em dispositivos NISQ (Noisy Intermediate-Scale Quantum).

• Benchmarking de Emaranhamento: O mecanismo de ET vibrônico depende de estados não-separáveis onde os sítios eletrônicos e os osciladores estão emaranhados. Simular com sucesso este processo até 10 sítios serve como uma métrica para a capacidade de um processador quântico de sustentar e escalar bipartições emaranhadas.
• Eficiência de Recursos: Ao tratar a dissipação do hardware como um recurso em vez de um defeito, o método evita as codificações bóson-para-qubit intensivas em recursos exigidas por outras abordagens.
• Escalabilidade: Os autores argumentam que, enquanto os métodos clássicos de redes de tensores (como o DAMPF) escalam mal com o tamanho do sistema devido às restrições de dimensão de ligação em regimes não-perturbativos, seu algoritmo quântico escala linearmente com o número de qubits e mantém uma profundidade de circuito constante por passo de Trotter (para topologias específicas).
• Limitações: Os autores notam modestamente que o gargalo atual é a disponibilidade limitada de qubits de alta fidelidade e longa coerência. Eles reconhecem que simular modelos alvo com taxas de amortecimento ainda menores (tempos de vida mais longos) requer melhorias de hardware além das capacidades atuais, mas os resultados atuais demonstram a viabilidade da abordagem para capturar características essenciais não-markovianas.