La reconstruction des sources de l'activité cérébrale à partir des mesures de potentiel fournies par un électroencéphalographie (EEG) nécessite de résoudre le problème connu sous le nom de « problème inverse de l'EEG ». La solution de ce problème dépend de la solution du « problème direct de l'EEG », qui fournit à partir de sources de courant connues, le potentiel mesuré au niveau des électrodes. Pour des modèles de tête réels, ce problème ne peut être résolut que de manière numérique. En particulier, les équations intégrales de surfaces requièrent uniquement la discrétisation des interfaces entre les différents compartiments constituant le milieu cérébral. Cependant, les formulations intégrales existant actuellement ne prennent pas en comptent l'anisotropie du milieu. Le travail présenté dans cette thèse introduit deux nouvelles formulations intégrales permettant de palier à cette faiblesse. Une formulation indirecte capable de prendre en compte l'anisotropie du cerveau est proposée. Elle est discrétisée à l'aide de fonctions conformes aux propriétés spectrales des opérateurs impliqués. L'effet de cette discrétisation de type mixe lors de la reconstruction des sources cérébrales est aussi étudié. La seconde formulation se concentre sur l'anisotropie due à la matière blanche. Calculer rapidement la solution du système numérique obtenu est aussi très désirable. Le travail est ainsi complémenté d'une preuve de l'applicabilité des stratégies de préconditionnement de type Calderon pour les milieux multicouches. Le théorème proposé est appliqué dans le contexte de la résolution du problème direct de l'EEG. Un préconditionneur de type Calderon est aussi introduit pour l'équation intégrale du champ électrique (EFIE) dans le cas de structures unidimensionnelles. Finalement, des résultats préliminaires sur l'impact d'un solveur rapide direct lors de la résolution rapide du problème direct de l'EEG sont présentés. / Electroencephalography (EEG) is a very useful tool for characterizing epileptic sources. Brain source imaging with EEG necessitates to solve the so-called EEG inverse problem. Its solution depends on the solution of the EEG forward problem that provides from known current sources the potential measured at the electrodes positions. For realistic head shapes, this problem can be solved with different numerical techniques. In particular surface integral equations necessitates to discretize only the interfaces between the brain compartments. However, the existing formulations do not take into account the anisotropy of the media. The work presented in this thesis introduces two new integral formulations to tackle this weakness. An indirect formulation that can handle brain anisotropies is proposed. It is discretized with basis functions conform to the mapping properties of the involved operators. The effect of this mixed discretization on brain source reconstruction is also studied. The second formulation focuses on the white matter fiber anisotropy. Obtaining the solution to the obtained numerical system rapidly is also highly desirable. The work is hence complemented with a proof of the preconditioning effect of Calderon strategies for multilayered media. The proposed theorem is applied in the context of solving the EEG forward problem. A Calderon preconditioner is also introduced for the wire electric field integral equation. Finally, preliminary results on the impact of a fast direct solver in solving the EEG forward problem are presented.

