Non Relativistic String
Non Relativistic String
Non-Relativistic Strings
87
88 CHAPTER 4. NON-RELATIVISTIC STRINGS
Figure 4.2: Left side: string with Dirichlet boundary conditions at the endpoints.
Right side: string with Neumann boundary conditions at the endpoints.
Let’ s see how we can solve the wave equation for a particular set of initial
conditions. The general solution of equation (4.1.6) is of the form
y(t, x) = h+ (x − v0 t) + h− (x + v0 t) , (4.2.3)
d2 y(x) µ0
2
+ ω2 y(x) = 0 . (4.3.2)
dx T0
This is an ordinary second-order differential equation for the profile y(x)
of the oscillations. The allowed frequencies are selected by this equation,
together with the boundary conditions. Since ω, µ0 and T0 are constants,
92 CHAPTER 4. NON-RELATIVISTIC STRINGS
These are the frequencies of oscillation for a Dirichlet string. The strings on
a violin are Dirichlet strings. To tune a violin to the right frequency one must
adjust the string tension. The higher the tension, the higher is the pitch, as
predicted by (4.3.4).
For the case of Neumann boundary conditions (4.2.2), we obtain the
spatial solutions
nπx
yn (x) = An cos n = 1, 2, . . . . (4.3.5)
a
This time the n = 0 case is a little less trivial: the string does not oscillate,
but it becomes rigidly translated to y(t, x) = A0 . The oscillation frequencies,
found by plugging (4.3.5) in (4.3.2), are the same as those in (4.3.4). There-
fore, the oscillation frequencies are the same in the Neumann and Dirichlet
problems. The Neumann case, however, admits one extra solution not in-
cluded in our oscillatory ansatz (4.3.1): the string can translate with constant
velocity. Indeed, y(t, x) = at + b, with a and b arbitrary constants, satisfies
both the boundary conditions and the original wave equation (4.1.7).
string that can be chosen to be sufficiently small so that the mass density is
approximately constant. We therefore get:
∂ 2 y µ(x) ∂ 2 y
− = 0. (4.4.1)
∂x2 T0 ∂t2
For normal oscillations, we use the ansatz in (4.3.1) and find
d2 y µ(x) 2
+ ω y(x) = 0 . (4.4.2)
dx2 T0
This equation is no longer simple to solve, and it can only be studied in detail
once the function µ(x) is specified. In Problems 4.4 and 4.5 you will consider
some specific mass distributions, and you will explore a variational approach
that gives an upper bound for the lowest oscillation frequency.
L=T −V , (4.5.1)
94 CHAPTER 4. NON-RELATIVISTIC STRINGS
where T is the kinetic energy of the system and V is the potential energy
of the system. For a point particle of mass m moving along the x axis
under the influence of a time-independent potential V (x), the non-relativistic
Lagrangian takes the form
1 dx(t)
L(t) = m [ ẋ(t)]2 − V [x(t)] , ẋ(t) ≡ . (4.5.2)
2 dt
We must emphasize that the above Lagrangian is implicitly a function of
time, but it has no explicit time dependence. All the time dependence arises
from the time dependence of the position x(t). The action S is defined as
S= L(t)dt , (4.5.3)
P
application, the input for the action functional is the function x(t) which
determines the path P. We can emphasize the argument of S by using the
notation S[x]. Here [x] represents the full function x(t). It is potentially
confusing to write S[x(t)], since it suggests that S is ultimately a function of
t, which it is not.
More explicitly, for any path x(t), the action is given by
tf
1 2
S[x] = m [ẋ(t)] − V [x(t)] dt . (4.5.4)
ti 2
It is very important to emphasize that the action S can be calculated for any
path x(t) and not only for paths that represent physically-realized motion.
It is because S can be calculated for all paths that it is a very powerful tool
to find the paths that can be physically realized.
Figure 4.4: A path x(t) and its variation x(t)+δx(t). This variation δx(t) vanishes
at t = ti and at t = tf .
Hamilton’s principle states that the path P which a system actually takes
is one for which the action S is stationary. More precisely, if this path P
is varied infinitesimally, the action does not change to first order in the
variation. In terms of the function x(t) specifying the path, the perturbed
path takes the form x(t) + δx(t), as shown in Figure 4.4. For any time t,
the variation δx(t) is the vertical distance between the original path and the
varied path. As in the figure, we consider variations where the initial and
final positions xi = x(ti ) and xf = x(tf ) are unchanged:
δx(ti ) = δx(tf ) = 0 . (4.5.5)
96 CHAPTER 4. NON-RELATIVISTIC STRINGS
We now calculate the action S[x + δx] for the perturbed path x(t) + δx(t):
tf
2
m d
S[x + δx] = (x(t) + δx(t)) − V [x(t) + δx(t)] dt (4.5.6)
ti 2 dt
tf
d
= S[x] + mẋ(t) δx(t) − V [x(t)]δx(t) dt + O((δx)2 ) ,
ti dt
where in the last equation we have expanded V in a Taylor series about x(t).
The terms of order (δx)2 and higher are unnecessary to determine whether
the action is stationary. We have thus left them undetermined and only
indicate them by O((δx)2 ). We can write the new action as S + δS, where
δS is linear in δx. From the equation above we see that δS is given by
tf
d
δS = mẋ(t) δx(t) − V [x(t)]δx(t) dt . (4.5.7)
ti dt
To find the equations of motion, the variation δS must be rewritten in the
form δS ∼ δx{. . .}. In particular, no derivatives must be acting on δx.
This can be achieved using integration by parts:
tf
d
δS = [mẋ(t)δx(t)] − mẍ(t)δx(t) − V [x(t)]δx(t) dt
ti dt
tf
= mẋ(tf )δx(tf ) − mẋ(ti )δx(ti ) + δx(t) {−mẍ(t) − V [x(t)]} dt .
ti
(4.5.8)
The action is stationary if δS vanishes for every variation δx(t). For this to
happen, the factor multiplying δx(t) in the integrand must vanish:
Suppose we have determined the path that the particle takes while going
from xi to xf . As we have seen, the action is then stationary under variations
that vanish at the initial and final times. Is the action also stationary under
variations that change the initial position at ti or the final position at tf ?
In general, the answer is no. This can be seen from equation (4.5.8). The
integral term vanishes by assumption, but if δx(tf ) = 0, the first term in the
right-hand side would not vanish unless mẋ(tf ), the final momentum of the
particle, happens to vanish. The situation is analogous for δx(ti ) = 0.
L(t) = µ0 − T0 dx ≡ L dx , (4.6.4)
0 2 ∂t 2 ∂x 0
98 CHAPTER 4. NON-RELATIVISTIC STRINGS
In this action the “path” is the function y(t, x). To find the equations of
motion, we must examine the variation of the action as we vary: y(t, x) →
y(t, x) + δy(t, x). Performing the variation as before, we get
tf a
∂y ∂(δy) ∂y ∂(δy)
δS = dt dx µ0 − T0 . (4.6.7)
ti 0 ∂t ∂t ∂x ∂x
−T0 δy + T0 2 δy . (4.6.8)
∂x ∂x ∂x
The total derivatives can be integrated. The total time derivative on the first
line reduces to evaluations at tf and ti , while the total space derivative on
the second line gives evaluations at the string endpoints:
a t=tf tf x=a
∂y ∂y
δS = µ0 δy dx − T0 δy dt
0 ∂t t=ti ti ∂x x=0
tf a
∂2y ∂2y
− dt dx µ0 2 − T0 2 δy . (4.6.9)
ti 0 ∂t ∂ x
Our final expression for δS contains three terms. Each one must vanish inde-
pendently. The third term, for example, is determined by the motion of the
string for x ∈ (0, a) and t ∈ (ti , tf ). The boundary conditions do not restrict
δy(t, x) here, so we set to zero the coefficient of δy, and recover our original
4.6. THE NON-RELATIVISTIC STRING LAGRANGIAN 99
∂y ∂y
(t, x = 0) = (t, x = a) = 0 , Neumann B.C. (4.6.10)
∂x ∂x
This will also cause the second term to vanish. Dirichlet boundary conditions
can be written in a form where the similarity to Neumann boundary condi-
tions is more apparent. If the string endpoints are fixed, the time derivatives
of the endpoint coordinates must vanish
∂y ∂y
(t, x = 0) = (t, x = a) = 0 , Dirichlet B.C. (4.6.11)
∂t ∂t
The similarity with (4.6.10) is quite striking. The only change is that spatial
derivatives were turned into time derivatives. If we write Dirichlet boundary
conditions in this form, we must still specify the values of the coordinates of
the fixed endpoints.
δS = dt dx δ ẏ + δy = dt dx P δ ẏ +P δy . (4.6.17)
t x
ti 0 ∂ ẏ ∂y ti 0
4.6. THE NON-RELATIVISTIC STRING LAGRANGIAN 101
∂ t ∂
δS = dt dx (P δy) + (P δy) − δy
x
+ . (4.6.18)
ti 0 ∂t ∂x ∂t ∂x
∂P t ∂P x
+ = 0, (4.6.19)
∂t ∂x
Using (4.6.16), one can see that this is the wave equation (4.1.6).
Note that P t , as given in (4.6.16), coincides with the momentum density
used before in equation (4.6.12). This is not a coincidence. In Lagrangian
mechanics, the derivative of the Lagrangian with respect to a velocity is
the conjugate momentum. For the string, ẏ plays the role of a velocity, and
therefore P t , a derivative of the Lagrangian density with respect to a velocity,
is a momentum density.
In addition, note that for free string endpoints, the vanishing of the varia-
tion δS requires that P x = 0. As we can see from (4.6.16), this is a Neumann
boundary condition. Furthermore, P t vanishes at the string endpoints when
impose a Dirichlet boundary condition ((4.6.11)). A more detailed analysis
of these ideas will be given in Chapter 8, where P t and P x will be shown to
have an interesting two-dimensional interpretation.
102 CHAPTER 4. NON-RELATIVISTIC STRINGS
Problems
Problem 4.1. Consistency of small transverse oscillations.
(a) Show that the above Dirichlet boundary conditions imply that in the
notation of equation (4.2.3)
Now consider an initial value problem for this string. At t = 0 the transverse
displacement is identically zero, and the velocity is
∂y x x
(0, x) = v0 1− , x ∈ (0, a) . (2)
∂t a a
(b) Calculate h+ (u) for u ∈ (−a, a). Does this define h+ (u) for all x?
Problems for Chapter 4 103
(d) At t = 0 the midpoint x = a/2 has the largest velocity of all points in
the string. Show that the velocity of the midpoint reaches the value
of zero at time t0 = a/(2v0 ) and that y(t0 , a/2) = a/12. This is the
maximum vertical displacement of the string.
(a) Set up a variational procedure that gives an upper bound on the lowest
frequency of oscillation ω0 . (This can be done as in quantum mechanics,
where the ground state energy E0 of a system with Hamiltonian H
satisfies E0 ≤ ψ,Hψ
ψ,ψ
). A useful first step may be to define an inner
product · , · such that ψi , ψj vanishes when ωi = ωj . Explain why
your variational procedure works.
(b) Consider the case µ(x) = µ0 xa . Use your variational principle to find
a simple bound on the lowest oscillation frequency. Compare with the
answer ω02 = (18.956) µT0 a0 2 obtained by a direct numerical solution of
the eigenvalue problem.