Ion Impact Energy Research Articles

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy (E in) and angles (A in) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in, and W impurity transport.

Two-dimensional collisional particle model of the divertor sheath with electron emissive walls

A novel two-dimensional particle-in-cell (PIC) code, named Divertor Edge Simulator of Plasma-wall Interaction with Consistent COllisions (DESPICCO) and developed at CNR-ISTP, is capable of simulating the thin plasma layer of several millimeters, adjacent to the divertor tiles of a Tokamak fusion reactor. Here, kinetic effects and non-neutral plasma physics in the Debye sheath can be self-consistently captured by the PIC approach. The code is firstly benchmarked against literature one-dimensional codes and additional theoretical predictions for a magnetized sheath. Then, it is applied to a realistic divertor scenario featuring an attached plasma with monoblocks (MBs) radial misalignment and gaps, to compute the energy flux amplification factor at the exposed MB edge. A non-ambipolar local current density close to the leading edge and an average sheath heat transmission coefficient larger than the one predicted by classical sheath theory, are found. The effects of electron wall emission and plasma-gas collisions on the ion Mach number and on particle and energy fluxes to the walls are finally estimated to determine future guidelines for simulations. Ion collisions with recycled neutrals and both secondary and thermionic electron emission from the wall are found to have a relevant impact, with the overall effect of reducing by 25% the average ion impact energy, and by 15%–20% the total heavy particles energy flux to the walls, with relevant implications on the divertor wall erosion.

Establishing a Best Practice for SDTrimSP Simulations of Solar Wind Ion Sputtering

Solar wind (SW) ion irradiation on airless bodies can play an important role in altering their surface properties and surrounding exosphere. Much of the ion sputtering data needed for exosphere studies come from binary collision approximation sputtering models such as TRansport of Ions in Matter and its more recent extension, SDTrimSP. These models predict the yield and energy distribution of sputtered atoms, along with the depth of deposition and damage of the substrate, all as a function of the incoming ion type, impact energy, and impact angle. Within SDTrimSP there are several user-specific inputs that have been applied differently in previous SW ion sputtering simulations. These parameters can influence the simulated behavior of both the target and sputtered atoms. Here, we have conducted a sensitivity study into the SDTrimSP parameters in order to determine a best practice for simulating SW ion impacts onto planetary surfaces. We demonstrate that ion sputtering behavior is highly sensitive to several important input parameters including the ion impact angle and energy distribution and the ejected atom surface binding energy. Furthermore, different parameters can still result in similarities in the total sputtering yield, potentially masking large differences in other sputtering-induced behaviors such as the elemental yield, surface concentration, and damage production. Therefore, it is important to consider more than just the overall sputtering behavior when quantifying the importance of different parameters. This study serves to establish a more consistent methodology for simulations of SW-induced ion sputtering on bodies such as Mercury and the Moon, allowing for more accurate comparisons between studies.

Effect of edge ion temperature on the divertor tungsten sputtering in WEST

The influence of upstream ion temperature in the scrape-off layer (SOL) on the tungsten (W) sputtering in the divertor is studied in the WEST tokamak. For an almost constant power into the SOL, the upstream ion temperature and its ratio over the electron temperature gradually increase with the decrease of electron density in the SOL. This increment is observed to enhance the energy transfer from ions to electrons. This increases the downstream electron temperature and by coupling of electrons and ions, the impact energy of ions causing W sputtering in the divertor. This enhancement mechanism may become crucial to sputtering the W material for high upstream T i/T e ratio since the impact energy of ions in the divertor would increase accordingly.

Influence of surface temperature, ion impact energy, and bulk tungsten content on the sputtering of steels: In situ observations from plasma exposure in PSI-2

• PSI-2 deuterium plasma exposures of 3 RAFM steels. • In-situ measurement of relative sputtering yields by optical emission spectroscopy. • In-situ measurement of the evolution of surface morphology by IR camera. • Post-mortem measurement of tungsten enrichment by RBS, SEM and EDX. • Independent measurements allow to disentangle the impact on the erosion of Fe. This study presents results from PSI-2 deuterium plasma exposures of reduced-activation ferritic martensitic steels with different tungsten content at a peak plasma flux of 9 · 10 20 D/m 2 s and 420–920 K sample temperature. Via optical emission spectroscopy the relative sputtering yield of Fe and Cr, and via an infrared camera the evolution of surface morphology is studied in-situ in 2D. After the exposure, RBS, SEM and EDX reveal tungsten enrichment and the formation of surface structures. Due to the independent measurement of surface morphology and tungsten enrichment, it is possible to disentangle the impact of both effects on the erosion of Fe. As part of this study, we investigate the effects of tungsten content, sample temperature and ion impact energy. The results show a reduction of the sputtering yield, depending on the exposure conditions, up to a factor of 2 at 130 eV ion impact energy and 720 K sample temperature.

