  Impurity temperature screening is a favorable neoclassical phenomenon
involving an outward radial flux of impurity ions from the core of fusion
devices. Quasisymmetric magnetic fields are intrinsically ambipolar and give
rise to temperature screening for low enough $\eta^{-1}\equiv d\ln n/d\ln T$.
In contrast, neoclassical fluxes in generic stellarators will depend on the
radial electric field, which is predicted to be inward for ion-root plasmas,
potentially leading to impurity accumulation. Here we examine the impurity
particle flux in a number of approximately quasisymmetric stellarator
configurations and parameter regimes while varying the amount of
symmetry-breaking in the magnetic field. The amount of symmetry-breaking (from
perfect quasisymmetry) required to lose temperature screening is dependent on
the particular configuration and flux surface. Neoclassical fluxes have been
obtained using the SFINCS drift-kinetic equation solver with the electrostatic
potential $\Phi=\Phi(\psi)$, where $\psi$ is a flux surface label. Results
indicate that achieving temperature screening is possible, but unlikely, at
reactor-relevant conditions in the core. Thus, the small departures from
symmetry in nominally quasisymmetric stellarators are large enough to
significantly alter the neoclassical impurity transport. A further look at
these fluxes when compared to a gyro-Bohm turbulence estimate suggests that
neoclassical fluxes are small in optimized configurations compared to turbulent
fluxes. Therefore, although neoclassical impurity accumulation is expected in
most situations, the strength of the turbulence may render this irrelevant.
