This is a simple code for calculating the valence bands of a nucleus, according to the Liquid Droplet Model in Nuclear Physics. We also calculate the binding energy of the nuclides and plot for certain values.
I developed this code with the help of AnthonyAposta for our Nuclear Physics Course at our home university (Universidade Federal de Santa Catarina - UFSC), if you like what we are doing, you should also check his profile.
We took our raw data (the isotopes descriptions) from the Nuclear Energy Agency website. The archive is originally a .txt (which is a poor choice of structured data). So we filtered the information we found relevant and built an .csv archive, under the name of isotopos.csv in this repository.
The binding energy can be calculated under the Bond.py code. If one wishes to calculate the plot for the valence band, the function is under the name "gota_liquida".
General information on specific Isotopes can be calculated in the main script of this repository "ModNucl.py". There is also an option for writing it directly as a LateX table, if necessary.
Binding Energy with deleted parameters:
Liquid Droplet Plots: With R=1,beta=1,gamma=1
Liquid Droplet Plots: With R=1,beta=4,gamma=30
Liquid Droplet Plots: With R=1,beta=8,gamma=60