• Support PF! Buy your school textbooks, materials and every day products Here!

2D Harmonic Oscillator

This might be another problem that our class hasn't covered material to answer yet - but I want to be sure.

The question is the following:
Find the eigenfunctions and eigenvalues of a two-dimensional isotropic harmonic oscillator.

Again, I need help simply starting.
 

Answers and Replies

jtbell
Mentor
15,395
3,181
Have you done the one-dimensional harmonic oscillator? Have you seen any other two- or three-dimensional setups yet? Have you used the method of separation of variables? (e.g. for the two-dimensional rectangular infinite square well)

If you've covered those topics, you should have all the tools you need. Write down the potential energy function for the two-dimensional oscillator, stick it into the two-dimensional Schrödinger equation, and separate the variables to get two one-dimensional equations.
 
We have covered the 1D harmonic oscillator, but we haven't seen any other higher dimensional setups yet. We have also used the seperation of variables so far, just not in regards to higher-dimensions.

Just as a general question - once the equation is broken down into two 1D equations, how are the eigenvalues and eigenfunctions obtained? Is it:
[tex]
H \Psi = E \Psi
[/tex]
with Psi being the eigenfunctions and E being the eigenvalues?
 
238
0
No, you first write the Hamiltonian operator (and use the SE) simply through operator of momentum and potential energy of the oscillator. Only now you will have a PDE, as both operators have two variables. Then use the separation of variables method.
 
Okay, I think I've got it then. Is this correct:
[tex]
\hat{H} = \frac{ (p_x)^2 }{2m} + \frac{ (p_y)^2 }{2m} + \frac{mw^2}{2} \left( x^2 + y^2 \right)
[/tex]
Which is broken up into components:
[tex]
\hat{H} = \hat{H_x} + \hat{H_y}
[/tex]
Noting the 1-D harmonic oscillator gives:
[tex]
E_x = \hbar w \left( n_x + \frac{1}{2} \right)
[/tex]
[tex]
E_y = \hbar w \left( n_y + \frac{1}{2} \right)
[/tex]
[tex]
\psi_n (x) = A_n (a_+)^n \psi_0 (x)
[/tex]
[tex]
\psi_n (y) = A_n (a_+)^n \psi_0 (y)
[/tex]
and noting the seperation of variables on the wavefunction:
[tex]
\Psi (x,y) = \Psi_x (x) \Psi_y (y)
[/tex]
Putting this together forms the results - first, using:
[tex]
\hat{H} \Psi = E \Psi
[/tex]
and plugging in the components gives:
[tex]
\left( H_x + H_y \right) \Psi_x \Psi_y = E \Psi_x \Psi_y
[/tex]
So the eigenvalues are:
[tex]
\hbar w \left( n_x + \frac{1}{2} \right) + \hbar w \left( n_y + \frac{1}{2} \right)
[/tex]
Which simplifies to:
[tex]
\hbar w \left( n + 1 \right)
[/tex]
And the eigenfunctions are simply:
[tex]
\psi_n (x) \psi_n (y) = (A_n)^2 (a_+)^{2n} \psi_0 (x) \psi_0 (y)
[/tex]
Is this right - or do I have something wrong here? Thanks for the help.
 
Last edited:
Gokul43201
Staff Emeritus
Science Advisor
Gold Member
6,987
14
Looks good too (though you haven't said anywhere what n is).
 
with n = 0,1,2......
correct?
 
Gokul43201
Staff Emeritus
Science Advisor
Gold Member
6,987
14
Yes.

Also [itex]p_x[/itex] should be [itex]\hat{p_x} = \frac {i \hbar}{\sqrt{2m}}\frac {\partial}{\partial x} [/itex]
 
185
3
So the eigenvalues are:
[tex]
\hbar w \left( n_x + \frac{1}{2} \right) + \hbar w \left( n_x + \frac{1}{2} \right)
[/tex]
Which simplifies to:
[tex]
\hbar w \left( n + 1 \right)
[/tex]
And the eigenfunctions are simply:
[tex]
\psi_n (x) \psi_n (y) = (A_n)^2 (a_+)^{2n} \psi_0 (x) \psi_0 (y)
[/tex]
This is not quite right.

The eigenvalues are right, but there's really no need to simplify
down to n. You have a set of eigenfunctions in x and y with independent
values, and two indices.

[tex] \Psi_{nm} = \psi_n(x) \psi_m(y) = C_{nm} (a_+^n \psi_0(x)) (b_+^m \psi_0(y)) [/tex]

where the "b" operators are exactly analagous to the "a" operators but operate
in y instead of x. So if we want the eigenfunctions that give [itex] p \hbar \omega [/itex]
we need combinations of n and m that add up to p. (That is, there isn't
one eigenfunction for a given "n".)
 
I see - so technically the eigenvalues are:
[tex]
\hbar w \left( n + m + 1 \right)
[/tex]
and the eigenfunctions are:
[tex]
\Psi_{nm} = C_{nm} (a_+^n \psi_0(x)) (b_+^m \psi_0(y))
[/tex]
with:
[tex]
C_{nm} = A_n A_m
[/tex]

right?
 
Last edited:
Gokul43201
Staff Emeritus
Science Advisor
Gold Member
6,987
14
Oops ! Didn't notice that glitch.

cyber : If you are using indexes n,m for the wavefunction, use the same indexes for the eigenvalues. n is your n_x and m is your n_y.

Just also noticed in #5 that an n_y changed into an n_x when writing the eigenvalue.
 
Last edited:
Corrected the problems you pointed out Gokul. Other than those minor issues though, my solution proposed in #10 is correct then?
 
Gokul43201
Staff Emeritus
Science Advisor
Gold Member
6,987
14
#10 is now correct.
 

Related Threads for: 2D Harmonic Oscillator

Replies
2
Views
2K
  • Last Post
Replies
0
Views
2K
  • Last Post
Replies
11
Views
5K
  • Last Post
Replies
8
Views
2K
  • Last Post
Replies
3
Views
2K
  • Last Post
Replies
3
Views
1K
  • Last Post
Replies
9
Views
2K
  • Last Post
Replies
9
Views
690
Top