7: Nuclear Structure

Example number 1, Page number 259

In [3]:
#importing modules
import math
from __future__ import division

#Variable declaration
A1=165;               #mass number
A2=4;                 #mass number
R1=7.731;             #radius(fermi)

#Calculation
R2=R1*(A2/A1)**(1/3);     #radius of He(fermi)

#Result
print "radius of He is",round(R2,4),"fermi"

radius of He is 2.2375 fermi

Example number 2, Page number 259

In [5]:
#importing modules
import math
from __future__ import division

#Variable declaration
n=1.008665;           #mass of neutron(amu)
p=1.007276;           #mass of proton(amu)
alpha=4.00150;        #mass of alpha particle(amu)
m=931;                

#Calculation
deltam=2*(p+n)-alpha;
BE=deltam*m;          #binding energy(MeV)
ABE=BE/4;             #average binding energy per nucleon(MeV)

#Result
print "average binding energy per nucleon is",round(ABE,2),"MeV"

average binding energy per nucleon is 7.07 MeV

Example number 3, Page number 260

In [7]:
#importing modules
import math
from __future__ import division

#Variable declaration
n=1.008665;           #mass of neutron(amu)
Li36=6.015125;        #mass of Li(amu)
Li37=7.016004;        #mass of Li(amu)
m=931;                

#Calculation
deltam=Li36+n-Li37;   
BE=deltam*m;          #binding energy of neutron(MeV)

#Result
print "binding energy of neutron is",round(BE,2),"MeV"

binding energy of neutron is 7.25 MeV

Example number 4, Page number 260

In [8]:
#importing modules
import math
from __future__ import division

#Variable declaration
BEHe=4*7.0;        #binding energy for He
BEH=2*1.1;         #binding energy for H

#Calculation
deltaE=BEHe-(2*BEH);     #energy released(MeV)   

#Result
print "energy released is",deltaE,"MeV"

energy released is 23.6 MeV

Example number 5, Page number 260

In [11]:
#importing modules
import math
from __future__ import division

#Variable declaration
n=1.008665;           #mass of neutron(amu)
p=1.007276;           #mass of proton(amu)
BE=160.647;           #binding energy(MeV)
m=931;                

#Calculation
Mx=10*(p+n)-(BE/m);   #mass(amu)

#Result
print "mass is",round(Mx,3),"amu"
print "answer given in the book is wrong"

mass is 19.987 amu
answer given in the book is wrong

Example number 6, Page number 261

In [14]:
#importing modules
import math
from __future__ import division

#Variable declaration
n=1.008665;           #mass of neutron(amu)
Ca41=40.962278;       #mass of Ca(amu)
Ca42=41.958622;       #mass of Ca(amu)
m=931;                

#Calculation
deltam=Ca41+n-Ca42;     
BE=deltam*m;          #binding energy of neutron(MeV)

#Result
print "binding energy of neutron is",round(BE,3),"MeV"
print "answer given in the book varies due to rounding off errors"

binding energy of neutron is 11.471 MeV
answer given in the book varies due to rounding off errors