We analyze the high-resolution emission spectrum of WASP-33b taken using the High Dispersion Spectrograph (R\,$\approx$\,165,000) on the 8.2-m Subaru telescope. The data cover $\lambda$\,$\approx$\,$6170$-$8817$\,\AA, divided over 30 spectral orders. The telluric and stellar lines are removed using a de-trending algorithm, {\sc SysRem}, before cross-correlating with planetary spectral templates. We calculate the templates assuming a 1-D plane-parallel hydrostatic atmosphere including continuum opacity of bound-free H$^{-}$ and Rayleigh scattering by H$_{2}$ with a range of constant abundances of Fe\,{\sc i}. Using a likelihood-mapping analysis, we detect an Fe\,{\sc i} emission signature at 6.4-$\sigma$ located at $K_{\mathrm{p}}$ of 226.0\,$^{+2.1}_{-2.3}$\,km\,s$^{-1}$and $v_{\mathrm{sys}}$ of -3.2\,$^{+2.1}_{-1.8}$\,km\,s$^{-1}$ -- consistent with the planet's expected velocity in the literature. We also confirm the existence of a thermal inversion in the day-side of the planet which is very likely to be caused by the presence of Fe\,{\sc i} and previously-detected TiO in the atmosphere. This makes WASP-33b one of the prime targets to study the relative contributions of both species to the energy budget of an ultra-hot Jupiter.

With a temperature akin to an M-dwarf, WASP-33b is among the hottest Jupiters known, making it an ideal target for high-resolution optical spectroscopy. By analyzing both transmission and emission spectra, we aim to substantiate previous reports of atmospheric TiO and a thermal inversion within the planet's atmosphere. We observed two transits and six arcs of the phase curve with ESPaDOns on the Canada-France-Hawaii Telescope and HIRES on the Keck telescope, which provide high spectral resolution and ample wavelength coverage. We employ the Doppler cross-correlation technique to search for the molecular signatures of TiO and H$_2$O in these spectra, using models based on the TiO line list of Plez (2012). Though we cannot exclude line-list-dependent effects, our data do not corroborate previous indications of a thermal inversion. Instead we place a $3\sigma$ upper limit of $10^{-9}$ on the volume mixing ratio of TiO for the T-P profile we consider. While we are unable to constrain the volume mixing ratio of water, our strongest constraint on TiO comes from day-side emission spectra. This apparent absence of a stratosphere sits in stark contrast to previous observations of WASP-33b as well as theoretical predictions for the atmospheres of highly irradiated planets. The discrepancy could be due to variances between line lists, and we stress that detection limits are only as good as the line list employed, and are only valid for the specific T-P profile considered due to the strong degeneracy between lapse rate ($dT/d\log P$) and molecular abundance.

Currently, we have only limited means to probe the presence of planets at large orbital separations. Foreman-Mackey et al. searched for long-period transiting planets in the Kepler light curves using an automated pipeline. Here, we apply their pipeline, with minor modifications, to a larger sample and use updated stellar parameters from Gaia DR2. The latter boosts the stellar radii for most of the planet candidates found by FM16, invalidating a number of them as false positives. We identify 15 candidates, including two new ones. All have sizes from 0.3 to 1 $R_{\rm J}$, and all but two have periods from 2 to 10 yr. We report two main findings based on this sample. First, the planet occurrence rate for the above size and period ranges is $0.70^{+0.40}_{-0.20}$ planets per Sun-like star, with the frequency of cold Jupiters agreeing with that from radial velocity surveys. Planet occurrence rises with decreasing planet size, roughly describable as $dN/d\log R \propto R^{\alpha}$ with $\alpha = -1.6^{+1.0}_{-0.9}$, i.e., Neptune-sized planets are some four times more common than Jupiter-sized ones. Second, five out of our 15 candidates orbit stars with known transiting planets at shorter periods, including one with five inner planets. We interpret this high incidence rate to mean: (1) almost all our candidates should be genuine; (2) across a large orbital range (from $\sim 0.05$ to a few astronomical units), mutual inclinations in these systems are at most a few degrees; and (3) large outer planets exist almost exclusively in systems with small inner planets.

Gravity darkening induced by rapid stellar rotation provides us with a unique opportunity to characterize the spin-orbit misalignment of a planetary system through analysis of its photometric transit. We use the gravity-darkened transit modeling code simuTrans to reproduce the transit light curve of Kepler-13Ab by separately analyzing phase-folded transits for 12 short-cadence Kepler quarters. We verify the temporal change in impact parameter indicative of spin-orbit precession identified by Szab\'o et al. (2012) and Masuda (2015), reporting a rate of change $db/dt = (-4.1 \pm 0.2) \times 10^{-5}$ day$^{-1}$. We further investigate the effect of light dilution on the fitted impact parameter and find that less than 1% of additional light is sufficient to explain the seasonal variation seen in the Kepler quarter data. We then extend our precession analysis to the phase curve data from which we report a rate of change $db/dt = (-3.2 \pm 1.3) \times 10^{-5}$ day$^{-1}$. This value is consistent with that of the transit data at a lower significance and provides the first evidence of spin-orbit precession based solely on the temporal variation of the secondary eclipse.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.