We present the detection and analysis of the phase curve of KELT-1b at optical wavelengths, analyzing data taken by the Transiting Exoplanet Survey Satellite (TESS). With a mass of ~27 M_J, KELT-1b is a low-mass brown dwarf. Due to the high mass and close proximity of its companion, the host star has a TESS light curve which shows clear ellipsoidal variations. We model the data with a six-component model: secondary eclipse, phase curve accounting for reflected light and thermal emission, Doppler beaming, ellipsoidal variations, stellar activity and the primary transit. We determine the secondary eclipse depth in the TESS bandpass to be 304 +/- 75 parts-per-million (ppm), the most accurate eclipse depth determined so far for KELT-1b. We measure the amplitude of the phase curve to be 128 +/- 27 ppm, with a corresponding eastward offset between the region of maximum brightness and the substellar point of 19.2 +/- 9.6 degrees, in good agreement with Spitzer measurements. We determine day and night brightness temperatures of 3201 +/- 147 K and 1484 +/- 110 K, respectively, slightly higher than those from Spitzer 3.6 and 4.5 micrometer data. A one-dimensional self-consistent atmospheric model can explain the TESS and Spitzer day side brightness temperatures with thermal emission alone and no reflected light. The night side data can be explained by a model with an internal temperature of ~1100 K, which may be related to the inflated radius. The difference between the TESS and Spitzer brightness temperatures can be explained by stronger molecular opacity in the Spitzer bands. On the night side, this opacity is due primarily to CH4 and CO while on the day side it is due to H2-H2 and H2-He collision-induced absorption.

We present the detection and characterization of the full-orbit phase curve and secondary eclipse of the ultra-hot Jupiter WASP-33b at optical wavelengths, along with the pulsation spectrum of the host star. We analyzed data collected by the Transiting Exoplanet Survey Satellite (TESS) in sector 18. WASP-33b belongs to a very short list of highly irradiated exoplanets that were discovered from the ground and were later visited by TESS. The host star of WASP-33b is of delta Scuti-type and shows nonradial pulsations in the millimagnitude regime, with periods comparable to the period of the primary transit. These completely deform the photometric light curve, which hinders our interpretations. By carrying out a detailed determination of the pulsation spectrum of the host star, we find 29 pulsation frequencies with a signal-to-noise ratio higher than 4. After cleaning the light curve from the stellar pulsations, we confidently report a secondary eclipse depth of 305.8 +/- 35.5 parts-per-million (ppm), along with an amplitude of the phase curve of 100.4 +/- 13.1 ppm and a corresponding westward offset between the region of maximum brightness and the substellar point of 28.7 +/- 7.1 degrees, making WASP-33b one of the few planets with such an offset found so far. Our derived Bond albedo, A_B = 0.369 +/- 0.050, and heat recirculation efficiency, epsilon = 0.189 +/- 0.014, confirm again that he behavior of WASP-33b is similar to that of other hot Jupiters, despite the high irradiation received from its host star. By connecting the amplitude of the phase curve to the primary transit and depths of the secondary eclipse, we determine that the day- and nightside brightness temperatures of WASP-33b are 3014 +/- 60 K and 1605 +/- 45 K, respectively. From the detection of photometric variations due to gravitational interactions, we estimate a planet mass of M_P = 2.81 +/- 0.53 M$_J.

HD 209458b was the first transiting planet discovered, and the first for which an atmosphere, in particular NaI, was detected. With time, it has become one of the most frequently studied planets, with a large diversity of atmospheric studies using low- and high-resolution spectroscopy. Here, we present transit spectroscopy observations of HD 209458b using the HARPS-N and CARMENES spectrographs. We fit the Rossiter-McLaughlin effect by combining radial velocity data from both instruments (nine transits in total), measuring a projected spin-orbit angle of $-1.6\pm0.3~{\rm deg}$. We also present the analysis of high-resolution transmission spectroscopy around the NaI region at $590~{\rm nm}$, using a total of five transit observations. In contrast to previous studies where atmospheric NaI absorption is detected, we find that for all of the nights, whether individually or combined, the transmission spectra can be explained by the combination of the centre-to-limb variation and the Rossiter-McLaughlin effect. This is also observed in the time-evolution maps and transmission light curves, but at lower signal-to-noise ratio. Other strong lines such as H$\alpha$, CaII IRT, the MgI triplet region, and KI D1 are analysed, and are also consistent with the modelled effects, without considering any contribution from the exoplanet atmosphere. Thus, the transmission spectrum reveals no detectableNaI absorption in HD 209458b. We discuss how previous pioneering studies of this benchmark object may have overlooked these effects. While for some star-planet systems these effects are small, for other planetary atmospheres the results reported in the literature may require revision.

