Observations of the X-ray Pulsars XTE J1946+274, 4U 0115+634, and GX 304−1 Der Naturwissenschaftlichen Fakultät der Friedrich-Alexander-Universität Erlangen-Nürnberg zur Erlangung des Doktorgrades Dr. rer. nat. vorgelegt von Sebastian Müller aus Fürth Als Dissertation genehmigt von der Naturwissenschaftlichen Fakultät der Friedrich-Alexander-Universität Erlangen-Nürnberg Tag der mündlichen Prüfung: 28. Oktober 2013 Vorsitzender der Prüfungskommission: Prof. Dr. Johannes Barth Erstberichterstatter: Prof. Dr. Jörn Wilms Zweitberichterstatter: Dr. José Miguel Torrejón 3 4 Zusammenfassung Akkretierende Röntgenpulsare sind Doppelsternsysteme, die aus einem Neutronen- stern und einem optischen Begleitstern bestehen. Mit Dichten ähnlich der von Atom- kernen und magnetischen Feldstärken in der Größenordung von etwa 1012 G gehören diese Neutronensterne zu den extremsten Objekten im Universum. Die optischen Be- gleitsterne der in dieser Arbeit betrachteten Systeme sind sogenannte „Be-Sterne“, welche relativ hohe Temperaturen und große Massen besitzen. Ein Hauptmerkmal eines Be-Sterns ist die in hohem Maße veränderliche, äquatoriale Scheibe, von der Materie auf den Neutronenstern übertragen werden kann. Die Wechselwirkungen zwischen dem akkretierten Material und dem Magnetfeld führen zur Bildung von Akkretionssäulen an den magnetischen Polen des Neutronensterns. Die dabei entste- henden Akkretionshügel am Fuß der Akkretionssäulen emittieren thermische Strah- lung, die wiederum mit dem einfallenden Material wechselwirken. Diese Wechsel- wirkungen, sowie Bremsstrahlungsprozesse innerhalb der Akkretionssäulen, führen zur Emission von Röntgenstrahlung, was zum Auftreten der hellsten Quellen dieser Strahlung am Himmel führen kann. Die Magnetfeldstärken in der näheren Umgebung von Neutronensternen über- steigen die von Menschen technisch realisierbaren Felstärken um viele Größenord- nungen. Aus diesem Grund bieten Röntgenpulsare die einzigartige Möglichkeit, die Physik von Materie, die solch starken Magnetfeldern ausgesetzt ist, zu studieren. So sagt etwa die Theorie eine Quantisierung der Bewegung der Elektronen senkrecht zu den Magnetfeldlinien voraus. Die entsprechenden Übergänge liegen, bei den vor- liegenden Magnetfeldstärken energetisch im Bereich der Röntgenstrahlung, was bei einigen Röntgenpulsaren in Absorption erscheinende Spektrallinien zur Folge hat. Diese Linien werden auch Zyklotronlinien genannt. Die beobachtete Abhängigkeit der Zentralenergie dieser Linien von der Röntgenleuchtkraft erlaubt es, Rückschlüs- se über die Emissionsregionen und somit auch über die physikalischen Prozesse in- nerhalb der Akkretionssäulen zu ziehen. Aus theoretischen Modellierungen dieser Säulen geht hervor, dass es unterschiedliche Leuchtkraftbereiche mit unterschied- lichen Vorhersagen für die Variation der Energie der Zyklotronlinien gibt. So sollte für relativ niedrige Leuchtkräfte eine positive Korrelation zwischen dem Röntgen- fluss und der Linienenergie auftreten, während für hohe Leuchtkräfte eine negative Korrelation dieser Größen erwartet wird. Eines der Hauptprobleme im Zusammenhang mit der Analyse von Röntgenspek- tren akkretierender Pulsare ist die Modellierung des Breitbandkontinuums. Vor al- lem die kaum verstandene Geometrie, aber auch die komplexen physikalischen Pro- zesse innerhalb der Akkretionssäulen machen die Modellierung sehr anspruchsvoll und zeitaufwendig. Folglich sind Modellspektren von Röntgenpulsaren, die auf phy- sikalischen Grundlagen basieren, noch in der Entwicklungsphase. Stattdessen wird bei den meisten Analysen auf empirische Breitbandmodelle zurückgegriffen. In dieser Arbeit beschreibe ich Beobachtungen der drei Röntgenpulsare XTE J1946+274, 4U 0115+634 und GX 304–1 mit verschiedenen Röntgenteleskopen und erörtere deren wissenschaftliche Analyse und Interpretation. Für das System XTE J1946+274 wurde 2001 die Entdeckung und Identifizierung einer Zyklotronlinie veröffentlicht (Heindl et al., 2001). Meine Analyse von RXTE-, INTEGRAL- und Swift-Daten zweier neuerer Ausbrüche dieser Quelle zeigt, dass diese Linie bei der damals bestimmten Energie nicht mehr vorhanden ist. Ferner finde ich schwache Hinweise darauf, dass die Linie sich stattdessen bei signifikant i niedrigeren Energien befindet. Typische Variabilitäten, hervorgerufen durch Ände- rungen des Röntgenflusses, können eine derart starke Änderung der Linienenergie nicht erklären. Dies würde bedeuten, dass XTE J1946+274 für verschiedene Ausbrü- che unterschiedlichen Ausbruchsmechanismen unterliegt. Mit fünf bestätigten Zyklotronlinien spielt 4U 0115+634, mit der höchsten Anzahl gesicherter Linien, eine herausragende Rolle unter allen Röntgenpulsaren. Praktisch alle bisherigen Analysen dieser Quelle resultierten in einer negativen Korrelation zwischen dem Röntgenfluss und der Energie der fundamentalen Zyklotronlinie. In dieser Arbeit kann ich dieses Ergebnis unter Verwendung empirischer spektraler Standardmodelle reproduzieren. Jedoch finde ich heraus, dass dieses Verhalten auf ein unphysikalisches Zusammenspiel zwischen den Kontinuums- und den Zyklotron- parametern zurückzuführen ist. Bei Anwendung eines alternativen Kontinuummo- dels auf den gleichen Datensatz verschwindet diese negative Korrelation und die Fundamentallinie zeigt keine erkennbare Variabilität der Energie für alle beobach- teten Röntgenflüsse. Dieses Ergebnis hat starke Auswirkungen auf die physikalische Interpretation dieses Paradebeispiels eines Röntgenpulsars. Dies ist das erste Mal, dass ein derart starker Einfluss der Wahl des Kontinuumsmodells auf die Parame- ter der Zyklotronlinie nachgewiesen wird. Ich bin in der Lage, dieses Ergebnis unter Verwendung von Suzaku-Daten eines späteren Ausbruchs zu bestätigen. Bei der Analyse des Systems GX 304–1 zeigt sich ein ungewöhnlich starker An- stieg der Wasserstoff-Säulendichte, die nach etwa drei Tagen wieder auf das ur- sprüngliche Niveau zurückgegangen ist. Mögliche Szenarien, die diese Beobachtung erklären könnten, werden diskutiert. Eine der vielversprechendsten Erklärungen ist die Bedeckung des Röntgenstrahlung emittierenden Neutronensterns durch die Be- Scheibe des optischen Begleitsterns. Das Hauptresultat dieser Arbeit ist die Abhängigkeit freier Fitparameter, vor al- lem der Zyklotronlinien, und deren Verhalten von der Wahl des Kontinuummodells. Im Rahmen dieser Arbeit wird die Stärke und Reichweite solcher Einflüsse verdeut- licht. Aus diesem Grund sind künftig Modelle, die auf physikalischen Prozessen fu- ßen, unabdingbar, um sichere Erkenntnisse über Röntgenpulsare zu erlangen. ii Abstract An accreting X-ray binary pulsar consists of a neutron star and an optical companion star. With densities of the order of atomic nuclei and magnetic field strengths of ∼1012 G, these neutron stars are among the most extreme objects in the universe. The optical companions of the systems treated in this work are so called “Be-type stars”. These relatively hot and massive stars host a highly variable circumstellar disk. Triggered by the periastron passage of the neutron star, by the expansion of the Be-disk, and/or by other processes, mass originating from the Be-disk can be transferred onto the neutron star. Interactions between the strong magnetic field and the accreted matter lead to the formation of columns at the magnetic poles of the neutron star. Thermal emission from accretion mounds, which are located at the polar caps of the neutron star, interacts with the infalling matter. These interactions and bremsstrahlung processes within the accretion columns lead to the emission of X-rays resulting in the brightest sources of this radiation in the sky. Since the magnetic field strengths in the vicinity of neutron stars exceed man- made field strengths by many orders of magnitude, observations of X-rays from these binary systems allow for studying and testing of physical theories for such extreme conditions, e.g., the quantization of the electrons’ motion perpendicular to the mag- netic field line into discrete Landau levels. In some accreting X-ray pulsars, tran- sitions of electrons at the respective energies imprint absorption features onto the X-ray spectra, the so called “cyclotron resonance scattering features”, or short “cy- clotron lines”. The behavior of the centroid cyclotron line energy depending on the X-ray luminosity allows conclusions about the variation of the X-ray emitting region. These results can be used to draw conclusions about the physical processes within the accretion column. Theory predicts several X-ray luminosity regimes with a different behavior of the cyclotron line, i.e., for relatively low luminosities a positive correla- tion between the X-ray flux and the cyclotron line energy is expected, while for high luminosities a negative correlation should be observed. One of the central problems concerning the analysis of X-ray spectra of these outbursts is the modeling of the broadband X-ray continua. The hardly known ge- ometries and the complex physical processes within the accretion columns make it very complicated to physically model the X-ray spectra. As a consequence, model spectra based on a physical description are still under development. Instead, there are several different semi-physical empirical models, which are used in most cases to describe the broadband continuum. In this thesis I report on X-ray observations of outbursts of the three neutron star binary pulsars XTE J1946+274, 4U 0115+634, and GX 304–1 with different tele- scopes and discuss the respective analyses and the resulting conclusions. For XTE J1946+274 a cyclotron