Spectropolarimetric modeling of interacting Type II supernovae. Application to early-time observations of SN1998S

High-cadence surveys of the sky are revealing that a large fraction of red-supergiant (RSG) stars, which are progenitors of Type II-Plateau (II-P) supernovae (SNe), explode within circumstellar material (CSM). Such SNe II-P/CSM exhibit considerable diversity, with interaction signatures lasting from hours to days, potentially merging with the Type IIn subclass for which longer-duration interaction typically occurs. To tackle this growing sample of transients and to understand the pre-SN mass loss histories of RSGs, we train on the highest quality, spectropolarimetric observations of a young Type IIn SN taken to date: Those of SN1998S at ~5d after explosion. We design an approach based on a combination of radiation hydrodynamics with HERACLES and polarized radiative transfer with CMFGEN and LONG_POL. The adopted asymmetries are based on a latitudinal, depth- and time-independent, scaling of the density of 1D models of SNe II-P/CSM (e.g., model r1w6b with a `wind' mass-loss rate of 0.01Msun/yr used for SN2023ixf). For a pole-to-equator density ratio of five, we find that the polarization reaches, and then remains for days, at a maximum value of 1.0, 1.4, and 1.8% as the CSM extent is changed from 6, to 8 and 10x10^14cm. The polarization is independent of wavelength away from funnel-shaped depolarizations within emission lines. Our models implicate a significant depolarization at line cores, which we use to constrain the interstellar polarization of SN1998S. Our 2D, prolate ejecta models with moderate asymmetry match well the spectropolarimetric observations of SN1998S at 5d, supporting a polarization level of about ~2%. This study provides a framework for interpreting future spectropolarimetric observations of SNe II-P/CSM and SNe IIn and fostering a better understanding of the origin of their pre-SN mass loss.