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

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

Wegen <math>\mathrm{d} \circ \mathrm{d} \equiv 0</math> 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 <math>U \subseteq \R^d</math> definierte geschlossene Differentialform exakt ist.

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

<math> \mathrm{H}_{\mathrm{dR}}^k(U) = 0 </math>

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 <math>\Phi(\mathbf r)</math> (<math>k=1</math>), ein quellfreies Vektorfeld auf einem konvexen Gebiet durch Rotation eines Vektorpotentials <math>\vec A(\mathbf r, t)</math> (<math>k=2</math>), und eine skalare Felddichte („Quellendichte“) als Divergenz eines Vektorfeldes (<math>k=3</math>) dargestellt werden können.

Beweis (konstruktiv)

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

<math>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}</math>.

(Das Dachsymbol in der <math>i_\alpha</math>-ten Spalte der rechten Seite bedeutet, dass das entsprechende Differential ausgelassen wird.)

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

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

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

Die angegebene Identität verallgemeinert zugleich das Poincarésche Lemma durch Zerlegung einer beliebigen Differentialform <math>\omega</math> 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 <math>P</math> eine kontrahierende Homotopie, die z. B. auf den zentralen Punkt des hier betrachteten sternförmigen Gebietes kontrahiert.

Umeichung

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

Der Zusatz <math>+\, \mathrm d \xi^{k-2}</math> wird auch als Eichtransformation bzw. Umeichung von <math>\eta_1^{k-1}</math> bezeichnet.

Anwendung in der Elektrodynamik

Aus der Elektrodynamik ist der Fall eines von einem stationären Strom erzeugten Magnetfeldes bekannt, mit dem sog. Vektorpotential <math>\vec A(\mathbf r )\,.</math> Dieser Fall entspricht <math>k=2</math>, wobei das sternförmige Gebiet der <math>\mathbb R^3</math> ist. Der Vektor der Stromdichte ist <math>\vec j</math> und entspricht der Stromform <math>\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\,.</math> Für das Magnetfeld <math>\vec B</math> gilt Analoges: es entspricht der Magnetflussform <math>\Phi_B:=B_1{\rm d}x_2\wedge {\rm d}x_3 +\dots </math> und lässt sich aus dem Vektorpotential ableiten: <math>\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</math>, oder <math>\Phi_B={\rm d}\mathbf A</math>. Dabei entspricht das Vektorpotential <math>\vec A</math> der Potentialform <math>\mathbf A:=A_1{\rm d}x_1+A_2{\rm d}x_2+A_3{\rm d}x_3\,.</math> Die Geschlossenheit der Magnetflussform entspricht der Quellenfreiheit des Magnetfeldes <math>(\operatorname{div} \vec B \equiv 0\,).</math>

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

<math>A_i(\vec r)

=\int \frac{\mu_0 j_i(\vec r^{\,'})\,\, dx_1'dx_2'dx_3'}{4\pi |\vec r -\vec r^{\,'}|}\,,</math>

dabei ist <math>\mu_0</math> 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 <math>\vec E</math> entspricht, dem Coulombpotential <math>\,\phi (x_1,x_2, x_3)</math> einer gegebenen Ladungsverteilung mit der Dichte <math>\rho (x_1,x_2,x_3)</math>. Man vermutet an dieser Stelle bereits, dass

<math>\vec E</math> und <math>\vec B</math> bzw.
<math>\rho </math> und <math>\vec j</math> sowie
<math>\,\phi</math> und <math>\vec A</math>

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 <math>A_i</math> zu den Raumkoordinaten das Zeitargument <math>t</math> hinzugefügt werden, während auf der rechten Seite in <math>j_i'</math> die sog. „retardierte Zeit“ <math>t':=t-\tfrac{|\vec r -\vec r^{\,'}|}{c}</math> zu ergänzen ist. Es wird dabei wie zuvor über die drei Raumkoordinaten <math> \vec r^{\,'}</math> integriert. Schließlich ist <math>c</math> 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:

:= margin-left: 1.5em;	::= margin-left: 3em;	:::= margin-left: 4.5em;	border:1px solid #BBBBBB; padding:2px;\| border-collapse:collapse;}}"
Formel		{{#if:\|{{{3}}}\|(3.5)}}
Formel		{{#if:\|{{{3}}}\|(3.5)}}

{\partial x_k}

=\vec{0} \quad\rightarrow\quad \exists\varphi\colon\vec{u}=\operatorname{grad}\varphi </math>|I}}

Der Operator „grad“ bildet den Gradient, die Vektoren <math>\hat{e}_{1,2,3}</math> sind die Standardbasis des kartesischen Koordinatensystems mit Koordinaten <math>x_{1,2,3}</math> und es wurde die einsteinsche Summenkonvention angewendet, der gemäß ü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 <math>\mathbf{T}=\hat{e}_i\otimes\vec{t}_i</math>, dessen Zeilenvektoren <math>\vec{t}_{1,2,3}</math> mit dem dyadischen Produkt „⊗“ zum Tensor zusammengefügt werden. Jeder Tensor zweiter Stufe kann in dieser Form dargestellt werden. Die Rotation des Tensors verschwinde

<math>\begin{align}

\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{align}</math>

so dass also jeder Zeilenvektor rotationsfrei ist. Dann gibt es für jeden Zeilenvektor ein Skalarfeld <math>u_i</math>, dessen Gradient er ist:

<math>

\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}\,, </math>

denn der Gradient des Vektors <math>\vec u:=u_i\hat{e}_i</math> bildet sich gemäß:

<math>

\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\,. </math>

Damit gilt die zweite Form des Lemmas:

:= margin-left: 1.5em;	::= margin-left: 3em;	:::= margin-left: 4.5em;	border:1px solid #BBBBBB; padding:2px;\| border-collapse:collapse;}}"
II		{{#if:\|{{{3}}}\|(3.5)}}
II		{{#if:\|{{{3}}}\|(3.5)}}

Wenn zusätzlich die Spur des Tensors verschwindet, dann ist das Vektorfeld divergenzfrei:

<math>

\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\,. </math>

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

<math>\begin{align}

\operatorname{rot}(\mathbf1\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{align}</math>

und der Tensor <math>\mathbf1\times\vec{u}</math> ist schiefsymmetrisch:

<math>

(\mathbf1\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 = -\mathbf1\times\vec{u} </math>

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

:= margin-left: 1.5em;	::= margin-left: 3em;	:::= margin-left: 4.5em;	border:1px solid #BBBBBB; padding:2px;\| border-collapse:collapse;}}"
III		{{#if:\|{{{3}}}\|(3.5)}}
III		{{#if:\|{{{3}}}\|(3.5)}}

oder mit <math>\mathbf{W}=\mathbf1\times\vec{u}</math> und dem Nabla-Operator

:= margin-left: 1.5em;	::= margin-left: 3em;	:::= margin-left: 4.5em;	border:1px solid #BBBBBB; padding:2px;\| border-collapse:collapse;}}"
III		{{#if:\|{{{3}}}\|(3.5)}}
III		{{#if:\|{{{3}}}\|(3.5)}}

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.