Home / Math and Physics Files / Physics / Physics Book Downloads / Math Methods in Physics Books / PDF Originals
Stone M. Methods of Mathematical Physics I (2002)(316s)-1
PDF · 316 pages · 1.8 MB
Open PDF file
Lecture notes by Michael Stone (University of Illinois) for a one-semester graduate mathematical methods course, stressing linear operators on function spaces as analogues of matrices. Chapters cover calculus of variations, function spaces, linear ODEs and differential operators, Green functions, PDEs, real waves and solitons, special functions, and integral equations, with a linear algebra appendix. This is a downloaded book in Phil's collection, not his own work.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
Metho dsofMathematical PhysicsI
Asetoflecture notes by
MichaelStone
PIMANDER-CASA UBON
AlexandriaFlorenceLondon
ii
Copyrightc
2001,2002 M.Stone.
Allrightsreserv ed.Nopartofthismaterial canbereproduced, stored or
transmitted without thewritten permission oftheauthor. Forinformation
contact:MichaelStone, LoomisLaboratory ofPhysics, UniversityofIllinois,
1110WestGreen Street, Urbana, IL61801, USA.
Preface
These notes wereprepared forPHYCS-498MMA, afairly traditional one-
semester mathematical metho dscourse forbegining graduate studen tsin
physics. Theemphasis isonlinear operators andstresses theanalogy between
suchoperators acting onfunction spaces andmatrices acting onnitedimen-
sional spaces. Theoperator language thenprovides aunied framew orkfor
investigating ordinary dieren tialequations, partial dieren tialequations,
andintegral equations.
Although thismathematics isapplicable toawiderange physical phenom-
ena,theillustrativ eexamples aremostly drawnfromclassical andquantum
mechanics. Classical mechanics isasubjectfamiliar toallphysicsstuden ts
andthepointbeingillustrated isimmediately understandable without any
further specialized knowledge. Similarly allphysics studen tshavestudied
quantummechanics, andherethematrix/dieren tial-op erator analogy lies
attheheart ofthesubject.
Themathematical prerequisites forthecourse areasound grasp ofun-
dergraduate calculus (including thevectorcalculus needed forelectricit yand
magnetism courses), linear algebra (themore thebetter), andcompetence
atcomplex arithmetic. Fourier sums andintegrals, aswellasbasic ordinary
dieren tialequation theory receiveaquickreview, butitwouldhelpifthe
reader hadsome prior experience tobuild on.Contourintegration isnot
required.
iii
iv PREF ACE
Contents
Preface iii
1Calculus ofVariations 1
1.1What isitgoodfor? .......................1
1.2Functionals ............................2
1.2.1 TheFunctional Derivative................2
1.2.2 Examples .........................3
1.2.3 FirstIntegral .......................8
1.3Lagrangian Mechanics ......................9
1.3.1 OneDegree ofFreedom ..................10
1.3.2 Noether's Theorem ....................14
1.3.3 ManyDegrees ofFreedom ................17
1.3.4 ContinuousSystems ...................17
1.4Variable EndPoints........................26
1.5Lagrange Multipliers .......................33
2Function Spaces 39
2.1Motiv ation .............................39
2.1.1 Functions asVectors ...................40
2.2Norms andInner Products ....................41
2.2.1 Norms andConvergence .................41
2.2.2 Norms fromIntegrals ...................43
2.2.3 HilbertSpace .......................45
2.2.4 Orthogonal Polynomials .................51
2.3Linear Operators andDistributions ...............55
2.3.1 Linear Operators .....................55
2.3.2 Distributions .......................57
2.4Fourier Series andIntegrals. ...................61
v
vi CONTENTS
2.4.1 Fourier Series .......................62
2.4.2 Fourier Integral Transforms ...............64
2.4.3 ThePoisson Summation Formula............66
3Linear Ordinary Dieren tialEquations 69
3.1Existence andUniqueness ofSolutions .............69
3.1.1 FlowsforFirst-Order Equations .............69
3.1.2 Linear Indep endence ...................71
3.1.3 TheWronskian ......................72
3.2Normal Form ...........................76
3.3Inhomogeneous Equations ....................77
3.3.1 Particular Integral andComplemen taryFunction ...77
3.3.2 Variation ofParameters .................78
3.4Singular Points..........................80
4Linear Dieren tialOperators 83
4.1Formal vs.Concrete Operators .................83
4.1.1 TheAlgebra ofFormal Operators ............83
4.1.2 Concrete Operators ....................85
4.2TheAdjoin tOperator ......................86
4.2.1 TheFormal Adjoin t....................86
4.2.2 ASimple EigenvalueProblem ..............90
4.2.3 Adjoin tBoundary Conditions ..............92
4.2.4 Self-adjoin tBoundary Conditions ............93
4.3Completeness ofEigenfunctions .................99
4.3.1 Discrete Spectrum ....................99
4.3.2 Continuousspectrum ...................104
5Green Functions 115
5.1Inhomogeneous Linear equations .................115
5.1.1 Fredholm Alternativ e...................115
5.2Constructing Green Functions ..................116
5.2.1 Sturm-Liouville equation .................117
5.2.2 Initial ValueProblems ..................119
5.2.3 Modied Green Functions ................124
5.3Applications ofLagrange's Identity...............126
5.3.1 Hermiticit yofGreen function ..............126
5.3.2 Inhomogeneous Boundary Conditions ..........127
CONTENTS vii
5.4Eigenfunction Expansions ....................129
5.5Analytic Properties ofGreen Functions .............130
5.5.1 Causalit yImplies Analyticit y..............130
5.5.2 Plemelj Formul.....................135
5.5.3 Resolv entOperator ....................137
5.6LocalityandtheGelfand-Dikii equation ............142
6Partial Dieren tialEquations 145
6.1Classication ofPDE's ......................145
6.1.1 CauchyData .......................147
6.1.2 Characteristics andrst-order equations ........149
6.2WaveEquation ..........................150
6.2.1 d'Alem bert'sSolution ...................150
6.2.2 Fourier's Solution .....................153
6.2.3 Causal Green Function ..................153
6.2.4 Oddvs.EvenDimensions ................158
6.3HeatEquation ...........................163
6.3.1 HeatKernel ........................164
6.3.2 Causal Green Function ..................165
6.3.3 Duhamel's Principle ...................167
6.4Laplace's Equation ........................169
6.4.1 Separation ofVariables ..................169
6.4.2 Green Functions .....................175
6.4.3 Metho dofImages .....................177
6.4.4 Kirchhovs.Huygens ..................179
7TheMathematics ofRealWaves 183
7.1Dispersivewaves.........................183
7.1.1 Ocean Waves.......................183
7.1.2 Group Velocity......................187
7.1.3 Wakes...........................190
7.1.4 Hamilton's Theory ofRays................193
7.2Making Waves...........................195
7.2.1 Rayleigh's Equation ...................195
7.3Non-linear Waves.........................199
7.3.1 Sound inAir.......................200
7.3.2 Shocks...........................202
7.3.3 WeakSolutions ......................208
viii CONTENTS
7.4Solitons ..............................209
8SpecialFunctions I 215
8.1Curvilinear Co-ordinates .....................215
8.1.1 Div,Grad andCurlinCurvilinear Co-ordinates ....218
8.1.2 TheLaplacian inCurvilinear Co-ordinates .......221
8.2Spherical Harmonics .......................221
8.2.1 Legendre Polynomials ..................222
8.2.2 Spherical Harmonics ...................227
8.3Bessel Functions .........................230
8.3.1 Cylindrical Bessel Functions ...............230
8.3.2 Orthogonalit yandCompleteness ............237
8.3.3 Modied Bessel Functions ................240
8.3.4 Spherical Bessel Functions ................243
8.4Singular Endpoints........................247
8.4.1 Weyl'sTheorem ......................247
9Integral Equations 255
9.1Illustrations ............................255
9.2Classication ofIntegral Equations ...............256
9.3Integral Transforms ........................257
9.3.1 Fourier Metho ds.....................258
9.3.2 Laplace Transform Metho ds...............260
9.4Separable Kernels .........................263
9.4.1 Eigenvalueproblem ....................263
9.4.2 Inhomogeneous problem .................264
9.5Singular Integral Equations ...................266
9.5.1 Solution viaTchebychefPolynomials ..........266
9.6Some Functional Analysis ....................269
9.6.1 Bounded andCompact Operators ............269
9.6.2 Closed Operators .....................272
9.7Series Solutions ..........................276
9.7.1 Neumann Series ......................276
9.7.2 Fredholm Series ......................276
AElemen taryLinear Algebra 281
A.1Vector Space ...........................281
A.1.1 Axioms ..........................281
CONTENTS ix
A.1.2 Bases andComp onents..................282
A.2Linear Maps ............................283
A.2.1 Range-Nullspace Theorem ................284
A.2.2 TheDualSpace ......................284
A.3Inner-Pro ductSpaces .......................286
A.3.1 Inner Products ......................286
A.3.2 Adjoin tOperators ....................288
A.4Inhomogeneous Linear Equations ................289
A.4.1 Fredholm Alternativ e...................291
A.5Determinan ts...........................292
A.5.1 Skew-symmetric n-linear Forms .............292
A.5.2 TheAdjugate Matrix ...................294
A.5.3 Dieren tiating Determinan ts...............296
A.6Diagonalization andCanonical Forms ..............297
A.6.1 Diagonalizing Linear Maps ................297
A.6.2 Quadratic Forms .....................301
A.6.3 Symplectic Forms .....................303
x CONTENTS
Chapter 1
Calculus ofVariations
Inthischapter wewillstudy whatiscalled thecalculus ofvariations .Many
physicsproblems canbeformulated inthelanguage ofthiscalculus, andonce
theyarethere areuseful toolstohand. Inthetextandassociated exercises
wewillmeetsomeoftheequations whose solution willoccupyusfortherest
ofthecourse.
1.1What isitgoodfor?
Theclassical problems ofthecalculus ofvariations include:
i)Dido's problem :InVirgil's Aeneid,Queen DidoofCarthage needs to
ndlargest areathatcanbeenclosed byacurve(astripofbull's hide)
ofxedlength.
ii)Plateau's problem :Findthesurface ofminim umareaforagivensetof
bounding curves.Asoaplmonawireframe willadopt thisminimal-
areaconguration.
iii)Johann Bernoulli's Brachistochrone:Abeadslides downacurvewith
xedends. Assuming thatthetotalenergy1
2mv2+V(x)isconstan t,
ndthecurvethatgivesthemostrapid descen t.
iv)Catenary :Findtheformofahanging heavychainofxedlength by
minimizing itspotentialenergy .
Allthese problems involvending maxima orminima, andhence equating
somesortofderivativetozero.Inthenextsection wewilldene thisderiva-
tive,andshowhowtocompute it.
1
2 CHAPTER 1.CALCULUS OFVARIA TIONS
1.2Functionals
Invariational problems weareprovided withanexpressionJ[y]that\eats"
whole functionsy(x)andreturns asingle number.Suchobjectsareoften
called functionals todistinguish them fromordinary functions. Anordinary
function isamapf:R!R.Afunctional, J,isamapJ:C1(R)!R
whereC1(R)isthespace ofsmooth(havingderivativesofallorders) func-
tions. Tondthefunctiony(x)thatmaximizes orminimizes agivenfunc-
tionalJ[y]weneedtodene, andevaluate, itsfunctional derivative.
1.2.1 TheFunctional Derivative
Wewillrestrict ourselv estoexpressions oftheform
J[y]=Zx2
x1f(x;y;y0;y00;y(n))dx; (1.1)
depending onthevalueofy(x)andonlynitely manyofitsderivatives.Such
functionals aresaidtobelocalinx.
Consider rstafunctional depending onlyonx,yandy0.Wevaryy(x)!
y(x)+(x)whereisanx-indep enden tconstan t,andwrite
J[y+] J[y]=Zx2
x1ff(x;y+;y0+0) f(x;y;y0)gdx
=Zx2
x1(
@f
@y+d
dx@f
@y0+O(2))
dx
="
@f
@y0#x2
x1+Zx2
x1((x))(@f
@y d
dx @f
@y0!)
dx+O(2):
Forthemomen tletusassume that(x1)=(x2)=0.Thatis,weareusing
\xed endpoint"variations. Inthiscasetheintegrated-out partvanishes,
and
J=Zx2
x1((x))(@f
@y d
dx @f
@y0!)
dx
=Zx2
x1y(x) J
y(x)!
dx: (1.2)
Herey(x)(x),andthequantity
J
y(x)@f
@y d
dx @f
@y0!
(1.3)
1.2.FUNCTIONALS 3
iscalled thefunctional (orFrechet)derivativeofJwithrespecttoy(x).We
canthink ofitasakindofgeneralization ofthenotion ofapartial derivative
@J=@yi,withthediscrete subscript \i"onybeingreplaced byacontinuous
label,\x".Thus
J=X
i@J
@yiyi!Zx2
x1dx J
y(x)!
y(x): (1.4)
Thecondition forthefunctional tobestationary under variationsy!
y+yis
J
y(x)=@f
@y d
dx @f
@y0!
=0; (1.5)
andthisisusually called theEuler-L agrangeequation.
Ifthefunctional dependsonmorethanonefunctiony,thenstationarit y
under allpossible variations requires oneequation
J
yi(x)=@f
@yi d
dx @f
@y0
i!
=0 (1.6)
foreachfunctionyi(x).
Ifthefunction dependsonhigher derivatives,y00,y(3),etc.,thenwehave
tointegrate byparts more times, andweendupwith
J
y(x)=@f
@y d
dx @f
@y0!
+d2
dx2 @f
@y00!
d3
dx3 @f
@y(3)!
+:(1.7)
1.2.2 Examples
Nowweapply ournewderivativetosolvesome simple problems.
Soaplmsupported byapairofcoaxial rings.
xy(x)
1 x x2
4 CHAPTER 1.CALCULUS OFVARIA TIONS
Herewewishtominimize thefreeenergy ofthelm,whichisequal totwice
(once foreachliquid-air interface) thesurface tensionofthesoapsolution
times theareaofthelm.Wetherefore needtominimize
J[y]=4Zx2
x1yq
1+y02dx: (1.8)
withy(x1)=y1andy(x2)=y2.Weformthepartial derivatives
@f
@y=4q
1+y02;@f
@y0=4yy0
q
1+y02(1.9)
andthuswrite downtheEuler-Lagrange equation
q
1+y02 d
dx0
@yy0
q
1+y021
A=0: (1.10)
Performing theindicated derivativewithrespecttoxgives
q
1+y02 (y0)2
q
1+y02 yy00
q
1+y02+y(y0)2y00
(1+y02)3=2=0: (1.11)
Collecting terms, thisis
1q
1+y02 yy00
(1+y02)3=2=0: (1.12)
Thisdieren tialequation looksatri
eintimidating. Tosimplify ,wemultiply
byy0toget
0=y0
q
1+y02 yy0y00
(1+y02)3=2
=d
dx0
@yq
1+y021
A: (1.13)
Thesolution totheminimization problem therefore reduces tosolving
yq
1+y02=; (1.14)
1.2.FUNCTIONALS 5
whereisanasyetundetermined integration constan t.Fortunately this
non-linear, rstorder, dieren tialequation iselemen tary.Wewrite itas
dy
dx=s
y2
2 1 (1.15)
andseparate variablesZ
dx=Zdyq
y2
2 1: (1.16)
Wenowmakethenatural substitution y=cosht,whence
Z
dx=Z
dt: (1.17)
Thuswendthatx+a=t,leading to
y=coshx+a
: (1.18)
Weselectandatottheendpointsy(x1)=y1andy(x2)=y2.
HeavyChain overPulleys. Wecannot yetconsider theformofahanging
chainofxedlength, butwecansolveasimpler problem ofaheavycable
drapedoverapairofpulleys located atx=L,y=h,andwiththeexcess
cable resting onahorizon talsurface.
h
−L +L
Thepotentialenergy ofthesystem is
P:E:=X
mgy=gZL
Lyq
1+(y0)2dx+const. (1.19)
6 CHAPTER 1.CALCULUS OFVARIA TIONS
Heretheconstan trefers totheunchanging potentialenergy ofthevertically
hanging cableandthecableonthehorizon talsurface. Notice thatthetension
inthecableisbeingtacitly determined bytheweightofthevertical segmen ts.
TheEuler-Lagrange equations coincide withthose ofthesoaplm,so
y=cosh(x+a)
(1.20)
where wehavetondanda.Wehave
h=cosh( L+a)=;
=cosh(L+a)=; (1.21)
soa=0andh=coshL=.Settingt=L=thisreduces to
h
L!
t=cosht: (1.22)
Byconsidering theintersection oftheliney=ht=Lwithy=coshtwesee
thatifh=Listoosmall there isnosolution (theweightofthesuspended
cable istoobigforthetension supplied bythedangling ends) andonceh=L
islargeenough there willbetwopossible solutions.
y
y= ht/Ly=cosht
t=L/κ
Intersection ofy=ht=Lwithy=cosht.
Further investigation willshowthatonlyoneofthese isstable.
Example: TheBrachistochrone.Thisproblem wasposedasachallenge by
Johann Bernoulli in1696. Heaskedwhatshapeshould awirewithendpoints
1.2.FUNCTIONALS 7
(0;0)and(a;b)takeinorder thatafrictionless beadwillslidefromrestdown
thewireintheshortest possible time(o&:shortest,oo&:time).
x
yg
(a,b)
When presen tedwithanostensibly anonymous solution, Johann made his
famous remark: Tanquam exunguem leonem1,|meaning thatherecognized
thattheauthor wasIsaac Newton.
Johann gaveasolution himself, butthatofhisbrother Jacob Bernoulli
wassuperiorandJohann triedtopassitoashis.Thiswasnotatypical.
Johann latermisrepresen tedthepublication dateofhisbookonhydraulics
tomakeitseem thathehadpriorit yinthiseldoverhisownson,Daniel
Bernoulli.
Webeginoursolution oftheproblem byobserving thatthetotalenergy
E=1
2m(_x2+_y2) mgy=1
2m_x2(1+y02) mgy; (1.23)
ofthebeadwillbeconstan t.Fromtheinitial condition weseethatthis
constan tiszero. Wetherefore wishtominimize
T=ZT
0dt=Za
01
_xdx=Za
0s
1+y02
2gydx (1.24)
soasndy(x),giventhaty(0)=0andy(a)=b.TheEuler-Lagrange
equation is
yy00+1
2(1+y02)=0: (1.25)
Again thislooksintimidating, butwecanusethesame trickofmultiplying
through byy0toget
y0
yy00+1
2(1+y02)
=1
2d
dxn
y(1+y02)o
=0: (1.26)
1Irecognize thelionbyhisclawmark.
8 CHAPTER 1.CALCULUS OFVARIA TIONS
Thus
2c=y(1+y02): (1.27)
Thishasaparametric solution
x=c( sin);
y=c(1 cos); (1.28)
(asyoushould verify) andthesolution isacycloid.
x
y(0,0)
(a,b)θ
θ(x,y)
Awheel rollsonthexaxis.Thedot,whichisxedtotherimofthewheel,
traces outacycloid.
Theparametercisdetermined byrequiring thatthecurvedoesinfactpass
through thepoint(a;b).
1.2.3 FirstIntegral
Howdidweknowthatwecould simplify boththesoap-lm problem and
thebrachistochrone bymultiplying theEuler equation byy0?Theanswer
isthatthere isageneral principle, closely related toenergy conserv ation in
mechanics, thattellsuswhen andhowwecanmakesuchasimplication. It
workswhen thefisoftheformf(y;y0),i.e.hasnoexplicit dependence on
x.Inthiscasethelasttermin
df
dx=y0@f
@y+y00@f
@y0+@f
@x(1.29)
isabsent,andwehave
d
dx
f y0@f
@y0!
=y0@f
@y+y00@f
@y0 y00@f
@y0 y0d
dx @f
@y0!
1.3.LAGRANGIAN MECHANICS 9
=y0 @f
@y d
dx @f
@y0!!
; (1.30)
andthisiszeroiftheEuler-Lagrange equation issatised. Thequantity
I=f y0@f
@y0(1.31)
isthusarstintegraloftheEuler-Lagrange equation. Inthesoap-lm case
f y0@f
@y0=yq
1+(y0)2 y(y0)2
q
1+(y0)2=yq
1+(y0)2: (1.32)
When there areanumberofdependen tvariableyi,sothatwehave
J[y1;y2;:::yn]=Z
dxf(y1;y2;:::yn;y0
1;y0
2;:::y0
n) (1.33)
thentherstintegral becomes
I=f X
iy0
i@f
@y0
i: (1.34)
Again
dI
dx=d
dx
f X
iy0@f
@y0
i!
=X
i
y0
i@f
@yi+y00
i@f
@y0
i y00
i@f
@y0
i y0
id
dx @f
@y0
i!!
=X
iy0
i @f
@yi d
dx @f
@y0
i!!
; (1.35)
andthiszeroiftheEuler-Lagrange equation issatised foreachyi.
Notethatthere isonlyonerstintegral, nomatter howmanyy'sthere
are.
1.3Lagrangian Mechanics
InhisMecanique Analytique (1788) Joseph-Louis deLaGrange, following
d'Alem bert(1742) andMaup ertuis (1744), showedthatmost ofclassical
10 CHAPTER 1.CALCULUS OFVARIA TIONS
mechanics canberecast asavariational principle: theprinciple ofleast
action .TheideaistointroducetheLagrangian functionL=T VwhereT
isthekinetic energy ofthesystem andVthepotentialenergy ,bothexpressed
interms ofgeneralizedcoordinatesqiandtheirtimederivatives_qi.Then
Lagrange showedthatthemultitude ofNewton's F=maequations, onefor
eachparticle inthesystem, could bereduced to
d
dt @L
@_qi!
@L
@qi=0; (1.36)
oneequation foreachgeneralized coordinateq.Quite remark ably|given
thatLagrange's derivation containsnomentionofmaxima orminima |we
observ ethatthisistheprecisely thecondition thattheaction integral
S=Ztfinal
tinitialL(qi;q0i)dt (1.37)
bestationary withrespecttovariations ofthetrajectoryqi(t)whichleavethe
initial andnalpointsxed. Thisfactsoimpressed itsdiscoverersthatthey
believedtheyhaduncoveredtheunifying principle oftheuniverse.Maup er-
tuis,forone,triedtobaseaproofoftheexistence ofGodonit.Todaythe
action integral, through itsstarring roleintheFeynman pathintegral for-
mulation ofquantummechanics, remains attheheart oftheoretical physics.
1.3.1 OneDegree ofFreedom
Wewillnotattempt toderiveLagrange from Newton andD'Alem bert's
extension oftheprinciple ofvirtual work{leavingthistasktoamechanics
course |butwillsatisfy ourselv eswithsome examples whichillustrate the
computational advantages ofLagrange's approac h,aswellasasubtle pitfall.
Example: Atwood'sMachine.Thisdevice, inventedin1784butstillafa-
miliar sightinundergraduate laboratories, isusedtodemonstrate Newton's
lawsofmotion andtomeasureg.Itconsists oftwoweightsconnected bya
lightstring whichpasses overalightandfrictionless pulley .
1.3.LAGRANGIAN MECHANICS 11
m21mT1x
Tx2g
Theelemen taryapproac histowrite anequation ofmotion foreachofthe
twoweights
m1x1=m1g T;
m2x2=m2g T: (1.38)
Wethentakeintoaccoun ttheconstrain t_x1= _x2toget
m1x1=m1g T;
m2x1=m2g T: (1.39)
Finally weeliminate theconstrain tforce, thetensionT,togettheaccelera-
tion
(m1+m2)x1=(m1 m2)g: (1.40)
TheLagrangian solution takestheconstrain tintoaccoun tfromthevery
beginning byintroducing asingle generalized coordinateq=x1= x2,and
writing
L=T V=1
2(m1+m2)_q2 (m2 m1)gq: (1.41)
Fromthisweobtain asingle equation ofmotion
d
dt @L
@_qi!
@L
@qi=0)(m1+m2)q=(m1 m2)g: (1.42)
12 CHAPTER 1.CALCULUS OFVARIA TIONS
TheadvantageofthetheLagrangian metho disthatconstrain tforces, which
dononetwork,neverappear.Thedisadv antageisexactly thesame: ifwe
needtondtheconstrain tforces {inthiscasethetension inthestring |
wecannot useLagrange alone.
Example: PolarCoordinates
ϑry
xaraϑ
Consider acentralforce problem withFr= @rV(r).TheNewtonian
metho dbegins bycomputing theacceleration inpolarcoordinates. This
ismosteasily donebysettingz=reianddieren tiating twice:
_z=(_r+ir_)ei;
z=(r r_2)ei+i(2_r_+r)ei: (1.43)
Reading othecomponentsparallel andperpendicular toeigivesforthe
acceleration
ar=r r_2;
a=r+2_r_; (1.44)
Newton's equations therefore become
m(r r_2)= @V
@r
m(r+2_r_)=0;)d
dt(mr2_)=0: (1.45)
Settingl=mr2_,theconserv edangular momen tum,andeliminating _gives
mr l2
mr3= @V
@r: (1.46)
1.3.LAGRANGIAN MECHANICS 13
(IfthiswereKepler's problem, whereV=GmM=r,wewouldnowproceed
tosimplify thisequation bysubstituting r=1=u,butthatisanother story.)
FollowingLagrange werstcompute thekinetic energy inpolarcoordi-
nates (thisrequires onelessderivativethancomputing theacceleration) and
set
L=T V=1
2m(_r2+r2_2) V(r): (1.47)
TheEuler-Lagrange equations arenow
d
dt @L
@_r!
@L
@r=0;)mr r2_2+@V
@r=0
d
dt @L
@_!
@L
@=0;)d
dt(mr2_)=0: (1.48)
Therstintegral forthisproblem is
E=_r@L
@_r+_@L
@_r L
=1
2m(_r2+r2_2)+V(r) (1.49)
whichisthetotalenergy .Thustheconstancy oftherstintegral states that
dE
dt=0; (1.50)
orthatenergy isconserv ed.
Warning :Wemightrealize, without havinggonetothetrouble ofderiving
itfromtheLagrange equations, thatrotational invariance guaran teesthat
theangular momen tuml=mr2_willbeaconstan t.Havingdoneso,itis
almost irresistible totrytoshort-circuit some ofthearithmetic byplugging
thispriorknowledge into
L=1
2m(_r2+r2_2) V(r) (1.51)
soastoeliminate thevariable _infavouroftheconstan tl.Ifwetrythiswe
get
L?!1
2m_r2+l2
mr2 V(r): (1.52)
14 CHAPTER 1.CALCULUS OFVARIA TIONS
Wecannowdirectly write downtheLagrange equationr,whichis
mr+l2
mr3?= @V
@r: (1.53)
Unfortunately thishasthewrong signbeforethel2=mr3term! Thelesson is
thatwemustbeverycareful inusing consequences ofavariational principle
tomodifytheprinciple. Itcanbedone, andinmechanics itleads tothe
Routhian or,inmore modernlanguage toHamiltonian reduction ,butit
requires using aLegendre transform. Thereader should consult abookon
mechanics fordetails.
1.3.2 Noether's Theorem
Thetime-indep endence oftherstintegral
d
dt(
_q@L
@_q L)
=0; (1.54)
andofangular momen tum
d
dtfmr2_g=0; (1.55)
areexamples ofconservation laws.Weobtained them bothbymanipulating
theEuler-Lagrange equations ofmotion, butalsoindicated thattheywere
insome wayconnected withsymmetries. Oneofthechiefadvantages ofa
variational formulation ofaphysical problem isthatthisconnection
Symmetry,Conserv ation Law
canbemade explicit byexploiting astrategy duetoEmmyNoether. She
showedhowtoproceeddirectly fromtheaction integral totheconserv ed
quantitywithout havingtoddle aboutwiththeequations ofmotion. We
beginbyillustrating hertechnique inthecaseofangular momen tum,whose
conserv ation isaconsequence therotational symmetry ofthecentralforce
problem. Theaction integral forthecentralforceproblem is
S=ZT
01
2m(_r2+r2_2) V(r)
dt: (1.56)
Noetherobserv esthattheintegrand isleftunchanged ifwemakethevariation
(t)!(t)+ (1.57)
1.3.LAGRANGIAN MECHANICS 15
whereisaxedangle andisasmall, time-indep enden t,parameter. This
invariance isthesymmetry weshallexploit. Itisamathematical identity:
itdoesnotrequire thatrandobeytheequations ofmotion. Shenext
observ esthatsincetheequations ofmotion areequivalenttothestatemen t
thatSisleftstationary under anyinnitesimal variations inrand,they
necessarily imply thatSisstationary under thespecic variation
(t)!(t)+(t) (1.58)
where nowisallowedtobetime-dep enden t.Thisstationarit yoftheaction
isnolonger amathematical identity,but,because itrequiresr,,toobey
theequations ofmotion, hasphysical content.Inserting=(t)intoour
expression forSgives
S=ZT
0n
r2_o
_dt: (1.59)
Notethatthisvariation dependsonlyonthetimederivativeof,andnot
itself. Thisisbecause oftheinvariance ofSunder time-indep enden trota-
tions. Wenowassume that(t)=0att=0andt=T,andintegrate by
parts totakethetimederivativeoandputitontherestoftheintegrand:
S= Z(d
dt(r2_))
(t)dt: (1.60)
Since theequations ofmotion saythatS=0under allinnitesimal varia-
tions, andinparticular those duetoanytimedependen trotation(t),we
deduce thattheequations ofmotion imply thatthecoecien tof(t)must
bezero,andso,providedr(t),(t),obeytheequations ofmotion, wehave
0=d
dt(r2_): (1.61)
Asasecond illustration wederiveenergy (rstintegral) conserv ation for
thecasethatthesystem isinvariantunder timetranslations |meaning
thatLdoesnotdependexplicitly ontime. Inthiscasetheaction integral
isinvariantunder constan ttimeshiftst!t+intheargumen tofthe
dynamical variable:
q(t)!q(t+)q(t)+_q: (1.62)
Theequations ofmotion tellusthatthattheaction willbestationary under
thevariation
q(t)=(t)_q; (1.63)
16 CHAPTER 1.CALCULUS OFVARIA TIONS
where again wenowpermit theparametertodependont.Weinsert this
variation into
S=ZT
0Ldt (1.64)
andnd
S=ZT
0(@L
@q_q+@L
@_q(q+_q_))
dt: (1.65)
Thisexpression containsundotted's.Because ofthisthechange inSisnot
obviously zerowhenistimeindependen t|buttheabsence ofanyexplicit
tdependence inLtellsusthat
dL
dt=(@L
@q_q+@L
@_qq)
: (1.66)
Asaconsequence, fortimeindependen t,wehave
S=ZT
0(
dL
dt)
dt=[L]T
0; (1.67)
showingthatthechange inScomes entirely fromtheendpointsofthetime
interval.These xedendpointsexplicitly break time-translation invariance,
butinatrivial manner. Forgeneral(t)wehave
S=ZT
0(
(t)dL
dt+@L
@_q_q_)
dt: (1.68)
Thisequation isanidentity.Itdoesnotrelyonqobeying theequation of
motion. After anintegration byparts, taking(t)tobezeroatt=0;T,it
isequivalentto
S=ZT
0(t)d
dt(
L @L
@_q_q)
dt: (1.69)
Nowweassume thatq(t)doesobeytheequations ofmotion. Thevariation
principle thensaysthatS=0forany(t),andwededuce thatforq(t)
satisfying theequations ofmotion wehave
d
dt(
L @L
@_q_q)
=0: (1.70)
Thegeneral strategy thatconstitutes \Noether's theorem" mustnowbe
obvious: welookforaninvariance oftheaction under asymmetry trans-
formation withatime-indep enden tparameter. Wethenobserv ethatifthe
1.3.LAGRANGIAN MECHANICS 17
dynamical variables obeytheequations ofmotion, thentheaction principle
tellsusthattheaction willremain stationary under suchavariation ofthe
dynamical variables evenaftertheparameter ispromoted tobeingtimede-
penden t.Theresultan tvariation ofScanonlydependontimederivativesof
theparameter. Weintegrate bypartssoastotakeallthetimederivativeso
it,andontotherestoftheintegrand. Since theparameter isarbitrary ,we
deduce thattheequations ofmotion tellusthatthatitscoecien tinthein-
tegral mustbezero.Since thiscoecien tisthetimederivativeofsomething,
thissomething isconserv ed.
1.3.3 ManyDegrees ofFreedom
Theextension oftheaction principle tomanydegrees offreedom isstraigh t-
forward.Asanexample consider thesmall oscillations aboutequilibrium of
asystem withNdegrees offreedom. Weparametrize thesystem interms of
deviations fromtheequilibrium position andexpand outtoquadratic order.
Weobtain aLagrangian
L=NX
i;j=11
2Mij_qi_qj 1
2Vijqiqj
; (1.71)
whereMijandVijareNNsymmetric matrices encodingtheinertial and
potentialenergy properties ofthesystem. Nowwehaveoneequation
0=d
dt @L
@_qi!
@L
@qi=NX
j=1
Mijqj+Vijqj
(1.72)
foreachi.
1.3.4 ContinuousSystems
Theaction principle canbeextended toeldtheories andtocontinuumme-
chanics. Hereonehasacontinuousinnit yofdynamical degrees offreedom,
either oneforeachpointinspace andtimeoroneforeachpointinthemate-
rial,buttheextension ofthevariational derivativetofunctions ofmorethan
onevariable should possess noconceptual diculties.
SupposewearegivenanactionSdepending onaeld'(x)anditsrst
derivatives
'@'
@x: (1.73)
18 CHAPTER 1.CALCULUS OFVARIA TIONS
Herex,=0;1;:::;d,arethecoordinates ofd+1dimensional space-time.
Itistraditional totakex0tandtheother coordinates spacelik e.Suppose
further that
S=Z
Ldt=Z
L(';')dd+1x; (1.74)
whereListheLagrangian density ,interms ofwhich
L=Z
Lddx; (1.75)
where theintegral isoverthespace coordinates. Now
S=Z(
'(x)@L
@'(x)+('(x))@L
@'(x))
dd+1x
=Z
'(x)(@L
@'(x) @
@x @L
@'(x)!)
dd+1x: (1.76)
Ingoing fromtherstlinetothesecond, wehaveobserv edthat
('(x))=@
@x'(x) (1.77)
andusedthedivergence theorem,
Z
@A
@x!
dn+1x=Z
@
AndS; (1.78)
where
issome space-time region and@
itsboundary ,tointegrate by
parts. HeredSistheelemen tofareaontheboundary ,andntheoutward
normal. Asbefore, wetake'tovanish ontheboundary ,andhence there
isnoboundary contribution tovariation ofS.Theresult isthat
S
'(x)=@L
@'(x) @
@x @L
@'(x)!
; (1.79)
andtheequation ofmotion comes fromsetting thistozero.Notethatasum
overtherepeated coordinate indexisimplied. Inpractice, however,itis
easier nottousethisformula,butinstead dothevariation explicitly asin
thefollowingexamples.
1.3.LAGRANGIAN MECHANICS 19
TheVibrating string
Thesimplest continuousdynamical system isthevibrating string. Wede-
scribethestring displacemen tbyy(x;t).
0 Ly(x,t)
Letussupposethatthestring hasxedends, amassperunitlength of,and
isunder tensionT.Ifweassume onlysmall displacemen tsfromequilibrium,
theLagrangian is
L=ZL
0dx1
2_y2 1
2Ty02
: (1.80)
Thevariation oftheaction is
S=ZZL
0dtdxf_y_y Ty0y0g
=ZZL
0dtdxfy(x;t)( y+Ty00)g: (1.81)
Toreachthesecond linewehaveintegrated byparts, and,because theends
arexed, andthereforey=0atx=0andL,there isnoboundary term.
Requiring thatS=0forallallowedvariationsythengivestheequation
ofmotion
@2y
@t2 T@2y
@x2=0: (1.82)
Thisisthewaveequation forwaveswithspeedc=q
T=.Observ ethat
from(1.81) wecanreadothefunctional derivativeofSwithrespecttothe
variabley(x;t)asbeing
S
y(x;t)= y(x;t)+Ty00(x;t): (1.83)
Inwriting downtherstintegral forthiscontinuoussystem, wemust
replace thesumoverdiscrete indices byanintegral:
E=X
i_qi@L
@_qi L!Z
dx(
_y(x)L
_y(x))
L: (1.84)
20 CHAPTER 1.CALCULUS OFVARIA TIONS
When computing L=_y(x)from
L=ZL
0dx1
2_y2 1
2Ty02
;
wemustremem berthatitisthecontinuousanalogue of@L=@_qi,andso,in
contrasttowhat wedowhen computating S=y(x),wemusttreat _y(x)as
avariable independen tofy(x).Wethenhave
L
_y(x)=_y(x); (1.85)
leading to
E=ZL
0dx1
2_y2+1
2Ty02
: (1.86)
This, asexpected, isthetotalenergy ,kinetic pluspotential,ofthestring.
Exercise :Consider anaction oftheform
S=Z
dd+1xL(';@') (1.87)
whichdoesnotdependexplicitly onx.Generalize theNoether derivation
oftheenergy conserv ation lawtooneexploiting variations oftheform
'=(x)@'; (1.88)
wheredependsonspace andtime, andhence showthat
@T
=0; (1.89)
where
T
=@L
@(@')@'