Yukawa Potential

ManyNames

Pertubation expansion is done by so:

[latex]E_n=\sum^{\infty}_{j=0} \lambda^{j}E^{j}_{n}[/latex]

Its wave function is expanded in terms given by:

[latex]|\psi>=\sum^{\infty}_{j=0} \lambda^{j}|\psi_{n}>[/latex]

This all works when dealing with Eigenvalues dealing with [latex]H=H_{o}+ \lambda H[/latex] where [latex]\lambda[/latex]is the expansion parameter.

nickjer

You are still leaving out constants, like charge. Also, you completely changed the solution. Not sure what happened but now things are worse.

Go back and carefully step through everything on paper that way you can keep track of everything.

Hart

These are my exact full calculations:

[tex]V_{\gamma}(r) = -\left(\frac{q^{2}}{4\pi \epsilon_{0}r}\right)e^{-\gamma r}[/tex]

[tex]\left<\psi(r)|V_{\gamma}(r)|\psi(r)\right> = \left<\psi(r)|-\left(\frac{q^{2}}{4\pi \epsilon_{0}r}\right)e^{-\gamma r}|\psi(r)\right> = 4\pi \int^{\infty}_{0}\psi(r)V_{\gamma}(r)\psi(r)r^{2}dr[/tex]

[tex]=\left(-\frac{4\pi q^{2}}{4\pi \epsilon_{0}}\right)\int^{\infty}_{0}r e^{-\gamma r}\psi(r)\psi(r) dr[/tex]

[tex]=\left(-\frac{q^{2}}{\epsilon_{0}\pi a^{3}}\right)\int^{\infty}_{0}r e^{-\left(\frac{2}{a}+\gamma\right)r} dr[/tex]

[tex]=\left(-\frac{q^{2}}{\epsilon_{0}\pi a^{3}}\right)\left(-\frac{4 \pi}{\frac{2}{a}+\gamma}\right)[/tex]

[tex]=\frac{4 q^{2}}{\epsilon_{0}a^{3}\left(\frac{2}{a}+\gamma\right)}[/tex]

.. I'm hoping this is along the right lines, though I can't see that I've missed anything else now, but I'm unsure what I need to do further now.

nickjer

The solution to the integral in the 2nd to last line is incorrect. Scroll back a page and see how we solved it.

Hart

.. that's what I've used, obviously with the value of [itex]\alpha[/itex] substituted.

.. OH wait, it's [itex]\alpha^{2}[/itex] not [itex]\alpha[/itex], so then the result should be:

[tex]=\frac{4 q^{2}}{\epsilon_{0}a^{3}\left(\frac{2}{a}+\gamma\right)^{2}}[/tex]

correct now?

nickjer

Then where did the 4*pi come from since you used it in line 2 already. And why isn't the denominator squared?

nickjer

Also, once you finished this, you will need to do the same for the coulomb potential. Since you need to do

[tex]<\psi_0|V_\lambda-V|\psi_0>=<\psi_0|V_\lambda|\psi_0>-<\psi_0|V|\psi_0>[/tex]

But that will be much easier to solve, since we solved it basically already. You would replace [tex]\alpha[/tex] from your previous integral with just [tex]2/a_0[/tex]

Hart

so it should be:

[tex]\left<\psi(r)|V_{\gamma}(r)|\psi(r)\right> = \left<\psi(r)|-\left(\frac{q^{2}}{4\pi \epsilon_{0}r}\right)e^{-\gamma r}|\psi(r)\right> =\left(-\frac{q^{2}}{\epsilon_{0}\pi a^{3}}\right)\left(-\frac{1}{\left(\frac{2}{a}+\gamma}\right)^{2}\right)[/tex]

better?

nickjer

Forgot to mention, the solution is positive. You subtract the last limit of r=0, making it positive.

Hart

erm.. yep. so should be:

[tex]\int_{0}^{\infty} r e^{-\left(\frac{2}{a}+\gamma\right)r}dr = \left(\left(\frac{r}{\left(\frac{2}{a}+\gamma\right)}-\frac{1}{\left(\frac{2}{a}+\gamma\right)^{2}}\right)e^{-\left(\frac{2}{a}+\gamma\right)r}\right)\right|^{\infty}_{0} = \frac{4\pi}{\left(\frac{2}{a}+\gamma\right)^{2}}[/tex]

right?

but that doesn't help simplify this:

[tex]\left<\psi(r)|V_{\gamma}(r)|\psi(r)\right> = \left<\psi(r)|-\left(\frac{q^{2}}{4\pi \epsilon_{0}r}\right)e^{-\gamma r}|\psi(r)\right> =\left(-\frac{q^{2}}{\epsilon_{0}\pi a^{3}}\right)\left(\frac{1}{\left(\frac{2}{a}+\gamma}\right)^{2}\right)[/tex]

does it?

nickjer

There is no 4*Pi in the first line, since you already pulled out the angular integrals. The last line is now correct. To simplify it, I would recommend pulling out the 2/a from the denominator so it more resembles the ground state of the coulomb potential.

Now go back and solve this again for just the coulomb potential. You will notice a similarity.

Hart

um, don't get how I can rearrange that how you said

.. also, so to get the final answer I need also caclulate this:

[tex]\int_{0}^{\infty} r e^{-\left(\frac{2}{a_{0}}\right)r}dr = \left(\left(-\frac{a_{0}r}{2}-\frac{a_{0}^{2}}{4}\right)e^{-\left(\frac{a_{0}}{2}\right)r}\right)\right|^{\infty}_{0} = \left(\frac{a_{0}r}{2}+\frac{a_{0}^{2}}{4}\right)[/tex]

Then take this away from what we're currently working out, so that the final answer is:

[tex]<\psi_0|V_\lambda-V|\psi_0>=<\psi_0|V_\lambda|\psi_0>-<\psi_0|V|\psi_0>[/tex]

correct??

nickjer

That is correct (in theory). But you solved the integral wrong. It is the exact same integral as before but without the gamma.

Hart

..huh?

nickjer

The same as this but now alpha = 2/a and not 2/a+gamma.

Hart

.. so you mean for the second part it will be this:

[tex]\int_{0}^{\infty} r e^{-\left(\frac{2}{a}\right)r}dr = \left(\frac{r}{\left(\frac{2}{a}\right)}-\frac{1}{\left(\frac{2}{a})^{2}}\right)\right|^{\infty}_{0} = \frac{4\pi}{\left(\frac{2}{a}\right)^{2}} = \pi a^{2}[/tex]

??

Going back to the first bit, how do I rearrange that result better then? I didn't get what you said about pulling the [itex]\frac{2}{a}[/itex] term out?!