snewpy

SNEWPY: Modelos de alerta precoce de neutrinos de supernova para Python

SNEWPY é um pacote Python para trabalhar com neutrinos de supernovas. Oferece…

• …uma interface simples e unificada para centenas de simulações de supernovas.
• … uma grande biblioteca de transformações de sabor que relacionam os fluxos de neutrinos produzidos na supernova com aqueles que atingem um detector na Terra.
• … e uma interface Python para SNOwGLoBES que permite estimar e traçar taxas de eventos em muitos detectores de neutrinos diferentes.

Instalação

Execute pip install snewpy para instalar o SNEWPY.

SNEWPY inclui um grande número de modelos de supernovas de diferentes grupos de simulação. Como esses modelos têm vários 100 MB de tamanho, eles não estão incluídos na instalação inicial, mas serão baixados automaticamente quando necessário. Alternativamente, você pode executar o seguinte comando para fazer download explicitamente dos modelos que deseja usar para um subdiretório chamado SNEWPY-models/<model_name>/ no diretório atual:

python -c 'import snewpy; snewpy.get_models()'

Uso e Documentação

SNEWPY oferece acesso fácil a centenas de simulações SN incluídas…

 import astropy . units as u
from snewpy . models . ccsn import Nakazato_2013 , Bollig_2016

# Initialise two SN models. This automatically downloads the required data files if necessary.
nakazato = Nakazato_2013 ( progenitor_mass = 20 * u . solMass , revival_time = 100 * u . ms , metallicity = 0.004 , eos = 'shen' )
bollig = Bollig_2016 ( progenitor_mass = 27 * u . solMass )

… e muitas transformações de sabor que os neutrinos poderiam experimentar no caminho para a Terra…

 from snewpy . flavor_transformation import AdiabaticMSW
from snewpy . neutrino import MassHierarchy

# Adiabatic MSW flavor transformation with normal mass ordering
msw_nmo = AdiabaticMSW ( mh = MassHierarchy . NORMAL )

… permitindo calcular rapidamente o fluxo de neutrinos que chega à Terra:

 times = [ 0.5 , 1 ] * u . s
energies = range ( 5 , 50 ) * u . MeV
# Assume a SN at the fiducial distance of 10 kpc and normal mass ordering.
flux = bollig . get_flux ( times , energies , distance = 10 * u . kpc , flavor_xform = msw_nmo )

Você também pode calcular a taxa de eventos observada em todos os detectores de neutrinos suportados pelo SNOwGLoBES, usar os modelos SN incluídos e transformações de sabor em código de terceiros (como geradores de eventos) e muito mais.

Os notebooks Jupyter apresentando os modelos de supernova para download disponíveis através do SNEWPY e muitas de suas funcionalidades estão disponíveis na subpasta doc/nb/ . Scripts de exemplo adicionais estão na subpasta python/snewpy/scripts/ .

Artigos que descrevem o SNEWPY e a física subjacente são publicados no Astrophysical Journal (DOI:10.3847/1538-4357/ac350f, arXiv:2109.08188) e no Journal of Open Source Software (DOI:10.21105/joss.03772).

Para obter mais informações, consulte a documentação completa em Leia os documentos.

Contribuindo

Suas contribuições para o SNEWPY são bem-vindas! Para pequenas alterações, basta enviar uma solicitação pull. Se você planeja mudanças maiores, provavelmente é uma boa ideia abrir um problema primeiro para coordenar nosso trabalho.

Usamos um fluxo de trabalho Fork & Pull Request, que é comum no GitHub. Consulte a página Contribuindo em nossa documentação completa para obter detalhes.