Benutzer:Tepask/Sandkiste HM

Die magnetische Helizität ist eine Größe in der Magnetohydrodynamik. Sie ist ein Maß, das zeigt, wie sehr die Feldlinien ineinander verketten und umeinander kreisen. Bei nicht vorhandenem elektrischem Widerstand im System ist sie eine Erhaltungsgröße, sie ist also eine "ideale Invariante". Wenn magnetische Helizität in einem Magnetfeld vorhanden ist, tendiert das System dazu von kleinskaligen magnetischen Strukturen immer größere Strukturen zu formen. Dies kann man auch als einen "inversen Transfer im Fourierraum" bezeichnen.

Diese letzte Eigenschaft macht die magnetische Helizität sehr besonders: turbulente dreidimensionale Strömungen, die in der Astrophysik besonders oft auftreten, tendieren dazu, Strukturen zu zerstören, indem große Wirbel in immer kleinere zerfallen (ein als “(direkte) Energiekaskade” bezeichnetes Phänomen, das durch Lewis Fry Richardson und Andrei Nikolajewitsch Kolmogorov beschrieben worden ist), wo sie letzten Endes durch Viskosität in Hitze verwandelt werden. Im Gegensatz dazu geschieht durch eine Art inverse magnetische Helizitätskaskade das Gegenteil: kleine helikale magnetische Strukturen (das heißt, Strukturen mit einer Helizität ungleich Null) führen zu großen Magnetfeldern.

Dies ist zum Beispiel durch die Heliosphärische Stromschicht - eine große magnetische Struktur in unserem Sonnensystem - sichtbar. Magnetische Helizität ist allgemein in vielen astrophysikalischen Systemen von großer Relevanz, da in diesen der elektrische Widerstand typischerweise sehr gering ist. Beispielsweise ist die Dynamik der magnetischen Helizität in Sonneneruptionen und koronale Massenauswürfen von Bedeutung. Der Sonnenwind enthält magnetische Helizität. Die Erhaltung der magnetischen Helizität ist auch in Dynamo-Prozessen von großer Relevanz. In der Fusionsforschung spielt magnetische Helizität beispielsweise in Reversed field pinch Experimenten eine wichtige Rolle.

Die magnetische Helizität eines Systems kann man jedoch nicht direkt messen. Unter bestimmten Annahmen und Bedingungen, kann man sie trotzdem aus der "Stromhelizität" schlussfolgern (siehe diesen Abschnitt).

Mathematische Formulierung

Ganz allgemein wird die Helizität $H^{\mathbf {f} }$ eines glatten dreidimensionallen Vektorfeldes $\mathbf {f}$ durch das Volumenintegral des Skalarproduktes von $\mathbf {f}$ und dessen Rotation $\nabla \times {\mathbf {f} }$definiert:

$H^{\mathbf {f} }=\int {\mathbf {f} }\cdot {\nabla \times {\mathbf {f} }}\,d^{3}{\mathbf {r} }$,

wobei $d^{3}{\mathbf {r} }$ das infinitesimale Volumenelement ist, und die Integration über das gesamte betrachtete Gebiet stattfindet.

Die magnetische Helizität $H^{\mathbf {M} }$ wird als die Helizität des magnetischen Vektorpotenzials ${\mathbf {A} }$ definiert, wobei ${\mathbf {B} }=\nabla \times {\mathbf {A} }$ dem Magnetfeld entspricht:

$H^{\mathbf {M} }=\int {\mathbf {A} }\cdot {\mathbf {B} }\,d^{3}{\mathbf {r} }$.

Die magnetische Helizität $\int {\mathbf {A} }\cdot {\mathbf {B} }\,d^{3}{\mathbf {r} }$ soll nicht mit der Helizität des Magnetfeldes $H^{\mathbf {J} }=\int {\mathbf {B} }\cdot {\mathbf {J} }\,d^{3}{\mathbf {r} }$, mit ${\mathbf {J} }=\nabla \times {\mathbf {B} }$ der elektrische Stromdichte, verwechselt werden. Diese Quantität wird Stromhelizität genannt (auf Englisch "current helicity"). Im Gegensatz zu der magnetischen Helizität ist die Stromhelizität keine ideale Invariante (selbst wenn der elektrische Widerstand Null ist, wird sie nicht erhalten).

