QM approximation (electron within nucleus, Griffiths 4.45b

Saraphim

1. The problem statement, all variables and given/known data
What is the probability that an electron in the ground state of hydrogen will be found inside the nucleus?

a) First calculate the exact answer, assuming the wave function $\psi(r,\theta,\phi) = \frac{1}{\sqrt{\pi a^3}} e^{-r/a}$ is correct all the way down to r=0. Let b be the radius of the nucleus

b) Expand your result as a power series in the small number $\epsilon = 2b/a$ and show that the lowest-order term is the cubic: $P \approx (4/3)(b/a)^3$. This should be a suitable appoximation, provided that b << a, which it is.

3. The attempt at a solution
Using partial integration a few times, I've gotten an answer for a) which I think is sufficient in order to solve b). $P (r<b) = 1-e^{-2b/a}(2b^2/a^2 + 2b/a + 1)$.

However, I'm completely lost on b). I'm looking for pointers here on how to even attack the problem. I can rewrite P so that I get $P = 1 - e^{-\epsilon}(\epsilon/2 + \epsilon + 1)$ but from there I'm stumped. I think the problem is that I've mostly forgotten how to do Taylor expansion, and a quick lookup does not help me in attacking this specific problem.

Related Introductory Physics Homework News on Phys.org

tiny-tim

Homework Helper
Hi Saraphim!
… show that the lowest-order term is the cubic: $P \approx (4/3)(b/a)^3$

$P (r<b) = 1-e^{-2b/a}(2b^2/a^2 + 2b/a + 1)$.
use the definition of ex (which is also its taylor series) …

ex = 1 + x + x2/2 + x3/6 + …​

then P = 1 - e-2ε(1 + 2ε + 2ε2)

= 1 - (1 - 2ε + 4e2/2 - 8ε3/6 + …)(1 + 2ε + 2ε2)

= … ?

Saraphim

Thanks for your response! I have a few problems, still. First, I feel the Taylor expansion should be $1 - \epsilon + \epsilon^2/2 - \epsilon^3/6 +\ ...$. Isn't that right?

Second, I feel that I should be able to realize that the terms lower than the cubic somehow get cancelled out, but I really don't see it. :(

tiny-tim

Homework Helper
Good morning!
Thanks for your response! I have a few problems, still. First, I feel the Taylor expansion should be $1 - \epsilon + \epsilon^2/2 - \epsilon^3/6 +\ ...$. Isn't that right?
that would be correct for e

the question gives you e-2ε
Second, I feel that I should be able to realize that the terms lower than the cubic somehow get cancelled out, but I really don't see it. :(
should cancel out now!

Saraphim

But how is that so? In the unexpanded P without substituting in $\epsilon = 2b/a$, I have $e^{-2b/a}$.

tiny-tim

Homework Helper
oh, i'm sorry, i was using ε = b/a

with ε = 2b/a,

then P = 1 - e(1 + ε + ε2/2)

= 1 - (1 - ε + e2/2 - ε3/6 + …)(1 + ε + ε2/2)

Saraphim

I think I understood what I was missing. I needed to set all terms in the product with exponents higher or equal to 4 to zero. Approximation is not my string suite. Thanks for the help!

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving