Solving the Spherical Oscillator Problem

T-7
Messages
62
Reaction score
0
Hi,

I'd really appreciate some speedy advice on this question. I need to know whether the expressions for energy and the angular momentum are correct, as I seem unable to produce the final expression desired from what I have derived.

Here's the problem:

Homework Statement



The Lagrangian is

L = \mu \nu^{2}/2 - kr^{2}/2

i. Calculate the angular momentum vector, L, proving that L is constant and that, unless L vanishes, the oscillation occurs in the x,y plane.

ii. Using x = x_{0}cos(\omega t), y = y_{0}sin(\omega t), show the trajectory is an ellipse. Express the energy, E, and then both the eccentricity and the parameter in terms of E and Lz.

iii. Representing the trajectory in polar coordinates by x = rcos(\phi), y = rsin(\phi), show that you get the simple relation

<br /> cos (2\phi) = \frac{1 - \nu^{2} / (2\xi^{2})}{\sqrt{1 - \nu^{2}}} ,<br /> r = \sqrt{ \frac{2E}{k} }\xi ,<br /> \nu = \frac{ L_{z}\omega }{E}<br />

The Attempt at a Solution



For i. I found that L = Lz = \mu \omega x_{0}y_{0}sin(\phi_{x}-\phi_{y}) (assuming L is not zero), by using three sine equations for the x,y and z components (with angle offsets) and taking L = m(r \times v); I chose my axis so that x and y could be rewritten as x = x_{0}cos(\omega t), y = y_{0}sin(\omega t), and then find I can write Lz = \mu \omega x_{0}y_{0}.

For ii., it seemed all that was required was to square both expressions and combine thus:

\frac{x}{x_{0}}^{2} + \frac{y}{y_{0}}^{2} = 1

Writing the energy E = \mu \nu^{2}/2 + kr^{2}/2 and substituting the expressions for x and y, I found the expression to boil down to E = \frac{(\mu \omega^{2})}{2} (x_{0}^{2}+y_{0}^{2}) where \omega^{2} = k/m.

I'm guessing I need to go ahead and solve this equation for x0 (assuming it's right) by using Lz to get rid of y0, and hope that something which looks like the eccentricity and the parameter will emerge (?).

Part iii. has me really stuck. I substitute x = rcos(\phi), y = rsin(\phi) into the ellipse equation, which looks like a promising start, use double angle results to turn the cos and the sine into cos (2\phi), and eventually obtain (after some simply rearranging):

cos (2\phi) = \frac{1-r^{2}/2(1/x_{0}^{2}+1/y_{0}^{2})}{r^{2}/2(1/x_{0}^{2}-1/y_{0}^{2})}

Multiplying numerator and denomiator by x_{0}y_{0} and using my expressions derived earlier for Energy and Lz I arrive at:

cos (2\phi) = \frac{L_{z}^{2}/\mu^{2}\omega^{2}-r^{2}/2(2E/\mu\omega^{2})}{r^{2}/2(y_{0}^{2}-x_{0}^{2})}

It looks promising, but once the first term on the numerator has been made one, it seems no manipulation can make this equation correspond to the one given. I can only presume my energy, angular moment or ellipse equations are wrong? But are they?

Many thanks for your insights!
 
Last edited:
Physics news on Phys.org
Hey, Can I just ask how you got the expression for E=\frac{(\mu\omega^2)}{2}(x_0^2+y^2_0)
 
I've been working on this problem as well, and I think that the expression for the energy should have a k term

E=\frac{(\mu\omega^2)}{2}(x_0^2+y^2_0)+\frac{k^2}{2}}(x_0^2+y^2_0)

Because if you put the terms for xo and yo in for r I can't see why the k term would disappear. Also for your expression for the ellipse the denominators should be squared.

I haven't gone any further with this yet, but the method seems like it should work from what you've explained.
 
Can I just ask if we are all handing this in tomorrow for a certain classical mechanics module at a certain Scottish University?

I tried using Maple to get the expression for E which comes out as:
E=\frac{k}{2}(x_0^2+y_0^2)

I can't seem to get very far with the next question and haven't even looked at the last one. Oh Classical Mechanics!

has anyone had any luck getting the eccentricity out?
 
Last edited:
Hi, I had an exam and I completely messed up a problem. Especially one part which was necessary for the rest of the problem. Basically, I have a wormhole metric: $$(ds)^2 = -(dt)^2 + (dr)^2 + (r^2 + b^2)( (d\theta)^2 + sin^2 \theta (d\phi)^2 )$$ Where ##b=1## with an orbit only in the equatorial plane. We also know from the question that the orbit must satisfy this relationship: $$\varepsilon = \frac{1}{2} (\frac{dr}{d\tau})^2 + V_{eff}(r)$$ Ultimately, I was tasked to find the initial...
The value of H equals ## 10^{3}## in natural units, According to : https://en.wikipedia.org/wiki/Natural_units, ## t \sim 10^{-21} sec = 10^{21} Hz ##, and since ## \text{GeV} \sim 10^{24} \text{Hz } ##, ## GeV \sim 10^{24} \times 10^{-21} = 10^3 ## in natural units. So is this conversion correct? Also in the above formula, can I convert H to that natural units , since it’s a constant, while keeping k in Hz ?
Back
Top