snewpy

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

SNEWPY es un paquete de Python para trabajar con neutrinos de supernova. Ofrece…

• … una interfaz sencilla y unificada para cientos de simulaciones de supernovas.
• … una gran biblioteca de transformaciones de sabor que relacionan los flujos de neutrinos producidos en la supernova con los que llegan a un detector en la Tierra.
• … y una interfaz Python para SNOwGLoBES que le permite estimar y trazar tasas de eventos en muchos detectores de neutrinos diferentes.

Instalación

Ejecute pip install snewpy para instalar SNEWPY.

SNEWPY incluye una gran cantidad de modelos de supernovas de diferentes grupos de simulación. Dado que estos modelos tienen un tamaño de varios 100 MB, no se incluyen en la instalación inicial, pero se descargarán automáticamente cuando sea necesario. Alternativamente, puede ejecutar el siguiente comando para descargar explícitamente los modelos que desea usar en un subdirectorio llamado SNEWPY-models/<model_name>/ en el directorio actual:

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

Uso y documentación

SNEWPY le brinda fácil acceso a cientos de simulaciones SN incluidas...

 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 )

…y muchas transformaciones de sabor que los neutrinos podrían experimentar en su camino a la Tierra…

 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 )

… permitiéndote calcular rápidamente el flujo de neutrinos que llega a la Tierra:

 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 )

También puede calcular la tasa de eventos observados en todos los detectores de neutrinos compatibles con SNOwGLoBES, utilizar los modelos SN incluidos y las transformaciones de tipo en código de terceros (como generadores de eventos) y mucho más.

Los cuadernos de Jupyter que muestran los modelos de supernova descargables disponibles a través de SNEWPY y gran parte de su funcionalidad están disponibles en la subcarpeta doc/nb/ . Se encuentran scripts de ejemplo adicionales en la subcarpeta python/snewpy/scripts/ .

Los artículos que describen SNEWPY y la física subyacente se publican en Astrophysical Journal (DOI:10.3847/1538-4357/ac350f, arXiv:2109.08188) y Journal of Open Source Software (DOI:10.21105/joss.03772).

Para obtener más información, consulte la documentación completa en Leer los documentos.

Contribuyendo

¡Tus contribuciones a SNEWPY son bienvenidas! Para cambios menores, simplemente envíe una solicitud de extracción. Si planea cambios mayores, probablemente sea una buena idea abrir una edición primero para coordinar nuestro trabajo.

Usamos un flujo de trabajo Fork & Pull Request, que es común en GitHub. Consulte la página de contribución en nuestra documentación completa para obtener más detalles.