See class o2scl::nucleus_rmf .
This example uses the NL3 EOS as implemented in o2scl::eos_had_rmf and computes the structure of using o2scl::nucleus_rmf. The results are stored in o2scl::table_units objects and output to HDF files. It also computes the energy of three unoccupied neutron and proton levels.
Typical output:
Proton orbitals: Index State Energy Nodes Deg. 0 ^{1}s_{1/2} -4.848748e+01 0 2 1 ^{1}p_{3/2} -4.288159e+01 0 4 2 ^{1}p_{1/2} -4.220057e+01 0 2 3 ^{1}d_{5/2} -3.586023e+01 0 6 4 ^{1}d_{3/2} -3.429919e+01 0 4 5 ^{2}s_{1/2} -3.054878e+01 1 2 6 ^{1}f_{7/2} -2.786126e+01 0 8 7 ^{1}f_{5/2} -2.508489e+01 0 6 8 ^{2}p_{3/2} -2.058053e+01 1 4 9 ^{2}p_{1/2} -1.951978e+01 1 2 10 ^{1}g_{9/2} -1.922203e+01 0 10 11 ^{1}g_{7/2} -1.501610e+01 0 8 12 ^{2}d_{5/2} -1.081390e+01 1 6 13 ^{1}h_{11/2} -1.020303e+01 0 12 14 ^{2}d_{3/2} -9.184272e+00 1 4 15 ^{3}s_{1/2} -8.054054e+00 2 2 Neutron orbitals: Index State Energy Nodes Deg. 16 ^{1}s_{1/2} -5.903999e+01 0 2 17 ^{1}p_{3/2} -5.283626e+01 0 4 18 ^{1}p_{1/2} -5.224434e+01 0 2 19 ^{1}d_{5/2} -4.536117e+01 0 6 20 ^{1}d_{3/2} -4.394391e+01 0 4 21 ^{2}s_{1/2} -4.083662e+01 1 2 22 ^{1}f_{7/2} -3.699136e+01 0 8 23 ^{1}f_{5/2} -3.438846e+01 0 6 24 ^{2}p_{3/2} -3.046551e+01 1 4 25 ^{2}p_{1/2} -2.940207e+01 1 2 26 ^{1}g_{9/2} -2.806148e+01 0 10 27 ^{1}g_{7/2} -2.402928e+01 0 8 28 ^{2}d_{5/2} -2.046789e+01 1 6 29 ^{1}h_{11/2} -1.884691e+01 0 12 30 ^{2}d_{3/2} -1.881771e+01 1 4 31 ^{3}s_{1/2} -1.807312e+01 2 2 32 ^{1}h_{9/2} -1.335125e+01 0 10 33 ^{2}f_{7/2} -1.100552e+01 1 8 34 ^{1}i_{13/2} -9.576804e+00 0 14 35 ^{2}f_{5/2} -8.991761e+00 1 6 36 ^{3}p_{3/2} -8.290647e+00 2 4 37 ^{3}p_{1/2} -7.515362e+00 2 2 Proton orbitals: Index State Energy Nodes Deg. 0 ^{3}s_{1/2} -8.057079e+00 2 2 1 ^{1}h_{9/2} -4.550103e+00 0 10 2 ^{2}f_{7/2} -1.365592e+00 1 8 Neutron orbitals: Index State Energy Nodes Deg. 3 ^{1}i_{11/2} -2.912219e+00 0 12 4 ^{1}j_{15/2} -5.639563e-01 0 16 5 ^{3}d_{5/2} -1.097961e+00 2 6 Proton RMS radius : 5.460268e+00 fm. Neutron RMS radius : 5.740404e+00 fm. Total energy per nucleon: -7.842549e+00 MeV. 1 tests performed. All tests passed.
Documentation generated with Doxygen. Provided under the
GNU Free Documentation License (see License Information).