Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

1 files changed, 174 insertions, 0 deletions

diff --git a/Modern_Physics_by_R_A_Serway/7-Tunnelling_phenomena.ipynb b/Modern_Physics_by_R_A_Serway/7-Tunnelling_phenomena.ipynb
new file mode 100644
index 0000000..05a54f5
--- /dev/null
+++ b/Modern_Physics_by_R_A_Serway/7-Tunnelling_phenomena.ipynb
@@ -0,0 +1,174 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 7: Tunnelling phenomena"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.1: Transmission_coefficient_for_an_oxide_layer.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"// Scilab code Ex7.1: Pg 235 (2005)\n",
+"clc; clear;\n",
+"c = 3e+08;    // Velocity of light, m/s\n",
+"m_e = 511e+03/(c^2);    // Mass of electron, eV\n",
+"U = 3.00;    // Ground state energy neglecting E, eV\n",
+"h_cross = (1.973e+03)/c;    // Reduced planck's constant, eV\n",
+"alpha = sqrt(2*m_e*U)/h_cross;\n",
+"L = 50;     // Thickness of the layer, angstrom\n",
+"T = 1/(1+1/4*10^2/(7*3)*sinh(alpha*L)^2);\n",
+"printf('\nThe transmission coefficient for the layer thickness of');\n",
+"printf('\n%2d angstrom = %5.3e', L, T);\n",
+"L = 10;    // // Thickness of the layer, angstrom\n",
+"T = 1/(1+1/4*10^2/(7*3)*sinh(alpha*L)^2);\n",
+"printf('\n%2d angstrom = %5.3e', L, T);\n",
+"\n",
+"// Result\n",
+"// The transmission coefficient for the layer thickness of\n",
+"// 50 angstrom = 9.628e-39\n",
+"// 10 angstrom = 6.573e-08 "
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.2: Tunnelling_current_through_an_oxide_layer.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"// Scilab code Ex7.2: Pg 236 (2005)\n",
+"clc; clear;\n",
+"e = 1.60e-19;    // Electrc charge, C\n",
+"i = 1.00e-03;      // Electron current, A\n",
+"N = i/e;       // Electrons per second\n",
+"T = 0.657e-07;     // Fraction of electrons transmitted\n",
+"T_e = N*T;    // Number of electrons transmitted per second\n",
+"T_i = T_e*e;    // Transmitted current, A\n",
+"printf('\nThe transmitted current through the oxide layer = %4.1f pA', T_i*1e+12);\n",
+"\n",
+"// Result\n",
+"// The transmitted current through the oxide layer = 65.7 pA "
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.5: Tunnelling_in_a_parallel_plate_capacitor.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"// Scilab code Ex7.5: Pg 241 (2005)\n",
+"clc; clear;\n",
+"epsilon_c = 5.5e+10;    // Characteristic field strength, V/m\n",
+"epsilon = 1.0e+09;    // Electric field, V/m\n",
+"f = 1.0e+30;     // Collision frequency, s(-1)cm(-2)\n",
+"lamda = f*exp(-epsilon_c/epsilon);   // Electron emission rate, electrons/sec\n",
+"e = 1.60e-19;    // Electrc charge, C\n",
+"I = lamda*e;    // Tunelling current, A\n",
+"printf('\nTunelling current in parallel plate capacitor = %4.2f pA', I/1e-12);\n",
+"printf('\n');\n",
+"\n",
+"// Result\n",
+"// Tunelling current in parallel plate capacitor = 0.21 pA"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.6: Estimating_halflives_of_Thorium_and_Polonium.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"// Scilab code Ex7.6: Pg 244 (2005)\n",
+"clc; clear;\n",
+"Z_T = 88;     // Atomic number of daughter nucleus\n",
+"E_T = 4.05e+06;     // Energy of ejected alphas, eV\n",
+"R = 9.00e-15;      // Nuclear radius, m\n",
+"r_o = 7.25e-15;     // Bohr radius, m\n",
+"E_o = 0.0993e+06;     // Energy analogous to the Rydberg in Atomic Physics \n",
+"T_T = exp(-4*%pi*Z_T*sqrt(E_o/E_T) + 8*sqrt((Z_T*R)/r_o));    // Transmission factor in case of Thorium\n",
+"f = 1e+21;    // Frequency of collisions, Hz\n",
+"lamda_T = f*T_T;    // Decay rate in case of Thorium, s^(-1)\n",
+"t_T = 0.693/lamda_T;     // Half-life time of Thorium, s\n",
+"Z_P = 82;    // Atomic number of daughter nucleus\n",
+"E_P = 8.95e+06;     // Energy of ejected alphas, eV\n",
+"R = 9.00e-15;     // Nuclear radius, m\n",
+"r_o = 7.25e-15;     // Bohr radius, m\n",
+"E_o = 0.0993e+06;      // Energy unit, eV\n",
+"T_P = exp(-4*%pi*Z_P*sqrt(E_o/E_P) + 8*sqrt((Z_P*R)/r_o));    // Transmission factor in case of Polonium\n",
+"f = 1e+21;    // Frequency of collisions, Hz\n",
+"lamda_P = f*T_P;      // Decay rate in case of Thorium, s^(-1)\n",
+"t_P = 0.693/lamda_P;    // Half-life time of Polonium, s\n",
+"\n",
+"printf('\nHalf-life time of Thorium = %3.1e s = %3.1e yrs', t_T, t_T/(365*24*60*60));\n",
+"printf('\nHalf-life time of Polonium = %3.1e s', t_P);\n",
+"\n",
+"// Result\n",
+"// Half-life time of Thorium = 5.3e+17 s = 1.7e+10 yrs\n",
+"// Half-life time of Polonium = 8.4e-10 s "
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}