  We study the polarization produced by scattering from dust in a bow
shock-shaped region of enhanced density surrounding a stellar source, using the
Monte Carlo radiative transfer code SLIP. Bow shocks are structures formed by
the interaction of the winds of fast-moving stars with the interstellar medium.
Our previous study focused on the polarization produced in these structures by
electron scattering; we showed that polarization is highly dependent on
inclination angle and that multiple scattering changes the shape and degree of
polarization. In contrast to electron scattering, dust scattering is
wavelength-dependent, which changes the polarization behaviour. Here we explore
different dust particle sizes and compositions and generate polarized spectral
energy distributions for each case. We find that the polarization SED behaviour
depends on the dust composition and grain size. Including dust emission leads
to polarization changes with temperature at higher optical depth in ways that
are sensitive to the orientation of the bow shock. In various scenarios and
under certain assumptions, our simulations can constrain the optical depth and
dust properties of resolved and unresolved bow shock-shaped scattering
regions.Constraints on optical depth can provide estimates of local ISM density
for observed bow shocks. We also study the impact of dust grains filling the
region between the star and bow shock. We see that as the density of dust
between the star and bow shock increases, the resulting polarization is
suppressed for all the optical depth regimes.
