(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 81514, 1904]*) (*NotebookOutlinePosition[ 82348, 1932]*) (* CellTagsIndexPosition[ 82304, 1928]*) (*WindowFrame->Normal*) Notebook[{ Cell["Hartree-Self Consistant Field (SCF)", "Subtitle", Evaluatable->False, FormatType->TextForm], Cell["\<\ Application for He ground state\ \>", "Subsubtitle", Evaluatable->False, FormatType->TextForm], Cell[CellGroupData[{ Cell["Coulomb potential", "Subsubsection", Evaluatable->False, FormatType->TextForm], Cell["\<\ kounter=1; pots=Table[i,{i,1,10}];\ \>", "Input"], Cell[CellGroupData[{ Cell["\<\ vc[r_]=-Z/r /. Z->2; veff[r_]=vc[r]+l(l+1)/2/r^2 /. l->0\ \>", "Input"], Cell[BoxData[ \(\(-\(2\/r\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Print[\"kounter=\",kounter]", "Input"], Cell[BoxData[ InterpretationBox[\("kounter="\[InvisibleSpace]1\), SequenceForm[ "kounter=", 1], Editable->False]], "Print"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ x0=10^(-12); xend=4.0; pots[[kounter]]=Plot[veff[r],{r,x0,xend},PlotRange->{-10,0}]\ \>", "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.238095 0.618034 0.0618034 [ [.2619 .60553 -3 -9 ] [.2619 .60553 3 0 ] [.5 .60553 -3 -9 ] [.5 .60553 3 0 ] [.7381 .60553 -3 -9 ] [.7381 .60553 3 0 ] [.97619 .60553 -3 -9 ] [.97619 .60553 3 0 ] [.01131 0 -18 -4.5 ] [.01131 0 0 4.5 ] [.01131 .12361 -12 -4.5 ] [.01131 .12361 0 4.5 ] [.01131 .24721 -12 -4.5 ] [.01131 .24721 0 4.5 ] [.01131 .37082 -12 -4.5 ] [.01131 .37082 0 4.5 ] [.01131 .49443 -12 -4.5 ] [.01131 .49443 0 4.5 ] [ 0 0 0 0 ] [ 1 .62428 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2619 .61803 m .2619 .62428 L s [(1)] .2619 .60553 0 1 Mshowa .5 .61803 m .5 .62428 L s [(2)] .5 .60553 0 1 Mshowa .7381 .61803 m .7381 .62428 L s [(3)] .7381 .60553 0 1 Mshowa .97619 .61803 m .97619 .62428 L s [(4)] .97619 .60553 0 1 Mshowa .125 Mabswid .07143 .61803 m .07143 .62178 L s .11905 .61803 m .11905 .62178 L s .16667 .61803 m .16667 .62178 L s .21429 .61803 m .21429 .62178 L s .30952 .61803 m .30952 .62178 L s .35714 .61803 m .35714 .62178 L s .40476 .61803 m .40476 .62178 L s .45238 .61803 m .45238 .62178 L s .54762 .61803 m .54762 .62178 L s .59524 .61803 m .59524 .62178 L s .64286 .61803 m .64286 .62178 L s .69048 .61803 m .69048 .62178 L s .78571 .61803 m .78571 .62178 L s .83333 .61803 m .83333 .62178 L s .88095 .61803 m .88095 .62178 L s .92857 .61803 m .92857 .62178 L s .25 Mabswid 0 .61803 m 1 .61803 L s .02381 0 m .03006 0 L s [(-10)] .01131 0 1 0 Mshowa .02381 .12361 m .03006 .12361 L s [(-8)] .01131 .12361 1 0 Mshowa .02381 .24721 m .03006 .24721 L s [(-6)] .01131 .24721 1 0 Mshowa .02381 .37082 m .03006 .37082 L s [(-4)] .01131 .37082 1 0 Mshowa .02381 .49443 m .03006 .49443 L s [(-2)] .01131 .49443 1 0 Mshowa .125 Mabswid .02381 .0309 m .02756 .0309 L s .02381 .0618 m .02756 .0618 L s .02381 .09271 m .02756 .09271 L s .02381 .15451 m .02756 .15451 L s .02381 .18541 m .02756 .18541 L s .02381 .21631 m .02756 .21631 L s .02381 .27812 m .02756 .27812 L s .02381 .30902 m .02756 .30902 L s .02381 .33992 m .02756 .33992 L s .02381 .40172 m .02756 .40172 L s .02381 .43262 m .02756 .43262 L s .02381 .46353 m .02756 .46353 L s .02381 .52533 m .02756 .52533 L s .02381 .55623 m .02756 .55623 L s .02381 .58713 m .02756 .58713 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .07155 0 m .07299 .01958 L .08269 .11816 L .09312 .19344 L .10458 .25366 L .11478 .29454 L .12409 .32455 L .145 .3752 L .16409 .40824 L .18485 .43528 L .20462 .45527 L .22563 .47221 L .26338 .49519 L .30359 .51284 L .34227 .52562 L .38341 .53619 L .42304 .54432 L .46115 .55074 L .50171 .55645 L .54075 .5611 L .58225 .56533 L .62223 .56885 L .6607 .57182 L .70162 .57461 L .74102 .577 L .77891 .57906 L .81925 .58104 L .85807 .58276 L .89935 .58442 L .93911 .58588 L .97619 .58713 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOooool00ol0003oooooooooo`3ooooooe;ooooo000Pooooo`03o`000?oooooooooo03[o oooo00?o0000ooooooooool0oooooomAooooo`008?ooool00ol0003oooooooooo`0kooooo`03o`00 0?oooooooooo0?ooooooD?ooool0023ooooo00?o0000ooooooooool0>oooool00ol0003ooooooooo o`3ooooooe3ooooo000Pooooo`03o`000?oooooooooo03cooooo00?o0000ooooooooool0oooooom? ooooo`008?ooool2o`0003kooooo00?o0000ooooooooool0oooooom>ooooo`008?ooool00ol0003o ooooooooo`0nooooo`03o`000?oooooooooo0?ooooooCOooool0023ooooo00?o0000ooooooooool0 ?oooool00ol0003oooooooooo`3oooooodcooooo000Pooooo`03o`000?oooooooooo043ooooo00?o 0000ooooooooool0oooooom;ooooo`008?ooool00ol0003oooooooooo`11ooooo`03o`000?oooooo oooo0?ooooooB_ooool0023ooooo00?o0000ooooooooool0@_ooool00ol0003oooooooooo`3ooooo odWooooo000Pooooo`03o`000?oooooooooo04?ooooo00?o0000ooooooooool0oooooom8ooooo`00 8?ooool00ol0003oooooooooo`14ooooo`03o`000?oooooooooo0?ooooooAoooool0023ooooo00?o 0000ooooooooool0AOooool00ol0003oooooooooo`3oooooodKooooo000Pooooo`03o`000?oooooo oooo04Kooooo0_l0003oooooodKooooo000Pooooo`03o`000?oooooooooo04Sooooo00?o0000oooo ooooool0oooooom3ooooo`008?ooool00ol0003oooooooooo`19ooooo`03o`000?oooooooooo0?oo oooo@_ooool0023ooooo00?o0000ooooooooool0B_ooool2o`000?oooooo@_ooool0023ooooo0_l0 001=ooooo`03o`000?oooooooooo0?oooooo?oooool0023ooooo00?o0000ooooooooool0COooool0 0ol0003oooooooooo`3oooooockooooo000Pooooo`03o`000?oooooooooo04kooooo00?o0000oooo ooooool0oooooolmooooo`008?ooool00ol0003oooooooooo`1?ooooo`;o0000oooooolmooooo`00 8?ooool00ol0003oooooooooo`1Aooooo`03o`000?oooooooooo0?oooooo>_ooool0023ooooo00?o 0000ooooooooool0D_ooool2o`000?oooooo>_ooool0023ooooo00?o0000ooooooooool0E?ooool2 o`000?oooooo>?ooool0023ooooo00?o0000ooooooooool0E_ooool00ol0003oooooooooo`3ooooo ocGooooo000Pooooo`03o`000?oooooooooo05Oooooo0_l0003oooooocGooooo000Pooooo`03o`00 0?oooooooooo05Wooooo00?o0000ooooooooool0oooooolbooooo`004oooool6o`0000Oooooo00?o 0000ooooooooool0F_ooool2o`000?oooooo<_ooool001?ooooo00?oV9RHo`000?oSooooo000Pooooo`03o`000?oooooooooo0:Gooooo1ol0 003Qooooo`008?ooool00ol0003oooooooooo`2/ooooo`Oo0000f_ooool0023ooooo00?o0000oooo ooooool0/oooool6o`000=Cooooo000Pooooo`03o`000?oooooooooo0;Wooooo1Ol0003?ooooo`00 8?ooool2o`000;oooooo1ol00038ooooo`008?ooool00ol0003oooooooooo`35ooooo`So0000`?oo ool0023ooooo00?o0000ooooooooool0cOooool=o`000;?ooooo000Pooooo`03o`000?oooooooooo 0=[ooooo3?l0002Wooooo`008?ooool00ol0003oooooooooo`3Vooooo`So0000Woooool0023ooooo 00?o0000ooooooooool0k_oooolooooo`008?ooool00ol0003ooooooooo o`3oooooohkooooo000Pooooo`03o`000?oooooooooo0?ooooooS_ooool0023ooooo00?o0000oooo ooooool0oooooon>ooooo`005_ooooooo`0009[o00000Oooool003?ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo013ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`08ooooo`00POooool00ol0003oooooooooo`1Nooooo`03o`000?oo oooooooo05kooooo00?o0000ooooooooool0Goooool00ol0003oooooooooo`08ooooo`00oooooona ooooo`00oooooonaooooo`00\ \>"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.33004, -10.3471, \ 0.0154179, 0.0593967}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ If[kounter==1,kounter2=kounter,kounter2=kounter-1]; Show[pots[[kounter]],pots[[kounter2]]] \ \>", "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.238095 0.618034 0.0618034 [ [.2619 .60553 -3 -9 ] [.2619 .60553 3 0 ] [.5 .60553 -3 -9 ] [.5 .60553 3 0 ] [.7381 .60553 -3 -9 ] [.7381 .60553 3 0 ] [.97619 .60553 -3 -9 ] [.97619 .60553 3 0 ] [.01131 0 -18 -4.5 ] [.01131 0 0 4.5 ] [.01131 .12361 -12 -4.5 ] [.01131 .12361 0 4.5 ] [.01131 .24721 -12 -4.5 ] [.01131 .24721 0 4.5 ] [.01131 .37082 -12 -4.5 ] [.01131 .37082 0 4.5 ] [.01131 .49443 -12 -4.5 ] [.01131 .49443 0 4.5 ] [ 0 0 0 0 ] [ 1 .62428 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2619 .61803 m .2619 .62428 L s [(1)] .2619 .60553 0 1 Mshowa .5 .61803 m .5 .62428 L s [(2)] .5 .60553 0 1 Mshowa .7381 .61803 m .7381 .62428 L s [(3)] .7381 .60553 0 1 Mshowa .97619 .61803 m .97619 .62428 L s [(4)] .97619 .60553 0 1 Mshowa .125 Mabswid .07143 .61803 m .07143 .62178 L s .11905 .61803 m .11905 .62178 L s .16667 .61803 m .16667 .62178 L s .21429 .61803 m .21429 .62178 L s .30952 .61803 m .30952 .62178 L s .35714 .61803 m .35714 .62178 L s .40476 .61803 m .40476 .62178 L s .45238 .61803 m .45238 .62178 L s .54762 .61803 m .54762 .62178 L s .59524 .61803 m .59524 .62178 L s .64286 .61803 m .64286 .62178 L s .69048 .61803 m .69048 .62178 L s .78571 .61803 m .78571 .62178 L s .83333 .61803 m .83333 .62178 L s .88095 .61803 m .88095 .62178 L s .92857 .61803 m .92857 .62178 L s .25 Mabswid 0 .61803 m 1 .61803 L s .02381 0 m .03006 0 L s [(-10)] .01131 0 1 0 Mshowa .02381 .12361 m .03006 .12361 L s [(-8)] .01131 .12361 1 0 Mshowa .02381 .24721 m .03006 .24721 L s [(-6)] .01131 .24721 1 0 Mshowa .02381 .37082 m .03006 .37082 L s [(-4)] .01131 .37082 1 0 Mshowa .02381 .49443 m .03006 .49443 L s [(-2)] .01131 .49443 1 0 Mshowa .125 Mabswid .02381 .0309 m .02756 .0309 L s .02381 .0618 m .02756 .0618 L s .02381 .09271 m .02756 .09271 L s .02381 .15451 m .02756 .15451 L s .02381 .18541 m .02756 .18541 L s .02381 .21631 m .02756 .21631 L s .02381 .27812 m .02756 .27812 L s .02381 .30902 m .02756 .30902 L s .02381 .33992 m .02756 .33992 L s .02381 .40172 m .02756 .40172 L s .02381 .43262 m .02756 .43262 L s .02381 .46353 m .02756 .46353 L s .02381 .52533 m .02756 .52533 L s .02381 .55623 m .02756 .55623 L s .02381 .58713 m .02756 .58713 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .07155 0 m .07299 .01958 L .08269 .11816 L .09312 .19344 L .10458 .25366 L .11478 .29454 L .12409 .32455 L .145 .3752 L .16409 .40824 L .18485 .43528 L .20462 .45527 L .22563 .47221 L .26338 .49519 L .30359 .51284 L .34227 .52562 L .38341 .53619 L .42304 .54432 L .46115 .55074 L .50171 .55645 L .54075 .5611 L .58225 .56533 L .62223 .56885 L .6607 .57182 L .70162 .57461 L .74102 .577 L .77891 .57906 L .81925 .58104 L .85807 .58276 L .89935 .58442 L .93911 .58588 L .97619 .58713 L s .07155 0 m .07299 .01958 L .08269 .11816 L .09312 .19344 L .10458 .25366 L .11478 .29454 L .12409 .32455 L .145 .3752 L .16409 .40824 L .18485 .43528 L .20462 .45527 L .22563 .47221 L .26338 .49519 L .30359 .51284 L .34227 .52562 L .38341 .53619 L .42304 .54432 L .46115 .55074 L .50171 .55645 L .54075 .5611 L .58225 .56533 L .62223 .56885 L .6607 .57182 L .70162 .57461 L .74102 .577 L .77891 .57906 L .81925 .58104 L .85807 .58276 L .89935 .58442 L .93911 .58588 L .97619 .58713 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOooool00ol0003oooooooooo`3ooooooe;ooooo000Pooooo`03o`000?oooooooooo03[o oooo00?o0000ooooooooool0oooooomAooooo`008?ooool00ol0003oooooooooo`0kooooo`03o`00 0?oooooooooo0?ooooooD?ooool0023ooooo00?o0000ooooooooool0>oooool00ol0003ooooooooo o`3ooooooe3ooooo000Pooooo`03o`000?oooooooooo03cooooo00?o0000ooooooooool0oooooom? ooooo`008?ooool2o`0003kooooo00?o0000ooooooooool0oooooom>ooooo`008?ooool00ol0003o ooooooooo`0nooooo`03o`000?oooooooooo0?ooooooCOooool0023ooooo00?o0000ooooooooool0 ?oooool00ol0003oooooooooo`3oooooodcooooo000Pooooo`03o`000?oooooooooo043ooooo00?o 0000ooooooooool0oooooom;ooooo`008?ooool00ol0003oooooooooo`11ooooo`03o`000?oooooo oooo0?ooooooB_ooool0023ooooo00?o0000ooooooooool0@_ooool00ol0003oooooooooo`3ooooo odWooooo000Pooooo`03o`000?oooooooooo04?ooooo00?o0000ooooooooool0oooooom8ooooo`00 8?ooool00ol0003oooooooooo`14ooooo`03o`000?oooooooooo0?ooooooAoooool0023ooooo00?o 0000ooooooooool0AOooool00ol0003oooooooooo`3oooooodKooooo000Pooooo`03o`000?oooooo oooo04Kooooo0_l0003oooooodKooooo000Pooooo`03o`000?oooooooooo04Sooooo00?o0000oooo ooooool0oooooom3ooooo`008?ooool00ol0003oooooooooo`19ooooo`03o`000?oooooooooo0?oo oooo@_ooool0023ooooo00?o0000ooooooooool0B_ooool2o`000?oooooo@_ooool0023ooooo0_l0 001=ooooo`03o`000?oooooooooo0?oooooo?oooool0023ooooo00?o0000ooooooooool0COooool0 0ol0003oooooooooo`3oooooockooooo000Pooooo`03o`000?oooooooooo04kooooo00?o0000oooo ooooool0oooooolmooooo`008?ooool00ol0003oooooooooo`1?ooooo`;o0000oooooolmooooo`00 8?ooool00ol0003oooooooooo`1Aooooo`03o`000?oooooooooo0?oooooo>_ooool0023ooooo00?o 0000ooooooooool0D_ooool2o`000?oooooo>_ooool0023ooooo00?o0000ooooooooool0E?ooool2 o`000?oooooo>?ooool0023ooooo00?o0000ooooooooool0E_ooool00ol0003oooooooooo`3ooooo ocGooooo000Pooooo`03o`000?oooooooooo05Oooooo0_l0003oooooocGooooo000Pooooo`03o`00 0?oooooooooo05Wooooo00?o0000ooooooooool0oooooolbooooo`004oooool6o`0000Oooooo00?o 0000ooooooooool0F_ooool2o`000?oooooo<_ooool001?ooooo00?oV9RHo`000?oSooooo000Pooooo`03o`000?oooooooooo0:Gooooo1ol0 003Qooooo`008?ooool00ol0003oooooooooo`2/ooooo`Oo0000f_ooool0023ooooo00?o0000oooo ooooool0/oooool6o`000=Cooooo000Pooooo`03o`000?oooooooooo0;Wooooo1Ol0003?ooooo`00 8?ooool2o`000;oooooo1ol00038ooooo`008?ooool00ol0003oooooooooo`35ooooo`So0000`?oo ool0023ooooo00?o0000ooooooooool0cOooool=o`000;?ooooo000Pooooo`03o`000?oooooooooo 0=[ooooo3?l0002Wooooo`008?ooool00ol0003oooooooooo`3Vooooo`So0000Woooool0023ooooo 00?o0000ooooooooool0k_oooolooooo`008?ooool00ol0003ooooooooo o`3oooooohkooooo000Pooooo`03o`000?oooooooooo0?ooooooS_ooool0023ooooo00?o0000oooo ooooool0oooooon>ooooo`005_ooooooo`0009[o00000Oooool003?ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo013ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0@ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo00?o0000ooooooooool0 4Oooool00ol0003oooooooooo`08ooooo`00POooool00ol0003oooooooooo`1Nooooo`03o`000?oo oooooooo05kooooo00?o0000ooooooooool0Goooool00ol0003oooooooooo`08ooooo`00oooooona ooooo`00oooooonaooooo`00\ \>"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.33004, -10.3471, \ 0.0154179, 0.0593967}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["eigcoulomb=-Z^2/n^2/2 /. {Z->2,n->1}", "Input"], Cell[BoxData[ \(\(-2\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ eig=eigcoulomb; eig=-2.0 Clear[y]; Print[\"Eigenguess=\",eig] wrule=NDSolve[{-1/2 y''[x]+ veff[x]y[x]- eig \ y[x]==0,y[x0]==0.0,y'[x0]==0.1},y,{x,x0,xend},WorkingPrecision->6,MaxSteps->\ 6000,Compiled->True][[1,1]] \"keep on guessing for the correct eigenvalue\";\ \>", "Input"], Cell[BoxData[ \(\(-2.`\)\)], "Output"], Cell[BoxData[ InterpretationBox[\("Eigenguess="\[InvisibleSpace]\(-2.`\)\), SequenceForm[ "Eigenguess=", -2.0], Editable->False]], "Print"], Cell[BoxData[ RowBox[{\(NDSolve::"precw"\), \(\(:\)\(\ \)\), "\<\"The precision of the \ differential equation (\\!\\({\\(\\(\\(\\(\\(\\(2.`\\\\ \ \\(\\(y[x]\\)\\)\\)\\) + \\(\\(\[LeftSkeleton] 1 \[RightSkeleton]\\)\\) - \\(\ \\(\\(\\*SuperscriptBox[\\\"y\\\", \\\"\[Prime]\[Prime]\\\", \ Rule[MultilineFunction, None]]\\)\\)[x]\\)\\/2\\)\\) \[Equal] 0\\)\\), \\(\\(\ \\(\\(y[\\(\\(1\\/1000000000000\\)\\)]\\)\\) \[Equal] 0.`\\)\\), \\(\\(\\(\\(\ \\(\\(\\*SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", Rule[MultilineFunction, \ None]]\\)\\)[\\(\\(1\\/1000000000000\\)\\)]\\)\\) \[Equal] 0.1`\\)\\)}\\)) is \ less than WorkingPrecision (\\!\\(6.`\\)). \\!\\(\\*ButtonBox[\\\"More\ \[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"NDSolve::precw\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{"y", "\[Rule]", TagBox[\(InterpolatingFunction[{{1.`5.999999999999995*^-12, 4.`}}, "<>"]\), False, Editable->False]}]], "Output"] }, Open ]], Cell["fun[x_]:=y[x] /. wrule ", "Input"], Cell[CellGroupData[{ Cell["Plot[fun[x],{x,x0,xend},PlotRange->{-0.5,0.5}]", "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.238095 0.309017 0.618034 [ [.2619 .29652 -3 -9 ] [.2619 .29652 3 0 ] [.5 .29652 -3 -9 ] [.5 .29652 3 0 ] [.7381 .29652 -3 -9 ] [.7381 .29652 3 0 ] [.97619 .29652 -3 -9 ] [.97619 .29652 3 0 ] [.01131 .0618 -24 -4.5 ] [.01131 .0618 0 4.5 ] [.01131 .18541 -24 -4.5 ] [.01131 .18541 0 4.5 ] [.01131 .43262 -18 -4.5 ] [.01131 .43262 0 4.5 ] [.01131 .55623 -18 -4.5 ] [.01131 .55623 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2619 .30902 m .2619 .31527 L s [(1)] .2619 .29652 0 1 Mshowa .5 .30902 m .5 .31527 L s [(2)] .5 .29652 0 1 Mshowa .7381 .30902 m .7381 .31527 L s [(3)] .7381 .29652 0 1 Mshowa .97619 .30902 m .97619 .31527 L s [(4)] .97619 .29652 0 1 Mshowa .125 Mabswid .07143 .30902 m .07143 .31277 L s .11905 .30902 m .11905 .31277 L s .16667 .30902 m .16667 .31277 L s .21429 .30902 m .21429 .31277 L s .30952 .30902 m .30952 .31277 L s .35714 .30902 m .35714 .31277 L s .40476 .30902 m .40476 .31277 L s .45238 .30902 m .45238 .31277 L s .54762 .30902 m .54762 .31277 L s .59524 .30902 m .59524 .31277 L s .64286 .30902 m .64286 .31277 L s .69048 .30902 m .69048 .31277 L s .78571 .30902 m .78571 .31277 L s .83333 .30902 m .83333 .31277 L s .88095 .30902 m .88095 .31277 L s .92857 .30902 m .92857 .31277 L s .25 Mabswid 0 .30902 m 1 .30902 L s .02381 .0618 m .03006 .0618 L s [(-0.4)] .01131 .0618 1 0 Mshowa .02381 .18541 m .03006 .18541 L s [(-0.2)] .01131 .18541 1 0 Mshowa .02381 .43262 m .03006 .43262 L s [(0.2)] .01131 .43262 1 0 Mshowa .02381 .55623 m .03006 .55623 L s [(0.4)] .01131 .55623 1 0 Mshowa .125 Mabswid .02381 .09271 m .02756 .09271 L s .02381 .12361 m .02756 .12361 L s .02381 .15451 m .02756 .15451 L s .02381 .21631 m .02756 .21631 L s .02381 .24721 m .02756 .24721 L s .02381 .27812 m .02756 .27812 L s .02381 .33992 m .02756 .33992 L s .02381 .37082 m .02756 .37082 L s .02381 .40172 m .02756 .40172 L s .02381 .46353 m .02756 .46353 L s .02381 .49443 m .02756 .49443 L s .02381 .52533 m .02756 .52533 L s .02381 .0309 m .02756 .0309 L s .02381 .58713 m .02756 .58713 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .30902 m .04262 .34765 L .06244 .37914 L .07286 .39182 L .08255 .40141 L .09388 .41025 L .10458 .41692 L .11437 .42179 L .1199 .4239 L .125 .4254 L .12929 .42633 L .13401 .42708 L .13895 .42765 L .14357 .42806 L .14597 .42823 L .14854 .42837 L .15001 .42844 L .15134 .42849 L .15199 .42851 L .15269 .42853 L .15395 .42855 L .15514 .42856 L .15645 .42856 L .15711 .42856 L .15782 .42855 L .15909 .42852 L .16036 .42848 L .16157 .42842 L .16382 .42828 L .16605 .4281 L .16843 .42785 L .17269 .42727 L .18234 .42556 L .2201 .41619 L .2603 .4021 L .29899 .3893 L .34013 .37719 L .37975 .36745 L .41786 .36001 L .43911 .35675 L .45842 .35434 L .46806 .35334 L .4785 .35242 L .4883 .3517 L .49746 .35116 L .50693 .35074 L .51195 .35058 L .51734 .35046 L .52004 .35041 L .52146 .35039 L Mistroke .52298 .35038 L .52444 .35036 L .52576 .35036 L .5264 .35035 L .52709 .35035 L .52834 .35035 L .52955 .35035 L .53024 .35035 L .53088 .35035 L .53218 .35036 L .53291 .35036 L .53357 .35037 L .53604 .35039 L .53838 .35043 L .54299 .35052 L .54806 .35066 L .55336 .35085 L .55833 .35108 L .56721 .35159 L .57677 .35233 L .59684 .35454 L .61796 .35778 L .65612 .36594 L .67577 .37176 L .69673 .37961 L .73828 .39942 L .77832 .42592 L .81683 .46199 L .8578 .5113 L .89726 .57693 L Mfstroke .89726 .57693 m .91606 .61803 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgooooo`03odA4A?oooooo`03 o`000?oooooooooo05oooooo0ol0002Kooooo`;o0000MOooool002Sooooo00?o0000ooooooooool0 3oooool00ol0003oooooooooo`1Kooooo`?o0000X?ooool2o`0007?ooooo000Xooooo`03o`000?oo oooooooo013ooooo00?o0000ooooooooool0Eoooool3o`000:Gooooo0ol0001`ooooo`00:?ooool0 0ol0003oooooooooo`0Aooooo`;o0000E?ooool3o`000:_ooooo0_l0001^ooooo`00:?ooool00ol0 003oooooooooo`0Cooooo`03o`000?oooooooooo04kooooo0ol0002`ooooo`;o0000K?ooool002So oooo00?o0000ooooooooool05?ooool00ol0003oooooooooo`1:ooooo`?o0000]Oooool2o`0006[o oooo000Xooooo`03o`000?oooooooooo01Gooooo00?o0000ooooooooool0A_ooool3o`000;[ooooo 0_l0001Xooooo`00:?ooool3o`0001Kooooo00?o0000ooooooooool0@Oooool4o`000;oooooo00?o 0000ooooooooool0IOooool002Sooooo00?o0000ooooooooool05oooool00ol0003oooooooooo`0m ooooo`?o0000a?ooool2o`0006Gooooo000Xooooo`03o`000?oooooooooo01Sooooo00?o0000oooo ooooool0>Oooool3o`000;ooooo0_l0001Kooooo`00:?ooool00ol0003o ooooooooo`0Rooooo`Co00006Oooool4o`000>Oooooo00?o0000ooooooooool0F?ooool000_ooooo 00Co^[Zjoa0@4?l0003oMgMg1_ooool00omEEEGo0000oeEEE@04ooooo`Ko00001_ooool00ol0003o ooooooooo`0Vooooo`Co00004Oooool4o`000>cooooo00?o0000ooooooooool0Eoooool000[ooooo 00Kocooooo`00:?oo ool00ol0003oooooooooo`3oooooocKooooo00?o0000ooooooooool0COooool002Sooooo00?o0000 ooooooooool0oooooolgooooo`03o`000?oooooooooo04cooooo000Xooooo`03o`000?oooooooooo 0?oooooo>?ooool00ol0003oooooooooo`1;ooooo`00:?ooool00ol0003oooooooooo`3oooooocWo oooo00?o0000ooooooooool0B_ooool002Sooooo00?o0000ooooooooool0ooooooljooooo`03o`00 0?oooooooooo04Wooooo000Xooooo`?o0000oooooolkooooo`03o`000?oooooooooo04Sooooo000X ooooo`03o`000?oooooooooo0?oooooo>oooool00ol0003oooooooooo`18ooooo`00:?ooool00ol0 003oooooooooo`3ooooooccooooo00?o0000ooooooooool0Aoooool002Sooooo00?o0000oooooooo ool0oooooolmooooo`03o`000?oooooooooo04Kooooo000Xooooo`03o`000?oooooooooo0?oooooo ?_ooool00ol0003oooooooooo`15ooooo`00:?ooool00ol0003oooooooooo`3oooooocoooooo00?o 0000ooooooooool0A?ooool002Sooooo00?o0000ooooooooool0oooooom0ooooo`03o`000?oooooo oooo04?ooooo000Xooooo`03o`000?oooooooooo0?oooooo@Oooool00ol0003oooooooooo`12oooo o`00:?ooool00ol0003oooooooooo`3ooooood;ooooo00?o0000ooooooooool0@Oooool002Sooooo 00?o0000ooooooooool0oooooom3ooooo`03o`000?oooooooooo043ooooo000Xooooo`03o`000?oo oooooooo0?ooooooA?ooool00ol0003oooooooooo`0oooooo`00:?ooool00ol0003oooooooooo`3o ooooodCooooo00?o0000ooooooooool0?oooool002Sooooo00?o0000ooooooooool0oooooom5oooo o`03o`000?oooooooooo03kooooo000Xooooo`?o0000oooooom6ooooo`03o`000?oooooooooo03go oooo000Xooooo`03o`000?oooooooooo0?ooooooAoooool00ol0003oooooooooo`0looooo`00:?oo ool00ol0003oooooooooo`3oooooodSooooo00?o0000ooooooooool0>oooool002Sooooo00?o0000 ooooooooool0oooooom9ooooo`03o`000?oooooooooo03[ooooo000Xooooo`03o`000?oooooooooo 0?ooooooB_ooool00ol0003oooooooooo`0iooooo`00:?ooool00ol0003oooooooooo`3ooooood_o oooo00?o0000ooooooooool0>?ooool002Sooooo00?o0000ooooooooool0oooooom"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.430978, -0.5384, \ 0.0157692, 0.006075}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell["\<\ If we are satisfied with our numerical radial wave function we \ normalize it\ \>", "Text", Evaluatable->False, FormatType->TextForm], Cell["\<\ norm1=NIntegrate[fun[x]^2,{x,x0,xend}]; nfun[x_]=fun[x]/Sqrt[norm1];\ \>", "Input"], Cell["\<\ Lets construct the Hartree Self-Consistant Field for the 1s orbital\ \ \>", "Text", Evaluatable->False, FormatType->TextForm], Cell["We now iterate", "Text", FormatType->TextForm], Cell[BoxData[ \(\(kounter = kounter + 1;\)\)], "Input"], Cell[BoxData[ \(\(\(va[r_] := If[r > x0, 1/r\ NIntegrate[\ nfun[z]^2, {z, x0, r}], 0] + \n\t If[r < xend, NIntegrate[\ nfun[z]^2/z, {z, r, xend}], 0]\)\(\n\) \)\)], "Input"], Cell[BoxData[ \(veff[r_] := vc[r] + \ va[r] + \ l \(\((l + 1)\)/2\)/r^2\ /. \ l -> 0\)], "Input"], Cell["\"We now have a new effective potential\";", "Text", FormatType->TextForm], Cell[CellGroupData[{ Cell[BoxData[ \(Print["\", kounter]\)], "Input"], Cell[BoxData[ InterpretationBox[\("kounter="\[InvisibleSpace]2\), SequenceForm[ "kounter=", 2], Editable->False]], "Print"] }, Open ]], Cell["\<\ Lets manually loop back to the beginning and calculate a new \ wavefunction & eigenvalue and iterate again... and so on\ \>", "Text", Evaluatable->False, FormatType->TextForm], Cell["Part II", "Subsubtitle", Evaluatable->False, FormatType->TextForm], Cell["\<\ The eigenvalue obtained above can be interpreted as a laboratory \ measurable quantity: It's negative value is equal to the IONIZATION POTENTIAL of the atom. It is \ the energy required to remove an electron from the atom (for example in hydrogen \ that is just equal to -13.6 eV). This fact is a consequence of what is called \ Koopman's Theorem.\.08\.08 \ \>", "Text", Evaluatable->False, FormatType->TextForm], Cell["\<\ In the ground state both electrons see the same Self consistant potential that we have evaluated above\ \>", "Text", Evaluatable->False, FormatType->TextForm], Cell["\<\ If we now use this potential as a mean-background field we have the following Hamiltonian\ \>", "Text", Evaluatable->False], Cell[BoxData[ FormBox[ FormBox[ RowBox[{"H", "=", RowBox[{ \(p\_1\%2\/\(2 m\)\), "+", \(p\_2\%2\/\(2 m\)\), "+", \(U(r\_1)\), "+", \(U(r\_2)\), "+", SuperscriptBox["H", "\[Prime]", MultilineFunction->None]}]}], "TraditionalForm"], TraditionalForm]], "DisplayFormula"], Cell["where U is the self consistant potential gotten above", "Text", Evaluatable->False], Cell[BoxData[ FormBox[ FormBox[ RowBox[{ SuperscriptBox["H", "\[Prime]", MultilineFunction->None], "=", \(e\^2\/r\_12 - U(r\_1) - U(r\_2) - Ze\^2\/r\_1 - Ze\^2\/r\_2\)}], "TraditionalForm"], TraditionalForm]], "DisplayFormula", Evaluatable->False], Cell[TextData[ "\[CapitalDelta]E=<\[CurlyPhi]|H'|\[CurlyPhi]>"], "DisplayFormula", Evaluatable->False], Cell["\<\ H.W. Show that this energy shift is equal to the following integral\ \ \>", "Text", Evaluatable->False, FormatType->TextForm], Cell[TextData[{ "We find\n\n\[CapitalDelta]E=-\[Integral]\[DifferentialD]", Cell[BoxData[ \(TraditionalForm\`r\_1\)]], "\[Integral]\[DifferentialD]", Cell[BoxData[ \(TraditionalForm\`r\_1\)]], Cell[BoxData[ FormBox[ RowBox[{" ", RowBox[{"|", RowBox[{"\[CurlyPhi]", "(", FormBox[\(r\_1\), "TraditionalForm"], ")"}], \( | \^2\)}]}], TraditionalForm]]], " ", Cell[BoxData[ FormBox[ RowBox[{"|", RowBox[{"\[CurlyPhi]", "(", FormBox[\(r\_2\), "TraditionalForm"], ")"}], \( | \^2\)}], TraditionalForm]]], " ", Cell[BoxData[ \(TraditionalForm\`e\^2\/r\_12\)]], "= -I(1s,1s)" }], "DisplayFormula"], Cell[TextData[{ "Note: This integral is just the direct integral we considered in the \ previous ", StyleBox["Mathematica", FontSlant->"Italic"], " notebook,\nexcept that the Coulomb radial wavefunctions are replaced by \ the wavefunctions obtained above by\nthe SCF procedure, thus " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(I1s1s = NIntegrate[nfun[x]^2\ va[x], {x, x0, xend}]\)], "Input"], Cell[BoxData[ \(0.2916192590080639`\)], "Output"] }, Open ]], Cell[TextData[{ "The total energy of the He(", Cell[BoxData[ \(TraditionalForm\`1 s\^2\)]], ") ground state is equal to\n\nE = \[Epsilon](1s)+\[Epsilon](1s)-I(1s,1s)" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Energy = 2\ eig\ - \ I1s1s\)], "Input"], Cell[BoxData[ \(\(-4.291619259008064`\)\)], "Output"] }, Open ]] }, Open ]] }, FrontEndVersion->"5.2 for Macintosh", ScreenRectangle->{{0, 1920}, {0, 1103}}, WindowSize->{1304, 865}, WindowMargins->{{Automatic, 274}, {Automatic, 0}}, PrintingPageRange->{Automatic, Automatic}, PrintingOptions->{"PaperSize"->{612, 792}, "PaperOrientation"->"Portrait", "Magnification"->1}, Magnification->1.5, StyleDefinitions -> "Classroom.nb" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1754, 51, 101, 2, 57, "Subtitle", Evaluatable->False], Cell[1858, 55, 109, 5, 106, "Subsubtitle", Evaluatable->False], Cell[CellGroupData[{ Cell[1992, 64, 89, 2, 64, "Subsubsection", Evaluatable->False], Cell[2084, 68, 59, 3, 108, "Input"], Cell[CellGroupData[{ Cell[2168, 75, 81, 3, 108, "Input"], Cell[2252, 80, 47, 1, 108, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2336, 86, 44, 0, 88, "Input"], Cell[2383, 88, 139, 3, 34, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[2559, 96, 108, 4, 128, "Input"], Cell[2670, 102, 23913, 489, 279, 3235, 229, "GraphicsData", "PostScript", \ "Graphics"], Cell[26586, 593, 130, 3, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26753, 601, 116, 4, 128, "Input"], Cell[26872, 607, 24399, 521, 279, 3721, 261, "GraphicsData", "PostScript", \ "Graphics"], Cell[51274, 1130, 130, 3, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51441, 1138, 53, 0, 88, "Input"], Cell[51497, 1140, 40, 1, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51574, 1146, 288, 9, 208, "Input"], Cell[51865, 1157, 42, 1, 90, "Output"], Cell[51910, 1160, 155, 3, 34, "Print"], Cell[52068, 1165, 817, 11, 124, "Message"], Cell[52888, 1178, 189, 5, 92, "Output"] }, Open ]], Cell[53092, 1186, 40, 0, 88, "Input"], Cell[CellGroupData[{ Cell[53157, 1190, 63, 0, 88, "Input"], Cell[53223, 1192, 23358, 517, 279, 3942, 273, "GraphicsData", "PostScript", \ "Graphics"], Cell[76584, 1711, 130, 3, 90, "Output"] }, Open ]], Cell[76729, 1717, 147, 5, 41, "Text", Evaluatable->False], Cell[76879, 1724, 93, 3, 108, "Input"], Cell[76975, 1729, 139, 5, 41, "Text", Evaluatable->False], Cell[77117, 1736, 56, 1, 41, "Text"], Cell[77176, 1739, 59, 1, 89, "Input"], Cell[77238, 1742, 196, 4, 135, "Input"], Cell[77437, 1748, 109, 2, 89, "Input"], Cell[77549, 1752, 82, 1, 41, "Text"], Cell[CellGroupData[{ Cell[77656, 1757, 63, 1, 89, "Input"], Cell[77722, 1760, 139, 3, 34, "Print"] }, Open ]], Cell[77876, 1766, 190, 6, 79, "Text", Evaluatable->False], Cell[78069, 1774, 76, 2, 83, "Subsubtitle", Evaluatable->False], Cell[78148, 1778, 430, 12, 193, "Text", Evaluatable->False], Cell[78581, 1792, 173, 6, 117, "Text", Evaluatable->False], Cell[78757, 1800, 135, 4, 79, "Text", Evaluatable->False], Cell[78895, 1806, 356, 9, 65, "DisplayFormula"], Cell[79254, 1817, 91, 1, 41, "Text", Evaluatable->False], Cell[79348, 1820, 308, 8, 69, "DisplayFormula", Evaluatable->False], Cell[79659, 1830, 104, 2, 39, "DisplayFormula", Evaluatable->False], Cell[79766, 1834, 139, 5, 40, "Text", Evaluatable->False], Cell[79908, 1841, 751, 26, 100, "DisplayFormula"], Cell[80662, 1869, 307, 7, 117, "Text"], Cell[CellGroupData[{ Cell[80994, 1880, 84, 1, 89, "Input"], Cell[81081, 1883, 53, 1, 90, "Output"] }, Open ]], Cell[81149, 1887, 192, 5, 117, "Text"], Cell[CellGroupData[{ Cell[81366, 1896, 60, 1, 89, "Input"], Cell[81429, 1899, 57, 1, 90, "Output"] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)