masstable 0.2.1

The Nuclear Mass Table Toolkit provides utilities to work with nuclear mass tables. At the moment the following tables are supported:

• AME2003: G. Audi, H. Wapstra, C. Thibault, Nucl. Phys. A 729 (2003) 337
• AME2003all: Same as above but including interpolated(“#”) values
• AME2012: G. Audi et al, Chinese Physics C. 36, No. 12(2012)
• AME2012all: Same as above but including interpolated(“#”) values
• AME1995: Audi, G., and A. H. Wapstra, Nuclear Physics A 595 (4) (1995): 409-480
• AME1995all: Same as above but including interpolated(“#”) values
• DUZU: J. Duflo, A.P. Zuker, Phys. Rev. C 52 (1995)
• FRDM95: Moller, P. et al., At. Data and Nuc. Data Tables 59 (1995) 185
• KTUY05: H. Koura, T.Tachibana, M. Uno, M. Yamada, Progr. Theor. Phys. 113 (2005) 305
• ETFSI12: Y. Aboussir et al., At. Data Nucl. Data Tables 61 (1995) 127
• MAJA88: Masson, P. J., and J. Jänecke, Atomic Data and Nuclear Data Tables 39.2 (1988): 273-280.
• GK88: Jänecke, J., and P. J. Masson, Atomic Data and Nuclear Data Tables 39.2 (1988): 265-271.
• HFB14: S. Goriely, M. Samyn, J.M. Pearson, Phys. Rev. C 75 064312 (2007)
• HFB26: S. Goriely, N. Chamel, J.M. Pearson, Phys. Rev. C 88 024308 (2013)
• WS32010: N. Wang, Z. Liang, M. Liu, X. Wu, Phys. Rev. C 82 044304 (2010)
• WS32011: M. Liu, N. Wang, Y. Deng, X. Wu, Phys. Rev. C 84 014333 (2011)

See the full documentation at http://elyase.github.io/masstable/.

Usage:

• Print first 5 elements from Audi 2003:

>>> from masstable import Table
>>> Table('AME2003').head()
Z  N
0  1     8.07132
1  0     7.28897
   1    13.13570
   2    14.94980
   3    25.90150

• Calculate the root mean squared error of Moller, et al. Atomic Data and Nuclear Data Tables, 59 (1995), 185-351.

>>> Table('FRDM95').rmse(relative_to='AME2003')
0.890859326191

• Calculate 2 neutron separation energies for even-even nuclei:

>>> table = Table('AME2012').even_even.s2n
Z  N
2  2           NaN
   4      0.975454
   6      2.125034
   8     -1.417666
4  2           NaN
       ...

• Select nuclei with Z,N > 28:

>>> condition = lambda Z,N: Z > 28 and N > 28
>>> table.select(condition)
30  30    28.016334
    32    23.136434
    34    20.978934
    36    19.037934
    38    17.250334
    40    15.700534
       ...

• Plot binding energies per nucleon:

>>> t = Table('AME2012')
>>> (t.binding_energy/t.A).plot()

Install

Just do:

pip install masstable

Requirements

• python >= 2.7
• pandas >= 0.11

Credits

Yaser Martinez