Poincaré-Lemma

Das Poincaré-Lemma ist ein Satz aus der Mathematik und wurde nach dem französischen Mathematiker Henri Poincaré benannt.

Exakte und geschlossene Differentialformen

• Eine Differentialform $\omega$ vom Grad $k$ heißt geschlossen, falls $\mathrm {d} \omega =0$ gilt. Dabei bezeichnet $\mathrm {d}$ die äußere Ableitung.
• Eine Differentialform $\omega$ vom Grad $k$ heißt exakt, falls es eine $(k-1)$-Differentialform $\nu$ gibt, so dass $\omega =\mathrm {d} \nu$ gilt. Die Form $\nu$ nennt man eine Potentialform von $\omega .$

Die Potentialform ist nicht eindeutig bestimmt, sondern nur "bis auf Umeichung" (siehe unten).

Wegen $\mathrm {d} \circ \mathrm {d} \equiv 0$ ist jede exakte Differentialform auch geschlossen. Das Poincaré-Lemma gibt Voraussetzungen an, unter denen auch die umgekehrte Aussage gilt. Beim Beweis ergibt sich darüber hinaus eine Verallgemeinerung des Lemmas: Von jeder Differentialform lässt sich „per Konstruktion“ ein exakter Anteil abspalten.

Aussage

Das Poincaré-Lemma besagt, dass jede auf einer sternförmigen offenen Menge $U\subseteq \mathbb {R} ^{d}$ definierte geschlossene Differentialform exakt ist.

Die Aussage lässt sich abstrakter auch so formulieren: Für eine sternförmige offene Menge $U\subseteq \mathbb {R} ^{d}$ verschwindet die $k$-te De-Rham-Kohomologie für alle $k>0$:

$\mathrm {H} _{\mathrm {dR} }^{k}(U)=0$

Im dreidimensionalen Spezialfall besagt das Poincaré-Lemma, in die Sprache der Vektoranalysis überführt, dass ein auf einem einfach-zusammenhängenden Gebiet definiertes wirbelfreies Vektorfeld als Gradient eines Potentialfeldes $\Phi (\mathbf {r} )$ ($k=1$), ein quellfreies Vektorfeld auf einem konvexen Gebiet durch Rotation eines Vektorpotentials ${\vec {A}}(\mathbf {r} ,t)$ ($k=2$), und eine skalare Felddichte („Quellendichte“) als Divergenz eines Vektorfeldes ($k=3$) dargestellt werden können.

Beweis (konstruktiv)

Sei $x_{0}\in U$ der Punkt, um welchen herum $U\subset \mathbb {R} ^{n}$ sternförmig ist. Das Poincaré-Lemma gibt explizit eine $(k-1)$-Form an, und zwar mit folgender Formel: Einer beliebigen $k$-Form $\textstyle \omega ^{k}=\sum \omega _{I}{\rm {d}}x_{I}$ lässt sich, Geschlossenheit nicht notwendig vorausgesetzt, eine $(k-1)$-Form $P^{k-1}(\omega ^{k})$ zuordnen, aus der sich bei Geschlossenheit die gesuchte Potentialform ergibt: Diese zugeordnete Form lässt sich durch folgende Abbildung definieren:

$P^{k-1}(\omega ^{k})(x)\,:=\,\sum _{i_{1}<\cdots <i_{k}}\sum _{\alpha =1}^{k}(-1)^{\alpha -1}{\Big (}\int _{0}^{1}t^{k-1}\omega _{i_{1}\cdots i_{k}}(x_{0}+tx)dt{\Big )}x^{i_{\alpha }}{\rm {d}}x^{i_{1}}\wedge \cdots \wedge {\widehat {{\rm {d}}x^{i_{\alpha }}}}\wedge \cdots \wedge {\rm {d}}x^{i_{k}}$.

(Das Dachsymbol in der $i_{\alpha }$-ten Spalte der rechten Seite bedeutet, dass das entsprechende Differential ausgelassen wird.)

Nun zeigt man direkt, dass folgende Identität gilt: $\omega ^{k}\equiv \mathrm {P} ^{k}({\rm {d}}\omega ^{k})+{\rm {d}}{\mathrm {P} ^{k-1}(}\omega ^{k}\mathrm {)} \,,$ was formal der Produktregel der Differentiation entspricht und die durch $\omega ^{k}$ repräsentierten Eigenschaften in zwei Anteile zerlegt, von denen der zweite die gesuchte Eigenschaft besitzt.

