Wavefunction for particle inside a disk

[lolgarithms]

I am trying to solve for the wavefunction for a particle inside a 2 dimensional disk of radius r_0.
The conditions are:
\int^{2\pi}_0 \int^{r_0}_0 |\psi|^2\, r \,dr\, d\theta = 1
Then I tried separations of variables (set psi = R(r)Theta(theta)) to solve the Schrödinger equattion. didnt' really understand how it worked, but tried it anyway.

I got R(r)=A_r r^{i\sqrt{2mE_r}\over\hbar}+B_r r^\frac{-i \sqrt{2mE_r}}{\hbar} and \Theta(\theta) = A_\theta \cos \frac{\sqrt{2mE_\theta}}{\hbar} \theta) +B_\theta \sin \frac{\sqrt{2mE_\theta}}{\hbar} \theta)
(substituting K_i = \frac{\sqrt{2mE_i}}{\hbar}) I plugged the wave function to the probability integral. After painstakingly doing that integral I got:
\left [ |A_\theta|^2 \left ( \pi+\frac {\sin 4\pi K_\theta \theta}{2K_\theta} \right ) +|B_\theta|^2 \left ( \pi-\frac{\sin 4\pi K_\theta \theta}{2K_\theta}\right ) \right ] |A_r^2+B_r^2| \ln r_0 \left ( \frac{(-1+2 \ln r_0)r_0^2}{2}\right ) \left ( \frac{r^{2iK_r+2}}{2iK_r+2}-\frac{r^{-2iK_r+2}}{2iK_r-2} \right ) = 1

Now i can't solve for the energy eigenvalues because they are stuck in the deep parts of the equation (I set E = E_r + E_theta). Also where is the quantization? I know I screwed up somewhere.

 

[PiratePhysicist]

Hi,
Just to make sure, what are you using for you Hamiltonian?

Also, what that condition tells you is that \Psi(r) [\tex] is zero at r_0 (because the particle has a 100% chance of being between 0 and r_0, so it must be zero outside r_0 and wave functions are assuming to be continuous). Your function though does not appear to hold this to be true. You should wait to integrate and set to zero until you've applied boundary conditions.

 

[lolgarithms]

my hamiltonian was -\frac{\hbar^2}{2m} \left ( \psi_{rr} +\frac{\psi_r}{r}+\frac{\psi_{\theta\theta}}{r^2}\right ). I will have to solve for the energy setting R(r0) = R'(r0) = 0.
with these conditions I get:
A_r r_0^{iK_r}+B_r r_0^{-iK_r}=0=A_r iK_r r_0^{iK_r-1}-B_r iK_r r_0^{-iK_r-1}
Nothing cancels. I can't frickin solve this.

 

[PiratePhysicist]

Why are you using something of this form?
<br /> A_r {r_0}^{i K_r}<br />
Try using a function of this form:
<br /> A e^{i x}<br />
(I just woke up, and I'm not generous enough in the mornings to offer hints about what x would be :P). Basically, I mean try something that would be sinusoidal over r.

 

[lolgarithms]

Why are you using something of this form?
<br /> A_r {r_0}^{i K_r}<br />
Try using a function of this form
<br /> A e^{i x}<br />

the laplacian in polar coordinates is (subscripts are derivatives)|\nabla|^2 f(r,\theta}) = f_{rr} + \frac{f_r}{r} + \frac{f_{\theta\theta}}{r^2}. (that's how i write the laplacian because it's an inner product of two dels = "absolute value of del squared". less of an abuse of notation than simply writing "del squared") Thus I will pick a function that will multiply r to the denom when it's differentiated, so that all of the terms in the laplacian are put on a common factor. Thus R(r) can't be an exponential. I guessed a function in the form of R(r)=r^c, c constant. I plugged r^c into the equation for R(r) and got my characteristic equation for c. I solved it and got c=\frac{\pm i\sqrt{2mE_r}}{\hbar}.

 

[PiratePhysicist]

Wait a sec, have you tried Bessel functions? It just occurred to me that if you have x=1/f, then the Schroedinger equation looks similar to Bessel's differential equation.

 

[lolgarithms]

yeah, i think i'll try bessel functions because my complex power function doesn't work for the equation. the dR/dr term doesn't cancel out.
don't spoil the answer for me however.

 

[Avodyne]

I will have to solve for the energy setting R(r0) = R'(r0) = 0.

No, R(r0)=0 only, not also R'(r0)=0. With both the only solution is R(r)=0 for all r. This is like the one-dimensional infinite square well with walls at x=+a and x=-a, where we require the wave function to vanish but not also its derivative.

 

