An Adaptive Machine Learning Decision System for Flexible Predictive Maintenance

Process monitoring and Predictive Maintenance (PdM) are gaining increasing attention in most manufacturing environments as a means of reducing maintenance related costs and downtime. This is especially true in industries that are data intensive such as semiconductor manufacturing. In this paper an adaptive PdM based flexible maintenance scheduling decision support system, which pays particular attention to associated opportunity and risk costs, is presented. The proposed system, which employs Machine Learning and regularized regression methods, exploits new information as it becomes available from newly processed components to refine remaining useful life estimates and associated costs and risks. The system has been validated on a real industrial dataset related to an Ion Beam Etching process for semiconductor manufacturing.

Quantum control of photodissociation using intense, femtosecond pulses shaped with third order dispersion

We demonstrate the ability to control the molecular dissociation rate using femtosecond pulses shaped with third-order dispersion (TOD). Explicitly, a significant 50% enhancement in the dissociation yield for the low lying vibrational levels (v ∼ 6) of an H+2 ion-beam target was measured as a function of TOD. The underlying mechanism responsible for this enhanced dissociation was theoretically identified as non-adiabatic alignment induced by the pre-pulses situated on the leading edge of pulses shaped with negative TOD. This control scheme is expected to work in other molecules as it does not rely on specific characteristics of our test-case H+2 molecule.

Alpha particles induce pan-nuclear phosphorylation of H2AX in primary human lymphocytes mediated through ATM

The use of high linear energy transfer radiations in the form of carbon ions in heavy ion beam lines or alpha particles in new radionuclide treatments has increased substantially over the past decade and will continue to do so due to the favourable dose distributions they can offer versus conventional therapies. Previously it has been shown that exposure to heavy ions induces pan-nuclear phosphorylation of several DNA repair proteins such as H2AX and ATM in vitro. Here we describe similar effects of alpha particles on ex vivo irradiated primary human peripheral blood lymphocytes. Following alpha particle irradiation pan-nuclear phosphorylation of H2AX and ATM, but not DNA-PK and 53BP1, was observed throughout the nucleus. Inhibition of ATM, but not DNA-PK, resulted in the loss of pan-nuclear phosphorylation of H2AX in alpha particle irradiated lymphocytes. Pan-nuclear gamma-H2AX signal was rapidly lost over 24h at a much greater rate than foci loss. Surprisingly, pan-nuclear gamma-H2AX intensity was not dependent on the number of alpha particle induced double strand breaks, rather the number of alpha particles which had traversed the cell nucleus. This distinct fluence dependent damage signature of particle radiation is important in both the fields of radioprotection and clinical oncology in determining radionuclide biological dosimetry and may be indicative of patient response to new radionuclide cancer therapies.

Observation of the inhomogeneous spatial distribution of MeV ions accelerated by the hydrodynamic ambipolar expansion of clusters

An inhomogeneous spatial distribution of laser accelerated carbon/oxygen ions produced via the hydrodynamic ambipolar expansion of CO₂ clusters has been measured by using CR-39 detectors. An inhomogeneous etch pits spatial distribution has appeared on the etched CR-39 detector installed on the laser propagation direction, while homogeneous ones are appeared on those installed at 45°and 90°from the laser propagation direction. From the range of ions in CR-39 obtained by using the multi-step etching technique, the averaged energies of carbon/oxygen ions for all directions are determined as 0.78 ± 0.09 MeV/n. The number of ions in the laser propagation direction is about 1.5 times larger than those in other directions. The inhomogeneous etch pits spatial distribution in the laser propagation direction could originate from an ion beam collimation and modulation by the effect of electromagnetic structures created in the laser plasma.

Ultrafast charge dynamics initiated by high-intensity, ultrashort laser-matter interaction

The interaction of high‐intensity laser pulses with matter releases instantaneously ultra‐large currents of highly energetic electrons, leading to the generation of highly‐transient, large‐amplitude electric and magnetic fields. We report results of recent experiment in which such charge dynamics have been studied by using proton probing techniques able to provide maps of the electrostatic fields with high spatial and temporal resolution. The dynamics of ponderomotive channelling in underdense plasmas have been studied in this way, as also the processes of Debye sheath formation and MeV ion front expansion at the rear of laser‐irradiated thin metallic foils. An application employing laser‐driven impulsive fields for energy‐selective ion beam focusing is also presented.

 Graphitization of Small Carbonate Samples for Paleoceanographic Research at the Godwin Radiocarbon Laboratory, University of Cambridge

A new radiocarbon preparation facility was set up in 2010 at the Godwin Laboratory for Palaeoclimate Research, at the University of Cambridge. Samples are graphitized via hydrogen reduction on an iron powder catalyst before being sent to the Chrono Centre, Belfast, or the Australian National University for accelerator mass spectrometry (AMS) analysis. The experimental setup and procedure have recently been developed to investigate the potential for running small samples of foraminiferal carbonate. By analyzing background values of samples ranging from 0.04 to 0.6 mg C along with similar sized secondary standards, the setup and experimental procedures were optimized for small samples. “Background” modern 14C contamination has been minimized through careful selection of iron powder, and graphitization has been optimized through the use of “small volume” reactors, allowing samples containing as little as 0.08 mg C to be graphitized and accurately dated. Graphitization efficiency/fractionation is found not to be the main limitation on the analysis of samples smaller than 0.07 mg C, which rather depends primarily on AMS ion beam optics, suggesting further improvements in small sample analysis might yet be achieved with our methodology.

Valence-shell photoionization of chlorinelike Ar+ ions

Absolute cross-section measurements for valence-shell photoionization of Ar + ions are reported for photon energies ranging from 27.4 eV to 60.0 eV. The data, taken by merging beams of ions and synchrotron radiation at a photon energy resolution of 10 meV, indicate that the primary ion beam was a statistically weighted mixture of the 2P o3/2 ground state and the 2P o1/2 metastable state of Ar +. Photoionization of this Cell-like ion is characterized by multiple Rydberg series of autoionizing resonances superimposed on a direct photoionization continuum. Observed resonance lineshapes indicate interference between indirect and direct photoionization channels. Resonance features are spectroscopically assigned and their energies and quantum defects are tabulated. The measurements are satisfactorily reproduced by theoretical calculations based on an intermediate coupling semi-relativistic Breit-Pauli approximation.

Laser Driven Neutron Sources: Characteristics, Applications and Prospects

The basics of laser driven neutron sources, properties and possible applications are discussed. We describe the laser driven nuclear processes which trigger neutron generation, namely, nuclear reactions induced by laser driven ion beam (ion n), thermonuclear fusion by implosion and photo-induced nuclear (gamma n) reactions. Based on their main properties, i.e. point source (< 100 μm) and short durations (< ns), different applications are described, such as radiography, time-resolved spectroscopy and pump-probe experiments. Prospects on the development of laser technology suggest that, as higher intensities and higher repetition rate lasers become available (for example, using DPSSL technology), laser driven methodologies may provide neutron fluxes comparable to that achieved by accelerator driven neutron sources in the near future.

Caracterização e modificação de heteroestruturas de nitretos do grupo III