Dans cette thèse, nous introduisons et développons une approche mathématiques originale des techniques d'imagerie biomédicales dites "hybrides". L'idée et d'appliquer une méthode d'imagerie mal posée, tout en perturbant le milieu à imager par des déplacements mécaniques. Ces déplacements provenant d'une équation de type onde élastique perturbent les mesures effectuées. En utilisant ces mesures perturbées, et profitant du caractère local des perturbations mécaniques, il est possible d'augmenter considérablement la résolution de la méthode de base. Le problème direct est donc un couplage d'une EDP décrivant la propagation utilisée pour la méthode de base et d'une seconde décrivant les champs de déplacement mécaniques. Dans toutes cette thèse, on fait l'hypothèse d'un milieu mécaniquement homogène afin d'assurer le contrôle et la géométrie des ondes perturbatrices utilisées. A partir des mesures perturbées, une étape d'interprétation permet de construire une donnée interne au domaine considéré. Cette étape nécessite en général l'inversion d'opérateurs géométriques intégraux de type Radon, afin d'utiliser le caractère localisant des perturbations utilisées. A partir de cette donnée interne, il est possible d'initier une procédure de reconstruction du paramètre physique recherché. Dans le chapitre 1, il est question d'un couplage entre micro-ondes et perturbations sphériques. Dans les chapitres 2, 3 et 4, nous étudions l'imagerie optique diffuse toujours couplée avec des perturbations sphériques. Enfin dans le chapitre cinq, nous donnons une méthode originale de reconstruction de la conductivité électrique par un couplage entre champs magnétique et perturbations acoustiques focalisées. / This thesis aims at developing an original mathematical approach for modeling hybrid biomedical imaging modalities. The core idea is to run an ill-posed imaging method while perturbing the medium using mechanical displacements. These displacements described by an elastic wave equation perturb the collected measurements. Using these perturbed measurements and taking advantage of the perturbation localizing e↵ect, it is possible to significantly overcome the resolution of the basic method. The direct problem here is a coupling between a PDE describing the propagation used for the basic method and a second one describing the mechanical displacements fields. In the whole thesis, we only consider mechanically homogeneous medium in order to assure the control and the geometry of the perturbing wavefronts. From these perturbed measurements, an interpretation step leads to an internal data map inside the considered medium. This step usually requires inversion of geometric integral operators such as Radon transform. This allows to use the geometrical localizing behavior of the perturbations. From this internal data, one can start a recovering procedure for the unknown physical parameter. This recovering step involves a new non physical PDE, non linearly coupled with the main modality equation. In the first chapter, we study a coupling between micro-waves and spherical perturbations. In chapter 2, 3 and 4, we propose a model for di↵use optical imaging coupled with spherical perturbations. In chapter 5, we introduce a new method for imaging the electric conductivity by a coupling between magnetic field and focused acoustic perturbations

Nous étudions des fluctuations qui sont omniprésentes dans des collisions entre particules aux hautes énergies. Ces fluctuations peuvent être de nature classique ou quantique et nous allons considérer ces deux cas. D'abord, nous étudions les fluctuations quantiques qui sont présentes dans des collisions entre protons. Celles-ci sont calculables en théorie quantique des champs, et nous allons nous concentrer sur une certaine classe de diagrammes dans ce cadre. Dans un second temps nous allons étudier des fluctuations qui sont présentes dans des collisions entre particules plus lourdes que le proton. Celles-ci sont décrites par les lois quantiques de la nature qui donnent les positions des nucléons dans le noyau, ou bien des fluctuations classiques, d'origine thermique, qui affectent l'évolution hydrodynamique du milieu produit dans une collision. Les fluctuations dans des collisions entre protons peuvent être calculées analytiquement jusqu'à un certain ordre en théorie quantique des champs. Nous allons nous concentrer sur des diagrammes à une boucle, d'une topologie donnée. Ces diagrammes aux boucles donnent des intégrales, qui typiquement sont difficiles à calculer. Nous allons démontrer comment des outils des mathématiques modernes peuvent être utilisés pour faciliter leur évaluation. En particulier, nous allons étudier des relations entre des coupures d'un diagramme, la discontinuité à travers d'un branchement et le coproduit. Nous allons démontrer comment l'intégrale originale peut être reconstruit à partir de l'information contenue dans le coproduit. Nous nous attendons à ce que ces méthodes seront utiles pour le calcul des diagrammes avec des topologies plus difficiles et ainsi aident au calcul des nouvelles amplitudes de diffusion. A la fin, nous étudions les deux types de fluctuations qui ont lieu dans des collisions entre ions lourds. Celles-ci sont liées soit à l'état initial de la matière, soit à l'état intermédiaire produit dans une telle collision. Les fluctuations de l'état initial ont été mesurées expérimentalement, et on voit qu'elles donnent lieu à des non-Gaussianités dans le spectre final de particules. Nous allons démontrer comment ces non-Gaussianités peuvent être comprises comme des positions et des énergies d'interaction aléatoires des 'sources' dans les noyaux entrant en collision. En plus, nous étudions le bruit hydrodynamique dans le milieu produit juste après une collision. Le comportement de ce milieu est celui d'un fluide à basse viscosité. / We study fluctuations that are omnipresent in high-energy particle collisions. These fluctuations can be either of either classical or quantum origin and we will study both. Firstly, we consider the type of quantum fluctuations that arise in proton-proton collisions. These are computable perturbatively in quantum field theory and we will focus on a specific class of diagrams in this set-up. Secondly, we will consider the fluctuations that are present in collisions between nuclei that can be heavier than protons. These are the quantum laws of nature that describe the positions of nucleons within a nucleus, but also the hydrodynamic fluctuations of classical, thermal origin that affect the evolution of the medium produced in heavy-ion collisions. The fluctuations arising in proton-proton collisions can be computed analytically up to a certain order in perturbative quantum field theory. We will focus on one-loop diagrams of a fixed topology. Loop diagrams give rise to integrals that typically are hard to evaluate. We show how modern mathematical methods can be used to ease their computation. We will study the relations among unitarity cuts of a diagram, the discontinuity across the corresponding branch cut and the coproduct. We show how the original integral corresponding to a given diagram can be reconstructed from the information contained in the coproduct. We expect that these methods can be applied to solve more complicated topologies and help in the computation of new amplitudes in the future. Finally, we study the two types of fluctuations arising in heavy-ion collisions. These are related either to the initial state or the intermediate state of matter produced in such collisions. The initial state fluctuations are experimentally observed to give rise to non-Gaussianities in the final-state spectra. We show how these non-Gaussianities can be explained by the random position and interaction energy of `sources' in the colliding nuclei. Furthermore, we investigate the effect of hydrodynamical noise in the evolution of the medium produced just after a collision. This medium behaves like a fluid with a very low viscosity, and so the corresponding evolution is hydrodynamical.

Cette thèse développe un cadre rigoureux qui permet de démontrer des représentations exactes associées à divers observables de la chaîne XXZ de Heisenberg de spin 1/2 à température finie. Il a était argumenté dans la littérature que l’énergie libre par site ou les longueurs de corrélations admettent des représentations intégrales où les intégrandes sont exprimées en termes de solutions d’équations intégrales non-linéaires. Les dérivations de ces représentations reposaient sur divers conjectures telles que l’existence d’une valeur propre de la matrice de transfert quantique, real, non-dégénérée, de module maximale, de l’échangeabilitée de la limite du volume infinie et du nombre de Trotter à l’infinie, de l’existence et de l’unicité des solutions des equation intégrales non-linéaires auxiliaires et finalement de l’identification des valeurs propers de la matrice de transfert quantiques avec les solutions de l’équations intégrales non-linéaires. Nous démontrons toutes ces conjectures dans le regime de haute température. Nôtre analyse nous permet aussi de démontrer que pour ces température suffisamment élevées, il est possible d’avoir une description d’un certain sous-ensemble de valeurs propres sous-dominante de la matrice de transfert quantique décrite en terme de solutions d’une chaîne de spin-1 de taille finie. / This thesis develops a rigorous framework allowing one to prove the exact representations for various observables in the XXZ Heisenberg spin-1/2 chain at finite temperature. Previously it has been argued in the literature that the per-site free energy or the correlation lengths admit integral representations whose integrands are expressed in terms of solutions of non-linear integral equations. The derivations of such representations relied on various conjectures such as the existence of a real, non-degenerate, maximal in modulus Eigenvalue of the quantum transfer matrix, the exchangeability of the infinite volume limit and the Trotter number limits, the existence and uniqueness of the solutions to the auxiliary non-linear integral equations and finally the identification of the quantum transfer matrix’s Eigenvalues with solutions to the non-linear integral equation. We rigorously prove all these conjectures in the high temperature regime. Our analysis also allows us to prove that for temperatures high enough, one may describe a certain subset of sub-dominant Eigenvalues of the quantum transfer matrix described in terms of solutions to a spin-1 chain of finite length.

Les propriétés de solvatation jouent un rôle important dans les problèmes chimiques et biochimiques. La théorie fonctionnelle de la densité moléculaire (MDFT) est l'une des méthodes frontières pour évaluer ces propriétés, dans laquelle une fonction d'énergie libre de solvatation est minimisée pour un soluté arbitraire dans une boîte de solvant cubique périodique. Dans cette thèse, nous travaillons sur l'évaluation du terme d'excès de la fonctionnelle d’énergie libre sous l’approximation du fluide de référence homogène (HRF), équivalent à l'approximation de la chaîne hypernettée (HNC) dans la théorie des équations intégrales. Deux algorithmes sont proposés: le premier est une extension d'un algorithme précédent, qui permet de traiter le cas d'un solvant moléculaire à trois dimensions (en fonction de trois angles d'Euler) au lieu d'un solvant linéaire (selon deux angles); L'autre est un nouvel algorithme qui intègre le traitement de la convolution angulaire de l'équation Ornstein-Zernike (OZ) moléculaire dans MDFT, et en fait développe la densité du solvant et le gradient fonctionnel en harmoniques sphériques généralisées (GSHs). On montre que le nouvel algorithme est beaucoup plus rapide que le précédent. Les deux algorithmes sont appropriés pour des solutés arbitraires tridimensionnel dans l'eau liquide, et pour prédire l'énergie libre et la structure de solvatation d'ions et de molécules. / Solvation properties play an important role in chemical and bio-chemical issues. The molecular density functional theory (MDFT) is one of the frontier numerical methods to evaluate these properties, in which the solvation free energy functional is minimized for an arbitrary solute in a periodic cubic solvent box. In this thesis, we work on the evaluation of the excess term of the free energy functional under the homogeneous reference fluid (HRF) approximation, which is equivalent to hypernetted-chain (HNC) approximation in integral equation theory. Two algorithms are proposed: the first one is an extension of a previously implemented algorithm, which makes it possible to handle full 3D molecular solvent (depending on three Euler angles) instead of linear solvent (depending on two angles); the other one is a new algorithm that integrates the molecular Ornstein-Zernike (OZ) equation treatment of angular convolution into MDFT, which in fact expands the solvent density and the functional gradient on generalized spherical harmonics (GSHs). It is shown that the new algorithm is much more rapid than the previous one. Both algorithms are suitable for arbitrary three-dimensional solute in liquid water, and are able to predict the solvation free energy and structure of ions and molecules.

Le Contrôle Non Destructif (CND) permet de sonder l’intérieur d’un milieu dans le but desurveiller son intégrité et son vieillissement. Assisté d’outils de simulation il permet de détecter, caractériseret localiser des défauts de structure du milieu inspecté mais sa fiabilité dépend de la précision des méthodesde simulation. Dans le cadre du CND ultrasonore, les méthodes usuelles (numériques et asymptotiques) sontbien souvent inadéquates pour simuler la diffraction par les défauts. On leur préfère des techniques hybrides.On propose dans cette thèse une nouvelle approche hybride pour la simulation numérique de la diffractionhaute fréquence en milieu étendu (configuration critique pour le CND). Combinant la méthode des équationsintégrales et la méthode des rayons, cette approche exploite le caractère multi-échelle du problème hautefréquence en proposant un modèle d’obstacle à deux échelles. Elle permet le calcul précis de la diffraction etla propagation rayon des champs. D’abord mise au point dans le cadre de la diffraction d’ondes acoustiquespar un obstacle de taille inférieure à la longueur d’onde (méthode barycentrique), l’approche est ensuiteétendue à des configurations de diffraction par des obstacles de l’ordre de la longueur d’onde grâce àl’introduction d’un partitionnement de l’unité de sa surface (méthode multi-centres). Pour accélérerl’approche hybride, on propose une procédure de résolution Online-Offline, basée sur un pré-calcul de lamatrice de diffraction associée à un ensemble réduit de directions d’incidence et d’observation et sur uneinterpolation polynomiale de ses vecteurs singuliers pour son évaluation dans des directions quelconquesd’émetteurs et de récepteurs. On étudie ensuite la stratégie dans le cadre de l’acoustique 3D puis on en faitune extension de principe à l’élastodynamique. On donne enfin un ensemble de pistes pour étendre l’approchehybride dans des cas de diffraction par un ou plusieurs obstacles pouvant être proches des bords du milieu. / Non Destructive Testing (NDT) aims at probing a medium to check its integrity and aging. Withthe help of simulation tools, it allows to detect, caracterize and locate flaws inside a material with a precisiondepending on that of the simulation methods. Usual numerical and asymptotic methods nevertheless often failat precisely computing diffraction for ultrasonic NDT. Hybrid approaches are thus prefered in this framework.In this thesis, we propose a new hybrid strategy combining the boundary integral equation method and raytracing to compute high frequency diffraction of an obstacle in a large medium (critical NDT configuration).This strategy allows to compute precisely the diffraction effects and to convert and propagate the diffractedfield as rays. The proposed strategy relies on a two-scale model of the diffracting obstacle. First developpedto simulate acoustic waves diffraction on an obstacle of size less than the wave length (barycentric method),the hybrid strategy is then extended to compute the diffraction by an obstacle of size some wave lengths(polycentric method) thanks to the introduction of a partition of unity of the obstacle surface. Besides, inorder to accelerate the hybrid approach, we propose an Online-Offline resolution procedure based on theOffline computation of the scattering matrix for a reduced set of incidence and observation directions and onthe use of a polynomial interpolation of its singular vectors for the Online evaluation of the scattering matrixfor any incidence and observation directions. We then study the possibility of extension of the hybrid strategyto 3D acoutics and elastodynamics. We finally give several perspectives for the adaptation of the approach todeal with diffraction by one or several obstacles potentially close to the propagating medium boundaries.

