axionbloch: an Open-Source Python Package for Simulating Axion-Induced Spin Dynamics

Resumo Técnico de axionbloch: um Pacote Python de Código Aberto para Simulação de Dinâmica de Spin Induzida por Áxions

Problema
A interação entre a matéria escura de bósons ultraleves (áxions e partículas do tipo áxion, ou ALPs) e a matéria ordinária pode ser modelada como um campo pseudo-magnético atuando sobre spins fermiônicos. Esta interação é descrita por um Hamiltoniano onde o gradiente do campo de áxion acopla-se ao operador de spin, induzindo dinâmica de spin detectável via experimentos de medição de precisão, tais como Ressonância Magnética Nuclear (RMN) e buscas de comagnetômetros. Embora estudos teóricos tenham derivado assinaturas de sinais de áxion e sensibilidades experimentais, houve uma carência de ferramentas numéricas dedicadas para verificar essas assinaturas ou simular sistematicamente a dinâmica de spin sob diversas condições físicas e configurações de campo de áxion (ex: fundos de áxion cósmico, áxions de grão fino ou áxions acoplados quadraticamente). Além disso, à medida que novos modelos de áxion e propostas experimentais emergem, é necessário um arcabouço flexível para fazer a ponte entre modelos teóricos e observáveis experimentais.

Metodologia
O artigo apresenta o axionbloch, um pacote Python de código aberto projetado para simular a dinâmica de spin induzida tanto por campos magnéticos padrão quanto por campos pseudo-magnéticos induzidos por áxions. O pacote emprega a seguinte abordagem técnica:

• Modelagem Física: A física central baseia-se nas equações de Bloch, que descrevem a evolução da magnetização ($M$) na presença de campos magnéticos efetivos ($B_{eff}$). O pacote foca na ressonância magnética de spin-1/2, tratando o gradiente do campo de áxion como um campo pseudo-magnético adicionado ao campo de polarização (bias field). Para lidar com a disparidade entre a rápida precessão de Larmor do campo de polarização e a dinâmica mais lenta induzida pelo áxion, a simulação utiliza um Referencial de Coordenadas Rotativas (RCF).
• Integração Numérica: A evolução temporal da magnetização é resolvida usando um método de Runge–Kutta de quarta ordem (RK4). O algoritmo é implementado em uma biblioteca C++ dedicada para eficiência computacional, enquanto o fluxo de trabalho de alto nível e o gerenciamento de parâmetros são tratados em Python. As duas linguagens são interconectadas via pybind11.
• Arquitetura: O software é construído utilizando programação orientada a objetos. Os componentes principais incluem:
  • Sample: Encapsula parâmetros experimentais como tempos de relaxação ($T_1, T_2, T_2^*$) e propriedades da amostra.
  • Magnet e MagField: Gerenciam campos de polarização, sua inhomogeneidade e a geração de campos pseudo-magnéticos.
  • MilkyWayAxionHalo (e classes extensíveis): Modela configurações específicas de campos de áxion, tratando o campo como um sinal oscilante com frequência determinada pela massa do áxion e largura de linha determinada pela dispersão de velocidade.
  • Simulation e Simulations: Gerenciam a execução de uma ou múltiplas execuções de simulação, integrando configurações experimentais com modelos de áxion.
• Unidades e Constantes: O pacote utiliza astropy.units para o tratamento de grandezas físicas, garantindo consistência dimensional e facilidade de conversão entre os sistemas SI e CGS. Constantes comuns (ex: razões giromagnéticas) estão pré-definidas dentro do pacote.

Principais Contribuições e Resultados

1. Arcabouço de Código Aberto: A principal contribuição é o lançamento do axionbloch (disponível em https://github.com/Yuzhe98/AxionBloch), uma ferramenta modular que permite aos pesquisadores simular sinais pseudo-magnéticos sob modelos de áxion e configurações experimentais específicos.
2. Calibração e Validação: Os autores realizaram uma extensa calibração contra expectativas teóricas para garantir a precisão. Três marcos específicos foram estabelecidos:
   • RMN Pulsada: Verificou o decaimento livre ($T_2^*$) e comportamentos de eco de spin ($T_2$) após pulsos de 90° e 180°.
   • RMN de Onda Contínua (CW): Validou o crescimento da magnetização transversal sob campos oscilantes fracos.
   • Amostras Hiperpolarizadas: Confirmou o decaimento exponencial da polarização devido ao relaxamento $T_1$.
   • Precisão Quantitativa: A discrepância entre a magnetização simulada e a esperada foi quantificada usando uma métrica $\chi^2$. Os resultados de calibração mostraram $\chi^2 < 10^{-4}$ para RMN pulsada e valores ainda menores para CW-RMN e relaxamento $T_1$, confirmando a confiabilidade da implementação numérica.
3. Demonstração de Capacidade: O artigo apresenta uma simulação de sinais induzidos por áxion de um halo de matéria escura de áxion local virializado (modelo da Via Láctea). A simulação gerou com sucesso campos pseudo-magnéticos e a magnetização transversal resultante, demonstrando a capacidade do pacote de lidar com campos de áxion estocásticos e inhomogeneidades de campo de polarização.
4. Desempenho: A arquitetura híbrida Python/C++ permite uma computação eficiente. Um caso de teste envolvendo $7 \times 10^8$ passos de integração foi concluído em 2,7 segundos em um computador pessoal padrão, permitindo a exploração de vários regimes de parâmetros.

Significância
O artigo posiciona o axionbloch como uma ferramenta necessária para a simulação sistemática da dinâmica de spin induzida por áxions, preenchendo uma lacuna no campo onde a verificação numérica de assinaturas de sinal era anteriormente inexistente. Ao fornecer uma plataforma de código aberto, bem documentada e testada, o pacote visa facilitar o design de buscas experimentais e análise de dados para detecção de áxions. Ele serve como uma ponte entre modelos teóricos de áxion e observáveis experimentais, permitindo que pesquisadores derivem assinaturas de sinal para configurações específicas. Adicionalmente, os autores observam o potencial de utilidade do pacote para fins educacionais dentro da comunidade de RMN, permitindo a visualização de sequências de pulsos e dinâmica de spin. A ferramenta foi projetada para ser extensível, acomodando futuros modelos de áxion e propostas experimentais, incluindo potenciais implementações futuras para RMN de campo zero a ultra-baixo (ZULF).