Cette thèse porte sur la modélisation théorique de l'équation d’état (EE) décrivant la matière nucléaire présente dans le coeur des étoiles à neutrons (EN), sous l'hypothèse qu'aucune transition de phase ne s'y produise. Nous utilisons un méta-modèle permettant i) d’incorporer directement les connaissances en physique nucléaire sous la forme de paramètres empiriques tels que la densité de saturation nucléaire, l’incompressibilité, l’énergie de symétrie; ii) de reproduire la plupart des modèles nucléoniques existants; et iii) d’explorer les régions inconnues à haute densité de façon la plus large possible. Pour chaque EE, nous déterminons un ensemble de solutions pour la masse et le rayon des EN, et nous effectuons une première sélection des EE compatibles avec la stabilité et la causalité de la matière nucléaire, ainsi que la masse maximale connues des EN. Nous confrontons ensuite ces EE aux observations d’émission thermique dans la gamme des rayons-X pour 7 EN soigneusement choisies. Pour la première fois, la modélisation théorique des EE est directement introduite dans l’analyse des données. Nous utilisons les dernières mesures effectuées par GAIA II pour fixer la distance des EN. Les paramètres du modèle d’émission thermique et de l’EE sont déterminés selon une méthode Bayésienne basée sur un algorithme Monte-Carlo par Chaîne de Markov. Nous déterminons ainsi la température de surface, la masse et le rayon des EN, ainsi que sur la valeur de certains paramètres empiriques tels que la dépendance en densité de l'énergie de symétrie (Lsym), la contribution isovectorielle au module d’incompressibilité (Ksym) ou encore le paramètre de distorsion isoscalaire (Qsat) / This thesis concerns theoretical modeling of the equation of state (EoS) describing nuclear matter in the core of neutron stars (NS), under the hypothesis that no phase transition occurs. We use a meta-model which is able to i) directly incorporate nuclear physics knowledge on the form of empirical parameters such as the nuclear saturation density, the incompressibility or the symmetry energy; ii) reproduce most of the existing models; iii) explore new behaviors at high densities in a very flexible way. For each EoS, we determine a set of solutions for the masses and radii of NS, and we make a first selection of the EoS that are compatible with the stability and causality constraints, as well as the maximum observed mass of NS. Then we confront these EoS to observational data coming from thermal emission in the soft X-ray domain, for 7 NS carefully chosen. For the first time, the theoretical modeling of the EoS is directly implemented in the data analysis. We use the recent measurments of GAIA II to constrain the distance to the NS. The parameters of the modeling of thermal emission as well as the empirical parameters entering in the EoS are determined by Bayesian methods using a Monte-Carlo by Morkov Chain algorithm. Therefore, we determine the surface effective tempreature, the masses and radii of NS, as well as some empirical nuclear parameters such as the density dependance of the symmetry energy (Lsym), the isovector incompressibility (Ksym), or the isoscalar squewness (Qsat)

Cette thèse est consacrée à l'étude des microquasars, des sources binaires de haute énergie appartenant à notre Galaxie, analogues morphologiquement aux quasars, mais sur des échelles de longueur et de masse considérablement inférieures. La première partie concerne l'un des meilleurs représentants des microquasars: GRS 1915+105, découvert par le télescope WATCH/GRANAT, et qui présente des mouvements en apparence superluminiques. Nous avons établi la nature galactique de cette source, située à une distance de 12.5 +/- 1.5 kpc, ainsi que son environnement, à partir des observations millimétriques. Nos observations de cette source dans le domaine infrarouge ont montré qu'elle subissait des variations importantes, de 1 à 2 magnitudes, sur des intervalles de temps variables. Il apparaît que GRS 1915+105 est un système binaire lumineux et de grande masse, dont l'objet compact, certainement un trou noir, est entouré d'un disque d'accrétion optiquement épais. Grâce à un ensemble de campagnes d'observations simultanées à plusieurs longueurs d'onde, nous avons mis en évidence le lien entre les phénomènes d'accrétion et d'éjection de matière se produisant autour de l'objet compact. Ces observations, en établissant le modèle du microquasar, montrent que les analogies morphologiques entre quasars et microquasars sont sous-tendues par des analogies dynamiques et physiques. Enfin, grâce à des observations millimétriques, nous avons analysé l'effet d'éventuelles interactions entre les éjections énergétiques de cette source et le milieu interstellaire environnant. Dans la deuxième partie de cette thèse, j'expose une étude traitant de plusieurs sources de haute énergie, pour lesquelles nous avons identifié et/ou étudié les contreparties infrarouge et radio. Il ressort de cette étude que ce sont les objets compacts, trous noirs ou étoiles à neutrons, dans des systèmes binaires de faible masse, qui subissent les plus grandes variations de luminosité infrarouge.

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

