Current neutrino detectors will observe hundreds to thousands of neutrinos from Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package that bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as the SuperNova Observatories with General Long Baseline Experiment Simulator (SNOwGLoBES). Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities.
In this work we present : A package dedicated to efficient computations of observables in the Early Universe with the focus on the cosmological era of Big Bang Nucleosynthesis (BBN). The code offers fast and precise evaluation of BBN light-element abundances together with the effective number of relativistic degrees of freedom, including non-instantaneous decoupling effects. is suitable for state-of-the-art analyses in the Standard Model as well as for general investigations into New Physics active during BBN. After reviewing the physics implemented in , we provide a short guide on how to use the code for applications in the Standard Model and beyond. The package is written in Python, but more advanced users can optionally take advantage of the open-source community for Julia. is publicly available on GitHub.
more » « less- Award ID(s):
- 2210283
- NSF-PAR ID:
- 10487782
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- The European Physical Journal C
- Volume:
- 84
- Issue:
- 1
- ISSN:
- 1434-6052
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Summary Computational cell-type deconvolution is an important analytic technique for modeling the compositional heterogeneity of bulk gene expression data. A conceptually new Bayesian approach to this problem, BayesPrism, has recently been proposed and has subsequently been shown to be superior in accuracy and robustness against model misspecifications by independent studies; however, given that BayesPrism relies on Gibbs sampling, it is orders of magnitude more computationally expensive than standard approaches. Here, we introduce the InstaPrism package which re-implements BayesPrism in a derandomized framework by replacing the time-consuming Gibbs sampling step with a fixed-point algorithm. We demonstrate that the new algorithm is effectively equivalent to BayesPrism while providing a considerable speed and memory advantage. Furthermore, the InstaPrism package is equipped with a precompiled, curated set of references tailored for a variety of cancer types, streamlining the deconvolution process.
Availability and implementation The package InstaPrism is freely available at: https://github.com/humengying0907/InstaPrism. The source code and evaluation pipeline used in this paper can be found at: https://github.com/humengying0907/InstaPrismSourceCode.
-
more » « less
We propose a non-exotic electromagnetic solution (within the standard model of particle physics) to the cosmological 7Li problem based upon a narrow 2 MeV photo-emission line from the decay of light glueballs (LGBs). These LGBs form within color superconducting quark clusters (SQCs), which are tens of Fermi in size, in the radiation-dominated post-BBN epoch. The mono-chromatic line from the LGB→γ+γ decay reduces Big Bang nucleosynthesis (BBN) 7Be by 2/3 without affecting other abundances or the cosmic microwave background (CMB) physics, provided the combined mass of the SQCs is greater than the total baryonic mass in the universe. Following the LGB emission, the in-SQC Quantum ChromoDynamics (QCD) vacuum becomes unstable and “leaks” (via quantum tunneling) into the external space-time (trivial) vacuum, inducing a decoupling of SQCs from hadrons. In seeking a solution to the 7Li problem, we uncovered a solution that also addresses the Dark Energy (DE) and dark matter (DM) problem, making these critical problems intertwined in our model. Being colorless, charge-neutral, optically thin, and transparent to hadrons, SQCs interact only gravitationally, making them a viable cold DM (CDM) candidate. The leakage (i.e., quantum tunneling) of the in-SQC QCD vacuum to the trivial vacuum offers an explanation of DE in our model and allows for a cosmology that evolves into a ΛCDM universe at a low redshift with a possible resolution of the Hubble tension. Our model distinguishes itself by proposing that the QCD vacuum within SQCs possesses the ability to tunnel into the exterior trivial vacuum, resulting in the generation of DE. This implies the possibility that DM and hadrons might represent distinct phases of quark matter within the framework of QCD, characterized by different vacuum properties. We discuss SQC formation in heavy-ion collision experiments at moderate temperatures and the possibility of detection of MeV photons from the LGB→γ+γ decay.
-
more » « less
SUMMARY PyLith, a community, open-source code for modelling quasi-static and dynamic crustal deformation with an emphasis on earthquake faulting, has recently been updated with a flexible multiphysics implementation. We demonstrate the versatility of the multiphysics implementation by extending the code to model fully coupled continuum poromechanics. We verify the newly incorporated physics using standard benchmarks for a porous medium saturated with a slightly compressible fluid. The benchmarks include the 1-D consolidation problem as outlined by Terzaghi, Mandel’s problem for the 2-D case, and Cryer’s problem for the 3-D case. All three benchmarks have been added to the PyLith continuous integration test suite. We compare the closed form analytical solution for each benchmark against solutions generated by our updated code, and lastly, demonstrate that the poroelastic material formulation may be used alongside the existing fault implementation in PyLith.
-
more » « less
Abstract We present a Python package geared toward the intuitive analysis and visualization of paleoclimate timeseries,
Pyleoclim . The code is open‐source, object‐oriented, and built upon the standard scientific Python stack, allowing users to take advantage of a large collection of existing and emerging techniques. We describe the code's philosophy, structure, and base functionalities and apply it to three paleoclimate problems: (a) orbital‐scale climate variability in a deep‐sea core, illustrating spectral, wavelet, and coherency analysis in the presence of age uncertainties; (b) correlating a high‐resolution speleothem to a climate field, illustrating correlation analysis in the presence of various statistical pitfalls (including age uncertainties); (c) model‐data confrontations in the frequency domain, illustrating the characterization of scaling behavior. We show how the package may be used for transparent and reproducible analysis of paleoclimate and paleoceanographic datasets, supporting Findable, Accessible, Interoperable, and Reusable software and an open science ethos. The package is supported by an extensive documentation and a growing library of tutorials shared publicly as videos and cloud‐executable Jupyter notebooks, to encourage adoption by new users.
