High-Precision Computation and PSLQ Identification of Stokes Multipliers for Anharmonic Oscillators

Título: Computação de Alta Precisão e Identificação PSLQ de Multiplicadores de Stokes para Osciladores Anarmônicos

Autor: Jian Zhou (Pesquisador Independente, XD.Inc, Xangai, China)
Data: Março de 2026

---

1. O Problema

O artigo aborda o cálculo dos multiplicadores de Stokes ($C_M$) para osciladores anarmônicos unidimensionais definidos pelo Hamiltoniano:
$$H = \frac{p^2}{2} + \frac{x^2}{2} + g x^{2M}$$
onde $M = 2, 3, \dots, 11$.

O problema central reside na natureza divergente da série de perturbação para a energia do estado fundamental. O comportamento de grande ordem dos coeficientes da série ($a_k$) é governado por soluções de instantons (soluções clássicas no tempo imaginário) e contém informações não perturbativas. Especificamente, os coeficientes crescem fatorialmente de acordo com:
$$a_k \sim C_M \cdot (-1)^{k+1} \cdot A_M^{-k} \cdot \Gamma((M-1)k + b + 1)$$
Enquanto a ação do instanton ($A_M$) e o deslocamento da função gama ($b$) são bem compreendidos teoricamente, expressões de forma fechada exata para o multiplicador de Stokes $C_M$ eram conhecidas apenas para o caso quartico ($M=2$). Para $M \ge 3$, a complexidade das soluções de instantons e dos determinantes de flutuação tornava o cálculo analítico extremamente difícil.

2. Metodologia

O autor desenvolveu um pipeline computacional de grande escala combinando aritmética de precisão arbitrária, aceleração de sequências e o algoritmo de relação inteira PSLQ. O processo divide-se em três etapas principais:

1. Recursão de Rayleigh-Schrödinger:

   • Cálculo de até 1200 coeficientes de perturbação ($a_k$) para cada $M$.
   • Utilização da biblioteca mpmath do Python com precisão de trabalho de até 300 dígitos (incluindo dígitos de guarda).
   • A matriz de elementos é construída no basis do oscilador harmônico, com complexidade $O(Mk)$ por ordem.

2. Extração de Parâmetros e Aceleração:

   • Ação do Instanton ($A_M$): Extraída via razão de coeficientes, acelerada por Extrapolação de Richardson (ordem 40–100).
   • Deslocamento Gama ($b$): Confirmado numericamente como $b = -1/2$ para todos os casos, consistente com argumentos de determinantes de flutuação de um loop.
   • Multiplicador de Stokes ($C_M$): A sequência convergente para $C_M$ é acelerada usando Extrapolação de Richardson de alta ordem para extrair 17 a 30 dígitos significativos confiáveis.

3. Identificação PSLQ:

   • O algoritmo PSLQ é aplicado para buscar relações inteiras entre $\ln|C_M|$ e uma base de constantes matemáticas (valores da função Gama em frações racionais, $\pi$, e números algébricos como 2 e 3).
   • A busca foi realizada com limites de coeficientes inteiros até 2000 e precisões variáveis (até 40 dígitos).

3. Principais Contribuições e Resultados

A. Novas Identidades Exatas

O estudo descobriu e verificou numericamente formas fechadas exatas para os multiplicadores de Stokes em casos que nunca foram publicados anteriormente:

• $M=3$ (Oscilador Sextico): $C_3^2 \pi^4 = 32$.
• $M=5$ (Oscilador Decico): $C_5^4 \Gamma(1/4)^4 \pi^5 = 2^{12} \cdot 3^2$.
• $M=7$ (Oscilador Tetradecico): $C_7^6 \Gamma(1/3)^9 \pi^6 = 2^{20} \cdot 3^3$.
• Confirmação do caso conhecido $M=2$: $C_2^2 \pi^3 = 6$.

B. O "Anomalia" $M=4$ (Oscilador Octico)

Um resultado negativo significativo foi encontrado para $M=4$:

• Após buscas exaustivas do PSLQ com 30 dígitos de precisão e limites de coeficientes até 2000, nenhuma forma fechada foi encontrada para $C_4$ na base padrão $\{\Gamma(1/3), \pi, 2, 3\}$.
• O valor numérico foi calculado com 30 dígitos confiáveis: $C_4 \approx -0.740051498259358506401511491622$.
• Isso sugere fortemente que o caso $x^8$ introduz um novo número transcendente genuíno, possivelmente relacionado ao determinante de flutuação que não se reduz a combinações algébricas de $\Gamma(1/3)$ e $\pi$.

C. Padrão Teórico-Numerário

O autor identifica uma correlação entre a existência de formas fechadas e a função totiente de Euler, $\phi(n)$.

• A existência de formas fechadas parece correlacionar-se com o número de transcendentais de Gama independentes, dado por $\phi(M-1)/2$.
• Para $\phi(M-1)/2 \le 1$ (casos $M \in \{2, 3, 4, 5, 7\}$), todos exceto $M=4$ possuem formas fechadas.
• Para $\phi(M-1)/2 \ge 2$ (casos $M \ge 6$), não foram encontradas formas fechadas simples, sugerindo que a complexidade algébrica aumenta com o número de transcendentais independentes.

4. Significado e Impacto

1. Validação da Teoria de Resurgência: Os resultados fornecem benchmarks de alta precisão para a análise de séries divergentes e a teoria de resurgência (resurgence), conectando setores perturbativos e não perturbativos.
2. Conexão com Teoria dos Números: O trabalho estabelece uma ponte intrigante entre a mecânica quântica (determinantes de flutuação de instantons) e a teoria dos números (propriedades da função Gama e a função totiente de Euler).
3. Mistério do Caso $M=4$: A ausência de uma forma fechada para o oscilador octico, enquanto o caso $M=7$ (que compartilha a mesma ação de instanton envolvendo $\Gamma(1/3)$) possui uma, indica que a complexidade reside especificamente no determinante de flutuação do operador de oscilação, e não na ação clássica.
4. Reprodutibilidade: Todo o código, dados e scripts de busca PSLQ foram disponibilizados publicamente, permitindo a verificação completa dos resultados.

Conclusão

Este artigo representa um avanço significativo na compreensão das propriedades assintóticas de osciladores anarmônicos. Ao combinar computação de precisão extrema com algoritmos de reconhecimento de constantes, o autor não apenas descobriu novas identidades matemáticas, mas também isolou um caso anômalo ($M=4$) que desafia a classificação atual baseada apenas na primalidade de $M$ ou na ação do instanton, apontando para uma estrutura algébrica mais profunda e complexa nos determinantes de flutuação quântica.