Investigation of Vacuum Arc-Deposited ta-C and ta-C:N Thin Films on Silicon and Stainless-Steel Foil Substrates Using Raman Spectroscopy

In this paper, we investigated the vacuum arc-deposited tetrahedral amorphous carbon (ta-C) and nitrogen-doped ta-C (ta-C:N) thin films grown on silicon and stainless-steel foil substrates by the visible Raman Spectroscopy. We found that carbon films grown on silicon surfaces prefer more sp3 hybridized carbon compared to the stainless-steel foil, which prefer more sp2 hybridized carbon even though the films are grown concurrently under the same conditions. The impact energy of plasma ions modifies the interfacial layer formation, giving a strong dependence of the sp2/sp3 ratio with the substrate bias, behaving differently for different substrates. However, nitrogen doping in amorphous carbon thin films helps to increase sp2 contents uniformly on both the surfaces for a fixed substrate bias. Understanding such interfaces are of interest for many electronic, optoelectronic, and energy storage devices as these interfaces get buried and can influence the properties significantly.

Data-driven surrogate modeling of hPIC ion energy-angle distributions for high-dimensional sensitivity analysis of plasma parameters' uncertainty

We present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.1

Ion beam figuring with focused anode layer thruster.

This work presents the peculiarities of cone ion beam formation with a focused thruster with anode layer (TAL) and its application to silicon carbide (SiC) ion beam figuring. Modeling results of Lorentz E × B force distribution in the discharge gap are presented. 3D particle tracing for keV Ar ions is carried out for the first time in the beam drift region of TAL with magnetic lens. Extracted ion beam full width at half maxima is about 2mm in the focal plane, where the SiC etching rate reaches 0.5 µm/min. The SiC sputter yields are measured as a function of the Ar ion impact energy and beam incidence angle. The maximum sputter yield of 2.8 atom/ion is observed at 45° of the beam-sample angle for the Si targets. Furthermore, the maximum sputter yield value of 1.7 atom/ion is measured at 30° of the beam-sample angle for the SiC targets. The novelty of present research is in the application of focused TAL keV Ar ion beam to the SiC ion beam figuring.

Inducing thermionic emission from lanthanum hexaboride probes in Magnum-PSI

For thermionic emission rates exceeding the incident plasma electron flux, recent theory proposes an inverse sheath regime, with promising properties for future application in fusion edge plasmas. With the aim of inducing thermionic emission in fusion-relevant plasma conditions, several lanthanum hexaboride probes were heated in the linear plasma generator Magnum-PSI. During exposures at low plasma power and additional pulsed laser heating, the probe’s floating potential was reduced by up to 12%, providing a possible indication of thermionic emission. However, these observations coincided with rapid erosion of probe material, attributed to enhanced lanthanum self-sputtering. During follow-up experiments with helium plasmas at electron temperatures around 1 eV, the lanthanum ion impact energy and sputtering yield were reduced, and rapid erosion was avoided, thus confirming the thesis of self-sputtering. A parameter scan of plasma power resulted in LaB 6 surface temperatures up to 2450 °C, exceeding the theoretical inverse sheath threshold temperature by over 300 °C. However, the probe’s floating potential did not deviate from reference measurements using a probe with high electronic work function, indicating absence of strong thermionic emission. This apparent discrepancy is attributed to the effects of probe surface modifications as observed during these experiments: impurity deposition, erosion and cavity formation. These modifications possibly affected the LaB 6 electronic work function, thereby keeping the inverse sheath threshold out of reach. In conclusion, although LaB 6 has one of the lowest work functions available, the inverse sheath threshold conditions could not be reached with the present setup in Magnum-PSI. Surface modifications thus do form a limiting factor for the application of LaB 6 in fusion-relevant plasma conditions. Moreover, the window of stable operation for LaB 6 in dense hydrogen plasmas is limited below ∼ 1.5 eV, and does not overlap with the conditions expected in the edge region of future fusion devices like ITER. • LaB 6 probes were heated up to 2450 °C under fusion-relevant plasma conditions. • Surface modifications limit the effective operation of powder-pressed LaB 6 . • Lanthanum self-sputtering in dense plasma conditions leads to rapid erosion. • The LaB 6 operation window does not overlap with ITER divertor conditions.

