The nucleardatapy toolkit for simple access to experimental nuclear data, astrophysical observations, and theoretical predictions

This paper introduces the Python toolkit \texttt{nucleardatapy} to facilitate transparent and unified access to experimental nuclear data, astrophysical observations, and theoretical predictions, while demonstrating its utility through a meta-analysis of dense matter properties that aligns nuclear-physics constraints with gravitational-wave findings.

The Gist

Imagine you are trying to build the ultimate LEGO castle, but the instructions are scattered across a billion different boxes, written in a dozen different languages, and some of the bricks are made of invisible material. This is the current state of nuclear physics and astrophysics. Scientists have data on how atoms stick together, how stars explode, and how neutron stars (the densest objects in the universe) behave, but it's all messy, disconnected, and hard to compare.

This paper introduces a solution: nucleardatapy, a free, open-source "Swiss Army Knife" for scientists. Think of it as a universal translator and a master librarian rolled into one.

Here is a breakdown of what this toolkit does, using simple analogies:

1. The Problem: A Library Without a Catalog

For years, if a scientist wanted to know how dense matter behaves inside a neutron star, they had to:

• Dig through hundreds of old research papers.
• Download data in different file formats (like trying to open a Word doc with a PDF reader).
• Re-calculate numbers because one author used "apples" and another used "oranges" to measure the same thing.
• Guess which theories were reliable.

It was like trying to cook a gourmet meal when your ingredients are in different jars, labeled in different languages, and you don't know which ones are fresh.

2. The Solution: The "Nuclear Data App"

The authors created a Python toolkit called nucleardatapy. Think of this as a centralized app store for nuclear physics.

• One-Stop Shop: Instead of hunting for data, you just open the app. It has everything: experimental measurements from labs, theoretical predictions from supercomputers, and observations from telescopes.
• Standardized Format: It takes all those messy, different formats and translates them into a single, clean language (Python). Now, comparing a theory from 1990 with a measurement from 2024 is as easy as comparing two numbers on a spreadsheet.
• The "Meta-Analysis" Magic: The toolkit doesn't just store data; it helps you mix it. It can take 50 different theories, average them out, and tell you, "Here is the most likely reality based on what we all know."

3. What's Inside the Toolkit? (The Ingredients)

The paper details four main "pantries" inside this app:

• The "Uniform Matter" Pantry: This is the recipe for the "soup" of the universe. It contains data on how protons and neutrons behave when they are packed together perfectly evenly (like in the core of a star) versus when they are packed with only neutrons (like in a neutron star). It includes both "microscopic" recipes (calculated from the bottom up) and "phenomenological" recipes (based on fitting the data).
• The "Finite Nuclei" Pantry: This is the catalog of actual atoms we can build on Earth. It has the weights, sizes, and stability of thousands of different elements. It's like a massive phonebook of every atom ever discovered, showing which ones are stable and which ones fall apart instantly.
• The "Exotic" Pantry: This section deals with "strange" matter. It includes Hypernuclei (atoms with extra, weird particles called hyperons) and data on the "crust" of neutron stars (the crunchy outer shell of a star).
• The "Astrophysics" Pantry: This connects the lab to the sky. It contains real observations from telescopes and gravitational wave detectors (like LIGO). It has the masses of heavy neutron stars and the "tidal deformability" (how squishy they are) from colliding stars.

4. The Big Discovery: Do the Stars Agree?

One of the coolest things the authors did with this toolkit was to play a game of "Cross-Check."

• Team A (Nuclear Physicists): They used the toolkit to calculate the pressure inside a neutron star based only on lab experiments and theory.
• Team B (Astrophysicists): They used gravitational wave data from colliding stars to guess the pressure.
• The Result: When they compared the two using the toolkit, they matched! The pressure predicted by nuclear physics fits perfectly with the pressure observed in the sky. It's like two detectives solving a crime from opposite sides of the world and realizing they are looking at the same suspect.

5. Why This Matters for Everyone

You might ask, "Why do I care about a Python script for nuclear physics?"

• It speeds up discovery: Scientists don't waste time formatting data; they spend time solving mysteries.
• It prevents errors: By having a community-verified database, we reduce the chance of using bad data.
• It connects the dots: It helps us understand the life cycle of stars, the origin of heavy elements (like gold), and the fundamental laws of the universe.

The Bottom Line

The nucleardatapy toolkit is a community-driven bridge. It takes the scattered, confusing pieces of the nuclear puzzle and snaps them together into a clear picture. It allows scientists to say, "We know what's happening inside a neutron star," with much more confidence than before.

The authors are essentially saying: "Here is the toolbox. It's free, it's open, and we invite everyone to add their own tools to it so we can build a better understanding of the universe together."

Technical Summary

1. Problem Statement

The field of nuclear physics and nuclear astrophysics faces a significant challenge in data accessibility and interoperability. Critical data regarding dense nuclear matter, finite nuclei, hypernuclei, and astrophysical observations (such as neutron star masses and radii) are dispersed across numerous sources, often stored in disparate formats (tables, raw data files, specific code outputs).

• Fragmentation: Researchers must manually reconstruct quantities, leading to potential inconsistencies in definitions (e.g., how symmetry energy or pressure is calculated) between different authors.
• Meta-analysis Limitations: Systematic comparisons (meta-analyses) across theoretical predictions and experimental constraints are hindered by the difficulty of aggregating large datasets.
• Lack of Standardization: There is no unified, community-driven repository that provides reconstructed quantities in a transparent, standardized format, making it difficult to assess model quality or perform Bayesian analyses efficiently.

