(*^ ::[ Information = "This is a Mathematica Notebook file. It contains ASCII text, and can be transferred by email, ftp, or other text-file transfer utility. It should be read or edited using a copy of Mathematica or MathReader. If you received this as email, use your mail application or copy/paste to save everything from the line containing (*^ down to the line containing ^*) into a plain text file. On some systems you may have to give the file a name ending with ".ma" to allow Mathematica to recognize it as a Notebook. The line below identifies what version of Mathematica created this file, but it can be opened using any other version as well."; FrontEndVersion = "NeXT Mathematica Notebook Front End Version 2.2"; NeXTStandardFontEncoding; fontset = title, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e8, 24, "Times"; ; fontset = subtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, bold, L1, e6, 18, "Times"; ; fontset = subsubtitle, inactive, noPageBreakBelow, noPageBreakInGroup, nohscroll, preserveAspect, groupLikeTitle, center, M7, italic, L1, e6, 14, "Times"; ; fontset = section, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, grayBox, M22, bold, L1, a20, 18, "Times"; ; fontset = subsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, blackBox, M19, bold, L1, a15, 14, "Times"; ; fontset = subsubsection, inactive, noPageBreakBelow, nohscroll, preserveAspect, groupLikeSection, whiteBox, M18, bold, L1, a12, 12, "Times"; ; fontset = text, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = smalltext, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 10, "Times"; ; fontset = input, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeInput, M42, N23, bold, L1, 12, "Courier"; ; fontset = output, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L-5, 12, "Courier"; ; fontset = message, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = print, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = info, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeOutput, M42, N23, L1, 12, "Courier"; ; fontset = postscript, PostScript, formatAsPostScript, output, inactive, noPageBreakInGroup, nowordwrap, preserveAspect, groupLikeGraphics, M7, l34, w282, h287, L1, 12, "Courier"; ; fontset = name, inactive, noPageBreakInGroup, nohscroll, preserveAspect, M7, italic, B65535, L1, 10, "Times"; ; fontset = header, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, italic, L1, 12, "Times"; ; fontset = leftheader, 12; fontset = footer, inactive, nohscroll, noKeepOnOnePage, preserveAspect, center, M7, italic, L1, 12, "Times"; ; fontset = leftfooter, 12; fontset = help, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = clipboard, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = completions, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12, "Courier"; ; fontset = special1, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special2, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special3, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special4, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; fontset = special5, inactive, nohscroll, noKeepOnOnePage, preserveAspect, M7, L1, 12; paletteColors = 128; automaticGrouping; remoteKernel = "Local"; ] :[font = title; inactive; preserveAspect; startGroup] Alpha Decay :[font = subtitle; inactive; preserveAspect; plain; fontName = "Times"] by Gerard P. Gilfoyle Department of Physics University of Richmond Richmond, VA 23173 gilfoyle@urvax.urich.edu :[font = subsubtitle; inactive; preserveAspect; plain; bold; italic; fontName = "Times"] Published in: Mathematica in Education and Research Vol.4 No.1 ;[s] 3:0,0;14,1;52,2;62,-1; 3:1,12,9,Times,1,14,0,0,0;1,13,10,Times,3,14,0,0,0;1,12,9,Times,1,14,0,0,0; :[font = subsubtitle; inactive; preserveAspect] Copyright 1995 TELOS/Springer-Verlag :[font = section; inactive; Cclosed; preserveAspect; startGroup] 1. Introduction :[font = text; inactive; preserveAspect; endGroup] One of the earliest and most convincing verifications of the validity of quantum mechanics was the resolution of the paradox of alpha decay observed in heavy nuclei. Gamow, Condon, and Gurney treated the problem as a barrier penetration problem and were able to provide a reasonable explanation of the decay rates for a broad range of alpha-emitting nuclei. In this exercise you will compile some of the existing data on alpha decay and make a comparison with a theoretical calculation based on the transfer matrix method. :[font = section; inactive; Cclosed; preserveAspect; startGroup] 2. Constants, Commands, and Functions :[font = text; inactive; preserveAspect] In this Section the ingredients necessary for the calculation of the transmission coefficient for the alpha potential barrier are generated. The commands in each Section should be executed since they will be used in the calculation of the transmission coefficient and they demonstrate the plotting of data with Mathematica. ;[s] 3:0,0;311,1;322,2;324,-1; 3:1,11,8,Times,0,12,0,0,0;1,10,8,Times,2,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] 2.1 Constants: :[font = text; inactive; preserveAspect] The mass and charge of the emitted particle are defined by "Aout" and "Zout" in units of atomic mass units electron charges respectively. The mass and charge of the parent nucleus is defined by "Anuc" and "Znuc". The other quantities below are necessary constants. To define the constants place the pointer in the cell, click, and then hit "Enter". The emitted particle is an alpha particle. :[font = text; inactive; preserveAspect] Note: Mathematica will warn you about "possible" spelling errors in some of the definitions below because of the similarity of some of the names used. You can usually, safely ignore this warning. List of constants. Aout - emitted particle mass Zout - emitted particle charge Anuc - atomic number of parent nucleus Mout - mass of parent Znuc - charge of parent R0 - radius of parent Rmax - classical turning point of emitted particle Ndiv - Number of segments step - width of an individual barrier Ve - array of (radius, potential energy) :[font = input; preserveAspect; endGroup] Aout = 4.0; Zout = 2.0; Anuc = 212; Znuc = 84; Ndiv = 10; hbarc = 197.0; (* h_bar times c *) e2 = hbarc/137; (* electron charge squared *) :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] 2.2 Necessary commands: :[font = text; inactive; preserveAspect] Some of the commands needed to calculate the transmission coefficient for a given emitted particle energy are demonstrated here. The technique for generating an array of numbers is shown first. You must determine what values R0, Rmax, and step should have in terms of the constants defined above. The values given below are for demonstration purposes only. :[font = input; preserveAspect] R0 = 8.0; Rmax = 30.0; step = 1.0; Ve = Table[{r, Zout*Znuc*e2/r}, {r, R0, Rmax-step,step}]; :[font = text; inactive; preserveAspect] Put zeros at each end of the array since we are treating the potential as if it is zero in those regions. :[font = input; preserveAspect] PrependTo[Ve, {R0-step, 0}]; AppendTo[Ve, {Rmax, 0}]; :[font = text; inactive; preserveAspect] Initialize the transfer matrix to the unit matrix. :[font = input; preserveAspect] test = { {1,0}, {0,1} }; :[font = text; inactive; preserveAspect] Use the Do command to loop over the segments and calculate the transfer matrix by multiplying successive matrices together. The commands for extracting elements from the array Ve and multiplying matrices are also demonstrated here within the loop. ;[s] 3:0,0;8,1;10,2;248,-1; 3:1,11,8,Times,0,12,0,0,0;1,10,8,Times,1,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = input; preserveAspect; endGroup] Do[ Vtest = Ve[[i,2]]; mat1 = {{ Vtest, 0}, { 0, Vtest} }; mat2 = {{-Vtest, 0}, { 0, Vtest} }; test = test.mat1.mat2;, {i,1,Ndiv+1} ]; :[font = subsection; inactive; Cclosed; preserveAspect; startGroup] 2.3 Plotting commands :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] 2.3.1 Plotting data. :[font = text; inactive; preserveAspect] Mathematica can be used plots of data with the following commands that define an array of x and y ordered pairs ("Stuff") that is then plotted with the command "ListPlot" as demonstrated below. ;[s] 2:0,0;11,1;196,-1; 2:1,10,8,Times,2,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = input; Cclosed; preserveAspect; startGroup] Stuff= { {1,2}, {2,4}, {3,7}, {4,6} } ListPlot[Stuff, PlotStyle->PointSize[0.025], AxesLabel->{FontForm["x",{"Times-Italic",12}], FontForm["y",{"Times-Italic",12}]}] :[font = output; output; inactive; preserveAspect] {{1, 2}, {2, 4}, {3, 7}, {4, 6}} ;[o] {{1, 2}, {2, 4}, {3, 7}, {4, 6}} :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 34; pictureWidth = 239; pictureHeight = 147] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Helvetica findfont 7 scalefont setfont % Scaling calculations -0.293651 0.31746 -0.220726 0.117721 [ [(1.5)] .18254 .01472 0 2 Msboxa [(2)] .34127 .01472 0 2 Msboxa [(2.5)] .5 .01472 0 2 Msboxa [(3)] .65873 .01472 0 2 Msboxa [(3.5)] .81746 .01472 0 2 Msboxa [(4)] .97619 .01472 0 2 Msboxa p /Times-Italic findfont 12 scalefont setfont [(x)] 1.025 .01472 -1 0 Msboxa P [(3)] .01131 .13244 1 0 Msboxa [(4)] .01131 .25016 1 0 Msboxa [(5)] .01131 .36788 1 0 Msboxa [(6)] .01131 .4856 1 0 Msboxa [(7)] .01131 .60332 1 0 Msboxa p /Times-Italic findfont 12 scalefont setfont [(y)] .02381 .61803 0 -4 Msboxa P [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .18254 .01472 m .18254 .02097 L s P [(1.5)] .18254 .01472 0 2 Mshowa p .002 w .34127 .01472 m .34127 .02097 L s P [(2)] .34127 .01472 0 2 Mshowa p .002 w .5 .01472 m .5 .02097 L s P [(2.5)] .5 .01472 0 2 Mshowa p .002 w .65873 .01472 m .65873 .02097 L s P [(3)] .65873 .01472 0 2 Mshowa p .002 w .81746 .01472 m .81746 .02097 L s P [(3.5)] .81746 .01472 0 2 Mshowa p .002 w .97619 .01472 m .97619 .02097 L s P [(4)] .97619 .01472 0 2 Mshowa p .001 w .05556 .01472 m .05556 .01847 L s P p .001 w .0873 .01472 m .0873 .01847 L s P p .001 w .11905 .01472 m .11905 .01847 L s P p .001 w .15079 .01472 m .15079 .01847 L s P p .001 w .21429 .01472 m .21429 .01847 L s P p .001 w .24603 .01472 m .24603 .01847 L s P p .001 w .27778 .01472 m .27778 .01847 L s P p .001 w .30952 .01472 m .30952 .01847 L s P p .001 w .37302 .01472 m .37302 .01847 L s P p .001 w .40476 .01472 m .40476 .01847 L s P p .001 w .43651 .01472 m .43651 .01847 L s P p .001 w .46825 .01472 m .46825 .01847 L s P p .001 w .53175 .01472 m .53175 .01847 L s P p .001 w .56349 .01472 m .56349 .01847 L s P p .001 w .59524 .01472 m .59524 .01847 L s P p .001 w .62698 .01472 m .62698 .01847 L s P p .001 w .69048 .01472 m .69048 .01847 L s P p .001 w .72222 .01472 m .72222 .01847 L s P p .001 w .75397 .01472 m .75397 .01847 L s P p .001 w .78571 .01472 m .78571 .01847 L s P p .001 w .84921 .01472 m .84921 .01847 L s P p .001 w .88095 .01472 m .88095 .01847 L s P p .001 w .9127 .01472 m .9127 .01847 L s P p .001 w .94444 .01472 m .94444 .01847 L s P p /Times-Italic findfont 12 scalefont setfont [(x)] 1.025 .01472 -1 0 Mshowa P p .002 w 0 .01472 m 1 .01472 L s P p .002 w .02381 .13244 m .03006 .13244 L s P [(3)] .01131 .13244 1 0 Mshowa p .002 w .02381 .25016 m .03006 .25016 L s P [(4)] .01131 .25016 1 0 Mshowa p .002 w .02381 .36788 m .03006 .36788 L s P [(5)] .01131 .36788 1 0 Mshowa p .002 w .02381 .4856 m .03006 .4856 L s P [(6)] .01131 .4856 1 0 Mshowa p .002 w .02381 .60332 m .03006 .60332 L s P [(7)] .01131 .60332 1 0 Mshowa p .001 w .02381 .03826 m .02756 .03826 L s P p .001 w .02381 .0618 m .02756 .0618 L s P p .001 w .02381 .08535 m .02756 .08535 L s P p .001 w .02381 .10889 m .02756 .10889 L s P p .001 w .02381 .15598 m .02756 .15598 L s P p .001 w .02381 .17952 m .02756 .17952 L s P p .001 w .02381 .20307 m .02756 .20307 L s P p .001 w .02381 .22661 m .02756 .22661 L s P p .001 w .02381 .2737 m .02756 .2737 L s P p .001 w .02381 .29724 m .02756 .29724 L s P p .001 w .02381 .32079 m .02756 .32079 L s P p .001 w .02381 .34433 m .02756 .34433 L s P p .001 w .02381 .39142 m .02756 .39142 L s P p .001 w .02381 .41497 m .02756 .41497 L s P p .001 w .02381 .43851 m .02756 .43851 L s P p .001 w .02381 .46205 m .02756 .46205 L s P p .001 w .02381 .50914 m .02756 .50914 L s P p .001 w .02381 .53269 m .02756 .53269 L s P p .001 w .02381 .55623 m .02756 .55623 L s P p .001 w .02381 .57977 m .02756 .57977 L s P p /Times-Italic findfont 12 scalefont setfont [(y)] .02381 .61803 0 -4 Mshowa P p .002 w .02381 0 m .02381 .61803 L s P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p .025 w .02381 .01472 Mdot .34127 .25016 Mdot .65873 .60332 Mdot .97619 .4856 Mdot P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup] The Unformatted text for this cell was not generated. Use options in the Actions Preferences dialog box to control when Unformatted text is generated. ;[o] -Graphics- :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] 2.3.2 Data entry for alpha decay. :[font = text; inactive; preserveAspect] Because of the broad range covered by the lifetime calculations it is best to plot the common log of the lifetime as a function of alpha energy. The Mathematica table below will hold the results from the theoretical calculations. The column on the left contains the alpha particle energy and the column on the right contains the lifetime calculation. Simply substitute your energies and lifetimes for those displayed here. Do NOT use the numbers displayed below. They are an example only. Energy Theoretical Lifetime ;[s] 3:0,0;149,1;160,2;590,-1; 3:1,11,8,Times,0,12,0,0,0;1,10,8,Times,2,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = input; Cclosed; preserveAspect; startGroup] theory = { { 8.78 , Log[10,5.9 10^-5] }, { 9.00 , Log[10,1.2 10^-12]}, { 12.0 , Log[10,2.0 10^-7 ] } } ;[s] 5:0,0;13,1;17,2;27,3;36,4;125,-1; 5:1,10,8,Courier,1,12,0,0,0;1,10,8,Courier,1,12,65535,0,0;1,10,8,Courier,1,12,0,0,0;1,10,8,Courier,1,12,65535,0,0;1,10,8,Courier,1,12,0,0,0; :[font = output; output; inactive; preserveAspect; endGroup] {{8.78, -4.229147988357856}, {9., -11.92081875395237}, {12., -6.698970004336018}} ;[o] {{8.78, -4.22915}, {9., -11.9208}, {12., -6.69897}} :[font = text; inactive; preserveAspect] The measured lifetimes taken from the table of isotopes are stored here in a manner similar to the table above. Substitute your results for the ones listed here. Energy Measured Lifetime :[font = input; Cclosed; preserveAspect; startGroup] data = { { 8.78 , Log[10,3.0 10^-7] }, { 9.00 , Log[10,2.0 10^-12]}, { 12.0 , Log[10,5.0 10^-8 ] } } ;[s] 5:0,0;13,1;17,2;27,3;36,4;125,-1; 5:1,10,8,Courier,1,12,0,0,0;1,10,8,Courier,1,12,65535,0,0;1,10,8,Courier,1,12,0,0,0;1,10,8,Courier,1,12,65535,0,0;1,10,8,Courier,1,12,0,0,0; :[font = output; output; inactive; preserveAspect; endGroup; endGroup] {{8.78, -6.522878745280337}, {9., -11.69897000433602}, {12., -7.30102999566398}} ;[o] {{8.78, -6.52288}, {9., -11.699}, {12., -7.30103}} :[font = subsubsection; inactive; Cclosed; preserveAspect; startGroup] 2.3.3 Plotting data and theory together. :[font = text; inactive; preserveAspect] The code in this cell will take the numbers you entered above and plot them together. It defines two graphical objects, "dataplot" and "theoryplot" and then combines them. Place the pointer anywhere in the cell, click, and hit "Enter". :[font = input; Cclosed; preserveAspect; startGroup] dataplot = ListPlot[ data, Frame->True, FrameLabel-> {FontForm["E(MeV)",{"Helvetica-Bold",10}], FontForm["Log(t)",{"Helvetica-Bold",10}], FontForm["Points->Data, Line->Theory",{"Helvetica-Bold",10}], " "}, Prolog->AbsolutePointSize[5], RotateLabel->False, DisplayFunction->Identity ] :[font = output; output; inactive; preserveAspect; endGroup] The Unformatted text for this cell was not generated. Use options in the Actions Preferences dialog box to control when Unformatted text is generated. ;[o] -Graphics- :[font = input; Cclosed; preserveAspect; startGroup] theoryplot = ListPlot[ Sort[theory], Frame->True, FrameLabel-> {FontForm["E(MeV)",{"Helvetica-Bold",10}], FontForm["Log(t)",{"Helvetica-Bold",10}], FontForm["Points->Data, Line->Theory",{"Helvetica-Bold",10}], " "}, PlotJoined->True, RotateLabel->False, PlotStyle->{RGBColor[1,0,0]}, DisplayFunction->Identity ] :[font = output; output; inactive; preserveAspect; endGroup] The Unformatted text for this cell was not generated. Use options in the Actions Preferences dialog box to control when Unformatted text is generated. ;[o] -Graphics- :[font = input; Cclosed; preserveAspect; startGroup] Show[{dataplot,theoryplot}, DisplayFunction->$DisplayFunction] :[font = postscript; PostScript; formatAsPostScript; output; inactive; preserveAspect; pictureLeft = 80; pictureWidth = 240; pictureHeight = 148] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Helvetica findfont 7 scalefont setfont % Scaling calculations -2.573055 0.29577 0.926954 0.076525 [ [(9)] .08888 0 0 2 0 Minner Mrotsboxa [(9.5)] .23676 0 0 2 0 Minner Mrotsboxa [(10)] .38465 0 0 2 0 Minner Mrotsboxa [(10.5)] .53253 0 0 2 0 Minner Mrotsboxa [(11)] .68042 0 0 2 0 Minner Mrotsboxa [(11.5)] .82831 0 0 2 0 Minner Mrotsboxa [(12)] .97619 0 0 2 0 Minner Mrotsboxa p /Helvetica-Bold findfont 10 scalefont setfont [(E\(MeV\))] .5 0 0 2 0 0 -1 Mouter Mrotsboxa P [(-12)] -0.0125 .00866 1 0 0 Minner Mrotsboxa [(-10)] -0.0125 .16171 1 0 0 Minner Mrotsboxa [(-8)] -0.0125 .31476 1 0 0 Minner Mrotsboxa [(-6)] -0.0125 .4678 1 0 0 Minner Mrotsboxa p /Helvetica-Bold findfont 10 scalefont setfont [(Log\(t\))] -0.0125 .30902 1 0 0 -1 0 Mouter Mrotsboxa P p /Helvetica-Bold findfont 10 scalefont setfont [(Points->Data, Line->Theory)] .5 .61803 0 -4 Msboxa P [( )] 1.0375 .30902 -1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .08888 0 m .08888 .00625 L s P [(9)] .08888 0 0 2 0 Minner Mrotshowa p .002 w .23676 0 m .23676 .00625 L s P [(9.5)] .23676 0 0 2 0 Minner Mrotshowa p .002 w .38465 0 m .38465 .00625 L s P [(10)] .38465 0 0 2 0 Minner Mrotshowa p .002 w .53253 0 m .53253 .00625 L s P [(10.5)] .53253 0 0 2 0 Minner Mrotshowa p .002 w .68042 0 m .68042 .00625 L s P [(11)] .68042 0 0 2 0 Minner Mrotshowa p .002 w .82831 0 m .82831 .00625 L s P [(11.5)] .82831 0 0 2 0 Minner Mrotshowa p .002 w .97619 0 m .97619 .00625 L s P [(12)] .97619 0 0 2 0 Minner Mrotshowa p .001 w .11846 0 m .11846 .00375 L s P p .001 w .14803 0 m .14803 .00375 L s P p .001 w .17761 0 m .17761 .00375 L s P p .001 w .20719 0 m .20719 .00375 L s P p .001 w .26634 0 m .26634 .00375 L s P p .001 w .29592 0 m .29592 .00375 L s P p .001 w .3255 0 m .3255 .00375 L s P p .001 w .35507 0 m .35507 .00375 L s P p .001 w .41423 0 m .41423 .00375 L s P p .001 w .4438 0 m .4438 .00375 L s P p .001 w .47338 0 m .47338 .00375 L s P p .001 w .50296 0 m .50296 .00375 L s P p .001 w .56211 0 m .56211 .00375 L s P p .001 w .59169 0 m .59169 .00375 L s P p .001 w .62127 0 m .62127 .00375 L s P p .001 w .65084 0 m .65084 .00375 L s P p .001 w .71 0 m .71 .00375 L s P p .001 w .73957 0 m .73957 .00375 L s P p .001 w .76915 0 m .76915 .00375 L s P p .001 w .79873 0 m .79873 .00375 L s P p .001 w .85788 0 m .85788 .00375 L s P p .001 w .88746 0 m .88746 .00375 L s P p .001 w .91704 0 m .91704 .00375 L s P p .001 w .94661 0 m .94661 .00375 L s P p .001 w .0593 0 m .0593 .00375 L s P p .001 w .02972 0 m .02972 .00375 L s P p .001 w .00015 0 m .00015 .00375 L s P p /Helvetica-Bold findfont 10 scalefont setfont [(E\(MeV\))] .5 0 0 2 0 0 -1 Mouter Mrotshowa P p .002 w 0 0 m 1 0 L s P p .002 w 0 .00866 m .00625 .00866 L s P [(-12)] -0.0125 .00866 1 0 0 Minner Mrotshowa p .002 w 0 .16171 m .00625 .16171 L s P [(-10)] -0.0125 .16171 1 0 0 Minner Mrotshowa p .002 w 0 .31476 m .00625 .31476 L s P [(-8)] -0.0125 .31476 1 0 0 Minner Mrotshowa p .002 w 0 .4678 m .00625 .4678 L s P [(-6)] -0.0125 .4678 1 0 0 Minner Mrotshowa p .001 w 0 .03927 m .00375 .03927 L s P p .001 w 0 .06988 m .00375 .06988 L s P p .001 w 0 .10049 m .00375 .10049 L s P p .001 w 0 .1311 m .00375 .1311 L s P p .001 w 0 .19232 m .00375 .19232 L s P p .001 w 0 .22293 m .00375 .22293 L s P p .001 w 0 .25354 m .00375 .25354 L s P p .001 w 0 .28415 m .00375 .28415 L s P p .001 w 0 .34537 m .00375 .34537 L s P p .001 w 0 .37597 m .00375 .37597 L s P p .001 w 0 .40658 m .00375 .40658 L s P p .001 w 0 .43719 m .00375 .43719 L s P p .001 w 0 .49841 m .00375 .49841 L s P p .001 w 0 .52902 m .00375 .52902 L s P p .001 w 0 .55963 m .00375 .55963 L s P p .001 w 0 .59024 m .00375 .59024 L s P p /Helvetica-Bold findfont 10 scalefont setfont [(Log\(t\))] -0.0125 .30902 1 0 0 -1 0 Mouter Mrotshowa P p .002 w 0 0 m 0 .61803 L s P P p p .002 w .08888 .61178 m .08888 .61803 L s P p .002 w .23676 .61178 m .23676 .61803 L s P p .002 w .38465 .61178 m .38465 .61803 L s P p .002 w .53253 .61178 m .53253 .61803 L s P p .002 w .68042 .61178 m .68042 .61803 L s P p .002 w .82831 .61178 m .82831 .61803 L s P p .002 w .97619 .61178 m .97619 .61803 L s P p .001 w .11846 .61428 m .11846 .61803 L s P p .001 w .14803 .61428 m .14803 .61803 L s P p .001 w .17761 .61428 m .17761 .61803 L s P p .001 w .20719 .61428 m .20719 .61803 L s P p .001 w .26634 .61428 m .26634 .61803 L s P p .001 w .29592 .61428 m .29592 .61803 L s P p .001 w .3255 .61428 m .3255 .61803 L s P p .001 w .35507 .61428 m .35507 .61803 L s P p .001 w .41423 .61428 m .41423 .61803 L s P p .001 w .4438 .61428 m .4438 .61803 L s P p .001 w .47338 .61428 m .47338 .61803 L s P p .001 w .50296 .61428 m .50296 .61803 L s P p .001 w .56211 .61428 m .56211 .61803 L s P p .001 w .59169 .61428 m .59169 .61803 L s P p .001 w .62127 .61428 m .62127 .61803 L s P p .001 w .65084 .61428 m .65084 .61803 L s P p .001 w .71 .61428 m .71 .61803 L s P p .001 w .73957 .61428 m .73957 .61803 L s P p .001 w .76915 .61428 m .76915 .61803 L s P p .001 w .79873 .61428 m .79873 .61803 L s P p .001 w .85788 .61428 m .85788 .61803 L s P p .001 w .88746 .61428 m .88746 .61803 L s P p .001 w .91704 .61428 m .91704 .61803 L s P p .001 w .94661 .61428 m .94661 .61803 L s P p .001 w .0593 .61428 m .0593 .61803 L s P p .001 w .02972 .61428 m .02972 .61803 L s P p .001 w .00015 .61428 m .00015 .61803 L s P p /Helvetica-Bold findfont 10 scalefont setfont [(Points->Data, Line->Theory)] .5 .61803 0 -4 Mshowa P p .002 w 0 .61803 m 1 .61803 L s P p .002 w .99375 .00866 m 1 .00866 L s P p .002 w .99375 .16171 m 1 .16171 L s P p .002 w .99375 .31476 m 1 .31476 L s P p .002 w .99375 .4678 m 1 .4678 L s P p .001 w .99625 .03927 m 1 .03927 L s P p .001 w .99625 .06988 m 1 .06988 L s P p .001 w .99625 .10049 m 1 .10049 L s P p .001 w .99625 .1311 m 1 .1311 L s P p .001 w .99625 .19232 m 1 .19232 L s P p .001 w .99625 .22293 m 1 .22293 L s P p .001 w .99625 .25354 m 1 .25354 L s P p .001 w .99625 .28415 m 1 .28415 L s P p .001 w .99625 .34537 m 1 .34537 L s P p .001 w .99625 .37597 m 1 .37597 L s P p .001 w .99625 .40658 m 1 .40658 L s P p .001 w .99625 .43719 m 1 .43719 L s P p .001 w .99625 .49841 m 1 .49841 L s P p .001 w .99625 .52902 m 1 .52902 L s P p .001 w .99625 .55963 m 1 .55963 L s P p .001 w .99625 .59024 m 1 .59024 L s P [( )] 1.0375 .30902 -1 0 Mshowa p .002 w 1 0 m 1 .61803 L s P P p P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath p p 5 Mabswid .02381 .42779 Mdot .08888 .03169 Mdot .97619 .36824 Mdot P p 1 0 0 r .004 w .02381 .60332 m .08888 .01472 L .97619 .41432 L s P P % End of Graphics MathPictureEnd :[font = output; output; inactive; preserveAspect; endGroup; endGroup; endGroup; endGroup] The Unformatted text for this cell was not generated. Use options in the Actions Preferences dialog box to control when Unformatted text is generated. ;[o] -Graphics- :[font = section; inactive; Cclosed; noPageBreak; preserveAspect; startGroup] 3. Procedure :[font = text; inactive; preserveAspect] 1. Clearly state the paradox associated with alpha decay. 2. Using the "Table of Isotopes" find seven examples of alpha emitters and make a table of nucleus, half-life, and alpha energy. See the accompanying sheet for examples of the data contained in the "Table of Isotopes". The nuclei that alpha decay are in the mass range A>200. It is best to find nuclei that emit alpha particles with primarily one kinetic energy. :[font = text; inactive; preserveAspect; endGroup] 3. The commands demonstrated in Section 2.2 can be used to calculate the transmission coefficient using the transfer matrix method and by breaking the potential barrier into a series of adjacent square barriers. The accuracy of the calculation depends on the number of barriers that the potential energy curve is divided into. It is important to "calibrate" the calculation by exploring the dependence of the result on the parameters of the calculation. For one of the entries in your table calculate the transmission coefficient using many values of the number of barriers. Does your result converge to a limiting value as the number of barriers increases? It is best to show your results with a plot. Section 2.3 demonstrates the procedure for making plots of data points with Mathematica. 4. Use the algorithm you programmed in part 3 to calculate the transmission coefficient for each of the entries in your table. Since the transmission coefficient represents the probability that a particle will penetrate the barrier during a collision, the inverse of it is the number of collisions that will occur before, on average, transmission will occur. You will use this fact in the next Section. ;[s] 3:0,0;780,1;791,2;1197,-1; 3:1,11,8,Times,0,12,0,0,0;1,10,8,Times,2,12,0,0,0;1,11,8,Times,0,12,0,0,0; :[font = section; inactive; Cclosed; preserveAspect; startGroup] 4. Analysis and Discussion :[font = text; inactive; preserveAspect; endGroup; endGroup] 1. If we know the number of collisions that must occur before a particle is transmitted (1/T), then the product of this number with the time between collisions will produce an estimate of the lifetime. Assume the kinetic energy of the alpha particle inside the nucleus is the same as when it is emitted. Show the time between collisions, Dt, is Dt = 2 R/v where R = 1.4A1/3 is the nuclear radius, v is the alpha particle speed inside the nucleus, and A is the mass number of the parent nucleus. Calculate the theoretical lifetimes for each of the entries in your table. Be sure to include one sample calculation in your report. 2. Make a plot comparing your theoretical and experimental results using the plotting commands defined in Section 2.4. How well does the theory describe the data? Does it reproduce the systematic trends of the data set? What attributes of the data does the theory fail to account for? Make at least one criticism of the procedure you followed to calculate the lifetime and why it might be wrong. Your report should include Purpose, Calculations, Analysis, and a Discussion. ;[s] 7:0,0;338,1;339,2;379,3;380,4;405,5;408,6;1139,-1; 7:1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,10,0,0,0;1,11,8,Times,0,12,0,0,0;1,0,0,Symbol,0,12,0,0,0;1,11,8,Times,0,12,0,0,0;1,8,6,Times,32,9,0,0,0;1,11,8,Times,0,12,0,0,0; ^*)