Refine
Has Fulltext
- yes (2)
Is part of the Bibliography
- yes (2)
Document Type
- Doctoral Thesis (2)
Language
- German (2) (remove)
Keywords
- Blazar (1)
- Computersimulation (1)
- Mathematisches Modell (1)
- Schockfront (1)
- Sonnenwind (1)
- Stoßfreies Plasma (1)
- Stoßwelle (1)
- Strahlung (1)
- Teilchenbeschleunigung (1)
- kinetische Plasmaphysik (1)
Die vorliegende Arbeit beschäftigt sich mit den Prozessen, die in einer Unterklasse der Aktiven Galaxienkerne, den Blazaren, das Emissionsspektrum dieser Objekte erzeugen. Dies beinhaltet insbesondere den Beschleunigungsprozess, der eine nichtthermische Teilchenverteilung erzeugt, sowie diverse Strahlungsprozesse. Das Spektrum dieser Quellen reicht dabei vom Radiobereich bis zu Energien im TeV-Bereich. Die Form des zeitlich gemittelten Spektrums kann durch Modelle bereits sehr gut beschrieben werden. Insbesondere die erste der beiden dominierenden Komponenten des Spektrums kann mit hoher Sicherheit mit Synchrotronemission einer Elektronenenergieverteilung in Form eines Potenzgesetzes identifiziert werden. Für den Ursprung der zweiten Komponente existieren jedoch verschiedene Erklärungsversuche. Dies sind im wesentlichen die inverse Compton-Streuung der internen oder externer Strahlung (leptonische Modelle) sowie die Emission und photohadronische Wechselwirkung einer hochenergetischen Verteilung von Protonen in der Quelle.
Eine räumliche Auflösung des Ursprungs der detektierten Strahlung ist mit den zur Verfügung stehenden Teleskopen nicht möglich. Einschränkungen für die Ausdehnung dieser Emissionszone ergeben sich lediglich aus der Variation des Emissionsspektrums. Eine Bestimmung der Morphologie ist jedoch im selbstabsorbierten Radiobereich des Spektrums durch die Ausnutzung von interferometrischen Beobachtungen möglich. Die resultierenden Längen, auf denen die im inneren der Quelle selbstabsorbierte Strahlung die Quelle schließlich verlässt, sind jedoch etwa zwei Größenordnungen oberhalb der aus den Variabilitätszeitskalen gefolgerten Limits.
Das im Rahmen dieser Arbeit entwickelte Modell soll dabei helfen, verschiedene Beobachtungen mit Hilfe eines quantitativen Modells zu beschreiben. Hier steht insbesondere die Korrelation zwischen den Verläufen der Hochenergie- und Radioemission im Vordergrund. Eine Aussage über die Existenz einer solchen Verbindung konnte aus den bisherigen Beobachtungen nicht getroffen werden.
Eine quantitative Modellierung könnte bei der Interpretation der bisher uneindeutigen Datenlage helfen. Eine weitere, durch Modelle bisher nicht beschreibbare, Beobachtungsevidenz sind extrem kurzzeitige Variationen des Flusszustands.
Die Lichtlaufzeit durch das für die Modellierung benötigte Raumgebiet ist zumeist größer als die beobachtete Zeitskala.
Zudem deuten die Beobachtungen darauf hin, dass manche dieser Flussausbrüche nicht zwischen den verschiedenen Bändern korreliert sind, wie es zumindest die leptonischen Modelle erwarten lassen würden.
Das hier beschriebene Modell verbindet eine räumliche Auflösung des Emissionsgebiets mit dem dominanten Beschleunigungsmechanismus. Hierdurch konnte zunächst gezeigt werden, dass die Beschreibung von Variabilität auch auf Skalen unterhalb der Lichtlaufzeit durch das modellierte Raumgebiet möglich ist. Zudem wurde ein Szenario quantifiziert, dass im leptonischen Fall unkorrelierte Ausbrüche vorhersagt.
\thispagestyle{empty}
Durch eine Erweiterung des Emissionsgebiets gegenüber anderen Blazar-Modellen um zwei Größenordnung konnte zudem eine Verknüpfung zwischen dem Hochenergie- und dem Radiobereich erfolgen. Die gefundene Morphologie des Einschlussgebiets der nichtthermischen Teilchenpopulation beinhaltet eine physikalisch sinnvolle Randbedingung für das Emissionsgebiet der Hochenergiestrahlung, die zudem den für die betrachtete Quelle korrekten Spektralindex im Radiobereich erzeugt.
Darüber hinaus wurden in das Modell sowohl leptonische als auch hadronische Prozesse integriert, die eine flexible und unvoreingenommene Modellierung potentieller Hybridquellen erlauben.
Mit dem entwickelten Modell ist es möglich, aus detailliert vermessenen Lichtkurven im Hochenergiebereich die zu erwartende Radioemission vorherzusagen. Die in diese Vorhersage eingehenden Parameter lassen sich aus der Modellierung des Gleichgewichtsspektrums bestimmen.
Das Magnetfeld der Sonne ist kein einfaches statisches Dipolfeld, sondern weist
wesentlich kompliziertere Strukturen auf. Wenn Rekonnexion die Topologie eines
Feldlinienbündels verändert, wird viel Energie frei, die zuvor im Magnetfeld
gespeichert war. Das abgetrennte Bündel wird mit dem damit verbundenen Plasma
mit großer Geschwindigkeit durch die Korona
von der Sonne weg bewegen. Dieser Vorgang wird als koronaler Massenauswurf
bezeichnet. Da diese Bewegung mit Geschwindigkeiten deutlich über der
Alfv\'en-Geschwindigkeit, der kritischen Geschwindigkeit im Sonnenwind,
erfolgen kann, bildet sich eine Schockfront, die durch den Sonnenwind
propagiert.
Satelliten, die die Bedingungen im Sonnenwind beobachten, detektieren beim
Auftreten solcher Schockfronten einen erhöhten Fluss von hochenergetischen
Teilchen. Mit Radioinstrumenten empfängt man zeitgleich elektromagnetische
Phänomene, die als Radiobursts bezeichnet werden, und ebenfalls für die
Anwesenheit energiereicher Teilchen sprechen. Daher, und aufgrund von
theoretischen Überlegungen liegt es nahe, anzunehmen, daß Teilchen an der
Schockfront beschleunigt werden können.
Die Untersuchung der Teilchenbeschleunigung an kollisionsfreien Schockfronten
ist aber noch aus einem zweiten Grund interessant. Die Erde wird kontinuierlich
von hochenergetischen Teilchen, die aus historischen Gründen als kosmische
Strahlung bezeichnet werden, erreicht. Die gängige Theorie für deren Herkunft
besagt, daß zumindest der galaktische Anteil durch die Beschleunigung an
Schockfronten, die durch Supernovae ausgelöst wurden, bis zu den beobachteten
hohen Energien gelangt sind. Das Problem bei der Untersuchung der Herkunft der
kosmischen Strahlung ist jedoch, daß die Schockfronten um Supernovaüberreste
aufgrund der großen Entfernung nicht direkt beobachtbar sind.
Es liegt dementsprechend nahe, die Schockbeschleunigung an den wesentlich
näheren und besser zu beobachtenden Schocks im Sonnensystem zu studieren, um so
Modelle und Simulationen entwickeln und testen zu können.
Die vorliegende Arbeit beschäftigt sich daher mit Simulationen von
Schockfronten mit Parametern, die etwa denen von CME getriebenen Schocks
entsprechen. Um die Entwicklung der Energieverteilung der Teilchen zu studieren,
ist ein kinetischer Ansatz nötig. Dementsprechend wurden die Simulationen mit
einem Particle-in-Cell Code durchgeführt. Die Herausforderung ist dabei die
große Spanne zwischen den mikrophysikalischen Zeit- und Längenskalen, die aus
Gründen der Genauigkeit und numerischen Stabilität aufgelöst werden müssen und
den wesentlich größeren Skalen, die die Schockfront umfasst und auf der
Teilchenbeschleunigung stattfindet.
Um die Stabilität und physikalische Aussagekraft der Simulationen
sicherzustellen, werden die numerischen Bausteine mittels Testfällen, deren
Verhalten bekannt ist, gründlich auf ihre Tauglichkeit und korrekte
Implementierung geprüft.
Bei den resultierenden Simulationen wird das Zutreffen von analytischen
Vorhersagen (etwa die Einhaltung der Sprungbedingungen) überprüft. Auch die
Vorhersagen einfacherer Plasmamodelle, etwa für das elektrostatischen
Potential an der Schockfront, das man auch aus einer Zwei-Fluid-Beschreibung
erhalten kann, folgen automatisch aus der selbstkonsistenten, kinetischen
Beschreibung. Zusätzlich erhält man Aussagen über das Spektrum und die Bahnen
der beschleunigten Teilchen.