Refine
Has Fulltext
- yes (3) (remove)
Is part of the Bibliography
- yes (3)
Document Type
- Doctoral Thesis (3)
Language
- German (3) (remove)
Keywords
- Astrophysik (3) (remove)
Institute
Die vorliegende Arbeit beschäftigt sich mit der Abstrahlung von Aktiven Galaxienkernen. Das erste Maximum der charakteristischen Doppelpeakstruktur des $\nu F_{\nu}$-Spektrums vom Blazaren ist zweifelsfrei Synchrotronstrahlung hochenergetischer Elektronen innerhalb des relativistischen Ausflusses des zugrundeliegenden Aktiven Galaxienkerns. Die zum zweiten (hochenergetischen) Maximum beitragenden Strahlungsprozesse und Teilchenspezies hingegen sind Gegenstand aktueller Diskussionen. In dieser Arbeit wir ein vollständig selbstkonsistentes und zeitabhängiges hybrides Emissionsmodell, welches auch Teilchenbeschleunigung berücksichtigt, entwickelt und auf verschiedene Blazar-Typen entlang der Blazar-Sequenz, von BL Lac Objekten mit verschiedenen Peakfrequenzen bis hin zu Flachspektrum-Radioquasaren, angewendet. Die spektrale Emission ersterer kann gut im rein leptonischen Grenzfall, d.h. der zweite $\nu F_{\nu}$-Peak kommt durch invers Compton-gestreute Synchrotronphotonen der abstrahlenden Elektronen selbst zustande, beschrieben werden. Zur Beschreibung letzterer muss man nicht-thermische Protonen innerhalb des Jets zulassen um die Dominanz des zweiten Maximums im Spektrum konsistent zu erklären. In diesem Fall besteht der zweite Peak aus Protonensynchrotronstrahlung und Kaskadenstrahlung der photohadronischen Prozesse. Mit dem entwickelten Modell ist es möglich auch die zeitliche Information, welche durch Ausbrüche von Blazaren bereitgestellt wird, auszunutzen um zum einen die freien Modellparameter weiter einzuschränken und -viel wichtiger- zum anderen leptonisch dominierte Blazare von hadronischen zu unterscheiden. Hierzu werden die typischen Zeitunterschiede in den Interbandlichtkurven als hadronischer Fingerabdruck benutzt.\\ Mit einer Stichprobe von 16 Spektren von zehn Blazaren entlang der Blazar-Sequenz, welche in unterschiedlichen Flusszuständen und mit starker Variabilität beobachtet wurden, ist es möglich die wichtigsten offenen Fragen der Physik relativistischer Ausbrüche in systematischer Art und Weise zu adressieren. Anhand der modellierten Ausbrüche kann man erkennen, dass sechs Quellen rein leptonisch dominiert sind, aber vier Protonen bis auf $\gamma \approx 10^{11}$ beschleunigen, was Auswirkungen auf die möglichen Quellen extragalaktischer kosmischer Strahlung unter den Blazaren hat. Darüber hinaus findet sich eine Abhängigkeit zwischen dem Magnetfeld der Emissionsregion und der injizierten Leuchtkraft, welche unabhängig von den zugrunde liegenden Teilchenpopulationen Gültigkeit besitzt. In diesem Zusammenhang lässt sich die Blazar-Sequenz als ein evolutionäres Szenario erklären: die Sequenz $FSRQ \rightarrow LBL/IBL \rightarrow HBL$ kommt aufgrund abnehmender Gasdichte der Hostgalaxie und damit einhergehender abnehmender Akkretionsrate zustande, dies wird durch weitere kosmologische Beobachtungen bestätigt. Eine abnehmende Materiedichte innerhalb des relativistischen Ausflusses wird von einem abnehmenden Magnetfeld begleitet, d.h. aber auch, dass Protonen weit vor den Elektronen nicht mehr im Strahlungsgebiet gehalten werden können. Die Blazar-Sequenz ist also ein Maß für die Hadronizität des Jets. Dies erklärt zudem die Dichotomie von FSRQs und BL Lac Objekten sowie die Zweiteilung in anderen Erscheinungsformen von AGN, z.B. FR-I und FR-II Radiogalaxien.\\ Während der Modellierung wird gezeigt, dass man Blazar-Spektren, speziell im hadronischen Fall, nicht mehr statisch betrachten kann, da es zu kumulierten Effekten aufgrund der langen Protonensynchrotronzeitskala kommt. Die niedrige Luminosität der Quellen und unterschiedlich lange Beobachtungszeiten verschiedener Experimente verlangen bei variablen Blazaren auch im leptonischen Fall eine zeitabhängige Betrachtung. Die Kurzzeitvariabilität scheint bei einzelnen Blazaren stets die selbe Ursache zu haben, unterscheidet sich aber bei der Betrachtung verschiedener Quellen. Zusätzlich wird für jeden Blazar, der in verschiedenen Flusszuständen beobachtet werden konnte, der Unterschied zwischen Lang- und Kurzzeitvariabilität, auch im Hinblick auf einen möglichen globalen Grundzustand hin, betrachtet.
Hochenergetische solare Teilchen werden bei ihrem Transport durch die Heliosphäre an turbulenten Magnetfeldern gestreut. Für das Verständnis dieses Streuprozesses ergeben sich aus heutiger Sicht zwei wesentliche Hindernisse: - Bei der Streuung hochenergetischer Teilchen an turbulenten Magnetfeldern handelt es sich um einen nichtlinearen Prozess, der durch analytische Theorien kaum zu beschreiben ist. - Der Streuprozess hängt stark von den tatsächlichen Magnetfeldern und somit auch von der Magnetfeldturbulenz ab. Unser bisheriges Verständnis der heliosphärischen Turbulenz ist leider aufgrund spärlicher experimenteller Daten deutlich eingeschränkt, was eine qualifizierte Umsetzung in analytischen und numerischen Ansätzen deutlich erschwert. Dies machte in der Vergangenheit künstliche Annahmen für die Modellerstellung notwendig. In dieser Arbeit wird der Teilchentransport mit Hilfe der Simulation von Testteilchen in einem turbulenten, magnetohydrodynamischen Plasma untersucht. Durch die Testteilchen werden auch die nichtlinearen Streuprozesse korrekt wiedergegeben, wodurch das erste hier genannte Hindernis überwunden wird. Dies wurde auch bereits in früheren numerischen Untersuchungen erfolgreich angewendet. Die Modellierung der Turbulenz für den Fall des Teilchentransports erfolgt in dieser Arbeit erstmalig auf Grundlage der magnetohydrodynamischen Gleichungen. Dabei handelt es sich um die mathematisch korrekte Wiedergabe der Magnetfeldturbulenz unterhalb der Ionen-Gyrofrequenz mit nur geringen numerischen Einschränkungen. Darüber hinaus erlaubt ein auf das physikalische Szenario anpassbarer Turbulenztreiber eine noch realistischere Simulation der Turbulenz. Durch diesen universell gültigen, numerischen Ansatz können für das zweite hier angegebene Hindernis jegliche künstlichen Annahmen vermieden werden. Die drei im Rahmen dieser Arbeit erstmals zusammengeführten Methoden (Testteilchen, magnetohydrodynamische Turbulenz, Turbulenztreiber) ermöglichen somit eine Untersuchung und Analyse von Transport- und Turbulenzphänomenen mit herausragender Qualität, die insbesondere für den Fall des Teilchentransports einen direkten Anschluss an experimentelle Ergebnisse ermöglichen. Wichtige Ergebnisse im Rahmen dieser Arbeit sind: - der Nachweis der Drei-Wellen-Wechselwirkung für schwache und einsetzende starke Turbulenz. - eine Analyse der Anisotropie der Turbulenz im Bezug auf das Hintergrundmagnetfeld in Abhängigkeit vom Treibmodell. Insbesondere die Anisotropie ist experimentell bislang kaum erfassbar. - eine Untersuchung der Auswirkung der Gyroresonanzen auf die Diffusionskoeffizienten hochenergetischer solarer Teilchen in allgemeiner Form. - die Simulation des Teilchentransports in der Heliosphäre auf Grundlage experimenteller Messdaten. Die genauere Analyse der Simulationsergebnisse ermöglicht insgesamt einen Zugang zum Verständnis des Transports, der durch experimentelle Untersuchungen nicht erfassbar ist. Bei der Simulation wurden lediglich die Magnetfeldstärke sowie die untersuchte Teilchenenergie vorgegeben. Aus der Analyse der Simulationsergebnisse ergibt sich dieselbe mittlere freie Weglänge, wie sie auch durch andere Verfahren direkt aus den Messergebnissen gewonnen werden konnte. Auch die vorwiegende Ausrichtung der hochenergetischen Teilchen parallel und antiparallel zum Hintergrundmagnetfeld in der Simulation entspricht experimentellen Untersuchungen. Es zeigt sich, dass diese allein aus den resonanten Streuprozessen der Teilchen mit den Magnetfeldern resultiert. Des Weiteren werden die Art der Diffusion, der Energieverlust der Teilchen während des Transportprozesses sowie die Gültigkeit der quasilinearen Theorie untersucht.
Im Rahmen dieser Arbeit wurde ein dreidimensionaler vollrelativistischer und parallelisierter Particle-in-Cell Code geschrieben, ausführlich getestet und angewandt. Der Code ACRONYM ist variabel einsetzbar und von der Genauigkeit und Stabilität her State-of-the-Art und somit konkurrenzfähig zu den sonstigen in der Astrophysik eingesetzten Codes anderer Gruppen. Die Energie bleibt bis auf einen Fehler von < 0.03% erhalten, die Divergenz des Magnetfeldes bleibt immer unter einem Wert von 10^{-12} und die Skalierung wurde mittlerweile bis zu einem Clustergröße von einigen 10000 CPUs getestet. In dieser Arbeit wurde dann, nach der Entwicklung des Codes, der Einfluss des fundamentalen Massenverhältnisses m_p/m_e auf die Teilchenbeschleunigung durch Plasmainstabilitäten untersucht. Dies ist relevant und wichtig, da in PiC-Simulationen in den allermeisten Fällen nicht mit dem realen Massenverhältnis gerechnet wird, da sonst viel zu viel Rechenleistung benötigt würde, um zu sehen, was mit den Protonen geschieht und was ihr Einfluss auf die leichten Teilchen wie Elektronen und Positronen ist. Zu diesem Zweck wurden Simulationen mit Massenverhältnissen zwischen m_p/m_e = 1.0 und 200.0 durchgeführt. Diese haben alle gemeinsam, dass periodische Randbedingungen verwendet wurden und das zur Verfügung stehende Simulationsgebiet mit jeweils zwei gegeneinander strömenden Plasmapopulationen vollständig gefüllt wurde, um jegliche Art von auftretenden Schocks auszuschließen. Die Rohdaten der einzelnen Simulationen wurden auf vielfältige Art und Weise analysiert, es wurden z.B. Schnitte durch die Teilchenverteilung erstellt, sowie ein- oder zweidimensionale Histogramme und Energieverläufe betrachtet. Dabei haben sich folgende Kernpunkte ergeben: Für Massenverhältnisse bis etwa m_p/m_e = 20 bildet sich die gesamte Zweistrom-Instabilität in nur einer Phase aus, das heißt, es bilden sich von ringförmigen Magnetfeldern umgebene Flussschläuche aus, die dann verschmelzen, bis nur noch zwei übrig sind und alle Teilchen werden über den gesamten Verlauf der Instabilität beschleunigt. Es ist damit zu folgern, dass die unterschiedlich schweren Teilchenspezies Protonen und Elektronen/Positronen durch die relativ nahe beieinander liegenden Massen noch so stark gekoppelt sind, dass sich nur eine Instabilität entwickeln kann. Bei großen Massenverhältnissen (m_p/m_e > 20) ist eine deutliche Trennung in zwei Phasen der Instabilität zu erkennen. Zuerst bilden sich wiederum Flussschläuche aus, diese verschmelzen miteinander (zu zweien oder mehr), bevor der erste Teil der Instabilität abflaut. Anschließend entstehen wieder ringförmige Magnetfelder und Flussschläuche, von denen einer meist deutlich stärker ist als all die anderen, das bedeutet, dass dieser von stärkeren Magnetfeldern umgeben ist und eine höhere Teilchendichte aufweist. Im Rahmen dieser zweigeteilten Instabilität werden die Elektronen und Positronen nur in der ersten Phase signifikant beschleunigt, die deutlich schwereren Protonen gewinnen über den gesamten Zeitraum Energie. Die höchstenergetischen Teilchen erreichen im Ruhesystem der jeweiligen Plasmapopulation Werte um gamma = 250. Man kann daraus für zukünftige Untersuchungen mit Hilfe von Particle-in-Cell Codes den Schluss ziehen, dass Rückschlüsse auf das tatsächliche Verhalten beim realen Massenverhältnis von m_p/m_e = 1836.2 nur aus den Simulationen mit m_p/m_e >> 20 gezogen werden können, da die starke Kopplung der leichten und schweren Teilchen bei kleineren Massenverhältnissen die Ergebnisse sehr stark beeinflusst. Es wurde anhand der gemessenen Zeitpunkte der Instabilitätsmaxima eine Extrapolation durchgeführt, die zeigt, dass die Instabilität beim realen Massenverhältnis etwa bei t = 1400 omega_{pe}^{-1} auftreten würde. Um dies wirklich zu simulieren müsste allerdings mehr als die 1000-fache Anzahl an CPU-Stunden aufgewandt werden. Des weiteren wurde eine Maxwell-Jüttner-Verteilung an die Teilchenverteilungen der einzelnen Simulationen auf dem Höhepunkt der Instabilität gefittet, um sowohl die neue Temperatur des Plasmas als auch die Beschleunigungseffizienz des Prozesses zu berechnen. Die Temperatur erhöht sich demnach durch die Instabilität von etwa 10^8K auf 10^{10} bis 10^{11}K, der Anteil suprathermischer Teilchen beträgt 2 bis 4%.