Near-Infrared Polarimetry and H$_2$ emission toward Massive Young Stars: Discovery of a Bipolar Outflow associated to S235 e2s3

We present a near-infrared $H$ band polarimetric study toward the S235 e2s3 protostar, obtained using the POLICAN instrument on the 2.1m OAGH telescope. The images reveal a bipolar outflow with a total length of about 0.5pc. The outflow nebulosity presents a high degree of linear polarization ($\sim80\%$) and reveals a centrosymmetric pattern with the polarization position angles. The polarization characteristics suggest their origin to be single scattering associated with dust in the outflow. Using multiwavelength archival data, we performed spectral energy distribution (SED) fitting based on radiative transfer models of turbulent core accretion theory. The best-fit SED model indicated that the protostar has a mass of $6.8\pm1.2\,M_\odot$, with a disk accretion rate of $3.6\pm1.2\times10^{-4}\,M_\odot\,yr^{-1}$ and a total bolometric luminosity of $9.63\pm2.1\times10^{3}\,L_\odot$. Narrowband H$_2$ ($2.12\,\mu$m) observations show shocked emission along the bipolar lobes tracing the jet's interaction with the surrounding medium. The estimated H$_2$ luminosity of the outflow is $2.3_{-1.3}^{+3.5}\,L_\odot$, which matched the known power-law correlation with the source bolometric luminosity, similar to other high-mass outflows. The orientation of the bipolar outflow was found to be parallel to the local magnetic field direction. The overall results assert the fact that the S235 e2s3 source is a massive young star driving a highly collimated bipolar outflow through disk accretion.