{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 9: Nuclear Reactions"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example number 1, Page number 299"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Q value in nuclear reaction is -1.1898 MeV\n"
     ]
    }
   ],
   "source": [
    "#importing modules\n",
    "import math\n",
    "from __future__ import division\n",
    "\n",
    "#Variable declaration\n",
    "N=14.003073;     #mass of N\n",
    "O=16.99913;      #mass of O\n",
    "H=1.007825;      #mass of H\n",
    "He=4.002604;     #mass of He\n",
    "m=931; \n",
    "\n",
    "#Calculation\n",
    "Q=(N+He-(O+H))*m;     #Q value in nuclear reaction(MeV)    \n",
    "\n",
    "#Result\n",
    "print \"Q value in nuclear reaction is\",round(Q,4),\"MeV\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example number 2, Page number 299"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "heat generated is 6.6 *10**6 KWH\n"
     ]
    }
   ],
   "source": [
    "#importing modules\n",
    "import math\n",
    "from __future__ import division\n",
    "\n",
    "#Variable declaration\n",
    "Li=7.01600;     #mass of Li\n",
    "H=1.007825;     #mass of H\n",
    "He=4.002604;    #mass of He\n",
    "m=931; \n",
    "e=1.6*10**-19;  #charge(coulomb)\n",
    "N=6.02*10**26;  #avagadro number\n",
    "M=0.1;          #mass(kg)\n",
    "x=1000*3600;\n",
    "\n",
    "#Calculation\n",
    "Q=(Li+H-(He+He))*m*10**6*e;     #heat generated by Li(J)\n",
    "mLi=Li/N;                       #mass of Li(kg) \n",
    "H=Q*M/(x*mLi);                  #heat generated(KWH)\n",
    "\n",
    "#Result\n",
    "print \"heat generated is\",round(H*10**-6,1),\"*10**6 KWH\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example number 3, Page number 300"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Q-value for the reaction is 5.485 MeV\n",
      "kinetic energy of Zn is 0.635 MeV\n"
     ]
    }
   ],
   "source": [
    "#importing modules\n",
    "import math\n",
    "from __future__ import division\n",
    "\n",
    "#Variable declaration\n",
    "Cu=62.929599;     #mass of Cu\n",
    "H=2.014102;       #mass of H(amu)\n",
    "n=1.008665;       #mass of n(amu)\n",
    "Zn=63.929145;     #mass of Zn(amu)\n",
    "m=931; \n",
    "Kx=12;          #energy of deuterons(MeV)\n",
    "Ky=16.85;       #kinetic energy of deuterons(MeV)\n",
    "\n",
    "#Calculation\n",
    "Q=(H+Cu-n-Zn)*m;     #Q-value for the reaction(MeV)\n",
    "K=Q+Kx-Ky;           #kinetic energy of Zn(MeV)\n",
    "\n",
    "#Result\n",
    "print \"Q-value for the reaction is\",round(Q,3),\"MeV\"\n",
    "print \"kinetic energy of Zn is\",round(K,3),\"MeV\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example number 4, Page number 301"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "threshold kinetic energy is 5.378 MeV\n",
      "answer given in the book is wrong\n"
     ]
    }
   ],
   "source": [
    "#importing modules\n",
    "import math\n",
    "from __future__ import division\n",
    "\n",
    "#Variable declaration\n",
    "P=1.007825;     #mass of P(amu)\n",
    "H2=2.014102;    #mass of H2(amu)\n",
    "H3=3.016049;    #mass of H3(amu)\n",
    "m=931; \n",
    "\n",
    "#Calculation\n",
    "Q=(P+H3-(2*H2))*m;         #Q-value(MeV)\n",
    "Kth=-Q*(1+(P/H3));          #threshold kinetic energy(MeV)\n",
    "\n",
    "#Result\n",
    "print \"threshold kinetic energy is\",round(Kth,3),\"MeV\"\n",
    "print \"answer given in the book is wrong\""
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}