1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Unidimensional potential and Schroedinger

  1. Nov 21, 2009 #1
    1. The problem statement, all variables and given/known data

    Consider a particle with charge $e$ which moves in $x$ axys (in particular in the positive region) with the potential $V(x)=-\alpha/x$ ($\alpha>0$). Find the ground state energy.


    2. Relevant equations

    I think that in some way I have to solve this ODE:

    $-\bar{h}^2/(2*m)*(d^2\psi)/(dx^2)+(\alpha/x+E)\psi=0


    3. The attempt at a solution
    I don't know how to solve this awful ODE. May you give me an hint?
     
  2. jcsd
  3. Nov 21, 2009 #2
    Welcome to Physics Forums Matt!

    [tex]
    -\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}=-\left(\alpha x+E\right)\psi
    [/tex]

    is an awful ODE? I think you might be confusing yourself. Is this an awful ODE:

    [tex]
    \frac{d^2y}{dx^2}=\kappa^2 y
    [/tex]

    EDIT: Actually, your ODE is wrong, it should be

    [tex]
    -\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}+\alpha x\psi=E\psi\Longrightarrow-\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2}+\left(\alpha x-E\right)\psi=0
    [/tex]
     
    Last edited: Nov 21, 2009
  4. Nov 21, 2009 #3
    I am sorry for my mistakes.
    The potential is in the form [tex]V(x)=\frac{-\alpha}{x}[/tex], where [tex]\alpha[/tex] is positive. The ODE should be:

    [tex]\frac{-\bar{h}^2}{2m}\frac{d^2\psi}{dx^2}-\frac{\alpha}{x}\psi=E\psi[/tex]

    I know that ODEs like [tex]y''(t)+ky(t)=0[/tex] are very simple, but this happens if [tex]k[/tex] is a real number.

    I think that I can solve this ODE in the same way as it is usually done for the hydrogen atom. Is this correct?
     
  5. Nov 21, 2009 #4
    Ah, that does certainly change the problem, but not by much. You are looking at a central potential (in 1D) and usually these potentials are solved as functions of [itex]V(x)[/itex]. So, from all the textbooks I've worked through on central potentials, you solve it as

    [tex]
    -\frac{\hbar^2}{2m}\,\frac{d^2\psi}{d x^2}+V(x)\psi(x)=E\psi(x)\longrightarrow\frac{d^2\psi}{dx^2}=\kappa^2\psi
    [/tex]

    where [tex]\kappa=\sqrt{2m\hbar^{-2}\left(V(x)-E\right)}[/tex].
     
  6. Nov 22, 2009 #5
    I am pretty sure that I have to find a set of solutions in function of [tex]x[/tex]. So the question still remains :biggrin:
     
  7. Nov 22, 2009 #6
    Well to find the energy eigenvalues (and thus the ground state energy), we'd need some boundary conditions. Where there any given in the problem? If there weren't any BC's, then you'll probably need the WKB approximation
     
  8. Nov 22, 2009 #7
    I don't think that this problem is so difficult (I don't even know what WKB approximation is :biggrin:).
    I finally managed to solve this ODE: I found the ground state autofunction

    [tex]\psi_1=\frac{2\sqrt{a}x}{a}\exp{-\frac{x}{a}}[/tex], where [tex]a=\frac{4\bar{h}^2}{me^2}[/tex] (I put [tex]\alpha=\frac{e^2}{4}[/tex], where [tex]e[/tex] is the charge of electron). Finding the fundamental energy and <x> is straightward.
     
  9. Nov 22, 2009 #8
    Actually, now that I've had a nights sleep on this problem, you probably could find (estimate) the ground state energy a lot more easily by using the Heisenberg uncertainty principle:

    Since your Hamiltonian takes the form

    [tex]
    H=\frac{p^2}{2m}-\frac{\alpha}{x}
    [/tex]

    then we can say that the ground state energy can be estimated via

    [tex]
    E_0=\langle\frac{p^2}{2m}\rangle-\langle\frac{\alpha}{x}\rangle
    [/tex]

    but through the uncertainty in [itex]p[/itex], we can say [itex]\langle p^2\rangle=(\Delta p)^2\geq\hbar^2/4a^2[/itex] where [itex]a[/itex] is the characteristic length of the atom, such that [itex]\langle \alpha/x\rangle=\alpha/a[/itex]. The ground state energy can then be written as

    [tex]
    E_0\geq\frac{\hbar^2}{8ma^2}-\frac{\alpha}{a}
    [/tex]

    to find [itex]a[/itex], we take the derivative of [itex]E_0[/itex] with respect to [itex]a[/itex] and set it equal to 0. You'll take this value of [itex]a[/itex], place it into the above equation and reduce it as much as possible and this would be your estimated ground state energy (usually close within a factor or two of [itex]\pi^2[/itex]). This type of ground state energy estimation can be done in this manner for any central potential.
     
  10. Nov 22, 2009 #9
    The way you solved the problem is pretty interesting.
    Thank you for the attention!
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook