Kontinuumsgrenze der Euler-Lagrange-Gleichung für die Lagrange-Dichte des harmonischen 1D-Gitters

Ich versuche, einer Ableitung der Euler-Lagrange-Gleichung an der Kontinuumsgrenze zu folgen, und finde einige Details schwer verständlich. Das 1D-Gitter hat eine monoatomare Basis mit Atomabstand A und Atommassenwesen M . Die relative Verschiebung des n-ten Atoms ist δ η N . Der Lagrange-Operator am diskreten Grenzwert wird geschrieben als

L = N = 1 N A 1 2 [ M A ( δ η ˙ N ) 2 κ A ( δ η N + 1 δ η N A ) 2 M A Ω 2 ( δ η N ) 2 ] =: N = 1 N A L N

Wo κ Und Ω 2 sind die Stärke der Atombindung bzw. des äußeren Feldes. An der Grenze von N , die zeitliche Ableitung der relativen Verschiebung δ η N zum Zeitpunkt T wurde durch den Wert der zeitlichen Ableitung ersetzt ( T φ ) an der Raum-Zeit-Koordinate ( X , T ) . Und wir haben eine periodische Randbedingung von

φ ( X + L , T ) = φ ( X , T ) , X ] 0 , L ] , T R .

Wenn wir lassen

μ := M A , Y := κ A

und ersetze die diskrete Summe N durch das Integral D X / A über das halboffene Intervall ] 0 , L ] , dann schreiben wir die Lagrange-Funktion als

L = 0 L D X 1 2 [ μ ( φ T ) 2 Y ( φ X ) 2 μ Ω 2 φ 2 ] = : 0 L D X L .

Wo L ist die "Lagrange-Dichte". Aus dem Buch von Christopher Mudry erhält man die Kontinuumsgrenze der Euler-Lagrange-Gleichungen als

(1) T δ L ( X , T ) δ ( T φ ) ( j , T ) + X δ L ( X , T ) δ ( X φ ) ( j , T ) = δ L ( X , T ) δ φ ( j , T ) .

Um (1) zu erhalten, erweitere ich zuerst L bezüglich φ , ( X φ ) , Und ( T φ ) geben

(2) δ L = L [ φ + δ φ , ( X φ ) + δ ( X φ ) , ( T φ ) + δ ( T φ ) ] L [ φ , ( X φ ) , ( T φ ) ] = L φ δ φ + L ( X φ ) δ ( X φ ) + L ( T φ ) δ ( T φ ) + .

Ignorieren Sie die Terme höherer Ordnung in (2) und beachten Sie das δ X φ = X δ φ , können wir (2) durch Teile integrieren, um zu geben

(3) δ L = 0 L D X δ L = 0 L D X { L φ δ φ + X ( L ( X φ ) δ φ ) + T ( L ( T φ ) δ φ ) δ φ X ( L ( X φ ) ) δ φ T ( L ( T φ ) ) } .

Um (1) aus (3) zu erhalten, muss das folgende Integral erfüllt sein:

(5) 0 L D X X ( L ( X φ ) δ φ ) = [ L ( X φ ) δ φ ] 0 L = 0 ,

(6) 0 L D X T ( L ( T φ ) δ φ ) = D D T 0 L D X ( L ( T φ ) δ φ ) = 0 .

Obwohl ich kein Problem mit (5) habe, da wir die Endpunkte nicht ändern, weiß ich nicht, warum Gl. (6) gilt . Ist es nur, weil wir Zeit behandeln T als fester Parameter in (6)?

Hallo Lonitch, ich studiere auch die stetige Grenze der Euler-Lagrange-Gleichungen und habe mich gefragt, welche Bibliographie Sie verwendet haben (beiseite von Christopher Mudrys Buch).
@JD_PM, ich habe Mudrys Buch verwendet, um mich mit einigen QFT-Problemen vertraut zu machen. Wenn Ihr Ziel darin besteht, die Euler-Lagrange-Gleichung und ihre Anwendung in der analytischen Mechanik zu lernen, würde ich als Einstieg „Die Variationsprinzipien der Mechanik“ von Cornelius Lanczos empfehlen. Wem das noch nicht reicht, der kann „Analytical Mechanics: An Introduction“ von Antonio Fasano und Stefano Marni ausprobieren. Ihr Buch hat eine sehr gute Berichterstattung über die Integrierbarkeit. Falls Sie an der Beziehung zwischen analytischer und Quantenmechanik interessiert sind, bietet das OGT-Buch von Oliver Johns eine gute Einführung.
Lonitch danke für die Empfehlung. Eigentlich interessiere ich mich besonders dafür, wie man die diskreten Euler-Lagrange-Gleichungen herleiten kann (mehr Details physicalforums.com/threads/… ). Ich habe kürzlich etwas über 3D-Gitterdiskretisierung gelernt und bin immer noch dabei, den Dreh raus zu haben. In meinen Grundkursen über Mechanik (wir konzentrierten uns auf Gregorys Buch über klassische Mechanik) haben wir eine solche Technik nicht behandelt, deshalb suche ich nach einer Bibliographie. Beschäftigt sich eine Ihrer empfohlenen Quellen mit 3D-Gitterdiskretisierung?
@JD_PM Es klingt für mich, als würden Sie versuchen, den Variationsintegrator auf einem 3D-Gitter zu lernen. Keines der Bücher gibt solche. Vielleicht möchten Sie einige Papiere dafür überprüfen. Zum Beispiel das zweite Kapitel der Doktorarbeit von Dr. Mathew West ( thesis.library.caltech.edu/2492/1/west_thesis.pdf )

Antworten (2)

OP hat Recht. Im Allgemeinen müssen wir, damit das Aktionsprinzip funktioniert, nur Randbedingungen (BCs) an der Grenze der Raumzeit auferlegen. Daher sollten wir Gl. (6) im Inneren der Raumzeit.

Mit anderen Worten, OP sollte seine obige Analyse vorzugsweise wiederholen, indem es die Aktion und nicht die Lagrange verwendet.

Danke für deinen Rat. Ich denke darüber nach, dass die Störung zu jedem Zeitpunkt in den Pfad eingeführt wird δ φ sollte nur eine Funktion von sein X . Und die Euler-Lagrange-Gleichung wird Gl. (1) ohne die T Begriff. Mein Versuch, dies zu korrigieren, wird in meiner Antwort gezeigt.
Beachten Sie, dass die infinitesimale Variation δ φ ( X , T ) kann auch abhängen T .

Mein Versuch, die obige Euler-Lagrange-Gleichung herzuleiten, geht von der Aktion aus,

(1) S = T 1 T 2 0 L L ( φ , X φ , T φ ) D X D T

und ich führe die Störung in den Pfad ein, δ φ ( X , T ) schreiben

S λ = T 1 T 2 0 L D X D T L ( φ + λ δ φ , X φ + λ X δ φ , T φ + λ T δ φ ) .

Dann

(2) S λ λ = T 1 T 2 0 L D X D T { δ φ L ( φ + λ δ φ ) + X ( δ φ ) L ( X φ + λ X ( δ φ ) ) + T ( δ φ ) L ( T φ + λ T ( δ φ ) ) }

Und

(3) S λ λ | λ = 0 = T 1 T 2 0 L D X D T { δ φ L φ + X ( δ φ ) L ( X φ ) + T ( δ φ ) L ( T φ ) } = 0 L [ δ φ L ( T φ ) ] T 1 T 2 D X + T 1 T 2 0 L D X D T { δ φ L φ δ φ T [ L ( T φ ) ] + X ( δ φ ) L ( X φ ) } = T 1 T 2 0 L D X D T { δ φ [ L φ T ( L ( T φ ) ) X ( L ( X φ ) ) ] } + T 1 T 2 [ δ φ L ( X φ ) ] 0 L D T = T 1 T 2 0 L D X D T { δ φ [ L φ T ( L ( T φ ) ) X ( L ( X φ ) ) ] } .

So

(4) S λ λ | λ = 0 = T 1 T 2 0 L D X D T { δ φ [ L φ T ( L ( T φ ) ) X ( L ( X φ ) ) ] } = 0 .

Aus (4) haben wir

(5) L φ = T ( L ( T φ ) ) + X ( L ( X φ ) ) .