Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Solving using schrodinger equation techniques

  1. Feb 20, 2009 #1
    Hello,
    I have problem I wish to solve, and I wonder if anyone already delt with it when solving the schrodinger 2D equation.

    say E(x,y) is a scalar field function that complies with

    ( [tex]\frac{d}{dx}[/tex]2+[tex]\frac{d}{dy}[/tex]2 ) *E(x,y)+k(x,y)*E(x,y)=k1*E(x,y)


    where k(x,y)={k2 for x2+y2<R2 and 0 otherwise}, i.e. a tube potential.
    All is known but E(x,y).
    I think it can be examined as a 2D schrodinger equation, even thow the eigenvalue k1 is known.

    How can I get to start finding the solutions of this equation?
    Can I expect to know how many are there? - one, two, many?
     
  2. jcsd
  3. Feb 20, 2009 #2

    StatusX

    User Avatar
    Homework Helper

    At least in 3D, a rotationally invariant system is best dealt with in spherical coordinates, expanding in eigenfunctions of angular momentum and then solving the radial equation, which will be an ODE. In 2D, the "spherical harmonics" are just the functions [itex]e^{im \phi}[/itex] where m is an integer. So try writing [itex]f_m(r,\phi) = u_m(r) e^{i m \phi} [/itex], which should give you an ordinary differential equation for [itex]u_m(r)[/itex]. Then the general solution will be a linear combination of the [itex]f_m[/itex], ie, you should expect a linearly independent solution for each integer m.

    This is a little messier in 2D than in 3D, and the equation will have a first derivative term that doesn't appear in the ordinary schrodinger equation. In fact, the solutions [itex]u_m(r,\phi)[/itex] are probably going to be Bessel functions.
     
  4. Feb 20, 2009 #3
    Numerically you could use pdetool within matlab (type: pdetool at the matlab prompt), where you could solve this eigen value problem in 2D (Its not that you know or could specify k1, you will get it from the numeric solution). Other better program is COMSOL multiphysics.

    This problem looks like cylinder symmetry, but depend on weather k(x,y)=f(r) or not. In that case you could obtain an effective 1D eigen value problem in Psi(r), using the angular quantum quantum number m as specified in the other post here.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Solving using schrodinger equation techniques
Loading...