Jlm

[beykyle]

• Add LDA and ILDA folding capabilities
• add JLM microscopic BHF nuclear matter self energy parameterization
• add JLMB energy renormalization
• add JLM and JLMB functions for computing finite nucleus self-energy given nuclear matter density
• add nuclear matter density tables and parameterizations from TALYS
• add extensize example demoing new capabilities

[beykyle]

@copilot Please review this pull request (You don't have to review the contents of @src/data/densities/, this is just data, just the source code, docs and example updates).

[Copilot]

@copilot Please review this pull request (You don't have to review the contents of @src/data/densities/, this is just data, just the source code, docs and example updates).

Reviewed the PR changes (excluding src/data/densities/). I found one medium-severity issue:

• src/jitr/folding/jlm.py (W0 around line 150, W1 around line 266): both use terms with 1 / (E - E_F) (or squared), so E == E_F can trigger divide-by-zero runtime warnings and unstable edge-case behavior near the Fermi surface. Suggest guarding with a small epsilon or safe np.where handling.

No other high-signal blockers stood out in the source/docs/examples changes.

[beykyle]

@copilot Please make a commit resolving the potential divide by 0 issue.

[Copilot]

@copilot Please make a commit resolving the potential divide by 0 issue.

Implemented in commit 2f03210.

W0 and W1 now use epsilon-stabilized denominators around E == E_F to prevent divide-by-zero warnings while keeping behavior finite near the Fermi surface. I also added regression coverage in tests/test_jlm_optical.py for both exact and near-Fermi energies.

[beykyle]

TODO:

• finish example_jlm.ipynb, including cross section calculation