nook.sources.ripl3.densities

The RIPL-3 densities segment: level-density parameters and HFB tables.

Analytic parameter files, one row per nuclide with a D0 fit:

• level-densities-egsm.dat – Enhanced Generalised Superfluid Model.

• level-densities-bfmeff.dat – Back-shifted Fermi gas (effective).

• level-densities-ctmeff.dat – Constant temperature (effective).

• level-densities-hfm.dat – Hartree-Fock microscopic normalisation.

The unit trap: EGSM lists D0 in keV, the other three in eV. Both normalise to keV here (spacing_kev), matching the resonances segment; the file’s own numbers stay in fields.

level-densities-hfb/z*.tab are the microscopic HFB-plus-combinatorial tables: per nuclide, one block per parity, rows of excitation energy against total and spin-resolved densities. Odd-A files label spin columns J=00, J=01... too – the column index means J = i + 1/2 there.

class nook.sources.ripl3.densities.LevelDensityParams(nuclide, model, target_spin, bn_mev, spacing_kev, fields)

Bases: object

One nuclide’s fit parameters for one analytic level-density model.

nuclide: Nuclide

model: str

target_spin: float | None

bn_mev: float | None

spacing_kev: Uncertain

D0 normalised to keV regardless of the file’s native unit.

fields: dict[str, float]

Every column as named in the file (D0 left in its native unit).

class nook.sources.ripl3.densities.DensityTable(nuclide, grids)

Bases: object

A spin- and parity-resolved HFB level-density table, rho in 1/MeV.

nuclide: Nuclide

grids: dict[int, tuple[tuple[float, ...], ...]]

parity (+1/-1) -> tuple of rows; each row is (U MeV, T MeV, cumulative, rho_observed, rho_total, spin-resolved…).

rho(e_mev, two_j=None, parity=None)

Nearest-grid-point level density. Total over J unless two_j given; summed over parity unless parity given.

The total is RHOOBS (the sum of the spin columns), so rho(e) == sum over J of rho(e, two_j). RHOTOT – the 2J+1-weighted state density – is a different quantity and stays available as row[4] of grids.

nook.sources.ripl3.densities.load_level_density_params(path, nuclide, model='egsm')

The analytic-model parameters for one nuclide.

nook.sources.ripl3.densities.load_hfb_density(path, nuclide)

The microscopic HFB level-density table for one nuclide.