We present a Python package LDTK that automates the calculation of custom stellar limb darkening (LD) profiles and model-specific limb darkening coefficients (LDC) using the library of PHOENIX-generated specific intensity spectra by Husser et al. (2013). The aim of the package is to facilitate analyses requiring custom generated limb darkening profiles, such as the studies of exoplanet transits-especially transmission spectroscopy, where the transit modelling is carried out for custom narrow passbands-eclipsing binaries (EBs), interferometry, and microlensing events. First, LDTK can be used to compute custom limb darkening profiles with uncertainties propagated from the uncertainties in the stellar parameter estimates. Second, LDTK can be used to estimate the limb-darkening-model specific coefficients with uncertainties for the most common limb-darkening models. Third, LDTK can be directly integrated into the log posterior computation of any pre-existing modelling code with minimal modifications. The last approach can be used to constrain the LD model parameter space directly by the LD profile, allowing for the marginalization over the LD parameter space without the need to approximate the constraint from the LD profile using a prior. Where µ = √ 1 − z 2 = cos γ, z is the normalized distance from the centre of the stellar disk, and γ is the foreshortening angle.

We present k2sc (K2 Systematics Correction), a Python pipeline to model instrumental systematics and astrophysical variability in light curves from the K2 mission. k2sc uses Gaussian process regression to model position-dependent systematics and time-dependent variability simultaneously, enabling the user to remove both (e.g., for transit searches) or to remove systematics while preserving variability (for variability studies). For periodic variables, k2sc automatically computes estimates of the period, amplitude and evolution timescale of the variability. We apply k2sc to publicly available K2 data from campaigns 3-5 showing that we obtain photometric precision approaching that of the original Kepler mission. We compare our results to other publicly available K2 pipelines, showing that we obtain similar or better results, on average. We use transit injection and recovery tests to evaluate the impact of k2sc on planetary transit searches in K2 pdc (Pre-search Data Conditioning) data, for planet-to-star radius ratio down R p /R = 0.01 and periods up to P = 40 d, and show that k2sc significantly improves the ability to distinguish between correct and false detections, particularly for small planets. k2sc can be run automatically on many light curves, or manually tailored for specific objects such as pulsating stars or large amplitude eclipsing binaries. It can be run on ASCII and FITS light curve files, regardless of their origin. Both the code and the processed light curves are publicly available, and we provide instructions for downloading and using them. The methodology used by k2sc will be applicable to future transit search missions such as TESS and PLATO.

We present a fast and user friendly exoplanet transit light curve modelling package PyTransit, implementing optimised versions of the Giménez (2006) and Mandel & Agol (2002) transit models. The package offers an object-oriented Python interface to access the two models implemented natively in Fortran with OpenMP parallelisation. A partial OpenCL version of the quadratic Mandel-Agol model is also included for GPU-accelerated computations. The aim of PyTransit is to facilitate the analysis of photometric time series of exoplanet transits consisting of hundreds of thousands of datapoints, and of multi-passband transit light curves from spectrophotometric observations, as a part of a researchers programming toolkit for building complex, problem-specific, analyses.

Astronomers have discovered thousands of planets outside the solar system 1 , most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star 2 , but more distant planets can survive this phase and remain in orbit around the white dwarf 3,4 . Some white dwarfs show evidence for rocky material floating in their atmospheres 5 , in warm debris disks [6][7][8][9] , or orbiting very closely [10][11][12] , which has been interpreted as the debris of rocky planets that were scattered inward and tidally disrupted 13 . Recently, the discovery of a gaseous debris disk with a composition similar to ice giant planets 14 demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether the planets can survive the journey. So far, the detection of intact planets in close orbits around white dwarfs has remained elusive. Here, we report the discovery of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95% confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red-giant phase and shrinks due to friction. In this case, though, the low mass and relatively long orbital period of the planet candidate make common-envelope evolution less likely. Instead, the WD 1856+534 system seems to demonstrate that giant planets can be scattered into tight orbits without being tidally disrupted, and motivates searches for smaller transiting planets around white dwarfs. WD 1856+534 (hereafter, WD 1856 for brevity) is located 25 parsecs away in a visual triple star system. It has an effective temperature of 4710 ± 60 Kelvin and became a white dwarf 5.9 ± 0.5 billion years ago, based on theoretical models for how white dwarfs cool over time. The total system age, including the star's main sequence lifetime, must be older. Table 1 gives the other key parameters of the star. WD 1856 is one of thousands of white dwarfs that was targeted for observations with NASA's Transiting Exoplanet Survey Satellite (TESS ), in order to search for any periodic dimming events caused by planetary transits. A statistically significant transit-like event was detected by the TESS Science Processing Operations Center (SPOC) pipeline based

