{
 "metadata": {
  "name": "Chapter13",
  "signature": "sha256:3612dc1ad0c2617d4c1e1da21ea7791334391562bfe93ddc9d0f98037d4fe15d"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 13:Nuclear Reaction and Application"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 13.1, Page 417"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#initiation of variable\n",
      "v=1*1.0*10**-6.0*10**2; p=7.9; m=p*v;Na=6.023*10**23     #given values and various constants in suitable units\n",
      "M=56.0;N=m*Na/M;                                 #number of atoms\n",
      "i=3.0*10**-6;\n",
      "q=1.6*10**-19;\n",
      "\n",
      "#calculation\n",
      "Io=i/q;                 #intensity\n",
      "s=0.6*10**-24;S=1;             #given values in suitable units\n",
      "R=N*s*Io/S;                     #rate of neutrons\n",
      "\n",
      "#result\n",
      "print\"The rate of neutrons emitted from the target in particles per second is %.1e\" %round(R,3);\n",
      "print\"Slight difference in answer due to approximation error\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The rate of neutrons emitted from the target in particles per second is 9.6e+07\n",
        "slight difference due to approximation error\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 13.2, Page 419"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#initiation of variable\n",
      "A=197.0; m=30*10**-3;phi=3.0*10**12;    #given values and various constants taken in suitable units\n",
      "Ar=99.0*10**-24; Na=6.023*10**23\n",
      "\n",
      "#calculation\n",
      "R=(phi*Na*Ar*m/A);        #rate or production of gold\n",
      "t=2.7*24*60         # time of decay\n",
      "Act=R*(0.693/t);     #activity /sec\n",
      "ActCi=Act/(3.7*10**4);     # in terms of curie(Ci)\n",
      "\n",
      "#result\n",
      "print\"The activity is found out to be %.1e\" %round(Act,3),\"/sec i.e \" ,round(ActCi,3),\"muCi\"\n",
      "print\"Slight difference in answer due to approximation error\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The activity is found out to be 4.9e+06 /sec i.e  131.229 muCi\n",
        "slight difference due to approximation error\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 13.3, Page 423"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#initiation of variable\n",
      "from math import exp\n",
      "v=1.5*1.5*2.5*(10**-6)*10**2; #volume in cm3\n",
      "p=8.9;         #density in g/cm3\n",
      "m=p*v;Na=6.023*10**23 #mass and Avagadro's number\n",
      "M=58.9;        #Given values\n",
      "\n",
      "#calculation\n",
      "N=m*Na/M;\n",
      "i=12*10**-6;       #thickness of beam\n",
      "q=1.6*10**-19;\n",
      "Io=i/(2*q);      #intensity\n",
      "s=0.64*10**-24;   #Given values\n",
      "S=1.5*1.5;\n",
      "R=N*s*Io/S;       #rate of production of 61Cu\n",
      "\n",
      "#result\n",
      "print \"The rate of neutrons emitted from the target in particles/second is %.1e\" %round(R,3);\n",
      "\n",
      "#part b\n",
      "act=R*(1-(exp((0.693)*(-2/3.41))));          #activity\n",
      "\n",
      "#result\n",
      "print\"The activity after 2.0h in /sec is %.1e\" %round(act,3),\"=4.9mCi\";\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The rate of neutrons emitted from the target in particles/second is 5.5e+08\n",
        "The activity after 2.0h in /sec is 1.8e+08 =4.9mCi\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 13.4, Page 425"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#initiation of variable\n",
      "m2H=2.014102;    #mass of various particles\n",
      "mn=1.008665;m63Cu=62.929599;\n",
      "m64Zn=63.929145;c2=931.5;     #c^2=931.5 MeV\n",
      "Q=(m2H+m63Cu-mn-m64Zn)*c2;       #Q of the reaction\n",
      "\n",
      "#result\n",
      "print\"The value of Q is in MeV\",round(Q,3);\n",
      "\n",
      "\n",
      "#part b\n",
      "Kx=12.00;Ky=16.85;\n",
      "Ky=Q+Kx-Ky              #kinetic energy of 64Zn\n",
      "\n",
      "#result\n",
      "print\"The value of Ky was found out to be in MeV\",round(Ky,3);"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The value of Q is in MeV 5.487\n",
        "The value of Ky was found out to be in MeV 0.637\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 13.5, Page 425"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#initiation of variable\n",
      "mp=1.007825;m3H=3.016049;    #mass of the particle\n",
      "m2H=2.014102;c2=931.5;      #constant\n",
      "Q=(mp+m3H-(2*m2H))*c2;     #Q of the reaction\n",
      "\n",
      "#result\n",
      "print\"The value of q was found out to be in MeV\",round(Q,3);\n",
      "\n",
      "#partb\n",
      "Kth1= -Q*(1+(mp/m3H));        #threshold energy of kinetic energy\n",
      "Kth2=-Q*(1+(m3H/mp));        #threshold kinetic energy in case2\n",
      "\n",
      "#result\n",
      "print\"The threshold kinetic energy in case-1 in MeV\",round(Kth1,3);\n",
      "print\"The threshold kinetic energy in case-2 in MeV\",round(Kth2,3);"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The value of q was found out to be in MeV -4.033\n",
        "The threshold kinetic energy in case-1 in MeV 5.381\n",
        "The threshold kinetic energy in case-2 in MeV 16.104\n"
       ]
      }
     ],
     "prompt_number": 6
    }
   ],
   "metadata": {}
  }
 ]
}