Phys410 - Classical Mechanics
University of Maryland, College Park
Fall 2011, Professor: Ted Jacobson
Homework
Problems from Taylor,
Classical Mechanics
HW0 - 7due at
beginning of class, Tuesday 9/06/11
HW12 -
due at beginning of class, Tuesday 12/13/11
16.31 (Earthquake
distance)
S12.1 (Fourier modes of
string in Lagrangian and Hamiltonian)
The Lagrangian for a string of length a, linear mass density µ, and tension
T, held fixed at the ends,
L = ∫ 1/2 [µ (∂y/∂t)2
- T (∂y/∂x)2]
dx
where y(x,t) describes the transverse displacement of the
string.
(a) Instead of using the function y(x) as the generalized
coordinate, rewrite the Lagrangian using the coefficients y_n, n
= 1,2,3,...
of the Fourier series y(x, t) = ∑ y_n(t) sin k_n x. (I'll let
you figure out how k_n depends on n and the string length.)
Show that the Lagrangian is a sum of harmonic oscillator
Lagrangians, one for each n. Find the frequency w_n of each
oscillator, and show that w_n = n w_1.
(b) Write out and solve Lagrange's equations for the system,
treating the y_n as the generalized coordinates.
(c) Find the canonical momentum conjugate to each y_n, and use
that find the Hamiltonian for the string using the phase space
of
these Fourier components and their conjugate momenta.
(d) You can "bend" a guitar note by pushing on the string to
increase the tension while the string is sounding. (The length
also changes,
but I think the effect of this is smaller.) (The note is a
superposition of all the harmonics, dominated by the
fundamental.)
(i) Explain why the variation of the tension is an adiabatic
change of the string Hamiltonian.
(ii) If the note is bent upwards by a whole step (2/12 of an
octave), by what percent does the energy of the string vibration
increase?
(Hint: The total energy
is the sum of the energies of all of the modes.)
(e) A quantum harmonic oscillator has a zero point energy hbar
w/2 in the ground state, and a string is a collection of
oscillators of
frequencies w_n = n w_1. Explain why the string doesn't have an
infinite zero point energy!
S12.2 (String with a
bending modulus)
In class we considered a stretched string whose only restoring
force is due to the tension. Suppose the string is somewhat
stiff and
has a potential energy associated with bending. The bending is
characterized by the curvature of y(x), i.e. ∂2y/∂x2.
The bending energy
has the form
∫ 1/2 ß (∂2y/∂x2)2
dx, where ß is
a bending modulus.
(a) Find the Lagrange equation for the stretched string
including this additional term, by imposing the condition that
the action is
stationary under all variations of y. It should be a modified
wave equation, with a fourth order spatial derivative term.
(b) The Fourier modes of the string are still harmonic
oscillators in the presence of the bending term.
Inspect your analysis for the Fourier series in S12.1 and read
off the frequencies w_n when the bending modulus is included.
You should find that the frequencies grow with n more quickly
than without the bending modulus. The string is no longer
"harmonic", in the sense that the frequencies w_n are not
integer multiples of the fundamental w_1.
S12.3 (Relativistic
string) A relativistic string has a rest
energy is equal to the string tension T times the proper length
of the string
(so its rest energy changes as it oscillates). The
action for a relativistic free particle is -mc^2 times the
proper time along the worldline.
If you replace mc^2 by the tension T times the proper length of
the string at each time, you get the action for the string: it
is -T times
the integral over proper length and proper time along the string
"worldsheet", i.e. -T times the proper spacetime area of the
"worldsheet".
Rather than ask you to work out the equations for a general
motion, I'll restrain myself and ask you only to consider
circular
configurations of a closed loop of string whose center is at
rest in some inertial frame, so the string is
characterized at a given time by
its radius R(t) alone, where t is the time coordinate in that
frame. The action discussed above, when restricted to such
configurations,
is -T ∫ dA = -2πT ∫ R [1 - (dR/dt)2/c2]1/2
dt. Let's choose units with c = 1 and
2πT = 1, so the action is just
S = - ∫ R [1 - (dR/dt)2]1/2
dt.
(a) Find the momentum p conjugate to R.
(b) Solve for dR/dt in terms of p and R.
(c) Find the Hamiltonian for the string, and compare it to the
Hamiltonian for a relativistic free particle.
(d) Find the equation of motion.
(e) (i) Explain how you know in advance that the Hamiltonian is
conserved, and (ii) use the equations of motion to show it
explicitly.
(f) Using (d) and (e), show that the string satisfies a harmonic
oscillator equation, and identify the frequency of the
oscillator.
(g) Show that when the string radius goes to zero, the string is
moving at the speed of light.
(h) A quantized string in its ground state has some zero point
motion. Estimate the size of the R fluctuations in the ground
state
by minimizing the energy subject to the uncertainty relation p R
≥ hbar. Use units with hbar = 1. Then find this zero point size
using
general units, i.e. restore the appropriate factors of hbar, c
and, 2πT
HW11
- pdf file - due at beginning of class, Tuesday
12/06/11
HW10 - due at
beginning of class, Tuesday 11/29/11
13.12 (bead on spinning
rod)
13.13 (particle on cylinder
with restoring force)
13.14 (particle on cone)
(Note: It should say
"for a given energy [and angular momentum], this occurs at
exactly two values of z,
[except for the special case when z is constant]." Also, note
that when p_z = 0 this Hamiltonian is equal to the effective
potential
for a fixed angular momentum.)
13.17 (nearly circular orbits
of particle on cone) (Suggestion:
Use the effective potential to do parts (a), (b), and (c).)
S10.1 (Hamiltonian for
relativistic particle and cyclotron motion) The action for a
relativistic particle is -m∫ dt (1-v2)1/2,
where v = dx/dt,
i.e. the Lagrangian is L = -m
(1-v2)1/2.
(a) Find the momentum p
conjugate to x. (b)
Show that p2
+ m2 = γ2m2,
where
γ = (1-v2)-1/2.
(c) Show that the Hamiltonian is H = (p2
+ m2)1/2
(in units with c = 1). (d) Now include an external
electromagnetic field, and show that for a charge e the
relativistic Hamiltonian is H = [(p - eA)2
+ m2]1/2+
eV,
where V and A are the
scalar and vector potentials. Note that in a magnetic field (with no
scalar potential),
the Hamiltonian is given by H = γm.
(e) Suppose there is a uniform magnetic field in the z direction,
which can be
derived from the vector potential A
= Bx yhat. Show that the
orbits in the xy plane are circles with angular velocity w = w_0/γ,
where w_0 = eB/m is the nonrelativistic cyclotron
frequency. In your analysis you may for convenience set
the conserved quantity
p_y to zero (after finding Hamilton's equations), since it just
determines the x coordinate of the center of the circular orbit.
(f) Eliminate the velocity using v = Rw, and show that w =
w_0/(1+R^2 w_0^2)^1/2. Note that this ensures that the speed is
never greater than the speed of light, and that as the radius
grows the speed approaches the speed of light.
S10.2 (Proper time and orbits)
Consider three worldlines that all begin and end at the same
radius above the earth's surface: (i) the circular orbit at that
radius,
(ii) the static worldline held in place by a long pole or a
rocket, and (iii) the radial up and down free-fall path (with
initial velocity
chosen so that the path will return at the same time as the
circular orbit does). Paths (i) and (iii) are both stationary
paths of the
proper time since they are freely falling, while path (ii) is
not. Show that the proper times on these three paths satisfy
tau_i < tau_ii < tau_iii. You need not evaluate the proper
times; rather, only show that they satisfy these inequalities.
For
(i) and (ii) you can easily show it directly (again, without
actually evaluating the times). For (iii) the most convenient
thing
is to argue that there must be a path of longer proper time than
(ii), and that path (iii) must be this path. (If you like a
challenge,
try to actually demonstrate that tau_iii is longer by using the
conserved energy and effective potential. I didn't succeed in
finding
a simple way --- I think it's rather tricky.)
S10.3 (Gravitational redshift)
Consider two light pulses sent radially outward from radius r_1
to radius r_2 in the Schwarzschild spacetime, with the emission
events separated
by a Schwarzschild coordinate time interval ∆t. (a) What is the
Schwarzschild coordinate time interval
between the arrival of the two pulses at r_2?
(b) What is the proper time ∆tau_1 between the emission of the
two pulses at r_1? (c) What is the
proper time ∆tau_2 between the reception of the
two pulses at r_2? (d) Suppose a photon is emitted radially from
r_1 with frequency w_1 measured by a static observer at r_1.
When the photon arrives
at r_2, what will be its frequency as measured by a static
observer at r_2? (e) As r_1 approaches the Schwarzschild
radius, what happens to the frequency
w_2, if w_1 is fixed?
HW9 - due at
beginning of class, Tuesday 11/15/11
S9.1 (Approximate
relativistic relation between energy and momentum)
Consider a particle of rest mass m, energy E and momentum p, and
first set c=1.
(a) Show that when p << m, we have E ≈
m + p2/2m.
(b) Show that when p >> m, we
have E ≈ p
+ m2/2p, and p ≈
E - m2/2E.
(c) Rewrite your answers for (a) and (b) including the
appropriate factors of c.
S9.2 (Inverse Compton
scattering)
If a high energy electron collides with a low energy photon,
much of the electron energy can be transferred to the outgoing
photon.
This is called "inverse Compton scattering", even though it's
the same as Compton scattering when viewed in a frame in which
the electron
is initially at rest. Let's consider the case in which the
incoming electron is highly relativistic (E0 >> m) and
moving in the +x direction,
and collides head on with a photon of frequency w0 moving along
the -x direction. Suppose that after the collision the outgoing
photon
is moving along the +x direction.
(a) Let p0, k0, p, and k, be the initial and final 4-momenta of
the electron and photons. Show that k0.p0 = k.p
(b) Using (a) show that w = w0(E0 +
p0)/(E - p) ≈ 2w0
E0/(E - p) (the last approximation follows from E0 >> m).
(c) Use conservation of energy and momentum to show that
E - p = E0 - p0 + 2w0.
(d) Combine the previous to results and S9.1 to show that w/E0 =
1/[1 + m^2/(4E0 w0)]. When the incoming electron energy is
sufficiently large, the outgoing photon therefore has almost all
of the energy.
(e) If w0 = 1 eV, what must be the incoming electron energy E0
if the final photon energy is to be approximately equal to E0/2?
S9.3 (Motion in the Schwarzschild spacetime)
S9.4 (Motion in a gravitational wave)
HW8 - due at
beginning of class, Tuesday 11/08/11
15.12 (Lorentz
contraction)
Answer only (b) How long is the pipe as "seen" (better word
would be "reckoned") by the pions, and how long does it take to
pass the pions?
15.56 (rest energy in chemical
reaction)
15.57 (rest energy in nuclear
reaction)
15.58 (rest energy vs. kinetic
energy)
15.71 (collision energy to
create a particle)
15.75 (parent particle
reconstruction)
15.87 (pion decay to photons:
pion velocity and photon angle)
15.90 (pair creation off
nucleus) Do part (a) by showing that the sum of two
future timelike 4-vectors cannot be lightlike.
15.92 (pion decay to lepton
and neutrino)
S8.1 (no vacuum Cerenkov
radiation) Show that the reaction e -> e + gamma
cannot satisfy both energy and momentum conservation.
You could do this at least three different ways: by (i) writing
out energy and momentum conservation and applying the mass
relation
in a general frame, (ii) doing the analysis in the rest frame of
the initial electron, or (iii) using 4-vectors and proving the
result for problem 15.53.
S8.2 (cosmic gamma ray cutoff)
Although gamma rays of energies up to 1020 eV are
created by ultra high energy cosmic rays throughout the
universe, they cannot travel all the way to the earth, because
on the way they collide with other photons and disappear,
creating electron-positron
pairs. The reaction is gamma1 + gamma2
-> e+ + e-. (a) If the two
photons collide head-on, show that for a given E2
the minimum energy E1
required to
create the pair is E1 = me2/E2.
(b) Evaluate this energy (in eV) assuming E2
is the energy of (i) a typical cosmic microwave background
photon, 3K = 0.0003 eV,
or (ii) an infrared photon of energy 0.5 eV. (The
actual cutoff of cosmic ray photons originating from more than 100
million light years (?) comes from
annihilation on the infrared background radiation. I'm not sure
about the numbers here...)
HW7 - due at
beginning of class, Tuesday 11/01/11
11.32 (CO2 vibrations) Do the eigenvalue and eigenvector
calculations by hand. (Suggestion:
Once you have worked out the equations
of motion, adopt units with k = m = 1.) Add to part (b):
Which mode has the higher frequency?
Try to explain "why", without reference to equations. Add parts:
(d) Show that the ratio of the two frequencies is √11/3, taking
into account
that the ratio of carbon to oxygen mass is 3/4. Show that the
observed frequency ratio is smaller by
about 12%.
(See www.phy.davidson.edu/StuHome/jimn/CO2/Pages/CO2Theory.htm
for information on carbon dioxide vibrations.*)
(e) Try to give a physical reason (or reasons) why the observed
ratio is smaller than that in our simple model. I don’t think
the difference
between quantum and classical mechanics is the most important
issue in this case.
(Hint: Think about the
physical nature of the “springs” in the molecule and how they
differ from those in our model.)
*The frequencies are given on this web page in units of cm−1.
This refers to the inverse wavelength of the photons that are
emitted in
transitions between the vibrational levels. The inverse
wavelength is proportional to the frequency of the photons,
which is proportional to the
frequency of the vibrations.
15.7 (muon lifetime in cosmic
ray showers)
S7.1 (masses suspended by
springs)
Consider a mass m suspended by a spring with spring constant k
from another mass m which is suspended from a fixed support by
another
spring with spring constant k. Let y_1 and y_2 denote the
displacements from their equilibrium positions of the top and
bottom masses
respectively, with the downward direction taken as positive.
Consider only vertical motion. (a) Write the Lagrangian for the
system, and find
the equations of motion. (Note the
equilibrium spring forces balance the gravitational forces,
which therefore drop out of the problem.)
(b) Determine (by hand) the normal mode frequencies and
displacement ratios y_2/y_1. (To check your result: The squared
frequencies are
(3± √5)/2 and the displacement ratios are y_2/y_1 = (1
-/+ √5)/2.) (c) Describe and indicate with arrows the nature of
the two normal mode motions,
showing both direction and approximate relative displacement of
each mass. Label with the frequency of each mode. Which is the
higher one?
Try to explain "why", without reference
to equations.
S7.2 (physical pendulum)
In class we saw the demonstration D2-13
RACING PENDULA. Show that the motion of a physical
pendulum consisting of a uniform rod of mass
M and length R is equivalent to that of a simple pendulum of
mass M and length 2/3 R.
S7.3 (physical penulum hanging
from a string)
Consider a
uniform rod of mass M and length R hanging by its end from a
massless string of length l. (a) Write the Lagrangian using the
angles of the string
and rod as generalized coordinates, and expand it to quadratic
order in the angles. (b) Find the frequencies of the normal mode
oscillations for motion in a plane.
(c) Evaluate the frequencies in three cases: l=R, l >> R,
and l << R.
(Notes: (i) When you
write the kinetic energy, decompose it into the center of mass
motion and the motion relative to the
center of mass, and for the latter make use of the moment of
inertia relative to the center of mass (not relative to the end
of the rod).
(ii) I suggest you adopt units with R=g=1 for part (b). (iii) Answers for part (c): case 1:
3.1√g/R, 0.8√g/R;
case 2: √6g/R, √g/l; case 3: √3g/2R,
2√g/l.)
S7.4 (relativistic
longitudinal Doppler effect) This problem refers to the
situation and notation with observers O1 and O2 described in the
notes of 10/20.
(a) Explain why the ratio t1/t0 is the Doppler factor for light,
i.e. the ratio of received frequency to transmitted frequency
for light sent from O1 to O2.
(b) Show that t1/t0 = √t1/t2. (c) The the
expression of t1 and t2 in terms of Dx and Dt to show that the
Doppler factor is √(1-v/c)/(1+v/c). (d) For v << c
it is useful to use an approximation. Expand the Doppler factor
to linear order in v/c to find this approximation.
HW6 - due at
beginning of class, Tuesday 10/18/11
9.2 (artifical gravity in
rotating space station)
9.7 (derivative of vector in
rotating frame)
9.8 (direction of centrifugal
and Coriolis forces)
9.11 (Lagrangian in rotating
frame and equations of motion) The textbook suggests to
do this problem using vector notation,
without coordinates. However in this case a judicious choice of
coordinates simplifies matters, or at least eliminates the need
for
fancy footwork with vector identities). So, instead of using
vectors, I suggest you adopt cylindrical coordinates, with the z
axis
along the angular velocity of the rotating frame. Then phi_0 =
phi + Omega t. Once you get the three Lagrange equations, you
can then show they are equivalent to (9.34).
9.22 (Larmor precession via
rotating frame trick)
9.26 (approximate trajectory
of falling particle on spinning earth)
S6.1 (tides for a constant
central force) It is often said - and I probably said
it - that the ocean tides can be traced to the fact that the
moon's gravitational force gets weaker with distance. But
actually we would have tides even if the moon's force on the
earth were
independent of distance, because of the direction dependence.
To see how this works explicitly, let's assume that the force
the moon
exerts is central and everywhere has the same magnitude as the
true gravitational force of the moon has at the center of the
earth.
That is, F =
-(GMm/d_0^2) dhat, where
M is the mass of the moon, m is the mass of the particle it is
acting on, d_0 is the distance from
he center of the moon to the center of the earth, and dhat is a unit vector
pointing from the moon to the particle. For this problem, you
will go through the same steps as in the textbook, but for the
modified force field. (a) Write down the tidal force field (as
in Eq. 9.12).
(b) Sketch the tidal force field at the surface of the earth (as
in Fig. 9.4). (c) Find the tidal potential (as in Eq. 9.13). (d)
Find the height
difference between high and low tides in the idealized model (as
in Eq. 9.18). How does your result compare with the true result
(9.18)?
HW5 - due at
beginning of class, Tuesday 10/11/11
8.2 (two bodies in an
external field) Add parts: (c) If the external field is
not uniform but is, say, a -1/r potential, then the
center of mass motion no longer separates completely from the
relative motion. For instance, this is the case for the
earth-moon system in the field of the sun. Write the Lagrangian
of the earth-moon system using the center of mass
and the earth-moon relative vector as your coordinates. (d) Make
a Taylor expansion in the small ratio of the earth-moon
distance divided by the distance from the center of mass to the
sun, to extract the leading order term in this ratio that
appears in the Lagrangian. (If you're ambitious (not required),
work out the next to leading order term as well.)
(e) Try to guess (after a bit of thought) what the effect of
this extra term is on the earth-moon system.
8.3 (two masses, spring and
gravity) Note:
The solution in the book is only valid if the table is removed
at t=0. You can assume
this is the case, to simplify the problem.
S5.1 (dark matter)
One of the first pieces of evidence for the existence of dark
matter was the "flat rotation curves" of galaxies, meaning the
fact that
at large radii the orbital velocity of stars (as measured by the
Doppler effect) approaches a constant, rather than falling off
with distance. (a) How does the velocity of a test body in a
circular Newtonian gravitational orbit around a central mass
scale
with the radius of the orbit? (b) Suppose the mass of a galaxy
is dominated by a spherical distribution of some kind of dark
matter
with a radial mass density profile rho(r). What form must the
function rho(r) take in order for the orbital velocities of
stars to be
independent of r?
S5.2 (innermost stable
circular orbit in GR)
In general relativity there is a correction to the Newtonian
orbit around a central gravitating body that can be approximated
for weak fields and low velocities by an additional, attractive
1/r^3 term in the potential energy. The total effective
potential for a
body of mass m then has the form U =
-a/r + b/r^2 - d/r^3. Here a = GMm, where M is the mass of
the central body, assumed much
larger than m, b = l^2/(2m), where l is the angular momentum of
the orbiting mass, and d = b r_g, where r_g = 2GM/c^2 is the
"gravitational radius" or "Schwarzschild radius"
(for M = M_sun the gravitational radius is 3 km), and c is the speed
of light.
We can simplify the algebra by choosing units with GM = c = m = 1,
in which case a = 1, b = l^2/2, and d = l^2.
Assume you can use Newtonian dynamics.
(a) Are there any unstable circular orbits? If so, over what range
of radii do they exist?
(b) Find the minimum value of r for which a stable circular orbit
exists, and find its orbital speed. (That orbit is called the ISCO,
"innermost stable circular orbit".) Use the above units in your
algebra, but once you have the result re-express the radius of the
ISCO
in terms of the gravitational radius and the orbital speed in terms
of the speed of light.
S5.3 (precession of perihelion of
Mercury in GR)
Using your results from S5.2, consider an elliptical
perturbation of any stable circular orbit, and compute the radial
oscillation frequency
w_r in terms of the radius r. (Eliminate the explicit l
dependence by solving for l in terms of r). The precession rate is
w_p = w_phi - w_r,
where w_phi is the angular velocity. Assuming that r >> r_g,
find the leading order precession rate by expanding in the small
number r_g/r
(which is 2/r in our nice units). Restore the dimensionful
quantities using dimensional analysis, and evaluate the precession
rate for the
perihelion of Mercury, expressing the answer in seconds of arc per
century. For the radius use the semimajor axis of Mercury's
elliptical
orbit. [Tip: When
restoring the dimensionful quantities, you need only multiply by
appropriate powers of r_g/2 and c (which are both equal to 1 in
these units) so that you wind up with something that has the
dimensions of inverse time.]
(Answer: w_p =
3r^(-5/2), 43"/century.)
HW4 - due at
beginning of class, Tuesday 10/04/11
7.49 (charged particle in
uniform magnetic field - polar gauge)
Add parts: (d) Show that the Hamiltonian is the kinetic energy.
(e) Find the conserved quantity conjugate to the angle
about the z-axis.
This is the z-component of angular
momentum, which is not the same as "mvr" (much as the linear
momentum is not mv.)
(f) For orbits of constant rho in the
xy plane, show that the energy is proportional to the
z-component of angular momentum,
and find
the coefficient of proportionality. (g) Show that dimensional analysis gives
the same result for the energy, up to a dimensionless constant
(which happens to be -1, and you could infer the sign from the
sign of the angular momentum). (h) Assuming the angular momentum
is
quantized in integer multiples of hbar, and using your classical
orbit results, find the quantized energies of a charge in a
planar orbit
in a uniform magnetic field. These
are the "Landau levels", except for a missing zero point energy.
(The quantum ground state has
zero angular momentum, but has a zero point energy hbar omega/2,
where omega is the orbital angular frequency.)
7.51 (pendulum with
constrained Cartesian coordinates)
S4.1 (ambiguity of Lagrangian)
Show explicitly that if you add a
total time derivative (d/dt)f(q^i,t) to the Lagrangian
all the added terms in Lagrange's equations cancel.
Allow for an arbitrary number of generalized coordinates q^i, and
make use of the Einstein summation convention when appropriate.
S4.2 (charged particle in uniform magnetic field - Cartesian gauge)
A uniform magnetic field of strength B in the zhat direction is
described by a vector potential A = Bx yhat.
Unlike the potential in problem
7.49 for the same magnetic field, this one is not obviously
rotation invariant, but it is y-translation invariant. (a) Write
the Lagrangian for a particle of
mass m and charge q using this gauge, and find the equations of
motion. (b) From now on, choose units with m = qB = 1, and show
that the solutions
that have no velocity in the z-direction are uniform circular
orbits with any center, with the same angular velocity you found
in problem 7.49.
(c) Show that the conserved quantity associated with the
y-translation symmetry is qB times the x-component of the center
of the circular orbit (!).
(d) Show that the Hamiltonian is the kinetic energy. (e) Show
directly that the equations of motion imply that the kinetic
energy is conserved.
S4.3 (constraint for bead on
spinning rod) Consider problem 7.21 again. (a)
This time include the angular coordinate as a degree of freedom,
and write the Lagrangian with a constraint with Lagrange
multiplier that imposes the condition that the rod spins with
constant angular velocity
omega. (b) Write the radial and angular equations of motion. (c)
Show that the Lagrange multiplier term in the angular equation
of motion yields
the torque on the bead.
S4.4 (bead on a
circular ring) Consider a bead of mass m sliding on a
frictionless ring of radius R in a fixed vertical plane.
(a) Write the Lagrangian and equations of motion using polar
coordinates based at the center of the ring, with angle theta =
0 at the top, including
the constraint r = R with a Lagrange multiplier. (b)
Suppose the bead starts from rest nearly at the top. Use the
Lagrange multiplier to determine
(i) at what angle the force of constraint goes to zero (this is
where a bead sliding on the surface of a sphere would leave the
sphere), and
(ii) the magnitude and direction of the force of constraint when
the bead reaches the bottom of the ring.
S5.5 (catenary) A
flexible rope or chain hanging from two points forms a shape
called a "catenary", which is actually a hyperbolic cosine!
You can show this directly using the condition that the net
force on an infinitesimal segment of rope vanishes, but you can
also show it using
variational calculus. The calculation is nearly identical to the
problem of the soap film bounded by two rings that you already
solved on HW1,
but you need to impose the constraint that the rope has a given
fixed length b. Suppose the rope hangs from two points (x1,
y1) and (x2,
y2),
and
express the gravitational potential energy of of rope as an
integral involving the function x(y). Using the method of
Lagrange multipliers,
find the equation on x(y) implied by the fact that the potential
is a minimum for all variations of x(y) that keep the rope
length fixed.
Then re-expressed this as an equation on y(x), and show that the
solutions are given by y = c + y0 cosh[(x-x0)/y0] (this should
be easy).
Compared to the soap film problem, there is an extra constant
here, c. This is necessary since there is one more condition to
be met:
besides the positions of the endpoints, the length of the rope
must be matched.
HW3 - due at
beginning of class, Tuesday 9/27/11
7.20 (particle on a helix)
(Note: this could equally well be solved by writing down the
total mechanical energy and setting its time derivative to
zero.)
7.21 (bead on a spinning rod)
Let's add two interesting parts (b): Suppose the rod has length
L from the pivot to one end, and the bead is released with zero
radial velocity at some initial radius r_0. In the limit that
r_0 goes to zero, what velocity does the bead have when it
reaches the end of the rod? Give the
components of the velocity along the rod and perpendicular to
the rod. Compare your result to what you would expect from
dimensional analysis.
(c) (i) Explain why the Hamiltonian is a conserved quantity, and
show that it is equal to E - Ω J,
where E and J are the kinetic energy and angular momentum
of the bead, and Ω is the angular velocity of the
rod (called lower case omega in the book). (ii) Use this
conservation law to evaluate the kinetic energy when
the bead leaves the rod, and show that it agrees with what you get
using your results from part (b).
7.22 (pendulum in an
accelerating elevator) Change the second part: instead
of finding the equation of motion, show that the Lagrangian
itself is equal,
up to addition of a total time derivative, to the Lagrangian for
a non-accelerating pendulum with g replaced by g+a. (The total
time derivative does not
affect the equations of motion.) (Hint: You'll need to do an integration by
parts on t d(cos theta)/dt.)
7.37 (two masses connected by
a string through a hole in a table) Add part (e) What
does dimensional analysis (and a bit of physics logic) tell you
about the oscillation frequency in part (d)? (The problem
involves m, L (length of the string), g, and r_0.)
S3.1 (double pendulum)
(a) Write the Lagrangian for the double pendulum illustrated in
Fig. 7.3, assuming the masses are equal. (b) Expand the
Lagrangian to quadratic order in the amplitude of the angles,
dropping the higher order terms. (This describes small
oscillations about the equilibrium.
Later we'll find the normal modes.)
S3.2 (conserved quantity
for bead on a spinning hoop) Consider Example
7.6, but call the angular velocity of the hoop Ω. (a)
Explain why neither the
total mechanical energy E nor the angular momentum J about the
vertical axis are conserved. (Hint:
What work and torque are exerted on the bead?)
(b) Show that the Hamiltonian is equal to E - Ω J. (c)
Explain why the system (bead on a spinning hoop)
is neither symmetric under time translation nor under rotation.
[N.B. I refer here to the system,
not to the Lagrangian
(d) On general grounds, the
Hamiltonian is conserved, since the Lagrangian has no explicit
t-dependence.
What is the symmetry whose corresponding conserved
quantity is the Hamiltonian, E - Ω J?
S3.3 (least action for
free-fall) Consider a particle of mass m in a uniform
gravitational potential U(y) = mgy. Let y represent the height
of the particle,
and consider the action for vertical paths y(t) over the time
interval [-T,T], with y(-T) = 0 = y(T). The action for the path
y(t)=0 is zero, of course.
(a) Show that the paths with uniform velocity on the way up and
the way down have minimum action if they reach the height h =
1/2 g T^2, and
show that the value of that minimum action is -1/4 m g^2 T^3.
(b) Show that the classical path reaches same same height as you
found in part (a),
and that the classical action on that path is -1/3 m g^2 T^3.
(c) Show that the m g^2 T^3 factor
follows from dimensional analysis.
HW2 - due at
beginning of class, Tuesday 9/20/11
6.6 (ds in various
settings)
6.16 (geodesics on a sphere)
6.19 (soap bubble)
Note: Assume y1 and y2
are both positive.
7.3 (2d oscillator using
Lagrangian) Modify this problem as follows: (a) Write
the Lagrangian using polar coordinates.
(b) Find the equations of motion for phi and r. (c) define the
angular momentum by J = mr2phidot, and use this to
eliminate
phidot from the equation for r. (d) Find the circular orbit
(motion with rdot=0) for a given J. (e) Describe the radial
motion
for orbits with nonzero rdot and nonzero J.
7.8 (CM
and relative motion of two bodies in 1d) Modify
this problem as follows: Instead of a spring, suppose the
potential
U is an arbitrary function of the separation of the two
particles, U(x), with x = x1 - x2. Replace
part (c) by "describe the resulting
dynamics". Add part (d) Now let the two masses be different.
Find the Lagrangian in terms of the center of mass position
X = (m1x1
+ m2x2)/M,
where
M
is the total mass M = m1 + m2, and x, and
again describe the resulting dynamics.
(The reduced mass µ = (m1m2)/M will
arise. Express your result using this quantity.)
S2 (motion about equilibria
for bead on spinning hoop)
Example 7.7 discusses motion near equilibria of a bead on a
spinning hoop. It does so by expanding the equation of motion
about the equilibrium. Instead, let's expand the Lagrangian
about the equilibrium. Hence: the Lagrangian (7.68) can be
written
in the form L = 1/2 m R^2 qdot^2 - U(q), where I'm using q to
stand for theta. At an equilibrium point q_0 we have U'(q_0)=0.
(a) Make a Taylor expansion of U(q) about an equilibrium point
q_0. Keep terms out to quadratic order in the displacement from
equilbrium. (b) Stability is determined by the sign of the
quadratic term in U(q): if it's positive the motion is a stable
oscillation, if
it's negative the motion runs away exponentially. (i) Show that
below the critical angular velocity, the equilibrium at theta=0
is
stable and that at theta=π is unstable. (ii) Show that above the
critical angular velocity the equilibria at theta = 0 and π are
both unstable,
while the two new equilibrium points are stable.
HW1 - due at
beginning of class, Tuesday 9/13/11
Read the Chapter summaries of Chapters 1-5, and browse the
chapters. If you see anything you're not comfortable
with read that part of the chapter.
1.32 - (Newton's 3rd law
failure for magnetic forces between point charges)
read, but do not do, this problem
1.33 (Newton's 3rd law for
magnetic forces between closed current loops) This problem
is intriguing.
A couple of hints: (i) Evaluate grad1(1/s),
and
use
that
to
show
that
a part of the integral vanishes
because it's the line integral of a gradient around a closed loop. (ii) To establish the 3rd law,
show that
the force of loop 2 on loop 1 is proportional to ∫∫ (dr1. dr2) s/s3, which
obviously changes sign when
considering the force of loop 1 on loop 2.
4.4 (Particle on table
attached to string through hole)
4.24 (Gravitational potential
of an infinite rod) Do only parts (a) and (d). [Once
you find a potential that
gives the force, you don't need to check that the curl of the
force vanishes, since curl grad U = 0 for any U.]
4.36 (ball-pulley-mass gizmo)
4.38 (pendulum period with
large amplitude) For part (b), do not use the "complete
elliptic integral" function.
Instead, just use whatever software you like to evaluate
numerically, by brute force, to sufficient accuracy, the
integral, and make the plot in question.
4.39 (pendulum period,
amplitude expansion) Add part (d): Plot the analytic
result you get for the approximate
period on the same graph as the numerical result from 4.38.
4.41 (virial theorem)
4.43 (central, spherically
symmetric forces) Add part (c) Better yet, show
directly by construction that one can
find a U such that F = -
grad U.