Presolar graphite grains carry the isotopic signatures of their parent stars. A significant fraction of presolar graphites shows isotopic abundance anomalies relative to solar for elements such as O, Si, Mg and Ca, which are compatible with nucleosynthesis in core-collapse supernovae (CCSNe). Therefore, they must have condensed from CCSN ejecta before the formation of the Sun. Their most puzzling abundance signature is the $^{22}$Ne-enriched component Ne-E(L), interpreted as the effect of the radioactive decay of $^{22}$Na ($T_{1/2}$ = 2.6 years). Previous works have shown that if H is ingested into the He shell and not fully destroyed before the explosion, the CCSN shock in the He shell material produces large $^{22}$Na amounts. Here we focus on such CCSN models, showing a radioactive $^{26}$Al production compatible with grains measurements, and analyze the conditions of $^{22}$Na nucleosynthesis. In these models, $^{22}$Na is mostly made in the He shell, with a total ejected mass varying between 2.6$\times$10$^{-3}$ M$_{\odot}$ and 1.9$\times$10$^{-6}$ M$_{\odot}$. We show that such $^{22}$Na may already impact the CCSN light curve 500 days after the explosion, and at later stages it can be the main source powering the CCSN light curve for up to a few years before the $^{44}$Ti decay becomes dominant. Based on the CCSN yields above, the 1274.53 keV $\gamma$-ray flux due to $^{22}$Na decay could be observable for years after the first CCSN light is detected, depending on the distance. This makes CCSNe possible sites to detect a $^{22}$Na $\gamma$-ray signature consistently with the Ne-E(L) component found in presolar graphites. Finally, we discuss the potential contribution from the $^{22}$Na decay to the galactic positron annihilation rate.

The enclosed data release consists of a subset of the calibration data from the Majorana Demonstrator experiment. Each Majorana event is accompanied by raw Germanium detector waveforms, pulse shape discrimination cuts, and calibrated final energies, all shared in an HDF5 file format along with relevant metadata. This release is specifically designed to support the training and testing of Artificial Intelligence (AI) and Machine Learning (ML) algorithms upon our data. This document is structured as follows. Section I provides an overview of the dataset's content and format; Section II outlines the location of this dataset and the method for accessing it; Section III presents the NPML Machine Learning Challenge associated with this dataset; Section IV contains a disclaimer from the Majorana collaboration regarding the use of this dataset; Appendix A contains technical details of this data release. Please direct questions about the material provided within this release to liaobo77@ucsd.edu (A. Li).

DMRadio-m$^3$ is an experiment that is designed to be sensitive to KSVZ and DFSZ QCD axion models in the 10--200\,MHz (41 neV$/c^2$ -- 0.83 $\mu$eV/$c^2$) range. The experiment uses a solenoidal dc magnetic field to convert an axion dark-matter signal to an ac electromagnetic response in a coaxial copper pickup. The current induced by this axion signal is measured by dc SQUIDs. In this work, we present the electromagnetic modeling of the response of the experiment to an axion signal over the full frequency range of DMRadio-m$^3$, which extends from the low-frequency, lumped-element limit to a regime where the axion Compton wavelength is only a factor of two larger than the detector size. With these results, we determine the live time and sensitivity of the experiment. The primary science goal of sensitivity to DFSZ axions across 30--200 MHz can be achieved with a $3\sigma$ live scan time of 2.9 years.

Sapphire has mechanical and electrical properties that are advantageous for the construction of internal components of radiation detectors such as time projection chambers and bolometers. However, it has proved difficult to assess its $\rm ^{232}Th$ and $\rm ^{238}U$ content down to the picogram per gram level. This work reports an experimental verification of a computational study that demonstrates $\gamma\gamma$ coincidence counting, coupled with neutron activation analysis (NAA), can reach ppt sensitivities. Combining results from $\gamma\gamma$ coincidence counting with those of earlier single-$\gamma$ counting based NAA shows that a sample of Saint Gobain sapphire has $\rm ^{232}Th$ and $\rm ^{238}U$ concentrations of <0.26 ppt and <2.3 ppt, respectively; the best constraints on the radiopurity of sapphire.

