L’étude de la dynamique d’un faisceau d’ions de haute intensité dans les lignes de basse énergie (LBE)représente l’un des défis majeurs de la science des accélérateurs. basse énergie, cette dynamique est dominée par le champ de charge d’espace induit par le faisceau lui-même, qui en général est non linéaire et peut entrainer des phénomènes de halo, de grossissement d’émittance et de pertes de faisceau. Toutefois, un faisceau à basse énergie se propageant dans une LBE induit l’ionisation du gaz résiduel présent dans la chambre. Les particules (ions et électrons) issues de l’ionisation sont repoussées ou attirées radialement par le champ de charge d’espace en fonction du signe de leur charge. D’autres réactions physiques ont lieu dans la ligne basse énergie, jouant ainsi un rôle dans la dynamique du faisceau et sur l’établissement du temps et du taux de compensation de charge d’espace. Afin d’obtenir des résultats prédictifs et fiables quantitativement, des simulations de transport de faisceau en régime de compensation de charge d’espace avec le code de calcul warp ont été réalisées en prenant également compte les réactions physiques les plus probables. On discutera ensuite des résultats de ces simulations en lien avec ceux issus de différentes compagnes de mesure réalisées auprès des LBE des projets MYRRHA et IFMIF. / The study of intense ion beam dynamics in lowenergy beam transport line (LEBT) representsone of the most important challenges inaccelerating sciences. At low energy, it isdominated by the space-charge field created bythe beam itself, which is generally non-linearand can induce halo, emittance growth and beamlosses. But, a ion beam at low energypropagating in a LEBT ionises the residual gas.The particles (ions and electrons) fromionisation are repelled or confined radially bythe space charge field according to their chargesign.Other interactions take place in the LEBT,modifying the beam dynamics and the space chargecompensation time and the space-chargecompensation yield. In order to obtain predictiveand precise results quantitatively, numericalsimulations of beam transport in space-chargecompensation regime with WARP code havebeen realized taking account the most probablephysical interactions. Then, we will discuss theresults with comparisons with experimental dataobtained on the MYRRHA and IFMIF LEBT’s

L’accélération de particules par interaction laser-plasma est une alternative prometteuse aux accélérateurs conventionnels qui permettrait de rendre plus compactes les machines du futur dédiées à la protonthérapie. Des champs électriques extrêmes de l’ordre du TV/m sont créés en focalisant une impulsion laser ultra-intense sur une cible solide mince de quelques micromètres d’épaisseur, ce qui produit un faisceau de particules de haute énergie. Ce dernier contient des protons ayant une énergie allant jusqu’à la dizaine de mégaélectron-volts, et est caractérisé par une forte divergence angulaire et un spectre en énergie très étendu.Le but de cette thèse est de caractériser parfaitement un accélérateur laser-plasma afin de produire un faisceau de protons stable, satisfaisant les critères d'énergie, de charge et d'homogénéité de surface requis pour son utilisation en radiobiologie. La conception, la réalisation et l’implémentation d’un système magnétique, constitué d'aimants permanents quadripolaires ont été optimisés au préalable avec des simulations numériques. Ce système permet d’obtenir un faisceau de protons ayant un spectre en énergie qui à été mise en forme, et dont le profil est uniforme sur une surface de taille adaptée aux échantillons biologiques.Une dosimétrie absolue et en ligne a également été établie, permettant le contrôle de la dose délivrée en sortie. Pour cela, une chambre d'ionisation à transmission, précédemment calibrée sur un accélérateur à usage médical de type cyclotron, a été mise en place sur le trajet du faisceau de protons. Des simulations Monte Carlo ont ensuite permis de calculer la dose déposée dans les échantillons. Ce système compact autorise maintenant de définir un protocole expérimental rigoureux pour la poursuite d’expériences in vitro de radiobiologie. De premières irradiations de cellules cancéreuses ont été ainsi réalisées in vitro, ouvrant la voie à l’exploration des effets de rayonnements ionisants pulsés à haut débit de dose sur les cellules vivantes. / Particle acceleration by laser-plasma interaction is a promising alternative to conventional accelerators that could make future devices dedicated to protontherapy more compact. Extreme electric fields in the order of TV/m are created when an ultra-intense laser pulse is focused on a thin solid target with a thickness of a few micrometers, which generates a beam of highly energetic particles. The latter includes protons with energies up to about ten megaelectron-volts and characterised by a wide angular divergence and a broad energy spectrum.The goal of this thesis is to fully characterise a laser-based accelerator in order to produce a stable proton beam meeting the energy, charge and surface homogeneity requirements for radiobiological experiments. The design, realisation and implementation of a magnetic system made of permanent magnet quadrupoles were optimised beforehand through numerical simulations. It enables to obtain a beam with a shaped energy spectrum and with a uniform profile over a surface with a size adapted to the biological samples.Deferred and online dosimetry was setup to monitor the delivered output dose. For that purpose, a transmission ionisation chamber, previously calibrated absolutely on a medical proton accelerator, was used. Monte Carlo simulations enabled to compute the dose deposited into the samples. This compact system allows now to define a rigorous experimental protocol for in vitro radiobiological experiments. First experiments of cancer cell irradiation have been carried out, paving the way for the exploration of the effects of pulsed ionizing radiations at extremely high dose rates on living cells.

