Robert Brose

• Supernova remnants are believed to be the source of cosmic rays with energies up to 10^15 eV that are produced within our Galaxy. The acceleration mechanism associated with the collision-less shocks in supernova remnants - diffusive shock acceleration - predicts a spectral index of the accelerated non-thermal particles of s = 2. However, measurements of non-thermal emission in radio, X-rays and gamma-rays reveal significant deviations of the particles spectral index from the canonical value of s = 2. The youngest Galactic supernova remnant G1.9+0.3 is an interesting target for next-generation gamma-ray observatories. So far, the remnant is only detected in the radio and the X-ray bands, but its young age of ≈100 yrs and inferred shock speed of ≈ 14, 000 km/s could make it an efficient particle accelerator. I performed spherical symmetric 1D simulations with the RATPaC code, in which I simultaneously solved the transport equation for cosmic rays, the transport equation for magnetic turbulence, and the hydro-dynamical equations forSupernova remnants are believed to be the source of cosmic rays with energies up to 10^15 eV that are produced within our Galaxy. The acceleration mechanism associated with the collision-less shocks in supernova remnants - diffusive shock acceleration - predicts a spectral index of the accelerated non-thermal particles of s = 2. However, measurements of non-thermal emission in radio, X-rays and gamma-rays reveal significant deviations of the particles spectral index from the canonical value of s = 2. The youngest Galactic supernova remnant G1.9+0.3 is an interesting target for next-generation gamma-ray observatories. So far, the remnant is only detected in the radio and the X-ray bands, but its young age of ≈100 yrs and inferred shock speed of ≈ 14, 000 km/s could make it an efficient particle accelerator. I performed spherical symmetric 1D simulations with the RATPaC code, in which I simultaneously solved the transport equation for cosmic rays, the transport equation for magnetic turbulence, and the hydro-dynamical equations for the gas flow. Separately computed distributions of the particles accelerated at the forward and the reverse shock were then used to calculate the spectra of synchrotron, inverse Compton, and Pion-decay radiation from the source. The emission from G1.9+0.3 can be self-consistently explained within the test-particle limit. I find that the X-ray flux is dominated by emission from the forward shock while most of the radio emission originates near the reverse shock, which makes G1.9+0.3 the first remnant with non-thermal radiation detected from the reverse shock. The flux of very-high-energy gamma-ray emission from G1.9+0.3 is expected to be close to the sensitivity threshold of the Cherenkov Telescope Array. The limited time available to grow large-scale turbulence limits the maximum energy of particles to values below 100 TeV, hence G1.9+0.3 is not a PeVatron. Although there are many models for the acceleration of cosmic rays in Supernova remnants, the escape of cosmic rays from these sources is yet understudied. I use our time-dependent acceleration code RATPaC to study the acceleration of cosmic rays and their escape in post-adiabatic Supernova remnants and calculate the subsequent gamma-ray emission from inverse-Compton scattering and Pion decay. My simulations span 100,000 years, thus covering the free-expansion, the Sedov-Taylor, and the beginning of the post-adiabatic phase of the remnant’s evolution. At later stages of the evolution cosmic rays over a wide range of energy can reside outside of the remnant, creating spectra that are softer than predicted by standard diffusive shock acceleration and feature breaks in the 10 - 100 GeV-range. The total spectrum of cosmic rays released into the interstellar medium has a spectral index of s ≈ 2.4 above roughly 10 GeV which is close to that required by Galactic propagation models. I further find the gamma-ray luminosity to peak around an age of 4,000 years for inverse-Compton-dominated high-energy emission. Remnants expanding in low-density media emit generally more inverse-Compton radiation matching the fact that the brightest known supernova remnants - RCW86, Vela Jr, HESSJ1721-347 and RXJ1713.7-3946 - are all expanding in low density environments. The importance of feedback from the cosmic-rays on the hydrodynamical evolution of the remnants is debated as a possibility to obtain soft cosmic-ray spectra at low energies. I performed spherically symmetric 1-D simulations with a modified version of the RATPaC code, in which I simultaneously solve the transport equation for cosmic rays and the hydrodynamical equations, including the back-reaction of the cosmic-ray pressure on the flow profiles. Besides the known modification of the flow profiles and the consequently curved cosmic-ray spectra, steady-state models for non-linear diffusive shock acceleration overpredict the total compression ratio that can be reached with cosmic-ray feedback, as there is limited time for building these modifications. Further, I find modifications to the downstream flow structure that change the evolutionary behavior of the remnant and trigger a cosmic-ray-induced instability close to the contact discontinuity, if and when the cosmic-ray pressure becomes dominant there.…
• Es wird vermutet das Supernovaüberreste die Quelle der galaktischen kosmischen Strahlung mit Energien bis zu 10^15eV sein können. Der Beschleunigungsprozess der mit den kollisionsfreien Schocks in Supernovaüberresten in Verbindung gebracht wird - diffuse Schockwellenbeschleunigung - sagt nicht-thermische Teilchenspektren mit einem Spektralindex von s=2 voraus. Messungen nicht-thermischer Strahlung im Radio-, Röntgen- und Gammastrahlenbereich zeigen teils deutliche Abweichungen von dieser Vorhersage. Der jüngste galaktische Supernovaüberrest G1.9+0.3 ist ein interessantes Ziel für zukünftige Gammastrahlenteleskope. Bis jetzt wurde der Überrest nur im Radio- und Röntgenband entdeckt aber sein geringes Alter von ~100 Jahren und die gemessenen hohen Schockgeschwindigkeiten von ~14,000km/s sollten Teilchenbeschleunigung auch bis zu sehr hohen Energien ermöglichen. In dieser Arbeit wurden 1D-Simulationen der Teilchenbeschleunigung in G1.9+0.3 mit Hilfe der RATPaC-Programmbibliothek durchgeführt, wobei das System der gekoppeltenEs wird vermutet das Supernovaüberreste die Quelle der galaktischen kosmischen Strahlung mit Energien bis zu 10^15eV sein können. Der Beschleunigungsprozess der mit den kollisionsfreien Schocks in Supernovaüberresten in Verbindung gebracht wird - diffuse Schockwellenbeschleunigung - sagt nicht-thermische Teilchenspektren mit einem Spektralindex von s=2 voraus. Messungen nicht-thermischer Strahlung im Radio-, Röntgen- und Gammastrahlenbereich zeigen teils deutliche Abweichungen von dieser Vorhersage. Der jüngste galaktische Supernovaüberrest G1.9+0.3 ist ein interessantes Ziel für zukünftige Gammastrahlenteleskope. Bis jetzt wurde der Überrest nur im Radio- und Röntgenband entdeckt aber sein geringes Alter von ~100 Jahren und die gemessenen hohen Schockgeschwindigkeiten von ~14,000km/s sollten Teilchenbeschleunigung auch bis zu sehr hohen Energien ermöglichen. In dieser Arbeit wurden 1D-Simulationen der Teilchenbeschleunigung in G1.9+0.3 mit Hilfe der RATPaC-Programmbibliothek durchgeführt, wobei das System der gekoppelten Differentialgleichungen für den Teilchentransport, den Transport der magnetischen Turbulenz und der Standardgasgleichungen gelöst wurde. Die separat berechneten Verteilungen der Teilchen an Vorwärts- und Rückwärtsschock wurden benutzt um die Emission des Überrests im Radio-, Röntgen und Gammastrahlungsbereich zu bestimmen. Die Emissionen von G1.9+0.3 können selbst konsistent in der Testteilchennäherung bestimmt werden. Die Röntgenemission wird vom Vorwärtsschock dominiert, während die Radioemissionen hauptsächlich vom Rückwätsschock stammen. Dies macht G1.9+0.3 zum ersten Überrest mit detektierter nich-thermischer Strahlung aus dem Bereich des Rückwärtsschocks. Die erwartet Gammastrahlungsemission ist nahe dem Detektionslimit des zukünftigen Cherenkov Telescope Arrays. Die geringe Alter von G1.9+0.3 begrenzt die Maximalenergie, die im Beschleunigungsprozess erreicht werden kann auf Werte unterhalb von 100TeV. Demnach ist G1.9+0.3 kein PeVatron. Auch wenn es zahlreiche Modelle zur Teilchenbeschleunigung in Supernovaüberresten gibt, ist das Entkommen der Teilchen aus den Überresten zur Zeit wenig erforscht. Mit Hilfe von RATPaC haben wir die Evolution und Teilchenbeschleunigung in Supernovaüberresten über 100,000 Jahre simuliert. Dieser Zeitraum deckt einen Großteil der Lebensspanne eines Supernovaüberrests ab und endet mit der letzten Teil der postadiabatischen Phase der Entwicklung des Überrests. In den späten Phasen der Entwicklung des Überrests können Teilchen in einem großen Energiebereich aus dem Überrest entweichen. Dies erzeugt Emissionspektren, die weicher sind als durch Fermibeschleunigung vorhergesagt, und die spektrale Brüche im Bereich von 10-100GeV aufweisen. Das Produktionsspektren der Überreste hat einen spektralen Index von s~2.4 oberhalb von 10GeV, was ungefähr mit den Spektren übereinstimmt, die Quellen in galaktischen Propagationsmodellen aufweisen müssen. Weiterhin erreichen die Überreste ihre größte Helligkeit im Gammabereich nach etwa 4000 Jahren wenn diese durch inverse Comptonstreuung erzeugt werden. Dabei erreichen Überreste in Medien mit geringer Umgebungsdichte größere Helligkeiten, was sich mit den Beobachtungen der hellsten Supernovaüberreste - RCW86, Vela Jr., HESSJ1721-347 und RXJ1713.7-3946 - deckt, die alle in Bereichen sehr geringer Umgebungsdichte expandieren. In der Literatur wird die Möglichkeit diskutiert, dass die Rückkopplung der Beschleunigten Teilchen auf die Struktur der Überreste für die beobachteten weichen Strahlungsspektren verantwortlich ist. Im Rahmen dieser Arbeit wurde eine modifizierte Version von RATPaC entwickelt, die diesen Prozess abbilden kann. Neben den bekannten Rückkopplungen zeigt sich, dass bisherige Modelle unter der Annahme eines Gleichgewichtszustandes für die beschleunigten Teilchen die erreicht maximal Kompression des Plasmas durch den Schock und damit die Härte der Teilchenspektren überschätzen. In unseren zeit aufgelösten Berechnungen ist die maximale Kompression durch die limitierte verfügbare Zeit begrenzt. Zusätzlich zeigt sich das Auftreten einer Instabilität die durch die Rückkopplung der kosmischen Strahlung nahe der Kontaktdiskontinuität hervorgerufen wird.…