We report the results of the second measurement campaign of the Karlsruhe Tritium Neutrino (KATRIN) experiment. KATRIN probes the effective electron anti-neutrino mass, $m_{\nu}$, via a high-precision measurement of the tritium $\beta$-decay spectrum close to its endpoint at $18.6\,\mathrm{keV}$. In the second physics run presented here, the source activity was increased by a factor of 3.8 and the background was reduced by $25\,\%$ with respect to the first campaign. A sensitivity on $m_{\nu}$ of $0.7\,\mathrm{eV/c^2}$ at $90\,\%$ confidence level (CL) was reached. This is the first sub-eV sensitivity from a direct neutrino-mass experiment. The best fit to the spectral data yields $m_{\nu}^2 = (0.26\pm0.34)\,\mathrm{eV^4/c^4}$, resulting in an upper limit of m_{\nu}<0.9\,\mathrm{eV/c^2} ($90\,\%$ CL). By combining this result with the first neutrino mass campaign, we find an upper limit of m_{\nu}<0.8\,\mathrm{eV/c^2} ($90\,\%$ CL).

The electro- and photo-production of $J/\psi$ meson near the threshold from the proton is relevant to the search of hidden charm pentaquark candidates reported by the LHCb collaboration, and the study of the QCD trace anomaly's contribution to the proton mass. It is also expected to be sensitive to the QCD van der Waals interaction, that is mediated by multi-gluon exchanges and expected to dominate the interaction between two hadrons with no common valence quarks. Subthreshold production of $J/\psi$ from a nuclear target is expected to enhance such attractive interaction, and also allows for a direct probe of short range correlations inside a nucleus. With the high luminosity capability of the 12-GeV CEBAF facility at Jefferson Lab, high-precision data on $J/\psi$ meson production from the proton is becoming available, providing also a reference for subthreshold $J/\psi$ production from the deuteron. Data from the deuteron will establish the baseline for subthreshold $J/\psi$ production from other nuclear targets. In this paper, we present our findings from a feasibility study of subthreshold $J/\psi$ production from the deuteron using the proposed Solenoidal Large Intensity Device (SoLID), and discuss the potential physics impact of such data.

The projected sensitivity of the effective electron neutrino-mass measurement with the KATRIN experiment is below 0.3 eV (90 % CL) after five years of data acquisition. The sensitivity is affected by the increased rate of the background electrons from KATRIN's main spectrometer. A special shifted-analysing-plane (SAP) configuration was developed to reduce this background by a factor of two. The complex layout of electromagnetic fields in the SAP configuration requires a robust method of estimating these fields. We present in this paper a dedicated calibration measurement of the fields using conversion electrons of gaseous $^\mathrm{83m}$Kr, which enables the neutrino-mass measurements in the SAP configuration.

Radioactive isotopes produced through cosmic muon spallation are a background for rare-event detection in $\nu$ detectors, double-$\beta$-decay experiments, and dark-matter searches. Understanding the nature of cosmogenic backgrounds is particularly important for future experiments aiming to determine the pep and CNO solar neutrino fluxes, for which the background is dominated by the spallation production of $^{11}$C. Data from the Kamioka liquid-scintillator antineutrino detector (KamLAND) provides valuable information for better understanding these backgrounds, especially in liquid scintillators, and for checking estimates from current simulations based upon MUSIC, FLUKA, and GEANT4. Using the time correlation between detected muons and neutron captures, the neutron production yield in the KamLAND liquid scintillator is measured to be $(2.8 \pm 0.3) \times 10^{-4} \mu^{-1} g^{-1} cm^{2}$. For other isotopes, the production yield is determined from the observed time correlation related to known isotope lifetimes. We find some yields are inconsistent with extrapolations based on an accelerator muon beam experiment.

The precision measurement of the tritium $\beta$-decay spectrum performed by the KATRIN experiment provides a unique way to search for general neutrino interactions (GNI). All theoretical allowed GNI terms involving neutrinos are incorporated into a low-energy effective field theory, and can be identified by specific signatures in the measured tritium $\beta$-spectrum. In this paper an effective description of the impact of GNI on the $\beta$-spectrum is formulated and the first constraints on the effective GNI parameters are derived based on the 4 million electrons collected in the second measurement campaign of KATRIN in 2019. In addition, constraints on selected types of interactions are investigated, thereby exploring the potential of KATRIN to search for more specific new physics cases, including a right-handed W boson, a charged Higgs or leptoquarks.

High-Purity Germanium (HPGe) detectors are a key technology for rare-event searches such as neutrinoless double-beta decay ($0\nu\beta\beta$) and dark matter experiments. Pulse shapes from these detectors vary with interaction topology and thus encode information critical for event classification. Pulse shape simulations (PSS) are essential for modeling analysis cuts that distinguish signal events from backgrounds and for generating reliable simulations of energy spectra. Traditional PSS methods rely on a series of first-principles corrections to replicate the effect of readout electronics, requiring challenging fits over large parameter spaces and often failing to accurately model the data. We present a neural network architecture, the Cyclic Positional U-Net (CPU-Net), that performs translations of simulated pulses so that they closely resemble measured detector signals. Using a Cycle Generative Adversarial Network (CycleGAN) framework, this Response Emulation Network (REN) learns a data-driven mapping between simulated and measured pulses with high fidelity, without requiring a predetermined response model. We use data from a High-Purity Germanium (HPGe) detector with an inverted-coaxial point contact (ICPC) geometry to show that CPU-Net effectively captures and reproduces critical pulse shape features, allowing more realistic simulations without detector-specific tuning. CPU-Net achieves up to a factor-of-four improvement in distribution-level agreement for pulse shape parameter reconstruction, while preserving the topology-dependent information required for pulse-shape discrimination.

We present a search for neutrinoless double-beta ($0\nu\beta\beta$) decay of $^{136}$Xe using the full KamLAND-Zen 800 dataset with 745 kg of enriched xenon, corresponding to an exposure of $2.097$ ton yr of $^{136}$Xe. This updated search benefits from a more than twofold increase in exposure, recovery of photo-sensor gain, and reduced background from muon-induced spallation of xenon. Combining with the search in the previous KamLAND-Zen phase, we obtain a lower limit for the $0\nu\beta\beta$ decay half-life of $T_{1/2}^{0\nu} > 3.8 \times 10^{26}$ yr at 90% C.L., a factor of 1.7 improvement over the previous limit. The corresponding upper limits on the effective Majorana neutrino mass are in the range 28-122 meV using phenomenological nuclear matrix element calculations.

The MAJORANA DEMONSTRATOR searched for neutrinoless double-$\beta$ decay ($0\nu\beta\beta$) of $^{76}$Ge using modular arrays of high-purity Ge detectors operated in vacuum cryostats in a low-background shield. The arrays operated with up to 40.4 kg of detectors (27.2 kg enriched to $\sim$88\% in $^{76}$Ge). From these measurements, the DEMONSTRATOR has accumulated 64.5 kg yr of enriched active exposure. With a world-leading energy resolution of 2.52 keV FWHM at the 2039 keV $Q_{\beta\beta}$ (0.12\%), we set a half-life limit of $0\nu\beta\beta$ in $^{76}$Ge at $T_{1/2}>8.3\times10^{25}$ yr (90\% C.L.). This provides a range of upper limits on $m_{\beta\beta}$ of $(113-269)$ meV (90\% C.L.), depending on the choice of nuclear matrix elements.

The authors review the evidence for the applicability of random--matrix theory to nuclear spectra. In analogy to systems with few degrees of freedom, one speaks of chaos (more accurately: quantum chaos) in nuclei whenever random--matrix predictions are fulfilled. An introduction into the basic concepts of random--matrix theory is followed by a survey over the extant experimental information on spectral fluctuations, including a discussion of the violation of a symmetry or invariance property. Chaos in nuclear models is discussed for the spherical shell model, for the deformed shell model, and for the interacting boson model. Evidence for chaos also comes from random--matrix ensembles patterned after the shell model such as the embedded two--body ensemble, the two--body random ensemble, and the constrained ensembles. All this evidence points to the fact that chaos is a generic property of nuclear spectra, except for the ground--state regions of strongly deformed nuclei.

Nucleons (protons and neutrons) are the building blocks of atomic nuclei, and are responsible for more than 99\% of the visible matter in the universe. Despite decades of efforts in studying its internal structure, there are still a number of puzzles surrounding the proton such as its spin, and charge radius. Accurate knowledge about the proton charge radius is not only essential for understanding how quantum chromodynamics (QCD) works in the non-perturbative region, but also important for bound state quantum electrodynamics (QED) calculations of atomic energy levels. It also has an impact on the Rydberg constant, one of the most precisely measured fundamental constants in nature. This article reviews the latest situation concerning the proton charge radius in light of the new experimental results from both atomic hydrogen spectroscopy and electron scattering measurements, with particular focus on the latter. We also present the related theoretical developments and backgrounds concerning the determination of the proton charge radius using different experimental techniques. We discuss upcoming experiments, and briefly mention the deuteron charge radius puzzle at the end.