Os nitretos binários semicondutores do grupo III, e respetivos compostos, são vastamente estudados devido à sua possível aplicabilidade em dispositivos optoeletrónicos, tais como díodos emissores de luz (LEDs) e LASERs, assim como dispositivos para a eletrónica de elevadas temperatura, potência e frequência. Enquanto se concretizou a comercialização na última década de LEDs e LASERs recorrendo ao ternário In1-yGayN, estudos das propriedades fundamentais estruturais e óticas, assim como de técnicas de processamento no desenvolvimento de novas aplicações de outros ternários do grupo III-N encontram-se na sua fase inicial. Esta tese apresenta a investigação experimental de filmes finos epitaxiais de Al1-xInxN crescidos sobre camadas tampão de GaN e de Al1-yGayN e o estudo do recozimento e implantação de super-redes (SL) compostas por pontos quânticos de GaN (QD) envolvidos por camadas de AlN. Apesar do hiato energético do Al1-xInxN poder variar entre os 0,7 eV e os 6,2 eV e, por isso, numa gama, consideravelmente superior à dos ternários Al1-yGayN e InyGa1-yN, o primeiro é o menos estudado devido a dificuldades no crescimento de filmes com elevada qualidade cristalina. É efetuada, nesta tese, uma caracterização estrutural e composicional de filmes finos de Al1-xInxN crescidos sobre camadas tampão de GaN e de Al1-yGayN usando técnicas de raios-X, feixe de iões e de microscopia. Mostra-se que o Al1-xInxN pode ser crescido com elevada qualidade cristalina quando a epitaxia do crescimento se aproxima da condição de rede combinada do Al1-xInxN e da camada tampão (GaN ou Al1-yGayN), isto é, com conteúdo de InN de ~18%, quando crescido sobre uma camada de GaN. Quando o conteúdo de InN é inferior/superior à condição de rede combinada, fenómenos de relaxação de tensão e deterioração do cristal tais como o aumento da rugosidade de superfície prejudicam a qualidade cristalina do filme de Al1-xInxN. Observou-se que a qualidade dos filmes de Al1-xInxN depende fortemente da qualidade cristalina da camada tampão e, em particular, da sua morfologia e densidade de deslocações. Verificou-se que, dentro da exatidão experimental, os parâmetros de rede do ternário seguem a lei empírica de Vegard, ou seja, variam linearmente com o conteúdo de InN. Contudo, em algumas amostras, a composição determinada via espetrometria de retrodispersão de Rutherford e difração e raios-X mostra valores discrepantes. Esta discrepância pode ser atribuída a defeitos ou impurezas capazes de alterar os parâmetros de rede do ternário. No que diz respeito às SL dos QD e camadas de AlN, estudos de recozimento mostraram elevada estabilidade térmica dos QD de GaN quando estes se encontram inseridos numa matriz de AlN. Por implantação iónica, incorporou-se európio nestas estruturas e, promoveu-se a ativação ótica dos iões de Eu3+ através de tratamentos térmicos. Foram investigados os efeitos da intermistura e da relaxação da tensão ocorridos durante o recozimento e implantação nas propriedades estruturais e óticas. Verificou-se que para fluências elevadas os defeitos gerados por implantação são de difícil remoção. Contudo, a implantação com baixa fluência de Eu, seguida de tratamento térmico, promove uma elevada eficiência e estabilidade térmica da emissão vermelha do ião lantanídeo incorporado nos QD de GaN. Estes resultados são, particularmente relevantes, pois, na região espetral indicada, a eficiência quântica dos LEDs convencionais de InGaN é baixa.

An experimental investigation of dechanneling in copper single crystals

This investigation comprises a comparison of experimental and theoretical dechanneling of MeV protons in copper single crystals. Dechanneling results when an ion's transverse energy increases to the value where the ion can undergo small impact parameter collisions with individual atoms. Depth dependent dechanneling rates were determined as functions of lattice temperature, ion beam energy and crystal axis orientation. Ion beam energies were IMeV and 2MeV,temperatures ranged from 35 K to 280 K and the experiment was carried out along both the (lOa) and <110) axes. Experimental data took the form of aligned and random Rutherford backscattered energy spectra. Dechanneling rates were extracted from these spectra using a single scattering theory that took explicit account of the different stopping powers experienced by channeled and dechanneled ions and also included a correction factor to take into account multiple scattering effects along the ion's trajectory. The assumption of statistical equilibrium and small angle scattering of the channeled ions allows a description of dechanneling in terms of the solution of a diffusion like equation which contains a so called diffusion function. The diffusion function is shown to be related to the increase in average transverse energy. Theoretical treatments of increase in average transverse energy due to collisions of projectiles with channel electrons and thermal perturbations in the lattice potential are reviewed. Using the diffusion equation and the electron density in the channel centre as a fitting parameter dechanneling rates are extracted. Excellent agreement between theory and experiment has been demonstrated. Electron densities determined in the fitting procedure appear to be realistic. The surface parameters show themselves to be good indicators of the quality of the crystal.

L’étude de l’InP et du GaP suite à l’implantation ionique de Mn et à un recuit thermique

