No description available.

Some of the brightest and most energetic events in the Universe are associated with the death of stars. These stellar deaths power transient electromagnetic emission which are routinely observed on Earth. This dissertation presents our research on various such transients. Its topics includes, supernova remnants, kilonovae, gamma-ray bursts (GRBs): The "long'' type produced from core-collapse supernovae and the "short'' type associated with neutron star merger events. It also focuses on the disruption of stars by the tidal forces of supermassive black holes i.e., tidal disruption events (TDEs). We model the emission from these transients and compare them to observations in order to draw a number of conclusions and make predictions for future detections. For example, we find that the non-thermal emission from supernovae and kilonovae associated with GRBs can produce long term emission which may be detected as a re-brightening in the overall emission. The sharp cut off observed in some TDE flares can be caused by a pre-existing accretion disk present around a supermassive black hole, which is expected in active galactic nuclei. Our work successfully predicted the nature of the very first electromagnetic detection from a neutron star merger, and was able to reproduce the emission that had been observed for more than one hundred days after the merger. This dissertation also provides frameworks on how the observable features of these transients can be leveraged to probe the properties of the progenitor system and their environment. <br>

X-ray binaries are binary stellar systems containing a compact object and a normal companion star which are gravitationally bound and rotate about a common center of mass. The compact object accretes matter from the companion star. The accreted matter may have a high angular momentum and hence follow a Keplarian orbit about the compact object. It slowly spirals inward as its angular momentum is redistributed via viscous forces and forms an accreting disk before being finally accreted onto the compact object. The compact object that is accreting matter may either be a neutron star or a black hole. X-ray binaries can be broadly classified into two classes depending on the mass of the companion star. Low Mass X-ray Binaries (LMXBs) have companion star masses and accrete mass via Roche lobe overflow of the companion star. High Mass X-ray Binaries (HMXBs) have companion star masses and in these systems the compact object accretes matter from the high velocity stellar winds of the companion star. For the work and results that are presented in the thesis we have studied the orbital evolution, apsidal motion and long term flux variations in High mass X-ray binaries which have a neutron star compact object with very high magnetic field of the order of B ~ 1012 G. Due to the high magnetic field, the accretion disk is disrupted at the Alfven radius where the magnetic field pressure equals the ram pressure of the infalling matter. From that boundary, the flow of the infalling matter will be guided by the magnetic field lines. The infalling matter will follow these lines, finally falling onto the magnetic poles with velocity nearly equal to the free fall velocity and form an accretion column over the magnetic poles. A hot spot is formed at both the magnetic poles and high energy photons are emitted from these regions. Inverse Compton scattering of these photons by high energy electrons in the accretion column can produce hard X-rays. If the optical depth of the accretion column is low, the radiation comes along the magnetic axis forming a pencil beam whereas if the optical depth is high, radiation escapes tangential to the accretion column forming a fan beam. Since the neutron star is rotating about its rotation axis, the radiation beam directed along magnetic axis non-aligned with the rotation axis will sweep across the sky. Whenever this beam of rotating radiation is aligned with the line of sight, a pulse of X-ray radiation is detected. Hence these systems are also called X-ray Binary Pulsars (XBP). These pulses are emitted at equal intervals of time, where the time between the emission of two pulses is the spin period of the neutron star. But since the neutron star is in a binary orbit, the arrival time of pulses as recorded by an observer will be delayed or advanced due to the motion of the neutron star. When the neutron star is moving towards the observer, the pulses arrive faster and when the neutron star is moving away from the observer, the pulses are delayed. These delays or advances of the arrival time of pulses can be measured accurately which allows us to measure the orbital elements (ax sin i, Porb, e, ω, Tω ) of the neutron star orbit. The neutron star orbit may evolve with time due to mass loss from the system, mass transfer from the companion star onto the neutron star and due to tidal interaction between the neutron star and the companion star. Gravitational wave radiation may also cause orbital evolution. However, in HMXBs this effect is likely to be much weaker compared to the effect of mass loss, mass exchange and tidal interaction. Rossi X-ray Timing Explorer (RXTE) is an X-ray astronomy satellite launched in 1995 by NASA. It has two pointed instruments, the Proportional Counter Array (PCA) and the High Energy X-ray Timing Experiment (HEXTE). PCA has a large effective area of 6500 sq cm and works in the energy range of 2-60 keV. It has a very good time resolution of 1 microsec. HEXTE observes in the energy range of 15-250 keV and has a time resolution of 8 microsec. RXTE also has an All Sky Monitor (ASM) which scans 80% of the sky in 90 minutes. We have used RXTE-PCA data for timing and spectral studies and ASM data for the long term flux variation studies of Cen X-3. The thesis presents details of our work, the analysis of the data, results of the analysis and our conclusions from these results. The first chapter of the thesis gives an overview of X-ray binaries, their orbital evolution and the instrument details of RXTE. In the second chapter we have presented our work of timing analysis of three persistent sources, namely Cen X-3, SMC X-1 and 4U 1538–52. For the SMC X-1 system, we have for the first time measured the eccentricity and the angle of periastron (ω). We found that the accuracy of pulse timing analysis is limited by the dependence of pulse profile on orbital phase. The new measurement of the orbit ephemeris of Cen X-3 when combined with the previous measurements of orbit ephemeris obtained by observations from other X-ray missions, gave an improved measurement of the rate of orbital decay P˙orb/Porb ~ -1.8 x 10−6yr−1 . A long observation of SMC X-1 made by RXTE in 2000 during the high state of SMC X-1 allowed us to measure the very small orbit eccentricity e ~ 0.00021 in this system. SMC X-1 was again observed for a long time by RXTE during 2003 during its low state. The SMC X-1 pulse fraction depends on the flux state of the source such that the pulse fraction decreases with decrease in the source flux. Thus the 2003 observations of SMC X-1 have higher error in measurement of pulse arrival times compared to the 2000 observations and could not be used to measure the eccentricity of the orbit. But combining the orbital ephemeris of SMC X-1 measured using the 2000 and the 2003 observation with the epoch history allowed us to improve the measurement of rate of orbit decay by an order of magnitude compared to previous observations P˙orb/Porb ~ - 3.4 x 10−6yr−1 . We observed 4U 1538–52 with RXTE under the guest observer program to measure the orbital evolution of this system. From observations of this system with BeppoSAX , a circular orbit similar to the SMC X-1 system was inferred. 4U1538–52 was observed with RXTE again in 1997 and analysis of this observation showed it to have eccentric orbit with a marginal evidence for an orbital decay. Our analysis carried out using the 2003 RXTE observation data confirmed that the orbit is eccentric with e ~ 0.18. But the new orbital ephemeris measured clearly shows that the orbit is not evolving with time as reported earlier. We have derived an upper limit on the rate of change of orbital period of this system to be P˙orb/Porb = 2.5 x 10−6yr−1 . 4U 1538–52 is similar to SMC X-1 in many respects, both have similar orbital period of Porb(SMC X -1) = 3.89 days and Porb(4U1538 - 52) = 3.72 days and companion star mass. But tidal interactions between the neutron star and the companion star have almost circularised the orbit of SMC X-1 where as the orbit of 4U 1538–52 is quite eccentric. Therefore we conclude that 4U 1538–52 is a young system and hence the orbit has not circularised by tidal interaction. The neutron star orbit also precesses due to tidal interaction and rotation of the companion star, which causes the longitude of periastron ω to change with time. The rate of change of ω can be measured by comparing the orbital elements of the neutron star orbit measured at different epochs of time. This rate of change of ω is directly related to the mass distribution of the companion star and hence the apsidal motion constant that are predicted by the theoretical models for stellar structure. Therefore measuring ˙ω will be a direct test for the stellar structure models. But ω can be measured only when the orbit is eccentric and for this purpose the Be-star/X-ray binary pulsars are the most suitable objects. The Be-star/X-ray binary pulsars are transient systems and have wide eccentric orbits of Porb > 10 days. The Be-stars are fast rotating stars with rotational velocity near to the break-up velocity. They eject matter along their equator in a circumstellar disk. When the neutron star intercepts this circumstellar disk during its periastron passage, the rate of mass accretion increases and the system becomes bright in X-rays. These short outburst are called the type-I X-ray bursts. The Be-star also has episodes of high mass ejection when the neutron star may accrete a larger amount of matter and can be seen over several binary orbits. These long duration outbursts are called type-II X-ray bursts. In the third chapter of the the thesis we have reported the analysis and results of three Be-/X-ray binary pulsars we have studied, namely 4U 0115+62, V0332+52 and 2S 1417-624 which were observed by RXTE during their respective type-II bursts. The X-ray pulse profiles of the Be-/X-ray systems evolve as a function of the source flux. Generally a simple single peaked pulse profile is seen during the onset of the outburst, which evolves into a more complex multiple peaked pulse profile as the source flux increases. The pulse profile again returns to the simple single peaked profile as the outburst fades off and the source flux decreases to persistent X-ray flux levels. Also due to varying mass accretion rate, the spin period evolves during the outburst. Both these factors together reduce the accuracy of measuring the arrival time of pulses. Hence we have used the instantaneous spin period measurements to deduce the orbital parameters of these system. The apparent spin period (Pspin) of the neutron star is modified by the radial velocity of the neutron star due to Doppler effect. The radial velocity of neutron star is dependant on the neutron star orbit and hence measurement of the spin period of the neutron star at different orbital phases allows us to determine the orbital elements. 4U 0115+63 was observed with the RXTE during two of its recent type II outbursts in 1999 and 2004. We measured the orbital parameters during both these outbursts independently. We combined the previous measurements of ω with our two measurements and measure the rate of apsidal motion of the system to be ˙ω ~ 0o .06 yr−1. V0332+52 was seen in outburst during 2004. During its previous outburst of 1983 only nine spin period measurements had been obtained and the orbital parameters measured from them were erroneous. We have measured the orbital parameters of this system accurately and determined the correct projected semi-major axis ax sin i and orbital period. The new orbit parameters can now be used to compare with future orbital element measurements to estimate any apsidal motion and/or orbital evolution in this system. We also used the 1999 outburst of 2S 1417–624 to accurately measure the orbital parameters of this system. We have also investigated the long term flux variations in the X-ray light curves of X-ray Binaries. Our studies on the flux variations observed in Cen X-3 are described in the fourth chapter of the thesis. Long term light curves of X-ray binaries show variations due to many reasons. Periodic variations of few milliseconds to a few hours in the light curve are seen due to spin of the neutron star. Light curves show variations due to motion of the neutron star in its orbit at timescales of few minutes to several days. Many sources also show quasi periodic variations in their X-ray light curves at timescales smaller than the neutron star orbital period which are believed to arise due to some material inhomogeneity orbiting the neutron star. These variations are called quasi periodic oscillations (QPOs). QPOs in X-ray binaries are observed between a frequency range of few millihertz to a few kilohertz. Long term X-ray light curves of many sources also reveal flux variations at time scales greater than the respective orbital period of the source. These variations are called superorbital variations. Systems like Her X-1, LMC X-4, 2S0114+650, SS 433, XTE J1716–389, 4U 1820–303 and Cyg X-1 show periodic superorbital variations whereas other systems like SMC X-1, GRS 1747-312, Cyg X2, LMC X-3 and the Rapid Burster show quasi periodic superorbital flux variations. These superorbital flux variations are understood as arising either due to a changing mass accretion rate which could be aperiodic in nature or as due to obscuration of the central X-ray source by a tilted, warped and precessing accretion disk. Many theoretical models have been proposed to explain the disk precession. The long term flux variations in the X-ray light curves of bright persistent X-ray binaries like Her X-1, SMC X-1 and LMC X-4 have been understood to be due to a periodic (in case of Her X-1 and LMC X-4) or a quasi periodic (for SMC X-1) precession of a warped accretion disk. We analysed the light curves of Cen X-3 obtained with the RXTE-ASM. The Cen X-3 light curves show aperiodic X-ray flux variations in all the energy bands of 1.5-3, 3-5 and 5-12 keV. The high and low states last for a few to upto a hundred days. The source also shows two spectral modes during the observations carried out with the ASM. The source was in a hard state during December 2000 to April 2004. At first look the aperiodic variations seen in Cen X-3 light curves seem to be arising due to a changing mass accretion rate. To investigate the cause of these aperiodic flux variations of Cen X-3 we studied the orbital modulation and the pulsed fraction as a function of source flux state. In the high state, the eclipse ingress and egress are found to be sharp whereas in the intermediate state, the transitions are more gradual. In the low state, instead of eclipse ingress and egress, the light curve shows a smooth intensity variation with orbital phase. The orbital modulation of the X-ray light curve in the low state shows that the X-ray emission observed in this state is from an extended object. The intensity dependent orbital modulations indicate that the different intensity states of Cen X-3 are primarily due to varying degree of obscuration. Measurements of the pulsed fraction in different intensity states are consistent with the X-ray emission of Cen X-3 having two components, one highly varying component with a constant pulsed fraction and a relatively stable component that is unpulsed and in the low state, the unpulsed component becomes dominant. The observed X-ray emission in the low state is likely to be due to scattering of X-rays from the stellar wind of the companion star. Though we can not ascertain the origin and nature of the obscuring material that causes the aperiodic long term intensity variation, we point out that a precessing accretion disk driven by radiative forces is a distinct possibility. We also studied the QPOs in Cen X-3 that are seen at 40 mHz. The QPOs are explained by the Beat Frequency Model (BFM) as arising due to the beat between the Keplarian frequency of the inner accretion disk and the spin of the neutron star. Thus when the mass accretion rate is high the inner disk radius decreases, increasing the Keplarian frequency and hence the observed QPO frequency and vise versa when the mass accretion rate decreases. Thus if the flux variations of Cen X-3 were due to a changing mass accretion rate then the observed QPO frequency should have a positive correlation with the observed X-ray flux of the source. But we find in our study that the QPO frequency does not have any correlation with the observed X-ray flux and the QPO frequencies does not follow the Frequency-Flux relation as expected in the Beat frequency model. Thus the QPO behaviour is in agreement that the observed X-ray flux does not indicate the true X-ray intensity state and hence the mass accretion rate in Cen X-3. Therefore, we conclude that X-ray variations of Cen X-3 are not due to changing mass accretion rate but due to varying obscuration of the central X-ray source, possibly by an accretion disk which precesses aperiodically. The conclusions from our studies presented in chapter 2, 3 and 4 of the thesis are summarised in the final chapter. The improved measurements of the rate of change of orbital periods from our work can now help us to detect any small departures from a constant period derivative in the persistent HMXB systems. The improved measurements of the orbital elements of Be-/X-ray binaries can now be used to study orbital evolution and apsidal motion in these system. New outbursts of the transient systems observed by future satellites providing good timing accuracy and large effective area, like LAXPC (Large Area X-ray Proportional Counter) of the ASTROSAT mission will facilitate such studies. The long term X-ray light curves study as done for Cen X-3 can be extended to other X-ray binary systems observed by All Sky Monitor. The method of source flux state dependent studies developed to study the Cen X-3 system can be easily extended to other systems that show long term superorbital flux variations. These kind of studies can be done by future proposed X-ray missions like ASTROSAT which will have a Sky Monitor similar to ASM dedicated to monitor X-ray sources. More sensitive measurements of long term X-ray light curves with the MAXI mission will allow similar studies of a large number of X-ray binaries and we will be able to see if aperiodically precessing accretion disk is present in many X-ray binaries.

