Algebraic and Arithmetic Attributes of Hypergeometric Functions in SageMath

Resumo Técnico: Atributos Algébricos e Aritméticos de Funções Hipergeométricas em SageMath

Enunciado do Problema
Funções hipergeométricas, definidas como séries de potências com parâmetros racionais, são objetos centrais no estudo de séries D-finitas. Embora suas propriedades algébricas e aritméticas — como álgebra, limitação global e comportamento sob redução de primos — sejam teoricamente classificadas, existe uma necessidade prática de ferramentas computacionais para formular e verificar conjecturas sobre essas propriedades. Especificamente, pesquisadores requerem software capaz de lidar com séries hipergeométricas sobre números racionais ($\mathbb{Q}$), corpos finitos ($\mathbb{F}_p$) e corpos $p$-ádicos ($\mathbb{Q}_p$) para testar hipóteses sobre grupos de Galois de reduções e outros fenômenos aritméticos. Os autores identificaram uma lacuna nos sistemas existentes de álgebra computacional quanto à implementação desses algoritmos específicos dentro de um framework unificado.

Metodologia
Os autores desenvolveram um novo pacote SageMath que implementa funções hipergeométricas como elementos de anéis simbólicos, permitindo a manipulação sobre anéis arbitrários (corpos racionais, anéis de polinômios, anéis de séries de potências e corpos $p$-ádicos). A implementação baseia-se em várias abordagens algorítmicas distintas adaptadas ao contexto aritmético específico:

1. Sobre $\mathbb{Q}$: O pacote implementa critérios para limitação global (global boundedness) e algebraicidade. Estes baseiam-se nas posições relativas das partes decimais dos parâmetros (critério de Christol) e extensões para casos envolvendo diferenças inteiras entre os parâmetros superiores e inferiores. O sistema também computa conjuntos de primos com boa redução analisando a valuação $p$-ádica dos coeficientes para pequenos primos e extrapolando com base em classes de congruência módulo o mínimo denominador comum dos parâmetros.
2. Sobre $\mathbb{F}_p$: Para séries com boa redução, o pacote utiliza operadores de seção para decompor funções hipergeométricas. Ele constrói relações de Dwork, expressando uma função como uma combinação linear polinomial de potências de ordem $p$ de outras funções hipergeométricas. Essas relações são usadas para computar polinômios anuladores (polinômios linearizados no operador de Frobenius). Além disso, o pacote inclui um algoritmo para verificar se duas séries hipergeométricas distintas são congruentes módulo $p$ comparando recursivamente suas seções e termos constantes. Para funções onde $m=n-1$, o pacote computa a $p$-curvatura do associado operador diferencial.
3. Sobre $\mathbb{Q}_p$: A implementação foca em convergência e valuação. Ela calcula o raio de convergência $p$-ádico logarítmico usando uma fórmula explícita baseada na partição dos parâmetros em inteiros $p$-ádicos e não-inteiros. Computa valuações $p$-ádicas com derivação de $\nu$ ($\nu$-drifted) usando uma recursão sobre o semiring tropical e o algoritmo de Floyd-Warshall para determinar o fechamento transitivo fraco de uma matriz tropical. Finalmente, computa polígonos de Newton e avalia as funções em pontos $p$-ádicos específicos dentro de seu raio de convergência.

Principais Contribuições

• Implementação Unificada: A criação de um pacote SageMath de código aberto que trata funções hipergeométricas sobre $\mathbb{Q}$, $\mathbb{F}_p$ e $\mathbb{Q}_p$ dentro de uma única interface.
• Realização Algorítmica: A tradução de classificações teóricas (ex: critérios de Christol, classificação de Beukers-Heckman) e propriedades aritméticas (relações de Dwork, $p$-curvatura) em código executável.
• Algoritmos Específicos:
  • is_globally_bounded() e is_algebraic() para parâmetros racionais.
  • good_reduction_primes() para determinar primos onde os denominadores dos coeficientes são coprimos a $p$.
  • section(), dwork_relation() e annihilating_ore_polynomial() para propriedades aritméticas em característica positiva.
  • is_equal_as_series() para detectar congruências entre diferentes séries hipergeométricas módulo $p$.
  • p_curvature() e p_curvature_coranks() para analisar operadores diferenciais em corpos finitos.
  • log_radius_of_convergence(), valuation() e newton_polygon() para análise $p$-ádica.
• Tratamento de Casos Limite: O sistema acomoda parâmetros gerais, incluindo inteiros não negativos nos parâmetros inferiores se compensados pelos parâmetros superiores, e lida com casos onde algoritmos podem teoricamente envolver sistemas infinitos (com um aviso de que a terminação não é estritamente provada para a implementação simplificada de polinômios anuladores).

Resultos
O artigo demonstra a funcionalidade do pacote através de vários exemplos:

• Caso Racional: Classifica com sucesso as funções $f$, $g$ e $h$ como limitadas globalmente ou algébricas e identifica classes de congruência específicas de primos para boa redução.
• Caso de Corpo Finito: Computa relações de Dwork para uma função módulo 19, derivando um dicionário de coeficientes polinomiais. Demonstra a computação de polinômios anuladores e verifica a congruência de duas séries hipergeométricas diferentes módulo 13. Também calcula a matriz de $p$-curvatura e seus corângs para vários primos.
• Caso $p$-ádico: Computa o raio de convergência logarítmico (ex: $-7/2$ para uma função específica em $p=5$), calcula valuações com derivação (retornando tuplas de valuação e índice mínimo) e constrói polígonos de Newton, lidando com sucesso com casos de inclinações infinitas por meio de truncamento.

Significância e Alegações
Os autores posicionam este trabalho como uma ferramenta prática para pesquisadores que trabalham em conjecturas relativas às propriedades algébricas e aritméticas de séries D-finitas, uma necessidade que surgiu durante seu trabalho conjunto sobre grupos de Galois de reduções módulo primos. O artigo afirma que o pacote fornece as funcionalidades necessárias para manipular funções hipergeométricas nesses diversos ambientes, facilitando o teste de previsões teóricas.

Os autores são modestos quanto à completude teórica da implementação; eles alertam explicitamente que o algoritmo para computar polinômios anuladores usa uma versão simplificada onde a finitude das séries hipergeométricas envolvidas não é provada, o que significa que as computações podem não terminar em todos os casos. No entanto, expressam confiança de que a implementação atual terminará na prática. O trabalho é apresentado como um passo em direção à integração em uma futura versão do SageMath, complementado por uma planilha interativa para exploração.