Nous avons étudié ici une nouvelle forme d'excitation, dite collective, dans la croûte interne des étoiles à neutrons. Dans cette partie de l'étoile, en se structurant, la matière prend des formes aux géométries originales appelées " phase pasta " et dont l'état est superfluide dans des conditions de température caractéristiques des étoiles à neutrons. L'objectif était de décrire des excitations dont les longueurs d'onde varient de la taille de la cellule de Wigner-Seitz à des très grandes valeurs. Ces dernières sont essentielles pour l'étude à basse énergie et représentent l'aspect macroscopique de ce milieu alors que les premières correspondent aux aspects microscopiques de la matière. Pour cela nous avons opté pour une approche fondée sur l'hydrodynamique des superfluides à plusieurs constituants. Ainsi, nous avons pris en compte les propriétés essentielles de la matière nucléaire telles que les vitesses du son. Le modèle que nous avons développé consiste essentiellement en l'étude de la propagation de perturbations à travers ces structures, par rapport à un équilibre hydrostatique. Nous nous sommes penchés sur le cas des structures en forme de plaques uniquement. Pour cela nous avons été conduits à définir des conditions aux bords, afin de modéliser les effets des interfaces. La périodicité du milieu a été prise en compte en utilisant le théorème de Floquet-Bloch. La résolution des équations nous a permis d'obtenir les spectres d'excitations de ces plaques. Nous en avons dégagé les propriétés physiques et nous avons utilisé ces spectres pour le calcul de chaleur spécifique. Ceci nous a permis d'étudier la contribution de ces excitations collectives en les comparant aux autres formes d'excitations possibles. Ainsi nous avons constaté que cette nouvelle contribution était loin d'être négligeable. Elle aura sa place à part entière dans les futurs calculs d'évolution thermique des étoiles à neutrons. Plus de détails sur http://lucdigallo.free.fr/spip.php?article14&lang=fr

L'expérience H.E.S.S. (High Energy Stereoscopic System), située en Namibie, est un réseau de télescopes basé sur la technique d'imagerie Tcherenkov atmosphérique en stéréoscopie et dédié à l'astronomie gamma entre 100 GeV et quelques dizaines de TeV. H.E.S.S. dispose d'une situation géographique idéale pour observer notre galaxie, la Voie Lactée, et a ainsi pu découvrir des dizaines de nouvelles sources gamma grâce à sa stratégie de relevé systématique du plan galactique. Après environ dix ans d'observations fructueuses, les limites de l'instrument doivent désormais être repoussées pour pouvoir détecter de nouvelles sources plus faibles. Le développement d'algorithmes avancés permettant d'améliorer la sensibilité et la résolution angulaire de l'instrument est au coeur de cette thèse. Un outil de recherche automatique de nouvelles sources ainsi qu'une amélioration de la soustraction du bruit de fond sont présentés dans le cadre de l'étude des sources très faibles, nécessitant une analyse particulièrement rigoureuse. Une méthode combinant différentes techniques de reconstruction des événements est développée et permet d'améliorer la résolution angulaire sans diminuer la statistique, cruciale pour l'étude des sources faibles. Ces techniques avancées sont appliquées à l'analyse d'une région complexe du plan galactique, autour du vestige de supernova G22.7-0.2, qui a conduit à la détection d'une nouvelle source, HESS J1832-092. Les contreparties multi-longueurs d'onde à cette source sont présentées et différents scénarios sont envisagés pour expliquer l'origine de l'émission gamma de cet objet astrophysique.

Cette thèse est consacrée à l'étude des oscillations périodiques détectées lors des sursauts X des étoiles à neutrons, dans des binaires X de faible masse. Ces oscillations offrent un moyen de sonder l'intérieur de ces objets, en mesurant notamment leur masse et leur rayon, pour ainsi contraindre l'équation d'état de la matière dense. J'ai développé des méthodes de détection et d'analyse de ces signaux, de leurs propriétés temporelles et leur dépendance en énergie. J'ai analysé les oscillations détectées dans tous les sursauts X de type 1 (ainsi qu'un super-sursaut) de 3 étoiles à neutrons observées avec l'instrument Rossi X-ray Timing Explorer/Proportional Counter Array . Sur les courbes de lumière des sursauts, j'ai sélectionné les segments donnant la meilleure signification statistique, pour construire un catalogue de profils moyens d'oscillations. La forme des profils varie grandement d'un sursaut à un autre, pour une même source. Un grand nombre de paramètres peuvent affecter les oscillations. J'ai élaboré un modèle de tache chaude à la surface de l'étoile en rotation rapide pour caractériser l'émission du sursaut X, dans un espace-temps relativiste. En utilisant les chaînes de Markov Monte Carlo pour explorer efficacement un espace des paramètres conséquent, les ajustements sur un échantillon de sursauts ont démontré l'applicabilité du modèle. Par contre, les contraintes obtenues sur la masse et le rayon de l'étoile sont limitées par la qualité des données de l'instrument utilisé. Enfin, des simulations révèlent que des mesures précises sur les paramètres sont possibles en augmentant la surface collectrice des détecteurs, comme le proposent les observatoires X du futur.