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!

Hydrogen atom:<K>,<V>,momentum distribution

  1. Mar 18, 2008 #1
    Hello to everybody! This is my first time in PF. I have problems with an QM exercise.

    1. The problem statement, all variables and given/known data
    About hydrogen atom in the ground state (n=1), evaluate:
    <K>, <V> (expectation value of the kinetic energy and potential energy) and the momentum distribution ([tex]\varphi[/tex](p)).



    3. The attempt at a solution

    for now I treat only <K>:
    see attachment 009 (<K>=<T>=kinetic energy)
    but i know the exact solution so that I evaluate is not right.
    where is my mistake??

    please help me.

    thanks
     

    Attached Files:

  2. jcsd
  3. Mar 19, 2008 #2
    okay. maybe my attachment is not readble.

    the wave function for the hydrogen atom in 1s state is:= [tex]\varphi(r)=\frac{exp(-r/a)}{\sqrt{\pi a^3}}[/tex]

    where:

    a is the Bohr radius := [tex]a=\frac{\hbar^2}{m e^2}[/tex] ;
    m is the electron mass ;
    [tex]e=\frac{q^2}{4\pi\epsilon}[/tex] ;
    q is the elementary charge ;

    now:

    to find <K> (expectation value of kinetic energy) I have to evaluate:

    [tex]\int \frac{p^2}{2m}\left|\varphi(r)\right|^2 r^2 sin\vartheta d\vartheta d\phi dr[/tex] from 0 to infinity

    i.e. = [tex]\frac{(-i\hbar)^2}{2m\pi a^3} \int \partial^{2}/\partial r^{2}[exp(-2r/a)] r^2 sin\vartheta dr d\vartheta d\phi = \frac{(-i\hbar)^2}{2m\pi a^3} 4 \pi \int \frac{4r^2}{a^2} exp(-2r/a) dr = \frac{(-i\hbar)^2}{2m\pi a^3} 4 \pi a = \frac{(-2i\hbar)^2}{m a^2} = \frac{-2 \hbar^2 m^2 e^4}{m \hbar^4} = \frac{-2me^4}{\hbar^2} [/tex]

    but the exact solutions found on Basdevant or Sakurai is:= [tex]K=\frac{+me^4}{2 \hbar^2} [/tex]

    so...where is my mistake??
    please, help me to find it!!
     
  4. Mar 19, 2008 #3

    clem

    User Avatar
    Science Advisor

    [tex]\nabla^2[/tex] is not [tex]\partial^2/\partial r^2[/tex] in spherical coordinates.
     
  5. Mar 19, 2008 #4
    ops....it's true...

    I tried with:

    [tex]<K>=\frac{1}{\pi a^3}\frac{-\hbar^2}{2m} \int \frac{1}{r^2} \partial/\partial r (r^2 \partial/\partial r) exp(\frac{-2r}{a}) r^2 sin\vartheta d \theta d \phi d r [/tex] from 0 to infinity

    .............................
    .....................................

    then

    [tex]\frac{-2 \hbar^2}{m a^3} \frac{-2}{a} \int exp(\frac{-2r}{a}) (2r-\frac{2r}{a}) dr =

    = \frac{-2 \hbar^2}{m a^3} \frac{-2}{a} [r^2 exp(\frac{-2r}{a})]^{inf}_{0} [/tex]

    and here I can't have a finite solution....
    is the integral correct?
    how to have a finite solution for [tex] [r^2 exp(\frac{-2r}{a})]^{inf}_{0} [/tex] ????

    thanks a lot
     
  6. Mar 19, 2008 #5

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    [tex][r^2 exp(\frac{-2r}{a})]^{\infty}_{0} = 0[/tex]

    and I dont think your interation is correctly.

    you should end up with an integral like:
    [tex] \int ^{\infty}_0(2r-2r^2/a)e^{-2r/a}dr [/tex]

    And use the standard integral:
    [tex] \int ^{\infty}_0r^ne^{-\alpha r}dr = (1/\alpha ^{n+1}) \Gamma (n+1) [/tex]


    Useful properties of the gamma function:
    [tex] \Gamma (m+1) = m\Gamma (m) [/tex]
    [tex]\Gamma (1) = 1 [/tex]
    [tex]\Gamma (1/2) = \sqrt{\pi} [/tex]
     
    Last edited: Mar 19, 2008
  7. Mar 19, 2008 #6
    I check the solution of the integral with the software "Derive" and it's correct!!!

    ....?????????????????
     
  8. Mar 19, 2008 #7

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    [tex]\int {r^2} \partial/\partial r (r^2 \partial/\partial r) exp(\frac{-2r}{a}) r^2 d r = (-2/a)\int (2r-2r^2/a)e^{-2r/a}dr[/tex]

    ??
     
    Last edited: Mar 19, 2008
  9. Mar 19, 2008 #8
    ...and I tried also to solve with Gamma function but the integral = 0 !!!!!
     
  10. Mar 19, 2008 #9
    yes...it's the same of mine!
     
  11. Mar 19, 2008 #10

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    ee no..

     
  12. Mar 19, 2008 #11

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    Strange.
     
  13. Mar 19, 2008 #12
    oh,iI'm sorry but it's only a writing-mistake....in my works I've what you've posted too.....
     
  14. Mar 19, 2008 #13

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    ok, but I think I know what the error is, you have done things in wrong order:

    [tex]\int \psi (r) ^*\Delta\psi (r)r^2dr [/tex]

    not

    [tex]\int \Delta|\psi (r)|^2r^2dr [/tex]

    if you know what I mean. This has to do with the definition of expectation value and general properties of derivative operators.

    \Delta is the Laplacian Operator, which acts on everything to the RIGHT.
     
    Last edited: Mar 19, 2008
  15. Mar 19, 2008 #14
    malawi_glen......you're great....i try to correct it right now!!!!
    thanks a lot.....
     
  16. Mar 19, 2008 #15

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    You shold not get 0 now, good luck!
     
  17. Mar 19, 2008 #16
    yes!!! the result is correct!! thanks a lot!!
     
  18. Mar 20, 2008 #17
    momentun distribution of the hydrogen atom in ground state

    ok, now I'm trying to evaluate the momentun distribution of the hydrogen atom in ground state.

    I should apply the Fourier Transform to the [tex]\psi (r)[/tex]

    i.e. [tex]\varphi ( \textbf{p} ) = \frac{1}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int^{-\infty}_{+\infty} exp(-\textbf{r} /a) exp(i \textbf{pr} ) d^{3}r [/tex]

    mmm..I'm not good in these -hard math works- ....

    my problem is also about the -morevariables- integral...

    I don't really know how to start....

    Since the wave function depends only from r, I've just tried to do:

    [tex]\varphi ( p ) = \frac{1}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int^{-\infty}_{+\infty} exp(-r /a) exp(ipr} ) dr =
    \frac{1}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int^{-\infty}_{+\infty} exp(\frac{-r}{a} + ipr) ) dr [/tex]

    ...but the result of the integral is something very strange...I think it's not the good way.....

    can you help me please???
     
  19. Mar 20, 2008 #18

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    dude [tex] d\mathbf{r} = r^2sin \theta d\theta drd\phi[/tex] or this is even better in this case:
    [tex] d\mathbf{r} = r^2 d(\cos \theta )drd\phi[/tex]

    so that:

    [tex]\varphi ( \textbf{p} ) = \frac{1}{(2 \pi)^{3/2}} \int \varphi (\textbf{r}) \exp(i \textbf{p} \cdot \textbf{r}) d^{3}r = [/tex]

    [tex]\frac{1}{(2 \pi)^{3/2}} \int \varphi (\textbf{r}) \exp (i \textbf{p} \cdot \textbf{r}) r^2d(\cos \theta ) drd\phi [/tex]

    And you integrate r from 0 to infinity, theta from 0 to pi (but you now have cos(theta) as integration variable, so be careful, and phi from 0 to 2*pi.
    Use the forumla for dot-product:
    [tex] \textbf{a} \cdot \textbf{b} = ab \cos \theta [/tex]

    But you must use the SPATIAL wave function, just not the radial part. i.e
    [tex] \varphi (\textbf{r}) = \varphi (r) \cdot Y_0^0 (\theta ,\phi )= \frac{\exp(-r/a)}{\sqrt{\pi a^3}}\cdot 1/(\sqrt{4 \pi}) [/tex]

    I hope you get the picture, I might miss some constants etc, since I have not done this QM for 1 year, but the procedures I remember ;-) I hope
     
    Last edited: Mar 20, 2008
  20. Mar 20, 2008 #19

    malawi_glenn

    User Avatar
    Science Advisor
    Homework Helper

    I did som crucial editing, my latex was wrong. Now it shold be more correct.
     
  21. Mar 20, 2008 #20
    my questions:

    1)
    if [tex] d cos (\theta) / d\theta = -sin(\theta)[/tex] ,
    so [tex] sin(\theta) d \theta = -cos(\theta) [/tex]....isn't true???

    2)
    I check for the spatial wave function and the term [tex]\frac{1}{\sqrt{4\pi}}[/tex] is included in the [tex]\varphi(r)=\frac{exp(-r/a)}{\sqrt{\pi a^3}}[/tex]....,no??

    then, I tried to cumpute:
    [tex]
    \varphi ( \textbf{p} ) = \frac{1}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int e^{-r/a -ipr cos \theta } r^2 sin \theta dr d\theta d\phi

    = \frac{1}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int^{\infty}_{0}\int^{\pi}_{0} \int ^{2\pi}_{0} e^{-r/a -ipr cos \theta } r^2 sin \theta dr d\theta d\phi =
    [/tex]

    [tex]\frac{2\pi}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int^{\infty}_{0}\int^{\pi}_{0} e^{-r/a -ipr cos \theta } r^2 sin \theta dr d\theta =
    [/tex]

    [tex]
    \frac{4\pi}{(2 \pi)^{3/2}} \frac{1}{ \sqrt {\pi a^3}} \int^{\infty}_{0} \frac{e^{-r/a } r sin (pr)}{p} dr =
    [/tex]

    but now I can't solve this.....

    how can I do??

    is what I wrote correct?
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?



Similar Discussions: Hydrogen atom:<K>,<V>,momentum distribution
Loading...