<p dir="ltr">Neutron star mergers generate powerful gravitational waves and various types of electromagnetic signals, including gamma-ray bursts (GRB), kilonovae, and their afterglows. Observing and modeling these signals help us understand the physical processes of the merger events. Radiation from mergers can also serve as probes to study nuclear physics and cosmology. In this report, I focus on two types of signals: the GRB afterglow and the kilonova. GRB afterglows are non-thermal radiation produced by the interaction of relativistic jets and circumburst material, where the jets are launched perpendicular to the merger plane. Kilonovae are the thermal radiation emitted from the hot materials ejected during the merger. Besides the modeling of these objects, I also investigate their application in multimessenger astrophysics, especially the constraint on the expansion rate of the Universe. </p><p dir="ltr">First, I developed a GRB afterglow model to account for the off-axis observation of a structured jet. Using a jet structure derived from a three-dimensional general relativistic magnetohydrodynamic simulation, we performed a joint analysis of the multimessenger data of the neutron star merger event GW170817, including the gravitational wave data and GRB afterglow data in the radio band. We have tightly constrained the observing angle of GW170817 and broken the degeneracy between the inclination angle and luminosity distance measured in gravitational waves. With a better constrained distance, we improved the standard siren measurement of the Hubble constant to $H_0 = 69.5\pm 4\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$. The error bar has been reduced by a factor of 2. This work demonstrates that the modeling of off-axis GRB afterglow can significantly improve the standard siren method, provided that we have a reliable jet structure.</p><p dir="ltr">Second, I upgrade the GRB afterglow model in the first work, extending it to the late time where lateral spreading of the GRB jet becomes important. In this model, the ultra-relativistic blastwave is approximated by an infinitely thin two-dimensional surface. With this approximation, the hydrodynamic equations can be analytically integrated over the radius. Further assuming axial symmetry, the three-dimensional hydrodynamic simulation can be reduced to one dimension, which significantly increases the computational efficiency. We have compared our method to full numerical simulations and existing GRB afterglow modeling tools. The comparison shows good agreement and verifies our approach. Compared to these tools, our model has better flexibility and is applicable in a broader context. This method has been developed into a numerical code, \texttt{jetsimpy}, which we have provided to the community. It will serve as a powerful tool in the era of multimessenger astrophysics.</p><p dir="ltr">Finally, I investigate the possibility of long-lived massive neutron stars as neutron star merger remnants. A long-lived massive neutron star can inject a significant amount of energy into the merger ejecta, boosting the luminosity of kilonova by several orders of magnitude. However, this type of event has not yet been observed in optical sky surveys. We developed a boosted kilonova model with a detailed calculation of the photoionization process to better describe the efficiency of energy injection from spin down power to the ejecta. Our study found that boosted kilonovae, if commonly occurring, they should have already been observed given the accumulated time in sky surveys. As a result, the absence of detection implies that long-lived massive neutron stars as neutron star merger remnants are likely to be rare in the Universe.</p>

