⚛️ quantum physics

High-order splitting of non-unitary operators on quantum computers

Este artigo apresenta um método de divisão de operadores de alta ordem com coeficientes complexos para simular dinâmicas dissipativas em processadores quânticos, demonstrando experimentalmente em um sistema de íons presos que essa abordagem supera métodos de baixa ordem em precisão, mesmo na presença de ruído.

Autores originais: Peter Brearley, Philipp Pfeffer

Publicado 2026-02-17

📖 4 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Peter Brearley, Philipp Pfeffer

Artigo original sob licença CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

Imagine que você tem um computador quântico. A grande vantagem dele é que ele é excelente para simular coisas que seguem as regras da mecânica quântica "perfeita", onde nada se perde e tudo é reversível, como um filme sendo passado para trás e para frente sem perder qualidade. Isso é chamado de dinâmica unitária.

No entanto, o mundo real é "sujo". Coisas esfriam, atrito para carros, ondas no mar perdem energia e o café esfria na xícara. Esses processos são chamados de dissipativos (eles perdem energia) e são irreversíveis (você não consegue fazer o café voltar a ficar quente sozinho). Simular isso em um computador quântico é muito difícil, porque a máquina "odeia" perder informação.

O Problema: O "Passo Atrás" Perigoso

Para simular coisas complexas, os cientistas usam uma técnica chamada "divisão de operadores". É como tentar desenhar uma curva complexa usando apenas retas. Você divide o tempo em pequenos passos.

Para fazer isso com alta precisão (alta ordem), a matemática tradicional exige que, em alguns momentos, você dê "passos para trás" no tempo (coeficientes negativos).

Para o mundo perfeito (unitário): Dar um passo para trás é fácil. É como tocar um vídeo ao contrário.
Para o mundo real (dissipativo): Dar um passo para trás é um pesadelo. Se você tentar simular o atrito "ao contrário", o computador tenta criar energia do nada, o que faz a simulação explodir numericamente (ficar instável e errada).

Por isso, até agora, os computadores quânticos só conseguiam simular esses processos "sujos" com baixa precisão (passos grandes e grosseiros).

A Solução: O "Truque do Tempo Imaginário"

Os autores deste artigo (Peter Brearley e Philipp Pfeffer) encontraram uma maneira genial de contornar isso. Eles usaram uma matemática especial com coeficientes complexos.

Pense nisso como uma viagem de duas dimensões:

Tempo Real: Onde a dissipação (perda de energia) acontece de verdade.
Tempo Imaginário: Um "mundo paralelo" matemático onde a dissipação se transforma em algo que o computador quântico adora: uma rotação unitária (como girar um cubo mágico).

A ideia é:

Eles dividem o problema em pequenas etapas.
Nas etapas de "perda de energia" (dissipação), eles usam números complexos. A parte real do número garante que a energia diminua (o que é o que queremos), mas a parte imaginária transforma essa perda em uma rotação segura que o computador pode fazer.
Nas etapas de "movimento" (como ondas), eles usam números reais normais.

É como se, para simular um carro freando, o computador dissesse: "Ok, vou girar o volante para a esquerda (tempo imaginário) e depois para a direita (tempo real) de uma forma tão específica que, no final, o carro parece que freou perfeitamente, sem nunca ter tentado acelerar para trás de forma perigosa."

O Experimento: Ondas que Morrem

Para provar que isso funciona, eles simularam um problema clássico: ondas mecânicas com atrito (como uma corda de violão que para de vibrar com o tempo).

A Simulação: Eles criaram um circuito quântico que divide o tempo em passos muito precisos (ordem 4 e ordem 6).
O Hardware: Eles rodaram isso em um computador quântico real da IonQ (com íons presos).
O Resultado:
- O método de baixa ordem (passos simples) foi impreciso.
- O método de alta ordem (passos complexos) foi muito mais preciso, mesmo que o circuito fosse maior e mais propenso a erros de ruído do computador.
- Curiosamente, o método de ordem 4 foi o "ponto ideal": mais preciso que os métodos simples, mas ainda estável o suficiente para funcionar no hardware atual.

