Is the 2s Hydrogen Atom Wavefunction Normalized?

[Reshma]

For this given wavefunction of a hydrogen atom in 2s state, verify if the function is normalized:
\psi_{200} = \frac{1}{\sqrt{32\pi a^3}}\left(2 - \frac{r}{a}\right)e^{\frac{-r}{2a}}

My work:
I have to verify:
\int_{all space} \psi_{200}^2 dV = 1
dV = 4\pi r^2dr

So,
\int_{all space} \psi_{200}^2 dV

= \frac{1}{8a^3}\int^{\infty}_{0} \left(2 - \frac{r}{a}\right) ^2 e^{\frac{-r}{2a}} r^2dr

This integral looks like a monster to evaluate . Someone help me out here!

 

[StatusX]

It's not that bad. Just integrate by parts. If you want a shortcut, look up the gamma function.

 

[Reshma]

Thanks, Gamma function reduced my work a lot. I used:
\int_0^{\infty} e^{-\lambda x} x^n = \frac{1}{\lambda^{n + 1}} \Gamma (n + 1)

I was able to get the value of the integral as 1.
Now, for the second part of the problem. I have to find the value of 'r' at which the probablity P of finding the electron is maximum.
P = 4\pi r^2 \psi^2

P = \frac{r^2}{8a^3} \left(2 - \frac{r}{a}\right)^2 e^{\frac{-r}{a}}

I am supposed to get r = 4a by solving \frac{dP}{dr} = 0

This is the expression I got for dP/dr:
\frac{8r^3}{a^2} - \frac{16r^2}{a} + 8r - \frac{r^4}{a^3} = 0

I am unable to solve this equation to get the value of "r". Someone please look into this and let me know where I am wrong.

 

[Physics Monkey]

Try plotting your function is as a function of r/a and you will clearly see that P(r) is not maximized at r/a = 4. I think somebody has made a mistake here, but the wavefunction you gave is correct for the 2s state. The supposed answer must be wrong. It possible to obtain a quadratic equation for the maximum if you use some information about the wavefunction to find one of the factors.

 

[siddharth]

Try plotting your function is as a function of r/a and you will clearly see that P(r) is not maximized at r/a = 4. I think somebody has made a mistake here, but the wavefunction you gave is correct for the 2s state.

Exactly.In fact, if you forget to multiply by r^2(a very common and easy to commit error), you get the incorrect answer r/a = 4.
Also, in
P = \frac{r^2}{8a^3} \left(2 - \frac{r}{a}\right)^2 e^{\frac{-r}{a}}
take the r^2 inside and then differentiate. You will get a quadratic after cancelling the common terms.

 

[Reshma]

Thanks Physics Monkey and Siddharth.
I did take the r² term inside and differentiated it, but I don't get a quadratic term. I have posted the equation I got in post #3. Anyway, I'll post again:
\frac{8r^3}{a^2} - \frac{16r^2}{a} + 8r - \frac{r^4}{a^3} = 0
Is it possible to solve this equation or am I going wrong somewhere?

 

[dextercioby]

One solution is r=0. Obviously non physical. Then you'd get a cubic equation which you can treat using Cardano's formulae.

Daniel.

 

[Physics Monkey]

As dextercioby already noted, there is a common factor of r that you can obviously factor out. A simple way to make progress on the cubic is to note the following: the probability density passes through zero at r = 2a, but it is also always greater than or equal to zero. What does this tell you about one of the roots of your equation?

 

[siddharth]

Wait a minute. You essentially have to differentiate
(2r - \frac{r^2}{a})^{2} e^{\frac{-r}{a}}
From that, won't you get
2(2r-\frac{r^2}{a}) (2-\frac{2r}{a}) e^{\frac{-r}{a}} + (\frac{-1}{a})(2r - \frac{r^2}{a})^2 e^{\frac{-r}{a}} =0And simplifying that,
r^2 -6ar + 4a^2 = 0

 

[Reshma]

As dextercioby already noted, there is a common factor of r that you can obviously factor out. A simple way to make progress on the cubic is to note the following: the probability density passes through zero at r = 2a, but it is also always greater than or equal to zero. What does this tell you about one of the roots of your equation?

Yes, I can factor out the common 'r'.
\frac{8r^2}{a^2} - \frac{16r}{a} - \frac{r^3}{a^3} +8 = 0

Probability density is greater than or equal to zero. Does this mean the roots are real?

 

[Physics Monkey]

The probability density is always greater than or equal to zero. Therefore, when the density equals zero that point will be a ... ?

 

[Reshma]

The probability density is always greater than or equal to zero. Therefore, when the density equals zero that point will be a ... ?

The point will be a minimum, right?

I need to clarify certain doubts or else this thread will become FUBAR.
Just as for particle in the box, there are some positions where the probability is zero. The radial probability distribution functions P(r) indicate the relative probability of finding the electron within a thin spherical shell of radius 'r'.
So if the P(r) = 0, the point is a minimum so at what value of 'r' is P(r) maximum?
How does this compare to 4a, the distance between the electron and the nucleus in the n=2 state of Bohr model?

 

[Gokul43201]

Your equation is basically :
x^3 - 8x^2 + 16x -8 = 0

Clearly, if there's a nice integer root, it would have to be an even integer. Inspection gives you a root immediately. That reduces the equation to a quadratic.