Le mécanisme permettant d’expliquer l’origine de la masse des neutrinos demeure, encore aujourd’hui, un mystère complet dont la résolution est susceptible de modifier considérablement la structure du modèle standard en physique des particules élémentaires. Dans la littérature, plusieurs candidats potentiels sont donc proposés afin de combler cette lacune et, ainsi, faire la lumière sur certaines des propriétés les plus étranges des neutrinos. Parmi ceux-ci, les mécanismes seesaw de type I, II et III constituent sans doute les approches les plus attrayantes et les plus étudiées. Cependant, bien que ces mécanismes offrent un cadre de travail simple et élégant pour expliquer la faible masse des neutrinos (l’ordre de grandeur), ceux-ci n’offrent aucune prédiction sur les paramètres fondamentaux caractérisant le phénomène d’oscillation, soit les angles de mélange, les phases complexes et la hiérarchie des masses. Afin d’obtenir des prédictions concrètes sur la phénoménologie des neutrinos, certaines hypothèses de travail supplémentaires doivent donc être formulées pour contraindre la structure des matrices de masse obtenue. Dans ce travail, l’hypothèse anarchique propre au secteur des neutrinos est adoptée. Les matrices de masse générées par les trois mécanismes seesaw dans la limite des basses énergies sont traitées dans le contexte d’une théorie de matrices aléatoires, ce qui permet de définir et d’analyser de nouveaux ensembles matriciels aléatoires appelés ensembles seesaw. Un cadre théorique unifié est donc présenté pour la construction de ces ensembles. Grâce au formalisme élaboré, qui repose sur les outils traditionnels relevant de la théorie des matrices aléatoires, les densités de probabilité jointes caractérisant ces ensembles sont obtenues de façon analytique. Une étude détaillée de leurs propriétés est alors réalisée, ce qui permet d’extraire les tendances dominantes propres à ces mécanismes de masse et d’analyser leurs conséquences pour le secteur des neutrinos du modèle standard étendu. En ce qui concerne le spectre de masse, les résultats obtenus indiquent que les mécanismes seesaw de type I et de type III sont plus adéquats pour reproduire les observations expérimentales. De plus, une forte préférence pour la différence de masses associée à la hiérarchie normale est observée. En contrepartie, il est également démontré que pour une différence de masses donnée entre les trois générations, toutes les permutations des masses sont équiprobables, ce qui rend hors de portée toute prédiction concernant la hiérarchie du spectre (normale ou inverse) sous l’hypothèse anarchique. En ce qui concerne les variables du groupe de symétrie (les angles de mélange et les phases complexes), on constate, d’une part, que la notion de mélange quasi-maximal est naturellement favorisée et, d’autre part, que la matrice PMNS peut être décrite comme une matrice unitaire générique tirée au hasard d’un ensemble matriciel caractérisé par la mesure de Haar du groupe de Lie correspondant. Par ailleurs, il est également démontré que ces conclusions sont indépendantes du mécanisme de masse considéré. / The neutrino mass generation mechanism remains, to this day, a complete mystery which is likely to play an important role in understanding the foundations of the Standard Model of particle physics. In an effort to fill this gap and, ultimately, shed some light on some of the most intriguing properties of neutrinos, many theoretical models are proposed in the literature. Among the many candidates, the type I, type II and type III seesaw mechanisms may very well be the most attractive and the most studied propositions. However, despite the fact that these mechanisms provide a simple and elegant framework for explaining the smallness of neutrino masses (the order of magnitude), no prediction can be made on the fundamental parameters governing neutrino oscillations (the mixing angles, the CP-violating phases and the mass differences). Thus, to obtain concrete results regarding neutrino phenomenology, additional working assumptions must be made in order to constrain the structure of the corresponding mass matrices. In this work, the anarchy hypothesis relevant to the neutrino sector is investigated. The mass matrices generated by the three seesaw mechanisms in the low-energy limit are studied within the framework of random matrix theory, which leads to the development and the analysis of the seesaw ensembles. A unified and precise theoretical formalism, based on the usual tools of random matrix theory, is presented for the construction of these new random matrix ensembles. Using this formalism, the joint probability density functions characterizing these ensembles are obtained analytically, thus paving the way for a detailed study of their properties. This study is then carried out, revealing the underlying trends in these ensembles and, thereby, offering a thorough analysis of their consequences for the neutrino sector of the seesaw-extended Standard Model. Regarding the mass spectrum, it is found that the type I and type III seesaw mechanisms are better suited to accommodate experimental data. Moreover, the results indicate a strong preference for the mass splitting associated to normal hierarchy. However, since all permutations of the masses are found to be equally probable for a particular mass splitting between the three generations, predictions concerning the hierarchy of the mass spectrum (normal or inverted) remains out of reach in the framework of anarchy. Regarding the group variables (the mixing angles an CP-violating phases), it is found that near-maximal mixing is naturally favored by these ensembles and, that the PMNS matrix can be described as a generic unitary matrix drawn at random from a matrix ensemble characterized by the Haar measure of the corresponding Lie group. Furthermore, these conclusions are found to be independent of the mass mechanism considered.