2. Methodology

The authors developed nucleardatapy, a Python-based open-source toolkit designed to centralize, standardize, and simplify access to nuclear data. The methodology involves:

• Unified Data Structure: Collecting data from microscopic calculations (e.g., Quantum Monte Carlo, Chiral EFT), phenomenological models (e.g., Skyrme, Relativistic Mean Field), experimental measurements, and astrophysical observations into a single repository.
• Object-Oriented Design: The toolkit is organized into modular Python classes (nuda) covering specific domains:
  • matter: Uniform nuclear matter (Symmetric, Neutron, Asymmetric) and Free Fermi Gas (FFG) calculations.
  • corr: Correlation diagrams between nuclear empirical parameters (NEPs).
  • nuc: Finite nuclei data (masses, charge radii, separation energies, odd-even staggering).
  • hnuc: Hypernuclei data (single/double $\Lambda$, $\Xi$).
  • crust: Neutron star crust models.
  • eos: Equations of State (EoS) combining nucleons and leptons (electrons/muons) in $\beta$-equilibrium.
  • astro: Astrophysical observations (pulsar masses, gravitational wave constraints, NICER mass-radius data).
• Reconstruction and Standardization: The toolkit reconstructs quantities (e.g., internal energy, pressure, symmetry energy) from raw data using consistent definitions, allowing for direct comparison between different models.
• Meta-analysis Tools: It includes functions to generate "reference bands" by averaging over multiple models with Gaussian weighting, allowing users to define uncertainty ranges based on community consensus.

3. Key Contributions

• The nucleardatapy Toolkit (v1.0): A publicly available Python library (GitHub) that serves as a central hub for nuclear data. It includes:
  • Microscopic Models: Variational (FP, APR), Brueckner-Hartree-Fock (BHF), Auxiliary Field Diffusion Monte Carlo (AFDMC), Quantum Monte Carlo (QMC), Many-Body Perturbation Theory (MBPT), and Nuclear Lattice EFT (NLEFT).
  • Phenomenological Models: Extensive lists of Energy Density Functionals (Skyrme, Gogny, Relativistic Hartree/Fock).
  • Experimental Data: AME mass tables (2012–2020), charge radii (Angeli & Marinova), neutron skin thickness, ISGMR energies, and hypernuclear binding energies.
  • Astrophysical Data: Pulsar masses (e.g., PSR J0348+0432, J0740+6620), gravitational wave events (GW170817, GW190814), and NICER mass-radius measurements.
• Meta-analysis Framework: The ability to construct reference bands for the Equation of State (EoS) and nuclear empirical parameters (NEPs) by statistically combining multiple models.
• EoS Construction Tools: Methods to construct the EoS for $\beta$-equilibrated matter, including lepton contributions (electrons and muons), and to connect crust and core EoS models seamlessly.
• Documentation and Tutorials: Comprehensive documentation and Jupyter notebook tutorials to facilitate adoption by the community.

4. Key Results and Illustrations

The paper demonstrates the toolkit's capabilities through several meta-analyses:

• Reference Bands: The authors constructed uncertainty bands for the internal energy of symmetric matter (SM) and neutron matter (NM) by averaging MBPT and QMC results. Most microscopic models fall within these bands, while many phenomenological models (especially older NLRH models) deviate significantly.
• Pressure at $2\rho_{sat}$: A critical finding is the prediction of pressure at twice the saturation energy density ($\approx 155$ MeV fm$^{-3}$).
  • The toolkit shows that nuclear physics predictions (restricted to models passing the reference band) are compatible with constraints derived from gravitational wave (GW) observations (GW170817).
  • However, nuclear models tend to favor the lower end of the GW uncertainty range.
• Correlation Diagrams: The toolkit visualizes correlations between nuclear parameters, such as $K_{sat}$ (incompressibility) vs. $Q_{sat}$ and $E_{sym}$ vs. $L_{sym}$. It highlights strong correlations in certain EDF families and the tension between different experimental constraints (e.g., PREX-II vs. CREX).
• Neutron Star Properties: By combining crust and core models, the toolkit generates complete EoSs. It calculates the probability distribution for the maximum mass of non-rotating neutron stars ($M_{TOV}$), showing that current data suggests $M_{TOV}$ is likely between 2.0 and 2.5 $M_{\odot}$.
• Hypernuclei: The toolkit aggregates data for single and double $\Lambda$ hypernuclei, noting that while the NAGARA event ($^6_{\Lambda\Lambda}$He) is unambiguous, other double-$\Lambda$ candidates remain uncertain.

5. Significance

• Democratization of Data: By providing a unified, open-source interface, nucleardatapy lowers the barrier to entry for researchers wishing to perform complex meta-analyses or Bayesian inference without needing to manually curate disparate datasets.
• Model Validation: It provides a standardized benchmark (reference bands) to test new theoretical models against the current state-of-the-art consensus.
• Bridging Disciplines: The toolkit effectively links nuclear physics (finite nuclei, uniform matter) with astrophysics (neutron stars, gravitational waves), facilitating a holistic understanding of dense matter.
• Community-Driven Evolution: The authors explicitly design the toolkit to be extensible, inviting the community to contribute new data and models, ensuring the database remains current with the rapidly evolving field of nuclear astrophysics.

In conclusion, nucleardatapy represents a significant infrastructure advancement for the nuclear physics community, transforming fragmented data into a coherent, accessible, and analyzable resource that supports rigorous testing of theoretical models against experimental and observational constraints.