Example IV. More accurate method for calculating bound states

<BODY LANG="EN"> <HR><center> <table border=1> <td><A NAME="tex2html42" HREF="node3.html"><img height=36 width=36 alt="PREVIOUS" src="/UCESIcons/tools/previous.gif"></A></td> <td><A NAME="tex2html48" HREF="Lab2.html"><img height=36 width=36 alt="UP" src="/UCESIcons/tools/up.gif"></A></td> <td><A NAME="tex2html50" HREF="node5.html"><img height=36 width=36 alt="NEXT" src="/UCESIcons/tools/next.gif"></A></td> <td><a href="contents-node4.html"><img height=36 width=36 alt="CONTENTS" src="/UCESIcons/tools/toc.gif"></a></td> <td><a href="/hamlet-bin/qhandler?exit" target="results"><img height=36 width=36 alt="RESET" src="/UCESIcons/tools/stop.gif"></a></td> <td><a href="/UCES/hamlet/"><img height=36 width=36 alt="HAMLET" src="/UCESIcons/tools/hamlet.gif"></a></td> <td><a href="copyright.html"><img height=36 width=36 alt="COPYRIGHT" src="/UCESIcons/tools/copyright.gif"></a></td> <td><a href="mailto:uces_info@krellinst.org"><img height=36 width=36 alt="FEEDBACK" src="/UCESIcons/tools/feedback.gif"></a></td> <td><a href="/UCES/hamlet/usage.html"><img height=36 width=36 alt="ABOUT" src="/UCESIcons/tools/about.gif"></a></td> </table> </center> <HR><H1><A NAME="SECTION00040000000000000000">Example IV. More accurate method for calculating bound states</A></H1> <P> To avoid the numerical problems in integrating <I>into</I> a classically forbidden region, which exist on each side of the potential well, the calculation can be broken up into two parts. First you generate a solution to the Schrödinger equation starting deep (i.e. `deep enough') inside the classically forbidden region on the left side of the well and go out to a point that is close to the middle of the well (the `joining point'). Then start on the other side of the well, deep inside the classically forbidden region and go in the opposite direction into the well up to the same point. You can make your two functions match at the joining point by renormalizing them, but the slope will not be continuous unless you happen to have chosen a correct bound state energy value in the numerical calculation. <P> A file called `boundstates.ma' located in the chem465 directory contains a Mathematica Notebook for doing this calculation. Make a copy of it in your directory by entering (in the UNIX window) (note the last "space dot") <P> <code>~chem465/boundstates.ma .</code> <P> The parameters chosen there for the Morse potential and the reduced mass are those for a <I>H</I>-<I>Cl</I> molecule. The starting guess for the energy is close to the ground state energy. <P> <I>Change the guess for the energy and observe how the slopes change at the joining point. Then implement the correction formula derived in the lecture notes to add a numerical estimate of how to modify the energy to approximate better the bound state energy. You will need to estimate the discontinuity of the slopes at the joining points (using finite differences) and you need to get the sum of the wavefunction squared.</I> <P> <HR><center> <table border=1> <td><A NAME="tex2html42" HREF="node3.html"><img height=36 width=36 alt="PREVIOUS" src="/UCESIcons/tools/previous.gif"></A></td> <td><A NAME="tex2html48" HREF="Lab2.html"><img height=36 width=36 alt="UP" src="/UCESIcons/tools/up.gif"></A></td> <td><A NAME="tex2html50" HREF="node5.html"><img height=36 width=36 alt="NEXT" src="/UCESIcons/tools/next.gif"></A></td> <td><a href="contents-node4.html"><img height=36 width=36 alt="CONTENTS" src="/UCESIcons/tools/toc.gif"></a></td> <td><a href="/hamlet-bin/qhandler?exit" target="results"><img height=36 width=36 alt="RESET" src="/UCESIcons/tools/stop.gif"></a></td> <td><a href="/UCES/hamlet/"><img height=36 width=36 alt="HAMLET" src="/UCESIcons/tools/hamlet.gif"></a></td> <td><a href="copyright.html"><img height=36 width=36 alt="COPYRIGHT" src="/UCESIcons/tools/copyright.gif"></a></td> <td><a href="mailto:uces_info@krellinst.org"><img height=36 width=36 alt="FEEDBACK" src="/UCESIcons/tools/feedback.gif"></a></td> <td><a href="/UCES/hamlet/usage.html"><img height=36 width=36 alt="ABOUT" src="/UCESIcons/tools/about.gif"></a></td> </table> </center> <HR><P><ADDRESS> <I>Hannes Jonsson<BR> Modified by Thomas L. Marchioro II<BR> and the <A NAME="tex2html1" HREF="/UCES/">Undergraduate Computational Engineering and Science project</A> </I> </ADDRESS> </BODY>