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

No solution to harmonic oscillator

  • Thread starter Dustinsfl
  • Start date
  • #1
699
5

Homework Statement


Given [itex](\mathcal{L} + k^2)y = \phi(x)[/itex] with homogeneous boundary conditions [itex]y(0) = y(\ell) = 0[/itex] where
\begin{align}
y(x) & = \frac{2}{\ell}\sum_{n = 1}^{\infty} \frac{\sin(k_nx)}{k^2 - k_n^2},\\
\phi(x) & = \frac{2}{\ell}\sum_{n = 1}^{\infty}\sin(k_nx),\\
u_n(x) &= \sqrt{\frac{2}{\ell}}\sum_{n = 1}^{\infty}\frac{\sin(k_nx)}{k^2 - k_n^2},
\end{align}
[itex]\mathcal{L} = \frac{d^2}{dx^2}[/itex], and [itex]k_n = \frac{n\pi}{\ell}[/itex].
If [itex]k = k_m[/itex], there is no solution unless [itex]\phi(x)[/itex] is orthogonal to [itex]u_m(x)[/itex].


Homework Equations





The Attempt at a Solution



Why is this?

Homework Statement





Homework Equations





The Attempt at a Solution

 
Last edited:

Answers and Replies

  • #2
Simon Bridge
Science Advisor
Homework Helper
17,848
1,645
Well crunch the numbers - what happens if the condition is not satisfied?
 
  • #3
699
5
I don't understand what numbers I need to crunch.
 
  • #4
699
5
For k_m = k_n,
$$
\sum_n\frac{0}{0}\sin(k_nx) = \sum_n\sin(k_nx)
$$
otherwise, k_m\neq k_n
$$
\sum_n\sin(k_nx) = \sum_n\sin(k_nx)
$$
so the DE holds.
We have an indeterminant form on the LHS for the first case and equality for the 2nd. How does orthogonality play a role?
 
  • #5
Simon Bridge
Science Advisor
Homework Helper
17,848
1,645
What does "orthogonal" mean?

What does ##k=k_m## mean when you only have ##k_n## in your expressions?
 
  • #6
699
5
What does "orthogonal" mean?

What does ##k=k_m## mean when you only have ##k_n## in your expressions?
##k=k_m = \frac{m\pi}{\ell}##

##k## is defined in the DE ##(\mathcal{L}+k^2)y = \phi(x)##.

Orthogonal is when the inner product is 0
 
  • #7
Simon Bridge
Science Advisor
Homework Helper
17,848
1,645
And ##k_n = \frac{n\pi}{l}## and ##k_{\text{carrot}}=\frac{\text{carrot}\pi}{l}##... I saw... but why the index "m" when index "n" is used in all the relations above?
What does that mean.......... use words not math?

Similarly for "orthogonal" - what is the relationship between ##u## and ##\phi## when their inner product is zero?
Why is that important to the solutions?

Also curious - how does ##u_n(x)## depend on ##n## when you have it as a sum over all ##n##?
 
Last edited:
  • #8
699
5
And ##k_n = \frac{n\pi}{l}## and ##k_{\text{carrot}}=\frac{\text{carrot}\pi}{l}##... I saw... but why the index "m" when index "n" is used in all the relations above?
What does that mean.......... use words not math?

Similarly for "orthogonal" - what is the relationship between ##u## and ##\phi## when their inner product is zero?
Why is that important to the solutions?
Because it says show that for the case when ##k=k_m## the DE has no solution unless ##\phi(x)## is orthogonal to ##u_m(x)##. The ##k_m## replaces all ##k## so there will be both ##k_m## and ##k_n##.

If there inner product is zero, one of them needs to be a cosine.
 
  • #9
Simon Bridge
Science Advisor
Homework Helper
17,848
1,645
I don't mean just say the math in words.
Try to get your head out of the math - this is a physical situation. What is the physical meaning?
I think you need to review the derivation from it's physical foundation. At the moment it seems to be an exercise in abstract mathematics for you.

What is special about ##k_m## is that it is a particular value of ##k##, which is one of the ##\{k_n\}## ... which means that one of the terms in the sums, the one where n=m, has a division-by-zero problem. That gives a non-physical result - therefore: no solution.

But isn't that a problem? Shouldn't there be some solutions to the situation - after all, you can physically set one up and watch it work away without blowing up or anything? Either the premise the equations were built on is just wrong or there's something else like............

........ what if that term where n=m did not exist for some reason? Something like that? Then the problem would never come up.

What roles do ##\phi## and ##u_n## play in the solution?
In what situation would you normally care about their orthogonality?
 
  • #10
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
14,623
1,257

Homework Statement


Given [itex](\mathcal{L} + k^2)y = \phi(x)[/itex] with homogeneous boundary conditions [itex]y(0) = y(\ell) = 0[/itex] where
\begin{align}
y(x) & = \frac{2}{\ell}\sum_{n = 1}^{\infty} \frac{\sin(k_nx)}{k^2 - k_n^2},\\
\phi(x) & = \frac{2}{\ell}\sum_{n = 1}^{\infty}\sin(k_nx),\\
u_n(x) &= \sqrt{\frac{2}{\ell}}\sum_{n = 1}^{\infty}\frac{\sin(k_nx)}{k^2 - k_n^2},
\end{align}
[itex]\mathcal{L} = \frac{d^2}{dx^2}[/itex], and [itex]k_n = \frac{n\pi}{\ell}[/itex].
If [itex]k = k_m[/itex], there is no solution unless [itex]\phi(x)[/itex] is orthogonal to [itex]u_m(x)[/itex].


Homework Equations





The Attempt at a Solution



Why is this?
Where did your expressions for ##u_n(x)##, ##\phi(x)##, and ##y(x)## come from? As Simon noted, your expression for ##u_n(x)## doesn't make sense. You probably meant, as you wrote in another thread, that ##u_n(x) = \sqrt{2/l} \sin k_n x##. The expression for ##\phi(x)## is clearly not orthogonal to ##u_m(x)##, so I don't see how it relates to your question.

Go back to the beginning. You're trying to solve the differential equation ##(\mathcal{L}+\lambda)y(x) = \phi(x)## where ##\mathcal{L}## is some Hermitian operator. You have solutions ##u_n(x)## to the homogeneous equation ##(\mathcal{L}+\lambda)y(x)=0## that satisfy the given boundary conditions. In the particular case of the linear oscillator, where ##\mathcal{L} = d^2/dx^2## and the boundary conditions are ##y(0)=y(l)=0##, you find that ##u_n(x) = \sqrt{2/l}\sin k_n x## where ##\lambda_n=k_n^2## and ##k_n = n\pi/l##.

One important property of Hermitian operators is that the eigenfunctions ##\{u_n\}## of such an operator form a complete set. This means you can express ##y## as a linear combination of the eigenfunctions, i.e.,
$$y(x) = \sum_n a_n u_n(x).$$ Try doing that for ##\phi(x)## as well and plugging the two expansions into the differential equation with ##k=k_m##.
 

Related Threads on No solution to harmonic oscillator

Replies
18
Views
2K
Replies
9
Views
2K
Replies
6
Views
846
Replies
2
Views
2K
  • Last Post
Replies
7
Views
4K
  • Last Post
Replies
1
Views
1K
  • Last Post
Replies
12
Views
2K
Replies
3
Views
1K
Top