Low-Energy Plasma Source for Clean Vacuum Environments: EUV Lithography and Optical Mirrors Cleaning

Plasma cleaning of extreme ultra-violet (EUV) optics for the semiconductor industry requires atomic-level precision. Low-energy ions and neutrals can be highly beneficial for this purpose. However, ion energies in many industrial capacitively or inductively coupled plasmas may be too high for atomic precision processing so that ions can cause sputtering and re-deposition of materials and produce vacuum system contamination. We discuss two sources that create low-energy ions: a capacitively coupled (CCP) plasma source operating at several tens of MHz and an electron beam gas ionization source. The goal here is to minimize the contamination by limiting the ion impact energy to a few tens of electron-volts. Ion energy and flux measurements on a grounded surface are characterized with a compact retarding field ion spectrometer. Plasma-induced contamination is quantified using X-ray photo-electron spectroscopy (XPS). Low ion energy plasma sources introducing little surface contamination may be interesting for cleaning and accelerated testing in EUV lithography (EUVL)-related research and for cleaning of front-end optical mirrors in fusion reactor diagnostics.

The influence of D2 pressure on D retention and release from Be co-deposits

Beryllium co-deposit thermal desorption studies are reported that focus on the effect of D 2 pressure during layer formation, as pressure has been linked in prior work to the formation of a high-retention-capacity sharp-release-feature. A pressure range from 0.3 to 13.3 Pa is explored. Additional insight into what drives the formation of T 0 is established through co-deposit formation parameter variations including deposition temperature spanning 373 K to 608 K, and ion impact energy in the range 10–100 eV. It is found that the formation of the sharp release feature exhibits a multi-parameter dependence, but is predominately favoured by low deposition temperature, higher D 2 pressure, and to a small degree, increased ion impact energy, during deposition. TESSIM simulations of the D 2 thermal release at the ITER bake temperatures of T b 1 (513 K) and T b 2 (623 K), show that the best efficacy for D removal (defined as the ratio of post to pre bake inventory) is found when the sharp release feature is present. However, when this feature is present, the total retained D inventory tends to be significantly increased. • Be co-deposit desorption studies are reported that focus on the role of D 2 pressure during layer formation. • A sharp release mechanism, T 0 , is found to be predominately inked to pressure. • TESSIM simulations show best ITER bake efficacy when a T 0 fraction is present in the co-deposit.

Electrical insulation characteristics of metal-insulator-metal structures using boron nitride dielectric films deposited with low-energy ion impact

Sp2-bonded boron nitride (BN) films are deposited on Ni and Si substrates with mean ion impact energies between a few eV and around 100 eV by surface-wave plasma enhanced chemical vapor deposition and their electrical insulation characteristics are examined for temperatures up to 300 °C. The overall crystallinity and order of sp2 structure in the films increase with decreasing ion energy, while the mass density becomes as high as sintered hexagonal BN only when the ion energy is above ∼50 eV. For a given temperature, the electrical resistivity of the films measured for metal-BN-metal structures increases with decreasing ion energy, and hence increasing crystallinity and order of sp2 structure, almost independent of mass density. A decrease in resistivity with temperature is small below 150 °C, then becomes large above 150 °C with a drastic increase in apparent activation energy for carriers, which is attributed to generation and transport of the carriers that get over higher barrier heights associated with stronger ion impact.

Mechanisms of crack generation in high-pure tungsten exposed to high power density plasma

The mechanisms of crack formation in the tungsten divertor plates irradiated by hydrogen plasma under ITER-type steady-state and transient events are investigated. The stresses and defects of samples‘ initial structure are annealed as a result of steady-state irradiation with Falcon ion source, which generates 2 keV hydrogen ion flux of ~ 1022 m-2s−1 and power density ~ 1.7 MW × m−2. Quasi-stationary plasma accelerator QSPA Kh-50 provides pulsed hydrogen plasma loads which operating mode simulates transient events (type I ELMs) in ITER: ion impact energy is about 400 eV, maximum plasma pressure is 0.32 MPa, pulse duration is 0.25 ms. The diffraction maxima profiles, positions, intensities, half-widths, and intensity distributions are analyzed. The origin of crack formation is found to be in the plastic deformation of surface layers by the twinning mechanism. The latter is reinforced by the interaction of twins with the hydrogen-filled micro-pores. The micro-pores appear as a result of hydrogen saturation of subtraction loops formed due to the coalescence of vacancy complexes.