Da das magnetische Vektorpotenzial nicht eichinvariant ist, ist die magnetische Helizität auch im allgemeinen Fall nicht eichinvariant (siehe diesen Abschnitt). Als Folge davon, kann man die magnetische Helizität nicht direkt messen. Man kann jedoch die Stromhelizität in einem System unter bestimmte Bedingungen und Annahmen bestimmen und daraus, unter weiteren Bedingungen und Annahmen, die magnetische Helizität schlussfolgern.

Topologische Interpretation

Abbildung 1: topologisch gesehen können Verwringung (auf Englisch "writhe") und Drehung (eng. "twist") ineinander umgewandelt werden. Dieses Bild wurde vom Wikipedia Artikel über "DNA supercoil" wiederbenutzt, wo andere topologische Aspekte detaillierter beschrieben werden.

Der Name "Helizität" beruht auf der Tatsache, dass in einer Strömung mit der Geschwindigkeit $\mathbf {v}$ und Wirbelstärke $\mathbf {\omega } =\nabla \times \mathbf {v}$, die Bewegung der Fluidpartikel in Gebieten in denen die kinetische Helizität $H^{K}=\int \mathbf {v} \cdot \mathbf {\omega } \neq 0$ eine Helix bildet. Für $H^{K}>0$ ist diese linkshändig, und für $H^{K}<0$ rechtshändig. Auf ähnliche Weise verhält es sich für die Form von Magnetfeldlinien.

Gebiete, in denen magnetische Helizität ungleich Null ist können zwar helikale Magnetfeldlinien enthalten, aber auch andere Arten von magnetischen Strukturen aufweisen. Magnetische Helizität ist nämlich eine Verallgemeinerung vom topologischen Konzept der Verknüpfungszahl (auf Englisch "linking number"). Die Verknüpfungszahl beschreibt wie sehr die Magnetfeldlinien einander verketten (für einen mathematischen Beweis dieser Eigenschaft, siehe ). Durch ein einfaches Experiment mittels einer Schere und einem Blatt Papier kann man zeigen, dass Magnetfeldinien, die umeinander kreisen, auch als verkettet angesehen werden können (siehe Abbild 5 in der Referenz). Deswegen kann man magnetische Helizität sowohl als helikalformige Magnetfeldlinien, verkettete magnetische Strukturen, aber auch als ein umeinander Kreisen der Magnetfeldlinien betrachten.

Umeinander kreisende Magnetfeldstrukturen können auch unterschiedliche Formen haben. Zum Beispiel betrachten wir eine gedrehte Ansammlung von benachbarten Magnetfeldlinien, die sogenannte geschlossene "magnetische Flußrohre" bilden (zur Veranschaulichung, siehe zum Beispiel die Rohre in Abbildung 1). Die Drehung bedeutet hier, dass sich die Rohre um ihre eigene Axe drehen (Figuren mit Twist=$\pm 1$). Topologisch gesehen, können Drehung und Verwringung ineinander umgewandelt werden, wobei die Verwringung bedeutet, dass die Achse der Rohre selbst Drehungen macht (Figuren mit Writhe=$\pm 1$). Man kann ebenfalls zeigen, dass Knoten auch äquivalent zu Drehungen und/oder Verwringungen betrachtet werden können.

Bildung größeren magnetischen Strukturen

Magnetische Helizität führt dazu, dass sich von kleinen helikalen Strukturen immer größere Strukturen ausbilden. Dies nennt man einen inversen Transfer im Fourierraum (im Gegensatz zu der (direkten) Energiekaskade in 3D-turbulente Strömungen). Die Möglichkeit von solch einem inversen Transfer wurde zuerst bei Uriel Frisch und dessen Mitarbeitern entdeckt und wurde durch mehrere numerische Experimente bestätigt. Deshalb ist magnetische Helizität ein potenzieller Kandidat, um die Existenz und die Erhaltung großer magnetischer Strukturen im Weltraum zu begründen.