Wegen der Voraussetzung ${\rm {d}}\omega ^{k}\equiv 0$ und wegen $\mathrm {d} \circ \mathrm {d} =0$ ergibt sich zunächst $0\equiv {\mathrm {d} }P^{k}({\mathrm {d} {\omega }^{k}}\to 0).$ Dies gilt ohne Einschränkung der Allgemeinheit auch ohne das vorderste $\mathrm {d}$ der rechten Seite, und zwar deshalb, weil durch die Forderung $\mathrm {d} \omega ^{k}\to 0$ die Form $\mathrm {d} P^{k}$ nur am Nullpunkt betrachtet wird, sodass wie beim Totalen Differential einer Funktion aus $\mathrm {d} P^{k}=0$ bis auf sog. Eichtransformationen (siehe unten) auch $\mathrm {P} ^{k}=0$ gefolgert werden kann.

Somit bleibt nur der letzte Term der obigen Identität, und es folgt die gesuchte Aussage: $\omega ^{k}\equiv \mathrm {d} \eta ^{k-1},$ mit $\eta ^{k-1}\,:=\,\mathrm {P} ^{k-1}(\omega ^{k})\,.$

Die angegebene Identität verallgemeinert zugleich das Poincarésche Lemma durch Zerlegung einer beliebigen Differentialform $\omega$ in einen nicht-exakten („anholonomen“) und einen exakten („holonomen“) Anteil (die eingeklammerten Bezeichnungen entsprechen den sog. Zwangskräften in der analytischen Mechanik). Es entspricht zugleich der Zerlegung eines beliebigen Vektorfeldes in einen Wirbel- und einen Quellen-Anteil.

In der Sprache der homologischen Algebra ist $P$ eine kontrahierende Homotopie, die z. B. auf den zentralen Punkt des hier betrachteten sternförmigen Gebietes kontrahiert.

Umeichung

Das so definierte $\eta ^{k-1}$ ist nicht die einzige $(k-1)$-Form, deren äußeres Differential ${\omega }^{k}$ ist. Alle anderen unterscheiden sich aber höchstens um das Differential einer $(k-2)$-Form voneinander: Sind $\eta _{2}^{k-1}$ und $\eta _{1}^{k-1}$ zwei solche $(k-1)$-Formen, so existiert eine $(k-2)$-Form $\xi ^{k-2}$ derart, dass $\eta _{2}^{k-1}=\eta _{1}^{k-1}+\mathrm {d} \xi ^{k-2}$ gilt.

Der Zusatz $+\,\mathrm {d} \xi ^{k-2}$ wird auch als Eichtransformation bzw. Umeichung von $\eta _{1}^{k-1}$ bezeichnet.

Anwendung in der Elektrodynamik

Aus der Elektrodynamik ist der Fall eines von einem stationären Strom erzeugten Magnetfeldes bekannt, mit dem sog. Vektorpotential ${\vec {A}}(\mathbf {r} )\,.$ Dieser Fall entspricht $k=2$, wobei das sternförmige Gebiet der $\mathbb {R} ^{3}$ ist. Der Vektor der Stromdichte ist ${\vec {j}}$ und entspricht der Stromform $\mathbf {I} :=j_{1}(x,y,z){\rm {d}}x_{2}\wedge {\rm {d}}x_{3}+j_{2}(x,y,z){\rm {d}}x_{3}\wedge {\rm {d}}x_{1}+j_{3}(x,y,z){\rm {d}}x_{1}\wedge {\rm {d}}x_{2}\,.$ Für das Magnetfeld ${\vec {B}}$ gilt Analoges: es entspricht der Magnetflussform $\Phi _{B}:=B_{1}{\rm {d}}x_{2}\wedge {\rm {d}}x_{3}+\dots$ und lässt sich aus dem Vektorpotential ableiten: $\textstyle {\vec {B}}=\operatorname {rot} {\vec {A}}=\left({\tfrac {\partial A_{3}}{\partial x_{2}}}-{\tfrac {\partial A_{2}}{\partial x_{3}}},{\tfrac {\partial A_{1}}{\partial x_{3}}}-{\tfrac {\partial A_{3}}{\partial x_{1}}},{\tfrac {\partial A_{2}}{\partial x_{1}}}-{\tfrac {\partial A_{1}}{\partial x_{2}}}\right)^{t}$, oder $\Phi _{B}={\rm {d}}\mathbf {A}$. Dabei entspricht das Vektorpotential ${\vec {A}}$ der Potentialform $\mathbf {A} :=A_{1}{\rm {d}}x_{1}+A_{2}{\rm {d}}x_{2}+A_{3}{\rm {d}}x_{3}\,.$ Die Geschlossenheit der Magnetflussform entspricht der Quellenfreiheit des Magnetfeldes $(\operatorname {div} {\vec {B}}\equiv 0\,).$

Unter Verwendung der Coulomb-Eichung $\operatorname {div} {\vec {A}}{\stackrel {!}{=}}0$ bzw. passend zu $\operatorname {div} {\vec {j}}{\stackrel {!}{=}}0$ gilt dann für i=1,2,3