Tungsten and Molybdenum Surfaces Exposed to Warm Deuterium Ion Plasma with Q-Nonextensive Distribution of Electrons

This work is concerned with the effect of non-Maxwellian electrons and obliqueness of magnetic field on magnetized plasma sheath characteristics, in which plasma interacts with tungsten (W) and molybdenum (Mo) surfaces via non-neutral plasma sheath using two fluids model. It is assumed that the singly charged positive ions are treated as warm fluid whereas the electrons obey q-nonextensive distribution. It is found that the q-nonextensive distributed electrons and the temperature of ions affect the entrance velocity of positive ions, which is a key parameter in the plasma sheath formation. Also, the nonextensive parameter q affects the distribution of ions and electrons in the sheath region and their distributions explicitly related with the electrostatic potential variation. The parallel and perpendicular components of ions velocity are affected by the obliqueness of magnetic field. As the nonextensivity of electrons increases, the gradient in electric potential increases towards the wall and hence the impact energy also increases. The obliqueness of magnetic field and impact energy of ions is a key factor that determines the physical sputtering rate, particle reflection and absorption from the target surface. Furthermore, the probability of particle reflection coefficient from the W-surface is higher than that of Mo-surface.

Effect of the boron-to-nitrogen ratio on leakage current characteristics of boron nitride films prepared by surface-wave plasma enhanced chemical vapor deposition

Sp2-bonded boron nitride (BN) films with different B-to-N ratios are deposited on Si and Ni substrates at an ion impact energy of around 50 eV by surface-wave plasma enhanced chemical vapor deposition. The overall crystallinity and order of sp2 bonding structure are almost retained with increasing BF3/N2 gas flow ratio in the plasma, while the B-to-N ratio in the film is varied from 0.89 to 0.96. The electrical resistivity of the films measured at room temperature for Ni-BN-Ni structures shows a large increase up to the order of 1010 Ω cm when the B-to-N ratio in the film is increased closer to the stoichiometry due to a lower density of defect states related to the boron vacancy. The characteristics of leakage current vs. applied bias voltage show Ohmic behavior at low fields and non-Ohmic behavior governed by the thermionic emission mechanism at high fields for any B-to-N ratio.

Possible mitigation of tokamak plasma–surface interactions using thermionic divertor plates with inverse sheaths

Recent studies demonstrate that when a plasma-facing surface emits a sufficient flux of electrons, it will form an inverse sheath. Here, we consider a possibility of using thermionic target plates with inverse sheaths as an innovative divertor operating scenario. We derive an electron heat flux boundary condition for inverse sheaths and show that for given power exhaust into a tokamak scrape-off-layer, an inverse sheath leads to a much lower target plasma electron temperature than a conventional sheath. Low enough target plasma temperatures for radiative divertor detachment could therefore be achieved using inverse sheaths instead of the usual need to inject neutral atoms that compromise the core plasma. Other advantages of inverse sheath detachment over conventional sheath operating scenarios include (a) ion impact energies are as low as possible, minimizing sputtering and tritium implantation, (b) surface recombination heat flux is reduced due to ion flux reduction, and (c) arcs are inhibited due to the sign of the surface electric field. This paper outlines the basic properties of inverse sheath detachment and considers the feasibility of implementation. We offer recommendations for future modeling efforts needed to better understand the effects of thermionic emission in tokamaks and whether inverse sheaths present a viable divertor solution.

Light-reflection-induced changes in the line shape of sputtered atoms

The erosion from plasma-facing components has to be monitored in many kind of laboratory and fusion plasmas. For this purpose, spectroscopy is an essential tool. Under certain conditions the particle flux can be calculated from the absolute line intensities of the sputtered material using so-called S/XB values. The impact of light reflection on the emission induced by sputtered particles at the mirror-grade polished surface of tungsten (W) and aluminum (Al) was investigated in a low-density (ne ≈ 2 × 1012 cm−3) and low-temperature (Te ≈ ) argon (Ar) plasma in the linear plasma device PSI-2 using high-resolution spectroscopy. Using the line shape affected by Doppler shift we show that the light reflection has a considerable impact on the number of measured photons and has to be taken into account for calculating particle fluxes. The Al target was sputtered by Ar ions at the incident ion energy of . The measured profile of the Al I line () was compared with a Doppler-shifted emission model based on the Thompson energy distribution function. In this new model, the instrumental broadening and the impact of the Zeeman effect were also taken into account. The parameter for the high-energy fall-off n of the energy distribution function (), the surface binding energy Eb and the surface reflectance were derived by comparing the experimental and the synthetic spectrum. The W target was sputtered by Ar ions at incident ion energies in the range of –. The influence of the ion impact energy on the energy distribution of the sputtered particles was demonstrated.

