{VERSION 4 0 "IBM INTEL NT" "4.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 256 "" 1 16 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 257 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 260 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 } {CSTYLE "" -1 261 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 265 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 266 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 268 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 271 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 273 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 1 3 1 } 1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE " " -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Map le Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT 256 12 "Atomic Model" }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 290 "A simple model fo r the structure of few-electron atoms is presented. A model potential \+ for the helium atom is introduced (it represents a good approximation \+ to the so-called Hartee-Fock potential for the helium atom). The eigen energies and eigenfunctions of this potential are determined by" }} {PARA 0 "" 0 "" {TEXT -1 96 "a) the variational method for the lowest- lying eigenstate for a given angular momentum symmetry;" }}{PARA 0 "" 0 "" {TEXT -1 62 "b) the numerical solution of the radial Schroedinger equation;" }}{PARA 0 "" 0 "" {TEXT -1 62 "c) a matrix diagonalization using a Slater-type orbital basis." }}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{PARA 0 "" 0 "" {TEXT -1 255 "A simple model potential for the ground state of helium can be obtained as follows: assume that in the ground state the wavefunction corresponds to an 1s^2 configuration (one spin -up, one spin-down electron). What potential should the electron exper ience?" }}{PARA 0 "" 0 "" {TEXT -1 96 "The nuclear attraction -2/r and a repulsive potential due to the presence of the other electron." }} {PARA 0 "" 0 "" {TEXT -1 133 "What is the repulsive potential due to t he other electron (which has the same 1s-wavefunction, it has just opp osite spin projection)?" }}{PARA 0 "" 0 "" {TEXT -1 426 "We have an ch icken-and-egg problem here: given a potential, we can determine the |1 s> state, and calculate the repulsive part of the potential correctly. This is the objective of a self-consistent field calculation. Then we would need to repeat the calculation of the new |1s> state using the \+ new potential. Repeating this procedure one can come up with a |1s> st ate and corresponding potential that provide the lowest energy." }} {PARA 0 "" 0 "" {TEXT -1 89 "We break the closed circuit by stating an approximate answer for the potential in helium:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "V:=r->-2/r+1/r*(1-exp(-3.36*r)*(1+1.665*r)); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"VGR6#%\"rG6\"6$%)operatorG%&ar rowGF(,&*&\"\"\"F.9$!\"\"!\"#*&,&F.F.*&-%$expG6#,$F/$!$O$F1F.,&F.F.*&$ \"%l;!\"$F.F/F.F.F.F0F.F/F0F.F(F(F(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 200 "This is the potential experienced by an 1s-electron in the hel ium atom: at short distances the electron just feels the full nuclear \+ attraction (-2/r) as the second expression approaches a constant as " }{TEXT 272 1 "r" }{TEXT -1 76 " goes to zero. At large distances the e lectronic repulsion term goes like 1/" }{TEXT 273 1 "r" }{TEXT -1 80 " and screens the nucleus by one unit so that the overall potential goe s like -1/" }{TEXT 274 1 "r" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 69 "We can check for the consistenc y by a simple variational calculation:" }}{PARA 0 "" 0 "" {TEXT -1 112 "We start with an unnormalized 1s state that depends on a 'charge' parameter (cf. the hydrogen-like wavefunction)" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 20 "chi:=r*exp(-beta*r);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$chiG*&%\"rG\"\"\"-%$expG6#,$*&%%betaGF'F&F'!\"\"F'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 134 "We are using the radial wavefu nction (i.e. r*R_nl(r)) so that the radial kinetic energy is just prop ortional to the second derivative." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "assume(beta>0);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "E1sT:=Int(chi*(-1/2*Diff(chi,r$2)+V(r)*chi),r=0..infi nity)/Int(chi^2,r=0..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%% E1sTG*&-%$IntG6$*(%\"rG\"\"\"-%$expG6#,$*&%&beta|irGF+F*F+!\"\"F+,&-%% DiffG6$*&F*F+F,F+-%\"$G6$F*\"\"##F2F;*(,&*&F+F+F*F2!\"#*&,&F+F+*&-F-6# ,$F*$!$O$F@F+,&F+F+*&$\"%l;!\"$F+F*F+F+F+F2F+F*F2F+F+F*F+F,F+F+F+/F*; \"\"!%)infinityGF+-F'6$*&)F*F;F+)F,F;F+FNF2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "E1sT:=value(E1sT);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%E1sTG,$*&*&,,$!'/Ff\"\"!\"\"\"*&$\"&+D'F+F,)%&beta|irG\"\"%F, F,*&$\"&+]'F+F,)F1\"\"$F,F,*&$\"'D*=&F+F,)F1\"\"#F,!\"\"*&$\"'[?wF+F,F 1F,F=F,F1F,F,,**$F6F,$\"&Dc\"F+*&$\"&](yF+F,F;F,F,*&$\"'+B8F+F,F1F,F,$ \"&)3uF+F,F=$\"++++]7!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "plot(E1sT,beta=1.55..1.8);" }}{PARA 13 "" 1 "" {GLPLOT2D 693 299 299 {PLOTDATA 2 "6%-%'CURVESG6$7S7$$\"3/++++++]:!#<$!3[dk:#[*y'*))!#=7$$\" 3km;a)G\\ab\"F*$!3%yPhZ#HS-*)F-7$$\"3G$ek`o!>g:F*$!3RVTChe52*)F-7$$\"3 mm\"z>)G_l:F*$!3CcW(p&H>7*)F-7$$\"3ymT&QU!*3d\"F*$!3*yQz\"Gw4<*)F-7$$ \"3S$eRZXKid\"F*$!3],9fLLw@*)F-7$$\"3i;z>,_=\"e\"F*$!3\")*)eW]l*e#*)F- 7$$\"3+]7G$[8je\"F*$!34\"R&4#z!)*H*)F-7$$\"3u;z%*frh\"f\"F*$!3W%>iNA&* R$*)F-7$$\"3%*\\ilFQ!pf\"F*$!3Y\"[shY%yP*)F-7$$\"3IL$3_\"=M-;F*$!3xx\" 3)p1YT*)F-7$$\"3omTg(fJrg\"F*$!36'>6GM7X%*)F-7$$\"3/+]7eP_7;F*$!3s68Kv 'Qx%*)F-7$$\"38+]Pf!Qzh\"F*$!3S4!)pk]v]*)F-7$$\"3++](=ubJi\"F*$!3[?(>! ))*\\M&*)F-7$$\"3o;zW(*Q*yi\"F*$!3kW n*)F-7$$\"3S$e*[.-dp;F*$!3g7^k^-Ao*)F-7$$\"3#o;/Egw[n\"F*$!3#\\\\I4h\\ )o*)F-7$$\"3i;z%*f%)Q!o\"F*$!3gJ<:X?Fp*)F-7$$\"3+](oza'=&o\"F*$!3I*)f \"yjZ%p*)F-7$$\"3um\"zWho.p\"F*$!3!fH7@9O%p*)F-7$$\"3'****\\i>Adp\"F*$ !3uT]#yh/#p*)F-7$$\"3'**\\i:jf4q\"F*$!3gf:M==wo*)F-7$$\"3(*\\7`>r-1n*)F-7$$\"3OLLeM%4nr\" F*$!3GcNs-&Qh'*)F-7$$\"3,+]P4v5A#o\")\\*)F-7$$\"3!o;aQz]Ow\"F*$!3)3CD.*evY*)F-7$$\"3K$ekG=4*o+N%*)F-7$$\"3)****\\i4TPx\"F*$!3:?u(z1:.%*)F-7$$\"3YL3F9!z#z'p\"ROk$*)F-7$$\"3ymm;%>KUy\"F*$!3?y&[)f0wK*)F-7$$\"31]7.EigW#*)F-7$$\"3/+++++++= F*$!3!)Q(*fdov>*)F--%'COLOURG6&%$RGBG$\"#5!\"\"$\"\"!F`[lF_[l-%+AXESLA BELSG6$Q&beta|ir6\"Q!6\"-%%VIEWG6$;$\"$b\"!\"#$\"#=F^[l%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 44.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "beta0:=fsolve(diff(E1sT,beta ),beta=1..2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&beta0G$\"+PG\\(o\" !\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "E0:=subs(beta=beta0 ,E1sT);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#E0G$!+K$o%p*)!#5" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 606 "Note that the beta-value is less \+ than 2. The value of beta=2 would be obtained if there was only one 1s -electron (the hydrogen-like solution for the ground state of the He+ \+ ion). The fact that the electron wants to have a slightly more diffuse wavefunction reflects the so-called inner screening: the electrostati c repulsion when combined with the simple model of two identical elect rons (apart from the spin projection) leads to two electrons which rep el each other on average (independent electron model). Both electrons \+ experience the same common central potential and are bound by the same eigenenergy." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 38 "Now we would lik e to check two things:" }}{PARA 0 "" 0 "" {TEXT -1 45 "1) how accurate is this variational solution?" }}{PARA 0 "" 0 "" {TEXT -1 85 "2) how \+ close is the used potential to the potential produced by this 1s-wavef unction?" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 98 "Part 1 can be answered by using dsolve[numeric]. For part 2 we nee d to solve the Poisson equation." }}{PARA 0 "" 0 "" {TEXT -1 119 "Let \+ us start with the first question: To solve the SE we start the integra tion not at zero, but at some small value of " }{TEXT 257 5 "r=eta" } {TEXT -1 26 " to avoid the singularity." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "eta:=10^(-8);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$e taG#\"\"\"\"*++++\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "IC:= phi(eta)=eta,D(phi)(eta)=1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#ICG6 $/-%$phiG6##\"\"\"\"*++++\"F*/--%\"DG6#F(F)F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "Et:=-0.9042;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%#EtG$!%U!*!\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "SE:=-1/ 2*diff(phi(r),r$2)+(V(r)-Et)*phi(r);" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%#SEG,&-%%diffG6$-%$phiG6#%\"rG-%\"$G6$F,\"\"##!\"\"F0*&,(*&\"\"\"F 6F,F2!\"#*&,&F6F6*&-%$expG6#,$F,$!$O$F7F6,&F6F6*&$\"%l;!\"$F6F,F6F6F6F 2F6F,F2F6$\"%U!*!\"%F6F6F)F6F6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "sol:=dsolve(\{SE,IC\},phi(r),numeric,output=listprocedure):" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "phir:=subs(sol,phi(r)):" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "plot(phir,0..10);" }} {PARA 13 "" 1 "" {GLPLOT2D 693 299 299 {PLOTDATA 2 "6%-%'CURVESG6$7go7 $$\"+v1h6o!#7$\"+pn]>nF*7$$\"+N@Ki8!#6$\"+[WvD8F07$$\"+-K[V?F0$\"+=h\" ='>F07$$\"+qUkCFF0$\"+$QK0e#F07$$\"+0k'p3%F0$\"+_dZnPF07$$\"+S&)G\\aF0 $\"+91g*)[F07$$F)F0$\"+j2!)\\fF07$$\"+5G$R<)F0$\"+Q:&3&pF07$$\"+X\\DO& *F0$\"+KOT&*yF07$$\"+3x&)*3\"!#5$\"+Fn.'y)F07$$F/FX$\"+5L^T5FX7$$\"+il yM;FX$\"+/kk&=\"FX7$$\"+*)4D2>FX$\"+&y\"o78FX7$$\"+;arz@FX$\"+:i:C9FX7 $$\"+)4bQl#FX$\"+0P:&e\"FX7$$\"+!y%*z7$FX$\"+#o&y4#>FX7$$\"+oMrU^F X$\"+?^^\\>FX7$$\"+*HB$4aFX$\"+E^1c>FX7$$\"+IJ$fn&FX$\"+2')Re>FX7$$\"+ hHaUfFX$\"+@Y#p&>FX7$$\"+\"z_\"4iFX$\"+Yt,_>FX7$$\"+m6m#G(FX$\"+f5//>F X7$$\"+S&phN)FX$\"+&o*\\B=FX7$$\"+:ddC%*FX$\"+?DgB/x\"F\\t$\"+\"o#zr) )F07$$\"+1J:w=F\\t$\"+$)f7?!)F07$$\"+3En$4#F\\t$\"+eA%HZ'F07$$\"+/RE&G #F\\t$\"+!y@eK&F07$$\"+D.&4]#F\\t$\"+L&y4D%F07$$\"+vB_$)=F*7$$\"+=>Y 2aF\\t$\"+'QJZ\\\"F*7$$\"+yXu9cF\\t$\"+mr+k6F*7$$\"+\\y))GeF\\t$\"+uB$ y)*)!#87$$\"+i_QQgF\\t$\"+L)z:)pFf\\l7$$\"+!y%3TiF\\t$\"+/o^uaFf\\l7$$ \"+O![hY'F\\t$\"+w_)>>%Ff\\l7$$\"+#Qx$omF\\t$\"++.%\\J$Ff\\l7$$\"+u.I% )oF\\t$\"+77G/EFf\\l7$$\"+(pe*zqF\\t$\"+I " 0 "" {MPLTEXT 1 0 8 "E 1s:=Et;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$E1sG$!%U!*!\"%" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 272 "The numerically exact solution fo r the radial SE yields a slightly lower eigenvalue of E_1s = -0.9042 a .u. (1 a.u. = 27.12 eV) compared to the variational result of E_1s_v = -0.897 a.u.. This means that the true eigenfunction is somewhat differ ent from the hydrogenic form." }}{PARA 0 "" 0 "" {TEXT -1 243 "The num erical eigenfunction is not normalized properly. To compare the graphs we simply change the normalization of the variational answer. In fact , the variational state was normalized such that the derivative of the function equals unity at " }{TEXT 258 1 "r" }{TEXT -1 3 "=0." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "eval(subs(r=0,diff(subs(beta =beta0,chi),r)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"\"\"\"\"!" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoords has been redefined\n " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "P1:=plot(phir,0..6,colo r=red):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "P2:=plot(subs(be ta=beta0,chi),r=0..6,color=blue):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "display([P1,P2]);" }}{PARA 13 "" 1 "" {GLPLOT2D 695 299 299 {PLOTDATA 2 "6&-%'CURVESG6$7do7$$\"+0k'p3%!#7$\"+r(4P0%F*7$$\" +5G$R<)F*$\"+Y6]T!)F*7$$\"+A**3E7!#6$\"+.fU'>\"F57$$\"+ilyM;F5$\"+vdG# e\"F57$$\"+V)z@X#F5$\"+vA5NBF57$$\"+DJdpKF5$\"+.#yK1$F57$$\"+(ofV!\\F5 $\"+#*4M[WF57$$\"+]i9RlF5$\"+E)zDu&F57$$\"+8G$R<)F5$\"+k:&3&pF57$$\"+v $>(3)*F5$\"+xOxx!)F57$$\"+%f]V9\"!#5$\"+)oLx7*F57$$\"+]#HyI\"Fhn$\"+c) *[55Fhn7$$\"+gIJ#f\"Fhn$\"+jPMk6Fhn7$$\"+pozw=Fhn$\"+8AG*H\"Fhn7$$\"+y 1Gh@Fhn$\"+0b2<9Fhn7$$\"+)[kdW#Fhn$\"+#pO$>:Fhn7$$\"+#3Gc3$Fhn$\"+?t(* *p\"Fhn7$$\"+v;\\DPFhn$\"+rgkB=Fhn7$$\"++nfpVFhn$\"+\"GjB!>Fhn7$$\"+D< q8]Fhn$\"+_.qW>Fhn7$$\"+^E'R<&Fhn$\"+L^^]>Fhn7$$\"+yNAM`Fhn$\"+a_ma>Fh n7$$\"+0X[%\\&Fhn$\"+2&[s&>Fhn7$$\"+JauacFhn$\"+rvNe>Fhn7$$\"+ej+:eFhn $\"+#H!3e>Fhn7$$\"+&Gn_(fFhn$\"+J**\\c>Fhn7$$\"+6#Gb8'Fhn$\"+7ap`>Fhn7 $$\"+Q\"*y&H'Fhn$\"+Z:u\\>Fhn7$$\"+8(=,*oFhn$\"+^8=E>Fhn7$$\"+)G[W[(Fh n$\"+1,%4*=Fhn7$$\"+))fB:()Fhn$\"+k4i\"z\"Fhn7$$\"+)Q=\"))**Fhn$\"+og2 m;Fhn7$$\"+k=pD6!\"*$\"+LRNI:Fhn7$$\"+lN?c7Fgu$\"+/&z$)Q\"Fhn7$$\"+U$e 6P\"Fgu$\"+A\"peE\"Fhn7$$\"+&>q0]\"Fgu$\"+-%))R8\"Fhn7$$\"+DM^I;Fgu$\" +Sk**45Fhn7$$\"+0ytb@Fgu$\"+yZs3jF57$$\"+4wY_AFgu$\"+#*ep3bF5 7$$\"+IOTqBFgu$\"+\"\\cc([F57$$\"+4\">)*\\#Fgu$\"+cw1cUF57$$\"+EP/BEFg u$\"+/j\"Ht$F57$$\"+)o:;v#Fgu$\"+'pb,D$F57$$\"+%)[opGFgu$\"+8)o\"eGF57 $$\"+i%Qq*HFgu$\"+kq)[[#F57$$\"+RIKHJFgu$\"+u\\!e9#F57$$\"+^rZWKFgu$\" +g\")f')=F57$$\"+[n%)oLFgu$\"+%*32S;F57$$\"+5FL(\\$Fgu$\"+Set<9F57$$\" +e6.BOFgu$\"+K#)GG7F57$$\"+p3lWPFgu$\"+.bDo5F57$$\"+A))ozQFgu$\"+rrxS \"*F*7$$\"+Ik-,SFgu$\"+E*y.%zF*7$$\"+D-eITFgu$\"+(p%HFoF*7$$\"+>_(zC%F gu$\"+OtZ]fF*7$$\"+b*=jP%Fgu$\"+hM2<^F*7$$\"+4/3(\\%Fgu$\"+(RgtV%F*7$$ \"+C4JBYFgu$\"+%fU6#QF*7$$\"+DVsYZFgu$\"+x>()*H$F*7$$\"+>n#f([Fgu$\"+Q E!)GGF*7$$\"+!)RO+]Fgu$\"+b9uPCF*7$$\"+_!>w7&Fgu$\"+]3%G4#F*7$$\"+*Q?Q D&Fgu$\"++.P)z\"F*7$$\"+5jyp`Fgu$\"+VU.k:F*7$$\"+Ujp-bFgu$\"+)RBCL\"F* 7$$\"+gEd@cFgu$\"+f%)Ga6F*7$$\"+4'>$[dFgu$\"+%ePU!**!#87$$\"+6EjpeFgu$ \"+*GDVb)Fcal7$\"\"'$\"+2?F4tFcal-%'COLOURG6&%$RGBG$\"*++++\"!\")$\"\" !FeblFdbl-F$6$7_o7$FdblFdbl7$$\"3$*****\\ilyM;!#>$\"31P>#[C/.f\"F]cl7$ $\"3')*****\\7t&pKF]cl$\"3Y_$)ypX1%4$F]cl7$$\"3z****\\(ofV!\\F]cl$\"3D dg9ZK\"[^%F]cl7$$\"3s******\\i9RlF]cl$\"3Wb?&*)Raf&eF]cl7$$\"3m****\\7 G$R<)F]cl$\"3Bd)fdNo27(F]cl7$$\"3e*****\\P>(3)*F]cl$\"30#p')H(eT7$)F]c l7$$\"3&****\\Pf]V9\"!#=$\"38,&z%4A%RV*F]cl7$$\"3%*******\\#HyI\"F\\el $\"3#Grls5G)[5F\\el7$$\"37+]PfIJ#f\"F\\el$\"35i8\"=J=r@\"F\\el7$$\"3-+ +voozw=F\\el$\"3/&z_\"3zKn8F\\el7$$\"3#***\\7y1Gh@F\\el$\"3TXE,@@y+:F \\el7$$\"33++]([kdW#F\\el$\"3cDksl;s=;F\\el7$$\"3K++D\"3Gc3$F\\el$\"3x 'R(H\\,>L=F\\el7$$\"3++++v;\\DPF\\el$\"3C/5+o*)z')>F\\el7$$\"3A++++nfp VF\\el$\"3UWKmDLH!4#F\\el7$$\"3W+++Dq0]\"Fa[ m$\"3Gl$RP#4r#>\"F\\el7$$\"3'******\\U80j\"Fa[m$\"3IR/r/6\"3/\"F\\el7$ $\"35+++0ytb@Fa[m$\"3>) 4J0fX,$fF]cl7$$\"3'****\\(3wY_AFa[m$\"3*pNw)>>zL]F]cl7$$\"3#)******HOT qBFa[m$\"3jJy#Hz?8M%F]cl7$$\"37++v3\">)*\\#Fa[m$\"379+1--=!o$F]cl7$$\" 3:++DEP/BEFa[m$\"3c^vlcdfOJF]cl7$$\"3=++](o:;v#Fa[m$\"3Y&3Sc/)e[EF]cl7 $$\"3=++v$)[opGFa[m$\"39A`bL_CjAF]cl7$$\"3%*****\\i%Qq*HFa[m$\"3C8Kub` e1>F]cl7$$\"3&****\\(QIKHJFa[m$\"3)znWsraCf\"F]cl7$$\"3#****\\7:xWC$Fa [m$\"3%yo\"yRmYf8F]cl7$$\"37++]Zn%)oLFa[m$\"3_L-+()*[V9\"F]cl7$$\"3y** ****4FL(\\$Fa[m$\"3kO$=j(zCk&*!#?7$$\"3#)****\\d6.BOFa[m$\"3OTx2UxH9!) Fcam7$$\"3(****\\(o3lWPFa[m$\"3'o=#=P5UYnFcam7$$\"3!*****\\A))ozQFa[m$ \"39P,h<\"o`c&Fcam7$$\"3e******Hk-,SFa[m$\"3K$e0i&puwYFcam7$$\"36+++D- eITFa[m$\"3=tGglA0!)QFcam7$$\"3u***\\(=_(zC%Fa[m$\"3I3FATh?tKFcam7$$\" 3M+++b*=jP%Fa[m$\"3%=n)yXjW:FFcam7$$\"3g***\\(3/3(\\%Fa[m$\"3;#>fT%*Qf F#Fcam7$$\"33++vB4JBYFa[m$\"3O)R$p#y<4*=Fcam7$$\"3u*****\\KCnu%Fa[m$\" 3g8Am'p/kd\"Fcam7$$\"3s***\\(=n#f([Fa[m$\"3_^BP1W4-8Fcam7$$\"3P+++!)RO +]Fa[m$\"3K!=EwM2C3\"Fcam7$$\"30++]_!>w7&Fa[m$\"3q.`p$QYX&*)!#@7$$\"3O ++v)Q?QD&Fa[m$\"3pf8#3:\\]T(Feem7$$\"3G+++5jyp`Fa[m$\"3tUCq$3@$[dFa[m$\"3'RP5Y/(*=_$Feem7$$\"37++D6EjpeFa[m$\"3 " 0 "" {MPLTEXT 1 0 14 "-Et*27.12*_eV;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%$_eVG$\")/>_C!\"'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 149 "For any atom the eigenenergy of the highest occupie d orbital should equal the negative of the ionization potential. Our r esult is indeed quite close." }}{PARA 0 "" 0 "" {TEXT -1 328 "There is another quantity that can be measured, namely the total energy of the atom (equal to the sum of ionization energies for both electrons). Th is is not simply twice the eigenenergy, since after ionizing on of the two He-electrons, the other is left in a hydrogen-like state (energy \+ of -2 a.u. due to Z=2, and E_1s=-Z^2/2)." }}{PARA 0 "" 0 "" {TEXT -1 97 "Combining this answer with the calculated 1s binding energy we hav e for the total binding energy:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "(Et-2)*27.12*_eV;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$%$_eVG$ !)/>wy!\"'" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 261 28 "The electrostatic \+ repulsion:" }}{PARA 0 "" 0 "" {TEXT -1 256 "Now we look at the questio n as to what potential is associated with the approximate wavefunction . For this we need a solution to the Poisson equation for a sphericall y symmetric charge distribution based on the properly normalized varia tional wavefunction." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "A1s :=1/sqrt(int(chi^2,r=0..infinity));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%$A1sG,$*$)%&beta|irG#\"\"$\"\"#\"\"\"F+" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 33 "rho:=subs(beta=beta,(A1s*chi)^2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$rhoG,$*()%&beta|irG\"\"$\"\"\")%\"rG\"\"#F*)-%$ex pG6#,$*&F(F*F,F*!\"\"F-F*\"\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "int(rho,r=0..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\" \"\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 116 "Note that the 4Pi from t he integration over theta and phi are cancelled by the square of the s pherical harmonic Y00!" }}{PARA 0 "" 0 "" {TEXT -1 91 "The solution to Poisson's equation in multipole expansion leads to the monopole expre ssion:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "V0:=unapply(simpl ify(int(rho,r=0..R)/R+int(rho/r,r=R..infinity)),R);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%#V0GR6#%\"RG6\"6$%)operatorG%&arrowGF(,$*&,(*(%&bet a|irG\"\"\"9$F1-%$expG6#,$*&F0F1F2F1!\"#F1F1F3F1F1!\"\"F1F2F9F9F(F(F( " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "plot([V(r)+2/r,subs(bet a=beta0,V0(r))],r=0..5,color=[red,blue],view=[0..5,0..2]);" }}{PARA 13 "" 1 "" {GLPLOT2D 693 299 299 {PLOTDATA 2 "6&-%'CURVESG6$7fn7$$\"3( *******R`!eS$!#?$\"35)QQy'Gz%p\"!#<7$$\"3&*******z1h6oF*$\"39[>G#f:Xp \"F-7$$\"3*******>gT<-\"!#>$\"3aPtP%QpTp\"F-7$$\"3++++O@Ki8F6$\"3q)\\q &Qav$p\"F-7$$\"3*******R?$[V?F6$\"3Ic')GH!HFp\"F-7$$\"3)******>FWYs#F6 $\"3GYO\"ffX9p\"F-7$$\"3)******zSmp3%F6$\"37-k1FL9)o\"F-7$$\"3/+++X&)G \\aF6$\"3g%f7!\\v\"Ro\"F-7$$\"3-+++3)QY*\\ \"F-7$$\"3G+++sZ3yTFZ$\"3yL:=uSZ'R\"F-7$$\"35+++]4\\Y_FZ$\"3k!oJK>)Q$H \"F-7$$\"3S+++U-/PiFZ$\"3#H)GXa1N,7F-7$$\"3s******empisFZ$\"3:6oHf]$=6 \"F-7$$\"3&******>*>VB$)FZ$\"3KE%[H**Ql-\"F-7$$\"3c*****R`l2Q*FZ$\"3a% Q2aRz?\\*FZ7$$\"3-+++0j$o/\"F-$\"3YS)*o&\\n\\x)FZ7$$\"3!******>&>jU6F- $\"3&)*RF^lO`?)FZ7$$\"3%******H;v/D\"F-$\"3,dzXMV&zi(FZ7$$\"3!******z= h(e8F-$\"3gLjWao\")4rFZ7$$\"35+++Q[6j9F-$\"3%4:xokREm'FZ7$$\"35+++\\z( yb\"F-$\"3gxq(>[QgH'FZ7$$\"3%******\\Xg0n\"F-$\"3[.QJD\")R.fFZ7$$\"3)* *****pJpW`(>F-$\"3r=e$)Q\"[R.&FZ7$$\"3#******4f#=$3#F-$\"39pu%4Cx2y%FZ7 $$\"3%)*****Hxpe=#F-$\"3MU_DAT7hXFZ7$$\"35+++uI,$H#F-$\"3Sgi3\"p-;N%FZ 7$$\"3=+++rSS\"R#F-$\"3I-%e2r'*[<%FZ7$$\"3-+++`?`(\\#F-$\"3))y)*>X&p#* *RFZ7$$\"3!********>pxg#F-$\"3E:V#)4#)[JQFZ7$$\"38+++g4t.FF-$\"31(*)f! )>&G'p$FZ7$$\"3*********Gst!GF-$\"3=6I\"orJ/c$FZ7$$\"38+++ERW9HF-$\"3E r$zURk+V$FZ7$$\"3@+++KE>>IF-$\"3$z_`NRf8J$FZ7$$\"3%******>RU07$F-$\"3I P?OOh,/KFZ7$$\"36+++?S2LKF-$\"3_Q7JdMl#4$FZ7$$\"3?+++$p)=MLF-$\"3-O[Yd C'*)*HFZ7$$\"3\"*******)=]@W$F-$\"3xJx*fMv\\!HFZ7$$\"3')******\\$z*RNF -$\"3x8INTAuCGFZ7$$\"3#)*****RYKpk$F-$\"3&f$3$e.Q>u#FZ7$$\"3))*****z+n vu$F-$\"3#*yD!4CK$oEFZ7$$\"3)******R5fF&QF-$\"3Oy*QCb'\\&f#FZ7$$\"3')* ****\\g.c&RF-$\"3\"yAT!pn-GDFZ7$$\"3K+++nAFjSF-$\"3MP,6G\"[5Y#FZ7$$\"3 q*****\\)*pp;%F-$\"3M#>;qZ4)*R#FZ7$$\"3#)*****z(e,tUF-$\"39?p,5kDSBFZ7 $$\"3G+++fO=yVF-$\"3;%yb1YWSG#FZ7$$\"3u*****f#>#[Z%F-$\"3C+l;A-sMAFZ7$ $\"3)******pG!e&e%F-$\"33hWY\\_u!=#FZ7$$\"3%)*****\\)Qk%o%F-$\"3v.a(\\ BJY8#FZ7$$\"3y*****>Mm-z%F-$\"3o**eryYc(3#FZ7$$\"3C+++60O\"*[F-$\"3!)f D:3'>W/#FZ7$$\"\"&\"\"!$\"3q.p&p0*****>FZ-%'COLOURG6&%$RGBG$\"*++++\"! \")$Fe]lFe]lF_^l-F$6$7fn7$F($\"33'=i**)eX(o\"F-7$F/$\"3w4_L%*eM(o\"F-7 $F4$\"3wuRH0T;(o\"F-7$F:$\"3]tQyg<\"po\"F-7$F?$\"3'*H=uTJm&*\\\"F-7$Fgo$\"3K4JW'y0sR\"F-7$F\\p$\"3AC!*)HU dVH\"F-7$Fap$\"3#ziOUd\\B?\"F-7$Ffp$\"34(G24_lF6\"F-7$F[q$\"3ul\"o/`^t -\"F-7$F`q$\"3i&o9>l)*))\\*FZ7$Feq$\"3It8\\V**\\!y)FZ7$Fjq$\"3A@'=u#o& )4#)FZ7$F_r$\"3*>5WKz*\\JwFZ7$Fdr$\"3af$zT7fD6(FZ7$Fir$\"3k\"30HOfZm'F Z7$F^s$\"3?2iCsWq(H'FZ7$Fcs$\"3!)))\\Dv+k/fFZ7$Fhs$\"3#\\'R#zhIVg&FZ7$ F]t$\"3W*4KG(R7)G&FZ7$Fbt$\"31f&Q2-#\\M]FZ7$Fgt$\"3W!*=xN)z6y%FZ7$F\\u $\"3kT?1%>u#FZ7$Fgy$\"3;!\\1Z5N$oEFZ7$F\\ z$\"3%)o*3Li)\\&f#FZ7$Faz$\"3xJ97w#G!GDFZ7$Ffz$\"3_*)zu2#\\5Y#FZ7$F[[l $\"3M&\\6%f-\")*R#FZ7$F`[l$\"3fe!oC(pDSBFZ7$Fe[l$\"3mxT\"e'[/%G#FZ7$Fj [l$\"3:&*Gr@0sMAFZ7$F_\\l$\"3wB#H7YX2=#FZ7$Fd\\l$\"3fF$y,RJY8#FZ7$Fi\\ l$\"3%*Qc3!zkv3#FZ7$F^]l$\"3D3><*o>W/#FZ7$Fc]l$\"3#Rs#[9\"*****>FZ-Fi] l6&F[^lF_^lF_^lF\\^l-%+AXESLABELSG6$Q\"r6\"Q!6\"-%%VIEWG6$;F_^lFc]l;F_ ^l$\"\"#Fe]l" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 300 "We recogn ize that the electronic repulsion in our helium atom of independent el ectrons is modelled after the simple hydrogen-like solution to the pro blem. The above potential shows how asymptotically the charge distribu tion of electron 1 screens one of the protons for electron 2 (located \+ at a large " }{TEXT 262 1 "r" }{TEXT -1 8 "-value)." }}{PARA 0 "" 0 " " {TEXT -1 216 "A sophisticated central-field or Hartree-Fock calculat ion take the electrostatic repulsion due to the numerically obtained c harge density and re-calculates the eigenenergy/eigenfunction until co nvergence is achieved." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 263 22 "Electronic excitations" }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 443 "We proceed to calculat e the energy levels for the 2s and 2p states. We simply assume that we can use the potential obtained for the ground state, and calculate th e energy spectrum for this potential. For the 2s-state we will not car ry out a variational calculation as the energy will not be guranteed t o be above the exact eigenenergy for the given potential. We repeat ou r trial-and-error procedure to find a radial function with one node at " }{TEXT 264 1 "r" }{TEXT -1 3 ">0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "IC:=phi(eta)=eta,D(phi)(eta)=1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#ICG6$/-%$phiG6##\"\"\"\"*++++\"F*/--%\"DG6#F(F)F+" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "Et:=-0.15768;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#EtG$!&od\"!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "SE:=-1/2*diff(phi(r),r$2)+(V(r)-Et)*phi(r);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#SEG,&-%%diffG6$-%$phiG6#%\"rG-%\"$G6$F,\" \"##!\"\"F0*&,(*&\"\"\"F6F,F2!\"#*&,&F6F6*&-%$expG6#,$F,$!$O$F7F6,&F6F 6*&$\"%l;!\"$F6F,F6F6F6F2F6F,F2F6$\"&od\"!\"&F6F6F)F6F6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "sol:=dsolve(\{SE,IC\},phi(r),numeri c,output=listprocedure):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "phir:=subs(sol,phi(r)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "plot(phir,0..20,-0.5..0.5);" }}{PARA 13 "" 1 "" {GLPLOT2D 694 299 299 {PLOTDATA 2 "6%-%'CURVESG6$7eo7$$\"+N@Ki8!#6$\"+hEpD8F*7$$\"+qUkCF F*$\"+Jp/!e#F*7$$\"+0k'p3%F*$\"+.l'ew$F*7$$\"+S&)G\\aF*$\"+yO&e)[F*7$$ \"+5G$R<)F*$\"+4lkQpF*7$$\"+3x&)*3\"!#5$\"+Hm\"FD7$$\"+#H[D:)FD$\"+\"z2*=8FD 7$$\"+;w)*=#*FD$\"+#>-q1\"FD7$$\"+%pU&G5!\"*$\"+g`a$*yF*7$$\"+Em=N6Fgo $\"+0&oT&\\F*7$$\"+e0$=C\"Fgo$\"+#\\tc#>F*7$$\"+'R\"=\\8Fgo$!+z>:\"F *7$$\"+LA`c9Fgo$!+Pw_=UF*7$$\"+qI)Qc\"Fgo$!+vhVIsF*7$$\"+3RBr;Fgo$!+Oz 0;5FD7$$\"+W^\"\\)=Fgo$!+LqRm:FD7$$\"+zjf)4#Fgo$!+!Q/Z#FD7$$\"+'4;[\\#Fgo$!+qGmBGFD7$$\"+!QZ**p#Fgo$!+7:nKJFD7$$ \"+j'y]!HFgo$!+*[UcQ$FD7$$\"+IdA+&Fgo$!+?!=.l$FD7$$\"+]Z/NaFgo$!+ny?3MFD7$$\"+]$fC&eFgo$!+rdQUJ FD7$$\"+'z6:B'Fgo$!+dj\"y)GFD7$$\"+<=C#o'Fgo$!+IMX#e#FD7$$\"+n#pS1(Fgo $!+7hvHBFD7$$\"+j`A3vFgo$!+/Fm\\?FD7$$\"+n(y8!zFgo$!+#G'*z\"=FD7$$\"+j .tK$)Fgo$!+\")Q)Qe\"FD7$$\"+)3zMu)Fgo$!+()=]\"Q\"FD7$$\"+#H_?<*Fgo$!+f ez\">\"FD7$$\"+!G;cc*Fgo$!+@qLO5FD7$$\"+4#G,***Fgo$!+ogDx))F*7$$\"+!o2 J/\"!\")$!+IzJGvF*7$$\"+%Q#\\\"3\"F\\[l$!+.C*>]'F*7$$\"+;*[H7\"F\\[l$! +m&GN`&F*7$$\"+qvxl6F\\[l$!+7e(*pYF*7$$\"+`qn27F\\[l$!+hzaWRF*7$$\"+cp @[7F\\[l$!+9NfTLF*7$$\"+3'HKH\"F\\[l$!+f?brFF*7$$\"+xanL8F\\[l$!+sv#oL #F*7$$\"+v+'oP\"F\\[l$!+!yKB%>F*7$$\"+S<*fT\"F\\[l$!+@?UQ;F*7$$\"+&)Hx e9F\\[l$!+>>+c8F*7$$\"+.o-*\\\"F\\[l$!+T*p38\"F*7$$\"+TO5T:F\\[l$!+4-e 6$*!#77$$\"+U9C#e\"F\\[l$!+c9zbwF[_l7$$\"+1*3`i\"F\\[l$!+X%[W='F[_l7$$ \"+$*zym;F\\[l$!+Oz)o(\\F[_l7$$\"+^j?4F\\[l$!+vwho>Fdal7$$\" +/Uac>F\\[l$\"+7fyxTFdal7$\"#?$\"+=AA&4\"F[_l-%'COLOURG6&%$RGBG$\"#5! \"\"$\"\"!F[clFjbl-%+AXESLABELSG6$Q!6\"F_cl-%%VIEWG6$;Fjbl$F`blF[cl;$! \"&Fibl$\"\"&Fibl" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "E 2s:=Et;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$E2sG$!&od\"!\"&" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 116 "Now let us carry out the calculat ion for the 2p state. First we carry out the hydrogen-like variational calculation." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "l:=1;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"lG\"\"\"" }}}{EXCHG {PARA 0 "" 0 " " {TEXT -1 146 "The trial function is hydrogen-like: the radial functi on picks up a factor of r^l, and we add the centrifugal potential to t he radial Hamiltonian:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "c hi:=r^(1+l)*exp(-beta*r);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$chiG*& )%\"rG\"\"#\"\"\"-%$expG6#,$*&%&beta|irGF)F'F)!\"\"F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 102 "E2pT:=Int(chi*(-1/2*Diff(chi,r$2)+ (V(r)+1/2*l*(l+1)/r^2)*chi),r=0..infinity)/Int(chi^2,r=0..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%E2pTG*&-%$IntG6$*()%\"rG\"\"#\"\" \"-%$expG6#,$*&%&beta|irGF-F+F-!\"\"F-,&-%%DiffG6$*&F*F-F.F--%\"$G6$F+ F,#F4F,*(,(*&F-F-F+F4!\"#*&,&F-F-*&-F/6#,$F+$!$O$FAF-,&F-F-*&$\"%l;!\" $F-F+F-F-F-F4F-F+F4F-*&F-F-*$F*F-F4F-F-F*F-F.F-F-F-/F+;\"\"!%)infinity GF--F'6$*&)F+\"\"%F-)F.F,F-FQF4" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "E2pT:=value(E2pT);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%E2pTG ,$*&*&,0*$)%&beta|irG\"\"'\"\"\"$\")+D1R\"\"!*&$\"*v=ny&F0F-)F+\"\"%F- F-*&$\"*+++]#F0F-)F+\"\"&F-F-*&$\"*++q\\(F0F-)F+\"\"$F-F-*&$\"+sI3L5F0 F-F+F-!\"\"*&$\"*+?N'HF0F-)F+\"\"#F-FC$\"*G\\wA&F0FCF-F+F-F-*$),&F+$\" #DF0$\"#UF0F-F:F-FC$\"++++]7!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "plot(E2pT,beta=0.45..0.6);" }}{PARA 13 "" 1 "" {GLPLOT2D 693 299 299 {PLOTDATA 2 "6%-%'CURVESG6$7S7$$\"35+++++++X!#=$!3+kqp*=!**[7F *7$$\"3#****\\7t&pKXF*$!3i;f!*=6\"4D\"F*7$$\"3;+v=7T9hXF*$!3K=SsM7]_7F *7$$\"3=+](=HPJf%F*$!3Wy2a+#*>a7F*7$$\"3/+]7VDMDYF*$!3s&omz&>\"eD\"F*7 $$\"3-+vVGZRdYF*$!3-5o$)G5Kd7F*7$$\"3')*\\(=276(o%F*$!3OUI4++D\"*309[F*$!3%G@'Q$y?LE\"F*7 $$\"3H+]i&e*yU[F*$!3Et)e`)Q]n7F*7$$\"3,+++&>0)H]F*$!3'4\\nG*f&yE\"F*7$$\"3++v=-p6j]F*$!3)f]LS Ys\"o7F*7$$\"3S++]2Mg#4&F*$!3[`O2IrOo7F*7$$\"3O+v=xZ&\\7&F*$!3g_vvT%)[ o7F*7$$\"3q*\\i:$4wb^F*$!3]\"G#3'G9&o7F*7$$\"3/+](=#R!z=&F*$!3QT++j![% o7F*7$$\"3Q+v$4A@u@&F*$!3Q5LG$[.$o7F*7$$\"3K**\\i:'f#\\_F*$!3eoZ]!pd!o 7F*7$$\"3?+vof2L#G&F*$!3/0IFAPqn7F*7$$\"3?+D\"yG>6J&F*$!3cCcg,QJn7F*7$ $\"3E+](oo6AM&F*$!33LWUDs!oE\"F*7$$\"3G++]xJLu`F*$!3aMo'Hv!>m7F*7$$\"3 t**\\P*yddS&F*$!3#Q8J28'\\l7F*7$$\"3W+v=xh7F*7$$\"3[*\\(o/Q*>c&F*$!3y3;['*4qg7F*7$$\"3k** *\\(Q(zSf&F*$!3m$eh%p,Wf7F*7$$\"3o*\\(=-,FCcF*$!3e\"*yCo!o\"e7F*7$$\"3 O*\\P4tFel&F*$!3#[yM4`\\nD\"F*7$$\"3;++D\"3\"o'o&F*$!3q*y-V$\\Fb7F*7$$ \"3Q+voz;)*=dF*$!3;E,#fIQOD\"F*7$$\"35+++&*44]dF*$!3kW/>)Hs>D\"F*7$$\" 3!)**\\7jZ!>y&F*$!31na " 0 " " {MPLTEXT 1 0 41 "beta0:=fsolve(diff(E2pT,beta),beta=0..2);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "E2pV:=subs(beta=beta0,E2pT);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&beta0G$\"+hBY\\^!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%E2pVG$!+9h^o7!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 66 "Now we verify the variational calculation by a numerical soluti on:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "IC:=phi(eta)=eta^(l+ 1),D(phi)(eta)=(l+1)*eta^l;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#ICG6 $/-%$phiG6##\"\"\"\"*++++\"#F+\"2++++++++\"/--%\"DG6#F(F)#F+\")+++]" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "Et:=-0.12699;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#EtG$!&*p7!\"&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "SE:=-1/2*diff(phi(r),r$2)+(V(r)+1/2*l*(l+1)/r^2-Et)*p hi(r);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#SEG,&-%%diffG6$-%$phiG6#% \"rG-%\"$G6$F,\"\"##!\"\"F0*&,**&\"\"\"F6F,F2!\"#*&,&F6F6*&-%$expG6#,$ F,$!$O$F7F6,&F6F6*&$\"%l;!\"$F6F,F6F6F6F2F6F,F2F6*&F6F6*$)F,F0F6F2F6$ \"&*p7!\"&F6F6F)F6F6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "sol :=dsolve(\{SE,IC\},phi(r),numeric,output=listprocedure):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "phir:=subs(sol,phi(r)):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "plot(phir,0..25);" }}{PARA 13 "" 1 "" {GLPLOT2D 695 243 243 {PLOTDATA 2 "6%-%'CURVESG6$7io7$$\"+pE!Hq\"!# 6$\"+E&)H$o#!#87$$\"+Q`!eS$F*$\"+[3Wb5!#77$$\"+2!3(3^F*$\"+`)[bL#F37$$ \"+v1h6oF*$\"+i>>\"3%F37$$\"+,;u@5!#5$\"+NKd&*))F37$$\"+N@Ki8FA$\"+P(* oI:F*7$$\"+-K[V?FA$\"+zznRKF*7$$\"+qUkCFFA$\"+z%*fAaF*7$$\"+0k'p3%FA$ \"+e@<\"4\"FA7$$\"+T&)G\\aFA$\"+Vk,YyFA$\"+i8KXIFA7$$\"+bhL0!*FA$\"+.lBOPFA7$$\"+O&o!>5!\"*$\"+pm5OW FA7$$\"+_MP_6Fio$\"+uL?L_FA7$$\"+n$ycG\"Fio$\"+$4c'>gFA7$$\"+#G$)*=9Fi o$\"+!Q![*y'FA7$$\"+)>)G_:Fio$\"+ONWHvFA7$$\"+WnZ'o\"Fio$\"+,]\\a#)FA7 $$\"+\"Hl1#=Fio$\"+_?xP*)FA7$$\"+QQ&[&>Fio$\"+]yu!e*FA7$$\"+&QU!*3#Fio $\"+U_G<5Fio7$$\"+HR9cBFio$\"+P@8Fio7$$\"+CU$\\P$Fio$\"+p oSc8Fio7$$\"+G$[8j$Fio$\"+p25v8Fio7$$\"+'GWwp$Fio$\"+hZ,z8Fio7$$\"+X-% Rw$Fio$\"+\"Qj5Q\"Fio7$$\"+.iBIQFio$\"+g^G#Q\"Fio7$$\"+i@`'*QFio$\"+?e q#Q\"Fio7$$\"+@\"GG'RFio$\"+k6N#Q\"Fio7$$\"+zS7HSFio$\"+4qC\"Q\"Fio7$$ \"+Q+U&4%Fio$\"+R\">%z8Fio7$$\"+&*frhTFio$\"+3(zeP\"Fio7$$\"+!Q\\gU%Fi o$\"+^gQg8Fio7$$\"+lFQ!p%Fio$\"+EOyL8Fio7$$\"+?:=M_Fio$\"+-cue7Fio7$$ \"+g(fJr&Fio$\"+G@D%FA7$$\"+&H\"fT5!\")$\"+)pg&*f$FA7$$ \"+')[$H4\"Fiz$\"+hhYcIFA7$$\"+Ol]Y6Fiz$\"+U#**\\c#FA7$$\"+N?q&>\"Fiz$ \"+oOFv@FA7$$\"+Egw[7Fiz$\"+s9<9=FA7$$\"+*f%)QI\"Fiz$\"+$41n\\\"FA7$$ \"+![l=N\"Fiz$\"+$pSAE\"FA7$$\"+Xho.9Fiz$\"+Kt4Z5FA7$$\"+i>Ad9Fiz$\"+* zgug)F*7$$\"+;jf4:Fiz$\"+iVQ'3(F*7$$\"+&>r-c\"Fiz$\"+]JbceF*7$$\"+4q`; ;Fiz$\"+FUAEZF*7$$\"+YV4n;Fiz$\"+2!3\"))QF*7$$\"+%4v5s\"Fiz$\"+x&**y9$ F*7$$\"+u'*)*pFiz$\"+#***Rh8F*7$$\"+-=!y(>Fiz$\"+Le$z3 \"F*7$$\"+LhjJ?Fiz$\"+cW=.&)F37$$\"+#*\\[$3#Fiz$\"+@KUBFiz$\"+0jZ(3\"F-7$$\"+qJ8&R#Fiz $!+<'*e)*)*F-7$$\"+b-oXCFiz$!+Mr'e3#F37$\"#D$!+#H=*pLF3-%'COLOURG6&%$R GBG$\"#5!\"\"$\"\"!F_dlF^dl-%+AXESLABELSG6$Q!6\"Fcdl-%%VIEWG6$;F^dl$Fd clF_dl%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "E2p:=Et;" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$E2pG$!&*p7!\"&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 115 "The answer is just slightly below the value obta ined by the variational calculation with a hydrogenic wavefunction." } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "(E2p-E1s)*27.12*_eV;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#,$%$_eVG$\"*_$z2@!\"(" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 122 "We learn that the model predicts an exci tation energy of slightly more than 21 eV. What does this have to do w ith reality?" }}{PARA 0 "" 0 "" {TEXT -1 318 "Unfortunately the connec tion is not very direct: inspection of a spectroscopic table for the H e atom reveals that it is complicated, and that many levels correspond to the primitive (1s 2p) configuration: this is the result of complic ations arising from coupling orbital and spin angular momenta for the \+ two electrons." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 418 "Nevertheless, if we are generous, we can compare the cal culated excitation energy with the average result which in some sense \+ ignores the spin-orbit interactions (these are also called magnetic or fine structure interactions). The model does provide a realistic asse ssment for the average transition energy (which is in the UV regime). \+ For an understanding of the visible spectrum one cannot ignore the fin e structure." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 98 " One important piece of information that does emerge from this \+ calculation is the lifting of the (" }{TEXT 265 2 "nl" }{TEXT -1 196 " )-degeneracy observed in the hydrogen atom spectrum: the 2s and 2p sta tes are no longer degenerate, and thus, there is the possibility of ob serving photons that correspond to 2p-2s de-excitations." }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 266 21 "Matrix re presentation" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 210 "We proceed with the calculation of approximate spectra in an a ngular momentum symmetry sector. For this purpose we first define ours elves a suitable basis. We defined so-called Slater-type orbitals for \+ a given " }{TEXT 267 1 "l" }{TEXT -1 66 "-symmetry, and use a Gram-Sch midt procedure to orthogonalize them." }}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 103 "Our radial functions are real-valued , and therefore we have the simple inner product and normalization:" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "IP:=(f,g)->int(f*g,r=0..in finity):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "NO:=xi->xi/sqrt (IP(xi,xi));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#NOGR6#%#xiG6\"6$%)o peratorG%&arrowGF(*&9$\"\"\"-%%sqrtG6#-%#IPG6$F-F-!\"\"F(F(F(" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "STO:=(n,l,beta)->r^(n+l)*exp (-beta*r);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$STOGR6%%\"nG%\"lG%%be taG6\"6$%)operatorG%&arrowGF**&)%\"rG,&9$\"\"\"9%F3F3-%$expG6#,$*&9&F3 F0F3!\"\"F3F*F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "L:=0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "N:=8;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"NG\" \")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "B1:=[seq(STO(n,L,17/ 10),n=1..N)];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#B1G7**&%\"rG\"\"\" -%$expG6#,$F'#!#<\"#5F(*&)F'\"\"#F(F)F(*&)F'\"\"$F(F)F(*&)F'\"\"%F(F)F (*&)F'\"\"&F(F)F(*&)F'\"\"'F(F)F(*&)F'\"\"(F(F)F(*&)F'\"\")F(F)F(" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 58 "These states are linearly independ ent, but not orthogonal!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "IP(B1[1],B1[2]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"%]P\"&@N)" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 119 "The Gram-Schmidt procedure takes \+ a list of functions (state vectors) and orthonormalizes them using the two procedures." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "GS:=pro c(vecs) local i,n,j,res,xi;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "n:=n ops(vecs); res:=[NO(vecs[1])];" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "f or i from 2 to n do:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "xi:=vecs[i] -add(IP(vecs[i],res[j])*res[j],j=1..i-1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "res:=[op(res),NO(xi)]; od: end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "B1ON:=GS(B1):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "IP(B1ON[5],B1ON[5]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "IP(B1ON[4],B1O N[5]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 8 "B1ON[2];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$* &,&*&)%\"rG\"\"#\"\"\"-%$expG6#,$F(#!#<\"#5F*F**&#\"#:\"# " 0 "" {MPLTEXT 1 0 65 "plot([seq(B1ON[i],i=1..4)],r=0..10,color=[red,blue,gr een,black]);" }}{PARA 13 "" 1 "" {GLPLOT2D 693 299 299 {PLOTDATA 2 "6( -%'CURVESG6$7do7$$\"\"!F)F(7$$\"3gmmTN@Ki8!#>$\"3$yBJ!)z$*4!fF-7$$\"3A LL$3FWYs#F-$\"3%3m.;0!=`6!#=7$$\"3%)***\\iSmp3%F-$\"3_$Hs'Q)p,p\"F57$$ \"3WmmmT&)G\\aF-$\"3`\\U'HCo>?#F57$$\"3m****\\7G$R<)F-$\"3gp4(\\9_M:$F 57$$\"3GLLL3x&)*3\"F5$\"3&4eUEf\"H9SF57$$\"3YmmTN@Ki8F5$\"3)*ocn&)Gu!z %F57$$\"3$*****\\ilyM;F5$\"3#p;PXn\"o)[&F57$$\"3RLLe*)4D2>F5$\"3miF^,Y i8hF57$$\"3emmm;arz@F5$\"3%4CKDw\\2n'F57$$\"3;L$e*)4bQl#F5$\"3G)[')[AP G\\(F57$$\"3v***\\7y%*z7$F5$\"3'GCeqU$fZ\")F57$$\"3Lm;ajW8-OF5$\"3u'[j Gy/gl)F57$$\"3[LL$e9ui2%F5$\"3v%[#4&z5o.*F57$$\"3\"pm\"H2Q\\4YF5$\"3MX l$ytYLL*F57$$\"3z***\\(oMrU^F5$\"30P0!R[?1^*F57$$\"3Cm\"z%*HB$4aF5$\"3 qv!\\Tt!Rg&*F57$$\"3oK$3-8Lfn&F5$\"3+&e!=Kc2(e*F57$$\"3B+v$4'HaUfF5$\" 3#3;1MOBEf*F57$$\"3nmmm\"z_\"4iF5$\"3K*>knya)y&*F57$$\"3Unmmm6m#G(F5$ \"33th0pc#3O*F57$$\"39ommT&phN)F5$\"3K\"GXZr\"**[*)F57$$\"3A,+v=ddC%*F 5$\"3$GR3(fE#oT)F57$$\"3KLLe*=)H\\5!#<$\"3N>!4X#[c9yF57$$\"3-++v=JN[6F \\s$\"3\")*yH[a\\oA(F57$$\"3smm\"z/3uC\"F\\s$\"3+Y'\\+r2Oj'F57$$\"3ULL e*ot*\\8F\\s$\"3?P)R^bY..'F57$$\"3!****\\7LRDX\"F\\s$\"3sU))G)H-.X&F57 $$\"3[LLekGhe:F\\s$\"3'G0F'p!QL)[F57$$\"3%om;zR'ok;F\\s$\"3K\"R(=]n3bV F57$$\"3OLL3_(>/x\"F\\s$\"3'\\hwd\"HqpQF57$$\"33++D1J:w=F\\s$\"3iq*31A _hU$F57$$\"3oLLL3En$4#F\\s$\"3T#\\D\"[1^TEF57$$\"3#pmmT!RE&G#F\\s$\"3i XfEmru\"3#F57$$\"3D+++D.&4]#F\\s$\"3]EW0XB!*y:F57$$\"3;+++vB_*Q'QF-7$$\"33++D\"o7Tv$F\\s$\"3rV$=,YWa\"GF-7$$\"3?LLL$Q*o ]RF\\s$\"3!G^$oP8>@@F-7$$\"3m++D\"=lj;%F\\s$\"3#[QK7W`.b\"F-7$$\"3S++v V&R$G0uSiFhy7$$\"3emm;/T1&*\\F\\s$\"3xt#R\"*eZ Ma%Fhy7$$\"3=nm\"zRQb@&F\\s$\"3'yQ,b(Q7hKFhy7$$\"3:++v=>Y2aF\\s$\"3W7z e]p&)RCFhy7$$\"3Znm;zXu9cF\\s$\"3!pjE*e-+\"y\"Fhy7$$\"34+++]y))GeF\\s$ \"3yhy@pou%G\"Fhy7$$\"3H++]i_QQgF\\s$\"3&y'R$yo)Q@$*!#@7$$\"3b++D\"y%3 TiF\\s$\"3*)y#)Qfl+EoF\\\\l7$$\"3+++]P![hY'F\\s$\"3\"\\*)Hnv%yB[F\\\\l 7$$\"3iKLL$Qx$omF\\s$\"37R4%3p8u_$F\\\\l7$$\"3Y+++v.I%)oF\\s$\"3aBVDnL yADF\\\\l7$$\"3?mm\"zpe*zqF\\s$\"3c]>9]mNg=F\\\\l7$$\"3;,++D\\'QH(F\\s $\"3]!pt^#=GK8F\\\\l7$$\"3%HL$e9S8&\\(F\\s$\"3zB0.dpXB(*!#A7$$\"3s++D1 #=bq(F\\s$\"3O]!4$*yI1*pF`^l7$$\"3\"HLL$3s?6zF\\s$\"31Hkx-gPf]F`^l7$$ \"3a***\\7`Wl7)F\\s$\"3l<\"4=vIRg$F`^l7$$\"3enmmm*RRL)F\\s$\"3'4\"eVh+ x(f#F`^l7$$\"3%zmmTvJga)F\\s$\"3C!*)>en$\\d=F`^l7$$\"3]MLe9tOc()F\\s$ \"3*Qstt![0J8F`^l7$$\"31,++]Qk\\*)F\\s$\"3!p\"R:&R$[%z*!#B7$$\"3![LL3d g6<*F\\s$\"3Cb'yhg\"Q()oFd`l7$$\"3%ymmmw(Gp$*F\\s$\"3X`@7Q$3T-&Fd`l7$$ \"3C++D\"oK0e*F\\s$\"3;k')fg,S(e$Fd`l7$$\"35,+v=5s#y*F\\s$\"3Qx,-v1h(f #Fd`l7$$\"#5F)$\"3k`6\\^#e_$=Fd`l-%'COLOURG6&%$RGBG$\"*++++\"!\")F(F(- F$6$7]pF'7$F+$!3Bs\"Qk^,j+\"F57$F1$!3;I[oe\"*oN>F57$F7$!3sCaV^J'=z#F57 $F<$!3Oa$olWy$yNF57$FA$!3uV(>grcf&\\F57$FF$!3;*QWo;ZT4'F57$FK$!3oCrKxr l;qF57$FP$!3U`8E9&=`u(F57$FU$!3uQEK3&=-I)F57$FZ$!3/%*HwYF#)*p)F57$$\"3 ()*\\7yD&y;CF5$!3&p$ozlkSM*)F57$Fin$!3qfp7Aheu!*F57$$\"3WmT5S\\#4*GF5$ !3+b/*Q]9*H\"*F57$F^o$!3\"*oa>9)R#4\"*F57$Fco$!3=m(zk6;?())F57$Fho$!3X :c&pPS7U)F57$F]p$!3JZe\"fX]1s(F57$Fbp$!3kt6>x:\"=(oF57$F\\q$!3q\"HJk?p N#fF57$Ffq$!3\"HN[Q$f'e\"\\F57$$\"3/nm;zp!fu'F5$!3/v%o-?jT(QF57$F[r$!3 *R/AUky8$GF57$$\"3ynm;a`T>yF5$!3u='R'yCI3=F57$F`r$!3J:,:HC,5#)F-7$$\"3 7M$3-js.*))F5$\"3!RQ>rQL39\"F-7$Fer$\"3!=]Zp3v0$**F-7$$\"3?n;H2)y(e**F 5$\"3x\"GzQ**H*4=F57$Fjr$\"3!*p2R9n#4c#F57$F`s$\"3[3NsfBatPF57$Fes$\"3 '*zunVBk`ZF57$Fjs$\"3!eQK!fyYNbF57$F_t$\"3oR3%))oj.5'F57$Fdt$\"3'\\A_1 k#f#['F57$Fit$\"3)o5gxPl\")o'F57$$\"3YL$ekt>6p\"F\\s$\"3RfZ3qrG:nF57$$ \"34+++vIb<F\\s$\"3<+;dI1$3_'F57$Fhu$\"33VY(*=\"45G'F57 $F]v$\"3_@5F#fEHt&F57$Fbv$\"3roTQEPj;]F57$Fgv$\"3]*o()3:NoF%F57$F\\w$ \"3SS>?G'Q00W#F57$F[x$\"3R%R G_aE(4?F57$F`x$\"3m&\\G(GJ8(e\"F57$Fex$\"3Q%=xuE=wF\"F57$Fjx$\"3E@1*G0 KV***F-7$F_y$\"3#GH)R:_**eyF-7$Fdy$\"35JNn_()*)ygF-7$Fjy$\"3Sn?cc#f#yZ F-7$F_z$\"35GM!QwA!oOF-7$Fdz$\"3_<\"Q/<8Rx#F-7$Fiz$\"3[R7I*Rjs;#F-7$F^ [l$\"3C8K#QW&[a;F-7$Fc[l$\"3-Iat1'*[Z7F-7$Fh[l$\"32(Q>s\\0WV*Fhy7$F^\\ l$\"3g^d-uJP!=(Fhy7$Fc\\l$\"3/zxFp5K(G&Fhy7$Fh\\l$\"3Hu[1GfS1SFhy7$F]] l$\"3Kt.(p!\\GsHFhy7$Fb]l$\"3/*4(=T3GjAFhy7$Fg]l$\"3[#ochBvnn\"Fhy7$F \\^l$\"3m#=mrL&=i7Fhy7$Fb^l$\"3#z+ke0>JO*F\\\\l7$Fg^l$\"3qX;&H.82)pF\\ \\l7$F\\_l$\"3*y/)og7*[7&F\\\\l7$Fa_l$\"3In?Ml:')*z$F\\\\l7$Ff_l$\"3Y$ HITojVz#F\\\\l7$F[`l$\"3#=&e0q$)F`^l7$F`al$\"3zF*)=\"o(GDh F`^l7$Feal$\"3*G\\G^Ou$QXF`^l7$Fjal$\"3)pyX6A>ZG$F`^l-F_bl6&FablF(F(Fb bl-F$6$7\\qF'7$F+$\"3y\"*z=`l1,9F57$F1$\"3YZ%)H&4tAl#F57$F7$\"3akg'pZ* =jPF57$F<$\"3;H\")z;K\"Hu%F57$$\"3MKL3x1h6oF-$\"3TcLIo%*3+cF57$FA$\"3M Z*fbAIHM'F57$$\"3±z%\\DO&*F-$\"3wcee%zK#zpF57$FF$\"31s0ZcvR;vF57$FK $\"3z'z#G!QK5K)F57$FP$\"34FNb4?t3))F57$$\"3!omTgx=5x\"F5$\"3'Q3E\"\\/[ [*)F57$FU$\"3Tg%oQhpg-*F57$$\"3)***\\7.K[V?F5$\"3!fy^qOkl/*F57$FZ$\"3' >fKl5_Z,*F57$Fin$\"3aj'Q/@<%e&)F57$F^o$\"3%3bi@*=[)o(F57$Fco$\"3y9#z]d A8a'F57$Fho$\"3:q`\\_rZE_F57$$\"3Z+DcwR)GM%F5$\"3G$eaFU**oW%F57$F]p$\" 3+\\]=)zB[l$F57$$\"3OL3-QO5w[F5$\"39;U;D6ChGF57$Fbp$\"3'y03&eDfv?F57$F gp$\"3])yP)4@#fI\"F57$F\\q$\"3-;,6#F57$F[r$!30S#Rtn#y0KF57$Fagl$!3$QA\"* H:<67%F57$F`r$!3+w>=ugX`[F57$Figl$!31x0!eq[SS&F57$Fer$!3WEJwC)yYy&F57$ $\"3wM3-jsn\"p*F5$!3^9O?UMh9fF57$Fahl$!3s\\d\\Y\\Q1gF57$$\"3'*\\i:N!)e A5F\\s$!3e>!G>\"*)phgF57$Fjr$!3)G3N/OTB3'F57$$\"3/+](=#>1u5F\\s$!3KJuN MK6sgF57$$\"3wmm;ac#))4\"F\\s$!3CVmc5#p^.'F57$$\"3HL$ekQ*eB6F\\s$!3%** )Q$e4GI(fF57$F`s$!3k*)3$Q'f@()eF57$$\"3CLLL$e!)y>\"F\\s$!3UJ,O``m]cF57 $Fes$!3?t9#)*Q\\vL&F57$Fjs$!3WAw5'G(R&\\%F57$F_t$!3'Gh'pM*p*pMF57$Fdt$ !3Sw&4MZ![/BF57$Fit$!3DN\"H)*=Sk5\"F57$F^u$\"3Ydw@c\\,ohFhy7$Fcu$\"3-t rB]GUg6F57$Fe[m$\"3)\\'Q;zBc'=#F57$Fhu$\"3as[U5j(o3$F57$$\"3_++Dc#o%*= #F\\s$\"3+@b[kj6nPF57$F]v$\"3A#\\@46)[RVF57$Fbv$\"3us&f$*3AmC&F57$$\"3 ?+++]jB4EF\\s$\"3%RdYToo)=bF57$Fgv$\"3y1o_GKg$o&F57$$\"3!)*******>*ppF F\\s$\"3q&=URuR$GdF57$$\"3!******\\-w=#GF\\s$\"3/r5U\"GIDv&F57$$\"3)** *****\\G0uGF\\s$\"3?1kiWtjddF57$F\\w$\"3#*3o_&f5^u&F57$$\"3pL$eky#*4-$ F\\s$\"37b1uq.)=o&F57$Faw$\"39_*)=naStbF57$$\"3>+]7.%Q%GKF\\s$\"3QJhW' *H@'R&F57$Ffw$\"3w+YecX\"z<&F57$F[x$\"3*G'3O<=aVZF57$F`x$\"3scE2$pUU=% F57$Fex$\"3o_PG/n/zOF57$Fjx$\"3rz`2wvwVJF57$F_y$\"3w,H)z'3]qEF57$Fdy$ \"3#>z%Ho:NDAF57$Fjy$\"3'3K**eBvU'=F57$F_z$\"3$o\"zBv5@E:F57$Fdz$\"3cF &=<4s(G7F57$Fiz$\"3-l3rnSq55F57$F^[l$\"3%)p`&f&R%H8)F-7$Fc[l$\"3EOgDRH 'pX'F-7$Fh[l$\"3-bvYTDAB^F-7$F^\\l$\"3sJ_7BH!e2%F-7$Fc\\l$\"3o!R;*pQ`X JF-7$Fh\\l$\"3=$)y\")\\\"o;[#F-7$F]]l$\"3;kGbi>$)=>F-7$Fb]l$\"3rOZ^\"F-7$Fg]l$\"39y:x#3$fl6F-7$F\\^l$\"3aw'))zP'>\"3*Fhy7$Fb^l$\"371]] C*\\Z(pFhy7$Fg^l$\"3mA5.'))*et`Fhy7$F\\_l$\"3QN+*G5d$ySFhy7$Fa_l$\"3oa `BF\\4>JFhy7$Ff_l$\"3'3$eSvTIlBFhy7$F[`l$\"3K0>I>hq$z\"Fhy7$F``l$\"3g! Gpn3p%)Q\"Fhy7$Ff`l$\"3=[EPuI7L5Fhy7$F[al$\"3$Hjs]G)\\;zF\\\\l7$F`al$ \"3!>F%\\K5Y\\fF\\\\l7$Feal$\"3*)F57$$\"3x*\\(=U#Q/V\"F5$!3e*3V2BBj (*)F57$$\"31L$e*[Vb)\\\"F5$!3Z&G.\\[/w+*F57$$\"3km\"HdXqmc\"F5$!3J&fe: Pd^,*F57$FP$!3pA$egHw-+*F57$FU$!3$4^tAMj2u)F57$FZ$!3SZ#**Q_\\Z@)F57$Fg dl$!3/(eA;p(y*e(F57$Fin$!3AXJn@MPZoF57$F_el$!357?]$\\-\">gF57$F^o$!3I= R_yISK^F57$$\"3/LeRAY1lLF5$!3%*)=PN2X4@%F57$Fco$!3Soz7Lc,vKF57$$\"3=+v o/V?RQF5$!3Bz8M&Qj=M#F57$Fho$!36?X0ZO,E9F57$Fbgm$!3gP(3Q?PnJ%F-7$F]p$ \"3Yyte3(Rh6&F-7$Fjgm$\"3_s@Zz_&HR\"F57$Fbp$\"3E*ob$[C6/AF57$F\\q$\"3] C%Qe#f9%f$F57$Ffq$\"3)>8I#Ha`cYF57$Fifl$\"31wQpA\\Z$R&F57$F[r$\"3p@\\o W#3B\"eF57$$\"3nm;a8(\\oT(F5$\"3A?y?M/IqeF57$$\"3/nmTg#Q5b(F5$\"3%HL\" =v#y0\"fF57$$\"3Sn;H2oA&o(F5$\"3k?;^LRnLfF57$Fagl$\"3B^`=y.9SfF57$$\"3 9o;/,Rg`zF5$\"3SiHn^#[0$fF57$$\"3_om\"zW#z(3)F5$\"38\"Hg*[H[0fF57$$\"3 ))o;z%*4)>A)F5$\"3![DHY*=aleF57$F`r$\"3cZ$py+K8\"eF57$Figl$\"3WV/J;=Dm aF57$Fer$\"30y\"finQR%\\F57$Fahl$\"3[3[e;.)HG%F57$Fjr$\"3\"\\*)fWRn,_$ F57$Fe[n$\"3if)HEa>F57$Fb\\n$\"3gDyFfZAO6F57$F es$\"3CpaI\\D'[I$F-7$$\"31++vo3p)H\"F\\s$!3]1P)z^Z,y%F-7$Fjs$!3!Q&H:Tu 7X7F57$$\"3xmmT5lD,9F\\s$!3A)fv$fjUe>F57$F_t$!3jjZK?.@3EF57$Fdt$!3]'=E #*zlAs$F57$Fit$!3V%y>k38:]%F57$F`jl$!3c)e+C?jHw%F57$F^u$!3yH`'4$R9U\\F 57$$\"3)***\\i!4`oz\"F\\s$!3Jw*yd.H=+&F57$F][m$!310(H'[y?U]F57$$\"3CL$ 3xw>(\\=F\\s$!3KP<\"\\8?Q1&F57$Fcu$!3'egN.X^s1&F57$$\"3l;/,WIM.>F\\s$! 3ap(R#))3Z_]F57$$\"3VL3x\")H`I>F\\s$!3nCwQyN$)>]F57$$\"3A]7`>Hsd>F\\s$ !3?a$>if*3q\\F57$Fe[m$!3PL_;4_,/\\F57$$\"37+D\"Gt#HR?F\\s$!3*=*[ov\"zg s%F57$Fhu$!3G)\\,&[Jc#\\%F57$Fb^n$!354X/y)zi'RF57$F]v$!3%HjG*Ra5CLF57$ $\"3eLLe9r5$R#F\\s$!3&p.B[/N^]#F57$Fbv$!3!QM[y=Ysi\"F57$F]_n$!3k)\\I>? FjD(F-7$Fgv$\"3\")eR%p/#>K;F-7$Fj_n$\"3jkjub'*)*G)*F-7$F\\w$\"3/XixT(Q yu\"F57$Fg`n$\"3=u&=Oq&y$Q#F57$Faw$\"3#*30.*pen&HF57$Ffw$\"3CsAW0pZYSF 57$F[x$\"3')3h!Qm6rm%F57$$\"3mm;HdO2VOF\\s$\"3v(f)ev%pz!\\F57$F`x$\"3Y $41>%z8o]F57$$\"39L3xcoD.QF\\s$\"3J+b+YT=:^F57$$\"3km;HK5S_QF\\s$\"3kO hPc@d[^F57$$\"39+D\"y?X:!RF\\s$\"3ar#=QS3!p^F57$Fex$\"3cSLW.)3s<&F57$$ \"39n;H#GF&eSF\\s$\"3G'zAOcZe:&F57$Fjx$\"3uO=\\MCf'3&F57$$\"33++]iB0pU F\\s$\"3AeJi(o^G)\\F57$F_y$\"3)Gc0TD*Q[[F57$Fdy$\"3]%H<^mgJ\\%F57$Fjy$ \"3')R0Y8k74TF57$F_z$\"39UC\")4B62l sI4s\"F57$F^\\l$\"3B#)QUeh*RW\"F57$Fc\\l$\"3g?2A%y?%y6F57$Fh\\l$\"39*f \")=$ex\\(*F-7$F]]l$\"3l\\())z2O9\"zF-7$Fb]l$\"3(*)yr\\*3]6lF-7$Fg]l$ \"3lP@G.KGM_F-7$F\\^l$\"3dnYs#>S]Fhy7$F`al$\"3%z`?-ZtP)QFhy7$Feal$\"3t?,5 ,y'3.$Fhy7$Fjal$\"3U?S(zx@gJ#Fhy-F_bl6&FablF)F)F)-%+AXESLABELSG6$Q\"r6 \"Q!6\"-%%VIEWG6$;F(Fjal%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 43.000000 43.000000 0 0 "Curve 1" "Curve 2" "Curve 3" "Curve 4" }}}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 162 "Now we recycle the code from Matr ixRep.mws to generate the Hamiltonian matrix: we define procedures to \+ calculate the kinetic and potential energy matrix elements:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "Tkin:=(phi,psi)->-1/2*int(ex pand(phi*diff(psi,r$2)),r=0..infinity);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%TkinGR6$%$phiG%$psiG6\"6$%)operatorG%&arrowGF),$-%$intG6$-%'e xpandG6#*&9$\"\"\"-%%diffG6$9%-%\"$G6$%\"rG\"\"#F6/F>;\"\"!%)infinityG #!\"\"F?F)F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 75 "Vpot:=(ph i,psi)->int(expand(phi*psi*(V(r)+L*(L+1)/(2*r^2))),r=0..infinity);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%%VpotGR6$%$phiG%$psiG6\"6$%)operator G%&arrowGF)-%$intG6$-%'expandG6#*(9$\"\"\"9%F5,&-%\"VG6#%\"rGF5*&*(#F5 \"\"#F5%\"LGF5,&F@F5F5F5F5F5*$)F;F?F5!\"\"F5F5/F;;\"\"!%)infinityGF)F) F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "with(LinearAlgebra): \+ Digits:=15:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "HM:=Matrix(N ):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 114 "We make use of the symmetr y of the hamiltonian matrix: it allows to save almost a factor of 2 in computation time." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 106 "for \+ i from 1 to N do: for j from 1 to i do: HM[i,j]:=Tkin(B1ON[i],B1ON[j]) +Vpot(B1ON[i],B1ON[j]); od: od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 71 "for i from 1 to N do: for j from i+1 to N do: HM[i,j] :=HM[j,i]: od: od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "HMf:= map(evalf,HM):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "evals:=Eigenvalue s(HMf, output='list'):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "ev_s:=sor t(map(Re,evals));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%ev_sG7*$!3NOe` Vu*>/*!#=$!33&\\3&3Ch^:F($\"3%emF;NE\"pG!#>$\"3E0m`%yZ/!QF($\"3mn;#y1I h4\"!#<$\"3zU-AM/B3FF2$\"3.12M.UyUtF2$\"3#)*=&)3zU>K$!#;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 134 "We see that the chosen basis set can rep roduce the 1s and 2s states, but that it does not have a prediction fo r a bound 3s eigenstate." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT 268 11 "Exercise 1:" }}{PARA 0 "" 0 "" {TEXT -1 121 "Incr ease the matrix size and observe the stability of the 1s and 2s eigenv alues. Do you find an acceptable 3s eigenvalue?" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 269 11 "Exercise 2:" }}{PARA 0 " " 0 "" {TEXT -1 158 "Change the value of the parameter that controls t he Slater type orbital (STO) basis from the chosen value of 17/10. Fin d the best basis for the 3s eigenvalue." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 123 "Now the eigenfunctions. We \+ need a sorting procedure to arrange the result from the eigenvector ca lculation in proper order." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "VE:=Eigenvectors(HMf,output='list'):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "Vp:=[seq([Re(VE[i][1]),VE[i][2],map(Re,VE[i][3])],i=1 ..nops(VE))]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 344 "Min:=proc (x,y); if type(x,numeric) and type(y,numeric) then if x<=y then RETURN (true): else RETURN(false): fi; elif type(x,list) and type(y,list) and type(x[1],numeric) and type(y[1],numeric) then if x[1]<=y[1] then RET URN(true): else RETURN(false): fi; elif convert(x,string)<=convert(y,s tring) then RETURN(true): else RETURN(false): fi: end:\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "VEs:=sort(Vp,Min):\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 252 "Suppose that we would like to see the ei genfunctions corresponding to the four lowest-lying eigenvalues. The e igenvector for a given eigenvalue contains the expansion coefficients \+ for the expansion of the eigenstate in terms of the chosen basis state s." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "for i from 1 to 4 do: " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "psi0:=add(VEs[i][3][j]*B1ON[j], j=1..N):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "No:=1/sqrt(int(expand(p si0^2),r=0..infinity));" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "phi_a[i] :=add(No*VEs[i][3][j]*B1ON[j],j=1..N): od:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 66 "plot([seq(phi_a[i],i=1..4)],r=0..15,color=[red,blue ,green,black]);" }}{PARA 13 "" 1 "" {GLPLOT2D 695 232 232 {PLOTDATA 2 "6(-%'CURVESG6$7eo7$$\"\"!F)F(7$$\"0]i:gT<-\"!#;$\"0'fg_I#[u%F-7$$\"0+ DJ?$[V?F-$\"0x?Ws\\uH*F-7$$\"0](o/[AlIF-$\"0'*3nkWkO\"!#:7$$\"0+]iSmp3 %F-$\"0!3sig@&y\"F:7$$\"0+v$4'\\/8'F-$\"0\">wPySrDF:7$$\"0++D\"G$R<)F- $\"0&=/jp)HH$F:7$$F,F:$\"0LHwf8V&RF:7$$\"0+v=#**3E7F:$\"0TxRTh%fXF:7$$ \"0](=U#Q/V\"F:$\"0,3(GjB7^F:7$$\"0++Dc'yM;F:$\"0g$QDS>;cF:7$$F2F:$\"0 4jyhE2\\'F:7$$\"0+]P%)z@X#F:$\"0ts\"oo42sF:7$$\"0+vV[w3'GF:$\"0kqB<$G' y(F:7$$\"0++]7t&pKF:$\"0e$GcNfY#)F:7$$\"0]P%[Ey!)RF:$\"0c#y!eg0\"))F:7 $$\"0+v=<#*>p%F:$\"0]d*G(>#*F:7$$\"0 ]7`p,KS&F:$\"0&H\"HVhlE*F:7$$\"0DJqX1)edF:$\"0p;J$fqx#*F:7$$\"0+](=7T9 hF:$\"0m#fn`rd#*F:7$$\"0+DJ?qSr(F:$\"0ztV_46)))F:7$$\"0++v=HPJ*F:$\"0a OZ$ol;#)F:7$$\"0++]<*R#4\"!#9$\"0)3u1&\\-U(F:7$$\"0+]7VDMD\"Fer$\"0;*Q yeD\"f'F:7$$\"0]7yN'o89Fer$\"0bpcQ()yy&F:7$$\"0+vVGZRd\"Fer$\"0;I`IKb. &F:7$$\"0+D\"y'HDs\"Fer$\"0Y$[1[R'R%F:7$$\"0+v=276(=Fer$\"0!QOedA=QF:7 $$\"0+vV`g\\-#Fer$\"0dI^8)=%G$F:7$$\"0+vo**3)y@Fer$\"0R&y$\\VM\"GF:7$$ \"0+voH>zL#Fer$\"0QM(38t)Q#F:7$$\"0+vofHq\\#Fer$\"09<#QPl@?F:7$$\"0+v$ f'HU\"GFer$\"028R?3\"Q9F:7$$\"0++D\"*309$Fer$\"07Uza!)R+\"F:7$$\"0+]i& e*yU$Fer$\"0m#Gm$>1F(F-7$$\"0++v[D9v$Fer$\"0(ecu+wD]F-7$$\"0++Dc$GwSFe r$\"05'G%\\G!\\MF-7$$\"0++D^W$*Q%Fer$\"0G@]S(o(Q#F-7$$\"0+vo%QjtYFer$ \"0LH+(f=.0)H&Fer$\" 0K7'Qq1C!)!#<7$$\"0+v=-p6j&Fer$\"0i*G3J1\\`Fbx7$$\"0++]2Mg#fFer$\"0XzT M:Lt$Fbx7$$\"0+v=xZ&\\iFer$\"0z?9bR)=DFbx7$$\"0+DcJ4wb'Fer$\"0wtaErot \"Fbx7$$\"0+](=#R!zoFer$\"0B.np5[=\"Fbx7$$\"0+v$4A@urFer$\"0?[\\\\\\[R )!#=7$$\"0+]i:'f#\\(Fer$\"02fSDg8%eFaz7$$\"0+vof2L#yFer$\"0QH[FGM0%Faz 7$$\"0+D\"yG>6\")Fer$\"0X'\\\\EXyHFaz7$$\"0+](oo6A%)Fer$\"0f5'yQdsY7$$\"0+]i0A#*p*Fer$\"0cKRcreC'Fe\\l7$$\"0 ++vgc-+\"!#8$\"0A4bD:hq%Fe\\l7$$\"0+]i0XE.\"F^]l$\"0xq8ON?[$Fe\\l7$$\" 0](o/Q*>1\"F^]l$\"0H#35['*[EFe\\l7$$\"0+](Q(zS4\"F^]l$\"0;L+cS5'>Fe\\l 7$$\"0](=-,FC6F^]l$\"0\\>')Q/SZ\"Fe\\l7$$\"0]P4tFe:\"F^]l$\"0?!f*Q**** 3\"Fe\\l7$$\"0+]73\"o'=\"F^]l$\"0sD.)3c$3)!#?7$$\"0](oz;)*=7F^]l$\"0n= r(R*\\)eF__l7$$\"0++]*44]7F^]l$\"0xR1#35:VF__l7$$\"0+DJw/>G\"F^]l$\"0x _AT$*o7$F__l7$$\"0](=(4bMJ\"F^]l$\"0e.LY05E#F__l7$$\"0++vdYCM\"F^]l$\" 0/>(*=%Qr;F__l7$$\"0+Dc3ucP\"F^]l$\"0um@`ej<\"F__l7$$\"0++];$R09F^]l$ \"0&RJO?wa&)!#@7$$\"0](=-*zqV\"F^]l$\"0Y43ZdV1'Fcal7$$\"0]7G:3uY\"F^]l $\"0$R!z'4eWVFcal7$$\"#:F)$\"0gz>0iD-$Fcal-%'COLOURG6&%$RGBG$\"*++++\" !\")F(F(-F$6$7doF'7$F<$!0skW++-$\\F-7$FF$!//#ewe+$*)F:7$FO$!0d4a`?'>7F :7$FY$!0H&yG^%p[\"F:7$F\\o$!0zcX(\\6\")=F:7$Ffo$!09x\\]uo7#F:7$$\"0v=n )y#F:7$F[p$!0qt(z'zDC#F:7$$\"0Dc,T(QOVF:$!0g4ig(eqAF: 7$F`p$!0]1y.W0G#F:7$Fep$!0#*>IILPF#F:7$Fjp$!0:f/#4S^AF:7$F_q$!0_`B))pZ @#F:7$Fdq$!0\">%>!o/l@F:7$Fiq$!0GwAO>2\"=F:7$F^r$!0'oxEp898F:7$Fcr$!0h (>,Im'[(F-7$Fir$!0rAD3!4u;F-7$F^s$\"0*Gig*4f,%F-7$Fcs$\"0')pt%H7=&*F-7 $Fhs$\"0YnzI#=V9F:7$F]t$\"0]Q**fgb\">F:7$Fbt$\"0nqR.+AQ#F:7$Fgt$\"0Knb \"\\T@GF:7$F\\u$\"0>;)z:YSKF:7$Fau$\"0;5!oh08A%F:7$F[v$\"0%4u8,nHYF:7$$\"0+vVQ-UG$Fer$\"0 #p'[Wkeu%F:7$F`v$\"0I`%)[(GG[F:7$$\"0+v=ng'*e$Fer$\"0-C#F:7$Fhz$\"0+>F:7$F][l$\"0)4 sW*pUi\"F:7$Fb[l$\"0Ia^nENN\"F:7$Fg[l$\"0]]c'fD16F:7$F\\\\l$\"0C^'QAAp *)F-7$Fa\\l$\"0F_m(HPTsF-7$Fg\\l$\"0q4acZ?k&F-7$F\\]l$\"0l2c*zGjWF-7$F b]l$\"0Zd-@?3W$F-7$Fg]l$\"0Dsx]#\\&p#F-7$F\\^l$\"0]LN[Bk/#F-7$Fa^l$\"0 #[4(ffpc\"F-7$Ff^l$\"0#*RKkdk<\"F-7$F[_l$\"0$H!f:'fE))Fbx7$Fa_l$\"0Lm( )f&R(['Fbx7$Ff_l$\"0OYVT#G\"z%Fbx7$F[`l$\"0=7v/tC\\$Fbx7$F``l$\"0/BWV3 t`#Fbx7$Fe`l$\"0Qo;q4C)=Fbx7$Fj`l$\"0Wp(>e_H8Fbx7$F_al$\"0#)QRd\"*[p*F az7$Feal$\"0*G+d?S!*oFaz7$Fjal$\"0')Q*G@kZ\\Faz7$F_bl$\"0&RN7**4]MFaz- Fdbl6&FfblF(F(Fgbl-F$6$7apF'7$F<$\"/<)>e:\"3ZF:7$FF$\"/Dmb#R+N*F:7$FO$ \"0L>:@uLN\"F:7$FY$\"0U;Zzdfq\"F:7$F\\o$\"0`8'Go.%>#F:7$Ffo$\"0MiGrUMU #F:7$F`dl$\"0,%zlV/gCF:7$F[p$\"0!y#HU]sY#F:7$Fhdl$\"0;2Ifu2X#F:7$F`p$ \"0`Mf$>l:CF:7$Fjp$\"0x?@)[q0BF:7$Fdq$\"0.&*)R[6i@F:7$Fiq$\"0&e0,tLt%fRF:7$$\"0]PfGpt(HFer$!0@5`A4(GRF:7$F[v $!00w%=DWmQF:7$F[il$!0cM0&ph$y$F:7$F`v$!0/8lCH!pOF:7$Fev$!0?Et\"pzdKF: 7$F[jl$!0x$QI\"*o`HF:7$Fjv$!0;SWLevd#F:7$Fcjl$!0#Q2!)HdZ@F:7$F_w$!01k! y!Hjl\"F:7$$\"0](oz\"*[JXFer$!0w:0w?X;\"F:7$Fdw$!/&=.jWKQ'F:7$$\"0](o/ wlU[Fer$\".8g*y#[w\"F:7$Fiw$\"/&43K@O)oF:7$$\"0+]ilm[:&Fer$\"0E.')=vND \"F:7$F^x$\"0%R)=)\\$R!=F:7$$\"0]Pf[5YY&Fer$\"0!3&)\\\"z:T#F:7$Fdx$\"0 /8$>piqHF:7$$\"0]P%[:gydFer$\"0tMS'Rh:MF:7$Fix$\"0KKRsau!QF:7$F^y$\"0S $\\2#eOY%F:7$$\"0+]PayNS'Fer$\"0&4=pY&[n%F:7$Fcy$\"07KjbZ?#[F:7$$\"0D1 9zmzj'Fer$\"01L39aV([F:7$$\"0](=nUK=nFer$\"0:=$)zj0\"\\F:7$$\"0voHu\"o )z'Fer$\"0-Ax%>IJ\\F:7$Fhy$\"0+Rn))es$\\F:7$$\"0D1kYKG&pFer$\"0]x^qY.$ \\F:7$$\"0]iSrDm-(Fer$\"0ZnT\\0A\"\\F:7$$\"0v=<'*=/5(Fer$\"0+xFJoM)[F: 7$F]z$\"0'\\()GiyW[F:7$$\"0]7G=/ML(Fer$\"0pMbe41t%F:7$Fcz$\"0wBNcm*zXF :7$Fhz$\"0%G:R@-zTF:7$F][l$\"0B.%f*3\"oPF:7$Fb[l$\"0.UIzrqH$F:7$Fg[l$ \"0w\\0QbB\"GF:7$F\\\\l$\"0Aa1W(4jBF:7$Fa\\l$\"0!>i,uzl>F:7$Fg\\l$\"0K qKfJkd\"F:7$F\\]l$\"0%Rwx*RcF\"F:7$Fb]l$\"0T(pDni/5F:7$Fg]l$\"0WimRcm+ )F-7$F\\^l$\"0vdqyA=='F-7$Fa^l$\"0d(Rm%=7![F-7$Ff^l$\"0\")G*GhP`OF-7$F [_l$\"0&ehZ;vtFF-7$Fa_l$\"0()GN/-=1#F-7$Ff_l$\"0/,SE3z`\"F-7$F[`l$\"0G .!QfXJ6F-7$F``l$\"0S!>gOu*G)Fbx7$Fe`l$\"0'y8wZG%>'Fbx7$Fj`l$\"0%[>rUE3 WFbx7$F_al$\"0B6Pr`ZB$Fbx7$Feal$\"0kL*[-Q8BFbx7$Fjal$\"0DTD_k.n\"Fbx7$ F_bl$\"0k^$*oy7<\"Fbx-Fdbl6&FfblF(FgblF(-F$6$7iqF'7$F<$!/S\"QX$QB))F:7 $FF$!0A?Q*4?([\"F:7$FO$!0)G-qSND>F:7$FY$!0Rp+!zPlAF:7$F\\o$!0n=34!)*)y #F:7$Ffo$!0EPR)QYzJF:7$F[p$!0js?oKjS$F:7$F`p$!0,h*o#3Q]$F:7$Fep$!06k2] 0&)\\$F:7$Fjp$!0\"[Z\\n.cMF:7$F_q$!0E4Q7'zwLF:7$Fdq$!0'Q]J&4@E$F:7$$\" 0]P4rSU\"pF:$!0,b(yLx')GF:7$Fiq$!0E1mv ;y\"F:7$F^r$!0Y:vp+%Q6F:7$$\"0+vo/')=,\"Fer$!/$R9/45$[F:7$Fcr$\"/2'\\' )zo\\\"F:7$$\"0+DJI7H<\"Fer$\"/'[(3%yUR(F:7$Fir$\"0%)QRt*)\\F\"F:7$F^s $\"0DDI:2H<#F:7$Fcs$\"0j=.+&H\")GF:7$Fhs$\"0zP1CZIU$F:7$F]t$\"0?#z5Z#R *QF:7$Fbt$\"0#3O0LM>VF:7$Fgt$\"0Jb()HrTm%F:7$$\"0+vo9k$eAFer$\"0EJk%3c )z%F:7$F\\u$\"0xlmsKY*[F:7$$\"0+v=(opxBFer$\"0dw[Mei#\\F:7$$\"0+voWuuT #Fer$\"0qK'HA!f%\\F:7$$\"0+v=-_sX#Fer$\"0?afdPG&\\F:7$Fau$\"03(RdgRY\\ F:7$$\"0](o/'zm`#Fer$\"0%y?xq1E\\F:7$F^fm$\"0s6!piL\"*[F:7$$\"0]7.izfh #Fer$\"0X3Ykt<%[F:7$F`hl$\"0&)[Y!4/xZF:7$Fffm$\"03\"Gxv>,YF:7$Ffu$\"0^ D92!3jVF:7$F^gm$\"0j@B+8Io$F:7$F[v$\"0_&z9OOtFF:7$$\"0]P%[cN7KFer$\"0b Eh*Hj7BF:7$F[il$\"0J]mdvF#=F:7$$\"0]7.7\\gN$Fer$\"0)*>VR;)48F:7$F`v$\" 0$4Obj)4!yF-7$$\"0]iSEy(3NFer$\"0&GobyO=QxY%F:7$F_w$!0G9wJnn,&F:7$F[im$!0JCh%QQJ`F:7$Fdw$!0%) 37Be9Z&F:7$$\"07Gjy*)er%Fer$!0=08$3Y![&F:7$$\"0D\"yDd9eZFer$!00uDJF]Z& F:7$$\"0PM_m,/![Fer$!0IQg8?bX&F:7$Fcim$!0+!G')4MAaF:7$$\"0v$f$[pr#\\Fe r$!0H+?s[nJ&F:7$Fiw$!0%fT!R;@;&F:7$F[jm$!0Vo;If.![F:7$F^x$!0)eU'phFL%F :7$Fcjm$!02Nv@$3(o$F:7$Fdx$!0PUKB*HnHF:7$F[[n$!0!z>qf_'H#F:7$Fix$!0-%Q +%Gsh\"F:7$$\"0]PM#4z(3'Fer$!/E!f6x*\\))F:7$F^y$!/hWZM^q=F:7$Ff[n$\"/x dSq$GI%F:7$Fcy$\"0\"*)3%Q\"R0**F-7$Fc\\n$\"0FIB/3k]\"F:7$Fhy$\"0TgF:7$F]z$\"0l*3tR\"*fDF:7$Fb^n$\"0>\"*GYEvy#F:7$Fcz$\"0%zHm2d[HF:7$$ \"0D1k,u_d(Fer$\"0#>DGNz2IF:7$$\"0]il(=&zl(Fer$\"0s1t!\\Y^IF:7$$\"0v=n tH1u(Fer$\"0pB%H*403$F:7$Fhz$\"0eFP.\")e4$F:7$$\"0](=<*y_*yFer$\"0!4%Q C$)))4$F:7$$\"0++vB]s'zFer$\"0v9Mz**G4$F:7$$\"0]7yb@#R!)Fer$\"0&HCY:dy IF:7$F][l$\"0IXPgLl0$F:7$$\"0]PM([lm#)Fer$\"0%\\%4&4c&)HF:7$Fb[l$\"0k8 OD%f()GF:7$Fg[l$\"0I$Q7z^BEF:7$F\\\\l$\"0?:T9g#=BF:7$Fa\\l$\"0S'e8bk3? F:7$Fg\\l$\"0$4UnHXt;F:7$F\\]l$\"0k*3V![WR\"F:7$Fb]l$\"02&oD?XG6F:7$Fg ]l$\"0cYCg^+>*F-7$F\\^l$\"0p@gqubC(F-7$Fa^l$\"0XO))Q9ps&F-7$Ff^l$\"0\\ 2#[)e*HWF-7$F[_l$\"0fA$\\bP7MF-7$Fa_l$\"0$yc8W\\rDF-7$Ff_l$\"0V8T%=@T> F-7$F[`l$\"0]orDNUW\"F-7$F``l$\"0/c^6S*o5F-7$Fe`l$\"0vp-&=Wc!)Fbx7$Fj` l$\"0Kek)e*fy&Fbx7$F_al$\"0G'R2*pyF%Fbx7$Feal$\"0())[Y(RB3$Fbx7$Fjal$ \"0gm#e,ZSAFbx7$F_bl$\"0D(R3Gd\"e\"Fbx-Fdbl6&FfblF)F)F)-%+AXESLABELSG6 $Q\"r6\"Q!6\"-%%VIEWG6$;F(F_bl%(DEFAULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" "Curve 3" "Curve \+ 4" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 53 "The basis functions do not \+ have a good span at large " }{TEXT 270 1 "r" }{TEXT -1 88 " where the \+ higher states wish to reside. A smaller value of beta should improve m atters." }}{PARA 0 "" 0 "" {TEXT -1 133 "For this reason one mixes usu ally STO's with different values of beta. The Gram-Schmidt process tak es care of the orthonormalization." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT 271 11 "Exercise 3:" }}{PARA 0 "" 0 "" {TEXT -1 300 "Calculate the eigenfunctions for a few low-lying eigenstates in t he L=2 and L=3 symmetry sectors. Observe how the eigenenergies approac h the hydrogenic result in this case of E_n=-1/(2n^2). Can you explain this behaviour? Hint: graph the effective potential, and compare with that of a hydrogen atom." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "0 0 0" 12 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }