Lors du traitement d'une tumeur par radiothérapie utilisant des ions comme le 12C, le faisceau peut subir des réactions nucléaires avec les noyaux des tissus menant à la création de particules secondaires plus légères. Cela a pour conséquences de délocaliser une partie de la dose dans les tissus sains et de créer un champ d'irradiation mixte. Afin d'améliorer les connaissances sur la fragmentation du projectile, une expérience a été réalisée par notre collaboration en Mars 2015 au GANIL avec un faisceau de 12C à 50 MeV/n. Cette expérience a permis de mesurer les sections efficaces doublement différentielles en angle et en énergie pour tous les isotopes créés, et ce, sur différentes cibles d'intérêt médical (C, H, O et Ca). Les distributions angulaires et les taux de production ont également été obtenus pour chaque fragment et pour chaque cible. Les résultats ont ensuite été comparés à différents modèles de réactions nucléaire disponibles dans les outils de simulation Monte Carlo GEANT4 et PHITS. Des écarts importants ont été observés avec les résultats expérimentaux, mais aussi entre les différents modèles. Ces conclusions avaient déjà poussé notre collaboration à créer un nouveau modèle, appelé SLIIPIE, spécialisé dans les réactions d'ions légers aux énergies intermédiaires. Les paramètres de ce modèle ont ainsi été ajustés de façon à reproduire au mieux les données expérimentales obtenues. / During the treatment of cancerous diseases using ions such as 12C, the beam will undergo nuclear reactions with the tissues leading to the production of lighter secondary fragments. This tends to delocate a part of the dose into healthy tissues and create a mixed radiation field. In order to improve the knowledge on the projectile fragmentation, an experiment was performed by our collaboration in March 2015 with a 50 MeV/n 12C beam at GANIL. This led to the measurement of double differential cross sections in angle and energy of every isotope created by the fragmentation of carbon on targets of medical interest (C, H, O and Ca). Energy and angular distributions, as well as production rates have also been obtained for every fragment and every target. Experimental data were then compared with nuclear reaction models included in the Monte-Carlo (MC) simulation codes, GEANT4 and PHITS. Important discrepancies were observed with our data, as well as between models. These observations had already led our collaboration in creating a new MC model called SLIIPIE, specialized in intermediate energies reactions with light ions. The parameters of this model were then adjusted in order to best fit the experimental data.

L’étude présentée porte sur la caractérisation et la comparaison des propriétés quasi-statiques et en fatigue de composites à fibres de lin et de verre avec une matrice époxy. Une importante campagne d’essais de fatigue a été réalisée à partir des propriétés quasi-statiques de traction et cisaillement plan mesurées. Il apparaît que les caractéristiques en statique et en fatigue du verre/époxy sont supérieures à celles du lin/époxy à taux de fibres et stratification identiques. Cependant, en raison de la faible densité des renforts de lin, les écarts entre les propriétés spécifiques des deux matériaux sont faibles. La mesure de l’évolution des propriétés en fatigue a mis à jour des comportements phénoménologiques particuliers tels que la diminution de l’amortissement et l’augmentation de la rigidité des lin/époxy  ayant des fibres parallèles à la direction du chargement, au cours de la vie des éprouvettes. Le suivi de l’endommagement par la mesure des densités de fissures montre que celle-ci augmente avec le chargement et le nombre de cycles appliqués. Les simulations par éléments finis de la durée de vie des composites lin/époxy, intégrant la variabilité des paramètres du modèle par la méthode Monté-Carlo, donnent des courbes de Wöhler conservatives par rapport aux données expérimentales. Conformément aux mesures expérimentales, la variabilité des durées de vie calculées diminue lorsque le chargement baisse / The present study focuses on the characterization and comparison of the quasi-static and fatigue properties of flax and glass fibre reinforced composites with an epoxy matrix. An extensive experimental campaign of fatigue tests have been performed using the measured quasi-static tension and in plane-shear properties. It appears that the static and fatigue characteristics of the glass/epoxy composites are higher than that of flax/epoxy ones with similar fibre volume fraction and lay-ups. However, due to the low density of flax reinforcements, the differences between the specific properties of the two materials are low. The measurement of the evolution of the fatigue properties has highlighted particular phenomenological behaviours such as the decrease of the damping and the increase of the stiffness of flax/epoxy composites having fibres parallel to the loading direction, during specimens’ life. The monitoring of material damage by measuring the crack density showed its increase with the loading level and the number of cycles applied. Finite elements simulations of the life cycle of flax/epoxy composites, incorporating the variability of the model’s parameters by Monte Carlo method, gave conservative Wöhler curves compared with experimental data. In agreement with the experimental measurements, the variability of lifetimes calculated decreases when the load decreases

L'étude des propriétés de transport dans les matériaux magnétiques a débuté dans les années 1950. Le magnétisme intrinsèque à ce type de solide représente une source de diffusion supplémentaire par rapport aux matériaux non-magnétiques qui vient s'ajouter aux phonons et aux impuretés ou défauts cristallins du réseau. L'étude de la diffusion magnétique est plus complexe que les deux précédentes sources diffusives énoncées. Car d'une part l'influence de la diffusion magnétique induit des comportements différents de résistivité dus au matériau, selon que ce dernier est un métal ou un semi-conducteur. D'autre part, le type de magnétisme porté par la structure est également à l'origine de différents comportements de la résistivité magnétique. Dès les années 1950 de nombreux auteurs se sont concentrés sur l'étude des métaux et semi-conducteurs magnétiques. Il résulte des différents travaux que chacune des théories s'applique à un cas particulier et ne peut pas rendre compte de l'ensemble des comportements susceptibles d'être observés expérimentalement. A titre d'exemple, la résistivité des métaux magnétiques présente une évolution monotone fonction de la température (seule la dérivée par rapport à la température présente une singularité), alors que les semi-conducteurs magnétiques présentent un pic à la température critique (Tc), température séparant les phases d'ordre et de désordre magnétique. Ajoutons que les différentes interprétations concernant la diffusion magnétique autour de la température critique est encore un point controversé. Aussi, les théories disponibles à ce jour concernant l'étude du transport dans les matériaux magnétiques doivent être choisies en fonction du matériau étudié (type de structure, métal ou semi-conducteur, type de magnétisme). Une méthode générale, capable de rendre compte de l'ensemble des comportements de la résistivité quelle que soit la nature du solide (métal ou semi-conducteur) et du type de magnétisme fait aujourd'hui cruellement défaut. C'est à ce manque que nous avons souhaité répondre en premier lieu dans le cadre de cette thèse. En effet, le développement spectaculaire de la spintronique repose sur l'utilisation des semi-conducteurs magnétiques qui ne peut faire fi d'une compréhension rigoureuse des mécanismes de diffusion autour de la température critique. Notre démarche a consisté à utiliser la méthode numérique Monte Carlo Métropolis afin d'étudier ce problème sous un angle nouveau. Nous avons ainsi construit un algorithme qui s'applique à un Hamiltonien assez général du système. Grâce à cet Hamiltonien et avec un jeu de paramètres d'entrées appropriées, nous sommes en mesure par le biais de la méthode Monte Carlo Métropolis de reproduire de façon générique les tendances expérimentales présentes dans la littérature pour des semi-conducteurs ferromagnétiques et antiferromagnétiques, frustrés et non-frustrés. Egalement, la méthode permet une étude systématique de matériaux caractérisés expérimentalement avec prise en compte des unités afin de réaliser des comparaisons directes entre nos résultats de simulations et les mesures expérimentales, avec un bon accord de forme et de grandeur, par exemple dans le cas du semi-conducteur magnétique MnTe. Si cette thèse se limite à l'étude de matériaux semi-conducteurs, nos perspectives vont consister à étendre l'étude aux systèmes dilués DMS, très étudiés aujourd'hui, ainsi qu'aux problèmes de diffusion aux interfaces dans les dispositifs GMR, problématique qui soulève actuellement de nombreuses questions. A ce jour notre démarche est unique et semble capable de rendre compte des comportements de résistivité expérimentale en particulier pour le cas des semi-conducteurs, de fournir une nouvelle explication quant à la controverse liée aux mécanismes de diffusion électronique autour de la température critique. Ma contribution à ce travail a consisté à développer, en collaboration avec mon directeur de thèse le Pr. Hung The Diep, un algorithme pour l'étude du transport de spin. Grâce à cet algorithme, j'ai obtenu de nombreux résultats qui ont permis l'interprétation des différents résultats expérimentaux présentés dans ma thèse. Ma contribution principale à ce travail est : - Mise au point d'une méthode Monte Carlo avec des techniques permettant la réduction des fluctuations statistiques des échantillons (multi-step averaging) - Etude des différents types de matériaux ferromagnétiques et antiferromagnétiques non frustrés. Observation des comportements très marquants autour de Tc et à basse T: la dépendance de la forme du pic à Tc dépend de variables physiques bien identifiées (interactions : portée, nature ; temps de relaxation, réseau, ...) - Etude de l'effet de la frustration sur la résistivité. Effet d'une transition du premier ordre sur la résistivité. - Interprétation du comportement de la résistivité utilisant les données numériques notamment en terme de clusters et de temps de relaxation. - Etude quantitative de MnTe et comparaison avec l'expérience : accord parfait. Développons ces points en quelques mots. Le premier point sur lequel je voudrais insister, est l'interprétation nouvelle en ce qui concerne les mécanismes de diffusion à proximité des températures de transition ordre/désordre pour les réseaux non-frustrés. En effet, j'ai mis en évidence lors de ma thèse que la forme du pic de résistivité autour de Tc est une conséquence directe de deux effets liés au magnétisme du réseau. Autour de Tc le réseau percole et présente de larges clusters de spins parallèles (antiparallèles) qui représentent autant de zones de basses (hautes) énergies pour les électrons qui s'y propagent. En addition les spins du réseau subissent un ralentissement critique qui rend l'évolution du paysage magnétique extrêmement lent. Le ralentissement critique des fluctuations des spins du réseau fige le paysage magnétique et permet aux spins itinérants de se localiser dans les clusters parallèles (basses énergies) et d'y être piégés, donnant naissance à un pic de résistivité à Tc. En conséquence nous soulignons que l'étude de transport dans les systèmes magnétiques doit tenir compte du paysage énergétique local du réseau magnétique, ainsi que de la dynamique des fluctuations des spins en fonction de la température (ce dont ne tenaient pas compte les travaux ultérieurs). C'est deux points nous permettent un traitement généralisé de la résistivité dans les systèmes non-frustrés. Quant aux systèmes antiferromagnétiques frustrés, ils se caractérisent par un grand nombre d'états dégénérés (souvent infini pour des spins d'Ising). Cependant, la structure en couches minces du système nous permet de rendre cette dégénérescence finie. Nous mettons premièrement en évidence que la dégénérescence du système induit une transition de résistivité du premier ordre, et que selon l'état dégénéré du système, la résistivité peut présenter deux types de transition : une transition allant des basses résistivités vers les hautes résistivités, dans le sens des températures croissantes, ou inversement. Soulignons que ce type de matériau peut présenter un intérêt dans l'élaboration de système tel que les rams assistées par la chaleur.

Le travail que nous présentons dans cette thèse s'articule autour du problème de la formation de coalitions entre des agents égoïstes qui planifient leurs activités, dans les systèmes multi-agents (SMA). Nous avons proposé, dans un premier temps, un mécanisme qui se fonde sur l’analyse des actions des agents dans leurs plans et le raisonnement sur les plans des autres, grâce notamment au calcul d’un degré de croyance sur les actions. Nous nous sommes, par ailleurs, intéressés au problème de la formation de coalitions avec des contraintes dynamiques et des préférences que les agents révèlent et communiquent aux autres lors de leurs négociations. Enfin, nous avons affiné notre mécanisme de formation des coalitions en permettant une recherche des coalitions guidée par la construction d'un arbre de contraintes et d'un arbre de coalitions, qui sont ensuite exploré par le biais de l'algorithme Monte-Carlo. / The work we present, in this thesis, focuses on the coalition formation problem for self-interested agents which plan their activities in multi-agents systems. As a first step, we have proposed, a mechanism that is based on the analysis of the agents' actions in their plans and reasoning about the plans of others. Additionally, we have addressed the problem of coalition formation with dynamic constraints and preferences that agents reveal and communicate to others during their negotiations. Finally, we have refined our coalition formation mechanism allowing a guided search of the coalitions by building a tree of constraints and a tree of coalitions. Each tree is explored by means of the Monte-Carlo algorithm.

Dans cette thèse, on propose des techniques d'inférence bayésienne dans les modèles où la vraisemblance possède une composante latente. La vraisemblance d'un jeu de données observé est l'intégrale de la vraisemblance dite complète sur l'espace de la variable latente. On s'intéresse aux cas où l'espace de la variable latente est de très grande dimension et comportes des directions de différentes natures (discrètes et continues), ce qui rend cette intégrale incalculable. Le champs d'application privilégié de cette thèse est l'inférence dans les modèles de génétique des populations. Pour mener leurs études, les généticiens des populations se basent sur l'information génétique extraite des populations du présent et représente la variable observée. L'information incluant l'histoire spatiale et temporelle de l'espèce considérée est inaccessible en général et représente la composante latente. Notre première contribution dans cette thèse suppose que la vraisemblance peut être évaluée via une approximation numériquement coûteuse. Le schéma d'échantillonnage préférentiel adaptatif et multiple (AMIS pour Adaptive Multiple Importance Sampling) de Cornuet et al. [2012] nécessite peu d'appels au calcul de la vraisemblance et recycle ces évaluations. Cet algorithme approche la loi a posteriori par un système de particules pondérées. Cette technique est conçue pour pouvoir recycler les simulations obtenues par le processus itératif (la construction séquentielle d'une suite de lois d'importance). Dans les nombreux tests numériques effectués sur des modèles de génétique des populations, l'algorithme AMIS a montré des performances numériques très prometteuses en terme de stabilité. Ces propriétés numériques sont particulièrement adéquates pour notre contexte. Toutefois, la question de la convergence des estimateurs obtenus parcette technique reste largement ouverte. Dans cette thèse, nous montrons des résultats de convergence d'une version légèrement modifiée de cet algorithme. Sur des simulations, nous montrons que ses qualités numériques sont identiques à celles du schéma original. Dans la deuxième contribution de cette thèse, on renonce à l'approximation de la vraisemblance et onsupposera seulement que la simulation suivant le modèle (suivant la vraisemblance) est possible. Notre apport est un algorithme ABC séquentiel (Approximate Bayesian Computation). Sur les modèles de la génétique des populations, cette méthode peut se révéler lente lorsqu'on vise uneapproximation précise de la loi a posteriori. L'algorithme que nous proposons est une amélioration de l'algorithme ABC-SMC de DelMoral et al. [2012] que nous optimisons en nombre d'appels aux simulations suivant la vraisemblance, et que nous munissons d'un mécanisme de choix de niveauxd'acceptations auto-calibré. Nous implémentons notre algorithme pour inférer les paramètres d'un scénario évolutif réel et complexe de génétique des populations. Nous montrons que pour la même qualité d'approximation, notre algorithme nécessite deux fois moins de simulations par rapport à laméthode ABC avec acceptation couramment utilisée. / This thesis consists of two parts which can be read independently.The first part is about the Adaptive Multiple Importance Sampling (AMIS) algorithm presented in Cornuet et al.(2012) provides a significant improvement in stability and Effective Sample Size due to the introduction of the recycling procedure. These numerical properties are particularly adapted to the Bayesian paradigm in population genetics where the modelization involves a large number of parameters. However, the consistency of the AMIS estimator remains largely open. In this work, we provide a novel Adaptive Multiple Importance Sampling scheme corresponding to a slight modification of Cornuet et al. (2012) proposition that preserves the above-mentioned improvements. Finally, using limit theorems on triangular arrays of conditionally independant random variables, we give a consistensy result for the final particle system returned by our new scheme.The second part of this thesis lies in ABC paradigm. Approximate Bayesian Computation has been successfully used in population genetics models to bypass the calculation of the likelihood. These algorithms provide an accurate estimator by comparing the observed dataset to a sample of datasets simulated from the model. Although parallelization is easily achieved, computation times for assuring a suitable approximation quality of the posterior distribution are still long. To alleviate this issue, we propose a sequential algorithm adapted fromDel Moral et al. (2012) which runs twice as fast as traditional ABC algorithms. Itsparameters are calibrated to minimize the number of simulations from the model.

Dans les collisions ion-solide, les états de l'ion projectile de grand moment angulaire l sont en moyenne beaucoup plus peuplés que lors des collisions ion-atome. L'utilisation d'ions projectiles Ar18+ d'énergie 13,6 MeV/u et de cibles solides de carbone nous a permis d'étudier les états excités de l'ion Ar17+ peuplés par capture. La gamme d'épaisseur de cibles choisie a permis d'effectuer cette étude de la condition de collision unique (3,5 µg/cm²) jusqu'à l'équilibre des populations (200 µg/cm²). Nous avons observé l'évolution des intensités des transitions Lyman en fonction du temps de transit de l'ion dans la cible (évolution sensible à la population des états de coeur) ainsi qu'en fonction du temps de vol de l'ion derrière la cible (évolution sensible à la population des états de Rydberg). Pour expliquer les résultats expérimentaux, une analyse complète du transport des états excités dans la cible a été réalisée. Nous avons utilisé deux modèles de type collisionnel. Le premier est un modèle d'équations d'évolution basé sur une description statistique des collisions binaires du projectile avec les atomes cibles. Le second est un modèle de transport classique basé sur l'équation décrivant le mouvement de l'électron projectile sur une orbite classique perturbé par une force stochastique. Une comparaison avec l'expérience montre que ces approches collisionnelles permettent de refléter assez bien la population en moment angulaire des états très excités mais nettement moins bien celle des états de coeur. Le mélange l observé pour ces états est beaucoup plus important que prévu par ces modèles. La polarisation du milieu induite par l'ion projectile (non prise en compte par les modèles) pourrait être responsable d'un tel mélange par l'intermédiaire de l'effet Stark.

L'étude présentée porte sur la caractérisation et la comparaison des propriétés quasi-statiques et en fatigue de composites à fibres de lin et de verre avec une matrice époxy. Une importante campagne d'essais de fatigue a été réalisée à partir des propriétés quasi-statiques de traction et cisaillement plan mesurées. Il apparaît que les caractéristiques en statique et en fatigue du verre/époxy sont supérieures à celles du lin/époxy à taux de fibres et stratification identiques. Cependant, en raison de la faible densité des renforts de lin, les écarts entre les propriétés spécifiques des deux matériaux sont faibles. La mesure de l'évolution des propriétés en fatigue a mis à jour des comportements phénoménologiques particuliers tels que la diminution de l'amortissement et l'augmentation de la rigidité des lin/époxy ayant des fibres parallèles à la direction du chargement, au cours de la vie des éprouvettes. Le suivi de l'endommagement par la mesure des densités de fissures montre que celle-ci augmente avec le chargement et le nombre de cycles appliqués. Les simulations par éléments finis de la durée de vie des composites lin/époxy, intégrant la variabilité des paramètres du modèle par la méthode Monté-Carlo, donnent des courbes de Wöhler conservatives par rapport aux données expérimentales. Conformément aux mesures expérimentales, la variabilité des durées de vie calculées diminue lorsque le chargement baisse