Στην παρούσα διατριβή υπολογίζονται σχετικιστικά πολυτροπικά μοντέλα περιστρεφομένων αστέρων νετρονίων, καθώς και μοντέλα που περιγράφονται από ρεαλιστικές καταστατικές εξισώσεις. Σκοπός αυτής της μελέτης είναι να υπολογιστούν σημαντικές φυσικές ποσότητες ενός αστέρα νετρονίων, στην περίπτωση της υδροστατικής ισορροπίας, της ομοιόμορφης αλλά και της διαφορικής περιστροφής, καθώς και στην περίπτωση που ο αστέρας έχει μαγνητικό πεδίο με πολοειδή και τοροειδή συνιστώσα. Μία σύντομη περιγραφή της αριθμητικής διαπραγμάτευσης έχει ως εξής. Καταρχάς, επιλύεται το σύστημα διαφορικών εξισώσεων Oppenheimer-Volkov (OV). Το σύστημα αυτό περιγράφει την υδροστατική ισορροπία μη περιστρεφομένων πολυτροπικών μοντέλων. Στη συνέχεια, θεωρείται η ομοιόμορφη περιστροφή ως διαταραχή, σύμφωνα με την «μέθοδο διαταραχής Hartle» και υπολογίζονται διορθώσεις στην μάζα και την ακτίνα, διορθώσεις που οφείλονται σε σφαιρικές και τετραπολικές παραμορφώσεις. Ακολούθως, εφαρμόζεται μία διαταρακτική προσέγγιση με όρους τρίτης τάξης στην γωνιακή ταχύτητα, Ω. Η στροφορμή, J, η ροπή αδράνειας, I, η περιστροφική κινητική ενέργεια, T, και η βαρυτική δυναμική ενέργεια, W, είναι ποσότητες που υφίστανται σημαντικές διορθώσεις από την προσέγγιση τρίτης τάξης. Η διαφορική περιστροφή ϑεωρείται ότι (i) υπακούει σε έναν συγκεκριμένο νόμο, ή (ii) επάγεται από το συνδυασμό ομοιόμορφης περιστροφής και ακτινικών ταλαντώσεων του αστέρα· ο στόχος είναι να υπολογισθεί η μεταβολή σημαντικών φυσικών ποσοτήτων που οφείλεται στη διαφορική περιστροφή. Στο δεύτερο μέρος, μελετάται η επίδραση του μαγνητικού πεδίου, το οποίο αποτελείται από πολοειδή και τοροειδή συνιστώσα, με τη «μέθοδο διαταραχής κατά Ioka-Sasaki» (IS). Στην παρούσα διαπραγμάτευση, το πρόβλημα περιγράφεται από μία «γενικευμένη διαφορική εξίσωση Grad-Shafranov» (GS),η επίλυση της οποίας δίνει τη συνάρτηση ροής (flux function), ψ. Μέσω αυτής της συνάρτησης υπολογίζονται οι συνιστώσες του μαγνητικού πεδίου και η γεωμετρική παραμόρφωση που υφίσταται ο αστέρας λόγω του μαγνητικού πεδίου. Η αντιμετώπιση του προβλήματος γίνεται και σε αυτήν την περίπτωση με τη ϑεωρία διαταραχών. ΄Εχοντας υπολογίσει μοντέλα περιστρεφομένων αστέρων νετρονίων και διάφορα μοντέλα με μαγνητικό πεδίο, μπορούμε να συνθέσουμε τα αποτελέσματά μας και να προσδιορίσουμε μοντέλα αστέρων νετρονίων μηδενικής φαινόμενης παραμόρφωσης (equalizers), δηλαδή αστέρων νετρονίων που η περιστροφή και το μαγνητικό πεδίο προκαλούν ίσες και αντίθετες γεωμετρικές παραμορφώσεις στο σχήμα του αστέρα. / We compute relativistic polytropic models as well as models obeying realistic equations of state, of rotating neutron stars. The purpose of this study is to calculate significant physical quantities of a neutron star, in the case of hydrostatic equilibrium, rigid and differential rotation, as well as in the case of a magnetic neutron star with both poloidal and toroidal components. A short description of the numerical treatment has as follows. First, we solve the Oppenheimer-Volkov system of differential equations. This system refers to hydrostatic equilibrium of non rotating polytropic models. Then, solid rotation is added as a perturbation, according to "Hartle’s perturbation method" and corrections to mass and radius are calculated, as also corrections due to spherical and quadrupole deformations. In addition a third order perturbation in angular velocity, Ω, is implemented. Angular momentum, J, moment of inertia, I, rotational kinetical energy, T, and gravitational potential energy, W, are quantites that are significally corrected by the third order approximation. Differential rotation is assumed that (i) obeys a specific law, or (ii) follows as a result of the solid rotation and radial oscillations combination; our purpose is the calculation of the main physical quantities that are altered by differential rotation. In the second part the effect of magnetic field is studied, which consists of a poloidal and a toroidal component. The "Ioka-Sasaki perturbation method" (IS) is implemented. This problem is described by the quantification of the flux function ψ, which comes as a solution of the "Grad-Shafranov" (GS) differential equation. Then the components of the magnetic field and the quadrupole deformation of the star are calculated. This method is also a perturbative method similar to "Hartle’s perturbation method". Having calculated models of rotating neutron stars, as also various models of magnetic fields, we can compose our results and determine models of neutron stars with zero deformation, the equalizers, these are neutron stars that are rotating and also have a magnetic field in a way that they, rotation and magnetic field, produce equal but opposite geometrical deformations in the shape of the star.