Impact of divertor material on neutral recycling and discharge fueling in DIII-D

Experiments with the lower divertor of DIII-D during the Metal Rings Campaign (MRC) show that the fraction F of atomic D in the total recycling flux is material-dependent and varies through the ELM cycle, which may affect divertor fueling. Between ELMs, FC ∼ 10% and FW ∼ 40%, consistent with expectations if all atomic recycling is due to reflections. During ELMs, FC increases to 50% and FW to 60%. In contrast, the total D recycling coefficient including atoms and molecules R stays close to unity near the strike point where the surface is saturated with D. During ELMs, R can deviate from unity, increasing during high energy ELM-ion deposition (net D release) and decreasing at the end of the ELM which leads to ability of the target to trap the ELM-deposited D. The increase of R > 1 in response to an increase in ion impact energy Ei has been studied with small divertor target samples using Divertor Materials Evaluation System (DiMES). An electrostatic bias was applied to DiMES to change Ei by 90 eV. On all studied materials including C, Mo, uncoated and W-coated TZM (>99% Mo, Ti, and Zr alloy), W, and W fuzz, an increase of Ei transiently increased the D yield (and R) by ∼10%. On C there was also an increase in the molecular D2 yield, probably due to ion-induced D2 desorption. Despite the measured increase in F on W compared to C, attached H-mode shots with OSP on W during MRC did not demonstrate a higher pedestal density. About 8% increase in the edge density could be seen only in attached L-mode scenarios. The difference can be explained by higher D trapping in the divertor and lower divertor fueling efficiency in H- versus L-mode.

Nonlinear dependence of plasma potential and ion impinging energy on energy loss of edge localized modes

Particle-in-cell (PIC) modelling has been performed to investigate the impact of energy loss during edge localized modes (ELMs) on the plasma potential and ion impinging energy on the divertor target. A double-peak structure of the ion impinging energy has been identified under JET-relevant ELM conditions. The ELM burst leads to a strong increase in the potential drop in front of the target plate, which accelerates the cold ions from the downstream divertor and accordingly causes a peak value of ion impinging energy. Moreover, the great potential drop helps confine the fast electrons and leads to a reduction in the potential drop and ion impinging energy. The arrival of the upstream hot ions results in the second peak value of ion impinging energy. The maximum potential drop and ion impact energy show a linear dependence on the pedestal temperature. Further, a nonlinear dependence of the peak potential drop and ion impact energy on the ELM energy loss can be ascertained based on the PIC simulations.

Estimation of tungsten production from the upper divertor in EAST during edge localized modes**Project supported by the National Key R&D Program of China (Grant Nos. 2017YFE0300400 and 2017YFE0301300), the National Natural Science Foundation of China (Grant Nos. 11475223 and 11775257), the National Magnetic Confinement Fusion Science Program of China (Grant No. 2015GB101003), and also partly supported by AHNSF of China

During edge localized modes (ELMs), the sheath evolution in front of the Experimental Advanced Superconducting Tokamak (EAST) upper divertor is studied to estimate the sputtered tungsten (W) atoms from the divertor target. A large potential drop across the sheath is formed during ELMs by compared with inter-ELMs, and the maximum of sheath potential drop can exceed one thousand of eV in current EAST operation. Due to the enhancement of the sheath potential drop during ELMs, the W physical sputtering yield from the deuterium (D) ions and the impurity ions on the upper divertor target is found to be significant. It is established that the sputtered W yield during ELMs is at least higher by an order of magnitude than inter-ELMs, and D ions and carbon (C) ions are the main ions governing the W production for the current H-mode with ELMs discharges. With increase in the pedestal electron temperature, the maximum of the D and C ion impact energy during ELMs shows a nearly linear increase, and the D ions have sufficient impact energy to cause the strong W physical sputtering. As a consequence, the D ions may dominate the sputtered W flux from the divertor target when the C concentration is controlled less than one percent for the higher heating power H-mode with ELM discharges in near future.