NuJuliet - Dev version
Many-body nuclear structure solver written in Julia [1]. Currently, the code features these solvers:
- Spherical Hartree-Fock (HF) solver, including leading order Many-body Perturbation Theory (MBPT) corrections.
- Spherical Hartree-Fock Tamm-Dancoff & Random-Phase Approximation (HF-TDA & HF-RPA) solver.
- Spherical HF-Renormalized RPA (HF-RRPA, also known as ERPA or IRPA) solver.
- Spherical HF-Bardeen-Cooper-Schrieffer (HF-BCS) solver.
- Spherical Hartree-Fock-Bogoliubov (HFB) solver, including leading order Many-body Perturbation Theory correlation energy calculation (BMBPT(2)).
- Spherical Quasiparticle-Tamm-Dancoff Approximation (QTDA) solver.
- Spherical Quasiparticle-Random-Phase Approximation (QRPA) solver.
With these one can study properties of closed-(sub)shell nuclei as well as basic properties of spherical open-shell nuclei.
The code accepts input NN & NNN interaction files in binary format generated by the NuHamil code [2].
To run the code, you have to download it, and have installed Julia [1] (v9.4 or newer should be sufficient, I suggest you do not use >v1.10. due to performance drop caused by the way compiler handles dictionaries used for interaction matrix element lookups, the lastest Long-Term Stable release is the best to use). Furthermore, the solver requires these Julia packages:
- LinearAlgebra
- DelimitedFiles
- BenchmarkTools
- Printf
- CGcoefficient
- Hungarian
- Mmap
Sample script "NuJuliet/Scripts/sample_script_run.jl" provides basic solver calls. You have to perform Julia calls of your scripts from "NuJuliet/" for proper linking. Output files are always saved under "IO/". Interaction input files for the interaction
More details how to run & how to set up the inputs can be found in pdf manual "NuJuliet_manual.pdf".
Pdf file with more detailed instructions to run the solver will be included ...
[1] Julia language: https://julialang.org/
[2] NuHamil code: https://github.com/Takayuki-Miyagi/NuHamil-public