Ein Argument für diesen inversen Transfer von Referenz ist im folgenden wiederholt. Dieser basiert auf der sogenannte "Realisierbarkeitsbedingung" ("realizability condition" auf Englisch) von der magnetischen Helizität im Fourierspektrum ${\hat {H}}_{\mathbf {k} }^{M}={\hat {\mathbf {A} }}_{\mathbf {k} }^{*}\cdot {\hat {\mathbf {B} }}_{\mathbf {k} }$ (wo ${\hat {\mathbf {B} }}_{\mathbf {k} }$ der Fourierkoeffizient vom Wellenvektor $\mathbf {k}$ des Magnetfeldes ${\mathbf {B} }$ ist, und in ähnlicher Weise für ${\hat {\mathbf {A} }}$). Der Stern * bedeutet das komplexe Konjugat. Die "Realisierbarkeitsbedingung" entspricht der Cauchy-Schwarzsche Ungleichung, angewandt an ${\hat {H}}_{\mathbf {k} }^{M}$, und lautet:

$|{\hat {H}}_{\mathbf {k} }^{M}|\leq {\frac {2{\hat {E}}_{\mathbf {k} }^{M}}{|{\mathbf {k} }|}}$,

mit $E_{\mathbf {k} }^{M}={\frac {1}{2}}{\hat {\mathbf {B} }}_{\mathbf {k} }^{*}\cdot {\hat {\mathbf {B} }}_{\mathbf {k} }$ der Spektrum der magnetischen Energie. Die Eigenschaft $|{\hat {\mathbf {B} }}_{\mathbf {k} }|=|{\mathbf {k} }||{\hat {\mathbf {A} }}_{\mathbf {k} }^{\perp }|$ (mit ${\hat {\mathbf {A} }}_{\mathbf {k} }^{\perp }$ der divergenzfreie Teil des Fourierspektrums des magnetischen Vektorpotenzialfeldes, orthogonal zum Wellenvektor $\mathbf {k}$ im Fourierraum) wurde benutzt, um diese Ungleichung zu erhalten, da ${\hat {\mathbf {B} }}_{\mathbf {k} }=i{\mathbf {k} }\times {\hat {\mathbf {A} }}_{\mathbf {k} }$ im Fourierraum (was ${\mathbf {B} }=\nabla \times {\mathbf {A} }$ entspricht). Im Artikel fehlt der Faktor 2, da in diesem eine alternative Definition der magnetischen Helizität, und zwar ${\frac {1}{2}}\int {\mathbf {A} }\cdot {\mathbf {B} }d^{3}{\mathbf {r} }$, benutzt wird.

Stellen wir uns eine Ausgangssituation ohne Geschwindigkeitsfeld und mit einem Magnetfeld vor, das aus einer Überlagerung von zwei Moden besteht, und zwar an die Wellenvektoren $\mathbf {p}$ und $\mathbf {q}$. Wir nehmen an, dass das Magnetfeld völlig helikal ist, das heißt, dass es die Realisierbarkeitsbedingung sättigt: $|{\hat {H}}_{\mathbf {p} }^{M}|={\frac {2E_{\mathbf {p} }^{M}}{|{\mathbf {p} }|}}$ und $|{\hat {H}}_{\mathbf {q} }^{M}|={\frac {2E_{\mathbf {q} }^{M}}{|{\mathbf {q} }|}}$. Angenommen, dass alle Energie- und magnetischen Helizitätstransfers zu einem dritten Wellenvektor $\mathbf {k}$ geschehen, liefern die Erhaltung der magnetischen Helizität und die Erhaltung der Gesamtenergie ${\hat {E}}^{G}={\hat {E}}^{M}+{\hat {E}}^{K}$ (die Summe von der (m)agnetischen und der (k)inetischen Energien) die Gleichungen:

${\hat {H}}_{\mathbf {k} }^{M}={\hat {H}}_{\mathbf {p} }^{M}+{\hat {H}}_{\mathbf {q} }^{M},$

${\hat {E}}_{\mathbf {k} }^{G}={\hat {E}}_{\mathbf {p} }^{G}+{\hat {E}}_{\mathbf {q} }^{G}={\hat {E}}_{\mathbf {p} }^{M}+{\hat {E}}_{\mathbf {q} }^{M}.$

Die zweite Gleichung für die Gesamtenergie beruht auf die Tatsache, dass wir eine Anfangsbedingung ohne Geschwindigkeitsfeld betrachten, sodass ${\hat {E}}_{\mathbf {p} }^{G}={\hat {E}}_{\mathbf {p} }^{M}$ und ${\hat {E}}_{\mathbf {q} }^{G}={\hat {E}}_{\mathbf {q} }^{M}$. Dann gilt zwangsläufig $|\mathbf {k} |\leq \max(|\mathbf {p} |,|\mathbf {q} |)$. Angenommen wir hätten $|\mathbf {k} |>\max(|\mathbf {p} |,|\mathbf {q} |)$, dann:

${\hat {H}}_{\mathbf {k} }^{M}={\hat {H}}_{\mathbf {p} }^{M}+{\hat {H}}_{\mathbf {q} }^{M}={\frac {2{\hat {E}}_{\mathbf {p} }^{M}}{|\mathbf {p} |}}+{\frac {2{\hat {E}}_{\mathbf {q} }^{M}}{|\mathbf {q} |}}>{\frac {2({\hat {E}}_{\mathbf {p} }^{M}+{\hat {E}}_{\mathbf {q} }^{M})}{|\mathbf {k} |}}={\frac {2{\hat {E}}_{\mathbf {k} }^{G}}{|\mathbf {k} |}}\geq {\frac {2{\hat {E}}_{\mathbf {k} }^{M}}{|\mathbf {k} |}},$

und dies würde die Realisierbarkeitsbedingung brechen. Dies impliziert, dass $|\mathbf {k} |\leq \max(|\mathbf {p} |,|\mathbf {q} |)$. Insbesondere, für $|{\mathbf {p} }|=|{\mathbf {q} }|$, wird die magnetische Helizität zu einem kleineren Wellenvektor transferiert, das heißt, zu größeren Skalen.

Ideale Invarianz

In den späten 1950er Jahren entdeckten Lodewijk Woltjer und Walter M. Elsasser unabhängig voneinander, dass die magnetische Helizität eine ideale Invariante ist. Dies bedeutet, dass sie in Systemen, die keinen elektrischen Widerstand aufweisen, eine Erhaltungsgröße ist. Systeme mit sehr geringem elektrischem Widerstand kann man oft in astrophysikalischen Systemen finden. Da die magnetische Helizität topologische Eigenschaften des Magnetfeldes beschreibt, führt diese Invarianz zu großen Beschränkungen in der Gesamttopologie des Magnetfeldes: wenn an einer Skala positive Helizität auftritt, muss auf einer anderen Skala im geschlossenen System negative Helizität auftreten.

Im folgenden wird Woltjers Beweis, gültig für ein geschlossenes System, wiedergegeben:

Die Gleichung zur zeitlichen Entwicklung des Magnetfeldes und des magnetischen Vektorpotenzials lauten im Fall eines idealen magnetohydrodynamischen Systems:

$\partial _{t}{\mathbf {B} }=\nabla \times ({\mathbf {v} }\times {\mathbf {B} }),$

$\partial _{t}{\mathbf {A} }={\mathbf {v} }\times {\mathbf {B} }+\nabla (\Phi )$.