We present TRICERATOPS, a new Bayesian tool that can be used to vet and validate TESS Objects of Interest (TOIs). We test the tool on 68 TOIs that have been previously confirmed as planets or rejected as astrophysical false positives. By looking in the false-positive probability (FPP)−nearby false-positive probability (NFPP) plane, we define criteria that TOIs must meet to be classified as validated planets (FPP < 0.015 and NFPP < 10−3), likely planets (FPP < 0.5 and NFPP < 10−3), and likely nearby false positives (NFPP > 10−1). We apply this procedure on 384 unclassified TOIs and statistically validate 12, classify 125 as likely planets, and classify 52 as likely nearby false positives. Of the 12 statistically validated planets, 9 are newly validated. TRICERATOPS is currently the only TESS vetting and validation tool that models transits from nearby contaminant stars in addition to the target star. We therefore encourage use of this tool to prioritize follow-up observations that confirm bona fide planets and identify false positives originating from nearby stars.

We introduce a new transit search and vetting pipeline for observations from the K2 mission, and present the candidate transiting planets identified by this pipeline out of the targets in Campaigns 5 and 6. Our pipeline uses the Gaussian Process-based k2sc code to correct for the K2 pointing systematics and simultaneously model stellar variability. The systematicscorrected, variability-detrended light curves are searched for transits with the Box Least Squares method, and a period-dependent detection threshold is used to generate a preliminary candidate list. Two or three individuals vet each candidate manually to produce the final candidate list, using a set of automatically-generated transit fits and assorted diagnostic tests to inform the vetting. We detect 147 single-planet system candidates and 5 multi-planet systems, independently recovering the previously-published hot Jupiters EPIC 212110888b, WASP55b (EPIC 212300977b) and Qatar-2b (EPIC 212756297b). We also report the outcome of reconnaissance spectroscopy carried out for all candidates with Kepler magnitude K p ≤ 13, identifying 12 targets as likely false positives. We compare our results to those of other K2 transit search pipelines, noting that ours performs particularly well for variable and/or active stars, but that the results are very similar overall. All the light curves and code used in the transit search and vetting process are publicly available, as are the follow-up spectra.

We present a detailed analysis of HARPS-N radial velocity observations of K2-100, a young and active star in the Praesepe cluster, which hosts a transiting planet with a period of 1.7 days. We model the activity-induced radial velocity variations of the host star with a multi-dimensional Gaussian Process framework and detect a planetary signal of 10.6 ± 3.0 m s −1 which matches the transit ephemeris, and translates to a planet mass of 21.8 ± 6.2 M ⊕ . We perform a suite of validation tests to confirm that our detected signal is genuine. This is the first mass measurement for a transiting planet in a young open cluster. The relatively low density of the planet, 2.04 +0.66 −0.61 g cm −3 , implies that K2-100b retains a significant volatile envelope. We estimate that the planet is losing its atmosphere at a rate of 10 11 − 10 12 g s −1 due to the high level of radiation it receives from its host star.

Ultra-hot Jupiters orbit very close to their host star and consequently receive strong irradiation, causing their atmospheric chemistry to be different from the common gas giants. Here, we have studied the atmosphere of one of these particular hot planets, MASCARA-2b/KELT-20b, using four transit observations with high resolution spectroscopy facilities. Three of these observations were performed with HARPS-N and one with CARMENES. Additionally, we simultaneously observed one of the transits with MuSCAT2 to monitor possible spots in the stellar surface. At high resolution, the transmission residuals show the effects of Rossiter-McLaughlin and centre-to-limb variations from the stellar lines profiles, which we have corrected to finally extract the transmission spectra of the planet. We clearly observe the absorption features of CaII, FeII, NaI, Hα, and Hβ in the atmosphere of MASCARA-2b, and indications of Hγ and MgI at low signal-to-noise ratio. In the case of NaI, the true absorption is difficult to disentangle from the strong telluric and interstellar contamination. The results obtained with CARMENES and HARPS-N are consistent, measuring an Hα absorption depth of 0.68 ± 0.05 and 0.59 ± 0.07%, and NaI absorption of 0.11 ± 0.04 and 0.09 ± 0.05% for a 0.75 Å passband, in the two instruments respectively. The Hα absorption corresponds to ~1.2 Rp, which implies an expanded atmosphere, as a result of the gas heating caused by the irradiation received from the host star. For Hβ and Hγ only HARPS-N covers this wavelength range, measuring an absorption depth of 0.28 ± 0.06 and 0.21 ± 0.07%, respectively. For CaII, only CARMENES covers this wavelength range measuring an absorption depth of 0.28 ± 0.05, 0.41 ± 0.05 and 0.27 ± 0.06% for CaII λ8498Å, λ8542Å and λ8662Å lines, respectively. Three additional absorption lines of FeII are observed in the transmission spectrum by HARPS-N (partially covered by CARMENES), measuring an average absorption depth of 0.08 ± 0.04% (0.75 Å passband). The results presented here are consistent with theoretical models of ultra-hot Jupiters atmospheres, suggesting the emergence of an ionised gas on the day-side of such planets. Calcium and iron, together with other elements, are expected to be singly ionised at these temperatures and be more numerous than its neutral state. The Calcium triplet lines are detected here for the first time in transmission in an exoplanet atmosphere.