Chapter4-Nuclear Energetics

Ex1-pg86

In [1]:
##Chapter 4, Example 4.1, Page 86
import math
##Binding energy
mn = 1.0078250
mp = 1.0086649
M = 4.0026032 ## mass of He
MD = 2.*mn+2.*mp-M ##Mass defect
BE = MD*931.5
print'%s %.2f %s'%("\n Mass defect = ",MD," u");
print'%s %.2f %s'%("\n Nuclear binding energy = ",BE," MeV"); ## answer provided in the textbook is wrong
##Answer may vary due to round off error
 Mass defect =  0.03  u

 Nuclear binding energy =  28.30  MeV

Ex2-pg89

In [2]:
##Chapter 4, Example 4.2, Page 89
import math
##Binding energy
O15 = 15.0030654 ## atomic mass of O15 isotope
mn = 1.00866492
O16 = 15.9949146 ## atomic mass of O16 isotope
c2 = 931.5 ## C^2 in MeV
S = (O15+mn-O16)*c2
print'%s %.2f %s'%("\n Binding energy = ",S," MeV"); 
##Answer may vary due to round off error
 Binding energy =  15.66  MeV

Ex3-pg94

In [3]:
##Chapter 4, Example 4.3, Page 94
import math
## Q value of an endothermic and exothermic reaction
##Exothermic reaction
Be = 9.012182 ##Reactants
He = 4.002603 ##Reactants
C12 = 12. ##Product
n = 1.008664 ##Product
C2 = 931.5 ## C^2 in MeV
Exo1 = Be+He
Exo2 = C12+n
Dif1 = Exo1-Exo2
Q1 = Dif1*C2
print'%s %.2f %s'%("\n Q of the exothermic reaction = ",Q1," MeV"); 
##Endothermic reaction
O = 15.994915 ##Reactants
n = 1.008664 ##Reactant
C13 = 13.003354 ##Product
He = 4.002603 ##product
End1 = O+n
End2 = C13+He
Dif2 = End1-End2
Q2 = Dif2*C2
print'%s %.2f %s'%("\n Q of the exothermic reaction = ",Q2," MeV"); 
##Answer may vary due to round off error 
 Q of the exothermic reaction =  5.70  MeV

 Q of the exothermic reaction =  -2.22  MeV

Ex4-pg95

In [4]:
##Chapter 4, Example 4.4, Page 95
import math
## Q value in a reaction
MH = 1.00782503
MD = 2.01410178
me = 0.00054858
C2 = 931.5
Q = (2*MH-MD-2*me)*C2
print'%s %.2f %s'%("\n Q of the reaction = ",Q," MeV");## Answer provided in the text is wrong 
##Answer may vary due to round off error 
 Q of the reaction =  0.42  MeV

Ex5-pg96

In [5]:
##Chapter 4, Example 4.5, Page 96
import math
## Q value of the reaction
mn = 1.0086649
MB = 10.0129370
MHe = 4.0026032
MLi = 7.0160040
C2 = 931.5
Q = (mn+MB-MHe-MLi)*C2 -0.48
print'%s %.2f %s'%("\n Q of the reaction = ",Q," MeV");
##Answer may vary due to round off error 
 Q of the reaction =  2.31  MeV