$A_{i}({\vec {r}})=\int {\frac {\mu _{0}j_{i}({\vec {r}}^{\,'})\,\,dx_{1}'dx_{2}'dx_{3}'}{4\pi |{\vec {r}}-{\vec {r}}^{\,'}|}}\,,$

dabei ist $\mu _{0}$ eine Naturkonstante, die sogenannte Magnetische Feldkonstante.

An dieser Gleichung ist u. a. bemerkenswert, dass sie vollständig einer bekannten Formel für das elektrische Feld ${\vec {E}}$ entspricht, dem Coulombpotential $\,\phi (x_{1},x_{2},x_{3})$ einer gegebenen Ladungsverteilung mit der Dichte $\rho (x_{1},x_{2},x_{3})$. Man vermutet an dieser Stelle bereits, dass

• ${\vec {E}}$ und ${\vec {B}}$ bzw.
• $\rho$ und ${\vec {j}}$ sowie
• $\,\phi$ und ${\vec {A}}$

zusammengefasst werden können und dass sich die relativistische Invarianz der Maxwellschen Elektrodynamik daraus ergibt, siehe dazu Elektrodynamik.

Wenn man die Bedingung der Stationarität aufgibt, muss auf der linken Seite der obigen Gleichung bei $A_{i}$ zu den Raumkoordinaten das Zeitargument $t$ hinzugefügt werden, während auf der rechten Seite in $j_{i}'$ die sog. „retardierte Zeit“ $t':=t-{\tfrac {|{\vec {r}}-{\vec {r}}^{\,'}|}{c}}$ zu ergänzen ist. Es wird dabei wie zuvor über die drei Raumkoordinaten ${\vec {r}}^{\,'}$ integriert. Schließlich ist $c$ die Lichtgeschwindigkeit im Vakuum.

Anwendung in der Kontinuumsmechanik

In der Kontinuumsmechanik wird das Lemma auf Tensoren angewendet, was z. B. für die Aufstellung der Kompatibilitätsbedingungen gebraucht wird. Ausgangspunkt ist das Lemma in der Formulierung:

$\operatorname {rot} {\vec {u}}={\hat {e}}_{k}\times {\frac {\partial {\vec {u}}}{\partial x_{k}}}={\vec {0}}\quad \rightarrow \quad \exists \varphi \colon {\vec {u}}=\operatorname {grad} \varphi$		(I)

Der Operator „grad“ bildet den Gradient, die Vektoren ${\hat {e}}_{1,2,3}$ sind die Standardbasis des kartesischen Koordinatensystems mit Koordinaten $x_{1,2,3}$ und es wurde die einsteinsche Summenkonvention angewendet, dergemäß über in einem Produkt doppelt vorkommende Indizes, hier k, von eins bis drei zu summieren ist, was auch im Folgenden praktiziert werden soll.

Gegeben sei nun ein Tensorfeld $\mathbf {T} ={\hat {e}}_{i}\otimes {\vec {t}}_{i}$, dessen Zeilenvektoren ${\vec {t}}_{1,2,3}$ mit dem dyadischen Produkt „⊗“ zum Tensor zusammengefügt werden. Jeder Tensor zweiter Stufe kann in dieser Form dargestellt werden. Die Rotation des Tensors verschwinde

${\begin{aligned}\operatorname {rot} (\mathbf {T} ):=\nabla \times (\mathbf {T} ^{\top })=&{\hat {e}}_{k}\times {\frac {\partial }{\partial x_{k}}}({\vec {t}}_{i}\otimes {\hat {e}}_{i})=\left({\hat {e}}_{k}\times {\frac {\partial {\vec {t}}_{i}}{\partial x_{k}}}\right)\otimes {\hat {e}}_{i}=\mathbf {0} \\&\rightarrow \quad {\hat {e}}_{k}\times {\frac {\partial {\vec {t}}_{i}}{\partial x_{k}}}={\vec {0}}\,,\quad i=1,2,3\end{aligned}}$

so dass also jeder Zeilenvektor rotationsfrei ist. Dann gibt es für jeden Zeilenvektor ein Skalarfeld $u_{i}$, dessen Gradient er ist:

${\vec {t}}_{i}=\operatorname {grad} u_{i}\quad \rightarrow \quad \mathbf {T} ={\hat {e}}_{i}\otimes {\vec {t}}_{i}={\hat {e}}_{i}\otimes \operatorname {grad} u_{i}=\operatorname {grad} {\vec {u}}\,,$

denn der Gradient des Vektors ${\vec {u}}:=u_{i}{\hat {e}}_{i}$ bildet sich gemäß:

$\operatorname {grad} {\vec {u}}:={\frac {\partial u_{i}}{\partial x_{k}}}{\hat {e}}_{i}\otimes {\hat {e}}_{k}={\hat {e}}_{i}\otimes {\frac {\partial u_{i}}{\partial x_{k}}}{\hat {e}}_{k}={\hat {e}}_{i}\otimes \operatorname {grad} u_{i}\,.$

Damit gilt die zweite Form des Lemmas:

$\operatorname {rot} (\mathbf {T} ):=\nabla \times (\mathbf {T} ^{\top })=\mathbf {0} \quad \rightarrow \quad \exists {\vec {u}}\colon \mathbf {T} =\operatorname {grad} {\vec {u}}$		(II)

Wenn zusätzlich die Spur des Tensors verschwindet, dann ist das Vektorfeld divergenz­frei:

$\operatorname {Sp} (\mathbf {T} )=\operatorname {Sp} ({\hat {e}}_{i}\otimes \operatorname {grad} u_{i})={\hat {e}}_{i}\cdot {\frac {\partial u_{i}}{\partial x_{k}}}{\hat {e}}_{k}={\frac {\partial u_{i}}{\partial x_{i}}}=\operatorname {div} {\vec {u}}=0\,.$

In diesem Fall berechnet sich mit dem Einheitstensor 1 = ê_j ⊗ ê_j:

${\begin{aligned}\operatorname {rot} (\mathbf {1} \times {\vec {u}})=&{\hat {e}}_{k}\times {\frac {\partial }{\partial x_{k}}}[({\hat {e}}_{j}\otimes {\hat {e}}_{j})\times u_{i}{\hat {e}}_{i}]^{\top }={\frac {\partial u_{i}}{\partial x_{k}}}[{\hat {e}}_{k}\times ({\hat {e}}_{j}\times {\hat {e}}_{i})]\otimes {\hat {e}}_{j}\\=&{\frac {\partial u_{i}}{\partial x_{k}}}(\delta _{ik}{\hat {e}}_{j}-\delta _{jk}{\hat {e}}_{i})\otimes {\hat {e}}_{j}={\frac {\partial u_{i}}{\partial x_{i}}}{\hat {e}}_{j}\otimes {\hat {e}}_{j}-{\frac {\partial u_{i}}{\partial x_{j}}}{\hat {e}}_{i}\otimes {\hat {e}}_{j}\\=&-\operatorname {grad} {\vec {u}}\end{aligned}}$

und der Tensor $\mathbf {1} \times {\vec {u}}$ ist schiefsymmetrisch:

$(\mathbf {1} \times {\vec {u}})^{\top }=({\hat {e}}_{i}\otimes {\hat {e}}_{i}\times u_{j}{\hat {e}}_{j})^{\top }=\epsilon _{ijk}u_{j}{\hat {e}}_{k}\otimes {\hat {e}}_{i}=-u_{j}{\hat {e}}_{k}\otimes {\hat {e}}_{k}\times {\hat {e}}_{j}=-\mathbf {1} \times {\vec {u}}$

Darin ist ϵ_ijk = (ê_i × ê_j) · ê_k das Permutationssymbol. Mit $\mathbf {W} =-\mathbf {1} \times {\vec {u}}$ folgt die dritte Form des Lemmas:

$\operatorname {rot} (\mathbf {T} )=\mathbf {0} \;{\text{und}}\;\operatorname {Sp} (\mathbf {T} )=0\;\rightarrow \quad \exists \mathbf {W} \colon \mathbf {T} =\operatorname {rot} \mathbf {W} \;{\text{mit}}\;\mathbf {W} =-\mathbf {W} ^{\top }$		(III)

oder mit $\mathbf {W} =\mathbf {1} \times {\vec {u}}$ und dem Nabla-Operator

$\nabla \times (\mathbf {T} ^{\top })=\mathbf {0} \;{\text{und}}\;\operatorname {Sp} (\mathbf {T} )=0\;\rightarrow \quad \exists \mathbf {W} \colon \mathbf {T} =\nabla \times \mathbf {W} \;{\text{mit}}\;\mathbf {W} =-\mathbf {W} ^{\top }$		(III)

Literatur

• Otto Forster: Analysis. Band 3: Integralrechnung im ℝⁿ mit Anwendungen. 4. Auflage. Vieweg + Teubner, Braunschweig u. a. 2007, ISBN 978-3-528-37252-1.
• John M. Lee: Introduction to Smooth Manifolds (= Graduate Texts in Mathematics 218). Springer-Verlag, New York NY u. a. 2003, ISBN 0-387-95448-1.
• C. Truesdell: Festkörpermechanik II in S. Flügge (Hrsg.): Handbuch der Physik, Band VIa/2. Springer-Verlag, 1972, ISBN 3-540-05535-5, ISBN 0-387-05535-5.