Les intégrales bi-électroniques jouent un rôle essentiel pour déterminer les propriétés électroniques des atomes et des molécules. Les oxydes de cuivre couvrent des domaines très divers et ont des applications extrêmement variés. Nous avons porté notre attention sur une étude de la contribution des intégrales bi-électroniques à la catalyse hétérogène, par une étude théorique de la réactivité sur la surface de cuivre, et en particulier l'oxydation ménagée des alcanes et alcènes sur la surface CuO(001). Ce travail de thèse a porté sur une étude théorique en deux parties principales : La première partie concernant l'étude des intégrales de Coulomb et d'échange effectuée par le logiciel STOP sur une base d'orbitales atomiques de Slater. L'équation de Poisson est résolue dans ce travail pour le potentiel de Coulomb moléculaire. La deuxième partie, est une étude des réactivités chimiques concernant l'oxydation ménagée des alcanes et alcènes sur les surfaces de métal de cuivre Cu(001) ou CuO(001). Cette étude est effectuée à l'aide du logiciel ABINIT DFT sur une base d'onde plane. Les résultats de ce travail ont validé une nouvelle approche des intégrales de Coulomb et d'échange. Ils ont abouti à la mise en évidence des mécanismes réactionnels de la réactivité d'oxydation ménagée sur la surface de cuivre CuO(001), ainsi que des produits sélectifs.

Nous présentons différents résultats liés par les outils et les structures qu'ils font intervenir (intégrales itérées, produits de mélange). Dans la première partie, nous considérons le calcul de certaines intégrales de type Selberg et leurs limites lorsque le nombre de variables tend vers l'infini. Dans le cas général, on montre que le résultat s'exprime comme un produit dont le nombre de facteurs ne dépend pas du nombre de variables (sous certaines conditions). Si la puissance du déterminant de Vandermonde vaut 2, il est possible de calculer la limite de ces intégrales lorsque le nombre de variables tend vers l'infini à l'aide d'opérateurs liés à l'interpolation de Newton. Dans la seconde partie, nous étudions les propriétés de dépendance linéaire de familles de fonctions obtenues par intégrales itérées et donnons un critère qui permet d'assurer l'indépendance linéaire d'une famille infinie de fonctions à partir de l'étude des relations entre les fonctions obtenues par intégrales simples. Nous montrons comment construire effectivement les corps de germes de fonctions analytiques nécessaires et en donnons quelques exemples qui permettent d'étendre les résultats connus sur les hyperlogarithmes. Ensuite, nous étudions certaines bases de l'algèbre libre dans le but d'appliquer la factorisation de Schützenberger. Nous rappelons quelques résultats classiques, puis nous intéressons à la famille obtenue à partir des mots de Lyndon. Celle-ci ne permet pas d'écrire la factorisation qui nous intéresse mais nous précisons les caractéristiques de sa famille duale. Enfin, nous donnons un critère relatif à deux familles en dualité assurant que l'on peut écrire cette factorisation.

Les calculs de précision dans les théories de jauge jouent un rôle très important pour l'étude de la physique du Modèle Standard et au-delà dans les super-collisionneurs de particules comme le LHC, TeVatron et ILC. Par conséquent, il est extrêmement important de fournir des outils du calcul d'amplitudes à une boucle stables, rapides, efficaces et hautement automatisés. Cette thèse a pour but de développer la bibliothèque d'intégrales Golem95. Cette bibliothèque est un programme écrit en Fortran95, qui contient tous les ingrédients nécessaires pour calculer une intégrale scalaire ou tensorielle à une boucle avec jusqu'à six pattes externes. Golem95 utilise une méthode traditionnelle de réduction (réduction à la Golem) qui réduit les facteurs de forme en des intégrales de base redondantes qui peuvent être scalaires (sans paramètres de Feynman au numérateur) ou tensorielles (avec des paramètres de Feynman au numérateur); ce formalisme permet d'éviter les problèmes de l'instabilité numérique engendrés par des singularités factices dues à l'annulation des déterminants de Gram. En plus, cette bibliothèque peut être interfacée avec des programmes du calcul automatique basés sur les méthodes d'unitarité comme GoSam par exemple. Les versions antérieures de Golem95 ont été conçues pour le calcul des amplitudes sans masses internes. Le but de ce travail de thèse est de généraliser cette bibliothèque pour les configurations les plus générales (les masses complexes sont incluses), et de fournir un calcul numériquement stable dans les régions problématique en donnant une représentation intégrale unidimensionnelle stable pour chaque intégrale de base de Golem95.