Measurements of the obliquities in exoplanet systems have revealed some remarkable architectures some of which are very different from the Solar System. Nearly 200 obliquity measurements have been obtained through observations of the RossiterMcLaughlin RM effect. Here we report on observations of 19 planetary systems that led to 17 clear detections of the RM effect and 2 less secure detections. After adding the new measurements to the tally we used the entire collection of RM measurements to investigate four issues that have arisen in the literature. i Does the obliquity distribution show a peak at approximately 90 We find tentative evidence that such a peak does exist when restricting attention to the sample of subSaturn planets and hot Jupiters orbiting F stars. ii Are high obliquities associated with high eccentricities We find the association to be weaker than previously reported and that a stronger association exists between obliquity and orbital separation possibly due to tidal obliquity damping at small separations. iii How low are the lowest known obliquities Among hot Jupiters around cool stars we find the dispersion to be 1.4  0.7 smaller than the 6 obliquity of the Sun which serves as additional evidence for tidal damping. iv What are the obliquities of stars with compact and flat systems of multiple planets We find that they generally have obliquities lower than 10 with several remarkable exceptions possibly caused by wideorbiting stellar or planetary companions.
['planets and satellites: dynamical evolution and stability', 'planets and satellites: formation', 'planet-star interactions', 'Astrophysics - Earth and Planetary Astrophysics']
Obliquities of exoplanet host stars Nineteen new and updated measurements and trends in the sample of 205 measurements
{'19 New and Updated Measurements, and Trends in the Sample of 205 Measurements': "E. Knudstrup 1 , 2 , S. H. Albrecht 2 , 3 , J. N. Winn 3 , D. Gandolfi 4 , J. J. Zanazzi 5 ⋆ , C. M. Persson 6 , M. Fridlund 7 , 6 , M. L. Marcussen 2 , A. Chontos 3 , 8 , M. A. F. Keniger 9 , N. L. Eisner 10 , A. Bieryla 11 , H. Isaacson 12 , A. W. Howard 13 , L. A. Hirsch 14 , F. Murgas 15 , 16 , N. Narita 17 , 18 , 15 , E. Palle 15 , 16 , Y. Kawai 19 , and D. Baker 20  \n- 1 Department of Space, Earth and Environment, Chalmers University of Technology, 412 93, Gothenburg, Sweden e-mail: emil.knudstrup@chalmers.se\n- 2 Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark\n- 3 Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA\n- 4 Dipartimento di Fisica, Universita degli Studi di Torino, via Pietro Giuria 1, I-10125, Torino, Italy\n- 5 Astronomy Department, Theoretical Astrophysics Center, and Center for Integrative Planetary Science, University of California, Berkeley, Berkeley, CA 94720, USA\n- 6 Chalmers University of Technology, Department of Space, Earth and Environment, Onsala Space Observatory, SE-439 92 Onsala, Sweden\n- 7 Leiden Observatory, University of Leiden, PO Box 9513, 2300 RA Leiden, The Netherlands\n- 8 Institute for Astronomy, University of Hawai'i, 2680 Woodlawn Drive, Honolulu, HI 96822, USA\n- 9 Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK\n- 10 Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010, USA\n- 11 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA\n- 12 501 Campbell Hall, University of California at Berkeley, Berkeley, CA 94720, USA\n- 13 Department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA\n- 14 Department of Chemical and Physical Sciences, University of Toronto Mississauga, Mississauga, ON, Canada\n- 15 Instituto de Astrofísica de Canarias (IAC), 38205 La Laguna, Tenerife, Spain\n- 16 Departamento de Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain\n- 17 Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan\n- 18 Astrobiology Center, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan\n- 19 Department of Multi-Disciplinary Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan\n- 20 Physics Department, Austin College, Sherman, TX 75090, USA  \nReceived May 8, 2024; accepted August 16, 2024", 'ABSTRACT': 'Measurements of the obliquities in exoplanet systems have revealed some remarkable architectures, some of which are very di ff erent from the Solar System. Nearly 200 obliquity measurements have been obtained through observations of the Rossiter-McLaughlin (RM) e ff ect. Here we report on observations of 19 planetary systems that led to 17 clear detections of the RM e ff ect and 2 less secure detections. After adding the new measurements to the tally, we use the entire collection of RM measurements to investigate four issues that have arisen in the literature. i) Does the obliquity distribution show a peak at approximately 90 · ? We find tentative evidence that such a peak does exist when restricting attention to the sample of sub-Saturn planets and hot Jupiters orbiting F stars. ii) Are high obliquities associated with high eccentricities? We find the association to be weaker than previously reported, and that a stronger association exists between obliquity and orbital separation, possibly due to tidal obliquity damping at small separations. iii) How low are the lowest known obliquities? Among hot Jupiters around cool stars, we find the dispersion to be 1 . 4 ± 0 . 7 · , smaller than the 6 · obliquity of the Sun, which serves as additional evidence for tidal damping. iv) What are the obliquities of stars with compact and flat systems of multiple planets? We find that they generally have obliquities lower than 10 · , with several remarkable exceptions possibly caused by wide-orbiting stellar or planetary companions.  \nKey words. Planet-star interactions - Planets and satellites: dynamical evolution and stability - Planets and satellites: formation', '1. Introduction': "Hot Jupiters, planets on extremely eccentric orbits, planets orbiting binary stars, planets orbiting white dwarfs and neutron stars -these and other types of exotic planetary systems have been  \ndiscovered over the last 30 years. A diversity of architectures has also been revealed by measurements of the host star's obliquity, i.e., the angle between the star's spin axis and a planet's orbital axis, or at least its sky projection. Thanks to many research groups pursuing these measurements, usually by means of the Rossiter-McLaughlin (RM) e ff ect, we now have a sample of nearly 200 measurements including prograde, polar, and retrograde orbits (see, e.g., Albrecht et al. 2022, for a review).  \nHere we present observations of 19 transiting exoplanets aiming to detect the RM e ff ect and measure the sky-projected stellar obliquities. Together with the sample drawn from the literature, we also wanted to investigate four issues regarding the distribution of obliquities.  \n- i) For most systems in the sample, the projected obliquity ( λ ) is known but not the obliquity itself ( ψ ) because of missing information about the inclination of the stellar rotation axis with respect to the line of sight ( i ⋆ ). Using the subset of 57 systems for which ψ has been measured, Albrecht et al. (2021) found evidence for a preponderance of perpendicular planets, i.e., a peak in the obliquity distribution near 90 · . Following up on this result, Siegel et al. (2023) and Dong & Foreman-Mackey (2023) applied more advanced statistical techniques that allowed the entire sample to be used, and did not find evidence for a peak. With an enlarged sample, we wanted to revisit this issue and see whether there is a particular category of planets for which a peak does exist.\n- ii) Three theories have been proposed to explain the existence of hot Jupiters: in-situ formation, disk-driven migration, and high-eccentricity or tidally-driven migration (see Dawson & Johnson 2018, for a review). The latter process would not only increase the orbit's eccentricity but might also raise its inclination relative to the star's equatorial plane. One might therefore expect high eccentricities and obliquities to be statistically associated, and evidence for such an association has been reported (Rice et al. 2022a). After noticing some errors in this earlier study, we decided to revisit the issue with an enlarged sample.\n- iii) Stars with T e ff ≲ 6 , 250 K that host hot Jupiters tend to have especially low obliquities. Indeed, the most precise such measurements ( σλ < 2 · ) show an obliquity dispersion less than one degree - smaller than 6.2 · obliquity of the Sun relative to the Solar System's invariable plane. Highlighting this result, Albrecht et al. (2022) interpreted the very low obliquities of hot Jupiters around cool stars as evidence for tidal obliquity damping. If this is true, then by continuing to perform precise measurements, we might gain a better understanding of the evolution of hot Jupiters and rates of tidal dissipation. With the enlarged sample, we wanted to see if this trend still holds and include the numerous but less precise obliquity measurements in the determination of the dispersion.\n- iv) The first five measurements of λ for stars with compact multiple-transiting planetary systems ('multis') were all consistent with low obliquities (Albrecht et al. 2013), following the blueprint of the Solar System. Since then, several exceptions have been discovered (Hjorth et al. 2021; Huber et al. 2013; Chaplin et al. 2013). The misalignments in these cases are suspected of being caused by the e ff ects of an outer massive companion. We wanted to assess the obliquity distribution of the multis, and see whether all of the misaligned multis have outer companions.  \nThe paper is structured as follows. Section 2 describes the general characteristics of the new data. Section 3 presents our  \ngeneral approach to analyzing the data, and Section 4 discusses the particulars for each of the 19 observed systems. Section 5 discusses the issues and questions posed above, and Section 6 summarizes our conclusions.", '2. Observations': 'To measure the projected obliquity of each host star, we performed high-resolution optical spectroscopy over a time range (typically 4-8 hours) spanning a planetary transit and tried to detect the RM e ff ect. Table 1 gives the names of the 19 systems, the names of the telescopes and instruments employed, and some key observational characteristics. The telescopes and instruments are further described in Table 2. Tables 3 and 4 give some of the basic parameters of the host stars and their planets. Table B.5 provides all of the apparent radial velocity (RV) measurements. These spectroscopic data were supplemented by the best available transit light curves, as described below.', '2.1. Spectroscopy': "Seven systems were observed with the Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO; Pepe et al. 2021) that can be fed by any of the four Very Large Telescopes at Paranal Observatory, Chile. Another seven systems were observed with the FIber-fed Echelle Spectrograph (FIES; Frandsen & Lindberg 1999; Telting et al. 2014) mounted on the Nordic Optical Telescope (NOT; Djupvik & Andersen 2010) located on the Roque de los Muchachos, La Palma, Spain. For five systems, we used the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N; Cosentino et al. 2014) mounted on the Telescopio Nazionale Galileo (TNG), also on the Roque de los Muchachos. For one system, we used the High Resolution Echelle Spectrometer (HIRES; Vogt et al. 1994) mounted on the 10-m Keck-1 telescope on Mauna Kea, Hawai'i, USA.  \nThe data from ESPRESSO and HARPS-N were reduced and radial velocities (RVs) were obtained using the standard data reduction software (DRS; Lovis & Pepe 2007; Pepe et al. 2021; Dumusque et al. 2021). The FIES data for all but two systems were reduced using FIEStool 1 , and RVs were extracted using FIESpipe 2 following an approach similar to that described by Zechmeister et al. (2018). The FIES observations of WASP-136 and WASP-186 were reduced and RVs extracted following the procedure described by Gandolfi et al. (2015). To account for any instrumental RV drift of the FIES spectrograph, we observed the spectrum of a Thorium-Argon (ThAr) lamp in between each science exposure. The HIRES spectra were obtained with the iodine cell in the light path, and the RV extraction followed the standard HIRES forward modelling pipeline of the California Planet Search (CPS; Howard et al. 2010) using the methodology of Butler et al. (1996).  \nFor some systems, in addition to the apparent RVs we modeled the shape of the Cross Correlation Function (CCF) of the absorption lines. For the ESPRESSO and HARPS-N data, we used the CCFs that are automatically produced by the standard reduction pipelines. For the FIES data, we derived CCFs using a template spectrum based on the most suitable model atmosphere from the library of Castelli & Kurucz (2003), as implemented in FIESpipe .  \nTable 1. Systems observed and spectroscopic observations.  \nIn the columns we list the Johnson V magnitudes, the instrument used, the program IDs, the dates of the observation nights, the number of exposures on a given night ( N ), the exposure duration ( t ), sampling time ( ∆ t ), and the total transit duration ( T 41).  \nTable 2. Telescopes and spectrographs used in this study.  \nR is the resolution and σ PSF the width of the point spread function in velocity units.", '2.2. Photometry': 'Light curves of all 19 of the systems are available from the databases of either the K2 mission (Howell et al. 2014) or the Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015). For the K2 data, we used the light curves from the EVEREST pipeline (Luger et al. 2016, 2018). The TESS data were downloaded and extracted using the lightkurve package (Lightkurve Collaboration et al. 2018), and the RegressionCorrector was applied to reduce the systematic variations due to scattered light. Our parametric model for each system was based on jointly fitting the photometric and spectroscopic data. In some cases, we supplemented the space-based photometry with ground-based photometry to refine the transit ephemerides.', '3.1. Photometric data': "The photometric data provide tight constraints on the orbital period ( P ), the time of conjunction of a transit chosen to be the  \nreference epoch, ( T 0), the planet-to-star radius ratio ( R p / R ⋆ ), the cosine of the orbital inclination (cos i ), and the ratio of the semimajor axis and stellar radius ( a / R ⋆ ). Uniform prior probability distributions were adopted for these parameters.  \nTo model the stellar limb darkening function, we used a standard quadratic law and allowed the sum of the two coe ffi cients (the center-to-limb intensity ratio) to be an adjustable parameter, while holding the di ff erence fixed at the value obtained from the tables of Claret et al. (2013); Claret (2018) as appropriate for the photometric bandpass and the star's e ff ective temperature ( T e ff ), surface gravity (log g ), and metallicity ([Fe / H]). See Table 5 and Table 7 for details. The sum of coe ffi cients was subjected to a Gaussian prior with a width of 0 . 1 and a central value from the aforementioned tables.  \nIn cases for which a nonzero orbital eccentricity has been reported (HD 118203 b, K2-261 b, K2-287 b, WASP-148 b, and WASP-186 b), we applied Gaussian priors on e and the argument of periastron ( ω ) using the values given in Table 3. For the multis, we applied Gaussian priors on the parameters of all planets apart from the transiting planet that was the target of our observations. The main e ff ect of these planets on our model is to change the slope of the RV time series on the night of the transit that was observed spectroscopically. When two planets are transiting simultaneously, the shape of the transit light curve is also a ff ected. Simultaneous transits occasionally occurred (for TOI-1130 and TOI-451A) within the long intervals spanned by the photometric data, but did not occur during any of our 19 spectroscopic observations.  \nWe modeled the light curves using the batman package (Kreidberg 2015), which is based on the equations of Mandel & Agol (2002). In addition to accounting for the transit and white noise in the measurements, we included a Gaussian process (GP) using celerite (Foreman-Mackey et al. 2017). We generally employed a Matérn-3 / 2 kernel, which is characterised by two hy-  \nTable 3. Key literature orbital and planetary parameters.  \nIn the period column, the number in parentheses is the uncertainty in the last digit of the period.  \nSources are given in Table 4 for the corresponding host.  \n(a) With 98% confidence.  \nperparameters: the amplitude ( A ) and the timescale ( τ ). For the time series with 30-minute sampling, evenly spaced model light curves were created and integrated over to mimic a cadence of 2 min.", '3.2. Stellar rotation': 'The analysis of the RM e ff ect is aided by knowledge of the stellar rotation velocity or period. We tried to use the available spectroscopic and photometric data to measure these quantities independently of the transit data, as describe below.', '3.2.1. Projected rotation velocity': 'The amplitude of the RM e ff ect scales with v sin i ⋆ , the product of the stellar rotation speed and the sine of the line-of-sight inclination of the stellar spin axis. Therefore, knowing v sin i ⋆ helps to plan observations and analyze the results. When the RM effect is detected with a high signal-to-noise ratio (SNR), v sin i ⋆ can be determined precisely from the RM data. Indeed, for stars with a high obliquity, the e ff ects of di ff erential rotation might be detectable (Gaudi & Winn 2007; Cegla et al. 2016a).  \nIn more typical cases in which the RM e ff ect is detected with a modest SNR, it is useful to determine v sin i ⋆ by modeling the spectral line broadening outside of transits, and thereby gain the option of applying a prior constraint on v sin i ⋆ when modeling the RM e ff ect. For this purpose, we used the broadening function (BF) method of Rucinski (1999) as implemented in FIESpipe , which is based on an appropriate model atmosphere from Castelli & Kurucz (2003). The BFs were created from high SNR orders from out-of-transit spectra, where we visually inspected the result to ensure that the BF created from a given order was well-behaved. This typically meant selecting  \ncentral orders typically with wavelengths ∼ 5000-6000 Å. Orders containing telluric lines were excluded. For a given epoch we stacked the BFs from the individual orders to create the highest possible SNR BF for that epoch. We applied this method not only to the FIES spectra but also the ESPRESSO and HARPSN spectra. To extract v sin i ⋆ we fitted the theoretical rotational BF from Kaluzny et al. (2006, Eq. (2)). The model includes the rotational broadening parameters along with a parameter representing non-rotational broadening (in most cases dominated by the finite instrumental resolution).  \nTable 6 gives the results. The tabulated v sin i ⋆ is the median of the best-fit results of analyzing spectra from di ff erent epochs, and the tabulated uncertainty is the standard deviation between those epochs, which in most cases are unrealistically small. This is especially true for observations where we only have few outof-transit observations. For systems where we only had two outof-transit spectra, we used the first ingress as a third spectrum. Our measurements are generally in agreement with previously published results. The most significant deviations are seen for the most slowly-rotating stars, where the details of extracting line profiles and modeling non-rotational broadening are most important. When modeling the RM e ff ect, we imposed Gaussian priors based on our v sin i ⋆ determinations, with a minimum width of σ = 0 . 5 km s -1 to account for possible systematic errors.  \nTo test if our choice of the width of the uncertainty interval in the v sin i ⋆ prior does significantly a ff ect our results for the projected obliquity, we selected four systems in our sample for which we expect the result of lambda to be most sensitive to v sin i ⋆ (low impact parameter, incomplete transit coverage) and performed the same analysis as described below. However, now with an uncertainty in v sin i ⋆ of 1.0 km s -1 . In all four cases we find that the highest probability value for λ changes  \nTable 4. Key literature stellar parameters.  \nThe sources refer to the following references: 1 Pepper et al. (2020), 2 Chontos et al. (2024), 3 Johnson et al. (2018), 4 Jordán et al. (2019), 5 Pepper et al. (2013), 6 Eastman et al. (2016), 7 Winters et al. (2022), 8 Newton et al. (2021), 9 Eisner et al. (2020), 10 Brahm et al. (2020), 11 Huang et al. (2020), 12 Gillon et al. (2011), 13 Hébrard et al. (2013), 14 Lam et al. (2017), 15 Hébrard et al. (2020), 16 Hellier et al. (2019a), 17 Hellier et al. (2019a), 18 Schanche et al. (2020), 19 Crouzet et al. (2020).  \nonly by a small fraction of its uncertainty interval, relative to the 0.5 km s -1 case.', '3.2.2. Photometric rotation period': "We attempted to measure the rotation period ( P rot) by seeking periodicities in the K2 and TESS light curves. The combination of the rotation period, v sin i ⋆ , and R ⋆ can be used to constrain the stellar inclination (Masuda & Winn 2020), which in turn can be used in combination with i orb and λ to determine the 'true' or three-dimensional obliquity.  \nWe used the method based on the auto-correlation function (ACF) described by McQuillan et al. (2014). After removing the data a ff ected by transits, we searched for peaks in the ACF as a function of lag. When multiple peaks were seen, we fitted a linear function between peak number and lag, the slope of which is an estimate of the rotation period (e.g., Hjorth et al. 2021).  \nTable 6 gives the results. Empty entries mark the cases for which no period could be determined. For two systems, TOI813 and TOI-892, peaks were seen in the ACF (corresponding to periods of 10 . 1 ± 0 . 6 d and 2 . 76 ± 0 . 11 d, respectively), but because the statistical significance was weak, the results do not appear in Table 6. For WASP-50, a series of strong peaks seemed to imply a period of 7 . 2 ± 1 . 8 d, which contrasts with the value of 16 . 3 ± 0 . 5 d reported by Gillon et al. (2011). They also noted that the true rotation period could be 32.6 days if there are similar spot patterns on opposite stellar hemispheres. We folded the light curve using trial periods of 16 . 3 d or 32 . 6 d but neither choice seemed compelling. Evidently, the variability of WASP-50 is complex, which is why Table 6 contains no entry for WASP-50.  \nWe note that the estimated rotation periods of WASP-136 and WASP-186, 5 . 4 ± 0 . 5 d and 4 . 8 ± 1 . 9 d, are both close to the orbital period of the transiting planets. These systems might have  \nbeen driven into spin-orbit synchronization by tidal interactions between the planet and star (e.g., Albrecht et al. 2012b; Brown 2014; Penev et al. 2018).", '3.3. The Rossiter-McLaughlin effect': "The RM e ff ect is the distortion of a star's absorption lines that arises from the combination of stellar rotation and a transiting planet. During a transit, the planet prevents the light from a portion of the photosphere from contributing to the disk-integrated spectral line. The missing spectral component causes the velocity profile of the spectral lines to have a slightly reduced flux at velocities centered on the rotational Doppler shift of the star at the 'subplanetary' point, i.e., the point on the the star directly behind the center of the planet. Observations of the RM e ff ect can therefore be used to determine the time series of subplanetary velocities, or equivalently, the trajectory of the transiting planet relative to the stellar rotation axis.  \nNeglecting di ff erential rotation, the RM e ff ect is mainly a function of λ , v sin i ⋆ , R p / R ⋆ , and the transit impact parameter b ≡ d / R ⋆ cos i , where d is the star-planet distance at the time of conjunction. For precise modeling it is also necessary to take into account limb darkening and non-rotational spectral line broadening. To model non-rotational broadening due to the star, we used the Gray (2005) model with parameters for the velocity spread generated by microturbulence ( ξ⋆ ) and macroturbulence ( ζ⋆ ). We included a separate term for instrumental broadening, another source of non-rotational broadening. We neglected the e ff ects of di ff erential rotation and the convective blueshift, which for solar-type stars are not expected to produce detectable e ff ects given the quality of our data (and indeed our modeling did not uncover any evidence for those e ff ects, with one possible exception; see Section 4.12.1). We did not model the e ff ects of  \nFig. 1. The Rossiter-McLaughlin e ff ect displayed in four ways, for the illustrative case of WASP-50. Top right: Schematic of the planet transiting the stellar disk. The wedge represents the best-fit λ and its uncertainty. The color scheme is based on the star's rotational Doppler shift; the same color scheme is used in the top left and lower left panels to convey the location of the planet on the star. Top left: Spectral lines. The thick black curve is the average out-of-transit cross-correlation function (CCF) and the superimposed white curve is a CCF observed near mid-transit. Middle left: Spectral line distortions. Shown are the results of subtracting five of the in-transit CCFs from the mean out-of-transit CCF. Solid curves show the data. Dashed curves are the best-fit Lorentzian profiles, with vertical lines marking their centers. The color coding is done according to the apparent velocity anomaly going from redshifted (in red) to blueshifted (in blue). Bottom left: The Doppler shadow. Each row shows the spectral line distortion at a particular time. The white contours are amplitude levels of the distortion. Each small circle marks the best-fit subplanetary velocity (the Doppler shift of the portion of the star directly behind the planet's center), with colors that convey the corresponding anomalous radial velocity shown in the lower right panel. Dark vertical lines mark the measured v sin i ⋆ , and the diagonal white line is the best-fit trajectory. Bottom right: Anomalous radial velocities. The dark shaded area is the range of times when the planet's silhouette is completely contained within the stellar disk.  \n<!-- image -->  \nstarspots, flares, or pulsations, since there was no evidence for such phenomena in our data.  \nThe RM e ff ect has been analyzed as a distortion of individual spectral lines (e.g., Albrecht et al. 2007), as distortion of the spectral cross-correlation function (e.g., Cegla et al. 2016a), or as an 'anomalous RV' obtained by an ordinary radial-velocity extraction code to a spectrum a ff ected by the RM e ff ect (e.g., Queloz et al. 2000). We used a combination of these approaches, depending on the system, as described below.  \nWe fitted the spectral CCFs with a parameterized model that incorporates the RM e ff ect, using a code first described by Knudstrup & Albrecht (2022) and Knudstrup et al. (2023a). In this model, the stellar disk is discretized (with a radius of 100 pixels) and a spectrum is assigned to each pixel according to the model parameters for limb darkening, stellar rotation, and non-  \nrotational broadening. A disk-integrated spectral line is created by summing the spectra of all pixels that are not concealed by the planet, and the disk-integrated line profile is compared to the observed spectral CCF.  \nAn alternative is to fit the time series of subplanetary velocities derived from the residuals between the out-of-transit CCF and the in-transit CCFs. The subplanetary velocity is calculated as  \nvp ( t ) = v sin i ⋆ x ( t ) , (1)  \nwhere x ( t ) is the distance from the projected rotation axis divided by the stellar radius. Once the orbital parameters are chosen, x ( t ) can be calculated. It is also instructive to calculate x at the times  \nTable 5. Limb-darkening coe ffi cients and velocity fields.  \nThe second, third, and fourth columns show the q 1 and q 2 quadratic limb darkening pairs we used as priors in the TESS, K2, and spectrograph bands. For the latter we simply assumed the V -filter. We obtained these using values for T e ff , log g , and [Fe / H] in tabulated in Table 4 and we queried the limbdarkening tables for the quadratic law by Claret et al. (2013, for K2 and V ) and Claret (2018, for TESS). The macro- ( ζ⋆ ) and micro-turbulence ( ξ⋆ ) priors were obtained from the relations in Doyle et al. (2014) and Bruntt et al. (2010), respectively, again using the values from Table 4, with the exception of WASP-172 and LTT 1445A (see Section 3).  \nof ingress ( x 1) and egress ( x 2) (Albrecht et al. 2011):  \nx 1 = √ 1 -b 2 cos λ -b sin λ, x 2 = √ 1 -b 2 cos λ + b sin λ.  \n(2)  \nTherefore, the very nearly linear trajectory of the planet can be calculated as a function of time for given values of t , T c, the transit duration, b , and λ .  \nThe simplest and most common approach is to fit a time series of anomalous RVs. In the model, the anomalous RV is computed as a function of vp by taking into account the size of the planet's shadow, the limb darkening function, and the response of the radial-velocity extraction code to the RM distortion; for this purpose we used the code described by Hirano et al. (2011).  \nThe di ff erent ways to sense the RM e ff ect are illustrated in Fig. 1 for the case of WASP-50. For all 19 systems in our sample, we fitted the time series of anomalous RVs (the 'RV-RM' method). For some systems, we also analyzed the distortions of the CCF, or fitted the time series of subplanetary velocities. Multiple methods were used to test for consistency in the results, and because for some systems it seemed plausible that the more sophisticated models would give more information about λ . In general, we expect line-profile modeling to be advantageous when  \nTable 6. Stellar rotation parameters.  \nMeasurements of v sin i ⋆ and P rot from the literature (see Table 4 for citations) and those derived in this work. The uncertainties reported for v sin i ⋆ generally refer to statistical uncertainties or the scatter in results of analyzing di ff erent spectra. Because they do not include systematic uncertainties, we imposed a minimum uncertainty of 0 . 5 km s -1 in our analysis.", '( R p / R ⋆ ) 2 v sin i ⋆ is larger than the non-rotational broadening (Albrecht et al. 2022).': 'We determined the best-fit parameters and their uncertainties via Markov Chain Monte Carlo (MCMC) sampling of the posterior probability distribution in parameter space. We used the emcee package (Foreman-Mackey et al. 2013). To ensure convergence we both inspected the chains visually and calculated the rank normalized R-hat 3 . As mentioned above, we always fitted the the transit light curve alongside the spectroscopic data.  \nWhen fitting the RVs and the line distortion directly, we used Gaussian priors for the macroturbulence and microturbulence parameters, and we adopted a fixed value for the instrumental broadening parameter of each spectrograph (see Table 2). The mean values for our Gaussian priors were obtained from the relationships presented by Doyle et al. (2014) for stars with T e ff ∈ [5200 K , 6400 K] and log g ∈ [4 . 0 , 4 . 6], and Bruntt et al. (2010) for stars with T e ff ∈ [5000 K , 6500 K] and log g > 4 . 0. The width of the prior in all cases was 1 km s -1 .  \nWe note that the properties of WASP-136, TOI-813, and HD 118203 are on the borderline of the range of properties for which the aforementioned relationships involving turbulent broadening are applicable; we applied them anyways. On the other hand, the properties of WASP-172 and LTT 1445A are far outside the applicable ranges. For WASP-172 ( T e ff = 6900 ± 150 K), we expect a high macroturbulent velocity. We chose 6 . 0 km s -1 and increased the width of the prior to 4 . 0 km s -1 to reflect our uncertainty. For the microturbulent velocity, we adopted ξ⋆ = 1 . 5 ± 1 . 5 km s -1 . In the other end of the spectrum,  \nwe expect the turbulent velocities of the M dwarf LTT 1445A to be smaller. We chose ζ = 1 . 0 ± 1 . 0 km s -1 and ξ⋆ = 1 . 0 ± 1 . 0 km s -1 . All of the priors are summarized in Table 5.  \nFor each system, we performed two RV-RM analyses, one in which we adopted a uniform prior on v sin i ⋆ , and one in which we applied a Gaussian prior on v sin i ⋆ based on the observed line broadening of the out-of-transit spectrum (Table 6) and a width of 0 . 5 km s -1 . When analyzing the line profiles or Doppler shadow, we always applied the Gaussian prior on v sin i ⋆ .  \nWhen fitting the time series of subplanetary velocities, we did not simultaneously fit the light curve. Instead, we applied Gaussian priors on P , T 0, R p / R ⋆ , i , and a / R ⋆ based on an external fit to the light curve, and adopted uniform priors on λ and v sin i .  \n⋆  \nIn some cases, to boost the SNR, we stacked the CCFs from di ff erent exposures in a way similar to Johnson et al. (2014), but with a code based more directly on the work of Hjorth et al. (2019). In these cases, we extracted the Doppler shadow as shown in Fig. 1, and then transformed the velocities by correcting for the slope in Eq. (2). The white line / slope in Fig. 1 is in this way used to bring each pixel directly beneath it in the Doppler shadow to be located at 0 km s -1 . We then collapse the slope corrected shadow by summing along the ordinate (time stamps) for a given geometry ( λ , v sin i ⋆ , b ). When the model parameters match reality, we expect to see a sharp peak at 0 km s -1 in the summed spectrum, and we found the height of the peak by fitting a Gaussian to it. Hjorth et al. (2019) did this for a dense grid in the parameter space of b , v sin i ⋆ , and λ . For our purposes, since b was always tightly constrained by the light curve, we neglected the uncertainty in b and simply adopted the best-fit value from the RV-RM method. We then searched the 2-d parameter space of v sin i ⋆ and λ for the strongest peak; in practice this was done by fitting a 2-d Gaussian function. We did not run an MCMC when stacking the CCFs, nor did we fit the light curve simultaneously. The quoted uncertainties are taken as the widths ( σ s) of the fitted 2D Gaussian and as such do not reflect the actual precision. The resulting values therefore only serves as an indication to whether the results in ( λ , v sin i ⋆ ) from the RV-RM runs are in (qualitative) agreement with that seen in the CCFs.', '4. Analysis of the individual systems': 'Below, we describe the particulars of each system, our analysis, and the results. The details for each system are also summarized in Table 1, Table 4, and Table 3.', '4.1. HD 118203': 'HD 118203 b was discovered by da Silva et al. (2006) through RV monitoring, and was later found to be a transiting planet by Pepper et al. (2020) using TESS data. TESS data are available for HD 118203 from Sectors 15, 16, 22, and 49 with a 2-minute cadence. The planet is a hot Jupiter on an eccentric orbit around a slightly evolved G-type star (Pepper et al. 2020). The TESS data centered on the mid-transit time are shown in Fig. A.1.  \nWe observed two transits of HD 118203 b with FIES. The first transit began on March 24, 2020 at around UT 21:00 and ended at 04:15. The observations were interrupted twice due to high humidity. The exposure time was 900 s and the time between exposures was 1171 s, with the time in between spent on ThAr calibration. Two exposures on this night were heavily a ff ected by weather and the resulting RVs deviate significantly from the rest (see the gray points in Fig. 2). These two data points were omitted from the analysis.  \nSince the first transit had been interrupted by the weather, we organized observations of a second transit on the night starting on June 18, 2022. Unfortunately, the first part of the transit was missed because the telescope was being used to observe a target-of-opportunity. Our transit observations started at around UT 22:30 and continued until around 03:00. As the FIES spectrograph was refurbished in the summer of 2021 resulting in an improved throughput, we opted for a shorter exposure time of 600 s for the second transit, which resulted in a sampling time of around 768 s. We note that the refurbishment introduced an o ff set in RV between observations. As such we allowed the systemic velocity parameter to di ff er between the two nights: γ FIES applies to the observations taken before summer 2021, and γ FIES + applies to subsequent observations.  \nAs mentioned above (Section 3) we performed two MCMC analyses, one with a uniform prior on v sin i ⋆ and one with a Gaussian prior. For this system, the results di ff ered substantially. The uniform prior in v sin i ⋆ yielded λ = -74 + 6 -9 · and the Gaussian prior yielded -23 + 25 -38 · . The data obtained from the first transit observations are obviously very noisy. This paired with the second transit observations only covering the second half of the transit meant that the amplitude of the signal was allowed to vary significantly, which made v sin i ⋆ wander o ff to large values when applying applying a uniform prior. In this case v sin i ⋆ ended up being 26 + 13 -23 km s -1 , which is both poorly constrained and at face value in disagreement with the value we derived and the value found in the literature. Therefore, our best estimate is -23 + 25 -38 · based on the analysis in which v sin i ⋆ was subjected to a Gaussian prior.', '4.2. HD 148193': "HD148193 is also known as TOI-1836. The planet HD 148193 b (TOI-1836.01) was confirmed and characterized by Chontos et al. (2024). It is a warm Saturn around a subgiant star in a system with a companion star. A second transiting planet candidate was announced later (TOI-1836.02), which is potentially a super-Earth-sized planet on a 1.7-day orbit. However, the signalto-noise ratio of the observed transits is modest and the planet remains unconfirmed.  \nTESS data are available with 30-min cadence from Sectors 16 and 25; with 2-min cadence in Sectors 23, 24, 49, 50, 51, and 52; and with 20-sec cadence in Sector 56. The light curve is shown in Fig. A.1.  \nWe used HARPS-N at the TNG to observe a transit of HD 148193 b. We started observations around UT 21:30 and continued until UT 05:45 on the night starting on June 15, 2021. The exposure time was 1140 s. After accounting for overheads, the sampling time was 1168 s. The transit observations are shown in Fig. 3.  \nThe combined fit to the light curve and RVs gave λ = 29 + 12 -13 · . However, we obtained only one pre-ingress RV and two postegress RVs. With such sparse out-of-transit data, any mismatch between the out-of-transit RV slope observed on that night and the calculated slope based on the planet's orbital parameters can lead to systematic errors. This problem, combined with the reasons discussed in Section 3.3 and the large projected rotation speed of this system ( ∼ 6 km s -1 ), led us to analyze the line distortions and Doppler shadow. We found λ = 14 + 7 -8 · , with results shown in Fig. 4. We adopted the value from the Doppler-shadow analysis as our best estimate for λ .  \n<!-- image -->  \nFig. 2. The RM e ff ect observed in HD 118203. Left: Time series and residuals of the FIES + RVs (orange) and FIES RVs (green). The black curve is the best-fit model and the gray shading shows the 1 and 2 σ confidence intervals created by randomly drawing values from the MCMC chains for b , v sin i ⋆ , R p / R ⋆ , and λ . The dashed curve is the best-fit model assuming λ = 0 · . Gray circles are the data that were obtained during poor atmospheric conditions and were excluded from the analysis. Right: Correlations between λ , v sin i ⋆ , and b . The red lines illustrate the prior that was applied (either uniform or Gaussian).  \n<!-- image -->  \n<!-- image -->  \nFig. 3. The RV-RM e ff ect for HD 148193. Same as for Fig. 2, but with RVs from HARPS-N.  \n<!-- image -->", '4.3. K2-261': 'K2-261 was observed in Campaign 14 of the K2 mission with 30-minute cadence. The planet was confirmed and characterized independently by Johnson et al. (2018), who named the planet K2-261 b, and Brahm et al. (2019), who named it K2-161 b. Here, we adopt the name K2-261 b, which is also the name chosen by Ikwut-Ukwa et al. (2020). The planet is a warm Saturnsized planet on an eccentric orbit around a slightly evolved Gtype star (Brahm et al. 2019). K2-261 was observed by TESS in Sectors 9, 35, 45, and 46. TESS data are available with 2-min cadence in Sectors 9, 45, and 46, and with 20-sec cadence in Sector 35. The light curve is shown in Fig. A.1.  \nFor K2-261 b we used HARPS-N at the TNG to observe a transit. The data are shown in Fig. 5. We carried out the observations on the night starting on March 24, 2022, from around UT  \n20:15 until 03:30. We used an exposure time of 600 s, resulting in a sampling time of around 625 s.  \nWe found a value for v sin i ⋆ for K2-261 that is somewhat smaller than the value of 3 . 70 ± 0 . 17 km s -1 reported by Johnson et al. (2018), both from our BF analysis (1 . 6 km s -1 ) and from our RM observations. The amplitude of the signal is not as strong as would be expected if v sin i ⋆ were really 3 . 70 ± 0 . 17 km s -1 . A uniform prior on v sin i ⋆ resulted in v sin i ⋆ = 1 . 0 + 0 . 5 -0 . 6 km s -1 and λ = 29 + 35 -56 · , and applying a Gaussian prior yielded a value for the projected obliquity of 32 + 48 -33 · . We adopt the latter as our best estimate for λ .', '4.4. K2-287': 'K2-287 was observed in Campaign 15 of the K2 mission. Jordán et al. (2019) discovered and characterized the system, finding the  \nFig. 4. Doppler shadow for HD 148193. Left: The observed Doppler shadow created by subtracting the in-transit CCFs from an average outof-transit CCF. Middle: Model Doppler shadow. Right: Residuals from subtracting the model from the observed shadow. The summed residuals are given in the adjacent panels. The color bar demarcate the contrast. Green horizontal lines denote points of first (solid), second (dashed), third (dashed), and fourth (solid) contact. The solid, orange, vertical lines denote ± v sin i ⋆ .  \n<!-- image -->  \n<!-- image -->  \nFig. 5. The RV-RM e ff ect for K2-261. Same as for Fig. 2, but with RVs from HARPS-N.  \n<!-- image -->  \nhost to be a G-dwarf and the planet to be a warm Saturn-sized planet on an eccentric orbit.  \nWe observed a transit of K2-287 b using ESPRESSO. Observations were carried out from UT 02:06 to 06:16 on the night July 6, 2021. The exposure time was 540 s, which including overhead comes out to a sampling time of roughly 600 s. The RVs derived from the spectra are shown in the left half of in Fig. 6.  \nGiven the 30-min time averaging of the K2 photometric time series, and the rather long orbital period of the planet, the K2 light curve is relatively poorly sampled. Furthermore, this star has evaded TESS observations to this point. As of the time of writing, there are no scheduled TESS observations that would cover this star.  \nThis situation, along with the lack of post-egress spectroscopoic observations, prompted us to search for a suitable archival ground-based light curve covering the egress. The one  \nwe found was obtained at the Adams Observatory at Austin College on the night starting June 30th, 2018 4 . More recently, the system was observed with the CHaracterising ExOPlanet Satellite (CHEOPS; Benz et al. 2021). Borsato et al. (2021) observed three separate transits with a 1-min cadence, and the data from all three provide good coverage of the entiure transit. In our analysis, we included the de-trended light curves presented by Borsato et al. (2021). The K2, CHEOPS, and ground-based light curves are shown in Fig. A.1.  \nThe two RV-RM analyses yielded similar results for λ , with 21 + 14 -11 · for the case of a uniform prior on v sin i ⋆ and 23 + 12 -13 · for the case of a Gaussian prior on v sin i ⋆ . We also investigated the in-transit subplanetary velocities, which are shown in Fig. 7 along with the best-fitting slope. In addition to the errors ob-  \n<!-- image -->  \nFig. 6. The RV-RM e ff ect for K2-287. Same as for Fig. 2, but with RVs from ESPRESSO.  \n<!-- image -->  \nen extracting the location of the subplanetary velocities, we add two di ff erent jitter terms in quadrature; one for the epochs taken when the planet is completely within the stellar disk (dark gray area in Fig. 7), and another when the planet is at the limb (light gray). This is because there generally (see, e.g, WASP-50 in Section 4.12) seems to be some additional scatter in the subplanetary velocities at the limb, which the errors from fitting a Gaussian to the residual peak do not fully capture.  \nFrom this we found values of v sin i ⋆ = 1 . 16 ± 0 . 19 km s -1 and λ = 37 + 8 -9 · , which are consistent with the values derived from the RV-RM e ff ect. Somewhat arbitrarily, we adopt λ = 23 + 12 -13 · as our best estimate; this derives from the RV-RM fit with a Gaussian prior on v sin i ⋆ . The posterior distributions are shown in Fig. 6.  \nFig. 7. Subplanetary velocities for K2-287. Subplanetary velocities for K2-287 with color-coding as described in Fig. 1. The dark shaded area corresponds to the time the planet is completely within the stellar disk and the light shaded area denote ingress and egress. The errorbars include a jitter term, where two separate terms have been added in quadrature; one for the epochs for the full transit and one for the epochs when the planet is partially outside of the limb of the star. The horizontal black lines in the top panel mark v sin i ⋆ from the fit.  \n<!-- image -->', '4.5. KELT-3': "Pepper et al. (2013) reported the discovery of the KELT-3 system, which consists of a hot Jupiter orbiting a late F star. TESS data with 2-min cadence are available from Sectors 21 and 48. The light curve is shown in Fig. A.1.  \nWe used FIES to observe a transit on the night starting April 8, 2022, with observations starting around UT 20:50 and continuing until UT 01:30 April 9, 2022. The exposure time was 900 s, which including overhead resulted in a cadence of 1080 s. The RM-RV e ff ect is shown on the left side of Fig. 8.  \nAs was the case for K2-261, the value for KELT-3's projected rotation velocity drawn from the literature (10 . 2 ± 0 . 5 km s -1 ; Pepper et al. 2013) is higher than what we obtained from the BF method (7 . 6 km s -1 ). The observed amplitude of the RM e ff ect suggests that the true v sin i ⋆ is even lower. Our RV-RM analysis gave 6 . 4 ± 0 . 4 km s -1 when using a uniform prior in v sin i ⋆ , and λ = -5 ± 4 · . When we applied a Gaussian prior on v sin i ⋆ , we obtained -6 ± 5 · .  \nGiven the discrepancy between the di ff erent measurements of v sin i ⋆ , we decided to look at the stacked CCFs resulting from our FIES spectra. In Fig. 9 we show the contours of the peak values obtained from stacking the CCFS on the 2D grid in the space of v sin i ⋆ and λ , compared to our posterior from the RVRMfit. From this analysis we found v sin i ⋆ = 5 . 96 ± 0 . 03 km s -1 and λ = 1 . 06 ± 0 . 16 · , implying that the system is well-aligned, and that v sin i ⋆ is closer to 6 km s -1 than the value of 10 . 2 ± 0 . 5 km s -1 reported by Pepper et al. (2013). Although evidently the uncertainties in some or all of the v sin i ⋆ determinations have been underestimated, we believe that the stacked-CCF method gives the most reliable result for v sin i ⋆ .  \nThis lends further credence to the RV-RM fit with a uniform prior on v sin i ⋆ , and as our observations of KELT-3 have a high SNR and cover the transit very well. We therefore adopt the result from this RV-RM fit, λ = -5 ± 4 · , as our best estimate for this system.", '4.6. KELT-4A': 'KELT-4 is a hierarchical triple star system discovered by Eastman et al. (2016). The primary star, KELT-4A, is an F star orbited by a hot Jupiter. The close binary system KELT-4BC is  \nFig. 8. KELT-3. Same as for Fig. 2. Here also with RVs from FIES.  \n<!-- image -->  \nFig. 9. Peak of the stacked CCFs of KELT-3. The peak of the stacked CCFs of KELT-3 created by assuming the best-fitting value of b from the RV-RM fit. The gray contours show a peak around ( v sin i ⋆ , λ ) = (6 . 0 km s -1 ,1 . 0 · ). The red contours are the same as created from the posterior shown to the right in Fig. 8.  \n<!-- image -->  \nengaged in a wide orbit with KELT-4A. TESS data are available from Sector 48 with 2-min cadence. The light curve is shown in Fig. A.1. Since the TESS cameras cannot resolve KELT-4A from KELT-4BC, we included a dilution factor in our analysis of δ M = 3 . 5 ± 0 . 1 (estimated from Table 3 in Eastman et al. 2016) to correct for the light from KELT-4BC.  \nOn the night starting January 27, 2021, we observed a transit of KELT-4Ab, starting the observations UT 01:40 and ending them UT 06:45 January 28, 2021. The exposure time was 900 s, meaning the e ff ective sampling was a spectrum every 1090 s.  \nThe observations are shown in Fig. 10 with our results from a run invoking a Gaussian prior on v sin i ⋆ . From this we got λ = 80 + 25 -22 · . Applying a uniform prior on v sin i ⋆ we found λ = 77 + 25 -23 · . Evidently, applying a Gaussian or a uniform prior on v sin i ⋆ makes little di ff erence, despite the somewhat noisy data, most likely because of the rather large impact parameter and large misalignment.  \n<!-- image -->  \nAs for KELT-3, we investigated the stacked CCFs for this system. The results coincided with the result from our RV-RM fit, and when fitting a 2D Gaussian to the grid of v sin i ⋆ and λ we obtained values of v sin i ⋆ = 6 . 0 ± 0 . 4 km s -1 and λ = 88 . 9 ± 0 . 5 · . The results are consistent with the polar configuration implied by the RV-RM fit. Our best estimate for λ for this system, 80 + 25 -22 · , is based on the RV-RM fit with a Gaussian prior on v sin i ⋆ .', '4.7. LTT 1445A': 'LTT 1445 is a triple star system consisting of three M-dwarfs. The primary star hosts an Earth-sized planet on a 5-day orbit, LTT 1445Ab, which was discovered by Winters et al. (2019). A second transiting planet, LTT 1445Ac, was reported by Winters et al. (2022), and an RV signal from a third planet was detected by Lavie et al. (2023). The orbit of LTT 1445Ab is sandwiched between the orbits of the other two planets.  \nLTT 1445A was observed by TESS in Sectors 4 and 31, and the data are available with 20-sec and 2-min cadence, respectively. The light curves show a photometric quasiperiodicity that probably arises from rotation of either star B or star C (Winters et al. 2019). To account for this rotational modulation, our model for the light curve including a GP with a kernel consisting of two stochastically driven damped harmonic oscillators. As the amplitude of the rotational modulation seemed to change between the two sectors (Winters et al. 2022), we did not require the hyperparameters to be the same in both sectors. The light curve is shown in Fig. A.2.  \nWe observed a spectroscopic transit of LTT 1445Ab employing HARPS-N. The observations were carried out on the night starting September 5th, 2020. We used an exposure time of 900 s, yielding a sampling time of 920 s. The first exposure was not useful because the telescope was erroneously aimed at LTT 1445BC. The remainder of the observations are shown in Fig. 11.  \nAs can be seen in Fig. 11 no RM e ff ect is apparent in the RV data. Formally, we found λ = 22 + 98 -83 · , implying a gentle preference for a prograde orbit. However, as described by Albrecht et al. (2011), such low SNR results are probably even more uncertain than the formal uncertainty suggests.  \nUsually, if the RV-RM e ff ect has a detectable amplitude, then the scatter of the RVs inside the transit should be larger than the  \n<!-- image -->  \nFig. 10. KELT-4A. Same as for Fig. 2. Here also with RVs from FIES.  \n<!-- image -->  \n<!-- image -->  \nFig. 11. LTT 1445A. Same as for Fig. 2, but with RVs from HARPS-N.  \n<!-- image -->  \nscatter of the RVs outside the transit. For LTT 1445A, the intransit scatter is σ = 5 . 9 m s -1 after subtracting the best-fit orbital solution. The out-of-transit scatter is σ = 7 . 1 m s -1 (dropping to σ = 6 . 0 m s -1 if the single pre-transit datum is excluded). The in-transit residuals relative to the best-fit RM model have a scatter of σ = 5 . 5 m s -1 . In light of these results, we do not consider the RM e ff ect to have been detected at all. Better results would probably require a bigger telescope.', '4.8. TOI-451A': 'TOI-451 is a young multiplanet system in the Pisces-Eridanus stream consisting of the G-dwarf planet host, TOI-451A, and a wide-orbiting (likely) M-dwarf binary (Newton et al. 2021). The planetary system around TOI-451A consists of three planets with radii 1.9, 3.1, and 4.1 R ⊕ and periods of 1.9, 9.2, and 16 days. TESS observed the system in Sector 4, 5, and 31 with a cadence of 2 min. The TESS observations of TOI-451A are shown in Fig. A.2.  \nWe observed a transit of TOI-451Ab, the innermost and smallest planet, using ESPRESSO. The observations were car-  \nried out on the night starting on September 24, 2022, from around UT 04:15 until 08:30 September 25, 2022. The exposure time was set to 720 s, resulting in a spectrum every 755 s.  \nAs reported by Newton et al. (2021), the light curve of TOI451 displays a rather large amplitude rotational modulation. We therefore adopted the same kernel as for LTT 1445 (described in the previous section) instead of the Matérn-3 / 2 kernel. We included all three planets in our transit modeling. We applied Gaussian priors with values drawn from Newton et al. (2021) for the transit parameters to the two, outer planets (c and d) for which we did not observe a spectroscopic transit. For planet b we followed our usual approach.  \nThe RV-RM analysis with a uniform prior on v sin i ⋆ gave λ = -27 + 64 -74 · and the analysis with a Gaussian prior gave -23 + 37 -40 · . All we can conclude is that the orbit is prograde.  \nWe investigated the RV scatter in the same manner we did for LTT 1445. We obtained σ = 2 . 8 m s -1 and σ = 2 . 4 m s -1 for the in-transit and out-of-transit scatter, respectively, and when subtracting the RM e ff ect from the in-transit data we found σ = 2 . 0 m s -1 . With these results we cannot be confident that the RM e ff ect was detected. For better results, one might ob-  \n<!-- image -->  \nFig. 12. The RV-RM e ff ect for TOI-451. Same as for Fig. 2, but with RVs from ESPRESSO.  \n<!-- image -->  \nve a second transit of TOI-451Ab or give up and choose to observe one of the larger planets (for which opportunities to observe transits are rarer).  \nWe did not use any results from LTT 1445A or TOI-451A in the ensemble analyses described in Section 5.2, Section 5.1, and Section 5.3. However, the results are displayed in Section 5.4, which focuses on multiplanet systems.', '4.9. TOI-813': "Eisner et al. (2020) used TESS data to discover a Saturn-sized planet on an 84-day orbit around this subgiant star. TESS data are available with 30-min cadence in Sector 2; 2-min cadence in Sectors 5, 8, 11, 30, 33, and 39; and 20-sec cadence in Sectors 27 and 61. Due to the planet's long orbital period, transits only occurred during Sectors 2, 27, and 61. Fig. A.2 shows the light curve.  \nWe observed a transit of TOI-813 b on the night January 20, 2023, using ESPRESSO. The exposure time was set to 800 s, yielding a sampling time of around 832 s. This transit was simultaneously observed with TESS. We observed the transit starting just before ingress and continued for for 6.5 hours until the target set. These observations cover roughly half of the 13.14 hr transit as shown in Fig. 13.  \nIn an attempt to cover as much of the transit as possible, the two last exposures were acquired through a large airmass ( > 2 . 2), and the RVs based on those exposures deviate from the trend established by the preceding RVs. Including the two highairmass RVs in our fit resulted in λ = -60 + 19 -21 · suggesting that the system is significantly misaligned. Omitting these two exposures yielded -32 ± 23 · . In both cases described above, a Gaussian prior on v sin i ⋆ was applied. Excluding the last two exposures while running an MCMC with a uniform prior in v sin i ⋆ , we obtained -53 + 22 -25 · for λ .  \nWe also analyzed the Doppler shadow for this system and found a value of -53 + 25 -28 · , which is more in-line with the results obtained using a uniform prior on v sin i ⋆ or including the two data points at high airmass. However, the lack of observations outside of the transit left us without a spectral CCF that is free of the RM e ff ect.  \nSince only half of the transit was covered, we found it necessary to apply a Gaussian prior on v sin i ⋆ to obtain reliable re-  \nlts. We also decided not to trust the two high-airmass RVs. We adopt the resulting value of -32 ± 23 · as our best estimate of λ for TOI-813. The RM e ff ect has a high enough amplitude to be detectable using smaller telescopes and less precise spectrographs. Another observation covering the second half of the transit is warranted, and would allow us to be judge more definitively if the system is aligned or not.", '4.10. TOI-892': 'TOI-892 b is a warm Jupiter on a circular orbit around an Ftype dwarf (Brahm et al. 2020). The system was first observed by TESS in Sector 6 with 30 min cadence, and then in Sector 33 with 2-min cadence. Fig. A.2 shows the light curve.  \nWe used HARPS-N to observe a transit of TOI-892 b. We started observations around UT 19:15 and continued until 03:15 on the night starting on January 1, 2022. The exposure time was set to 960 s, which including overhead comes out to a sampling time of around 987 s. The observations were interrupted due to high humidity.  \nThe results with and without the Gaussian prior on v sin i ⋆ are in agreement, giving λ = ! 10 + 39 -35 · and 9 + 30 -23 · , respectively. Since about half of the transit was not observed, the run including a Gaussian prior on v sin i ⋆ yielded a better constraint on λ , as one would expect.  \nWe analyzed the Doppler shadow and found λ = -12 + 18 -16 · , which di ff ers from the RV-RM results more than one would expect for two analyses of the same data. Since the the shadow shown in Fig. 15 looks convincing, we decided to adopt the value of -12 + 18 -16 · as our best estimate of λ for this system.', '4.11. TOI-1130': 'TOI-1130 is a K dwarf with two transiting planets, an inner Neptune and an outer Jupiter, discovered by Huang et al. (2020). We used ESPRESSO to observe a transit of the larger planet, TOI1130 c, on the night starting on June 1, 2021. The observations were carried out from around UT 01:30 to 06:00 on June 2, 2021. Weopted for an exposure time of 660 s, yielding a sampling time of 701 s. TOI-1130 was observed by TESS in Sector 13 in 30minute cadence and in Sectors 27 and 67 in 20-second cadence. The TESS light curve for planet c is shown in Fig. A.2.  \n<!-- image -->  \nFig. 13. The RV-RM e ff ect for TOI-813. Same as for Fig. 2, but with RVs from ESPRESSO.  \n<!-- image -->  \n<!-- image -->  \nFig. 14. The RV-RM e ff ect for TOI-892. Same as for Fig. 2, but with RVs from HARPS-N.  \n<!-- image -->  \nBecause the planets in TOI-1130 are in or near a 2:1 meanmotion resonance, the planets show relatively large transit timing variations (TTVs). The TTV amplitudes are on the order of a few hours for planet b and 15 min for planet c (Korth et al. 2023). Based on the comprehensive model of the system developed by Korth et al. (2023), the mid-transit time of planet c on June 2, 2021 should have occurred at BJD 2459367 . 6304 ± 0 . 0003 d (J. Korth, private communication). We used this result as a Gaussian prior on T 0 in our analysis.  \nThe TTVs also a ff ect the phase-folding process that is part of the construction of a high-SNR light curve. To deal with the TTVs, we identified the transits of both planet b and c in the TESS data, and fitted a model the TESS data alone, holding fixed all of the parameters except for the individual midtransit times. Then, in our RV-RM analysis, we fitted the TESS data jointly with the RV time series while applying Gaussian priors on the mid-transit times with a width of 0 . 005 d. We also applied a Gaussian prior on the stellar mean density of ρ⋆ = 2 . 98 ± 0 . 18 g cm -3 (Korth et al. 2023). Most of the transits of planet b were irrelevant, but there was one instance of overlapping transits of planets b and c.  \nOur BF analysis gave a lower projected rotation velocity for TOI-1130 (1 . 3 km s -1 ) than we found in the literature (4 . 0 ± 0 . 5 km s -1 , Huang et al. 2020). Our result is more in line with the upper limit of 3 km s -1 reported by Korth et al. (2023). From the RV-RM analysis, we obtained v sin i ⋆ = 1 . 3 ± 0 . 5 km s -1 and λ = 4 + 5 -6 · when applying a uniform prior on v sin i ⋆ . With a Gaussian prior on v sin i ⋆ , we obtained λ = 4 + 5 -6 · , which we regard as the best estimate. We need to stress that the accuracy of the result does depend on the prior on T 0 that was obtained from the model of (Korth et al. 2023).  \nWe investigated the subplanetary velocities for TOI-1130 to see if we could properly trace the deformation of the stellar lines for this grazing transit, where (unlike for K2-287) the planet only obscures part of the stellar limb. It proved more di ffi cult to locate the distortion in the first and last few in-transit spectra, as shown in Fig. 17. We also found it necessary to adopt a Gaussian prior on v sin i .  \n⋆  \nWhen fitting the subplanetary velocities using all of the data, the posterior for λ had a peak centered on 0 · , and a secondary peak near λ = 180 · . We found the secondary peak to be a consequence of including the noisiest data points from the first and  \nFig. 15. Doppler shadow for TOI-892. Same as in Fig. 4, but for TOI-892.  \n<!-- image -->  \n<!-- image -->  \nFig. 16. The RV-RM e ff ect for TOI-1130. Same as for Fig. 2, but with RVs from ESPRESSO.  \n<!-- image -->  \nlast few in-transit spectra. Out of concern over systematic errors in our model for limb-grazing transits, we repeated the fit after omitting the in-transit data that were obtained when the loss of light was smaller than 0.75% (the square data points in Fig. 17). The result was λ = 6 + 10 -9 · . Although excluding data without a firm justification is undesirable, we judge the evidence as indicating that any misalignment is smaller than about 20 · .', '4.12. WASP-50': "WASP-50 is a G9 dwarf hosting a hot Jupiter on a 1 . 9 d orbit discovered by Gillon et al. (2011). TESS data are available from observations in Sector 4 and 31 with 2 min and 20 sec cadence, respectively. Fig. A.2 shows the TESS transit light curve.  \nOnthe night starting on August 24, 2022, we observed a transit of WASP-50 b with ESPRESSO. The observations started at around UT 05:10 and continued until UT 09:40 on August 25, 2022. An exposure time of 240 s ensured a cadence of 270 s.  \nThe time series is well-sampled and the SNR is very high, exactly the case where we do not expect the results to be sensitive  \nto the prior on v sin i ⋆ . Indeed, for this system we found -2 . 9 ± 1 . 2 · when applying a uniform prior and -2 . 9 + 1 . 2 -1 . 3 · for a Gaussian prior.  \nThe high-quality data prompted us to investigate the RM effect further. We examined the subplanetary velocitiesm, which were shown in Fig. 1, but in Fig. 19 we display them in a slightly di ff erent manner (similar to the format of Fig. 7 for K2-287). As was the case for K2-287, we included two jitter terms, one for when the planet's silhouette is on or near the limb, and another when it is completely within the stellar disk. The result of our analysis using the subplanetary velocities was v sin i ⋆ = 2 . 77 ± 0 . 14 km s -1 and λ = -1 . 8 ± 1 . 0 · . The agreement with the RV-RM result is reassuring, but for our best estimate we chose the result -2 . 9 ± 1 . 2 · from the RV-RM analysis.", '4.12.1. Convective blueshift': "This term 'convective blueshift' refers to the net blueshift of the disk-integrated stellar spectrum caused by the imbalance in brightness between di ff erent components of convective cells.  \nFig. 17. Subplanetary velocities for TOI-1130. In the top panel the distortion from subtracting the in-transit CCF from the out-of-transit are shown as colored (wiggly) lines. Lines of corresponding colors with a black outline are (Gaussian) fits to the distortion in order to extract the position of the distortion denoted by the dashed vertical line. The bottom panel shows the location of the distortion plotted around the midtransit time. The squares denote epochs for which the decrease in flux from the transiting planet is less than 0.75%. The solid white (black dashed) line shows a model for an (anti-)aligned configuration.  \n<!-- image -->  \nThe upwelling gas is hotter, brighter, and blueshifted, while the sinking gas is cooler, fainter, and redshifted. Shporer & Brown (2011) showed that this e ff ect contributes to the anomalous RVs with amplitudes of a few m s -1 , and variability on a timescale of several minutes. Within our sample, WASP-50 is the only system for which the quality of the data is high enough to contemplate detecting and modeling the convective blueshift.  \nWe modelled the convective blueshift as described by Shporer & Brown (2011). In the model, the stellar disk is pixelated, and each pixel is assigned an intensity according to a limb-darkening law and a net convective velocity ( V CB) directed toward the center of the star. The radial component of each blueshifted surface element is therefore V CB cos θ where θ is the angle between the normal to the stellar surface and our line of sight. The value of V CB is expected to range from about -200 m s -1 for K-type stars to -1000 m s -1 for F-type stars (Dravins & Nordlund 1990a,b; Dravins 1999). The solar value is about -300 m s -1 (Dravins 1987). Since WASP-50 is similar to the Sun, with T e ff = 5400 ± 100 K, a reasonable expectation is V CB ≈ -300 m s -1 .  \nSince the anomalous RV caused by the convective blueshift is smallest near the limb, and WASP-50 b has a high impact parameter ( b ≈ 0 . 7), the overall e ff ect of the convective blueshift is expected to be rather weak, as shown in Fig. 20. Nevertheless, we tested whether the inclusion of the convective blueshift in our RV-RM model leads to a substantial improvement in the quality of the fit to the data,. We created a grid of λ values between -45 · and + 45 · in steps of 0 . 25 · . For each choice of λ , we computed the Bayesian information criterion (BIC; Schwarz 1978) for RV-RM models with and without the convective blueshift. The standard model has a minimum BIC for λ = -2 . 9 · (as expected given our result of -2 . 9 ± 1 . 2 · ), whereas the model with the convective blueshift has a minimum BIC for λ = -3 . 9 · . The di ff erence in the optimized BIC values is 4.5, favoring the model without the convective blueshift. We calculated the di ff erence in BIC ( ∆ BIC) between the two models at each λ . ∆ BIC changes  \nsign around λ = -4 . 9 · meaning that the net e ff ect of including convective blueshift results in slightly lower (more misaligned) values for λ .  \nEven though we did not find that including the convective blueshift significantly improves the fit, we decided to repeat our RV-RM analysis after including the convective blueshift. For convenience, we did not jointly fit the TESS light curve; we simply held b , R p / R ⋆ , P , and T 0 fixed at the best-fit values from our previous run above. We applied a Gaussian prior on V CB centered on -400 m s -1 with a width of 100 m s -1 . From this run we obtained λ = -3 . 5 ± 1 . 3 · again suggesting a slightly more misaligned system. A smaller absolute value for V CB (probably more in-line with the expected value for WASP-50) would make the di ff erence in λ from the two runs even smaller as such we do not expect convective blueshift to significantly a ff ect the results. Thus, when attempting to measure λ with a precision at the level of one degree, the results depend on how much one trusts the model for the convective blueshift and its assumed amplitude.", '4.13. WASP-59': "The WASP-59 system was discovered by Hébrard et al. (2013). The system consists of a warm Jupiter, WASP-59 b, orbiting a K-dwarf. We observed a transit of WASP-59 b with HIRES. The observations were carried out on November 23, 2016, starting at UT 05:19 to 09:46 with an exposure time of 787 s. Including overhead the sampling time came out to 835 s. The series of 20 observations spanned a transit.  \nTESS observed the system in Sector 56 with a cadence of with a cadence of 2 min. Given the long time interval that elapsed between our HIRES observations and the TESS observations, one might be concerned about the extra uncertainty in the transit ephemerides. Fortunately, we also arranged for groundbased photometric observations of a transit a few weeks after our spectroscopic observation. We used KeplerCam, which is mounted on the 1.2 m telescope at the Fred L. Whipple Observatory on Mount Hopkins, Arizona (Szentgyorgyi et al. 2005). We observed with an r ' filter on the night starting that began on the night of December 8, 2016. The light curves are shown in Fig. A.3.  \nBecause the HIRES spectra were imprinted with iodine absorption lines, we did not derive a value for v sin i ⋆ with the BF method. Instead, we adopted the value of 2 . 3 ± 1 . 2 km s -1 reported by Hébrard et al. (2013). When using this result as a Gaussian prior, we obtained λ = 16 + 17 -16 · and v sin i ⋆ = 14 . 8 ± 0 . 5 km s -1 , clearly suggesting that the Gaussian prior is not particularly informative. Indeed assuming a uniform prior on v sin i ⋆ we got v sin i ⋆ = 11 + 4 -11 km s -1 and λ = 16 + 18 -16 · , which we adopt as a out final value for λ . The resulting RM curve and posteriors are shown in Fig. 21.", '4.14. WASP-136': 'WASP-136 b is an inflated hot Jupiter orbiting an F9 subgiant star. The system was discovered by Lam et al. (2017). TESS data are available with 2-min cadence in Sector 29 and 20-sec cadence in Sector 41. The light curve is shown in Fig. A.3.  \nWe observed a transit with FIES on the night starting August 30, 2021, with observations running from UT 23:00 to UT 05:50 August 31, 2021. We opted for an exposure time of 1140 s, yielding a sampling time of 1320 s.  \nBoth of the RV-RM fits resulted in a significant misalignment for the WASP-136 system. We obtained λ = -57 + 9 -14 · when  \n<!-- image -->  \nFig. 18. The RM-RV e ff ect for WASP-50. Same as for Fig. 2, but with RVs from ESPRESSO.  \n<!-- image -->  \nFig. 19. Subplanetary velocities for WASP-50. Subplanetary velocities for WASP-50 with color-coding as described in Fig. 1. The dark shaded area correspond to the time the planet is completely within the stellar disk and the light shaded area denote ingress and egress. The errorbars include a jitter term, where two separate terms have been added in quadrature; one for the epochs for the full transit and one for the epochs when the planet is partially outside of the limb of the star. The horizontal black lines in the top panel mark v sin i ⋆ from the fit.  \n<!-- image -->  \napplying a uniform prior on v sin i ⋆ , and -38 + 16 -18 · when applying a Gaussian prior. We adopt the latter as our best estimate of the projected obliquity.', '4.15. WASP-148': "The WASP-148 system was discovered by Hébrard et al. (2020). The system consists of two giant planets orbiting a late-G dwarf with strong TTVs. The inner planet, WASP-148 b, is a transiting hot Jupiter with a period of 8.8 days. The star's projected obliquity with respect to the orbit of planet b was measured by Wang et al. (2022) to be λ = -8 + 9 -10 · .  \nWe observed a transit of WASP-148 b using HARPS-N on the night starting March 21, 2021, with observations starting around UT 01:45 to 06:15 March 22, 2021. The exposure time was set to 900 s, resulting in a sampling time of 920 s. Be-  \nFig. 20. Convective blueshift in WASP-50. Same as the left panel of Fig. 18, but also showing a model as a red solid line that includes the e ff ect of convective blueshift ( V CB = -300 m s -1 ).  \n<!-- image -->  \ncause the strong TTVs lead to extra uncertainty in the transit ephemerides, we also arranged for simultaneous photometric observations of the transit with the Multicolor Simultaneous Camera for studying Atmospheres of Transiting exoplanets (MuSCAT-2; Narita et al. 2015) mounted on the Telescopio Carlos Sánches, Tenerife, Spain. MuSCAT-2 allows for simultaneous observations in four filters; we used a griz filter set. The light curves, after detrending, are shown in Fig. A.4. When modeling these light curves, we did not use GP regression. The limbdarkening coe ffi cients tabulated by Claret et al. (2013) for these filters and for the appropriate stellar parameters are reproduced in Table 7.  \nAt first, we fitted the RM e ff ect for this system jointly with the TESS and MuSCAT-2 light curves, allowing the mid-transit time of each transit to be a freely adjustable parameter (as we did for TOI-1130, Section 4.11). We also tried fitting only the MuSCAT-2 light curves and not the TESS light curve. In both cases, we found λ = 16 ± 40 · , which is consistent with, but less precise than, the measurement reported by Wang et al. (2022).  \nAlmenara et al. (2022) performed a 'photodynamical' analysis of the system, in which the light curve is fitted with a model  \n<!-- image -->  \nFig. 21. The RM-RV e ff ect for WASP-59. Same as for Fig. 2, but with RVs from HIRES.  \n<!-- image -->  \n<!-- image -->  \nFig. 22. The RV-RM e ff ect for WASP-136. Same as for Fig. 2. Here also with RVs from FIES.  \n<!-- image -->  \nin which all the mutual gravitational interactions between the star and both planets are taken into account. They used the data presented by Hébrard et al. (2020) along with additional groundbased data and TESS data. Given their comprehensive modeling e ff ort, we decided to use the parameters from Almenara et al. (2022) as priors, and fitted only the HARPS-N RV measurements and the MuSCAT-2 light curve. The free parameters relevant to the light curve were the mid-transit time, a / R ⋆ , cos i , R p / R ⋆ , and the sum-of-limb-darkening coe ffi cients for each filter. The period was held fixed. In addition, we included a Gaussian prior on the projected obliquity using the value of λ = -8 + 9 -10 · as to incorporate all the information available for the system. Again, we performed two MCMCs obtaining a value of -2 ± 8 · for both the run with a uniform and the one with a Gaussian prior on v sin i ⋆ , respectively.  \nAs this is a dynamically interesting system, for which a large collection of transit times over a long interval can help to constrain the parameters of the planetary system, we measured the mid-transit time based on a fit solely to the MuSCAT-2 light curves. The results are given in Table 7.", '4.16. WASP-172': "WASP-172 is an F star with a bloated hot Jupiter discovered by Hellier et al. (2019a). Data from TESS are available with 30-min cadence from Sector 11 and 20-sec cadence from Sectors 37 and 38. The TESS light curve is shown in Fig. A.3.  \nOn the night starting June 1, 2022, we observed a transit using ESPRESSO with observations being carried out from UT 23:00 to UT 06:40. We opted for an exposure time of 900 s, resulting in a sampling time of 940 s. The RVs are shown in Fig. 24, where three data points are marked with gray points as they were excluded from the analysis, since they were taken at high airmass.  \nBy fitting the RV-RM e ff ect and the TESS light curve, we found λ = 118 ± 6 · and 115 ± 6 · for λ when applying a uniform prior or a Gaussian prior on v sin i ⋆ , respectively. Given the high quality of the data, and the relatively rapid rotation of the star ( v sin i ⋆ = 13 . 0 + 0 . 4 -0 . 5 km s -1 ), we decided to analyze the Doppler shadow. The results are displayed in Fig. 25. The Doppler-shadow analysis gave λ = 121 + 11 -13 · , consistent with the  \n<!-- image -->  \nFig. 23. WASP-148. Same as for Fig. 2, but with RVs from HARPS-N.  \n<!-- image -->  \n<!-- image -->  \nFig. 24. WASP-172. Same as for Fig. 2, but with RVs from ESPRESSO. The gray data points were excluded in the fit due to high airmass.  \n<!-- image -->  \nTable 7. WASP-148, photometry only.  \nTop: Results from an MCMC using only the MuSCAT-2 photometry for the transit of WASP-148 b. Bottom: The limbdarkening coe ffi cients for each filter using Claret et al. (2013).  \nRV-RM results. We adopted 115 ± 6 · as our best estimate for λ , based on the RV-RM fit with a Gaussian prior on v sin i ⋆ .  \nThe ESPRESSO transit spectroscopy presented here was also used by Seidel et al. (2023) in an investigation of the planet's atmosphere. They reported the detection of sodium and hydrogen absorption, as well as a tentative detection of iron absorption. In their study, a value of λ = 121 ± 13 · was adopted to correct for the RM e ff ect, which was taken from our Dopplershadow analysis prior to the completion of this paper.", '4.17. WASP-173A': 'WASP-173Ab is a hot Jupiter discovered by Hellier et al. (2019a). The host star is of spectral type G3V and has a wideorbiting stellar companion. The planet was also dubbed KELT22Ab by Labadie-Bartz et al. (2019), who reported an independent discovery of the transits using data from the KELT data. Here, we use the name WASP-173Ab.  \nTESS observed the system twice. Data with 2-min cadence are available from Sector 2, and with 20-sec cadence from Sector 29. The TESS light curve is shown in Fig. A.3. We observed a transit on the night starting July 23, 2022, with observations  \nFig. 25. Doppler shadow for WASP-172. Same as in Fig. 4, but for WASP-172.  \n<!-- image -->  \nbeginning at UT 05:40 until UT 10:00. The exposure time was set to 555 s, yielding a sampling time of 590 s. The RVs are shown to the left in Fig. 26.  \nOur analysis followed similar steps as for WASP-172 (see the previous section). We performed the usual RV-RM analyses and also an analysis of the Doppler shadow, for which the results are shown in Fig. 27. The RV-RM results were λ = 16 + 38 -31 · when applying a uniform prior on v sin i ⋆ and 11 + 32 -20 · for a Gaussian prior, while the shadow run resulted in a value for λ of 25 + 12 -11 · .  \nThe v sin i ⋆ derived from the Doppler-shadow analysis was 7 . 2 + 0 . 4 -0 . 4 km s -1 , which is a bit higher than the results obtained from the RV-RM anlayses (6 . 6 + 0 . 6 -0 . 8 km s -1 ) and a Gaussian prior (6 . 5 + 0 . 3 -0 . 4 kms -1 ). The RV-RM-based values are closer to the value derived from the BF of 6.6 km s -1 . As our best estimate of λ , we therefore adopt the value of 11 + 32 -20 · from our run applying a Gaussian prior on v sin i ⋆ .', '4.18. WASP-186': 'Schanche et al. (2020) discovered the WASP-186 b, a massive hot Jupiter with an eccentric orbit around a mid-F star. As noted by Schanche et al. (2020), WASP-186 was observed by TESS in Sector 17 and produced data with a 30-min cadence. In addition, TESS observed this system in Sectors 42 and 57 and the data are available with 2-min and 20-sec cadence, respectively. The light curves are shown in Fig. A.3.  \nUsing FIES, we observed a transit on the night October 11, 2021, with observations starting UT 22:30 and continuing until UT 02:50 October 12, 2021. The exposure time was set to 900 s, and the time between mid-exposures was 1080 s. The transit observations are shown in Fig. 28.  \nOur two RV-RM fits gave λ = -10 + 29 -32 · and λ = -10 + 18 -16 · for the run with a uniform and the one with a Gaussian prior on v sin i ⋆ , respectively. The values are in agreement and the prior on v sin i ⋆ leads to a more precise result for λ . Therefore, we chose λ = -10 + 18 -16 · as our best estimate for WASP-186.  \n4.19. XO-7  \nXO-7 b was discovered by Crouzet et al. (2020) as part of the XO project (McCullough et al. 2005). The transiting planet is a hot Jupiter. There is also a massive companion on an orbit with a period of at least 2 years around the G2V host star (Crouzet et al. 2020). TESS data are available with 30-min cadence in Sectors 18, 19, and 20; with 2-min cadence in Sectors 40, 47, 53, 53; and with 20-sec cadence in Sectors 59 and 60. The TESS light curves are shown in Fig. A.3.  \nWeobserved a transit of XO-7 b with FIES on the night starting on August 27, 2022, with observations starting around UT 00:15 and continuing until 04:00 August 28, 2022. We opted for an exposure time of 840 s, resulting in a sampling time of 1000 s. The RVs from the transit night are shown in Fig. 29. The egress was not completely covered, and no observations were made after the egress.  \nThe RV-RM results for λ are -4 + 27 -23 · and 0 + 20 -14 · for the runs with a uniform and a Gaussian prior on v sin i ⋆ , respectively, with the posterior distribution for the latter case shown in Fig. 29. Despite the incomplete coverage of the transit, the high impact parameter of b ∼ 0 . 7 enabled decent constraints to be placed upon λ .  \nThe posterior for λ , shown on the right side of Fig. 29, has a tail on the left side extending toward -90 · , allowing for the possibility that system is strongly misaligned. This tail is probably caused by the lack of egress / post-egress observations and also shows up in the 2 σ contours in the left panel of Fig. 29. We wondered if this unusual feature is related to the lack of egress and post-egress observations. We investigated the issue by examining the stacked spectral CCFs. The stacked CCFs showed contours that largely agree with the posteriors from our RV-RM. However, a peak is seen at around ( v sin i ⋆ , λ ) = (4 km s -1 , -60 · ) that is more consistent with the negative tail of the posterior from the RV-RM analysis, favoring a more misaligned system. A second and more complete transit observation seems warranted. For now, we adopted 0 + 20 -14 · as our best estimate, based on the run with a Gaussian prior on v sin i ⋆ .  \n<!-- image -->  \nFig. 26. The RV-RM e ff ect for WASP-173A. Same as for Fig. 2, but with RVs from ESPRESSO.  \n<!-- image -->  \nFig. 27. Doppler shadow for WASP-173A. Same as in Fig. 4, but for WASP-173A.  \n<!-- image -->', '4.20. WASP-26': 'WASP-26 b is a hot Jupiter orbiting an early G-type star, discovered by Smalley et al. (2010). Albrecht et al. (2012b) analyzed the RV-RM e ff ect based on data obtained with the Planet Finder Spectrograph (PFS; Crane et al. 2010) mounted at the 6.5 m Magellan II telescope. The RV time series is shown in the left half of Fig. 30. They reported λ = -34 + 36 -26 · .  \nThis case is unlike all of other cases presented in this paper in that we did not obtain any new data. Rather, we wanted to take advantage of TESS data that did not exist at the time of the prior study, and correct an error that was made in the prior study. The error was the adoption of an incorrect prior constraint on the transit duration ( T 12, in the terminology of Albrecht et al. 2012b), which can be seen in their Table 3. This error resulted in an incorrect value of b that a ff ected the results for λ .  \nRather than placing a prior on T 21, we followed our usual procedure of fitting the TESS light curve jointly with the RVs. We used TESS 20-sec data from Sector 29, shown in Fig. A.5). As in Albrecht et al. (2012b), we applied a Gaussian prior on  \nv sin i ⋆ of 2 . 4 ± 1 . 5 km s -1 , and as usual, we applied Gaussian priors on the macro- and microturbulent velocities of 4 . 0 ± 1 . 0 km s -1 and 1 . 1 ± 1 . 0 km s -1 , respectively. We also applied a Gaussian prior on T 0 taken from Anderson et al. (2011), as did Albrecht et al. (2012b). The reason for this choice is that about 10 years elapsed between the spectroscopic transit observations and the TESS observations; small systematic uncertainties in the TESS-based ephemeris might propagate into large uncertainties in the relevant transit time. We found a projected obliquity of -16 + 14 -10 · . Fig. 30 shows the RV-RM curve and posteriors.', '4.21. Results': 'Table 8 gives our best estimates of λ and v sin i ⋆ for each system. For cases in which the rotation period could be measured, the table also includes the stellar inclination, i ⋆ , calculated as per Masuda & Winn (2020), and the true obliquity. Other key parameters are given in Table B.1.  \n)  \n1  \n-  \ns  \n(m  \nV  \nR  \n)  \n1  \n-  \ns  \n(m  \nResiduals  \n100  \n50  \n0  \n50  \n-  \n-  \n-  \n100  \n100  \n0  \n100  \n-  \n<!-- image -->  \n-  \n1  \n0  \n1  \n2  \nHours From Midtransit  \nFig. 28. The RV-RM e ff ect for WASP-186. Same as for Fig. 2. Here also with RVs from FIES.  \n<!-- image -->  \nFig. 29. The RV-RM e ff ect for XO-7. Same as for Fig. 2. Here also with RVs from FIES.  \n<!-- image -->', '5. Discussion': 'We now return to the four questions posed in the introduction about the distribution of obliquities of stars with transiting planets. We combined the λ measurements presented in this paper (and listed in Table 8) with those measured previously. Specifically, we used the catalog assembled by Albrecht et al. (2022) with the help of the TEPcat catalog (Southworth 2011), and several other measurements that have appeared in the more recent literature (up until February 26th, 2024). Table B.2 gives the key parameters for these additional systems, along with the systems presented in this study. Fig. 31 shows plots of λ and ψ as a function of T e ff , a / R ⋆ and M p / M ⋆ .  \nAs mentioned we excluded our results for LTT 1445Ab and TOI-451Ab in the ensemble analyses below. Furthermore as this study builds upon the sample presented in Albrecht et al. (2022) systems excluded in there are also excluded here (based on the criteria outlined in their Appendix A).', '5.1. A Preference for Polar Planets?': "Albrecht et al. (2021) investigated a sample of 57 systems for which ψ has been determined by combining measurements of the projected obliquity ( λ ) and the stellar inclination angle ( i ⋆ ). Of the 19 values of ψ that are inconsistent with 0 · , 18 are between 80 · and 125 · . They hypothesized that the misaligned systems show a preference for approximately polar orbits. Siegel et al. (2023) and Dong & Foreman-Mackey (2023) followed up by replicating the preceding work and analyzing the larger sample of systems for which only λ has been measured but not i ⋆ . The larger sample did not show a significant preference for polar orbits, suggesting that the sample analyzed by Albrecht et al. (2021) was somehow biased toward polar orbits.  \nWe revisited the issue with a sample that is about 40% larger. The current tally of systems with λ measurements is 205, of which ψ is known for 87. In the analysis that follows, when considering the ψ measurements, we decided to omit the 6 systems for which the gravity-darkening technique was employed, because parameter degeneracies make it di ffi cult to measure ψ when it is near 0 · , 90 · , or 180 · . When considering only the  \n2  \n<!-- image -->  \nFig. 30. The RV-RM e ff ect for WASP-26. Same as for Fig. 2, but with RVs from PFS.  \n<!-- image -->  \nFig. 31. Obliquities and projected obliquities as a function of T e ff (left), a / R ⋆ (center), and M p / M ⋆ (right). Triangles indicate the a measurement whose value is outside of the plotted range. In all panels, the color of each data point conveys T e ff using the color scale shown in the lower left panel. The contours are KDEs illustrating the density of measurements in a given parameter space.  \n<!-- image -->  \nλ measurements, we retained the 6 gravity darkening measurements.  \nFirst, we repeated the statistical tests performed by Albrecht et al. (2021) using the enlarged set of ψ measurements. We selected the 27 systems for which cos ψ < 0 . 75 and computed two statistics that quantify clustering: the dispersion around cos ψ = 0 and the standard deviation relative to the mean. We then used Monte Carlo simulations to find the probability of obtaining clustering statistics at least as large as those that were observed if the obliquities were drawn from an isotropic distribu-  \nn (i.e., a uniform distribution in cos ψ ). A subtlety of these calculations is that for given values of the observables | cos i | (from the transit impact parameter), | sin i ⋆ | (from rotational broadening) and λ (from the RM e ff ect), there are two closely-spaced solutions for ψ :  \ncos ψ = | sin i ⋆ sin i | cos λ ± | cos i ⋆ cos i | . (3)  \nFollowing Albrecht et al. (2021), our Monte Carlo simulations take this discrete degeneracy into account by randomly choosing one of the two solutions in each iteration. The results were  \nTable 8. Final results for v sin i ⋆ , λ , i ⋆ , and ψ .  \np = 9 . 4 × 10 -4 for the dispersion and p = 3 . 8 × 10 -4 for the standard deviation. Applying the same test to the sample of 14 systems that were available to Albrecht et al. (2021) and were not based on gravity darkening gives p = 1 . 4 × 10 -2 and 2 . 1 × 10 -3 , respectively. Thus, the inclusion of 13 new data points has reduced the p -values by factors of 15 and 5.5, allowing a firmer rejection of the null hypothesis.  \nNext, we used the 'DFM' (Dong & Foreman-Mackey) code 5 provided by Dong & Foreman-Mackey (2023) to fit a model to the obliquity data consisting of the sum of two Beta distributions, representing the aligned and misaligned populations. 6 Each Beta distribution is characterized by three hyperparameters: a weight w , a mean µ , and an inverse variance κ . We adopted the same priors on the hyperparameters as outlined in Dong & ForemanMackey (2023). Table B.4 gives the results for the hyperparameters when the DFM code is applied to di ff erent samples of the data, and Fig. 32 shows the posterior probability density for cos ψ .  \nThe left panel of Fig. 32 shows the results of analyzing the complete sample of projected obliquity measurements. There is a strong peak at cos ψ = 1 representing the aligned systems, and a much broader distribution representing the misaligned systems. The distribution falls o ff toward retrograde systems  \n(cos ψ < 0 . 75) but otherwise there is little or no evidence for a concentration near cos ψ = 0, in agreement with the findings by Dong & Foreman-Mackey (2023) and Siegel et al. (2023).  \nThe middle panel of Fig. 32 shows the results of analyzing only those systems for which constraints on i ⋆ are available, but without using the i ⋆ information to perform the inference. The inferred probability density of cos ψ does not look too di ff erent from the case described above.  \nThe right panel of Fig. 32 shows the results of analyzing only those systems for which constraints on i ⋆ are available, and using the i ⋆ information to perform the inference. This is the type of sample that led Albrecht et al. (2021) to propose that there is a peak near cos ψ = 0, and indeed, we find a similar result.  \nThe preceding results suggest that the obliquity distribution of the systems for which measurements of both λ and i ⋆ have been reported is di ff erent from that of the systems for which λ was measured and no constraints are available on i ⋆ . This is true even though we have excluded inclination measurements obtained via gravity darkening, thereby eliminating the bias discussed by Siegel et al. (2023).  \nIt is probably relevant that the λ & i ⋆ sample is enriched in the types of systems that are known to show frequent misalignments: sub-Saturns (0 . 025 MJ < M p < 0 . 2 MJ), and hot Jupiters around hot stars ( a / R ⋆ < 7, M p > 0 . 3 MJ, T e ff ≥ 6250 K). Stellar inclination measurements in those systems almost all come from the v sin i ⋆ method - the combination of a photometric rotation period, rotational Doppler broadening of the spectral absorption lines, and the stellar radius. In the entire sample, there are 17 sub-Saturns, of which inclination measurements are available for 11, of which 5 have cos ψ < 0 . 75. The entire sample also contains 41 hot Jupiters around hot stars, of which 15 have inclination measurements. After excluding the gravity-darkening measurements, we are left with 10 hot Jupiters, all of which have cos ψ < 0 . 75.  \nBy combining the 5 misaligned hot Neptunes and the 10 misaligned hot Jupiters, we constructed a sample of 15 systems and subjected it to the same Monte Carlo experiments described above to test the null hypothesis that they are drawn from an isotropic distribution. For the dispersion-around-zero and standard deviation, we found p = 3 . 1 × 10 -3 and p = 3 . 5 × 10 -3 , respectively. We also applied the DFM code to a sample consisting only of sub-Saturns, and a sample consisting only of hot Jupiters around hot stars, with results that are given in Table B.4 and displayed in Fig. 33, in the same format as in Fig. 32. For the sub-Saturns, a central peak appears to be favored by the data, although the uncertainty in the distribution is large due to the small number of systems (17). A central peak is even more pronounced in the sample of hot Jupiters around hot stars (41 systems including the measurements employing the gravity darkening method). For the subset of systems with inclination measurements from this hot Jupiter sample, since all of the systems are misaligned, we modeled the distribution with a single Beta distribution rather than using two components. There is a clear peak at cos ψ ≈ -0 . 2 or ψ ≈ 100 · . Indeed, all of these 10 systems are consistent with ψ = 90 · with varying degrees of uncertainty, and a simple kernel-density estimate of the distribution shows a peak at cos ψ = 0. Both distributions are shown in Fig. 33.  \nTo investigate the importance of the priors assumed for the hyperparameters in the 'DFM' code for these subpopulations, we tried adopting less informative priors on the variance, κ or rather log κ , hyperparameter. Here we tried a uniform prior instead of the Gaussian prior used in Dong & Foreman-Mackey (2023, i.e., U ( -4 , 10) instead of N (0 , 3)). The resulting hyperparameters are tabulated in Table B.4 and distributions are dis-  \nFig. 32. Hierarchical Bayesian inference of the obliquity distribution using the code of Dong & Foreman-Mackey (2023). The left panel is for all systems for which λ was measured, the middle panel is for systems for which both λ and i ⋆ were measured but without using the i ⋆ information, and the right panel is for systems for which both λ and i ⋆ were measured and making use of the i ⋆ information.  \n<!-- image -->  \nayed in Fig. 34. As the variance ( κ ) is now poorly constrained the distributions do look somewhat di ff erent, and the peaks for the misaligned populations, especially for the hot Jupiter sample, are not as sharp. At present there are probably too few systems or too little information to properly infer the distributions for these two subpopulations, without assuming some value for the variance to describe the morphology. When applying the information on i ⋆ we have for the subset of 10 systems, the resulting distribution looks very similar to that in Fig. 33.  \nAlthough the results are too prior-dependent to be sure, these results lead us to hypothesize that the preponderance of polar planets is a phenomenon specific to hot Jupiters around hot stars and possibly also sub-Saturns. The hot stars, in particular, are more amenable to the determination of i ⋆ via the v sin i method, given their higher rotation velocities, possibly explaining their over-representation in the λ & i ⋆ sample.  \nIf there truly is a preponderance of systems with nearly polar orbits, it would seem to be a clue about the history of these systems. For example, for the sub-Saturn population, Petrovich et al. (2020) proposed that nearly polar orbits of Neptune-mass planets can result from a secular resonance between the nodal precession frequencies induced by the protoplanetary disk and by an outer companion. One might test this theory by searching for distant massive companions to the polar-orbiting planets. In any case, Petrovich et al. (2020) note that their theory is not applicable to hot Jupiters.  \nTidal evolution might lead to a pile-up of polar orbits of hot Jupiters around hot stars. Lai (2012) demonstrated that orbital orientations can linger near 90 · when damping is from the dissipation of inertial waves in the convective envelope of a star. In general, this theory requires the planet's orbital angular momentum to be larger than the star's spin angular momentum to avoid producing more systems near 180 · than are observed (e.g. Li & Winn 2016).", '5.2. Obliquity - Eccentricity correlation?': "Explaining the existence of giant planets on orbits well within the 'ice line' of the protoplanetary disk is one of the oldest unsolved problems in exoplanetary science. Disk-driven migration is one possibility, and would generically yield low eccentricities and obliquities (e.g., Lin et al. 1996; Dawson & Johnson  \n2018). In this scenario, the observed cases of high eccentricity and obliquity arise from other processes. An alternative idea is to combine eccentricity-raising dynamical interactions - such as planet-planet scattering or Kozai-Lidov cycles - and tidal dissipation (e.g., Fabrycky & Tremaine 2007; Chatterjee et al. 2008; Nagasawa et al. 2008). If the eccentricity reaches a high enough value, or more to the point, if the periastron distance reaches a low enough value, the damping of tidal oscillations in either the planet or star can convert orbital energy into heat and thereby shrink the semi-major axis.  \nIn many of the proposed scenarios for raising of the eccentricity, the orbital inclination relative to the initial plane is also increased, which could lead to a high stellar obliquity. Therefore, if the dynamical scenarios are correct, one might expect to see a correlation between e and λ for the a ff ected planets. An important caveat is that tidal damping tends to reduce both the eccentricity and stellar obliquity, but on timescales that might di ff er by orders of magnitude. This is not only because eccentricity and obliquity tides excite di ff erent types of perturbations, but also because eccentricity damping can take place inside either body while significant obliquity damping requires the energy to be dissipated inside the star (Ogilvie 2014). Thus, tides might confound the interpretation of any correlation. For example, if eccentricity damping is much faster than obliquity damping, then any initial correlation might be diminished on astronomically short timescales. Or, tidal dissipation might increase the correlation by boosting the number of lowe , lowλ systems even if there were little or no correlation before tidal evolution.  \nRegarding the latter point, a / R ⋆ and λ are known to be positively correlated for cool stars (Albrecht et al. 2022). This has been interpreted as evidence for tidal dissipation, since dissipation rates should decline sharply with a / R ⋆ . Furthermore, a / R ⋆ is known to be positively correlated with e , presumably for the same reason (although the causal variable might be a / R p in that case). Taken together, these two correlations might cause e and λ to be correlated without invoking any new explanation.  \nNevertheless, it is well worth testing for a correlation between e and λ in the current sample and Rice et al. (2022a) investigated both of these parameters. Specifically they found at the Kraft break, the cumulative sum of λ for systems with eccentric orbits di ff ers by 6 . 5 σ from the λ distribution of systems with circular orbits. For systems with giant-planet they obtain an  \nFig. 33. Hierarchical Bayesian inference of the obliquity distribution for subsets of the sample that are suspected of being especially prone to misalignment. The left panel is for all sub-Saturns with λ measurements, without using any i ⋆ information even when it is available. The middle panel is for all hot Jupiters around hot stars with λ measurements, without using any i ⋆ information even when it is available. The right panel is for hot Jupiters around hot stars for which both λ and i ⋆ are known and utilized. In the right panel, the black upward tick marks on the horizontal axis are the individual values of cos ψ , and the horizontal bands surrounding the tick marks convey the uncertainties. The solid gray curve is a KDE estimated from these measurements.  \n<!-- image -->  \nFig. 34. Hierarchical Bayesian inference of the obliquity distribution for the same subsets as in Fig. 33, but here a uniform prior, U ( -4 , 10), was applied to log κ in the 'DFM' code.  \n<!-- image -->  \n8 . 7 σ result. Their procedure was as follows, as we understand it. They selected systems with periastron distances ≤ 0 . 1 AU for which a value of λ was reported in the TEPCat obliquities table, looked up e in the NASA Exoplanet Archive (NEA) Confirmed Planets list, and separated the systems into two samples: an eccentric sample ( e ≥ 0 . 1), and a circular sample ( e = 0). They omitted borderline systems with 0 < e < 0 . 1. They arranged the planets in order of the star's e ff ective temperature, giving each planet a position i in a line of N planets. They examined the cumulative sums of | λ | for the eccentric sample, with the j th cumulative sum defined as the sum of all | λ i | values for i ≤ j . Since the eccentric sample contained fewer systems and had a di ff erent distribution of stellar types than the circular sample, Rice et al. (2022a) used a Monte Carlo procedure to make 'matching' eccentric samples for comparison. They created 5000 matching circular samples by drawing randomly (with replacement) from the circular sample the same number of stars with T e ff < 6100 K and the same number of stars with T e ff > 6100 K as in the eccentric sample. They then compared the N th cumulative sums  \nof the eccentric sample and the ensemble of matching circular samples, and found them to di ff er by 6 . 5 σ , where σ was determined from the spread amongst the 5000 trials. For convenience we will denote the 'number-of-sigma' by the symbol ∆ .  \nThey also tried some variations on the samples. When performing the same steps but with the samples restricted to giant planets, defined as M p > 0 . 3 MJ, they found ∆ = 8 . 7. To test the notion that a / R ⋆ is a 'hidden variable' responsible for the e /λ correlation, they created samples that were not divided by eccentricity, but rather on whether a / R ⋆ is greater or smaller than 12. In this case, they found ∆ = 4 . 9, and since this is smaller than 6.5 they concluded that the variation in a / R ⋆ cannot explain the entire e ff ect they observed between e and λ .  \nWith the enlarged sample, we decided to revisit these intriguing results. After implementing the method of Rice et al. (2022a) to the best of our ability, we found σ = 2 . 1 for the whole sample shown in Fig. 35 and σ = 4 . 1 for planets with M p > 0 . 3 MJ. To check on our code, we tried reproducing their results using the systems known at that time and were unsuccessful. Further in-  \nstigation showed that the main reason for the discrepancy can be traced to 6 systems in the eccentric sample that should have been in the circular sample. These 6 systems (with details in the footnotes) are HAT-P-23 b and HAT-P-32 b 7 KELT-6 b 8 .; and HATS-14 b, HATS-70 b, & Kepler-63 b 9 Taken together, these errors caused σ to appear larger than in reality.  \nWe tried some variations on the samples. When including systems with 0 < e < 0 . 1 in the circular group, instead of omitting them, we found ∆ = 1 . 8. When focusing on giant planets ( M p > 0 . 3 MJ), we obtained ∆ = 4 . 6. We note that there are only 10 known eccentric giant planets around cool stars, making the results sensitive to an error for an individual system. When we tried dividing the sample according to whether a / R ⋆ is smaller or greater than 12, instead of by eccentricity, we found ∆ = 8 . 1 as shown in Fig. 36. Thus, in our case, the σ statistic is larger when testing the e ff ect of a / R ⋆ than when testing for the e ff ect of eccentricity.", '5.3. Tidal realignment': "The main evidence that tidal obliquity damping occurs on astrophysically relevant timescales, despite many researchers' early expectations (see, e.g., Queloz et al. 2000; Winn et al. 2005), is that the types of systems where one might expect tidal dissipation rates to be fastest - with massive, close-orbiting planets and host stars with thick convective envelopes - do indeed have an obliquity distribution that is concentrated around 0 · . And conversely, systems for which misalignment is proposed to be part of the formation process but for which tidal dissipation rates are expected to be much slower - somewhat wider-orbiting giant planets, and stars with radiative envelopes - have a broad obliquity distribution. For some details on these findings, see Schlaufman (2010); Winn et al. (2010a), and Hébrard et al. (2011). To demonstrate this point more explicitly, Albrecht et al. (2012b) constructed a crude 'tidal dissipation figure of merit' for each system, meant to be a proxy for the tidal dissipation rate, and showed that the obliquity distribution does indeed narrow for the systems where the figure of merit predicts more rapid dissipation.  \nHowever, the precise mechanisms of tidal dissipation are not well understood and are probably quite complex, depending on forcing frequency, stellar structure, rotation rate, etc. Another problem with tidal damping is that the expected endpoint of tidal evolution for most of the known close-orbiting giant planets is tidal orbital decay and destruction (Brown et al. 2011). Thus, an observed planet around a star that experienced tidal damping must somehow have avoided orbital decay before the star was aligned to the observed degree.  \nSeveral theories exist that might prevent or at least delay planet destruction for long enough to be consistent with the data. The orbital angular momentum of some lucky planets may exceed the host's spin angular momentum, causing the star to align before its orbital decays (see, e.g., Hansen 2012; Valsecchi &  \nRasio 2014; Dawson 2014). If tidal damping is primarily driven by the dissipation of inertial waves in a star's convective zone, then alignment could be faster than orbital decay (Lai 2012; Lin &Ogilvie 2017; Damiani & Mathis 2018). The damping of gravity waves in radiative cores might similarly cause alignment before tidal in-spiral (Zanazzi et al. 2024). Another idea is that if the e ffi ciency of tidal dissipation drops sharply with increasing forcing frequency (shorter orbital periods), then alignment can be arranged to occur before tidal evolution slows down dramatically, allowing the orbit to survive (Barker & Ogilvie 2010; Penev et al. 2018; Anderson et al. 2021; Ma & Fuller 2021). Alternatively, and more speculatively, tides from the planet might only a ff ect a star's outer convective zone, allowing the planet to donate a smaller amount of its orbital angular momentum to spin up the star and save itself (Winn et al. 2010a). With all of this theoretical uncertainty, is there a relatively model-independent way to seek corroborating evidence for tidal obliquity damping?  \nThe Sun's obliquity is 7 . 155 · relative to the ecliptic (Beck & Giles 2005) and 6 . 2 · relative to the Solar System's invariable plane (see, e.g., Souami & Souchay 2012; Gomes et al. 2017). If we assume this angle to be representative of the degree to which stars and their planetary systems are initially aligned - admittedly, a big assumption - then we would expect some systems that have experienced strong tidal damping to have obliquities much smaller than 6 · . This can be tested with high-precision RM observations of the types of systems where strong tidal damping is expected, and comparisons to systems where it is not expected.  \nFocusing on only the most precise measurements ( σλ < 2 · ) in cool hot Jupiter systems, Albrecht et al. (2022, Section 3.1.8) found the dispersion of λ measurements to be 0 . 91 · . Given an average measurement uncertainty of 0 . 82 · , this should probably be regarded as an upper limit on the true dispersion of a few degrees. There are drawbacks to this simple approach. First, one must decide what constitutes a 'precise' measurement. Somewhat arbitrarily, Albrecht et al. (2022) chose a threshold precision of 2 · . Second, it ignores the information contained in the much more numerous and less precise measurements.  \nWe tried overcoming these problems with a Hierarchical Bayesian Model (HBM), following up on the work by Siegel et al. (2023), who carried out similar tests to those described below. (To eliminate any suspense, we obtained similar results, too.) The hyper-parameter of interest is the dispersion of the underlying λ distribution around 0 · , independent of measurement uncertainties. To identify systems for which rapid tidal dissipation is expected, we required T e ff < 6100 K, a / R ⋆ < 10, and M p > 0 . 3 MJ. We also required σλ < 50 · . This led to a sample of 60 systems, only one of which stood out: TOI-858Bb (Hagelberg et al. 2023a), which is located in a wide binary system and has a retrograde orbit. Since tidal dissipation is evidently not the only important obliquity-altering process in this system - perhaps due to the stellar companion - we decided to omit it from consideration, leaving a sample of 59 systems. Fig. 37 shows the key characteristics of those 59 systems (as well as some systems with higher values of a / R ⋆ ). The data points with a precision in λ better than 2 · are highlighted; the measurement of WASP-50 of -2 . 9 ± 1 . 2 · presented in this paper is included in this sample.  \nThe result of the HBM was a dispersion in λ of 1 . 4 ± 0 . 7 · , depicted in Fig. 37 as a Gaussian function with a width of 1 . 4 · . For prograde systems the expected endpoint of the alignment process is λ = 0 · , which is why our HBM did not include a hyperparameter for the mean of the obliquity distribution. However, as an experiment to check for surprises or systematic errors (such as the neglect of the convective blueshift) we repeated the analysis after allowing for the mean to be a second hyperparameter. The  \nFig. 35. Obliquity distribution for stars hosting planets on eccentric and circular orbits. Left: The cumulative sum of | λ | for eccentric planets shown in red, sorted by T e ff . The black line is the average of 5000 randomly sampled sets from the circular population with the gray shaded area denoting the 1 σ and 2 σ intervals. The vertical blue line denotes the Kraft break ( T e ff = 6100 K). Right: Histogram showing the distribution of the circular population and the value of the eccentric population at the Kraft break. The Gaussian is a fit to the histogram. The circular and eccentric samples only di ff er by 2 . 1 σ , meaning the eccentric population does not appear to be significantly more misaligned. Adapted from Rice et al. (2022a).  \n<!-- image -->  \nFig. 36. Obliquity distribution for stars hosting planets on close-in and wide orbits. Same as Fig. 35, but for systems with planets on close-in or wide orbits with the division given for a / R ⋆ = 12. The systems with planets on wider orbits appear to be significantly more misaligned compared to those with planets closer in.  \n<!-- image -->  \nresults were a a dispersion of 1 . 5 ± 0 . 8 · and a mean of -0 . 2 ± 0 . 7 · , giving no clear evidence for surprises or systematic errors.  \nIn the light of these results, our measurement of WASP-50 ( -2 . 9 ± 1 . 2 · ) is more interesting than we originally anticipated. According to the rough criteria described above, the tidal dissipation timescale for this system should be faster than in most of the other systems in the sample. Gillon et al. (2011) reported  \nan age of 7 . 0 ± 3 . 5 Gyr based on comparing the usual photometric and spectroscopic observables to stellar-evolutionary models. However, the estimated age based on the chromospheric activity and rotational period was only 0 . 8 ± 0 . 4 Gyr (which could also be a result of tidal spin-up, e.g., Tejada Arevalo et al. 2021). If the old age is correct, the system is more likely to have had enough time to achieve good alignment than if the young age is correct.  \n/star  \nFig. 37. Precise projected obliquity measurements. Here we highlight the most precise ( σλ ≤ 2 · ) measurements of λ for systems with cool stellar hosts ( T e ff < 6100 K). WASP-50 is highlighted with a black center. The color-coding is done according to T e ff and the marker sizes scale with the tidal alignment timescale (as given by Zahn 1977). The gray shaded area shows the Gaussian distribution inferred from cool hosts ( T e ff < 6100 K) with close-in ( a / R ⋆ < 10), massive ( M p < 0 . 3 MJ) planets. The horizontal lines demarcate the Solar System obliquity of 7 . 155 · (Beck & Giles 2005). The light gray error bars show the less precise measurements. Adapted from Albrecht et al. (2022, Figure 10).  \n<!-- image -->  \nAs discussed in Section 3.2.2, Gillon et al. (2011) reported a rotation period for WASP-50 of 16 . 3 ± 0 . 5 d or potentially twice this value. Using these values ψ comes out to 34 + 12 -9 and 5 . 0 + 1 . 9 -3 . 0 , respectively. The value of the rotation period has dramatic consequences for the inferred orientation of the system; the shorter period we identified in the TESS light curve results implies an even more misaligned system, which seems worth checking on, given that so many other similar systems are well-aligned. Future TESS observations might be able to resolve the issue.  \nConsidering planets on wider orbits, with 10 ≤ a / R ⋆ < 20, the HBM yielded a dispersion of 8 ± 2 · . The spread in λ is therefore seen to increase with orbital separation (as also discussed in Section 5.2). When restricted to planets with 20 ≤ a / R ⋆ < 30 (a sample of only 7 systems), the HBM yielded a dispersion of 4 ± 3 · , which is consistent with the 10 ≤ a / R ⋆ < 20 results. It would be interesting to obtain more measurements for planets with large separations to bolster these results and possibly even investigate the separation-dependence of tidal dissipation rates.", '5.4. Compact multi-transiting systems, with and without outer companions': "At a time when only five λ measurements had been made for stars with more than one known transiting planets, Albrecht et al. (2013) noted that they were all consistent with good alignment. The number is now 22, including four systems we have added to the tally in this work. For WASP-148, λ had already been measured by Wang et al. (2022), and our result is consistent with theirs but with lower precision. For TOI-1130, we found λ = 4 + 5 -6 · , while for two other systems, LTT 1445A and TOI451A, our formal results were that the orbits are prograde, but given the low SNR of our measurements we do not include them  \nin the following discussion. We do include an additional 4 systems for which there are good constraints on the stellar inclination even though λ is not known. Our discussion is therefore based on the properties of 24 systems.  \nFig. 38 is a family portrait of these systems. The Solar System is shown at the top for scale and reference. There are 19 systems for which the data are consistent with zero obliquity. There are 5 systems with large obliquities relative to at least one transiting planet: HD 3167, Kepler -56, HIP 41378, Kepler -129, and K2-290. HD 3167 has a very unusual architecture in which the planets appear to be on perpendicular orbits (Bourrier et al. 2021a). This case is illustrated in Fig. 38 as the outer planet floating over Kepler -129. HIP 41378 is a five-planet system in which Grou ff al et al. (2022a) measured the projected obliquity with respect to planet d to be 57 + 26 · -18 . However, due to the long duration ( ∼ 12 . 5 hr and P ∼ 278 d), a full transit was not observed, preventing us from having high confidence in the result.  \nIn Fig. 38 Kepler -30 is grouped with the aligned systems, but there are conflicting indications about the obliquity. By using two di ff erent methods Sanchis-Ojeda et al. (2012) found λ = -1 ± 10 · and λ = 4 ± 10 · , indicating a well-aligned system. However, Morgan et al. (2023) found i ⋆ = 43 + 15 · -9 , indicating a large misalignment, based on the combination of the rotation period, v sin i ⋆ , and R ⋆ . However, this result depends critically on the value v sin i ⋆ = 2 . 0 ± 0 . 2 km / s reported by Fabrycky et al. (2012), and the measurement uncertainty of 0.2 km / s is likely to have been underestimated. Systematic errors of at least 0.5 km / s are typical, due to e ff ects such as turbulent and instrumental broadening. A value of 3 . 0 ± 0 . 5 km s -1 would allow for consistency with alignment ( i ⋆ ≈ 90 · ).  \nIn contrast to hot Jupiter systems, where misalignments might have occurred in the tumultuous process of higheccentricity migration after disk dispersal, the misalignments in multiplanet systems might have occurred while the disk was present. Gravitational interactions between a planet-forming disk, a binary companion, and a spinning host star can be enhanced when the disk looses mass, magnifying the host star's tilt through a secular resonance (e.g. Batygin & Adams 2013). For a binary companion to generate spin-orbit misalignments, the stellar spin must be gravitationally decoupled from the planet forming in the disk, otherwise the star's spin axis and planet's orbit normal precess in-unison (eq. 80 of Zanazzi & Lai 2018). The spin is decoupled from the orbits of LTT 1445Ab, and K2290Ab,c, but not for TOI-451Ab, when torqued by their binary companions.  \nSimilarly, sweeping secular resonances from a dispersing disk can enhance gravitational interactions between forming planets, magnifying their inclinations (e.g. Ward 1981). However, for planetary companions to generate misalignments, they must also supply the system with a large angular momentum deficit to misalign an inner planet. Based on the work by Petrovich et al. (2020); Zanazzi & Chiang (2024), we assume a planetary companion (mass m 2, semi-major axis a 2) could secularly misalign an inner planet (mass m 1, semi-major axis a 1) if they are well separated ( a 2 / a 1 > 2), and they have a large angular momenta ratio (( m 2 a 1 / 2 2 ) / ( m 1 a 1 / 2 1 ) > 10, see e.g. eq. 22 of Zanazzi & Chiang 2024). The planet pairs Kepler -56 c,d and Kepler -129 c,d meet this criteria, while Kepler -25 c,d, WASP 47 b,c, and TOI-1130 c,d do not, although we note that the companion TOI-1130 d is as of yet poorly constrained (Korth et al. 2023). LTT 144A might be aligned when a misalignment is predicted by Zanazzi & Lai (2018), while the other systems with distant, massive companions appear consistent with spin-orbit  \nmisalignments being generated (or not) while gas disk is present (Fig. 38).  \nSecular interactions from companions exciting an inner planet's inclination during the disk-hosting phase is not the only mechanism which can explain misalignments in multiplanet systems. For instance, Gratia & Fabrycky (2017) postulate a dynamical instability between three outer planets could have misaligned Kepler -56 b and c, while Best & Petrovich (2022) find secular chaos induced by the two binary companions of K2290A is su ffi cient to misalign K2-290Ab and c, with both mechanisms taking place after disk dispersal. Additional stellar obliquity measurements in multiplanet systems with companions (e.g. LTT 1445A), and further constraints on the presence or absence of companions in inclined multiplanet systems (e.g. HD 3167, HIP 41378), could constrain which mechanisms are responsible for generating misalignments.  \nFor the multiplanet systems which have obliquities consistent with alignment (again barring TOI-451A and LTT 1445A), we applied the HBM framework described in the previous section to infer the mean µ and dispersion σ of the intrinsic distribution of obliquities. The results were µ = 0 . 0 ± 1 . 7 · and σ = 2 . 0 ± 1 . 5 · . When we removed Kepler -30 c from the sample, out of the concern raised above, we found µ = 0 . 0 ± 1 . 8 · and σ = 2 . 1 ± 1 . 6 · . Thus, the obliquity dispersion of the host stars of multis is similar to that of cool stars with hot Jupiters. However, tidal obliquity damping is expected to be much weaker for the multis than for the hot Jupiters, because of the planets' lower masses and relatively wider orbits. The low obliquities of the multis might therefore be primordial, i.e., due to the flatness of the protoplanetary disk and its strong coupling to the accreting young star. From this point of view, the solar obliquity of 6 · may be unusually high - and if the high solar obliquity is simply a peculiarity of the Solar System, then the evidence for tidal obliquity damping presented in the previous section would be undermined.  \nWe finally note that an additional aligned system could potentially be the HD 148193 system ( λ = 14 + 7 -8 · , Section 4.2), if the inner transiting candidate is confirmed as a bona fide planet at some point.", '6. Conclusions': "We presented new observations of the RM e ff ect for 19 stars with transiting planets. We combined these measurements with results from the literature to arrive at a sample of 205 systems for which λ has been determined with reasonable confidence. With this sample, we revisited four questions:  \n- i) Is there a preponderance of perpendicular planets? First raised as a possibility by Albrecht et al. (2021), and then questioned Siegel et al. (2023) and Dong & ForemanMackey (2023), our updated sample does not settle the matter. It remains the case that a 'polar peak' is seen in the obliquity distribution of the sample of planets for which λ and i ⋆ have been measured, but no such peak is seen clearly when analyzing the larger sample of systems for which λ has been measured but not necessarily i ⋆ . We showed that the two samples might di ff er from each other because certain kinds of systems - namely, hot stars with hot Jupiters -are more amenable to i ⋆ measurements than others and also more likely to have near-polar orbits. The subsample of sub-Saturns around cooler stars also seems to have a preference for polar orbits, although the sample needs to be in-\n- to be sure. Thus, the resolution of this issue still seems likely to be astrophysically interesting.\n- ii) Is there a correlation between orbital eccentricity and stellar obliquity for close-orbiting giant planets? We concluded that the correlation is weaker than had been previously reported by Rice et al. (2022a). Furthermore, we argued that the correlation exists because both eccentricity and obliquity are correlated with a third variable, orbital separation. The causal mechanism is likely to be tidal dissipation, which lowers both eccentricity and obliquity at short orbital separations.\n- iii) How low are the obliquities of stars suspected of having been subject to tidal obliquity damping? Measurements of λ of cool stars with close-orbiting giant planets ( T e ff < 6100 K, a / R ⋆ < 10, M p > 0 . 3 MJ) have an intrinsic dispersion of 1 . 4 ± 0 . 7 · in λ , according to our Hierarchical Bayesian Model (HBM). Such fine alignment can be considered as supporting evidence for tidal obliquity damping, if the solar obliquity of 6 · is taken as representative of a typical un-damped 'primordial' obliquity. More precise measurements of λ in systems hosting massive planets around cool stars might even help to quantify the dependence of tidal dissipation rates on orbital separation.\n- iv) How well aligned are the central stars of compact systems of multiple transiting planets? Except for a few notable exceptions, the λ measurements in compact multiplanet systems are consistent with low obliquities. For this sample, we found an intrinsic dispersion in λ of 2 . 1 ± 1 . 6 · . In this sense, the multis are better aligned than the Solar System, indicating that compact multis are born unusually flat or that the Sun is unusually tilted, complicating the conclusions drawn above about tidal obliquity damping. The few examples of misaligned stars with compact systems of multiple transiting planets might be explained by interactions between the protoplanetary disk and a stellar companion, or by interactions with outer giant planets after the epoch of planet formation.  \nAcknowledgements. The authors thank the anonymous referee for helpful and insightful comments that greatly improved the quality of the manuscript. The authors would like to thank R. I. Dawson, J. Siegel, J. Dong, D. Foreman-Mackey, M. Rice, and S. Wang for helpful comments and discussions. Funding for the Stellar Astrophysics Centre is provided by The Danish National Research Foundation (Grant agreement no.: DNRF106). E.K. and S.H.A. acknowledge the support from the Danish Council for Independent Research through grant No.203200230B. Parts of the numerical results presented in this work were obtained at the Centre for Scientific Computing, Aarhus, https://phys.au.dk/forskning/ faciliteter/cscaa/ . This research has made use of data obtained from or tools provided by the portal https://exoplanet.eu/home/ of The Extrasolar Planets Encyclopaedia. This research has made use of NASA's Astrophysics Data System. This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. This work is based on observations made with the Very Large Telescope (VLT) at Cerro Paranal, Chile, the Telescopio Nazionale Galileo (TNG) at Roque de los Muchachos, La Palma, Spain, and the Nordic Optical Telescope (NOT) also at Roque de los Muchachos. NOT is owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofísica de Canarias. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. This paper includes data collected by the  \nKepler mission and obtained from the MAST data archive at the Space Telescope Science Institute (STScI). Funding for the Kepler mission is provided by the NASA Science Mission Directorate. STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. This paper includes data collected with the TESS mission, obtained from the MAST data archive at the Space Telescope Science Institute (STScI). Funding for the TESS mission is provided by the NASA Explorer Program. STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. The K2 light curves were extracted through the EVEREST pipeline (Luger et al. 2016, 2018) hosted on the Mikulski Archive for Space Telescopes (MAST). The TESS light curves were extracted using Lightkurve (Lightkurve Collaboration et al. 2018). EsoReflex (Freudling et al. 2013) was used to extract CCFs from the ESPRESSO data. This work made use of the following Python packages; NumPy (Harris et al. 2020), SciPy (Virtanen et al. 2020), Matplotlib (Hunter 2007), corner (Foreman-Mackey 2016), PyMC (Wiecki et al. 2023), and Astropy (Astropy Collaboration et al. 2013, 2018, 2022).", 'References': '- Ahlers, J. P., Kruse, E., Colón, K. D., et al. 2020, ApJ, 888, 63\n- Albrecht, S., Re ff ert, S., Snellen, I., Quirrenbach, A., & Mitchell, D. S. 2007,\n- A&A, 474, 565\n- Albrecht, S., Winn, J. N., Butler, R. P., et al. 2012a, ApJ, 744, 189\n- Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2011, ApJ, 738, 50\n- Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2012b, ApJ, 757, 18\n- Albrecht, S., Winn, J. N., Johnson, J. A., et al. 2012c, ApJ, 757, 18\n- Albrecht, S., Winn, J. N., Marcy, G. W., et al. 2013, ApJ, 771, 11\n- Albrecht, S. H., Dawson, R. I., & Winn, J. N. 2022, PASP, 134, 082001\n- Albrecht, S. H., Marcussen, M. L., Winn, J. N., Dawson, R. I., & Knudstrup, E. 2021, ApJ, 916, L1\n- Almenara, J. M., Hébrard, G., Díaz, R. F., et al. 2022, A&A, 663, A134\n- Alsubai, K., Tsvetanov, Z. I., Latham, D. W., et al. 2018, AJ, 155, 52\n- Anderson, D. R., Collier Cameron, A., Gillon, M., et al. 2011, A&A, 534, A16 Anderson, K. R., Winn, J. N., & Penev, K. 2021, ApJ, 914, 56\n- Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167\n- Astropy Collaboration, Price-Whelan, A. M., Sip"ocz, B. M., et al. 2018, AJ, 156, 123\n- Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33\n- Bakos, G. Á., Torres, G., Pál, A., et al. 2010, ApJ, 710, 1724\n- Barker, A. J. & Ogilvie, G. I. 2010, MNRAS, 404, 1849\n- Batygin, K. & Adams, F. C. 2013, ApJ, 778, 169\n- Beatty, T. G., Stevens, D. J., Collins, K. A., et al. 2017, AJ, 154, 25\n- Beck, J. G. & Giles, P. 2005, ApJ, 621, L153\n- Benz, W., Broeg, C., Fortier, A., et al. 2021, Experimental Astronomy, 51, 109\n- Best, S. & Petrovich, C. 2022, ApJ, 925, L5\n- Bieryla, A., Hartman, J. D., Bakos, G. Á., et al. 2014, AJ, 147, 84\n- Bonomo, A. S., Desidera, S., Benatti, S., et al. 2017a, A&A, 602, A107\n- Bonomo, A. S., Desidera, S., Benatti, S., et al. 2017b, A&A, 602, A107\n- Borsato, L., Piotto, G., Gandolfi, D., et al. 2021, MNRAS, 506, 3810\n- Bourrier, V., Attia, O., Mallonn, M., et al. 2023, A&A, 669, A63\n- Bourrier, V., Ehrenreich, D., Lendl, M., et al. 2020, A&A, 635, A205\n- Bourrier, V., Lovis, C., Cretignier, M., et al. 2021a, A&A, 654, A152\n- Bourrier, V., Lovis, C., Cretignier, M., et al. 2021b, A&A, 654, A152\n- Brahm, R., Espinoza, N., Rabus, M., et al. 2019, MNRAS, 483, 1970\n- Brahm, R., Nielsen, L. D., Wittenmyer, R. A., et al. 2020, AJ, 160, 235\n- Brown, D. J. A. 2014, MNRAS, 442, 1844\n- Brown, D. J. A., Anderson, D. R., Doyle, A. P., et al. 2014, arXiv e-prints, arXiv:1412.7761\n- Brown, D. J. A., Collier Cameron, A., Hall, C., Hebb, L., & Smalley, B. 2011, MNRAS, 415, 605\n- Brown, D. J. A., Triaud, A. H. M. J., Doyle, A. P., et al. 2017, MNRAS, 464, 810 Bruntt, H., Bedding, T. R., Quirion, P. O., et al. 2010, MNRAS, 405, 1907  \nButler, R. P., Marcy, G. W., Williams, E., et al. 1996, PASP, 108, 500  \n- Castelli, F. & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, Vol. 210, A20\n- Castro-González, A., Díez Alonso, E., Menéndez Blanco, J., et al. 2022, MNRAS, 509, 1075\n- Cegla, H. M., Lovis, C., Bourrier, V., et al. 2016a, A&A, 588, A127\n- Cegla, H. M., Lovis, C., Bourrier, V., et al. 2016b, A&A, 588, A127 Cegla, H. M., Roguet-Kern, N., Lendl, M., et al. 2023, A&A, 674, A174 Chaplin, W. J., Sanchis-Ojeda, R., Campante, T. L., et al. 2013, ApJ, 766, 101 Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580 Chontos, A., Huber, D., Grunblatt, S. K., et al. 2024, arXiv e-prints, arXiv:2402.07893\n- Christiansen, J. L., Vanderburg, A., Burt, J., et al. 2017, AJ, 154, 122  \n- Claret, A. 2018, A&A, 618, A20\n- Claret, A. 2021, Research Notes of the American Astronomical Society, 5, 13\n- Claret, A., Hauschildt, P. H., & Witte, S. 2012, A&A, 546, A14\n- Claret, A., Hauschildt, P. H., & Witte, S. 2013, A&A, 552, A16\n- Cosentino, R., Lovis, C., Pepe, F., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9147, Ground-based and Airborne Instrumentation for Astronomy V, ed. S. K. Ramsay, I. S. McLean, &H. Takami, 91478C\n- Crane, J. D., Shectman, S. A., Butler, R. P., et al. 2010, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 7735, Ground-based and Airborne Instrumentation for Astronomy III, ed. I. S. McLean, S. K. Ramsay, & H. Takami, 773553\n- Cristo, E., Esparza Borges, E., Santos, N. C., et al. 2024, A&A, 682, A28\n- Crossfield, I. J. M., Ciardi, D. R., Isaacson, H., et al. 2017, AJ, 153, 255 Crouzet, N., Healy, B. F., Hébrard, G., et al. 2020, AJ, 159, 44\n- Crouzet, N., McCullough, P. R., Long, D., et al. 2017, AJ, 153, 94  \nCzesla, S., Schneider, P. C., & Hatzes, A. 2023, A&A, 677, L12  \n- da Silva, R., Udry, S., Bouchy, F., et al. 2006, A&A, 446, 717\n- Dai, F., Masuda, K., Beard, C., et al. 2023, AJ, 165, 33\n- Dalal, S., Hébrard, G., Lecavelier des Étangs, A., et al. 2019, A&A, 631, A28\n- Damasso, M., Esposito, M., Nascimbeni, V., et al. 2015, A&A, 581, L6\n- Damiani, C. & Mathis, S. 2018, A&A, 618, A90\n- Dawson, R. I. 2014, ApJ, 790, L31\n- Dawson, R. I. & Johnson, J. A. 2018, ARA&A, 56, 175\n- Demangeon, O. D. S., Faedi, F., Hébrard, G., et al. 2018, A&A, 610, A63\n- Djupvik, A. A. & Andersen, J. 2010, in Astrophysics and Space Science Proceedings, Vol. 14, Highlights of Spanish Astrophysics V, 211\n- Dong, J. & Foreman-Mackey, D. 2023, AJ, 166, 112\n- Dorval, P., Talens, G. J. J., Otten, G. P. P. L., et al. 2020, A&A, 635, A60\n- Doyle, A. P., Davies, G. R., Smalley, B., Chaplin, W. J., & Elsworth, Y. 2014, MNRAS, 444, 3592\n- Doyle, L., Cegla, H. M., Anderson, D. R., et al. 2023, MNRAS, 522, 4499 Dravins, D. 1987, A&A, 172, 200\n- Dravins, D. 1999, in Astronomical Society of the Pacific Conference Series, Vol. 185, IAU Colloq. 170: Precise Stellar Radial Velocities, ed. J. B. Hearnshaw &C. D. Scarfe, 268\n- Dravins, D. & Nordlund, A. 1990a, A&A, 228, 184\n- Dravins, D. & Nordlund, A. 1990b, A&A, 228, 203\n- Dumusque, X., Cretignier, M., Sosnowska, D., et al. 2021, A&A, 648, A103\n- Eastman, J. D., Beatty, T. G., Siverd, R. J., et al. 2016, AJ, 151, 45\n- Ehrenreich, D., Lovis, C., Allart, R., et al. 2020, Nature, 580, 597  \nEisner, N. L., Barragán, O., Aigrain, S., et al. 2020, MNRAS, 494, 750  \n- Fabrycky, D. & Tremaine, S. 2007, ApJ, 669, 1298\n- Fabrycky, D. C., Ford, E. B., Ste ff en, J. H., et al. 2012, ApJ, 750, 114  \nForeman-Mackey, D. 2016, The Journal of Open Source Software, 1, 24  \n- Foreman-Mackey, D., Agol, E., Ambikasaran, S., & Angus, R. 2017, AJ, 154, 220\n- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306\n- Frandsen, S. & Lindberg, B. 1999, in Astrophysics with the NOT, ed. H. Karttunen & V. Piirola, 71\n- Frazier, R. C., Stefánsson, G., Mahadevan, S., et al. 2023, ApJ, 944, L41\n- Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96  \nGandolfi, D., Parviainen, H., Deeg, H. J., et al. 2015, A&A, 576, A11  \n- Garhart, E., Deming, D., Mandell, A., et al. 2020, AJ, 159, 137\n- Gaudi, B. S. & Winn, J. N. 2007, ApJ, 655, 550\n- Gillon, M., Doyle, A. P., Lendl, M., et al. 2011, A&A, 533, A88\n- Gomes, R., Deienno, R., & Morbidelli, A. 2017, AJ, 153, 27\n- Gratia, P. & Fabrycky, D. 2017, MNRAS, 464, 1709\n- Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres\n- Grou ff al, S., Santerne, A., Bourrier, V., et al. 2022a, A&A, 668, A172\n- Grou ff al, S., Santerne, A., Bourrier, V., et al. 2022b, A&A, 668, A172\n- Hagelberg, J., Nielsen, L. D., Attia, O., et al. 2023a, A&A, 679, A70\n- Hagelberg, J., Nielsen, L. D., Attia, O., et al. 2023b, A&A, 679, A70\n- Hansen, B. M. S. 2012, ApJ, 757, 6\n- Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357\n- Hébrard, G., Collier Cameron, A., Brown, D. J. A., et al. 2013, A&A, 549, A134\n- Hébrard, G., Díaz, R. F., Correia, A. C. M., et al. 2020, A&A, 640, A32\n- Hébrard, G., Ehrenreich, D., Bouchy, F., et al. 2011, A&A, 527, L11\n- Hellier, C., Anderson, D. R., Bouchy, F., et al. 2019a, MNRAS, 482, 1379\n- Hellier, C., Anderson, D. R., Triaud, A. H. M. J., et al. 2019b, MNRAS, 488, 3067  \nHirano, T., Suto, Y., Winn, J. N., et al. 2011, ApJ, 742, 69  \n- Hixenbaugh, K., Wang, X.-Y., Rice, M., & Wang, S. 2023, ApJ, 949, L35 Hjorth, M., Albrecht, S., Hirano, T., et al. 2021, Proceedings of the National Academy of Science, 118, e2017418118\n- Hjorth, M., Albrecht, S., Talens, G. J. J., et al. 2019, A&A, 631, A76\n- Howard, A. W., Johnson, J. A., Marcy, G. W., et al. 2010, ApJ, 721, 1467\n- Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398\n- Huang, C. X., Quinn, S. N., Vanderburg, A., et al. 2020, ApJ, 892, L7  \nHuber, D., Carter, J. A., Barbieri, M., et al. 2013, Science, 342, 331 Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90 Ikwut-Ukwa, M., Rodriguez, J. E., Bieryla, A., et al. 2020, AJ, 160, 209 Johnson, J. A., Winn, J. N., Bakos, G. Á., et al. 2011, ApJ, 735, 24 Johnson, M. C., Cochran, W. D., Albrecht, S., et al. 2014, ApJ, 790, 30 Johnson, M. C., Cochran, W. D., Collier Cameron, A., & Bayliss, D. 2015, ApJ, 810, L23 Johnson, M. C., Dai, F., Justesen, A. B., et al. 2018, MNRAS, 481, 596 Jordán, A., Brahm, R., Espinoza, N., et al. 2019, AJ, 157, 100 Kaluzny, J., Pych, W., Rucinski, S. M., & Thompson, I. B. 2006, Acta Astron., 56, 237 Kawai, Y., Narita, N., Fukui, A., Watanabe, N., & Inaba, S. 2024, MNRAS, 528, 270 Knudstrup, E. & Albrecht, S. H. 2022, A&A, 660, A99 Knudstrup, E., Albrecht, S. H., Gandolfi, D., et al. 2023a, A&A, 671, A164 Knudstrup, E., Albrecht, S. H., Gandolfi, D., et al. 2023b, A&A, 671, A164 Knudstrup, E., Serrano, L. M., Gandolfi, D., et al. 2022, A&A, 667, A22 Korth, J., Gandolfi, D., Šubjak, J., et al. 2023, A&A, 675, A115 Kreidberg, L. 2015, PASP, 127, 1161 Kuhn, R. B., Rodriguez, J. E., Collins, K. A., et al. 2016, MNRAS, 459, 4281 Labadie-Bartz, J., Rodriguez, J. E., Stassun, K. G., et al. 2019, ApJS, 240, 13 Lai, D. 2012, MNRAS, 423, 486 Lam, K. W. F., Faedi, F., Brown, D. J. A., et al. 2017, A&A, 599, A3 Lavie, B., Bouchy, F., Lovis, C., et al. 2023, A&A, 673, A69 Li, G. & Winn, J. N. 2016, ApJ, 818, 5 Libby-Roberts, J. E., Schutte, M., Hebb, L., et al. 2023, AJ, 165, 249 Lightkurve Collaboration, Cardoso, J. V. d. M., Hedges, C., et al. 2018, Lightkurve: Kepler and TESS time series analysis in Python, Astrophysics Source Code Library Lin, D. N. C., Bodenheimer, P., & Richardson, D. C. 1996, Nature, 380, 606 Lin, Y. & Ogilvie, G. I. 2017, MNRAS, 468, 1387 Lovis, C. & Pepe, F. 2007, A&A, 468, 1115 Lubin, J., Wang, X.-Y., Rice, M., et al. 2023, ApJ, 959, L5 Luger, R., Agol, E., Kruse, E., et al. 2016, AJ, 152, 100 Luger, R., Kruse, E., Foreman-Mackey, D., Agol, E., & Saunders, N. 2018, AJ, 156, 99 Lund, M. B., Rodriguez, J. E., Zhou, G., et al. 2017, AJ, 154, 194 Lund, M. N., Knudstrup, E., Silva Aguirre, V., et al. 2019, AJ, 158, 248 Ma, L. & Fuller, J. 2021, ApJ, 918, 16 Mancini, L., Esposito, M., Covino, E., et al. 2018, A&A, 613, A41 Mancini, L., Hartman, J. D., Penev, K., et al. 2015, A&A, 580, A63 Mandel, K. & Agol, E. 2002, ApJ, 580, L171 Masuda, K. & Winn, J. N. 2020, AJ, 159, 81 McCullough, P. R., Stys, J. E., Valenti, J. A., et al. 2005, PASP, 117, 783 McQuillan, A., Mazeh, T., & Aigrain, S. 2014, ApJS, 211, 24 Morgan, M., Bowler, B. P., Tran, Q. H., et al. 2023, arXiv e-prints, arXiv:2310.18445 Mounzer, D., Lovis, C., Seidel, J. V., et al. 2022, A&A, 668, A1 Moˇcnik, T., Clark, B. J. M., Anderson, D. R., Hellier, C., & Brown, D. J. A. 2016, AJ, 151, 150 Moya, A., Bouy, H., Marchis, F., Vicente, B., & Barrado, D. 2011, A&A, 535, A110 Nagasawa, M., Ida, S., & Bessho, T. 2008, ApJ, 678, 498 Narita, N., Fukui, A., Kusakabe, N., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 045001 Narita, N., Hirano, T., Fukui, A., et al. 2017, PASJ, 69, 29 Neveu-VanMalle, M., Queloz, D., Anderson, D. R., et al. 2014, A&A, 572, A49 Newton, E. R., Mann, A. W., Kraus, A. L., et al. 2021, AJ, 161, 65 Ogilvie, G. I. 2014, ARA&A, 52, 171 Osborn, H. P., Bonfanti, A., Gandolfi, D., et al. 2022, A&A, 664, A156 Penev, K., Bouma, L. G., Winn, J. N., & Hartman, J. D. 2018, AJ, 155, 165 Pepe, F., Cristiani, S., Rebolo, R., et al. 2021, A&A, 645, A96 Pepper, J., Kane, S. R., Rodriguez, J. E., et al. 2020, AJ, 159, 243 Pepper, J., Siverd, R. J., Beatty, T. G., et al. 2013, ApJ, 773, 64 Petrovich, C., Muñoz, D. J., Kratter, K. M., & Malhotra, R. 2020, ApJ, 902, L5 Piaulet, C., Benneke, B., Rubenzahl, R. A., et al. 2021, AJ, 161, 70 Queloz, D., Eggenberger, A., Mayor, M., et al. 2000, A&A, 359, L13 Rice, M., Wang, S., Gerbig, K., et al. 2023a, AJ, 165, 65 Rice, M., Wang, S., & Laughlin, G. 2022a, ApJ, 926, L17 Rice, M., Wang, S., Wang, X.-Y., et al. 2022b, AJ, 164, 104 Rice, M., Wang, X.-Y., Wang, S., et al. 2023b, AJ, 166, 266 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 014003 Rodriguez, J. E., Quinn, S. N., Zhou, G., et al. 2021, AJ, 161, 194 Rosenthal, L. J., Fulton, B. J., Hirsch, L. A., et al. 2021, ApJS, 255, 8 Rucinski, S. 1999, in Astronomical Society of the Pacific Conference Series, Vol. 185, IAU Colloq. 170: Precise Stellar Radial Velocities, ed. J. B. Hearnshaw &C. D. Scarfe, 82 Saha, S., Chakrabarty, A., & Sengupta, S. 2021, AJ, 162, 18  \nSanchis-Ojeda, R., Fabrycky, D. C., Winn, J. N., et al. 2012, Nature, 487, 449  \n- Sanchis-Ojeda, R. & Winn, J. N. 2011, ApJ, 743, 61\n- Sanchis-Ojeda, R., Winn, J. N., Dai, F., et al. 2015, ApJ, 812, L11\n- Sanchis-Ojeda, R., Winn, J. N., Marcy, G. W., et al. 2013a, ApJ, 775, 54\n- Sanchis-Ojeda, R., Winn, J. N., Marcy, G. W., et al. 2013b, ApJ, 775, 54\n- Santerne, A., Hébrard, G., Lillo-Box, J., et al. 2016, ApJ, 824, 55\n- Santerne, A., Malavolta, L., Kosiarek, M. R., et al. 2019, arXiv e-prints, arXiv:1911.07355\n- Schanche, N., Hébrard, G., Collier Cameron, A., et al. 2020, MNRAS, 499, 428 Schlaufman, K. C. 2010, ApJ, 719, 602\n- Schwarz, G. 1978, The Annals of Statistics, 6, 461\n- Sebastian, D., Guenther, E. W., Deleuil, M., et al. 2022, MNRAS, 516, 636\n- Sedaghati, E., Jordán, A., Brahm, R., et al. 2023, AJ, 166, 130\n- Seidel, J. V., Prinoth, B., Knudstrup, E., et al. 2023, A&A, 678, A150\n- Shporer, A. & Brown, T. 2011, ApJ, 733, 30\n- Siegel, J. C., Winn, J. N., & Albrecht, S. H. 2023, ApJ, 950, L2\n- Singh, V., Scandariato, G., Smith, A. M. S., et al. 2023, arXiv e-prints, arXiv:2311.03264\n- Siverd, R. J., Collins, K. A., Zhou, G., et al. 2018, AJ, 155, 35\n- Smalley, B., Anderson, D. R., Collier Cameron, A., et al. 2010, A&A, 520, A56 Souami, D. & Souchay, J. 2012, A&A, 543, A133\n- Southworth, J. 2011, MNRAS, 417, 2166\n- Southworth, J., Dominik, M., Jørgensen, U. G., et al. 2011, A&A, 527, A8  \nSteiner, M., Attia, O., Ehrenreich, D., et al. 2023, A&A, 672, A134  \n- Szentgyorgyi, A. H., Geary, J. G., Latham, D. W., et al. 2005, in American Astronomical Society Meeting Abstracts, Vol. 207, American Astronomical Society Meeting Abstracts, 110.10\n- Tejada Arevalo, R. A., Winn, J. N., & Anderson, K. R. 2021, ApJ, 919, 138 Telting, J. H., Avila, G., Buchhave, L., et al. 2014, Astronomische Nachrichten, 335, 41\n- Temple, L. Y., Hellier, C., Albrow, M. D., et al. 2017, MNRAS, 471, 2743 Valsecchi, F. & Rasio, F. A. 2014, ApJ, 786, 102\n- Vanderburg, A., Becker, J. C., Buchhave, L. A., et al. 2017, AJ, 154, 237  \nVirtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261  \n- Vogt, S. S., Allen, S. L., Bigelow, B. C., et al. 1994, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 2198, Instrumentation in Astronomy VIII, ed. D. L. Crawford & E. R. Craine, 362\n- von Braun, K., Boyajian, T. S., ten Brummelaar, T. A., et al. 2011, ApJ, 740, 49 Wang, X.-Y., Rice, M., Wang, S., et al. 2022, ApJ, 926, L8\n- Wang, Y.-H., Wang, S., Liu, H.-G., et al. 2017, AJ, 154, 49\n- Ward, W. R. 1981, Icarus, 47, 234\n- Wiecki, T., Salvatier, J., Vieira, R., et al. 2023, pymc-devs / pymc: v5.8.2\n- Winn, J. N., Fabrycky, D., Albrecht, S., & Johnson, J. A. 2010a, ApJ, 718, L145\n- Winn, J. N., Johnson, J. A., Howard, A. W., et al. 2010b, ApJ, 723, L223\n- Winn, J. N., Noyes, R. W., Holman, M. J., et al. 2005, ApJ, 631, 1215\n- Winters, J. G., Cloutier, R., Medina, A. A., et al. 2022, AJ, 163, 168\n- Winters, J. G., Medina, A. A., Irwin, J. M., et al. 2019, AJ, 158, 152\n- Wright, J., Rice, M., Wang, X.-Y., Hixenbaugh, K., & Wang, S. 2023, AJ, 166, 217\n- Yee, S. W., Winn, J. N., Hartman, J. D., et al. 2023, ApJS, 265, 1\n- Zahn, J. P. 1977, A&A, 57, 383\n- Zanazzi, J. J. & Chiang, E. 2024, MNRAS, 527, 7203\n- Zanazzi, J. J., Dewberry, J., & Chiang, E. 2024, arXiv e-prints, arXiv:2403.05616  \nZanazzi, J. J. & Lai, D. 2018, MNRAS, 477, 5207  \n- Zechmeister, M., Reiners, A., Amado, P. J., et al. 2018, A&A, 609, A12  \nZhao, L. L., Kunovac, V., Brewer, J. M., et al. 2023, Nature Astronomy, 7, 198  \n- Zhou, G., Bakos, G. Á., Bayliss, D., et al. 2019, AJ, 157, 31\n- Zhou, G., Rodriguez, J. E., Vanderburg, A., et al. 2018, AJ, 156, 93  \nFig. 38. Obliquities of stars with multiple transiting planets, as constrained by either the RM e ff ect or asteroseismology. Each horizontal line shows the a / R ⋆ values of the planets in a given system, with the solar system on top (with gaps to allow Jupiter and Neptune to be shown). For the planets, the size of each circle conveys the planet\'s radius relative to Jupiter. For the stars, symbol size conveys the radius relative to the Sun, with a color conveying the e ff ective temperature. Gray circles are planets for which the RM e ff ect has been measured - except for β Pic b where spectro-interferometry was used. The wedge indicates the projected obliquity for a given system, based on the most precise measurement for any planet in the system, with the exception of HD 3167 for which two planets appear to have very di ff erent inclinations. The asteroseismic measurements constrain i ⋆ and not λ ; in those cases the wedge should therefore be read as aligned ( Kepler -50 and Kepler -65) or misaligned ( Kepler -56, Kepler -129). Planetary (tan circle) or stellar companions (red star, M-dwarfs) to a given system are shown to the right, where √ ( ÷ ) denotes that the companion could (not) have influenced the obliquity of the inner planetary system. The asterisks for LTT 1445A and TOI-451A denote that these measurements should be taken with a grain of salt. Adapted from Wang et al. (2022, Figure 3).  \n<!-- image -->  \n/star', 'Appendix A: Additional figures': 'Knudstrup et al.: Obliquities of Exoplanet Host Stars  \nFig. A.1. Phasefolded light curves. The best-fitting models are shown as the white line in the top panel and the residuals are given in the bottom panel. (a) HD 118203 b: TESS 2 min. cadence data shown as gray points with error bars. Black markers are binned data in ∼ 10 min. intervals. (b) HD 149193 b: TESS light curves 20 sec., 2 min., and 30 min. cadence data are shown in light gray, gray, and black, respectively. (c) K2-261 b: K2 light curve in black, TESS 20 sec. and 2 min. in light gray and gray, respectively. A light curve for each passband is shown. (d) K2-287 b: K2 30 min cadence shown in orange with CHEOPS 1 min. and ground-based 2 min. observations in light gray and gray, respectively. A light curve for each passband is also shown with a matching color. (e) KELT-3 b: TESS 2 min. cadence data shown in gray. (f) KELT-4Ab: TESS 2 min. cadence data shown in gray.  \n<!-- image -->  \nFig. A.2. Phasefolded light curves. Same as in Fig. A.1. (a) LTT 1445Ab: Black TESS 2 min. and gray TESS 20 sec. (b) TOI-451Ab: Unbinned 2 min. cadence data are shown in gray and in black the binned ( ∼ 6 min.) data are shown. (c) TOI-813 b: TESS 20 sec., 2 min., and 30 min. cadence data are shown in light gray, gray, and black, respectively. (d) TOI-892 b: TESS 2 min. and 30 min. cadence data are shown in gray and black, respectively. (e) TOI-1130 c: TESS 20 sec. and 30 min. cadence data shown in gray and black, respectively. (f) WASP-50 b: TESS 20 sec. and 2 min cadence data shown in gray and black, respectively.  \n<!-- image -->  \nFig. A.3. Phasefolded light curves. Same as in Fig. A.1. (a) WASP-59 b: TESS 2 min. cadence data shown in gray and ground-based 2 min. KeplerCam photometry in orange. (b) WASP-136 b: TESS 20 sec. and 2 min. cadence data shown in gray and black, respectively. (c) WASP172 b: TESS 20 sec. and 30 min. cadence data shown in gray and black, respectively. (d) WASP-173Ab: TESS 20 sec. and 2 min. cadence data shown in gray and black, respectively. (e) WASP-186 b: TESS 20 sec., 2 min., and 30 min. cadence data are shown in light gray, gray, and black, respectively. (f) XO-7 b: TESS 20 sec., 2 min., and 30 min. cadence data are shown in light gray, gray, and black, respectively.  \n<!-- image -->  \n<!-- image -->', '(a) WASP-59 b': '(b) WASP-136 b  \n<!-- image -->  \n(c) WASP-172 b  \n<!-- image -->  \n(d) WASP-173Ab  \n<!-- image -->  \n<!-- image -->  \n(e) WASP-186 b  \nFig. A.4. Phasefolded light curves of WASP-148 b. The MuSCAT-2 photometry obtained simultaneous with our spectroscopic transit observations. The observations in the di ff erent griz filters are shown as blue, orange, green, and red markers, respectively.  \n<!-- image -->  \nFig. A.5. Phasefolded light curve of WASP-26 b. Unbinned 20 sec. cadence data are shown in gray and in black the binned ( ∼ 6 min.) data are shown.  \n<!-- image -->  \nTable B.1. MCMCposteriors.  \nFor WASP-148 b we were stepping in i , where we got i = 86 . 57 + 0 . 16 · -0 . 16 .', 'Appendix B: Additional tables': 'Table B.2. Extension to Table A1 of Albrecht et al. (2022).  \nKnudstrup et al.: Obliquities of Exoplanet Host Stars  \nTable B.2 - continued from previous page  \n.  \nValues for the new literature systems are drawn from the following references (with bold denoting the source for λ ): 55 Cnc e: 1 Rosenthal et al. (2021), 2 Zhao et al. (2023). CoRoT-36 b: 1 Sebastian et al. (2022). HAT-P-11 b: 1 Sanchis-Ojeda & Winn (2011), 2 Winn et al. (2010b). HAT-P-3 b: 1 Bonomo et al. (2017b) , 2 Bourrier et al. (2023), 3 Mancini et al. (2018). HAT-P-30 b: 1 Bonomo et al. (2017b), 2 Cegla et al. (2023). HAT-P-32 b: 1 Bonomo et al. (2017b), 2 Czesla et al. (2023), 3 Albrecht et al. (2012c). HAT-P-33b: 1 Wang et al. (2017), 2 Bourrier et al. (2023). HAT-P-49 b: 1 Bieryla et al. (2014), 2 Bourrier et al. (2023). HD 3167 c: 1 Bourrier et al. (2021b). HD 106315 c: 1 Zhou et al. (2018), 2 Bourrier et al. (2023). HD 189733 b: 1 Bonomo et al. (2017b), 2 Cristo et al. (2024). HIP 41378 d: 1 Lund et al. (2019), 2 Santerne et al. (2019), 3 Grou ff al et al. (2022b). K2-29 b: 1 Santerne et al. (2016). K2-105 b: 1 Narita et al. (2017), 2 Castro-González et al. (2022), 3 Bourrier et al. (2023). KELT-10 b: 1 Kuhn et al. (2016), 2 Steiner et al. (2023). KELT-11 b: 1 Beatty et al. (2017), 2 Mounzer et al. (2022). KELT-19 b: 1 Kawai et al. (2024), 2 Siverd et al. (2018). Kepler-63 b: 1 Sanchis-Ojeda et al. (2013b). MASCARA-2 b: 1 Singh et al. (2023). MASCARA-4 b: 1 Ahlers et al. (2020), 2 Dorval et al. (2020). Qatar-6 b: 1 Rice et al. (2023a). TOI-640 b: 1 Knudstrup et al. (2023b). TOI-677 b: 1 Sedaghati et al. (2023). TOI-1670 c: 1 Lubin et al. (2023). TOI-858Bb: 1 Hagelberg et al. (2023b). TOI-1136 d: 1 Dai et al. (2023). TOI-1478 b: 1 Rodriguez et al. (2021), 2 Rice et al. (2022b). TOI-1842 b: 1 Hixenbaugh et al. (2023). TOI-1937 b: 1 Yee et al. (2023). TOI-2025 b: 1 Knudstrup et al. (2022). TOI-2076 b: 1 Osborn et al. (2022), 2 Frazier et al. (2023). TOI-2202 b: 1 Rice et al. (2023b). TOI-3884 b: 1 Libby-Roberts et al. (2023). WASP-7 b: 1 Bonomo et al. (2017b), 2 Albrecht et al. (2012a). WASP-12 b: 1 Bonomo et al. (2017b), 2 Albrecht et al. (2012c). WASP-33 b: 1 Johnson et al. (2015). WASP-47 b: 1 Bourrier et al. (2023), 2 Sanchis-Ojeda et al. (2015). WASP-52 b: 1 Bonomo et al. (2017b), 2 Cegla et al. (2023). WASP-62 b: 1 Garhart et al. (2020), 2 Brown et al. (2017). WASP-76 b: 1 Ehrenreich et al. (2020). WASP-85 b: 1 Moˇcnik et al. (2016). WASP-94Ab: 1 NeveuVanMalle et al. (2014). WASP-106 b: 1 Wright et al. (2023). WASP-107 b: 1 Piaulet et al. (2021), 2 Bourrier et al. (2023). WASP-121 b: 1 Bourrier et al. (2020). WASP-131 b: 1 Doyle et al. (2023). WASP-156 b: 1 Demangeon et al. (2018), 2 Saha et al. (2021), 3 Bourrier et al. (2023). WASP-166 b: 1 Bourrier et al. (2023). WASP-167 b: 1 Temple et al. (2017). XO-6 b: 1 Crouzet et al. (2017).  \nTable B.3. Extension to Table A2 of Albrecht et al. (2022).  \nKnudstrup et al.: Obliquities of Exoplanet Host Stars  \nTable B.3 - continued from previous page  \nValues for the new literature systems are drawn from the following references: 55 Cnc: 1 von Braun et al. (2011), 2 Zhao et al. (2023). CoRoT-36: 1 Sebastian et al. (2022). HAT-P-11: 1 Bakos et al. (2010), 2 Winn et al. (2010b). HAT-P-3: 1 Bourrier et al. (2023), 2 Mancini et al. (2018). HAT-P-30: 1 Johnson et al. (2011), 2 Cegla et al. (2023). HAT-P-32: 1 Bonomo et al. (2017b), 2 Albrecht et al. (2012c). HAT-P-33: 1 Wang et al. (2017), 2 Bourrier et al. (2023). HAT-P-49: 1 Bieryla et al. (2014), 2 Bourrier et al. (2023). HD 3167: 1 Christiansen et al. (2017), 2 Dalal et al. (2019), 3 Bourrier et al. (2021b). HD 106315: 1 Crossfield et al. (2017), 2 Bourrier et al. (2023). HD 189733: 1 Cegla et al. (2016b). HIP 41378: 1 Lund et al. (2019), 2 Grou ff al et al. (2022b). K2-29b: 1 Santerne et al. (2016). K2-105: 1 Castro-González et al. (2022), 2 Bourrier et al. (2023). KELT10: 1 Kuhn et al. (2016), 2 Steiner et al. (2023). KELT-11: 1 Beatty et al. (2017), 2 Mounzer et al. (2022). KELT-19: 1 Siverd et al. (2018). Kepler-63: 1 Sanchis-Ojeda et al. (2013b). MASCARA-2: 1 Lund et al. (2017), 2 Singh et al. (2023). MASCARA-4: 1 Dorval et al. (2020). Qatar-6: 1 Alsubai et al. (2018), 2 Rice et al. (2023a). TOI-640: 1 Rodriguez et al. (2021), 2 Knudstrup et al. (2023b). TOI-677: 1 Sedaghati et al. (2023). TOI-1670: 1 Lubin et al. (2023). TOI-858B: 1 Hagelberg et al. (2023b). TOI-1136: 1 Dai et al. (2023). TOI-1478: 1 Rodriguez et al. (2021), 2 Rice et al. (2022b). TOI-1842: 1 Hixenbaugh et al. (2023). TOI-1937: 1 Yee et al. (2023). TOI-2025: 1 Knudstrup et al. (2022). TOI-2076: 1 Frazier et al. (2023). TOI-2202: 1 Rice et al. (2023b). TOI-3884: 1 LibbyRoberts et al. (2023). WASP-7: 1 Southworth et al. (2011), 2 Albrecht et al. (2012a). WASP-12: 1 Bonomo et al. (2017b), 2 Albrecht et al. (2012c). WASP-33: 1 Moya et al. (2011), 2 Johnson et al. (2015). WASP-47: 1 Vanderburg et al. (2017), 2 Sanchis-Ojeda et al. (2015). WASP-52: 1 Bonomo et al. (2017b), 2 Cegla et al. (2023). WASP-62: 1 Brown et al. (2017). WASP-76: 1 Ehrenreich et al. (2020). WASP-85: 1 Brown et al. (2014), 2 Moˇcnik et al. (2016). WASP-94A: 1 Bonomo et al. (2017b), 2 Neveu-VanMalle et al. (2014). WASP-106: 1 Wright et al. (2023). WASP-107: 1 Piaulet et al. (2021), 2 Bourrier et al. (2023). WASP-121: 1 Bourrier et al. (2020). WASP-131: 1 Doyle et al. (2023). WASP-156: 1 Demangeon et al. (2018), 2 Bourrier et al. (2023). WASP-166: 1 Hellier et al. (2019b), 2 Bourrier et al. (2023). WASP-167: 1 Temple et al. (2017). XO-6: 1 Crouzet et al. (2017).  \nTable B.4. Hyperparameters for the beta-distribution.  \nSubscript 0 denotes the parameters for the misaligned component, while subscript 1 is for the aligned component. The true mean of cos ψ is given as 2 µ -1. a ) Results from Dong & Foreman-Mackey (2023). A Gaussian prior was applied to log κ ( N (0 , 3) as in Dong & Foreman-Mackey 2023) for the runs above the dashed line, whereas a uniform prior of U ( -4 , 10) for the runs below.  \nKnudstrup et al.: Obliquities of Exoplanet Host Stars  \nTable B.5. Radial velocity measurements for all targets.  \nThis table is available online in its entirety in machine-readable form.'}