We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on Carbon Oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M$_\odot$ to 1.35 M$_\odot$) and the composition of the accreted material. Our simulations are guided by the results of multi-dimensional studies of TNRs in WDs that find sufficient mixing with WD core material occurs after the TNR is well underway, reaching levels of enrichment that agree with observations of CN ejecta abundances. We use NOVA (our 1-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial $^{12}$C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of $^7$Be in the ejected gases implying that CO CNe may be responsible for a significant fraction ($\sim$ 100 M$_\odot$) of the $^7$Li in the galaxy ($\sim$1000 M$_\odot$). In addition, although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion-outburst-accretion cycle and CO CNe may be an important channel of Supernova Ia progenitors.

The PRad experiment has credibly demonstrated the advantages of the calorimetric method in e-p scattering experiments to measure the proton root-mean-square (RMS) charge radius with high accuracy. The PRad result, within its experimental uncertainties, is in agreement with the small radius measured in muonic hydrogen spectroscopy experiments and it was a critical input in the recent revision of the CODATA recommendation for the proton charge radius. Consequently, the PRad result is in direct conflict with all modern electron scattering experiments. Most importantly, it is 5.8% smaller than the value from the most precise electron scattering experiment to date, and this difference is about three standard deviations given the precision of the PRad experiment. As the first experiment of its kind, PRad did not reach the highest precision allowed by the calorimetric technique. Here we propose a new (and) upgraded experiment -- PRad-II, which will reduce the overall experimental uncertainties by a factor of 3.8 compared to PRad and address this as yet unsettled controversy in subatomic physics. In addition, PRad-II will be the first lepton scattering experiment to reach the Q^2 range of 10^{-5} GeV^2 allowing a more accurate and robust extraction of the proton charge radius. The muonic hydrogen result with its unprecedented precision (~0.05%) determines the CODATA value of the proton charge radius, hence, it is critical to evaluate possible systematic uncertainties of those experiments, such as the laser frequency calibration that was raised in recent review articles. The PRad-II experiment with its projected total uncertainty of 0.43% could demonstrate whether there is any systematic difference between $e-p$ scattering and muonic hydrogen results. PRad-II will establish a new precision frontier in electron scattering and open doors for future physics opportunities.

A method based on Monte Carlo techniques is presented for evaluating thermonuclear reaction rates. We begin by reviewing commonly applied procedures and point out that reaction rates that have been reported up to now in the literature have no rigorous statistical meaning. Subsequently, we associate each nuclear physics quantity entering in the calculation of reaction rates with a specific probability density function, including Gaussian, lognormal and chi-squared distributions. Based on these probability density functions the total reaction rate is randomly sampled many times until the required statistical precision is achieved. This procedure results in a median (Monte Carlo) rate which agrees under certain conditions with the commonly reported recommended "classical" rate. In addition, we present at each temperature a low rate and a high rate, corresponding to the 0.16 and 0.84 quantiles of the cumulative reaction rate distribution. These quantities are in general different from the statistically meaningless "minimum" (or "lower limit") and "maximum" (or "upper limit") reaction rates which are commonly reported. Furthermore, we approximate the output reaction rate probability density function by a lognormal distribution and present, at each temperature, the lognormal parameters miu and sigma. The values of these quantities will be crucial for future Monte Carlo nucleosynthesis studies. Our new reaction rates, appropriate for bare nuclei in the laboratory, are tabulated in the second paper of this series (Paper II). The nuclear physics input used to derive our reaction rates is presented in the third paper of this series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.

Globular clusters are of paramount importance for testing theories of stellar evolution and early galaxy formation. Strong evidence for multiple populations of stars in globular clusters derives from observed abundance anomalies. A puzzling example is the recently detected Mg-K anticorrelation in NGC 2419. We perform Monte Carlo nuclear reaction network calculations to constrain the temperature-density conditions that gave rise to the elemental abundances observed in this elusive cluster. We find a correlation between stellar temperature and density values that provide a satisfactory match between simulated and observed abundances in NGC 2419 for all relevant elements (Mg, Si, K, Ca, Sc, Ti, and V). Except at the highest densities ($\rho \gtrsim 10^8$~g/cm$^3$), the acceptable conditions range from $\approx$ $100$~MK at $\approx$ $10^8$~g/cm$^3$ to $\approx$ $200$~MK at $\approx$ $10^{-4}$~g/cm$^3$. This result accounts for uncertainties in nuclear reaction rates and variations in the assumed initial composition. We review hydrogen burning sites and find that low-mass stars, AGB stars, massive stars, or supermassive stars cannot account for the observed abundance anomalies in NGC 2419. Super-AGB stars could be viable candidates for the polluter stars if stellar model parameters can be fine-tuned to produce higher temperatures. Novae, either involving CO or ONe white dwarfs, could be interesting polluter candidates, but a current lack of low-metallicity nova models precludes firmer conclusions. We also discuss if additional constraints for the first-generation polluters can be obtained by future measurements of oxygen, or by evolving models of second-generation low-mass stars with a non-canonical initial composition.