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
How do we get an equation for H(p(Ф), p(z), Ф, z)?
Substitute for Ф’ and z’ using equations found before for p(Ф) and p(z)
What do we do after finding H(p(Ф), p(z), Ф, z)? What does this mean?
-Apply Hamiltons equations: p(Ф)’ = -dH/dФ = 0, so p(Ф) = const
p(z)’ = -dH/dz = -kz
-Since p(Ф) = const, if it is rotating, it will keep moving (p(Ф) is conserved)