[lolgarithms]

but which bessel function should i use? according to my knowledge there are several types of bessel functions

 

[Avodyne]

First you have to separate variables properly to get two equations, one for R(r) alone and one for Theta(theta) alone. Then you should find that your radial equation has two solutions (that is, two different Bessel functions solve it, before imposing boundary conditions), but that only one of these can be normalized; the other blows up at r=0. Finally, imposing R(r0)=0 will give you a condition that determines the energy eigenvalues in terms of the zeroes of the Bessel function.

 

[lolgarithms]

since i assumed \psi(r,\theta) = R(r) \Theta(\theta)
for the radial part of the wave equation i got r^2 R'' + r R' - (2mE_r / hbar^2)R = 0. this is not a bessel equation, it's an euler equation. so the general solution is not a bessle function. should be
R(r) = A r^(-1 + sqrt(1 - (2mE_r / hbar^2))) + B r^(-1 - sqrt(1 - (2mE_r / hbar^2))).

now, i don't want R to bolw up at 0, so sqrt(1 - (2mE_r / hbar^2)) can't be ±1.

 

[Avodyne]

Your radial wave equation is wrong in two respects. First, the energy term should have a factor of r^2 (and what you call E_r is just the energy E). Second, you are missing the angular-momentum term that comes from introducing a separation constant when separating variables. Your equation (after fixing the first mistake) assumes Theta(theta)=constant.

 

[lolgarithms]

i did not assume the angular component was constant. infact you can see i tried to solve for it:
\Theta(\theta) = A_\theta \cos \frac{\sqrt{2mE_\theta}}{\hbar} \theta +B_\theta \sin \frac{\sqrt{2mE_\theta}}{\hbar} \theta

First, the energy term should have a factor of r^2 (and what you call E_r is just the energy E).

oh really? where does the r^2 come from?

 

[Avodyne]

Inside the disk, the time-independent Schrödinger equation is

-{\hbar^2\over 2m}\nabla^2\psi(r,\theta) = E\psi(r,\theta).

Define k^2 \equiv 2mE/\hbar^2. Then the Schrödinger equation becomes

(\nabla^2+k^2)\psi(r,\theta) = 0.

In polar coordinates,

\nabla^2 = {\partial^2\over\partial r^2} + {1\over r}{\partial\over\partial r} + {1\over r^2}{\partial^2\over\partial\theta^2}.

So if we assume \psi(r,\theta) = R(r)\Theta(\theta), then we have

0 = (\nabla^2+k^2)\psi(r,\theta)

{} = \left[{\partial^2\over\partial r^2} + {1\over r}{\partial\over\partial r} + {1\over r^2}{\partial^2\over\partial\theta^2}+k^2\right]\!R(r)\Theta(\theta)

{}=\Theta(\theta)\!\left[{\partial^2\over\partial r^2} + {1\over r}{\partial\over\partial r}\right]\!R(r) + R(r)\!\left[{1\over r^2}{\partial^2\over\partial\theta^2}\right]\!\Theta(\theta) + k^2 R(r)\Theta(\theta)

{}=\Theta(\theta)[R&#039;&#039;(r)+r^{-1}R&#039;(r)+k^2 R(r)] + r^{-2}R(r)\Theta&#039;&#039;(\theta).

Now divide both sides by r^{-2}R(r)\Theta(\theta); we get

0={r^2[R&#039;&#039;(r)+r^{-1}R&#039;(r)+k^2 R(r)]\over R(r)} + {\Theta&#039;&#039;(\theta)\over\Theta(\theta).

So we have two terms that sum to zero. One is a function of r only, the other is a function of \theta only. If we vary r while holding \theta fixed, the first term can vary, but the second cannot. But since they sum to zero, the first term cannot vary either; it must equal a constant. Call the constant \ell^2. Then we have

{r^2[R&#039;&#039;(r)+r^{-1}R&#039;(r)+k^2 R(r)]\over R(r)} = \ell^2

and

{\Theta&#039;&#039;(\theta)\over\Theta(\theta)}=-\ell^2.

The two linearly independent solutions of the second equation are \Theta(\theta)=\exp(+i\ell\theta) and \Theta(\theta)=\exp(-i\ell\theta). We can treat both at once by taking \Theta(\theta)=\exp(+i\ell\theta) but letting \ell be positive or negative. The angle \theta is periodic with period 2\pi, and \Theta(\theta) must be single valued, so \ell must be an integer (positive, negative, or zero).

The two linearly independent solutions of the first equation are then the Bessel functions J_\ell(kr) and Y_\ell(kr).