Activity Induced Variation in Spin-Orbit Angles As Derived from Rossiter–Mclaughlin Measurements M
Total Page:16
File Type:pdf, Size:1020Kb
A&A 619, A150 (2018) https://doi.org/10.1051/0004-6361/201833709 Astronomy & © ESO 2018 Astrophysics Activity induced variation in spin-orbit angles as derived from Rossiter–McLaughlin measurements M. Oshagh1,2, A. H. M. J. Triaud3, A. Burdanov4, P. Figueira2,5, A. Reiners1, N. C. Santos2,6, J. Faria2, G. Boue7, R. F. Díaz8,9, S. Dreizler1, S. Boldt1, L. Delrez10, E. Ducrot4, M. Gillon4, A. Guzman Mesa1, E. Jehin4, S. Khalafinejad11,12, S. Kohl11, L. Serrano2, and S. Udry13 1 Institut für Astrophysik, Georg-August Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany e-mail: [email protected] 2 Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal 3 University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK 4 Space Sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, Allée du 6 Août 17, Bat. B5C, 4000 Liège, Belgium 5 European Southern Observatory, Alonso de Cordova 3107, Vitacura Casilla 19001, Santiago 19, Chile 6 Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, 4169-007 Porto, Portugal 7 IMCCE, Observatoire de Paris, UPMC University Paris 6, PSL Research University, Paris, France 8 Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Buenos Aires 1428, Argentina 9 CONICET – Universidad de Buenos Aires, Instituto de Astronomía y Física del Espacio (IAFE), Buenos Aires 1428, Argentina 10 Astrophysics Group, Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, UK 11 Hamburg Observatory, Hamburg University, Gojenbergsweg 112, 21029 Hamburg, Germany 12 Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany 13 Observatoire de Genève, Université de Genève, 51 chemin des Maillettes, 1290 Versoix, Switzerland Received 25 June 2018 / Accepted 17 August 2018 ABSTRACT 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 present 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◦) 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. Key words. planets and satellites: fundamental parameters – planets and satellites: detection – stars: activity – starspots – techniques: radial velocities – techniques: photometric 1. Introduction formation, planets should form on aligned orbits λ ∼ 0 (Pollack et al. 1996; and see a comprehensive review in Dawson & Understanding the orbital architecture of exoplanetary systems Johnson 2018, and references therein), as observed in the case is one of the main objectives of exoplanetary surveys, in order of the solar system. Planetary evolution mechanisms such as to provide new strong constraints on the physical mechanisms disk migration predict that planets should maintain their pri- behind planetary formation and evolution. The sky-projected 1 mordial spin-orbit angle, but the dynamical interactions with a spin-orbit angle (hereafter spin-orbit angle λ) is one of the third body and tidal migration could alter the spin-orbit angle most important parameters of the orbital architecture of plane- (see a comprehensive review in Baruteau et al. 2014, and ref- tary systems. According to the core-accretion paradigm of planet erences therein). Moreover, several studies have shown that the misalignment of a planet could be primordial considering the 1 The angle between the stellar spin axis and normal to the planetary disks being naturally and initially misaligned (e.g., Batygin 2012; orbital plane. Crida & Batygin 2014). In this context, the measurement of Article published by EDP Sciences A150, page 1 of 17 A&A 619, A150 (2018) spin-orbit angles under a wide range of host and planet condi- In this paper we assess the impact of stellar activity on the tions (including multiple transiting planets) is extremely impor- spin-orbit angle estimation by observing several RM observa- tant in order to obtain feedback on modeling. tions of known transiting planets which all transit extremely As a star rotates, the part of its surface that rotates toward active stars. We organize our paper as follows. We describe our the observer is blueshifted and the part that rotates away is targets’ selection and their observations in Sect.2. In Sect.3 we redshifted. During the transit of a planet, the corresponding rota- define our fitting procedure to estimate the spin-orbit angles, and tional velocity of the portion of the stellar disk that is blocked also present its result. We investigate the possibility of mitigating by the planet is removed from the integration of the velocity the influence of stellar activity on RM observations by combin- over the entire star, creating a radial velocity (RV) signal known ing several RM observations in Sect.4. In Sect.5 we present as the Rossiter–McLaughlin (RM) effect (Holt 1893; Rossiter the simultaneous photometric transit observation with our RM 1924; McLaughlin 1924). The RM observation is the most pow- observation and we discuss the advantages of having simultane- erful and efficient technique for estimating the spin-orbit angle λ ous photometric transit observation in eliminating stellar activity of exoplanetary systems or eclipsing binary systems (e.g., Winn noise in RM observations. In Sect.6, we discuss possible con- et al. 2005; Hébrard et al. 2008; Triaud et al. 2010; Albrecht tamination from other sources of stellar noise that affect our et al. 2012, and a comprehensive review in Triaud 2017, and refer- observations, and conclude in Sect7. ences therein). The number of spin-orbit angle λ measurements has dramatically increased over the past decade (200 systems to date)2, and have revealed a large number of unexpected systems, 2. Observations and data reduction which include planets on highly misaligned (e.g., Triaud et al. 2.1. Target selection 2010; Albrecht et al. 2012; Santerne et al. 2014), polar (e.g., Our stars were selected following four main criteria: Addison et al. 2013), and even retrograde orbits (e.g., Bayliss – the star should have a known transiting planet; et al. 2010; Hébrard et al. 2011). – the planet host star should be very active (e.g., the published The contrast of active regions (i.e., stellar spots, plages, photometric transit light curves clearly show large stellar and faculae) present on the stellar surface during the transit of spots’ occultation anomalies); an exoplanet alter the high-precision photometric transit light – at least one RM measurement of the planet has been curve (both outside and inside transit) and lead to an inaccurate observed (preferably with the HARPS spectrograph as we estimate of the planetary parameters (e.g., Czesla et al. 2009; use HARPS observations in our study) to ensure the feasi- Sanchis-Ojeda & Winn 2011; Sanchis-Ojeda et al. 2013; Oshagh bility of detecting the RM signal; et al. 2013b, 2015a,b; Barros et al. 2013; Ioannidis et al. 2016). – the transiting planet should have several observable tran- Since the physics and geometry behind the transit light curve sits from the southern hemisphere during the period from and the RM measurement are the same, the RM observations are March 2017 to October 2017 (this criterion is to ensure expected to be affected by the presence of stellar active regions observability). in a similar way. Furthermore, the inhibition of the convective Our criteria resulted in the selection of six transiting exo- blueshift inside the active regions can induce extra RV noise on planets, namely, WASP-6b, WASP-19b, WASP-41b, WASP-52b, the RM observations, which is absent in the photometric transit CoRoT-2b, and Qatar-2b. Their published transit light curves observation. have indicated the presence of stellar spots on their surfaces Boldt et al. (in prep.) demonstrate, using simulations, that with a filling factor of 6, 8, 3, 15, 16, and 4%, respectively unocculted stellar active regions during the planetary transit (Nutzman et al. 2011; Tregloan-Reed et al. 2015; Sedaghati et al. influence the RV slope of out-of-transit in RM observations (due 2015; Southworth et al. 2016; Mancini et al. 2017; Dai et al. to a combination of contrast and stellar rotation) and leads to 2017). Most of them also have at least one observed RM with an inaccurate estimate of the planetary radius; however, their HARPS (Bouchy et al. 2008; Gillon et al. 2009; Hellier et al. influence on the estimated λ was negligible. On the other hand, 2011; Neveu-VanMalle et al. 2016), except for Qatar-2, which has Oshagh et al.(2016) demonstrated that the occulted stellar active one observed RM with the HARPS-N spectrograph (Esposito regions can lead to quite significant inaccuracies in the estima- et al.