Die zweite Gleichung erhält man bei der „Entrotationierung“ (manchmal wird der Begriff „uncurling“ auf Englisch benutzt) der ersten Gleichung. Der Term $\Phi$ ist ein Skalarpotential, definiert durch die gewählte Eichung (siehe den Abschnitt über die Eichung). Wenn die Eichung so gewählt wird, dass das entsprechende Gradientenfeld verschwindet ($\nabla (\Phi )$=0), wird die Zeitentwicklung der magnetischen Helizität durch die folgende Gleichung beschrieben:

$\partial _{t}H^{\mathbf {M} }=\int (\partial _{t}{\mathbf {A} })\cdot {\mathbf {B} }+{\mathbf {A} }\cdot \partial _{t}{\mathbf {B} }d^{3}{\mathbf {r} }=\int ({\mathbf {v} }\times {\mathbf {B} })\cdot {\mathbf {B} }d^{3}{\mathbf {r} }+\int {\mathbf {A} }\cdot (\nabla \times \partial _{t}{\mathbf {A} })d^{3}{\mathbf {r} }$.

Das erste Integral ist Null, da ${\mathbf {B} }$ orthogonal zum Kreuzprodukt ${\mathbf {v} }\times {\mathbf {B} }$ ist. Das zweite Integral kann partiell integriert werden, sodass:

$\partial _{t}H^{\mathbf {M} }=\int _{V}\nabla \times {\mathbf {A} }\cdot (\partial _{t}{\mathbf {A} })d^{3}{\mathbf {r} }+\int _{S}{\mathbf {A} }\times (\partial _{t}{\mathbf {A} })d^{2}{\mathbf {r} }=0.$

Die erste Integration geschieht über das Gesamtvolumen und ist Null, da $\nabla \times {\mathbf {A} }={\mathbf {B} }\perp \partial _{t}{\mathbf {A} }$, wie oben bereits erwähnt. Das zweite Integral ist das Oberflächenintegral über die Grenzen des geschlossenen Systems $S$. Dieses Integral ist Null, da die Bewegungen im geschlossenen System nicht das magnetische Vektorpotenzial außerhalb beeinflussen können, sodass auf der Systemgrenze $\partial _{t}A=0$, da das magnetische Vektorpotenzial eine kontinuierliche Funktion ist.

Wenn die magnetische Helizität eichinvariant ist, bleibt deshalb die magnetische Helizität selbst dann ideal erhalten, wenn die besondere Eichung $\nabla (\Phi )=0$ nicht erfüllt ist.

Die magnetische Helizität bleibt auch in Situationen mit einem kleinen, endlichen elektrischen Widerstand näherungsweise erhalten. Hierbei tritt magnetische Rekonnexion auf, wodurch magnetische Energie in kinetische Energie umgewandelt wird, jedoch die magnetische Helizität ungefähr erhalten bleibt.

Über die Eichung

Da das magnetische Vektorpotenzial $\mathbf {A}$ nicht eichinvariant ist (es kann umdefiniert werden, indem man einen Gradienten addiert, ohne das Magnetfeld zu ändern. Die Rotation eines Gradientenfeldes ist nämlich Null), ist auch die magnetische Helizität in der Regel nicht eichinvariant. Für perfekt-leitende Systemgrenzen oder periodische Systeme ohne magnetischen Fluss durch die Grenzen des physikalischen Systems ist die magnetische Helizität jedoch eichinvariant. Das heißt, sie ist dann unabhängig von der gewählte Eichung. Periodische Systeme ohne magnetischen Fluss durch die Grenzen des Systems werden deshalb oft bei numerischen Simulationen für die Untersuchung der Dynamik der magnetischen Helizität benutzt (zum Beispiel).

Im allgemeinen Fall kann man aber eine sogenannte eichinvariante "relative Helizität" definieren. "Relativ" bedeutet hier im Vergleich zu einem "gut gewählten" magnetischen Vektorpotenzial und dem entsprechenden Magnetfeld, siehe Literatur. Dies ist dann besonders wichtig, wenn man realistischere Systeme betrachtet: zum Beispiel sind astrophysikalische Systeme in Wirklichkeit nicht periodisch und meistens nicht geschlossen.