Resonant false vacuum decay in two dimensions on a 4000-qubit quantum annealer

Enunciado do Problema
Estados metaestáveis em física de muitos corpos, como o falso vácuo em teoria quântica de campos e cosmologia, tipicamente decaem através da nucleação de bolhas localizadas do vácuo verdadeiro, seguidas por sua expansão. A teoria convencional trata a nucleação e o crescimento como estágios separados, sendo a nucleação a etapa limitante do processo, governada pela teoria da bolha crítica. No entanto, a dinâmica do crescimento pós-nucleação, particularmente em duas dimensões, permanece menos compreendida. Especificamente, não está claro se o próprio crescimento pode se tornar o regime dinâmico dominante e não trivial, controlado por processos ressonantes locais na fronteira da bolha, em vez de um movimento de interface suave impulsionado por pressão. Embora o decaimento do falso vácuo (FVD) tenha sido explorado em sistemas unidimensionais e diversas plataformas experimentais, o crescimento real de bolhas de falso vácuo ressonantes em duas dimensões permaneceu amplamente inexplorado.

Metodologia
Os autores investigam este problema usando um recozimento quântico programável de larga escala (sistemas D-Wave Advantage2 com mais de 4000 qubits físicos) para realizar um modelo Ising quântico de campo misto bidimensional em uma topologia de toro. O modelo é definido pelo Hamiltoniano:
$$\hat{H} = -J \sum_{\langle i,j \rangle} \hat{\sigma}^z_i \hat{\sigma}^z_j - h_x \sum_i \hat{\sigma}^x_i - h_z \sum_i \hat{\sigma}^z_i$$
onde $J \gg h_x$. O sistema é preparado em um estado metaestável de "falso vácuo" (todos os spins para cima) com um único "seed" (semente) de spin invertido para o estado de "vácuo verdadeiro" (para baixo) no centro.

O estudo emprega três abordagens complementares:

1. Simulação Quântica: Usando um protocolo de recozimento reverso (reverse-annealing), os autores controlam o campo transversal dependente do tempo $h_x(t)$ e o campo longitudinal $h_z(t)$. Eles sintonizam especificamente $h_z$ para a condição de ressonância $|h_z| = 2J$. Neste valor, inversões de spin único na fronteira de uma bolha de vácuo verdadeiro existente tornam-se conservadoras de energia (ressonantes), permitindo uma expansão com restrição cinética.
2. Simulações de Redes de Tensores: Redes de Tensores de Árvore (TTNs) são usadas para simular a dinâmica quântica completa de sistemas menores ($16 \times 16$ spins) para capturar o emaranhamento e verificar efeitos quânticos de tempo inicial.
3. Modelagem de Circuito Estocástico (SC): Um modelo estocástico clássico efetivo é desenvolvido para simular a dinâmica de longo prazo de sistemas muito maiores (centenas de milhares de spins). Este modelo impõe a condição de ressonância local como uma restrição nas inversões de spin, permitindo a análise estatística do crescimento da bolha ao longo de escalas de tempo estendidas.

Principais Contribuições e Resultados

• Regime de Crescimento Ressonante: Os autores identificam um regime específico onde o crescimento de um domínio de vácuo verdadeiro supera a nucleação espontânea em aproximadamente três ordens de magnitude. Na ressonância $|h_z| = 2J$, o custo energético de criar duas paredes de domínio adicionais é exatamente compensado pelo ganho de energia do campo longitudinal ao inverter um único spin na fronteira da bolha. Isso faz com que o processo de crescimento seja um efeito de primeira ordem em $h_x$, enquanto a nucleação espontânea de novas bolhas permanece um processo de ordem superior (quarta ordem).
• Expansão Balística: Dados experimentais e simulações revelam que a bolha semeada expande-se com dinâmica quase balística. A frente de propagação segue uma estrutura de cone de luz com um expoente $p \approx 1$ (onde o raio $r \sim \tau^p$), indicando que a interface se move a uma velocidade constante determinada pela ressonância.
• Alargamento de Interface KPZ: Embora a expansão seja balística, a interface não é rígida. O estudo observa um alargamento sub-balístico da frente de onda. Ao analisar o escalonamento da largura da interface, os autores descobrem que a dinâmica pertence à classe de universalidade de Kardar–Parisi–Zhang (KPZ), caracterizada por um expoente de escala de $2/3$. Isso se mantém válido através das simulações quânticas, cálculos de TTN e o modelo de circuito estocástico.
• Coerência Quântica: Os autores fornecem evidências de que a dinâmica observada é genuinamente quântica. O tempo gasto no regime ressonante (dezenas de nanossegundos) é consistente com escalas de tempo de operação coerente. Além disso, o escalonamento observado da magnetização $M(t) \sim h_x^2 \tau$ decorre da natureza quântica da ressonância (largura de ressonância $\delta h_z \propto h_x$), um recurso reproduzido por TTNs, mas distinto da ativação térmica clássica.
• Isotropia Emergente: Apesar da simetria $D_4$ subjacente da rede quadrada, o crescimento médio da bolha exibe uma notável invariância esférica, sugerindo que as regras de ressonância local levam a uma propagação isotrópica emergente.

Significância
O artigo estabelece um regime de "crescimento dominado" do decaimento do falso vácuo, onde a expansão de domínios de vácuo verdadeiro é governada por processos de ressonância local, em vez do padrão semiclássico de nucleação e crescimento. Os resultados demonstram que a simulação quântica de larga escala pode acessar a dinâmica fora do equilíbrio relevante para a teoria quântica de campos e cosmologia. Especificamente, o trabalho mostra que o hardware quântico programável pode isolar e estudar canais específicos de crescimento pós-nucleação, revelando que as propriedades estatísticas assintóticas de tal crescimento (expansão balística com alargamento KPZ) podem ser capturadas por modelos clássicos efetivos, embora a dinâmica microscópica subjacente gere emaranhamento substancial. Isso sugere que o comportamento estatístico de larga escala do crescimento de bolhas ressonantes pode ser robusto contra detalhes quânticos microscópicos, ao mesmo tempo em que destaca o potencial dos recedores quânticos para explorar dinâmicas de muitos corpos complexas além do alcance dos métodos numéricos clássicos atuais.