9: Nuclear Reactions

Example number 1, Page number 299

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

#Variable declaration
N=14.003073;     #mass of N
O=16.99913;      #mass of O
H=1.007825;      #mass of H
He=4.002604;     #mass of He
m=931; 

#Calculation
Q=(N+He-(O+H))*m;     #Q value in nuclear reaction(MeV)    

#Result
print "Q value in nuclear reaction is",round(Q,4),"MeV"
Q value in nuclear reaction is -1.1898 MeV

Example number 2, Page number 299

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

#Variable declaration
Li=7.01600;     #mass of Li
H=1.007825;     #mass of H
He=4.002604;    #mass of He
m=931; 
e=1.6*10**-19;  #charge(coulomb)
N=6.02*10**26;  #avagadro number
M=0.1;          #mass(kg)
x=1000*3600;

#Calculation
Q=(Li+H-(He+He))*m*10**6*e;     #heat generated by Li(J)
mLi=Li/N;                       #mass of Li(kg) 
H=Q*M/(x*mLi);                  #heat generated(KWH)

#Result
print "heat generated is",round(H*10**-6,1),"*10**6 KWH"
heat generated is 6.6 *10**6 KWH

Example number 3, Page number 300

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

#Variable declaration
Cu=62.929599;     #mass of Cu
H=2.014102;       #mass of H(amu)
n=1.008665;       #mass of n(amu)
Zn=63.929145;     #mass of Zn(amu)
m=931; 
Kx=12;          #energy of deuterons(MeV)
Ky=16.85;       #kinetic energy of deuterons(MeV)

#Calculation
Q=(H+Cu-n-Zn)*m;     #Q-value for the reaction(MeV)
K=Q+Kx-Ky;           #kinetic energy of Zn(MeV)

#Result
print "Q-value for the reaction is",round(Q,3),"MeV"
print "kinetic energy of Zn is",round(K,3),"MeV"
Q-value for the reaction is 5.485 MeV
kinetic energy of Zn is 0.635 MeV

Example number 4, Page number 301

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

#Variable declaration
P=1.007825;     #mass of P(amu)
H2=2.014102;    #mass of H2(amu)
H3=3.016049;    #mass of H3(amu)
m=931; 

#Calculation
Q=(P+H3-(2*H2))*m;         #Q-value(MeV)
Kth=-Q*(1+(P/H3));          #threshold kinetic energy(MeV)

#Result
print "threshold kinetic energy is",round(Kth,3),"MeV"
print "answer given in the book is wrong"
threshold kinetic energy is 5.378 MeV
answer given in the book is wrong