nuclear_mass

nuclear_mass

This project allows users to work with their own nuclear masses in NucNet Tools/Projects XML files. The default mass formula is that of Myers-Swiateck.

Installing the project.

If you have not previously done so, install the necessary [libraries]. Next, if you have not previously done so, create a projects directory:

** mkdir projects **

Change into the projects directory:

** cd projects **

Now you can checkout the project:

** svn checkout svn://svn.code.sf.net/p/nucnet-projects/code/trunk/projects/nuclear_mass**

Compiling the executables.

Change into the project directory:

** cd nuclear_mass **

and compile the codes:

** ./project_make **

The project compiled successfully if you see Make succeeded. Otherwise, please submit a Ticket or make a post to the Discussion board (see links above).

Compute nuclear mass terms from the Myers-Swiateck mass formula.

You will need a nuclear data file to specify the nuclei. Download the JINA reaclib database by typing

** make data**

Now print out the mass terms by typing

** ./mass_terms ../nucnet-tools-code/data_pub/my_net.xml **

This will print out the species name, Z, A, volume energy, surface energy, Coulomb energy, pairing energy, and shell energy. To select out only certain nuclei, use an XPath expression. For example, to print out only terms for n, p, and iron isotopes, type

** ./mass_terms ../nucnet-tools-code/data_pub/my_net.xml "[a = 1 or z = 26]" **

Note that you must include the neutron and proton (a = 1) because the code needs the mass for those species.

Update the nuclear data in an XML file with values from the Myers-Swiatecki mass formula.

To update the nuclear data with that from the Myers-Swiateck mass formula, type

** ./update_net_xml ../nucnet-tools-code/data_pub/my_net.xml new.xml **

The newly created file new.xml will have mass excesses for all species other than n and p updated with the Myers-Swiatecki values. It will also have simplified partition function data.

To select a subset of nuclei to be update, use an XPath expression. For example, to update only iron, cobalt, and nickel isotopes, type

** ./update_net_xml ../nucnet-tools-code/data_pub/my_net.xml new.xml "[z >= 26 and z <= 28]" **

Compute the mass term energies per nucleon in a network calculation.

You can compute the mass term energies per nucleon from a network calculation. To do so, you need the output XML. As an example, one might use my_alpha_output.xml, as computed in this blog post. Assuming you have moved my_alpha_output.xml to the nuclear_mass directory, type

** ./compute_mass_terms_in_zones my_alpha_output.xml my_new_alpha_output.xml **

The newly created file my_new_alpha_output.xml has the energies per nucleon added as optional properties. To extract them, use the NucNet Tools print_properties code. Type

** cd ../nucnet-tools-code/examples/analysis **

** make all_analysis **

** cd ../../../nuclear_mass **

** ../nucnet-tools-code/examples/analysis/print_properties my_new_alpha_output.xml time "Volume term" "Surface term" "Coulomb term" "Pairing term" "Shell term" **

This prints out, for each timestep in the calculation, the timestep number, the time in seconds, and then each of the energy terms per nucleon.

Project Members:

• Bradley S. Meyer (admin)