Dans ce projet, nous présentons une méthode nouvelle et précise d’estimation de la direction et des délais d’arrivée dans un environnement à trajets multiples, à des fins d’estimation de canal. Récemment, les méthodes de super-résolution ont été largement utilisées pour l’estimation à haute-résolution de la direction d’arrivée (DOA) ou de la différence de temps d’arrivée (TDOA). L’algorithme proposé dans ce travail est applicable à l’estimation d’un canal espace-temps pour des systèmes de traitement spatio-temporel qui emploient la technologie hybride DOA / TDOA. L’estimateur est basé sur l’algorithme MUSIC classique pour trouver la DOA et en profitant d’un simple corrélateur, il est possible de trouver le retard de chaque arrivée. Il est pertinent d’associer chaque angle à son propre retard pour être capable d’estimer les caractéristiques du canal quand nous ne connaissons pas la séquence transmise par l’émetteur. Pour ce faire, nous proposons une formation de faisceaux (voix) très simple et optimale par l’application du MVDR (Maximum Variance Distortion-less Response). Cette formation de faisceaux maximise le signal desiré par rapport aux autres signaux. Après détermination de l’angle d’arrivée par l’algorithme MUSIC, nous appliquons l’algorithme de formation de faisceaux MVDR pour obtenir le signal qui est reçu par le réseau d’antennes pour une direction. Ce signal est corrélé avec les autres signaux correspondants aux autres directions d’arrivée. Les pics dans les figures ainsi obtenues montrent le décalage temporel de chaque source par rapport à celle obtenue par la formation de faisceaux MVDR. La soustraction du plus petit décalage, correspondant au premier signal reçu à chaque décalage temporel, nous donne le temps d’arrivée de chaque source. Pour être plus précis, nous pouvons choisir la moyenne des vecteurs des délais estimés, chacun étant obtenu à partir d’une angle pour l’algorithme MVDR. / In this project, we present a novel and precise way of estimating the direction and delay of arrivals in multipath environment for channel estimation purposes. Recently, super-resolution methods have been widely used for high resolution Direction Of Arrival (DOA) or Time Difference Of Arrival (TDOA) estimation. The proposed algorithm in this work is applicable to space-time channel estimation for space-time processing systems that employ hybrid DOA/TDOA technology. The estimator is based on the conventional MUSIC algorithm to find the DOA and by using a simple correlator it is possible to find the delay of each arrival. It is of interest to associate each angle to its proper delay to be able to estimate the characteristics of the channel when we have no knowledge about the transmitted sequence. To do this, we suggest a very simple and optimal beamforming method by performing Maximum Variance Distortion-less Response (MVDR). This beamforming maximizes the desired signal in the desired direction compare to the other signals that come from other directions. After finding the DOAs by MUSIC algorithm and selecting our desired direction, we obtain the signal from this direction by applying MVDR beamforming. Then, we perform a correlation between this signal and the others incoming signals from other directions. The peaks in the simulation figures illustrate the delay between each source with the obtained signal from MVDR. If we subtract the delay of the first arrival (the smallest delay in time), from the delays indicated in the figures, we can obtain the delay of each arrival. To be more precise, the mean of these estimated TOAs vector follows the exact TOA of each source.

La très haute luminosité visée par les futurs collisionneurs linéaires nécessite une très forte focalisation finale des faisceaux au point d’interaction jusqu’à des dimensions transverses nanométriques. Dans le cadre de ce travail, la ligne de haute énergie de l’ILC délivrant le faisceau au point d’interaction a d’abord été optimisée pour permettre le ‘push-pull’ des détecteurs, ainsi que pour étudier l’impact d’une réduction d’une centaine de mètre de la longueur totale de la ligne. L’objet du travail a ensuite consisté à optimiser la région d’interaction pour conserver la luminosité en présence du détecteur de l’expérience contenant un solénoïde et un dipôle. Dans ce but, un modèle de la région d’interaction a été établi, afin d’être en mesure de simuler le transport du faisceau dans l’ensemble de la ligne de haute énergie en intégrant les éléments non coaxiaux du détecteur. Cette modélisation inclut pour la première fois tous les éléments électromagnétiques de la région d’interaction (cavité crabe, quadripôles et sextupôles du système de focalisation, solénoïde, dipôle intégré au détecteur). Elle a permis l’optimisation de l’anti-solénoïde, élément essentiel du système de correction des effets du solénoïde de l’expérience. Pour mesurer les performances de la machine après compensation totale des effets du solénoïde, un outil de calcul de la luminosité apte à utiliser des distributions quelconques a été développé. On montre alors que l’acceptance en moment de la ligne est réduite en présence du solénoïde compensé. Il a de plus été mis en évidence que l’insertion du solénoïde induit le transfert de l’effet de la cavité crabe du plan horizontal vers le plan vertical, ce qui provoque une nouvelle perte de luminosité. Enfin la dernière partie de ce travail de thèse est consacrée à l’étude de l’application d’une configuration à grand angle de Piwinski aux collisionneurs linéaires. Pour cela les paramètres des effets faisceau-faisceau en présence d’un angle de croisement ont été évalués. Il est possible de réduire la disruption du faisceau après collision en conservant la luminosité, en revanche réduire le paramètre de beamstrahlung est moins aisé en raison de la déviation horizontale de la trajectoire centrale. / Strong focalisation of the beam is mandatory at the interaction point of the future linear collider in order to reach very high luminosity. In this work, the ILC (International Linear Collider) beam delivery system has been re-optimised, first to take the ‘push-pull’ of the two detectors into account, then to evaluate the influence of a reduction of the total length by hundred meters approximately. In the following part, the interaction region has been optimised to restore the nominal luminosity in the presence of the detector, containing a solenoid and a magnetic dipole. Due to the crossing angle of the beams, these elements are not coaxial and a model for the interaction region had to be developed. This model enables to track the beam in the entire beam delivery system, from the end of the linac to the interaction point. The simulation includes for the first time all the electromagnetic elements of the interaction region (crab cavity, final focusing system quadrupoles and sextupoles, solenoid, detector integrated dipole). Thanks to this model, the weak anti-solenoid could be added and optimised as the main corrector of the solenoid effects on the beam. To study the new performances of the collider after full compensation of these effects, a luminosity calculation tool has been developed. It is shown that the momentum acceptance is reduced after compensation of the solenoid effects. Moreover, transverse coupling induces the transfer of the crab cavity horizontal kick to the vertical plane, implying a significant luminosity loss. Finally the last part of the thesis concerns the application of a large Piwinski’s angle to the linear colliders. The calculation of the beam-beam interaction parameters in the presence of a crossing angle is studied. Due to the crossing angle the central trajectory is deviated in the horizontal plane, preventing the beamstrahlung to be reduced at constant luminosity. However the disruption could be made significantly smaller.

De nombreux efforts sont actuellement portés sur le développement de sources laser à fibre optique, qui ont des applications dans une variété de domaines scientifiques et industriels. Les propriétés intrinsèques de la fibre dopée sont favorables à l'obtention de fortes puissances moyennes, grâce à une excellente capacité de dissipation thermique. Cependant, le confinement de la lumière dans le coeur de la fibre et les grandes longueurs d'interaction sont à l'origine d'effets non linéaires, comme la diffusion Brillouin stimulée dans le cas d'une source de faible largeur spectrale. L'objectif de ce travail de thèse est de développer une source laser fibrée de forte énergie alliant bonne qualité spatiale et finesse spectrale, dans le but de répondre à un besoin de LIDAR à détection cohérente. L'approche suivie consiste à séparer l'obtention de l'énergie de celle de la qualité de faisceau. Avec une chaîne amplificatrice dont le dernier étage est une fibre multimode co-dopée Erbium-Ytterbium, des impulsions laser à la longueur d'onde de 1,5 μm avec une énergie supérieure à 1 mJ ont été obtenues. Le très large coeur de la fibre évite l'apparition de l'effet Brillouin, mais le profil transverse du faisceau est fortement multimode. Pour retrouver une bonne qualité de faisceau, un étage de remise en forme du profil est mis en place, mettant à profit la propriété de nettoyage de faisceau par effet Brillouin dans une fibre multimode àgradient d'indice. Par ailleurs, un schéma original d'amplification Brillouin est proposé et validé expérimentalement, permettant de conserver la largeur spectrale de la source.

En radiothérapie, la tomodensitométrie (CT) fournit l’information anatomique du patient utile au calcul de dose durant la planification de traitement. Afin de considérer la composition hétérogène des tissus, des techniques de calcul telles que la méthode Monte Carlo sont nécessaires pour calculer la dose de manière exacte. L’importation des images CT dans un tel calcul exige que chaque voxel exprimé en unité Hounsfield (HU) soit converti en une valeur physique telle que la densité électronique (ED). Cette conversion est habituellement effectuée à l’aide d’une courbe d’étalonnage HU-ED. Une anomalie ou artefact qui apparaît dans une image CT avant l’étalonnage est susceptible d’assigner un mauvais tissu à un voxel. Ces erreurs peuvent causer une perte cruciale de fiabilité du calcul de dose. Ce travail vise à attribuer une valeur exacte aux voxels d’images CT afin d’assurer la fiabilité des calculs de dose durant la planification de traitement en radiothérapie. Pour y parvenir, une étude est réalisée sur les artefacts qui sont reproduits par simulation Monte Carlo. Pour réduire le temps de calcul, les simulations sont parallélisées et transposées sur un superordinateur. Une étude de sensibilité des nombres HU en présence d’artefacts est ensuite réalisée par une analyse statistique des histogrammes. À l’origine de nombreux artefacts, le durcissement de faisceau est étudié davantage. Une revue sur l’état de l’art en matière de correction du durcissement de faisceau est présentée suivi d’une démonstration explicite d’une correction empirique. / Computed tomography (CT) is widely used in radiotherapy to acquire patient-specific data for an accurate dose calculation in radiotherapy treatment planning. To consider the composition of heterogeneous tissues, calculation techniques such as Monte Carlo method are needed to compute an exact dose distribution. To use CT images with dose calculation algorithms, all voxel values, expressed in Hounsfield unit (HU), must be converted into relevant physical parameters such as the electron density (ED). This conversion is typically accomplished by means of a HU-ED calibration curve. Any discrepancy (or artifact) that appears in the reconstructed CT image prior to calibration is susceptible to yield wrongly-assigned tissues. Such tissue misassignment may crucially decrease the reliability of dose calculation. The aim of this work is to assign exact physical values to CT image voxels to insure the reliability of dose calculation in radiotherapy treatment planning. To achieve this, origins of CT artifacts are first studied using Monte Carlo simulations. Such simulations require a lot of computational time and were parallelized to run efficiently on a supercomputer. An sensitivity study on HU uncertainties due to CT artifacts is then performed using statistical analysis of the image histograms. Beam hardening effect appears to be the origin of several artifacts and is specifically addressed. Finally, a review on the state of the art in beam hardening correction is presented and an empirical correction is exposed in detail.

Le but de ce travail de thèse est d’étudier les moyens de lithographie haute résolution pour la fabrication de résonateurs à ondes élastiques de surface, et de l’illustrer à travers la réalisation de sources de fréquences fonctionnant au-delà du GigaHertz. Dans un premier temps nous abordons les différents dispositifs fondés sur les ondes élastiques de surface puis les sources de fréquence (instabilités caractéristiques) et fixons les objectifs de l’étude au travers notamment d’un état de l’art. Dans un second temps, nous présentons les moyens de lithographie étudiés dans le cadre de ce travail, à savoir la lithographie électronique, la gravure par faisceau d’ions focalisés, la lithographie UV par projection et la lithographie par nano-impression. Pour chacune d’entre elles, nous détaillons le principe de fonctionnement et montrons, notamment grâce à des simulations, leur intérêt et leurs limitations. Ensuite, nous présentons la fabrication et la caractérisation de résonateurs sur différents types de substrats ayant des propriétés innovantes par rapport à nos applications. Le PZT élaboré par épitaxie, présentant des coefficients de couplage élevés (plusieurs pourcents) couplés à une granulométrie fine et une orientation cristalline selon l’axe 001. Le diamant, qui permet d’atteindre des vitesses de phase de l’ordre de 10000 m.s-1, soit une vitesse deux fois supérieure à celles des ondes transverses sur substrat de quartz, quartz que nous avons également étudié afin de rechercher de nouveaux points de fonctionnement à haute fréquence. Pour chaque matériau, nous identifions un ou plusieurs moyen(s) de lithographie qui nous permettent de fabriquer les résonateurs. Les étapes de conception, de fabrication et de caractérisation sont décrites en détail. La dernière partie du manuscrit consiste à exposer les caractéristiques des oscillateurs fondés sur les résonateurs à haut produit Qf ainsi fabriqués (Qf > 5.1012). Nous reportons les résultats obtenus à des fréquences de 1,5 GHz (sur quartz) et à 3 GHz (sur diamant nanocristallin). Le bruit de phase à 10 kHz de la porteuse est compris entre -100 et -110 dBc.Hz-1, et le bruit plancher est de -160 dBc.Hz-1. Nous concluons en donnant des pistes afin d’améliorer ces caractéristiques / The purpose of this thesis is to study high-resolution lithography for the fabrication of surface acoustic wave resonators, and to illustrate this technology through the realization of frequency sources operating beyond the Gigahertz. At first, we detail several devices based on surface acoustic waves and frequency sources (instability characteristics) and set the goals of the study in particular through a state of the art. In a second step, we present the lithography methods studied in this work : electron beam lithography, focused ion beam, UV lithography (stepper) and nano-imprint lithography. For each, we detail the operating principle and show, in particular through simulations, their interest and limitations. Then, we present the fabrication and characterization of resonators on different types of substrates with innovative properties compared to our applications. The epitaxial PZT exhibits high coupling coefficients (several percent) coupled with a fine particle size and crystal orientation along the axis 001. The diamond, which achieves phase velocities of about 10000 m.s-1, twice higher than those of STW waves on quartz substrate, quartz that we have also studied in order to search new operating points at high frequency. For each material, we identify one or more lithography methods that allow to manufacture the resonators. Design, fabrication and characterization steps are described in detail. The last part of the manuscript describes the characteristics of oscillators based on high Qf resonators (Qf > 5.1012). We report the results for operating frequencies of 1.5 GHz (on quartz) and 3 GHz (on nanocrystalline diamond). The phase noise at 10 kHz from the carrier is between -100 and -110 dBc.Hz-1, and the noise floor is -160 dBc.Hz-1. We conclude by giving ideas to improve these characteristics

Les nanofils (NFs), de par leur propriétés opto et nanoélectroniques sont devenus des éléments indispensable à la fabrication des dispositifs de la nanoélectronique. Le problème principal reste la reproductibilité en terme de densité de NFs, de diamètre... Cette thèse a pour but de développer grâce à la technologie FIB, un procédé permettant l'élaboration de NFs organisés et homogène en taille. / Due to their ease of fabrication and unique physical properties, semiconductor nanowires (NWs) have been proposed as building blocks for new nanoelectronic and photonic devices. Various processes have been developed to obtain large density of ultra-small NWs but naturally forms nanowires often lack reproducibility. We propose to develop a bottom-up (B-U)processes which is based on naturally formed NWs grown on a patterned substrate resulting from self-assembly of metallic clusters or exposition to a focused ion beam (FIB). The major goal consist to obtain organized and homogeneous NWs.

