Green function for forced harmonic oscillator

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
3 replies · 3K views
Judas503
Messages
22
Reaction score
0

Homework Statement


The problem requires to solve the integration to find ## G(t) ## after ##G(\omega)## is found via Fourier transform. We have [tex]G(\omega)= \frac{1}{2\pi}\frac{1}{\omega _{0}^2 - \omega ^2}[/tex]

Homework Equations


As mentioned previously, the question asks to find ##G(t)##

The Attempt at a Solution


It is obvious that calculus of residues is required. To account for causality (## G(t<0)=0 ##), the poles at ## \omega=\pm \omega_{0} ## are shifted to the lower half plane by ## i\epsilon ## and integrated along the contour in the lower half plane. Then,
[tex]G(t)=-\frac{1}{2\pi}\int_{-\infty}^{+\infty}\frac{e^{-i\omega t}}{(\omega +i\epsilon)^2 -\omega _{0} ^2}d\omega[/tex]
After calculating the residues at ## \omega =\pm \omega _{0} - i\epsilon ##, I found
[tex]G(t)=\frac{i}{2\omega_{0}}\sin \omega_{0}t[/tex]

Is my answer correct?
 
Physics news on Phys.org
The presence of ##i## suggests it's not. G(t) should be real, shouldn't it?

The sine makes sense. The system is at rest and then you impart an impulse, causing it to oscillate. It has to be sine because it starts from x=0.
 
It looks like the i and the 2 might both get absorbed into the ##\sin \omega_0 t =\frac{e^{i \omega_0 t }-e^{-i \omega_0 t }}{2i}## term.
 
Yes, sorry! I forgot to put the "i" in the exponential form of sine. That should clear the problem.