This meeting will evaluate and recommend fundamental data concerning tungsten, hydrogen, their ions and molecules in the edge plasma region of experimental nuclear fusion devices with a view to quantifying and reducing the uncertainties in the modelling of its collisional, radiative and plasma-material interaction properties.
Research scientists related nuclear fusion, at any stage of their career, who are involved in the quantitative and qualitative study, analysis, modelling, simulation or prediction of collisional-radiative properties or plasma-wall interactions of tungsten and hydrogen in fusion devices, using or producing fundamental data on these processes. Experimental and theoretical approaches are welcomed.
The IAEA is committed to the promotion of gender equality in all its activities and qualified women scientists are particularly encouraged to participate in the meeting.
All persons wishing to participate the event must be designated by an IAEA Member State or should be a member of an organization that has been invited to attend.
For registration, please see information in Participation and Registration.
Quick guide can be found here.
Participants wishing to give an oral presentation or a poster presentation, are expected to submit an abstract through Call for Abstracts.
Guidelines for preparing and submitting slides and posters are given in Presentation Guidelines and Poster Guidelines.
IAEA does not share meeting participant information with third parties for accommodation booking or for any other purposes. If participants ever get any email messages that make it seem like their information is being shared, they should
Experiments 1: tokamaks, linear plasmas, stellarators
Erosion of the ITER main chamber first wall (FW) beryllium (Be) armour is expected to affect the lifetime of FW panel, dust formation, and release of Be impurities potentially leading to enhanced sputtering of the W divertor and tritium retention due to co-deposition. The ITER FW panels are shaped to protect leading edges and misalignments, which leads to magnetically shadowed regions, where impurity re-deposition and fuel co-deposition can occur. In addition to plasma ions, charge exchange (CX) neutrals in ITER and future reactors will play an important role on the first wall erosion and overall fuel retention, but the extent to which they contribute is still unknown.
A low-energy neutral particle analyzer (LENPA) based on the time of- flight method has been developed on EAST to measure the flux and energy of neutral particles to the first wall. The LENPA works in the pulse-counting mode and the energy distribution of neutrals is obtained based on the flight-time of particles. A 220 mm diameter chopper disc with 32 equally spaced slots has a maximum rotating speed of 300 Hz, resulting in a LENPA detecting period of 104 μs and a slit half-open time of 1.07 μs. After freely streaming 4170 mm, the chopped neutral particles reach the detector located at the end of LENPA. In the LENPA detection range, more than 85 % of neutral particles are in the energy range of 20–1000 eV. The integrated neutral flux in the energy range of 20–1000 eV increases with line-averaged density in ohmic discharges due to the increased CX reaction rates in a higher density plasma. Compare to the discharges fuelled with SMBI, the discharges without SMBI have more lower energy neutrals at the similar line-averaged density due to higher edge neutral density and proximity of fuelling. Due to the deeper penetration depth by SMBI fuelling, the neutral particles are generated closer to the core plasma and have higher energy. It was found that the neutral flux increases with heating power in all energy range. The mean energy of neutral particles increases significantly compared to the ohmic discharges. The higher ion temperature and edge density in auxiliary heated discharges result in higher neutral flux in all the energy range.
A quartz crystal microbalance (QMB) was installed together with the LENPA system on the equatorial port of EAST [2]. The neutral-induced material erosion rates and the neutral energy spectrum were measured simultaneously by the two real-time and in-situ diagnostics. The neutral-induced Al erosion rates for the 11 long-pulse full discharges are given by experimental measurement from the QMB and the theoretical calculations according to the neutral energy spectrum from the LENPA, which are consistent with each other. It is proved that higher density and heating power can increase the flux and energy of neutral particles, which results in stronger neutral-induced material erosion. The real-time Li powder injection during discharges can reduce the erosion rate due to the lower recycling and the resulting lower neutral flux.
Special session for Collisional Radiative Models (CR models)
Collisional radiative (CR) models represent the most versatile type of population models used for determining plasma parameters or for predicting for known plasma parameters the plasma behavior. In plasmas used in fusion research the particle temperatures and densities can cover a wide range. Typically, the edge plasma of fusion devices is much cooler than the core plasma, resulting in the presence of atoms and molecules and thus a high complexity of the reaction kinetics. Together with the presence of strong density gradients, processes like Molecular Assisted Recombination or the transition from an ionizing to a recombining plasma regime can play a crucial role. The situation is similar in the negative hydrogen ion sources for neutral beam injection (NBI) at ITER. An ionizing plasma (T$_\mathrm{e}$>10 eV) in the driver region is cooled down by a magnetic filter field to a recombining plasma with T$_\mathrm{e}$≈1 eV close to the extraction system. Consequently, CR models for atomic and molecular hydrogen are needed, being precise over a broad parameter range.
While the development of models for atomic hydrogen is mostly finished and well benchmarked models exist, the situation for molecular hydrogen is less definite. The reason is that due to the presence of vibrational and rotational sublevels a huge number of input parameters (reaction probabilities for collisional and radiative processes) are needed. Up to a few years ago mostly models resolving only the electronic states were applied due to the lack of appropriate reaction probabilities even for electron collision excitation. Even within the input used for these models large inconsistencies were present. Recently, some gaps in the available set of electronically and vibrationally resolved excitation cross sections for molecular hydrogen have been filled, shifting the main point of interest from electronically resolved models towards vibrationally or even ro-vibrationally resolved models. The presence of deuterium and tritium in fusion-relevant plasmas makes the development of molecular CR models desirable also for isotopomeres of H$_2$, further increasing the data needs.
Starting point for the present investigations are the well-known Yacora CR models for molecular hydrogen. A combined approach is followed, using three models with different levels of detail: an electronically resolved model for basic investigations on the kinetics of excitation and de- excitation of electronic molecular states, a vibrationally resolved model mainly for the reaction kinetics within the manifold of vibrational states within the electronic ground state X$^1$ and their impact on electronically excited states and a ro-vibrationally resolved Corona model for the Fulcher band transition (d$^3$→a$^3$). Final aim are well-benchmarked ro-vibrationally resolved models that can be applied to selected molecular optical emission bands over a broad range of plasma regimes. The current status of the models and the process of benchmarking them versus experimental data is introduced and current data needs are pointed out and discussed.
Tungsten used as a plasma-facing material for divertor in fusion devices is sputtered and transported into plasmas. The tungsten behaviors in divertor, edge, and core plasmas should be examined by spectroscopic methods, and the spectroscopic models and atomic data for a wide charge state range of tungsten ions are required to analyze tungsten spectra. Many studies have been done on tungsten ions experimentally and theoretically. Extreme ultraviolet (EUV) spectra of tungsten ions have been examined by comparison of measured spectra in fusion devices and electron beam ion traps (e.g. [1,2]) and calculated spectra by collisional-radiative (CR) models (e.g.[3,4]). Wide peaked spectral feature, so-called the unresolved transition array (UTA) at 4.5-7nm wavelength region is found in plasma with electron temperature ~1keV and is known as numerous overlapped 4d-4f and 4p-4d transitions of tungsten ions. Many little-wide peaks at 2-4 nm are produced mainly by 4g-5f and 4g-6f transitions of Wq+ with q=22-30. These peaks are useful for estimating the charge state distribution and behaviors of these ions in plasmas [4,5]. For ions with q<22, no peaks are found in the 2-4 nm region, and we need to find some spectral peaks for ions with q<22.
We extend our study to EUV spectra at 10-30 nm where n=5-5 transitions of mid-charged tungsten ions are found. We have performed plasma experiments to measure tungsten spectra by pellet injection into Large Helical Device plasmas for this wavelength region. We also try to extend our CR model for tungsten ions down to q=20. Details of the comparison will be presented at the conference.
[1] R. Neu et al., Plasma Phys. Control. Fusion 38, A165 (1996)
[2] H. A. Sakaue et al., AIP Conf. Proc. 1438, 91 (2012)
[3] T. Putterich, R. Neu, R. Dux et al., Plasma Phys. Control. Fusion 50, 085016. (2008)
[4] I. Murakami, H. A. Sakaue, C. Suzuki et al., Nucl. Fusion 55, 093016 (2015)
[5] I. Murakami, D. Kato, T. Oishi et al., Nucl. Mater. Energy 26, 100923 (2021)
Experiments 1: tokamaks, linear plasmas, stellarators
Fulcher Bands of the hydrogenic molecular spectra are complex triplet bands in the visible region (590-630 nm). The electronic transition is identified as d (3p$^3$Π$_u$) at approximately 14 eV level energy to a (2s$^3$Σ$^+_g$) at 12 eV, but this is not the sole electronic transition in this region, as the g→c (3d $^3$Σ$^+_g$ → 2p$^3$Π$_u$) transition is present in this wavelength region as well. Lines from both of those bands have been identified e. g. in the experimental measurements from divertor regions of the JET tokamak. The d→a transition is widely used to learn about molecular behavior in the divertor plasma, but the spectral analysis of such a complex spectrum is challenging. There are large catalogues of experimentally verified wavelength of the rotationally resolved wavelengths of the electronic-vibrational-rotational transitions in the homonuclear H$_2$ [1] and D$_2$ [2] molecular spectra (also some for the heteronuclear HD). They include also those of the g→c interfering band. Unfortunately, experimental data even for main FB transition are sparse for tritium-containing molecules, and most of the diatomic constants are either calculated using isotope effect or from relatively restricted and old research published by Dieke and Tomkins in 1950. Even for homonuclear tritium the data are restricted to first several (at most ten) rotational transitions for each vibrational transition. In the TEXTOR [3] and JET-C deuterium [4] and tritium/DT [5] experiments the higher rotational transitions were not observed, but in JET-ILW we see a change of ro-vibrational population depending on local injection and recycling at the target places which was not seen in previous experiments. Those conditions result in a rotational overpopulation of the first main-diagonal band, where in the case of high molecular density the observed rotational temperature can even reach the electron temperature. In such cases lines up to J= 17 can be easily observed (e.g. [6] in D$_2$, [7] in H$_2$), which adds to the complexity of the analysis, both because the wavelengths for such lines may not be following low-order diatomic expansion values, and because the spectra from different vibrational contributions do overlap. For heteronuclear molecules the analysis even more challenging, as in this case there is an overlap between spectra from at least one homonuclear molecule and a heteronuclear one (see e.g. [8]).
In this presentation will be shown experimental spectra from JET tokamak divertor measurements, containing T$_2$, HT and DT molecular contributions. The existing molecular data and resulting line identifications and estimations of energy distribution functions will be discussed and compared with available results from other isotope spectra.
[1] H.M. Crosswhite, G.H. Dieke, eds., The Hydrogen Molecule Wavelength Tables of Gerhard Heinrich Dieke, Wiley-Interscience, New York, 1972.
[2] B.P. Lavrov, I.S. Umrikhin, JQSRT 182 (2016) 180–192.
[3] S. Brezinsek et al, Journal of Nuclear Materials 313–316 (2003) 967–971.
[4] A. Pospieszczyk et al , Journal of Nuclear Materials 337–339 (2005) 500–504.
[5] A. Pospieszczyk et al, Journal of Nuclear Materials 363–365 (2007) 811–815.
[6] G. Sergienko et al, Journal of Nuclear Materials 438 (2013) S1100–S1103.
[7] E. Pawelec et al, Europhysics Conference Abstracts, 46A, (2022) O2.J501
[8] S. Brezinsek et al, Physica Scripta T103 (2003) 63.
Modelling 1: tokamaks, linear plasmas, stellarators
A summary on the numerical assessment of detailed distribution of neutral particles impinging specific diagnostic surfaces in ITER (i.e. mirrors on upper-outer limiter or midplane) was given. Demonstrated are the EIRENE post-processing of ITER SOLPS4.3 simulations including an extension of plasma on an artificial grids up to first wall (c.f. based on previous simulations by A. Khan NME 2019). Some modifications in EIRENE were required to also have both energy and angular distributions were implemented (for the poloidal angle $\alpha$ on a given surface). As a main result it was found that detailed distributions give $2-3$ larger D-on-W physical sputter yields $⟨𝒀(𝑬,\cos(\alpha))⟩$ compared to standard estimates $𝒀(⟨𝑬⟩)$. The actual factor depends on far-SOL assumptions used for the extension of the grid to the wall, as well on H/L-mode. The $\cos(\alpha)$-dependence gives a factor $2$. In H-mode the main contribution comes from the tail of the energy distribution (physical sputtering is suppressed for particle below the threshold $~200$ eV).
As a next step Ne-on-W calculation should be pursued and compared w.r.t its relevance to D-on-W. Also, refined SOLPS-ITER with wide-grid option should provide a better picture. So far only uncorrelated energy and angular distributions collected and an extension to multi-variate distribution functions possible $𝑓=(𝐸, \cos(\alpha))$ should be pursued (requiring longer EIRENE run-times for improved statistics and requires large memory). Data compression through MaxEnt regularization techniques are proposed, also to find a better figure-of-merit when mixed (i.e. multi-component) energy and angular distributions are saturated.
Measurements in JET ITER-like Wall low-confinement mode plasmas showed that the scrape-off layer at the low-field side midplane was $30 \%$ broader and the low-field side divertor plasma a factor of $2-3$ more strongly detached in pure-tritium than in pure-protium plasmas for otherwise similar core plasma conditions. The conditions in pure-deuterium and in a $40 \% - 60 \%$ deuterium-tritium plasmas were found to be between those in protium and tritium plasmas. These results are statistically significant and consistent with the lower triton than proton velocities and shorter ionisation mean free path for tritium than for protium atoms, as predicted by the EDGE2D-EIRENE code. For all four isotope species and mixes, the onset of detachment was observed within $10 \%$ of the same electron density at the edge of the core plasma. The formation of detachment with increasing edge density is identical between the isotope species and mixes: the electron temperature at and adjacent to the low-field side target plate decreased from $20$ eV to $1$ eV in non-linear fashion over a narrow edge density range, and the onset of detachment was achieved when the electron temperature at the low-field side plate reached $2$ eV. The high electron density (pressure) region, formed at the low-field side plate in the near-separatrix SOL when the electron temperature was $0.7-1.0$ eV, was observed to transition non-linearly from the target plate to the low-field side X-point region. The formation of the high-density/pressure region in the low-field side divertor was independent of the molecular gas pressure in the sub-divertor, and thus the fuelling and pumping rates of the experimental setup. The density limit was $30 \%$ lower in tritium than in protium plasmas. The potential impact of hydrogenic collisional-radiative models including the isotope mass and Lyman-alpha opacity will be discussed.
Special session for Collisional Radiative Models (CR models)
The need for quantitative spectroscopy of molecular hydrogen has experienced substantial demand, leading to the accumulation of diverse elementary-processes data encompassing radiative transitions, electron-impact transitions, and predissociations. In this study, we attempt an experimental validation of this dataset by comparing the vibronic populations across multiple molecular hydrogen levels: EF1Σ+g, H1Σ+g, D1Π±u, GK1Σ+g, I1Π±g , J1∆±g , h3Σ+g, e3Σ+u, d3Π±u, g3Σ+g , i3Π±g, and j3∆±g , measured from thousands of emission lines observed with a custom-made Echelle spectrometer. Our analysis incorporates a collisional-radiative model (CRM) that relies on the most up-to-date dataset and on spectrally observed populations from Large Helical Device (LHD) plasmas. Remarkably, we find that predissociation, which is the spontaneous dissociation process of excited molecules, is crucial for accurately replicating the observed outcomes, although many of the previously-reported CRM for molecular hydrogen have neglected this effect. We also demonstrate the possibility of a new divertor diagnostic technique utilizing the measured molecular spectrum interpreted through the CRM with predissociation included. This research also highlights that incorporation of predissociation data for hydrogen molecular isotopologues ($\mathrm{H}_2$, $\mathrm{H}_2$, $\mathrm{T}_2$, HD, HT, DT) is essential for the interpretation of future burning-fusion plasma spectroscopic diagnostics.
Modelling 1: tokamaks, linear plasmas, stellarators
Kinetic sputtering of materials during bombardment with energetic particles is a well understood process for the case of perfectly flat surfaces. Both analytical models and numerical simulation codes have been established in the last decades showing good agreement for quantities accessible by experiments like the sputtering yield. Still, for many applications, e.g., plasma-facing materials in a fusion device, the assumption of a perfectly flat surface may fail to account for the important role of surface roughness.
A range of effects can be introduced by the geometry of a rough surface, e.g., redeposition of sputtered particles, shadowing, or the influence of locally tilted surface segments on the incidence angle of incoming projectiles. Thus, the assessment of sputtering yields for rough surfaces remains a difficult task.
During the last years, the effect of roughness on W sputtering has been comprehensively investigated at TU Wien by conducting laboratory experiments with a quasi-non-invasive low-flux ion source and a high precision quartz crystal microbalance. Furthermore, the development of a binary collision approximation - based raytracing code called SPRAY allowed for quick calculation of sputtering yield data, where large-scale AFM images of rough W surfaces could be considered as input. It was possible to identify the mean inclination angle δm as a scale-independent roughness parameter that allows consistent prediction of sputtering yields for rough W surfaces [1,2].
Furthermore, oriented and periodically patterned topographies like nano-columnar W structures have also been investigated. In comparison to flat W surfaces, nano-columnar W shows a substantial reduction in the sputtering yield, which appears interesting for applications in nuclear fusion [3,4,5].
In this invited talk, an overview and explanation of roughness effects on W sputtering will be presented, including results for conventionally rough and nano-columnar W surfaces.
[1] C. Cupak, et al., Appl. Surf. Sci. 570 (2021) 151204
[2] P.S. Szabo et al., Surf. Interfaces. 30 (2022) 101924
[3] A. Lopez-Cazalilla, C. Cupak et al., Phys. Rev. Mat. 6 (2022) 075402
[4] C. Cupak et al., Phys. Rev. Mat. 7 (2023) 065406
[5] J. Brötzner et al., Nucl. Mater. Energy; in press (2023), 101507
(https://doi.org/10.1016/j.nme.2023.101507)
Molecular Dynamics simulations is a powerful tool in investigating the effects on surfaces under various conditions. It is especially effective for studying the interaction between energetic particles and surfaces. Due to relevancy to fusion, tungsten surfaces have been studied under various conditions. Both pristine surfaces and atomistically rough surfaces have been studied under irradiation by various ion species and energies [1,2]. Simulations have in the last years revealed a plethora of insight into how different surface features will affect the sputtering and evolve under continuous irradiation [3,4]. It has been shown that the sputtering at high energies is directly related to the channeling in a certain direction of the incoming ion, which renders some of the simulation tools assuming amorphous materials questionable [5].
Another important factor under fusion conditions is the effect of hydrogen or its isotopes that are impacting the surface, causing sputtering, or get implantated into the surface layers. We have found that the sputtering by light elements, like hydrogen or its isotopes, is more complicated than previously thought. We found that the sputtering can happen due to several different mechanisms, which are not seen for heavier elements. We have also studied the effect of having the surface saturated with deuterium on the sputtering by both light and heavy ions. Again, we see that having a saturated surface versus a pristine one will affect the results.
All obtained results can be used as input in higher scale models in order to further predict their impact on material lifetimes and plasma properties.
[1] Jussila, J., Granberg, F., and Nordlund, K. Nuclear Materials and Energy. 17 (2018) 113-122
[2] Lopez-Cazalilla, A., Jussila, J., Nordlund, K., and Granberg, F. Computational Materials Science 216 (2023) 111876
[3] Lopez-Cazalilla, A., Cupak, C., Fellinger, M., Granberg, F., Szabo, P. S., Mutzke, A., Nordlund, K., Aumayr, F., and Gonzalez-Arrabal, R. Physical Review Materials 6 (2022) 075402
[4] Cupak, C., Biber, H., Brötzner, J., Fellinger, M., Brandstätter, F., Aumayr, F., Lopez-Cazalilla, A., Granberg, F., Nordlund, K., Szabo, P. S., Mutzke, A., and González-Arrabal, R. Physical Review Materials 7 (2023) 065406
[5] Schlueter, K., Nordlund, K., Hobler, G., Balden, M., Granberg, F., Flinck, O., da Silva, T. F. and Neu, R. Physical Review Letters 125 (2020) 225502
Modelling 2: A+M data, cross-sections, processes
In this talk will be presented our last results on electron-H$_2$ and -D$_2$ collisions. We will focus in particular on the vibrational-excitation, dissociative-attachment, dissociative-recombination and dissociative-excitation processes rotationally and vibrationally resolved.
Experiments 1: tokamaks, linear plasmas, stellarators
Upper and lower graphite divertors in EAST tokamak were replaced by tungsten divertors in 2014 and 2021, respectively, to improve plasma performance in long pulse discharges and accumulate knowledges on the tungsten divertor operation. Studies on the tungsten behavior in edge and core plasmas are then crucially important for improving the plasma performance. For the purpose four fast-time-response [1-2] and four space-resolved [3-4] extreme ultraviolet (EUV) spectrometers have been installed on EAST to observe line emissions from tungsten ions and their radial intensity profiles in wide wavelength range of 5-520Å.
Photon emission coefficient (PEC) data for W43+ at 61.334 and 126.29Å and W45+ at 62.336 and 126.998Å have been used to estimate the density profile of W43+ and W45+ions in the core plasma [3,5]. Tungsten unresolved transition arrays (W-UTA) emitted in long wavelength ranges of 168-225Å, 225-268Å and 278-332Å are analyzed for the study of tungsten behaviors in the edge plasma. As a result, three lines of 186.28 Å, 190.48 Å and 192.02 Å with relatively strong intensities were identified as line emissions from W8+ ions by comparing the time behavior of line emissions from W6+ at 216.219 and 261.387 Å [6,7] and W7+ at 200.367 Å and 200.483 Å [8] [1-4], which are well known as the EUV line emission from low-ionized tungsten ions. The number of ionization events per photon (or inverse photon efficiency), S/XB, for line emissions from low-ionized ions are strongly required to evaluate the tungsten influx rate at plasma edge. A visible spectrometer aimed at spatial profile measurement covering the whole poloidal cross section has been newly installed on EAST tokamak for investigation of atomic tungsten and the M1 transition from tungsten ions, which are observed for W26+-W28+ ions in LHD [9] and W8+-W12+ ions in EBIT [10]. Analyses of full radial density profiles of highly ionized ions tungsten ions and influx rate of low-ionized tungsten ions are attempted based on the PEC and S/XB data.
Thailand is one of the countries that recognizes the importance of fusion technology and has a national plan for its implementation. The Thailand Institute of Nuclear Technology (TINT) has been engaged in the development of a tokamak device known as Thailand Tokamak-1 (TT-1) since 2015. This device is derived from the HT-6M, a previous tokamak device developed by the Chinese Academy of Science (ASIPP) in China. The TT-1 tokamak is scheduled to be the inaugural fusion device to be installed at the Thailand Institute of Nuclear Technology (TINT), located at the Ongkharak location in Nakhonnayok province. The objective of the TT-1 is to provide researchers from Thailand and surrounding areas an opportunity to access a facility that facilitates scientific investigation of fusion technology and design. Additionally, it provides a platform for human resource training in plasma control, experimental planning, tokamak operation, and the conduct of fusion-related research. The TT-1 exhibits a compact form, characterized by a major radius (R) of 0.65 m and a minor radius (a) of 0.20 m. The plasma current (Ip), and the toroidal magnetic field strength (Bt) have the potential to be enhanced to values of 150 kA and 1.52 T, respectively. The main parts of the device were assembled at ASIPP before being shipped and installed at Thailand Institute of Nuclear Technology (TINT). In terms of hydrogen plasma discharge, it will be initiated with ohmic heating in the first phase of plasma operation, which is expected to start in 2024. In the next operating phase, the TT-1 will be equipped with the neutral beam injection (NBI) system in order to generate plasmas characterized by high temperatures and density. Based on the on-site discussion of the device design, the diagnostics that will be developed and installed as planned for the first phase are introduced in this paper. In addition, a schematic overview of the sites for installing the diagnostics is reported.
Modelling 2: A+M data, cross-sections, processes
W generation and transport study represents one of the hottest topics for next generation tokamaks with full metal plasma-facing components [1-4]. The existing W sputtering models are based on low density plasma sheath approximation [5-7], which might be not applicable for high density divertor plasmas expected in future machines, e.g. like ITER and DEMO.
In the given presentation we discuss results of recent full kinetic modelling of W sputtering and re-deposition in a high-density plasma sheath. The simulations indicate that
i) using atomic data based on the coronal approximation results at negligible prompt re-deposition of the W ions. This might have important consequences, as the corresponding net sputtering rates (i.e. gross sputtering rates minus prompts re-deposition) might become significant and cannot be neglected as it was assumed up to now.
ii) the coronal approximation is not sufficient for such study and multiple atomic transitions have to be included in the model. We identified the corresponding atomic processes which have to be included into the W sputtering models in a high density plasma.
Based on above-listed results, we present a list of W-relevant data necessary for realistic estimation of the W net sputtering rated in future fusion machines.
[1] T. Pütterich, et al., Nucl. Fusion 50 (2010) 025012
[2] R. Dux, et al., Nucl. Fusion, 51 (2011) 053002.
[3] S. Brezinsek, et al., Phys. Scr., T145 (2011) 014016.
[4] R.A. Pitts, et al., Nucl. Mater. and Energy, 20 (2019) 100696.
[5] G.J. van Rooij, et al., J. Nucl. Mater., 438 (2013) S42–S47.
[6] D. Tskhakaya, et al., J. of Nucl. Mater., 463 (2015) 624-628.
[7] A.V. Chankin, et al., Plasma Phys. Control. Fusion, 56 (2014) 025003.
Electron-molecule cross sections resolved in the rovibrational levels of the target are required for determining the properties and dynamics of many low-temperature plasmas. We have applied the Molecular Convergent Close-Coupling (MCCC) method to produce a comprehensive set of vibrationally-resolved cross sections for electron collisions with molecular hydrogen and its isotopologues comprised of more than 60,000 entries. This complete collision data set is available to the research community via the dedicated MCCC database (mccc-db.org). For H$_2$ the data set includes transitions from all bound vibrational levels of the ground electronic state to all vibrational levels of 18 excited electronic states (all states in the $n=2,3$ shells) with similar datasets for each of the isotopologues. We are working on producing a vibrationally resolved set of cross sections for scattering on the electronically excited $c~^3\Pi_u$, $a~^3\Sigma_g^+$, $B~^1\Sigma_u^+$, $C~^1\Pi_u$, and $EF~^1\Sigma_g^+$ states of all H$_2$ isotopologues with first results (not vibrationally-resolved) already available. Recent calculations have also been performed for vibrationally-resolved ionisation of the $n=1$-2 states, as well as fully rovibrationally-resolved transitions within the ground state and between the ground and $n=2$ states. An overview of the available MCCC dataset will be presented, and examples of its utilisation in fusion and astrophysical plasma modelling will be given.
An understanding of the description of the discharge fuel, either D or He, is essential in order to model processes in the plasma edge and divertor, these processes including transport simulations of the plasma exhaust and the determination of plasma-facing surface erosion. The emission from the fuel species is described by collisional-radiative models. However, discrepancies have been found between the models and observations made in the JET tokamak. Lawson et al. (2022 and 2023a) give detailed measurements of He II (He+) Lyman series line intensities and their ratios and also make comparisons with He II Balmer and Paschen series members. A collisional-radiative model involving routinely occurring opacity is developed in order to explain the observations. Attention is given both to the equilibrium discharge phase when constant line intensity ratios indicate very similar temperatures for the plasma regions emitting the hydrogen-like radiation and exceptional phases in which the line intensity ratios deviate from these near-constant values. In particular, the variation of the electron temperature throughout discharges is illustrated. Increases in temperature during the non-equilibrium, exceptional phases are understood in terms of transient effects of recycling particles. Pulses in which additional opacity beyond that routinely occurring are investigated. Both the routinely occurring opacity and the limited temperature range of the emitting plasmas are expected to have implications for the interpretation of various analyses.
Lawson K D et al., 2022, Report ‘He II Lyman series line intensity measurements in the JET tokamak’
Lawson K D et al., 2023a, Report ‘He II VUV and visible line intensity measurements in the JET tokamak’
*See author list of ‘Overview of JET results for optimizing ITER operation’ by J Mailloux et al., 2022, Nucl. Fusion, 62, 042026
Experiments 2: A+M data, cross-sections, processes
Sputtering is one of the fundamental processes at plasma-facing components in fusion experiments as the lifetime of the components strongly depends on the sputtering rates. Therefore, it is of paramount importance to understand the processes underlying the sputtering, including the energy, angular and atomic level distribution functions of the sputtered tungsten (W). These distribution functions are further essential input parameter for erosion codes.
One aspect of the presented research deals with the investigation of the initial atomic level distribution of sputtered W. This distribution needs to be taken into account for the interpretation of spectroscopic data and further the determination of the net erosion via so-called S/XB values in fusion devices [1]. However, this initial atomic level population within the fivefold ground term $^5$D$_{0-4}$ and the metastable $^7$S$_3$ level of sputtered W remains an open question. On the one hand, experiments in the tokamak TEXTOR suggest a nonphysical effective temperature to describe the atomic level population distribution via a Boltzmann distribution [2]. This results in a strong population also of other levels and not only of the ground state. On the other hand, ion beam experiments for different materials show a strong population of the ground state only [3]. Using a new approach, we studied this open question in the low density and temperature plasma of the linear plasma device PSI-2. Via an imaging spectrometer with a high spatial resolution (50 µm/pixel) the temporal line intensity development as a function of distance to the target was studied for different transitions of W I. In the experiments, the W atoms were sputtered by mono-energetic Ar$^+$ ions at 80 eV. The target was cooled using water-cooling and different target temperatures were investigated.
A second aspect of the presented research deals with investigating the energy and angular distribution functions. These strongly influence the transport into the plasma and, moreover, the re-deposition of sputtered material on the inner vessel. Via ion beam experiments, these distribution functions are well accessible for high impact energies and can be described by a cosine angular distribution and the so-called Thompson energy distribution [4]. However, for the impact energies relevant in fusion research, of up to several hundred electron volts, the data is only rare and deviations from the distribution functions are reported [2]. To investigate these deviations and close the remaining gap of information, we carried out experiments and subsequently modeled the line shape emitted by sputtered W atoms via a Doppler-shifted emission model. The experiments were carried out in the linear plasma device PSI-2, where W samples were exposed to low density ($\text{n}_\text{e} ≈ 2 \times 10^{12} \text{cm}^{-3}$) and temperature ($\text{T}_\text{e} ≈ 3 \text{eV}$) argon plasmas. The ions hitting the sample were accelerated due to biasing the sample to mono-energetic impact energies in the order of 100 eV. The light emitted by the sputtered atoms was detected via a high-resolution spectrometer with a resolving power of $ \lambda/\Delta\lambda \approx7 \times 10^5$. Until now, the finite size of the targets has not been taken into account in the analysis of the observed spectral lines during sputtering. The standard emission model used so far, which assumes an infinite target or point source, fails for observations with two lines of sight. This gap could be removed due to the new finite-size model.
[1] S. Brezinsek et al., Physica Scripta, T170, 014 052, 2017.
[2] I. Beigman, et al., Plasma Phys. Control. Fusion, 49, 1833-1847, 2007.
[3] A. P. Yalin et al., Applied Optics, 44, 6496-6505, 2005.
[4] R. Behrisch and W. Eckstein, Sputtering by Particle Bombardment, ser. Topics in
Applied Physics. Berlin, Heidelberg: Springer Berlin Heidelberg, 110, 2007.
[5] A. Goehlich et al., Journal of Nuclear Materials, 266, 501-506, 1999.
Tungsten is being considered as a plasma-facing material in magnetically confined fusion devices, such as ITER. Electron collision ionization is a dominant atomic process in fusion plasma which determines the ionization balance of the non-local thermal equilibrium plasmas. Despite great effort have been pain the experimental measurements and theoretical calculations, however, the effect of long-lived excited states in low charged ionic stages need to be investigated[1]. Moreover, reliable EISI data are not available for many tungsten ions[2]. Last but not least, suprathermal electron influence and non-Maxwellian rate coefficients of high charged W ions remain unclear.
Therefore, we investigate the contributions of metastable states and non-Maxwellian EEDF to electron-impact ionization of tungsten ions, and aim to provied accurate electron-impact ionization rate coefficients of tungsten ions. Comparison between the previous experimental measurement results and present calculation show a prominent contribution of metastable states in low charged states such as W7+-W10+ ions[3-4]. Moreover, we performed calculations of detailed electron-impact single ionization cross sections for tungsten ions, spanning charge states W38+-W45+[5]. We demonstrate the influence of non-Maxwellian distribution on the rate coefficient of the W46+-W55+[6]. The data obtained are expected to be useful for modelling plasmas for fusion applications.
[1] 1st IAEA Technical Meeting on Tungsten and Hydrogen in Edge Plasmas meeting, https://amdis.iaea.org/meetings/tm-tungsten-hydrogen/
[2] I. Murakami, et al., NIFS Database, https://dbshino.nifs.ac.jp/nifsdb/.
[3] L. Chen, et al., J Quant. Spectrosc. Radiat. Tran. 285, 108179 (2022).
[4] R. Bao, et al., Chin. Phys. B 32, 063401 (2023).
[5] R. Bao, et al., Atoms, 10, 92 (2022).
[6] R. Bao, et al., Plasma Phys. Control. Fusion 65, 105004 (2023).
Modelling 1: tokamaks, linear plasmas, stellarators
The reflection and sputtering properties of surfaces under energetic particle impact have been studied theoretically and experimentally since several decades and for many systems (especially planar, mono-elemental and amorphous target-systems) a comprehensive level of knowledge does exist. The situation is considerably less satisfying once the target system exhibits more complex features. Here we outline the present state of the modelling capabilities with a focus on SDTrimSP, their comparison and validation with experimental data and on some unexpected features of Fe-W-systems under irradiation.
Dedicated computer codes are used to address the physics of plasma-wall interaction and related issues of wall lifetime and long-term tritium retention in fusion devices. The Monte-Carlo code ERO2.0 is such an example, which allows numerical modeling of realistic three-dimensional device geometries based on its massive parallelization [1]. In addition to (net) erosion and re-deposition maps on plasma-facing components, the code yields volumetric information about impurity concentrations in the plasma. A fully kinetic approach is used to trace impurities throughout the background plasma, thereby also addressing atomic processes like ionization, recombination, and collisions along the particle trajectory. The code was validated against experimental data from the largest European fusion device JET in the past [2], and predictions were made regarding the wall lifetime for the next step international experimental fusion reactor ITER [3].
In the present contribution, recent code developments, such as the contribution of spatially inhomogeneous multi-species plasma backgrounds to gross erosion calculations will be presented. These developments will then be applied for the first time to modeling the full tungsten European demonstration reactor DEMO, which is currently in the design phase. The ERO2.0 simulations are based on the currently only available DEMO-specific plasma solution [4], which was generated on basis of the 2017 Baseline equilibrium. Due to gaps in the range of (15 - 80) cm between the last simulated plasma grid-point and the first wall, special emphasis will be put on extrapolation assumptions of plasma parameters to the wall and their impact on the transport of tungsten in the scrape-off-layer plasma. The simulations also reveal charge-exchange neutrals to be the driven source of main chamber erosion, while the divertor erosion is mainly driven by seeded Argon ions and tungsten self-sputtering. Beyond, the importance of energy- and angular-resolved charge-exchange spectral characteristics for main chamber erosion calculations will be discussed in detail.
[1] J. Romazanov et al., Physica Scripta T170, 014018 (2017)
[2] J. Romazanov et al., Nuclear Materials and Energy 18, 331-338 (2019)
[3] J. Romazanov et al., Contributions to Plasma Physics 60, e201900149 (2019)
[4] F. Subba et al., Nuclear Fusion 61, 106013 (2021)
Introduction to open issues in the new ITER baseline with a W wall for Q = 10 operation that require experimental assessment.
Experiments 1: tokamaks, linear plasmas, stellarators
Modelling 2: A+M data, cross-sections, processes
Special session for Collisional Radiative Models (CR models)
Modelling 1: tokamaks, linear plasmas, stellarators
Experiments 2: A+M data, cross-sections, processes