Weizäcker's mass formula

Chart files for Isotope Explorer

Weizäcker's mass formula

File name: $kranem.chr

Weizäcker's mass formula compared with experimental data. The parameters are from Krane, Nuclear Physics. The chart is color-coded according to the difference between the calculated and experimental binding energy per nucleon.

Requirements

None.

Build section


// Weizacker's mass formula with parameters from Krane, Nuclear Physics

// define A, Z and N:

Def(Z=nuclZ())  

Def(A=nuclA())  

Def(N=nuclN())  

// define constants:

Def(av=15500)   //volume term constant

Def(as=16800)   //surface term constant

Def(ac=720)     //coulomb term constant

Def(aa=23000)   //asymmetry term constant

Def(ap=34)      //pairing term constant

// pairing term:

Def(Delta=ap/A^(3/4)*If(Even(A),If(Even(Z),1,-1),0)) 

Def(beE=be/A)   //experimental binding energy per nucleon

// calculated binding energy per nucleon:

Def(beW=(av-as*A^(-1/3)-ac*Z*(Z-1)/A^(4/3)-aa*(A-2*Z)^2/A^2+Delta/A))

// difference between experimental and calculated value:

Def(diff=beW-beE) 

ColorAccordingTo(diff)  // color scale on difference

// load data:

FindOut(beW)

FindOut(beE)

FindOut(diff)

FindOut(Delta)