Mostly multiband photometric transit observations have been used so far to retrieve broadband transmission spectra of transiting exoplanets in order to study their atmosphere. An alternative method has been proposed and has only been used once to recover transmission spectra using chromatic Rossiter-McLaughlin observations. Stellar activity has been shown to potentially imitate narrow and broadband features in the transmission spectra retrieved from multiband photometric observations; however, there has been no study regarding the influence of stellar activity on the retrieved transmission spectra through chromatic Rossiter-McLaughlin. In this study with the modified SOAP3.0 tool, we consider different types of stellar activity features (spots and plages), and we generated a large number of realistic chromatic Rossiter-McLaughlin curves for different types of planets and stars. We were then able to retrieve their transmission spectra to evaluate the impact of stellar activity on them. We find that chromatic Rossiter-McLaughlin observations are also not immune to stellar activity, which can mimic broadband features, such as Rayleigh scattering slope, in their retrieved transmission spectra. We also find that the influence is independent of the planet radius, orbital orientations, orbital period, and stellar rotation rate. However, more general simulations demonstrate that the probability of mimicking strong broadband features is lower than 25% and that can be mitigated by combining several Rossiter-McLaughlin observations obtained during several transits.

One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit. Stellar activity can alter the shape of photometric transit light curves and in a similar way they can deform the RM signal. In this paper we discuss several RM observations, obtained using the HARPS spectrograph, of known transiting planets that all transit extremely active stars, and by analyzing them individually we assess the variation in the estimated spin-orbit angle. Our results reveal that the estimated spin-orbit angle can vary significantly (up to 42 degrees) from transit to transit, due to variation in the configuration of stellar active regions over different nights. This finding is almost two times larger than the expected variation predicted from simulations. We could not identify any meaningful correlation between the variation of estimated spin-orbit angles and the stellar magnetic activity indicators. We also investigated two possible approaches to mitigate the stellar activity influence on RM observations. The first strategy was based on obtaining several RM observations and folding them to reduce the stellar activity noise. Our results demonstrated that this is a feasible and robust way to overcome this issue. The second approach is based on acquiring simultaneous high-precision short-cadence photometric transit light curves using TRAPPIST/SPECULOOS telescopes, which provide more information about the stellar active region's properties and allow a better RM modeling.