Solve Schroedinger Equation with Mathematica DSolve for given potential

donquixote17
Messages
9
Reaction score
0

Homework Statement


I need to solve the Schroedinger equation (Using DSolve in Mathematica) for a potential that is infinite below z=0 and V=mgz for positive z.


Homework Equations


TISE:
\psi \text{''}[z]+\frac{2 m}{\hbar }(\text{En}-V[z])\psi [z]==0



The Attempt at a Solution


First I defined the potential as
V[z]=\text{Piecewise}[\{\{\infty ,z<0\},\{m g z,z\geq 0\}\}];

Then I told Mathematica to solve the Diffeq
\text{DSolve}\left[\left\{\psi \text{''}[z]+\frac{2 m}{\hbar }(\text{En}-V[z])\psi [z]==0,\psi [0]==0,\psi [\infty ]==0\right\},\psi [z],z\right]

The error mathematica gave me was
InverseFunction::ifun: Inverse functions are being used. Values may be lost for multivalued inverses. >>

I'm not really sure how to get mathematica to give me the solution. I know I need 2 conditions besides the differential equation in order to solve. The only boundary condition that I know to use is that the wavefunction must be zero at z=0 since the potential is infinite there.

Any ideas on what second boundary condition I should use or any ideas on where I'm going wrong? Thanks
 
Last edited:
Physics news on Phys.org
It might be choking, in part, on the infinite potential. You know the wavefunction vanishes there, so try using V(z)=mgz and solve for just z>0.
 
Thread 'Need help understanding this figure on energy levels'
This figure is from "Introduction to Quantum Mechanics" by Griffiths (3rd edition). It is available to download. It is from page 142. I am hoping the usual people on this site will give me a hand understanding what is going on in the figure. After the equation (4.50) it says "It is customary to introduce the principal quantum number, ##n##, which simply orders the allowed energies, starting with 1 for the ground state. (see the figure)" I still don't understand the figure :( Here is...
Thread 'Understanding how to "tack on" the time wiggle factor'
The last problem I posted on QM made it into advanced homework help, that is why I am putting it here. I am sorry for any hassle imposed on the moderators by myself. Part (a) is quite easy. We get $$\sigma_1 = 2\lambda, \mathbf{v}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} \sigma_2 = \lambda, \mathbf{v}_2 = \begin{pmatrix} 1/\sqrt{2} \\ 1/\sqrt{2} \\ 0 \end{pmatrix} \sigma_3 = -\lambda, \mathbf{v}_3 = \begin{pmatrix} 1/\sqrt{2} \\ -1/\sqrt{2} \\ 0 \end{pmatrix} $$ There are two ways...
Back
Top