In this thesis, we construct a new hydrodynamical algorithm able of handling general compressible reactive flow problems, based on a finite-volume method inspired by the original MUSCL scheme of van Leer (1979). The algorithm is of second-order in the smooth part of the flow and avoids dimensional splitting. It uses MPI to achieve parallelism, and includes an astrophysical equation of state and a nuclear reaction network. It proves to be robust to tests cases. In particular it reproduces quite well the reactive and non-reactive results obtained with two different numerical methods (Fryxell & al. 1989, Busegnies & al. 2007). Moreover the time-dependent results are in agreement with the corresponding steady state solution. This gives us confidence in applying it to an astrophysical situation which has never been studied, the propagation of a detonation in conditions relevant to superbursts. The algorithm is described in (Noel & al. 2007).<p><p>In a firt step we obtain the detonation profiles in pure carbon and in a mixture of carbon and iron. In both cases we underline the large difference between the total reaction length and the length on which some species burn. This difference leads to enormous numerical difficulties because all the length scales cannot be resolved at the same time in a single simulation. We show that the carbon detonation might be studied in a partial resolution approach like the one of Gamezo & al. (1999).<p><p>In a second step we construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegrations of heavy elements, like ruthenium, which are thought to occur during superbursts in mixed H/He accreting systems. Using this new nuclear network we simulate detonations in mixture of carbon and ruthenium. An interesting feature is that, in this case, all the reaction lengths can be resolved in the same simulation. This makes the C/Ru detonations easier to study in future multi-dimensional simulations than the pure carbon ones (Noel & al. 2007b).<p><p>Finally we perform some numerical experiments which show that our algorithm is able to deal with initially inhomogeneous medium, and that the multi-dimensional simulations are attainable even if they are quite computational time consuming.<p><p>- B. Van Leer, J. Comp. Phys. 21, 101, 1979<p>- Fryxell, B.A. Muller, E. and Arnett, W.D. Technical report MPA 449, 1989<p>- Busegnies, Y. Francois, J. and Paulus, G. Shock Waves, 11, 2007<p>- Gamezo, V.N. Wheeler, J.C. Khokhlov, A.M. and Oran, E.S. ApJ, 512, 827, 1999<p>- Noël, C. Busegnies, Y. Papalexandris, M.V. & al. A&A, 470, 653, 2007<p>- Noël, C. Goriely, S. Busegnies, Y. & Papalexandris, M.V. submitted to A&A, 2007b<p><p>/<p><p>Un algorithme parallèle basé sur une méthode aux volumes finis inspirée du schéma MUSCL de Van Leer (1979) a été construit. Il a été développé sur base de la méthode de Lappas & al. (1999) qui permet de résoudre simultanément toutes les dimensions spatiales. Cette méthode se base sur la construction de surfaces appropriées dans l'espace-temps, le long desquelles les équations de bilan se découplent en équations plus simples à intégrer. Cet algorithme est actuellement le seul à éviter le "splitting" des dimensions spatiales. Dans les modèles conventionnels (PPM, FCT, etc.), l'intégration spatiale des équations est réalisée de manière unidimensionnelle pour chaque direction. <p>Un réseau de réactions nucléaires ainsi qu'une équation d'état astrophysique ont été inclus dans l'algorithme et celui-ci a ensuite été soumis à une grande variété de cas tests réactifs et non réactifs. Il a été comparé à d'autres codes généralement utilisés en astrophysique (Fryxell & al. 1989, Fryxell & al. 2000, Busegnies & al. 2007) et il reproduit correctement leurs résultats. L'algorithme est décrit dans Noël & al. (2007).<p><p>Sur base de cet algorithme, les premières simulations de détonation dans des conditions thermodynamiques représentatives des Superbursts ont été réalisées. Différentes compositions du milieu ont été envisagées (carbone pur, mélange de carbone et de fer, mélange de carbone et de cendres du processus rp). Dans la plupart des systèmes où des Superbursts ont été observés, la matière accrétée est un mélange d'hydrogène et d'hélium. Dans ce cas, des phases de combustion précédant le Superburst produisent des nucléides plus lourd que le fer (Schatz & al. 2003). Ces nucléides peuvent être photodésintégrés durant le Superburst. Pour prendre en compte ces réactions endothermiques de photodésintégration, nous avons construit un nouveau réseau réduit de réactions nucléaires qui a été incorporé dans l'algorithme hydrodynamique (Noël & al. 2007b). Ce réseau réduit reproduit globalement l'énergétique d'un réseau complet et a permis de faire la première simulation numérique de détonation dans des conditions caractéristiques de systèmes accréteurs d'un mélange hydrogène-hélium. <p>Finallement quelques simulations multidimensionelles préliminaires ont éte réalisées.<p><p>- Busegnies, Y. Francois, J. and Paulus, G. Shock Waves, 11, 2007<p>- Fryxell, B.A. Muller, E. and Arnett, W.D. Technical report MPA 449, 1989<p>- Fryxell, B.A. Olson, K. Ricker, P. & al. ApJS, 131, 273, 2000<p>- Lappas, T. Leonard, A. and Dimotakis, P.E. SIAM J. Sci. Comput. 20, 1064, 1999<p>- Noël, C. Busegnies, Y. Papalexandris, M.V. & al. A&A, 470, 653, 2007<p>- Noël, C. Goriely, S. Busegnies, Y. & Papalexandris, M.V. submitted to A&A, 2007b<p>- Röpke, F. K. PhD thesis, Technischen Universitat Munchen, 2003<p>- Schatz, H. Bildsten, L. Cumming, A. and Ouellette, M. Nuclear Physics A, 718, 247, 2003<p>- Van Leer, B. Comp. Phys. 21, 101, 1979<p>- Weinberg, N.N. and Bildsten, L. ArXiv e-prints, 0706.3062, 2007 / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Cherenkov-Teleskope sind in der Lage, das schwache Cherenkovlicht aus Teilchenschauern zu detektieren, die von kosmischen Teilchen mit Energien von ca. 100 GeV bis 100 TeV in der Erdatmosphäre initiiert werden. Das Ziel ist die Detektion von Cherenkovlicht aus Schauern, die von Gammastrahlen erzeugt wurden, der größte Teil der Schauer stammt jedoch von geladenen Teilchen. Im Jahr 2012 wurde das H.E.S.S.-Observatorium in Namibia, bis dahin bestehend aus vier Teleskopen mit 100 m²-Spiegeln, um ein fünftes Teleskop mit einer Spiegelfläche von ca. 600 m² ergänzt. Aufgrund der großen Spiegelfläche besitzt dieses Teleskop die niedrigste Energieschwelle aller Teleskope dieser Art. In dieser Dissertation wird ein schneller Algorithmus namens MonoReco präsentiert, der grundlegende Eigenschaften der Gammastrahlen wie ihre Energien und Richtungen rekonstruieren kann. Dieser Algorithmus kann weiterhin unterscheiden, ob Schauer von Gammastrahlen oder von geladenen Teilchen der kosmischen Strahlung initiiert wurden. Diese Aufgaben werden mit mithilfe von künstlichen neuronalen Netzwerken erfüllt, welche ausschließlich die Momente der Intensitätsverteilungen in der Kamera des neuen Teleskops analysieren. Eine Energieschwelle von 59 GeV und Richtungsauflösungen von 0.1°-0.3° werden erreicht. Das Energiebias liegt bei wenigen Prozent, die Energieauflösung bei 20-30%. Unter anderem mit dem MonoReco-Algorithmus wurden Daten, die in der Zeit um das Periastron des Binärsystems PSR B1259-63/LS 2883 im Jahre 2014 genommen wurden, analysiert. Es handelt sich hierbei um einen Neutronenstern, der sich in einem 3,4-Jahres-Orbit um einen massereichen Stern mit einer den Stern umgebenden Scheibe aus Gas und Plasmen befindet. Zum ersten Mal konnte H.E.S.S. das Gammastrahlenspektrum dieses Systems bei Energien unterhalb von 200 GeV messen. Weiterhin wurde bei erstmaligen Beobachtungen zur Zeit des Periastrons ein lokales Flussminimum gemessen. Sowohl vor dem ersten als auch nach dem zweiten Transit des Neutronensterns durch die Scheibe wurden hohe Flüsse gemessen. Im zweiten Fall wurden Beobachtungen erstmals zeitgleich mit dem Fermi-LAT-Experiment durchgeführt, das wiederholt sehr hohe Flüsse in diesem Teil des Orbits messen konnte. Ein Vergleich der gemessenen Flüsse mit Vorhersagen eines leptonischen Modells zeigt gute Übereinstimmungen. / Cherenkov telescopes can detect the faint Cherenkov light emitted by air showers that were initiated by cosmic particles with energies between approximately 100 GeV and 100 TeV in the Earth's atmosphere. Aiming for the detection of Cherenkov light emitted by gamma ray-initiated air showers, the vast majority of all detected showers are initiated by charged cosmic rays. In 2012 the H.E.S.S. observatory, until then comprising four telescopes with 100 m² mirrors each, was extended by adding a much larger fifth telescope with a very large mirror area of 600 m². Due to the large mirror area, this telescope has the lowest energy threshold of all telescopes of this kind. In this dissertation, a fast algorithm called MonoReco is presented that can reconstruct fundamental properties of the primary gamma rays like their direction or their energy. Furthermore, this algorithm can distinguish between air showers initiated either by gamma rays or by charged cosmic rays. Those tasks are accomplished with the help of artificial neural networks, which analyse moments of the intensity distributions in the camera of the new telescope exclusively. The energy threshold is 59 GeV and angular resolutions of 0.1°-0.3° are achieved. The energy reconstruction bias is at the level of a few percent, the energy resolution is at the level of 20-30%. Data taken around the 2014 periastron passage of the gamma-ray binary PSR B1259-63/LS 2883 were analysed with, among others, the MonoReco algorithm. This binary system comprises a neutron star in a 3.4 year orbit around a massive star with a circumstellar disk consisting of gas and plasma. For the first time the gamma-ray spectrum of this system could be measured by H.E.S.S. down to below 200 GeV. Furthermore, a local flux minimum could be measured during unprecedented measurements at the time of periastron. High fluxes were measured both before the first and after the second transit of the neutron star through the disk. In the second case measurements could be performed for the first time contemporaneously with the Fermi-LAT experiment, which has repeatedly detected very high fluxes at this part of the orbit. A good agreement between measured fluxes and predictions of a leptonic model is found.