Por que isso é importante?

Antes, para simular coisas do mundo real (como clima, fluidos, reações químicas) em um computador quântico, tínhamos que usar métodos "toscos" e imprecisos.

Agora, temos uma ferramenta que permite usar métodos sophisticados e precisos (alta ordem) mesmo para processos que perdem energia. Isso significa que, em breve, poderemos usar computadores quânticos para resolver problemas reais de engenharia e ciência que envolvem atrito, calor e difusão, com uma precisão que antes era impossível.

Resumo da Ópera:
Os autores criaram uma "ponte matemática" que permite ao computador quântico simular o mundo "sujo" e irreversível da vida real, transformando a perda de energia em movimentos seguros e controlados, permitindo simulações muito mais precisas do que nunca antes foi possível.

1. O Problema

A maioria dos processos físicos relevantes na ciência computacional, como fricção, viscosidade, difusão e fluxos de fronteira, é dissipativa e irreversível. Essas dinâmicas são descritas por operadores não unitários.

Desafio Quântico: Os computadores quânticos são nativamente adequados para simular dinâmicas unitárias (evolução temporal reversível de Hamiltonianos locais). Simular dinâmicas dissipativas é desafiador porque a evolução reversa (necessária para métodos de alta ordem) é numericamente instável para operadores dissipativos e não pode ser facilmente codificada em blocos unitários sem redimensionamento.
Limitação Atual: As fórmulas de produto de alta ordem (acima de 2ª ordem), que garantem alta precisão na simulação de Hamiltonianos, exigem coeficientes negativos (passos de tempo negativos). Para dinâmicas dissipativas, isso leva a instabilidades numéricas. Até este trabalho, a simulação quântica de dissipação estava limitada a métodos de baixa ordem (1ª e 2ª ordem), que exigem circuitos profundos para atingir tolerâncias de erro razoáveis.

2. Metodologia

Os autores propõem um algoritmo quântico estável para integrar dinâmicas não unitárias utilizando fórmulas de produto com coeficientes complexos.

Decomposição do Gerador: O gerador da evolução $M$ $M$ (onde $d\vec{\psi}/dt = M\vec{\psi}$ $d ψ / d t = M ψ$ ) é decomposto em componentes Hermitianos ( $H_1$ $H_{1}$ ) e anti-Hermitianos ( $iH_2$ $i H_{2}$ ), de modo que $M = H_1 + iH_2$ $M = H_{1} + i H_{2}$ .
- $H_1$ representa a parte dissipativa (autovalores reais e não positivos).
- $H_2$ representa a parte unitária (autovalores imaginários).
Fórmulas de Produto Complexas: Em vez de usar coeficientes reais negativos, o método utiliza coeficientes complexos $a_j$ $a_{j}$ para a evolução dissipativa e coeficientes reais positivos $b_j$ $b_{j}$ para a evolução unitária.
- Condição de Estabilidade: Os coeficientes complexos $a_j$ são escolhidos de modo que sua parte real seja positiva ( $\Re(a_j) > 0$ ). Isso garante que a evolução no "tempo imaginário" (dissipativa) permaneça contrativa (estável), enquanto a parte imaginária de $a_j$ introduz uma evolução unitária adicional.
- Estrutura do Circuito: A dinâmica total é decomposta em uma sequência de evoluções Hamiltonianas simples em tempo real e tempo imaginário.
  - Subetapas unitárias: $e^{i b_j H_2 \Delta t}$ (coeficientes reais positivos).
  - Subetapas dissipativas: $e^{a_j H_1 \Delta t}$ (coeficientes complexos com parte real positiva).
