Investigating the impact of the molecular charge-exchange rate on detached SOLPS-ITER simulations

TL;DR

This study demonstrates that disabling isotope mass rescaling in SOLPS-ITER simulations significantly increases molecular ion content and improves agreement with experimental observations of detachment phenomena in tokamak plasmas.

Contribution

The paper introduces a modified simulation approach that disables isotope mass rescaling, revealing its impact on molecular ion dynamics and plasma detachment modeling.

Findings

Increased D2+ content by over 100 times with modified rate setup.

Enhanced molecular activated recombination and dissociation effects.

Better agreement with experimental spectroscopic measurements.

Abstract

Plasma-molecular interactions generate molecular ions which react with the plasma and contribute to detachment through molecular activated recombination (MAR), reducing the ion target flux, and molecular activated dissociation (MAD), both of which create excited atoms. Hydrogenic emission from these atoms have been detected experimentally in detached TCV, JET and MAST-U deuterium plasmas. The TCV findings, however, were in disagreement with SOLPS-ITER simulations for deuterium indicating a molecular ion density ($D_2^+$) that was insufficient to lead to significant hydrogenic emission, which was attributed to underestimates of the molecular charge exchange rate ($D_2 + D^+ \rightarrow D_2^+ + D$) for deuterium (obtained by rescaling the hydrogen rates by their isotope mass). In this work, we have performed new SOLPS-ITER simulations with the default rate setup and a modified rate setup where ion isotope mass rescaling was disabled. This increased the $D_2^+$ content by $> \times 100$. By disabling ion isotope mass rescaling: 1) the total ion sinks are more than doubled due to the inclusion of MAR; 2) the additional MAR causes the ion target flux to roll-over during detachment; 3) the total $D\alpha$ emission in the divertor increases during deep detachment by roughly a factor four; 4) the neutral atom density in the divertor is doubled due to MAD, leading to a 50\% increase in neutral pressure; 5) total hydrogenic power loss is increased by up to 60\% due to MAD. These differences result in an improved agreement between the experiment and the simulations in terms of spectroscopic measurements, ion source/sink inferences and the occurrence of an ion target flux roll-over.