Βασικός άξονας της παρούσας διατριβής είναι οι σχετικιστικοί υπολογισμοί πολυτροπικών μοντέλων περιστρεϕόμενων αστέρων νετρονίων. Επειδή δεν υπάρχει ακριβής αναλυτική λύση των εξισώσεων του Einstein για το ϐαρυτικό πεδίο ενός περιστρεϕόμενου αστέρα νετρονίων, επιχειρούμε την αϱιθμητική επίλυση στο μιγαδικό επίπεδο όλων των διαϕορικών εξισώσεων, που εμπεριέχονται στην διαταρακτική μέθοδο του Hartle. Δίνουμε έμϕαση στον υπολογισμό φυσικών ποσοτήτων, που περιγράϕουν την γεωμετρία ταχέως περιστρεϕόμενων μοντέλων. Συγκρίνοντας τα αριθμητικά αποτελέσματα που ϐρίσκουμε με ορισμένες πολύπλοκες επαναληπτικές μεθόδους, ελέγχουμε την αξιόλογη ϐελτίωση των αποτελεσμάτων μας, έναντι εκείνων που δίνονται από το κλασσικό διαταρακτικό σχήμα του Hartle. Η παρούσα διατριβή χωρίζεται σε τέσσερα μέρη, που αναπτύσσονται στα κεϕάλαια 1, 2, 3 και 4. Στο πρώτο κεϕάλαιο, ϑα εστιάσουμε την προσοχή μας στο σύστημα διαφορικών εξισώσεων Oppenheimer − Volkov, που εξάγονται από τις εξισώσεις πεδίου του Einstein. Σε συνδυασμό με μια καταστατική εξίσωση περιγράφουμε σχετικιστικά πολυτροπικά μοντέλα μη περιστρεϕόμενων αστέρων νετρονίων σε υδροστατική ισορροπία. Ακολούθως, περιγράϕουμε ένα καθαϱά σχετικιστικό φαινόμενο, τον συρμό των αδρανειακών συστημάτων λόγω της περιστροϕής του αστέρα. Στην συνέχεια, χρησιμοποιούμε την μέθοδο διαταραχής του Hartle, σύμϕωνα με την οποία δεχόμαστε ότι ο στατικός αστέρας είναι το αδιατάρακτο σύστημα, πάνω στο οποίο εϕαρμόζουμε μικρές διαταραχές (ϑεωρώντας την ομοιόμορϕη περιστροϕή ως διαταραχή) και έτσι υπολογίζουμε τις διορθώσεις στην μάζα και την ακτίνα, λόγω των σϕαιρικών και τετραπολικών παραμορϕώσεων. Τέλος, εϕαρμόζουμε μία διαταρακτική προσέγγιση με όρους τρίτης τάξης στην γωνιακή ταχύτητα. Στο δεύτερο κεϕάλαιο, ϑα κάνουμε μια εκτενή περιγραϕή της στρατηγικής του μιγαδικού επιπέδου (Complex-Plane Strategy, εν συντομία CPS). Σύμϕωνα με αυτή την μέθοδο, η αριθμητική ολοκλήρωση των διαϕορικών εξισώσεων γίνεται στο μιγαδικό επίπεδο και όλες οι εμπλεκόμενες συναρτήσεις του προβλήματός μας είναι μιγαδικές, μιγαδικής μεταβλητής. Συνεπώς, για την αποϕυγή διαϕόρων ιδιομορϕιών ή και απροσδιόριστων μορϕών, που προκύπτουν από τις οριακές συνθήκες του προβλήματος, κυρίως στο κέντρο και στην επιϕάνεια του αστέρα, μας δίνεται η δυνατότητα να επιλέξουμε ένα κατάλληλο μιγαδικό μονοπάτι για την εκτέλεση πάνω σ΄ αυτό της αριθμητικής ολοκλήρωσης των διαϕορικών εξισώσεων. Επιπλέον, οι αριθμητικές ολοκληϱώσεις όλων των διαϕορικών εξισώσεων του προβλήματος συνεχίζονται πολύ πέραν της επιϕάνειας του αδιατάρακτου μοντέλου, με αποτέλεσμα η ακτίνα υπολογίζεται εύκολα ως η ϱίζα του πραγματικού μέρους της συνάρτησης της πυκνότητας (χωρίς να είμαστε αναγκασμένοι να εκτελέσουμε οποιεσδήποτε αριθμητικές προεκβολές, που είναι γνωστό ότι επιϕέρουν σημαντικά σϕάλματα). Στο τρίτο κεϕάλαιο, υπολογίζουμε σημαντικές φυσικές ποσότητες που αφορούν τον αστέρα νετρονίων, ολοκληρώνοντας αριθμητικά ένα σύστημα διαφορικών εξισώσεων πρώτης τάξης. Ιδιαίτερα, υπολογίζουμε το σύνορο της περιστρεϕόμενης αστρικής δομής με δύο τρόπους. Ο πρώτος είναι με ϐάση την κλασική διαπραγμάτευση της διαταρακτικής μεθόδου του Hartle και ο δεύτερος με τον αλγόριθμο λεπτής ϱύθμισης που αναπτύσσουμε με την ϐοήθεια του οποίου παίρνουμε αξιόλογα αριθμητικά αποτελέσματα. Στην συνέχεια περιγράϕουμε το λογισμικό πακέτο ATOMFT System. Ακολούθως, με την ϐοήθεια των λύσεων των διαϕορικών εξισώσεων τρίτης τάξης ως προς την γωνιακή ταχύτητα, υπολογίζουμε τις διορθώσεις στην στροϕορμή, την ϱοπή αδράνειας, την περιστροϕική κινητική ενέργεια και την ϐαρυτική δυναμική ενέργεια του αστέρα. Εϕαρμόζοντας τέλος μια κατάλληλη μέθοδο, υπολογίζουμε το όριο της μάζας διαϕυγής. Στο τέταρτο κεϕάλαιο, εκθέτουμε πίνακες αποτελεσμάτων και κάποιες σημαντικές γραϕικές παραστάσεις. Δίνουμε επίσης ορισμένες λεπτομέρειες της εϕαρμογής του προγράμματός μας. Επιπλέον, δίνουμε έμϕαση στο γνωστό «παράδοξο» που αϕορά την μέθοδο διαταραχών του Hartle,σύμϕωνα με την οποία αυτή η μέθοδος αν και αντιπροσωπεύει μια προσέγγιση αργής πεϱιστροϕής ενός αστέρα νετρονίων, δίνει αξιόλογα αποτελέσματα ακόμη και όταν εϕαρμόζεται σε ταχέως περιστρεϕόμενα μοντέλα. Στην παρούσα έρευνα αϕαιρέσαμε τον κρίσιμο περιορισμό του τερματισμού των αριθμητικών ολοκληρώσεων λίγο πριν από την επιϕάνεια του μη περιστρεϕόμενου αστέρα, συνεχίζοντας την ολοκλήρωση αρκετά πέραν του συνόρου του. Αυτό σημαίνει ότι η CPS ¨γνωρίζει¨ την παραμόρϕωση που προκαλείται από την περιστροϕή για ένα αρκετά εκτεταμένο διάστημα που περιβάλλει την αρχικά σϕαιρική μορϕή του αστέρα. Συνεπώς, για τους υπολογισμούς που απαιτούνται για τον περιστρεϕόμενο αστέρα, η CPS δεν προεκβάλλει ποτέ, με αποτέλεσμα τα σϕάλματα των υπολογισμών είναι πολύ μικρά. Τέλος, λαμβάνοντας υπόψη κατάλληλα στους υπολογισμούς μας ένα ορισμένο αριθμό συνθηκών, συνδυάζοντας την κλασική διαπραγμάτευση του διαταρακτικού σχήματος του Hartle και τις σχέσεις που απορρέουν από την δομή της στρατηγικής του μιγαδικού επιπέδου, οδηγηθήκαμε τελικά στην επινόηση του αλγόριθμου λεπτής ϱύθμισης, αποτέλεσμα του οποίου είναι η σημαντική ϐελτίωση της ακρίβειας των αριθμητικών αποτελεσμάτων που αϕορούν την γεωμετρία του συνόρου του αστέρα νετρονίων. ΄Αμεση συνέπεια όλων αυτών είναι ο υπολογισμός με ικανοποιητική ακρίβεια του ορίου της μάζας διαϕυγής, εϕαρμόζοντας μια κατάλληλη μέθοδο. / In the present dissertation we solve numerically in the complex plane all the differential equations involved in Hartle’s perturbation method for computing general-relativistic polytropic models of rotating neutron stars. We give emphasis on computing quantities describing the geometry of models in rapid rotation. Compared to numerical results obtained by certain sophisticated iterative methods, we verify appreciable improvement of our results vs to those given by the classical Hartle’s perturbative scheme. The description of the present investigation is constituted by four parts and has as follows. In the first chapter, we start to describe the nonrotating neutron star model. Then, according to "Hartle’s perturbation method", the solid rotation is added as a perturbation. So, the equations of structure for uniformly rotating stars are given up to second order in the angular velocity and the distortions to mass and radius are calculated as corrections owing to spherical and quadrupole deformations. Subsequently, the equations are given up to third order in the angular velocity. In the second chapter, we describe extensively the numerical method called Complex-Plane Strategy (abbreviated CPS). According to this method, we solve numerically in the complex plane all the differential equations involved in Hartle’s perturbation method. Any function of our problem is interpreted as a complex-valued function of a complex variable. CPS offers an alternative for avoiding any singularities and/or indeterminate forms, especially near the center and the surface of the nonrotating star, by performing numerical integration along a proper complex path. Moreover, the numerical integrations of all the differential equations governing the problem are continued well beyond the surface of the nonrotating star, thus, the radius is readily calculated as root of the density function (without been forced to perform any numerical extrapolations). In the third chapter, we solve numerically in the complex plane the system of first-order differential equations resulting from Hartle’s perturbation method. We give emphasis on computing the boundary of the rotating configuration by the so-called fine tuning algorithm which gives appreciably improved results. Then, we describe the software systems that we use in our investigation, with emphasis on the ATOMFT System. Finally, we compute the third order corrections in the uniform angular velocity for the angular momentum, moment of inertia, rotational kinetical energy and gravitational potential energy. Furthermore, we describe a method for computing the mass-shedding limit. In the fourth chapter, we present several numerical results and some significant graphical representations. We also give certain details of our program implementation. Concluding, we emphasize on the well-known "paradox" concerning Hartle’s perturbation method, according to which this method, although representing a slow-rotation approximation, gives remarkably accurate results even when applied to rapidly rotating models. In the present work, we have removed the certain critical limitations of terminating integrations below the radius of the star. Instead, the numerical integration of our problem continues well beyond the boundary of the star. This means that CPS knows the distortion to be caused by rotation over a sufficiently extended space surrounding the initially spherical configuration. So, to the computation of a particular rotating configuration, CPS never extrapolates beyond the end of the function tables computed by such extended numerical integrations. It is exactly the avoidance of any extrapolation which keeps the error in the computations appreciably small. Finally, we have properly taken into account certain conditions matching Hartle’s perturbative scheme and the relations arising in the framework of the Complex-Plane Strategy. This treatment has led to the fine tuning algorithm which, in turn, has improved appreciably the accuracy of our numerical results related to the geometry of the star’s boundary. Consequently, the mass-shedding limit can be calculated using a proper procedure which gives remarkably accurate results.