Ordens Alcançadas: O trabalho implementa e valida fórmulas de ordem 4 (baseada em Castella et al.) e ordem 6 (baseada em Bernier et al.), que satisfazem as restrições de estabilidade.

3. Contribuições Principais

Algoritmo Estável de Alta Ordem: Demonstração de que é possível realizar divisão de operadores de ordem superior (4 e 6) para dinâmicas dissipativas sem instabilidade numérica, utilizando coeficientes complexos.
Codificação Direta: O método produz uma codificação em bloco (block encoding) direta e não escalada do operador de evolução dissipativa até uma alta ordem de precisão.
Derivação de Circuitos Quânticos: Desenvolvimento de circuitos específicos para a equação de onda amortecida clássica, incluindo a decomposição de operadores de propagação de ondas e amortecimento exponencial usando portas quânticas padrão (como rotações $R_Y$ e portas de fase).
Execução em Hardware Real: Implementação e teste bem-sucedido do algoritmo no processador quântico de íons aprisionados IonQ Forte 1 (36 qubits).

4. Resultados

Os resultados foram validados através de simulações de vetores de estado e execução em hardware real:

Simulações de Estado (Statevector):
- As fórmulas de ordem superior (4 e 6) atingiram a precisão teórica esperada ( $O(\Delta t^p)$ ).
- Para tolerâncias de erro baixas, as fórmulas de alta ordem tornaram-se significativamente mais eficientes em termos de contagem de portas CNOT do que as de baixa ordem, apesar de terem circuitos mais profundos por passo.
- A ordem 6 alcançou a precisão de máquina clássica com cerca de $10^4$ portas CNOT.
Execução em Hardware (IonQ Forte 1):
- O problema simulado foi a propagação de ondas mecânicas com perdas (equação de onda amortecida) em 16 pontos de grade.
- Comparação de Ordens:
  - Ordem 1: Inacurada.
  - Ordem 2 e 4: Ambas mostraram alta precisão. A Ordem 4 foi ligeiramente mais precisa que a Ordem 2, mesmo com maior profundidade de circuito e maior suscetibilidade ao ruído.
  - Ordem 6: A precisão degradou-se devido ao ruído do hardware (profundidade de circuito excessiva para o dispositivo atual), tornando-se menos precisa que a Ordem 1 neste cenário específico.
- Probabilidade de Sucesso: As probabilidades de pós-seleção (sucesso da medição do ancilla) para as ordens 2 e 4 foram próximas das probabilidades ideais, indicando que a evolução contrativa foi preservada corretamente.

5. Significado e Conclusão

Viabilidade Prática: O trabalho demonstra que a divisão de operadores de alta ordem é uma abordagem precisa e prática para simular dinâmicas dissipativas em processadores quânticos de curto prazo (NISQ), superando as limitações dos métodos de baixa ordem.
Eficiência de Recursos: Embora circuitos de alta ordem sejam mais profundos, eles reduzem drasticamente o número total de portas necessárias para atingir alta precisão, tornando-os vantajosos mesmo em regimes de baixa precisão.
Escalabilidade: A complexidade de portas CNOT escala favoravelmente com o tempo de simulação e a precisão desejada. Para a ordem 6, a dependência do tempo de simulação é quase linear ( $O(t^{7/6})$ ) e a dependência do erro é muito fraca ( $O(\epsilon^{-1/6})$ ), uma melhoria substancial sobre a abordagem de ordem 2 existente.
Futuro: O método abre caminho para simulações mais complexas em ciência computacional (como dinâmica de fluidos e transporte de escalares) em hardware quântico. Trabalhos futuros focarão em encontrar fórmulas de ordem ainda superior (8+) e otimizar a compressão de qubits ancilla para reduzir a sobrecarga de recursos.

Em resumo, o artigo resolve um obstáculo fundamental na simulação quântica de sistemas abertos, permitindo a aplicação de técnicas de alta precisão a problemas dissipativos reais em hardware quântico atual.