Cette thèse est dédiée à l’étude des matériaux InMnP et GaMnP fabriqués par implantation ionique et recuit thermique. Plus précisément nous avons investigué la possibilité de former par implantation ionique des matériaux homogènes (alliages) de InMnP et GaMnP contenant de 1 à 5 % atomiques de Mn qui seraient en état ferromagnétique, pour des possibles applications dans la spintronique. Dans un premier chapitre introductif nous donnons les motivations de cette recherche et faisons une revue de la littérature sur ce sujet. Le deuxième chapitre décrit les principes de l’implantation ionique, qui est la technique utilisée pour la fabrication des échantillons. Les effets de l’énergie, fluence et direction du faisceau ionique sur le profil d’implantation et la formation des dommages seront mis en évidence. Aussi dans ce chapitre nous allons trouver des informations sur les substrats utilisés pour l’implantation. Les techniques expérimentales utilisées pour la caractérisation structurale, chimique et magnétique des échantillons, ainsi que leurs limitations sont présentées dans le troisième chapitre. Quelques principes théoriques du magnétisme nécessaires pour la compréhension des mesures magnétiques se retrouvent dans le chapitre 4. Le cinquième chapitre est dédié à l’étude de la morphologie et des propriétés magnétiques des substrats utilisés pour implantation et le sixième chapitre, à l’étude des échantillons implantés au Mn sans avoir subi un recuit thermique. Notamment nous allons voir dans ce chapitre que l’implantation de Mn à plus que 1016 ions/cm2 amorphise la partie implantée du matériau et le Mn implanté se dispose en profondeur sur un profil gaussien. De point de vue magnétique les atomes implantés se trouvent dans un état paramagnétique entre 5 et 300 K ayant le spin 5/2. Dans le chapitre 7 nous présentons les propriétés des échantillons recuits à basses températures. Nous allons voir que dans ces échantillons la couche implantée est polycristalline et les atomes de Mn sont toujours dans un état paramagnétique. Dans les chapitres 8 et 9, qui sont les plus volumineux, nous présentons les résultats des mesures sur les échantillons recuits à hautes températures : il s’agit d’InP et du GaP implantés au Mn, dans le chapitre 8 et d’InP co-implanté au Mn et au P, dans le chapitre 9. D’abord, dans le chapitre 8 nous allons voir que le recuit à hautes températures mène à une recristallisation épitaxiale du InMnP et du GaMnP; aussi la majorité des atomes de Mn se déplacent vers la surface à cause d’un effet de ségrégation. Dans les régions de la surface, concentrés en Mn, les mesures XRD et TEM identifient la formation de MnP et d’In cristallin. Les mesures magnétiques identifient aussi la présence de MnP ferromagnétique. De plus dans ces mesures on trouve qu’environ 60 % du Mn implanté est en état paramagnétique avec la valeur du spin réduite par rapport à celle trouvée dans les échantillons non-recuits. Dans les échantillons InP co-implantés au Mn et au P la recristallisation est seulement partielle mais l’effet de ségrégation du Mn à la surface est beaucoup réduit. Dans ce cas plus que 50 % du Mn forme des particules MnP et le restant est en état paramagnétique au spin 5/2, dilué dans la matrice de l’InP. Finalement dans le dernier chapitre, 10, nous présentons les conclusions principales auxquels nous sommes arrivés et discutons les résultats et leurs implications.

Mécanismes de déformation de nanoparticules d’Au par irradiation ionique