Les télécommunications laser en espace libre sont limitées en portée par la turbulence atmosphérique. L'optique adaptative, par la correction de la phase turbulente à l'émission du faisceau, a permis d'étendre leur domaine d'exploitation.Toutefois, sur de longues distances de propagation, cette correction n'est plus suffisante et il faut également précompenser l'amplitude du faisceau émis. De premières études numériques ont montré que le principe de retournement temporel, ou plus exactement de conjugaison de phase bidirectionnelle itérative, conduirait à des conditions satisfaisantes de focalisation du faisceau laser en fort régime de turbulence.Le principe de conjugaison de phase n'ayant été étudié que théoriquement jusque-là, mon travail de thèse s'est attaché à mettre en oeuvre un démonstrateur expérimental pour quantifier les performances de cette technique dans des conditions maîtrisées. En parallèle, une simulation de bout en bout de l'expérience a permis d'évaluer l'influence d'erreurs d'étalonnage sur les performances finales de la correction et de valider les résultats expérimentaux obtenus. Les points durs de la mise en oeuvre d'un système de télécommunications laser en espace libre ont ainsi été identifiés.L'ensemble de ces travaux constitue la toute première démonstration expérimentale du principe de retournement temporel optique. D'autres domaines d'application comme les lasers de puissance ou la propagation à travers des milieux biologiques très diffusants, nécessitant également de corriger le faisceau à l'émission, sont concernés. / Free Space Optical communications (FSO) are range limited due to atmospheric turbulence. Adaptive optics can mitigate turbulence effects by adding a phase modulation on the emitted laser beam. However, both phase and amplitude modulation are needed to perform long range FSO. Previous numerical studies have shown that iterative phase conjugation is an efficient modulation technique for lasercom systems.This PhD thesis is dedicated to the development and the realization of the first experimental demonstration of the iterative phase conjugation principle in a controlled turbulence environment. An optical bench, representative of a long range propagation through strong turbulence, has been scaled down to few-meters propagation in visible.Several methods for complex field measurement and modulation are numerically studied. Selected methods are implemented and tested, such as a novel focal plane technique for complex field measurement. Finally, iterative phase conjugation is performed and results cross-correlated with an end-to-end model representative of the optical bench.This work is the first experimental demonstration of the optical phase conjugation principle. Applications can be found in other fields than lasercoms, such as high power lasers or propagation through highly diffusing biological tissues, both in need of laser emission modulation.

PIXE, Particle Induced X-ray Emission, est une méthode d’analyse multiélémentaire, rapide, non destructive, basée sur la détection des rayons X caractéristiques émis suite à l’interaction de particules chargées avec la matière. Cette méthode est usuellement utilisée avec des protons accélérés à une énergie de l’ordre de quelques MeV dans des domaines d’applications variés, atteignant une limite de détection de l’ordre de quelques μg/g (ppm). Cependant, la profondeur d’analyse est relativement limitée. Grâce au cyclotron ARRONAX, nous pouvons utiliser des protons ou des particules alpha jusqu’à une énergie de 70 MeV pour mettre en œuvre la technique PIXE à haute énergie. Avec de telles énergies, nous excitons préférentiellement les raies X K, plus énergétiques que les raies L utilisées dans la PIXE classique pour l’analyse des éléments lourds. L’analyse d’échantillons épais, en profondeur, est ainsi accessible. Pour l’analyse des éléments légers, nous pouvons utiliser la détection de rayons gamma émis pas les noyaux excités en combinant les méthodes PIGE et PIXE. Nous allons tout d’abord présenter les caractéristiques et les principes d’analyse de la méthode PIXE à haute énergie que nous avons développée à ARRONAX. Nous détaillerons ensuite les performances atteintes, notamment en termes de limite de détection dans différentes conditions expérimentales. Enfin, nous présenterons les résultats obtenus pour l’analyse d’échantillons multicouches et la quantification d’éléments traces dans des échantillons épais. / Particle Induced X-ray Emission (PIXE) is a fast, nondestructive, multi-elemental analysis technique. It is based on the detection of characteristic X-rays due to the interaction of accelerated charged particles with matter. This method is successfully used in various application fields using low energy protons (energies around few MeV), reaching a limit of detection of the order the μg/g (ppm). At this low energy, the depth of analysis is limited. At the ARRONAX cyclotron, protons and alpha particles are delivered with energy up to 70 MeV, allowing the development of the High Energy PIXE technique. Thanks to these beams, we mainly produce KX-rays, more energetic than the LX-rays used with standard PIXE for the heavy elements analysis. Thus, in depth analysis in thick materials is achievable. For light element analysis, the PIGE technique, based on the detection of gamma rays emitted by excited nuclei, may be used in combination with PIXE. First of all, we will introduce the characteristics and principles of high energy PIXE analysis that we have developed at ARRONAX. Then we will detail the performance achieved, particularly in terms of detection limit in various experimental conditions. Finally, we present the results obtained for the analysis of multilayer samples and quantification of trace elements in thick samples.