Units of radial wavefunction

I am solving the Hydrogen wavefunction using FORTRAN.
Now using the Euler method, I am given a solution to match which is given by u10(r) = 1.06r*exp(-3.74r) (where unl(r) = rRnl(r) in general) which says it has a normalisation chosen to match what i should get from my code.

Then I use the backwards Euler method getting a more accurate solution with a correct shape and a much larger amplitude. I then normalise this and plot it against the analytic solution given by
u = 1/pi^1/2*(1/ao)^3/2*r*exp(-r/a0) which is apparently also normalised.

I am working in natural units so h-bar = c = 1 and me = 511.7keV the electron mass, which implies that my Bohr radius = 0.26737keV^-1 and I am working in units of r given by keV^-1 also.

so this would imply that the units of u for the analytic solution must be keV^(1/2) however looking at the original equation shouldn't the units be keV^-1 as u = rR(r) where r has units keV^-1 and R(r) is dimensionless, I'm confused. Is my normalisation unit calculated in fortran supposed to have units like the analytic normalisation constant?

On my graphs should I just put no units? :S
Any suggestions would be appreciated, units of wavefunctions has always confused me.
Last edited:


The radial part of the wave function is normalized according to
\int_0^\infty R_{nl}(r) r^2 \, dr = 1
which means that ##R_{nl}## has the same units as ##r^{-3}##.

Want to reply to this thread?

"Units of radial wavefunction" You must log in or register to reply here.

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

Hot Threads