Ultra-hot Jupiters are emerging as a new class of exoplanets. Studying their chemical compositions and temperature structures will improve the understanding of their mass loss rate as well as their formation and evolution. We present the detection of ionized calcium in the two hottest giant exoplanets - KELT-9b and WASP-33b. By utilizing transit datasets from CARMENES and HARPS-N observations, we achieved high confidence level detections of Ca II using the cross-correlation method. We further obtain the transmission spectra around the individual lines of the Ca II H&K doublet and the near-infrared triplet, and measure their line profiles. The Ca II H&K lines have an average line depth of 2.02 $\pm$ 0.17 % (effective radius of 1.56 Rp) for WASP-33b and an average line depth of 0.78 $\pm$ 0.04 % (effective radius of 1.47 Rp) for KELT-9b, which indicates that the absorptions are from very high upper atmosphere layers close to the planetary Roche lobes. The observed Ca II lines are significantly deeper than the predicted values from the hydrostatic models. Such a discrepancy is probably a result of hydrodynamic outflow that transports a significant amount of Ca II into the upper atmosphere. The prominent Ca II detection with the lack of significant Ca I detection implies that calcium is mostly ionized in the upper atmospheres of the two planets.

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

High-resolution transmission spectroscopy is a method for understanding the chemical and physical properties of upper exoplanetary atmospheres. Due to large absorption cross-sections, resonance lines of atomic sodium D-lines (at 5889.95 $\AA$ and 5895.92 $\AA$) produce large transmission signals. Our aim is to unveil the physical properties of WASP-17b through an accurate measurement of the sodium absorption in the transmission spectrum. We analyze 37 high-resolution spectra observed during a single transit of WASP-17b with the MIKE instrument on the 6.5 meter Magellan Telescopes. We exclude stellar flaring activity during the observations by analyzing the temporal variations of H$_{\alpha}$ and Ca II infra-red triplet (IRT) lines. Then we obtain the excess absorption light curves in wavelength bands of 0.75, 1, 1.5 and 3 $\AA$ around the center of each sodium line (i.e., the light curve approach). We model the effects of differential limb-darkening, and the changing planetary radial velocity on the light curves. We also analyze the sodium absorption directly in the transmission spectrum, which is obtained through dividing in-transit by out-of-transit spectra (i.e., the division approach). We then compare our measurements with a radiative transfer atmospheric model. Our analysis results in a tentative detection of exoplanetary sodium: we measure the width and amplitude of the exoplanetary sodium feature to be $\sigma_{\mathrm{Na}}$ = (0.128 $\pm$ 0.078) $\AA$ and A$_{\mathrm{Na}}$ = (1.7 $\pm$ 0.9)% in the excess light curve approach and $\sigma_{\mathrm{Na}}$ = (0.850 $\pm$ 0.034) $\AA$ and A$_{\mathrm{Na}}$ = (1.3 $\pm$ 0.6)% in the division approach. By comparing our measurements with a simple atmospheric model, we retrieve an atmospheric temperature of 1550 $^{+170} _{-200}$ K and radius (at 0.1 bar) of 1.81 $\pm$ 0.02 R$_{\rm Jup}$ for WASP-17b.

The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of these systems. To this end, we have developed a photodynamical model, and we demonstrate its performance using the Kepler-9 system as an example. Our comprehensive analysis combines the numerical integration of the system's dynamics over the time span of the observations along with the transit light curve model. This model is coupled with a Markov chain Monte Carlo algorithm, allowing the exploration of the model parameter space. Applied to the Kepler-9 long cadence data, short cadence data and 13 new transit observations collected by KOINet between the years 2014 to 2017, our modelling provides well constrained predictions for the next transits and the system's parameters. We have determined the densities of the planets Kepler-9b and 9c to the very precise values of rho_b = 0.439 +/-0.023 g/cm3 and rho_c = 0.322 +/- 0.017 g/cm3. Our analysis reveals that Kepler-9c will stop transiting in about 30 years. This results from strong dynamical interactions between Kepler-9b and 9c, near 2:1 resonance, that leads to a periodic change in inclination. Over the next 30 years the inclination of Kepler-9c (-9b) will decrease (increase) slowly. This should be measurable by a substantial decrease (increase) in the transit duration, in as soon as a few years' time. Observations that contradict this prediction might indicate the presence of additional objects. If this prediction proves true, this behaviour opens up a unique chance to scan the different latitudes of a star.

During its four years of photometric observations, the Kepler space telescope detected thousands of exoplanets and exoplanet candidates. One of Kepler's greatest heritages has been the confirmation and characterization of hundreds of multi-planet systems via Transit Timing Variations (TTVs). However, there are many interesting candidate systems displaying TTVs on such long time scales that the existing Kepler observations are of insufficient length to confirm and characterize them by means of this technique. To continue with Kepler's unique work we have organized the "Kepler Object of Interest Network" (KOINet). The goals of KOINet are, among others, to complete the TTV curves of systems where Kepler did not cover the interaction timescales well. KOINet has been operational since March, 2014. Here we show some promising first results obtained from analyzing seven primary transits of KOI-0410.01, KOI-0525.01, KOI-0760.01, and KOI-0902.01 in addition to Kepler data, acquired during the first and second observing seasons of KOINet. While carefully choosing the targets we set demanding constraints about timing precision (at least 1 minute) and photometric precision (as good as 1 part per thousand) that were achieved by means of our observing strategies and data analysis techniques. For KOI-0410.01, new transit data revealed a turn-over of its TTVs. We carried out an in-depth study of the system, that is identified in the NASA's Data Validation Report as false positive. Among others, we investigated a gravitationally-bound hierarchical triple star system, and a planet-star system. While the simultaneous transit fitting of ground and space-based data allowed for a planet solution, we could not fully reject the three-star scenario. New data, already scheduled in the upcoming 2018 observing season, will set tighter constraints on the nature of the system.

HD 189733 is an exoplanetary system consisting of a transiting hot Jupiter and an active K2V-type main sequence star. We aim to use VLT/UVES high resolution echelle spectra to study a stellar flare. We have performed simultaneous analyses of the temporal evolution in several chromospheric stellar lines, namely, the Ca II H and K lines, Halpha, Hbeta, Hgamma, Hdelta, Hepsilon, the Ca II infrared triplet line, and He I D3. Observations were carried out with a time resolution of approximately 1 min for a duration of four hours, including a complete planetary transit. We determine the energy released during the flare in all studied chromospheric lines combined to be about 8.7e31 erg, which puts this event at the upper end of flare energies observed on the Sun. Our analysis does not reveal any significant delay of the flare peak observed in the Balmer and Ca II H and K lines, although we find a clear difference in the temporal evolution of these lines. The He I D3 shows additional absorption possibly related to the flare event. Based on the flux released in Ca II H and K lines during the flare, we estimate the soft X-ray flux emission to be 7e30 erg. The observed flare can be ranked as a moderate flare on a K-type star and confirms a rather high activity level of HD 189733 host star. The cores of the studied chromospheric lines demonstrate the same behavior and let us study the flare evolution. We demonstrate that the activity of an exoplanet host star can play an important role in the detection of exoplanet atmospheres, since these are frequently discovered as an additional absorption in the line cores. A possible star-planet interaction responsible for a flare occurrence during a transit can neither be confirmed nor ruled out.

During primary transits, the spectral signatures of an exoplanet atmosphere can be measured using transmission spectroscopy. The goal of this work is to accurately measure the atomspheric sodium absorption light curve in HD189733b, correcting for the effects of stellar differential limb-darkening, stellar activity and a "bump" caused by the changing radial velocity of the exoplanet. In fact, due to the high cadence and quality of our data, it is the first time that the last feature can be detected even by visual inspection. We use 244 high-resolution optical spectra taken by the UVES instrument mounted at the VLT. Our observations cover a full transit of HD 189733b, with a cadence of 45 seconds. To probe the transmission spectrum of sodium we produce excess light- curves integrating the stellar flux in passbands of 1 \AA, 1.5 \AA, and 3 \AA inside the core of each sodium D-line. We model the effects of external sources on the excess light-curves, which correspond to an observed stellar flare beginning close to mid-transit time and the wavelength dependent limb-darkening effects. In addition, by characterizing the effect of the changing radial velocity and Doppler shifts of the planetary sodium lines inside the stellar sodium lines, we estimate the depth and width of the exoplanetary sodium feature. We estimate the shape of the planetary sodium line by a Gaussian profile, with an equivalent width of 0.0023 \pm 0.0010 \AA, thereby confirming the presence of sodium in the atmosphere of HD189733b with excess absorption levels of 0.72 \pm 0.25%, 0.34 \pm 0.11%, and 0.20 \pm 0.06% for the integration bands of 1 \AA, 1.5 \AA, and 3 \AA, respectively. From these, we produce a first order estimate of the number density of sodium in the exoplanet atmosphere.

The formation of dust gaps in protoplanetary disks is one of the most important signposts of disk evolution and possibly the formation of planets. We aim to characterize the 'flaring' disk structure around the Herbig Ae/Be stars HD 100453 and HD 34282. Their spectral energy distributions (SEDs) show an emission excess between 15-40{\mu}m, but very weak (HD 100453) and no (HD 34282) signs of the 10 and 20 {\mu}m amorphous silicate features. We investigate whether this implies the presence of large dust gaps. In this work, spatially resolved mid-infrared Q-band images taken with Gemini North/MICHELLE are investigated. We perform radiative transfer modeling and examine the radial distribution of dust. We simultaneously fit the Q-band images and SEDs of HD 100453 and HD 34282. Our solutions require that the inner-halos and outer-disks are likely separated by large dust gaps that are depleted wih respect to the outer disk by a factor of 1000 or more. The inner edges of the outer disks of HD 100453 and HD 34282 have temperatures of about $160 \pm 10$ K and $60 \pm 5$ K respectively. Because of the high surface brightnesses of these walls, they dominate the emission in the Q-band. Their radii are constrained at 20+2 AU and 92+31 AU, respectively. We conclude that, HD 100453 and HD 34282 likely have disk dust gaps and the upper limit on the dust mass in each gap is estimated to be about $10^{-7}$M$_{\odot}$. We find that the locations and sizes of disk dust gaps are connected to the SED, as traced by the mid-infrared flux ratio F30/F13.5. We propose a new classification scheme for the Meeus groups (Meeus et al. 2001) based on the F30/F13.5 ratio. The absence of amorphous silicate features in the observed SEDs is caused by the depletion of small (smaller than 1 {\mu}m) silicate dust at temperatures above 160 K, which could be related to the presence of a dust gap in that region of the disk.

The center-to-limb variation (CLV) describes the brightness of the stellar disk as a function of the limb angle. Across strong absorption lines, the CLV can vary quite significantly. We obtained a densely sampled time series of high-resolution transit spectra of the active planet host star HD 189733 with UVES. Using the passing planetary disk of the hot Jupiter HD 189733 b as a probe, we study the CLV in the wings of the Ca II H and K and Na I D1 and D2 Fraunhofer lines, which are not strongly affected by activity-induced variability. In agreement with model predictions, our analysis shows that the wings of the studied Fraunhofer lines are limb brightened with respect to the (quasi-)continuum. The strength of the CLV-induced effect can be on the same order as signals found for hot Jupiter atmospheres. Therefore, a careful treatment of the wavelength dependence of the stellar CLV in strong absorption lines is highly relevant in the interpretation of planetary transit spectroscopy.