Résumé Dans la présente thèse, nous avons étudié la déformation anisotrope par bombardement ionique de nanoparticules d'or intégrées dans une matrice de silice amorphe ou d'arséniure d’aluminium cristallin. On s’est intéressé à la compréhension du mécanisme responsable de cette déformation pour lever toute ambigüité quant à l’explication de ce phénomène et pour avoir une interprétation consistante et unique. Un procédé hybride combinant la pulvérisation et le dépôt chimique en phase vapeur assisté par plasma a été utilisé pour la fabrication de couches nanocomposites Au/SiO2 sur des substrats de silice fondue. Des structures à couches simples et multiples ont été obtenues. Le chauffage pendant ou après le dépôt active l’agglomération des atomes d’Au et par conséquent favorise la croissance des nanoparticules. Les nanocomposites Au/AlAs ont été obtenus par implantation ionique de couches d’AlAs suivie de recuit thermique rapide. Les échantillons des deux nanocomposites refroidis avec de l’azote liquide ont été irradiés avec des faisceaux de Cu, de Si, d’Au ou d’In d’énergie allant de 2 à 40 MeV, aux fluences s'étendant de 1×1013 à 4×1015 ions/cm2, en utilisant le Tandem ou le Tandetron. Les propriétés structurales et morphologiques du nanocomposite Au/SiO2 sont extraites en utilisant des techniques optiques car la fréquence et la largeur de la résonance plasmon de surface dépendent de la forme et de la taille des nanoparticules, de leur concentration et de la distance qui les séparent ainsi que des propriétés diélectriques du matériau dans lequel les particules sont intégrées. La cristallinité de l’arséniure d’aluminium est étudiée par deux techniques: spectroscopie Raman et spectrométrie de rétrodiffusion Rutherford en mode canalisation (RBS/canalisation). La quantité d’Au dans les couches nanocomposites est déduite des résultats RBS. La distribution de taille et l’étude de la transformation de forme des nanoparticules métalliques dans les deux nanocomposites sont déterminées par microscopie électronique en transmission. Les résultats obtenus dans le cadre de ce travail ont fait l’objet de trois articles de revue. La première publication montre la possibilité de manipuler la position spectrale et la largeur de la bande d’absorption des nanoparticules d’or dans les nanocomposites Au/SiO2 en modifiant leur structure (forme, taille et distance entre particules). Les nanoparticules d’Au obtenues sont presque sphériques. La bande d’absorption plasmon de surface (PS) correspondante aux particules distantes est située à 520 nm. Lorsque la distance entre les particules est réduite, l’interaction dipolaire augmente ce qui élargit la bande de PS et la déplace vers le rouge (602 nm). Après irradiation ionique, les nanoparticules sphériques se transforment en ellipsoïdes alignés suivant la direction du faisceau. La bande d’absorption se divise en deux bandes : transversale et longitudinale. La bande correspondante au petit axe (transversale) est décalée vers le bleu et celle correspondante au grand axe (longitudinale) est décalée vers le rouge indiquant l’élongation des particules d’Au dans la direction du faisceau. Le deuxième article est consacré au rôle crucial de la déformation plastique de la matrice et à l’importance de la mobilité des atomes métalliques dans la déformation anisotrope des nanoparticules d’Au dans les nanocomposites Au/SiO2. Nos mesures montrent qu'une valeur seuil de 2 keV/nm (dans le pouvoir d'arrêt électronique) est nécessaire pour la déformation des nanoparticules d'or. Cette valeur est proche de celle requise pour la déformation de la silice. La mobilité des atomes d’Au lors du passage d’ions est confirmée par le calcul de la température dans les traces ioniques. Le troisième papier traite la tentative de formation et de déformation des nanoparticules d’Au dans une matrice d’arséniure d’aluminium cristallin connue pour sa haute résistance à l’amorphisation et à la déformation sous bombardement ionique. Le résultat principal de ce dernier article confirme le rôle essentiel de la matrice. Il s'avère que la déformation anisotrope du matériau environnant est indispensable pour la déformation des nanoparticules d’or. Les résultats expérimentaux mentionnés ci-haut et les calculs de températures dans les traces ioniques nous ont permis de proposer le scénario de déformation anisotrope des nanoparticules d’Au dans le nanocomposite Au/SiO2 suivant: - Chaque ion traversant la silice fait fondre brièvement un cylindre étroit autour de sa trajectoire formant ainsi une trace latente. Ceci a été confirmé par la valeur seuil du pouvoir d’arrêt électronique. - L’effet cumulatif des impacts de plusieurs ions conduit à la croissance anisotrope de la silice qui se contracte dans la direction du faisceau et s’allonge dans la direction perpendiculaire. Le modèle de chevauchement des traces ioniques (overlap en anglais) a été utilisé pour valider ce phénomène. - La déformation de la silice génère des contraintes qui agissent sur les nanoparticules dans les plans perpendiculaires à la trajectoire de l’ion. Afin d’accommoder ces contraintes les nanoparticules d’Au se déforment dans la direction du faisceau. - La déformation de l’or se produit lorsqu’il est traversé par un ion induisant la fusion d’un cylindre autour de sa trajectoire. La mobilité des atomes d’or a été confirmée par le calcul de la température équivalente à l’énergie déposée dans le matériau par les ions incidents. Le scénario ci-haut est compatible avec nos données expérimentales obtenues dans le cas du nanocomposite Au/SiO2. Il est appuyé par le fait que les nanoparticules d’Au ne se déforment pas lorsqu’elles sont intégrées dans l’AlAs résistant à la déformation.

Préparation d'échantillons pour l'étude par GISAXS des mécanismes de déformation des matériaux par faisceaux d'ions lourds de haute énergie.

Le mécanisme menant à des déformations structurales suivant le bombardement d'échantillons de a-Si d'un faisceau d'ions lourds et rapides est sujet de controverses. Nous nous sommes penchés sur l'hypothèse de la formation d'une zone liquide causée par la déposition d'énergie des ions incidents dans le contexte de la théorie du pic thermique. Des échantillons de silicium amorphe furent préparés dans le but d'observer les indices d'une transition de phase l-Si/a-Si suivant la déposition locale d'énergie sur le parcours d'un ion lourd énergétique dans le a-Si. Les échantillons furent implantés d'impuretés de Cu ou d'Ag avant d'être exposés à un faisceau d'ions Ag12+ de 70 MeV. L'utilisation de l'analyse GISAXS est projetée afin d'observer une concentration locale d'impuretés suivant leur ségrégation sur la trace de l'ion. Des masques d'implantation nanométriques d'oxide d'aluminium ont été fabriqués afin d'augmenter la sensibilité de l'analyse GISAXS et une méthode d'alignement de ces masques selon la direction du faisceau fut développée. Le bombardement d'échantillons au travers de ces masques a donné lieu à un réseau de sites d'impacts isolés presque équidistants.

É́tude sur la cinétique des défauts structuraux dans le silicium amorphe.

Cette thèse présente à la fois des résultats de simulations numériques en plus de ré- sultats expérimentaux obtenus en laboratoire sur le rôle joué par les défauts de structure dans le silicium amorphe. Nos travaux de simulation numérique furent réalisés avec une nouvelle méthode de simulation Monte-Carlo cinétique pour décrire l’évolution tempo- relle de modèles de silicium amorphe endommagés sur plusieurs échelles de temps jus- qu’à une seconde à la température pièce. Ces simulations montrent que les lacunes dans le silicium amorphe sont instables et ne diffusent pas sans être détruites. Nous montrons également que l’évolution d’un modèle de silicium amorphe endommagé par une colli- sion ionique lors d’un recuit peut être divisée en deux phases : la première est dominée exclusivement par la diffusion et la création/destruction de défauts de liaison, alors que la deuxième voit les créations/destructions de liens remplacées par des échanges de liens entre atomes parfaitement coordonnés. Les défauts ont aussi un effet sur la viscosité du silicium amorphe. Afin d’approfondir cette question, nous avons mesuré la viscosité du silicium amorphe et du silicium amorphe hydrogéné sous l’effet d’un faisceau d’ions. Nous montrons que la variation de la viscosité dans les deux matériaux est différente : le silicium amorphe hydrogéné a une viscosité constante en fonction de la fluence des ions alors que le silicium amorphe pur a une viscosité qui augmente de façon linéaire. Pour de faibles fluences, la viscosité du silicium hydrogéné est plus grande que la viscosité sans hydrogène. La présence d’hydrogène diminue également l’amplitude de la variation logarithmique de la contrainte observée lors de la relaxation à la température de la pièce.

Spectrum simulation of rough and nanostructured targets from their 2D and 3D image by Monte Carlo methods

Corteo is a program that implements Monte Carlo (MC) method to simulate ion beam analysis (IBA) spectra of several techniques by following the ions trajectory until a sufficiently large fraction of them reach the detector to generate a spectrum. Hence, it fully accounts for effects such as multiple scattering (MS). Here, a version of Corteo is presented where the target can be a 2D or 3D image. This image can be derived from micrographs where the different compounds are identified, therefore bringing extra information into the solution of an IBA spectrum, and potentially significantly constraining the solution. The image intrinsically includes many details such as the actual surface or interfacial roughness, or actual nanostructures shape and distribution. This can for example lead to the unambiguous identification of structures stoichiometry in a layer, or at least to better constraints on their composition. Because MC computes in details the trajectory of the ions, it simulates accurately many of its aspects such as ions coming back into the target after leaving it (re-entry), as well as going through a variety of nanostructures shapes and orientations. We show how, for example, as the ions angle of incidence becomes shallower than the inclination distribution of a rough surface, this process tends to make the effective roughness smaller in a comparable 1D simulation (i.e. narrower thickness distribution in a comparable slab simulation). Also, in ordered nanostructures, target re-entry can lead to replications of a peak in a spectrum. In addition, bitmap description of the target can be used to simulate depth profiles such as those resulting from ion implantation, diffusion, and intermixing. Other improvements to Corteo include the possibility to interpolate the cross-section in angle-energy tables, and the generation of energy-depth maps.