Hamilton's Equations and more examples Flashcards
What are Hamiltons equations?
qi’ = dH/dpi, pi’ = -dH/dqi
What do we start with for the change in Lagrangian dL?
Sum over i of dL/dqi dqi + dL/dqi’ dqi’ = sum over i of dpi/dt dqi + pi dqi’
What do we do with this version of the Lagrangian?
- d(sum over i of pi*qi’) = sum over i of dpi qi’ + sum over i of pi dqi’
- sum over i of pi dqi’ = dL - sum over i of dpi/dt dqi
What is the third step after rearranging the Lagrangian again?
- dL = sum over i of (pi’ dqi 0 qi’ dpi) + d(sum over i of pi qi’)
- Rearrange so d(sum over i of pi qi’ - L) = sum over i of (qi’ dpi - pi’ dqi), so dH = sum over i of qi’ dpi - pi’ dqi
- Need H(pi, qi), so dh = sum over i of dH/dpi dpi + sum over i of dH/dqi dqi
- Compare thee two equations for dH
For the harmonic oscillator example, how to we convert H(x, v) to H(px, x)?
- use v = px/m to re-write H
- apply hamiltons equations
- solve the equations we get from these
What is phase space?
- For n qi’s, imagine space that has 2n dimensions
- Position in phase space is (p1, p2……pn; q1, q1,……..qn)
- A point in phase space gives us the current state of the system
- Path/line tells us the evolution of the system
For the 1D harmonic oscillator example, how can we use the Hamiltonian to plot px against x?
H = px^2/2m + 1/2kx^2 = const, which is equation for an ellipse, so can plot an ellipse with known axis-crossing points
What is the equation for the Lagrangian for a pendulum?
L = T-V = 1/2ml^2 *θ’^2 - mgl(1-cosθ)
What does p(θ) equal and what can we use this for?
- p(θ) = dL/dθ’ = ml^2θ’
- Can use this in the Hamiltonian: H = θ’ *dL/dθ’ - L = 1/2 ml^2 *θ’^2 + mgl(1-cosθ)
- Then sub in p(θ) to get H(p(θ), θ) = E
What do we do after finding H(p(θ), θ)?
- Use Hamilton’s equations: p(θ)’ = -dH/dθ = -mglsinθ, and θ’ = dH/dp(θ) = p(θ)/ml^2
- Can use this to plot phase space picture
What do we find in terms of the phase space for small values of E?
- Small oscillations so sinθ ~ θ, cosθ ~ 1-θ^2/2
- Sub in θ for sinθ in hamiltons equations, and then set p(θ) = 0 for one solution and θ=0 for the other solutions
- Plot θ against p(θ) for phase space graph.
What do we find in terms of the phase space for very large values of E? (enough for pendulum to rotate completely)
- sinθ = 0 when θ=0, π, 2π etc
- sketch p(θ) - θ curve by considering θ=0 in Hamiltons equations and θ=π in Hamiltons equations
- Get graph symmetric in θ axis looks like very shallow sine curve going from +π to -π
For a particle, mass m, moving on surface of a cylinder radius R, with F = -kr, what does V equal?
V = 1/2k(x^2 +y^2 + z^2) = 1/2k(R^2 +z^2)
For a particle, mass m, moving on surface of a cylinder radius R, with F = -kr, what does T equal?
T = 1/2m(R^2 *Ф’^2 + z’^2), where R term in KE in x-y and z term in kinetic energy in z
What do we do after determining the Lagrangian for the particle moving on the surface of a cylinder problem?
- Get p(Ф) and p(z): p(Ф) = dL/dФ’, p(z) = dL/dz’
- H = sum over i of qi’*pi - L = Ф’ p(Ф) + z’ p(z) - L