The aim of this thesis is to create collections of examples for the subject Astronomy taught for students of pedagogical faculties, studying this discipline as a part of physics courses. Due to very different mathematical knowledge of students I have chosen typical and not much difficult examples oriented to several branches of astronomy. Each part of examples begins with a self-contained theoretical introduction. The difficulty rises gradually from trivial to more complicated examples. The examples are mainly focused on motions in radial gravitational fields.

This thesis concerns with the study of X-ray binaries which are gravitationally bound systems consisting of a compact object (either a neutron star or a black hole) and usually a non degenerate companion star, both rotating around the common centre of mass. The compact star shines brightly in the X-ray regime. Emission from these systems are powered by accretion which is the most radioactively efficient mechanism known in the universe by the release of gravitational potential energy when matter from the companion star falls on the compact object. Accretion onto high magnetic field neutron stars are special as the magnetic field plays a crucial role in governing the dynamics of gas flow and the flow of the matter close to the compact object. The radiation emitted from these systems are anisotropic and for a distant observer, the intensity is modulated at the spin period of the neutron star, hence these objects are called accretion powered pulsars. The angular pattern of the emitted radiation is also highly anisotropic and depends on the mass accreted and hence the luminosity. The beaming pattern commonly known as the pulse profiles exhibit a wide variety in the pulse shape and pulse fraction and vary with energy as well as intensity. They also exhibit cyclotron absorption features in their energy spectrum which are a direct probe to the magnetic field geometry of these systems. This thesis is dedicated to the study of the magnetic field and emission geometry of accretion powered pulsars through the pulse phase resolved studies of the cyclotron absorption features which are a direct probe of the magnetized plasma. In order to study these features in detail broadband continuum modeling of the energy spectrum is done, taking care of all other factors which may smear the pulse phase dependence. Another prerequisite for detailed continuum modeling is accounting for the low absorption dips in the pulse profiles of many these sources. The dips are presumably formed by phase locked accretion stream causing partial covering absorption when the stream is along our line of sight towards the emission region. Studying the pulse phase dependence of this partial covering absorber also provides us with important clues on the local environment of the neutron star and the structure of the accretion stream. All of these studies are performed with data from the broadband and most sensitive instruments onboard the Japanese satellite Suzuki. Lastly we provide estimates of the polarization expected to be detected from these sources by a Thomson scattering polarimeter being developed to observe the polarization of X-rays in the energy range of 5--30 keV. Along with the X-ray pulsars, we also make an estimate of the likelihood of detection of X-ray polarization from black hole X-ray binaries in different spectral states. This is a particularly interesting topic as it will play a crucial role in providing additional handles on the magnetic field geometry in accretion powered pulsars as well as constrain the fundamental parameters of a black hole like its spin.