Search

Stellar active regions like spots and faculae can distort the shapes of spectral lines, inducing variations in the radial velocities that are often orders of magnitude larger than the signals from Earth-like planets. Efforts to mitigate these activity signals have hitherto focused on either the time or the velocity (wavelength) domains. We present a physics-driven Gaussian process (GP) framework to model activity signals directly in time series of line profiles or Cross-Correlation Functions (CCFs). Unlike existing methods which correct activity signals in line profile time series, our approach exploits the time correlation between velocity (wavelength) bins in the line profile variations, and is based on a simplified but physically motivated model for the origin of these variations. When tested on both synthetic and real data sets with signal-to-noise ratios down to $\sim$ 100, our method was able to separate the planetary signal from the activity signal, even when their periods were identical. We also conducted injection/recovery tests using two years of realistically sampled HARPS-N solar data, demonstrating the ability of the method to accurately recover a signal induced by a 1.5-Earth mass planet with a semi-amplitude of 0.3 m/s and a period of 33 days during high solar activity., Comment: Accepted for publication in MNRAS

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

Stellar magnetic activity induces both distortions and Doppler-shifts in the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with the high-precision spectrograph HARPS-N. We show that the RV signal arising from line-shape variations on time-scales associated with the solar rotation and activity cycle can be robustly extracted from the data, reducing the RV dispersion by half. Once these have been corrected, activity-induced Doppler-shifts remain, that are modulated at the solar rotation period, and that are most effectively modelled in the time domain, using Gaussian Processes (GPs). Planet signatures are still best retrieved with multi-dimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca II H and K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6-0.8 m/s, likely to be dominated by signals induced by super-granulation. Finally, we find that the statistical properties of the RVs evolve significantly over time, and that this evolution is primarily driven by sunspots, which control the smoothness of the signal. Such evolution, which reduces the sensitivity to long-period planet signatures, is no longer seen in the activity-induced Doppler-shifts, which is promising for long term RV monitoring surveys such as the Terra Hunting Experiment or the PLATO follow-up campaign., Comment: Accepted by MNRAS on the 2024 May 20. 17 pages, 14 figures (plus 8 pages of Appendix, 9 Appendix figures)

We report on the discovery and validation of a transiting long-period mini-Neptune orbiting a bright (V = 9.0 mag) G dwarf (TOI 4633; R = 1.05 RSun, M = 1.10 MSun). The planet was identified in data from the Transiting Exoplanet Survey Satellite by citizen scientists taking part in the Planet Hunters TESS project. Modeling of the transit events yields an orbital period of 271.9445 +/- 0.0040 days and radius of 3.2 +/- 0.20 REarth. The Earth-like orbital period and an incident flux of 1.56 +/- 0.2 places it in the optimistic habitable zone around the star. Doppler spectroscopy of the system allowed us to place an upper mass limit on the transiting planet and revealed a non-transiting planet candidate in the system with a period of 34.15 +/- 0.15 days. Furthermore, the combination of archival data dating back to 1905 with new high angular resolution imaging revealed a stellar companion orbiting the primary star with an orbital period of around 230 years and an eccentricity of about 0.9. The long period of the transiting planet, combined with the high eccentricity and close approach of the companion star makes this a valuable system for testing the formation and stability of planets in binary systems., Comment: 24 pages, 16 figures, 4 tables

We present an exhaustive photometric and spectroscopic analysis of TOI-837, a F9/G0 35 Myr young star, hosting a transiting exoplanet, TOI-837b, with an orbital period of 8.32d. Utilising data from TESS and ground-based observations, we determine a planetary radius of 0.82 R_J for TOI-837b. Through detailed HARPS spectroscopic time series analysis, we derive a Doppler semi-amplitude of 35 m/s, corresponding to a planetary mass of 0.39 M_J. The derived planetary properties suggest a substantial core of approximately 70 M_E, constituting about 60% of the planet's total mass. This finding poses a significant challenge to existing theoretical models of core formation. We propose that future atmospheric observations with JWST could provide insights into resolving ambiguities of TOI-837b, offering new perspectives on its composition, formation, and evolution., Comment: MNRAS published version, MNRAS, 531, 4275

As the hunt for an Earth-like exoplanets has intensified in recent years, so has the effort to characterise and model the stellar signals that can hide or mimic small planetary signals. Stellar variability arises from a number of sources, including granulation, supergranulation, oscillations and activity, all of which result in quasi-periodic or stochastic behaviour in photometric and/or radial velocity observations. Traditionally, the characterisation of these signals has mostly been done in the frequency domain. However, the recent development of scalable Gaussian process regression methods makes direct time-domain modelling of stochastic processes a feasible and arguably preferable alternative, obviating the need to estimate the power spectral density of the data before modelling it. In this paper, we compare the two approaches using a series of experiments on simulated data. We show that frequency domain modelling can lead to inaccurate results, especially when the time sampling is irregular. By contrast, Gaussian process regression results are often more precise, and systematically more accurate, in both the regular and irregular time sampling regimes. While this work was motivated by the analysis of radial velocity and photometry observations of main sequence stars in the context of planet searches, we note that our results may also have applications for the study of other types of astrophysical variability such as quasi-periodic oscillations in X-ray binaries and active galactic nuclei variability., Comment: Accepted for publication in MNRAS

Space-based photometry missions produce exquisite light curves that contain a wealth of stellar variability on a wide range of timescales. Light curves also typically contain significant instrumental systematics -- spurious, non-astrophysical trends that are common, in varying degrees, to many light curves. Empirical systematics-correction approaches using the information in the light curves themselves have been very successful, but tend to suppress astrophysical signals, particularly on longer timescales. Unlike its predecessors, the PLATO mission will use multiple cameras to monitor the same stars. We present REPUBLIC, a novel systematics-correction algorithm which exploits this multi-camera configuration to correct systematics that differ between cameras, while preserving the component of each star's signal that is common to all cameras, regardless of timescale. Through simulations with astrophysical signals (star spots and planetary transits), Kepler-like errors, and white noise, we demonstrate REPUBLIC's ability to preserve long-term astrophysical signals usually lost in standard correction techniques. We also explore REPUBLIC's performance with different number of cameras and systematic properties. We conclude that REPUBLIC should be considered a potential complement to existing strategies for systematic correction in multi-camera surveys, with its utility contingent upon further validation and adaptation to the specific characteristics of the PLATO mission data, Comment: 12 pages, 6 figures, accepted for publication in RAS Techniques and Instruments

The detection of planetary transits in the light curves of active stars, featuring correlated noise in the form of stellar variability, remains a challenge. Depending on the noise characteristics, we show that the traditional technique that consists of detrending a light curve before searching for transits alters their signal-to-noise ratio, and hinders our capability to discover exoplanets transiting rapidly-rotating active stars. We present nuance, an algorithm to search for transits in light curves while simultaneously accounting for the presence of correlated noise, such as stellar variability and instrumental signals. We assess the performance of nuance on simulated light curves as well as on the TESS light curves of 438 rapidly-rotating M dwarfs. For each dataset, we compare our method to 5 commonly-used detrending techniques followed by a search with the Box-Least-Square algorithm. Overall, we demonstrate that nuance is the most performant method in 93% of cases, leading to both the highest number of true positives and the lowest number of false positive detections. Although simultaneously searching for transits while modeling correlated noise is expected to be computationally expensive, we make our algorithm tractable and available as the JAX-powered Python package nuance, allowing its use on distributed environments and GPU devices. Finally, we explore the prospects offered by the nuance formalism, and its use to advance our knowledge of planetary systems around active stars, both using space-based surveys and sparse ground-based observations.

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S+LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to a) measure rotation periods for 49 slow rotators in our sample, b) explore the impact of these results on the spin-down of middle-aged late F, G & K stars, and c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations program, and help fine-tune current stellar structure and evolution models., Comment: Accepted for publication in MNRAS

We present ground and space-based photometric observations of TOI-270 (L231-32), a system of three transiting planets consisting of one super-Earth and two sub-Neptunes discovered by TESS around a bright (K-mag=8.25) M3V dwarf. The planets orbit near low-order mean-motion resonances (5:3 and 2:1), and are thus expected to exhibit large transit timing variations (TTVs). Following an extensive observing campaign using 8 different observatories between 2018 and 2020, we now report a clear detection of TTVs for planets c and d, with amplitudes of $\sim$10 minutes and a super-period of $\sim$3 years, as well as significantly refined estimates of the radii and mean orbital periods of all three planets. Dynamical modeling of the TTVs alone puts strong constraints on the mass ratio of planets c and d and on their eccentricities. When incorporating recently published constraints from radial velocity observations, we obtain masses of $M_{\mathrm{b}}=1.48\pm0.18\,M_\oplus$, $M_{c}=6.20\pm0.31\,M_\oplus$ and $M_{\mathrm{d}}=4.20\pm0.16\,M_\oplus$ for planets b, c and d, respectively. We also detect small, but significant eccentricities for all three planets : $e_\mathrm{b} =0.0167\pm0.0084$, $e_{c} =0.0044\pm0.0006$ and $e_{d} = 0.0066\pm0.0020$. Our findings imply an Earth-like rocky composition for the inner planet, and Earth-like cores with an additional He/H$_2$O atmosphere for the outer two. TOI-270 is now one of the best-constrained systems of small transiting planets, and it remains an excellent target for atmospheric characterization., Comment: 22 pages, 7 figures

Detecting planetary signatures in radial velocity time-series of young stars is challenging due to their inherently strong stellar activity. However, it is possible to learn information about the properties of the stellar signal by using activity indicators measured from the same stellar spectra used to extract radial velocities. In this manuscript, we present a reanalysis of spectroscopic HARPS data of the young star K2-233, which hosts three transiting planets. We perform a multidimensional Gaussian Process regression on the radial velocity and the activity indicators to characterise the planetary Doppler signals. We demonstrate, for the first time on a real dataset, that the use of a multidimensional Gaussian Process can boost the precision with which we measure the planetary signals compared to a one-dimensional Gaussian Process applied to the radial velocities alone. We measure the semi-amplitudes of K2-233 b, c, and d as 1.31(-0.74)(+0.81), 1.81(-0.67)(+0.71), and 2.72(-0.70)(+0.66) m/s, which translates into planetary masses of 2.4(-1.3)(+1.5), 4.6(-1.7)(+1.8), and 10.3(-2.6)(+2.4), respectively. These new mass measurements make K2-233 d a valuable target for transmission spectroscopy observations with JWST. K2-233 is the only young system with two detected inner planets below the radius valley and a third outer planet above it. This makes it an excellent target to perform comparative studies, to inform our theories of planet evolution, formation, migration, and atmospheric evolution., Comment: Accepted for publication in MNRAS

The last two decades have seen a major expansion in the availability, size, and precision of time-domain datasets in astronomy. Owing to their unique combination of flexibility, mathematical simplicity and comparative robustness, Gaussian Processes (GPs) have emerged recently as the solution of choice to model stochastic signals in such datasets. In this review we provide a brief introduction to the emergence of GPs in astronomy, present the underlying mathematical theory, and give practical advice considering the key modelling choices involved in GP regression. We then review applications of GPs to time-domain datasets in the astrophysical literature so far, from exoplanets to active galactic nuclei, showcasing the power and flexibility of the method. We provide worked examples using simulated data, with links to the source code, discuss the problem of computational cost and scalability, and give a snapshot of the current ecosystem of open source GP software packages. Driven by further algorithmic and conceptual advances, we expect that GPs will continue to be an important tool for robust and interpretable time domain astronomy for many years to come., Comment: Updated in response to community feedback; source code: https://github.com/dfm/araa-gps

In recent years, Gaussian Process (GP) regression has become widely used to analyse stellar and exoplanet time-series data sets. For spotted stars, the most popular GP covariance function is the quasi-periodic (QP) kernel, whose the hyperparameters of the GP have a plausible interpretation in terms of physical properties of the star and spots. In this paper, we test the reliability of this interpretation by modelling data simulated using a spot model using a QP GP, and the recently proposed quasi-periodic plus cosine (QPC) GP, comparing the posterior distributions of the GP hyperparameters to the input parameters of the spot model. We find excellent agreement between the input stellar rotation period and the QP and QPC GP period, and very good agreement between the spot decay timescale and the length scale of the squared exponential term. We also compare the hyperparameters derived from light and radial velocity (RV) curves for a given star, finding that the period and evolution timescales are in good agreement. However, the harmonic complexity of the GP, while displaying no clear correlation with the spot properties in our simulations, is systematically higher for the RV than for the light curve data. Finally, for the QP kernel, we investigate the impact of noise and time-sampling on the hyperparameters in the case of RVs. Our results indicate that good coverage of rotation period and spot evolution time-scales is more important than the total number of points, and noise characteristics govern the harmonic complexity., Comment: 18 pages, 24 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Observations of exoplanet atmospheres in high resolution have the potential to resolve individual planetary absorption lines, despite the issues associated with ground-based observations. The removal of contaminating stellar and telluric absorption features is one of the most sensitive steps required to reveal the planetary spectrum and, while many different detrending methods exist, it remains difficult to directly compare the performance and efficacy of these methods. Additionally, though the standard cross-correlation method enables robust detection of specific atmospheric species, it only probes for features that are expected a priori. Here we present a novel methodology using Gaussian process (GP) regression to directly model the components of high-resolution spectra, which partially addresses these issues. We use two archival CRIRES/VLT data sets as test cases, observations of the hot Jupiters HD 189733 b and 51 Pegasi b, recovering injected signals with average line contrast ratios of $\sim 4.37 \times 10^{-3}$ and $\sim 1.39 \times 10^{-3}$, and planet radial velocities $\Delta K_\mathrm{p} =1.45 \pm 1.53\,\mathrm{km\,s^{-1}}$ and $\Delta K_\mathrm{p}=0.12\pm0.12\,\mathrm{km\,s^{-1}}$ from the injection velocities respectively. In addition, we demonstrate an application of the GP method to assess the impact of the detrending process on the planetary spectrum, by implementing injection-recovery tests. We show that standard detrending methods used in the literature negatively affect the amplitudes of absorption features in particular, which has the potential to render retrieval analyses inaccurate. Finally, we discuss possible limiting factors for the non-detections using this method, likely to be remedied by higher signal-to-noise data., Comment: 14 pages, 7 figures, 3 tables. Accepted for publication in MNRAS

We present a spectroscopic analysis of a 1-year intensive monitoring campaign of the 22-Myr old planet-hosting M dwarf AU Mic using the HARPS spectrograph. In a companion paper, we reported detections of the planet radial velocity (RV) signatures of the two close-in transiting planets of the system, with respective semi-amplitudes of 5.8 $\pm$ 2.5 m/s and 8.5 $\pm$ 2.5 m/s for AU Mic b and AU Mic c. Here, we perform an independent measurement of the RV semi-amplitude of AU Mic c using Doppler imaging to simultaneously model the activity-induced distortions and the planet-induced shifts in the line profiles. The resulting semi-amplitude of 13.3 $\pm$ 4.1 m/s for AU Mic c reinforces the idea that the planet features a surprisingly large inner density, in tension with current standard models of core accretion. Our brightness maps feature significantly higher spot coverage and lower level of differential rotation than the brightness maps obtained in late 2019 with the SPIRou spectropolarimeter, suggesting that the stellar magnetic activity has evolved dramatically over a $\sim$1-yr time span. Additionally, we report a 3-$\sigma$ detection of a modulation at 8.33 $\pm$ 0.04 d of the He I D3 (587.562 nm) emission flux, close to the 8.46-d orbital period of AU Mic b. The power of this emission (a few 10$^{17}$ W) is consistent with 3D magnetohydrodynamical simulations of the interaction between stellar wind and the close-in planet if the latter hosts a magnetic field of $\sim$10 G. Spectropolarimetric observations of the star are needed to firmly elucidate the origin of the observed chromospheric variability., Comment: Accepted for publication in MNRAS (15/03/2022); 19 pages including 6 pages of appendices; 14 figures including 5 in appendices

The system of two transiting Neptune-sized planets around the bright, young M-dwarf AU Mic provides a unique opportunity to test models of planet formation, early evolution, and star-planet interaction. However, the intense magnetic activity of the host star makes measuring the masses of the planets via the radial velocity (RV) method very challenging. We report on a 1-year, intensive monitoring campaign of the system using 91 observations with the HARPS spectrograph, allowing for detailed modelling of the $\sim 600\, \rm{m\,s^{-1}}$ peak-to-peak activity-induced RV variations. We used a multidimensional Gaussian Process framework to model these and the planetary signals simultaneously. We detect the latter with semi-amplitudes of $\rm{K_b} = 5.8 \pm 2.5\, \rm{m\,s^{-1}}$ and $\rm{K_c} = 8.5 \pm 2.5\, \rm{m\, s^{-1}}$, respectively. The resulting mass estimates, $\rm{M_b} = 11.7 \pm 5.0\, \rm{M_{\rm \oplus}}$ and $\rm{M_c} = 22.2 \pm 6.7\, \rm{M_{\rm \oplus}}$, suggest that planet b might be less dense, and planet c considerably denser than previously thought. These results are in tension with the current standard models of core-accretion. They suggest that both planets accreted a H/He envelope that is smaller than expected, and the trend between the two planets' envelope fractions is the opposite of what is predicted by theory.

We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14-17 M$_{\odot}$) is more massive than the combined mass of the inner binary (10.9-13.2 M$_{\odot}$). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne-Zytkow object. Further observational characterisation of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products., Comment: Accepted for publication in MNRAS (11 pages, 8 figures)

Measured spectral shifts due to intrinsic stellar variability (e.g., pulsations, granulation) and activity (e.g., spots, plages) are the largest source of error for extreme precision radial velocity (EPRV) exoplanet detection. Several methods are designed to disentangle stellar signals from true center-of-mass shifts due to planets. The EXPRES Stellar Signals Project (ESSP) presents a self-consistent comparison of 22 different methods tested on the same extreme-precision spectroscopic data from EXPRES. Methods derived new activity indicators, constructed models for mapping an indicator to the needed RV correction, or separated out shape- and shift-driven RV components. Since no ground truth is known when using real data, relative method performance is assessed using the total and nightly scatter of returned RVs and agreement between the results of different methods. Nearly all submitted methods return a lower RV RMS than classic linear decorrelation, but no method is yet consistently reducing the RV RMS to sub-meter-per-second levels. There is a concerning lack of agreement between the RVs returned by different methods. These results suggest that continued progress in this field necessitates increased interpretability of methods, high-cadence data to capture stellar signals at all timescales, and continued tests like the ESSP using consistent data sets with more advanced metrics for method performance. Future comparisons should make use of various well-characterized data sets -- such as solar data or data with known injected planetary and/or stellar signals -- to better understand method performance and whether planetary signals are preserved., Comment: 33 pages (+12 pages of Appendix), 10 figures, 8 tables, accepted for publication in AJ

The two most successful methods for exoplanet detection rely on the detection of planetary signals in photometric and radial velocity time-series. This depends on numerical techniques that exploit the synergy between data and theory to estimate planetary, orbital, and/or stellar parameters. In this work, we present a new version of the exoplanet modelling code pyaneti. This new release has a special emphasis on the modelling of stellar signals in radial velocity time-series. The code has a built-in multidimensional Gaussian process approach to modelling radial velocity and activity indicator time-series with different underlying covariance functions. This new version of the code also allows multi-band and single transit modelling; it runs on Python 3, and features overall improvements in performance. We describe the new implementation and provide tests to validate the new routines that have direct application to validate the new routines that have direct application to exoplanet detection and characterisation. We have made the code public and freely available at https://github.com/oscaribv/pyaneti. We also present the codes citlalicue and citlalatonac that allow one to create synthetic photometric and spectroscopic time-series, respectively, with planetary and stellar-like signals., Comment: This version looks closer to the published one: MNRAS, 509, 866

We report on the discovery and validation of a two-planet system around a bright (V = 8.85 mag) early G dwarf (1.43 $R_{\odot}$, 1.15 $M_{\odot}$, TOI 2319) using data from NASA's Transiting Exoplanet Survey Satellite (TESS). Three transit events from two planets were detected by citizen scientists in the month-long TESS light curve (sector 25), as part of the Planet Hunters TESS project. Modelling of the transits yields an orbital period of \Pb\ and radius of $3.41 _{ - 0.12 } ^ { + 0.14 }$ $R_{\oplus}$ for the inner planet, and a period in the range 19.26-35 days and a radius of $5.83 _{ - 0.14 } ^ { + 0.14 }$ $R_{\oplus}$ for the outer planet, which was only seen to transit once. Each signal was independently statistically validated, taking into consideration the TESS light curve as well as the ground-based spectroscopic follow-up observations. Radial velocities from HARPS-N and EXPRES yield a tentative detection of planet b, whose mass we estimate to be $11.56 _{ - 6.14 } ^ { + 6.58 }$ $M_{\oplus}$, and allow us to place an upper limit of $27.5$ $M_{\oplus}$ (99 per cent confidence) on the mass of planet c. Due to the brightness of the host star and the strong likelihood of an extended H/He atmosphere on both planets, this system offers excellent prospects for atmospheric characterisation and comparative planetology., Comment: Accepted for publication in MNRAS (15 pages, 10 figures, 3 tables)

We present a new analysis of the light curve of the young planet-hosting star TOI 451 in the light of new observations from TESS Cycle 3. Our joint analysis of the transits of all three planets, using all available TESS data, results in an improved ephemeris for TOI 451 b and TOI 451 c, which will help to plan follow-up observations. The updated mid-transit times are $\textrm{BJD}-2,457\,000=$ $1410.9896_{ - 0.0029 }^{ + 0.0032 }$, $1411.7982_{-0.0020}^{+0.0022}$, and $1416.63407_{-0.00100}^{+0.00096}$ for TOI 451 b, c, and d, respectively, and the periods are $1.8587028_{-10e-06}^{+08e-06}$, $9.192453_{-3.3e-05}^{+4.1e-05}$, and $16.364932_{-3.5e-05}^{+3.6e-05 }$ days. We also model the out-of-transit light curve using a Gaussian Process with a quasi-periodic kernel and infer a change in the properties of the active regions on the surface of TOI 451 between TESS Cycles 1 and 3., Comment: To appear as an RNAAS, 5 pages, 1 Figure

Binary supermassive black hole (SMBH) systems result from galaxy mergers, and will eventually coalesce due to gravitational wave (GW) emission if the binary separation can be reduced to $\lesssim 0.1$ pc by other mechanisms. Here, we explore a gravitational self-lensing binary SMBH model for the sharp (duration $\sim 1$ hr), quasi-regular X-ray flares -- dubbed quasi-periodic eruptions -- recently observed from two low mass active galactic nuclei: GSN 069 and RX J1301.9+2747. In our model, the binary is observed $\sim$edge-on, such that each SMBH gravitationally lenses light from the accretion disc surrounding the other SMBH twice per orbital period. The model can reproduce the flare spacings if the current eccentricity of RX J1301.9+2747 is $\epsilon_0 \gtrsim 0.16$, implying a merger within $\sim 1000$ yrs. However, we cannot reproduce the observed flare profiles with our current calculations. Model flares with the correct amplitude are $\sim 2/5$ the observed duration, and model flares with the correct duration are $\sim 2/5$ the observed amplitude. Our modelling yields three distinct behaviours of self-lensing binary systems that can be searched for in current and future X-ray and optical time-domain surveys: i) periodic lensing flares, ii) partial eclipses (caused by occultation of the background mini-disc by the foreground mini-disc), and iii) partial eclipses with a very sharp in-eclipse lensing flare. Discovery of such features would constitute very strong evidence for the presence of a supermassive binary, and monitoring of the flare spacings will provide a measurement of periastron precession., Comment: Accepted for publication in MNRAS

We present the results from the first two years of the Planet Hunters TESS citizen science project, which identifies planet candidates in the TESS data by engaging members of the general public. Over 22,000 citizen scientists from around the world visually inspected the first 26 Sectors of TESS data in order to help identify transit-like signals. We use a clustering algorithm to combine these classifications into a ranked list of events for each sector, the top 500 of which are then visually vetted by the science team. We assess the detection efficiency of this methodology by comparing our results to the list of TESS Objects of Interest (TOIs) and show that we recover 85 % of the TOIs with radii greater than 4 Earth radii and 51 % of those with radii between 3 and 4 Earth radii. Additionally, we present our 90 most promising planet candidates that had not previously been identified by other teams, 73 of which exhibit only a single transit event in the TESS light curve, and outline our efforts to follow these candidates up using ground-based observatories. Finally, we present noteworthy stellar systems that were identified through the Planet Hunters TESS project., Comment: Accepted for publication in MNRAS (22 pages, 12 figures, 3 tables)

Observational studies of exoplanets are suggestive of a ubiquitous presence of clouds. The current modelling techniques used in emission to account for the clouds tend to require prior knowledge of the cloud condensing species and often do not consider the scattering effects of the cloud. We explore the effects that thermal scattering has on the emission spectra by modelling a suite of hot Jupiter atmospheres with varying cloud single-scattering albedos (SSAs) and temperature profiles. We examine cases ranging from simple isothermal conditions to more complex structures and physically driven cloud modelling. We show that scattering from nightside clouds would lead to brightness temperatures that are cooler than the real atmospheric temperature if scattering is unaccounted for. We show that scattering can produce spectral signatures in the emission spectrum even for isothermal atmospheres. We identify the retrieval degeneracies and biases that arise in the context of simulated JWST spectra when the scattering from the clouds dominates the spectral shape. Finally, we propose a novel method of fitting the SSA spectrum of the cloud in emission retrievals, using a technique that does not require any prior knowledge of the cloud chemical or physical properties., Comment: Accepted to MNRAS

Oppenheimer et al. (1997) discovered two M5 dwarfs in the Pleiades with nearly primordial lithium. These stars are not low enough in mass to represent the leading edge of the lithium depletion boundary at Pleiades age (~125 Myr). A possible explanation for the enhanced lithium in these stars is that they are actually not members of the Pleiades but instead are members of a younger moving group seen in projection towards the Pleiades. We have used data from Gaia DR2 to confirm that these two stars, HHJ 339 and HHJ 430, are indeed not members of the Pleiades. Based on their space motions, parallaxes and positions in a Gaia-based CMD, it is probable that these two stars are about 40 parsecs foreground to the Pleiades and have ages of ~25 Myr. Kinematically they are best matched to the 32 Ori moving group., Comment: Accepted by AJ

We present Mon-735, a detached double-lined eclipsing binary (EB) member of the $\sim$3 Myr old NGC 2264 star forming region, detected by Spitzer. We simultaneously model the Spitzer light curves, follow-up Keck/HIRES radial velocities, and the system's spectral energy distribution to determine self-consistent masses, radii and effective temperatures for both stars. We find that Mon-735 comprises two pre-main sequence M dwarfs with component masses of $M = 0.2918 \pm 0.0099$ and $0.2661 \pm 0.0095$ $\rm{M}_{\odot}$, radii of $R = 0.762 \pm 0.022$ and $0.748 \pm 0.023$ $\rm{R}_{\odot}$, and effective temperatures of $T_{\rm eff} = 3260 \pm 73$ and $3213 \pm 73$ $\rm{K}$. The two stars travel on circular orbits around their common centre of mass in $P = 1.9751388 \pm 0.0000050$ days. We compare our results for Mon-735, along with another EB in NGC 2264 (CoRoT 223992193), to the predictions of five stellar evolution models. These suggest that the lower mass EB system Mon-735 is older than CoRoT 223992193 in the mass-radius diagram (MRD) and, to a lesser extent, in the Hertzsprung-Russell diagram (HRD). The MRD ages of Mon-735 and CoRoT 223992193 are $\sim$7-9 and 4-6 Myr, respectively, with the two components in each EB system possessing consistent ages., Comment: 18 pages, 14 figures, 5 tables, accepted for publication in MNRAS

Exoplanet emission spectra are often modelled assuming that the hemisphere observed is well represented by a horizontally homogenised atmosphere. However this approximation will likely fail for planets with a large temperature contrast in the James Webb Space Telescope (JWST) era, potentially leading to erroneous interpretations of spectra. We first develop an analytic formulation to quantify the signal-to-noise ratio and wavelength coverage necessary to disentangle temperature inhomogeneities from a hemispherically averaged spectrum. We find that for a given signal-to-noise ratio, observations at shorter wavelengths are better at detecting the presence of inhomogeneities. We then determine why the presence of an inhomogeneous thermal structure can lead to spurious molecular detections when assuming a fully homogenised planet in the retrieval process. Finally, we quantify more precisely the potential biases by modelling a suite of hot Jupiter spectra, varying the spatial contributions of a hot and a cold region, as would be observed by the different instruments of JWST/NIRSpec. We then retrieve the abundances and temperature profiles from the synthetic observations. We find that in most cases, assuming a homogeneous thermal structure when retrieving the atmospheric chemistry leads to biased results, and spurious molecular detection. Explicitly modelling the data using two profiles avoids these biases, and is statistically supported provided the wavelength coverage is wide enough, and crucially also spanning shorter wavelengths. For the high contrast used here, a single profile with a dilution factor performs as well as the two-profile case, with only one additional parameter compared to the 1-D approach., Comment: Accepted for publication by MNRAS

Doppler spectroscopy is a powerful tool for discovering and characterizing exoplanets. For decades, the standard approach to extracting radial velocities (RVs) has been to cross-correlate observed spectra with a weighted template mask. While still widely used, this approach is known to suffer numerous drawbacks, and so in recent years increasing attention has been paid to developing new and improved ways of extracting RVs. In this proof-of-concept paper we present a simple yet powerful approach to RV extraction. We use Gaussian processes to model and align all pairs of spectra with each other; without constructing a template, we combine pairwise RV shifts to produce accurate differential stellar RVs. Doing this on a highly-localized basis enables a data-driven approach to identifying and mitigating spectral contamination, even without the input of any prior astrophysical knowledge. We show that a crude implementation of this method applied to an inactive standard star yields RVs with comparable precision to and significantly lower rms variation than RVs from industry-standard pipelines. Though amenable to numerous improvements, even in its basic form presented here our method could facilitate the study of smaller planets around a wider variety of stars than has previously been possible., Comment: 27 pages, 16 figures. Accepted for publication in MNRAS

We present a novel retrieval technique that attempts to model phase curve observations of exoplanets more realistically and reliably, which we call the 2.5-dimension (2.5-D) approach. In our 2.5-D approach we retrieve the vertical temperature profile and mean gaseous abundance of a planet at all longitudes and latitudes \textbf{simultaneously}, assuming that the temperature or composition, $x$, at a particular longitude and latitude $(\Lambda,\Phi)$ is given by $x(\Lambda,\Phi) = \bar{x} + (x(\Lambda,0) - \bar{x})\cos^n\Phi$, where $\bar{x}$ is the mean of the morning and evening terminator values of $x(\Lambda,0)$, and $n$ is an assumed coefficient. We compare our new 2.5-D scheme with the more traditional 1-D approach, which assumes the same temperature profile and gaseous abundances at all points on the visible disc of a planet for each individual phase observation, using a set of synthetic phase curves generated from a GCM-based simulation. We find that our 2.5-D model fits these data more realistically than the 1-D approach, confining the hotter regions of the planet more closely to the dayside. We then apply both models to WASP-43b phase curve observations of HST/WFC3 and Spitzer/IRAC. We find that the dayside of WASP-43b is apparently much hotter than the nightside and show that this could be explained by the presence of a thick cloud on the nightside with a cloud top at pressure $< 0.2$ bar. We further show that while the mole fraction of water vapour is reasonably well constrained to $(1-10)\times10^{-4}$, the abundance of CO is very difficult to constrain with these data since it is degenerate with temperature and prone to possible systematic radiometric differences between the HST/WFC3 and Spitzer/IRAC observations. Hence, it is difficult to reliably constrain C/O., Comment: Submitted to M.N.R.A.S. Revised Submission 9th December 2019

Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) with M_* = 0.892 +/- (0.051, 0.035) M_sun and R_* = 0.8717 +/- (0.0064, 0.0061) R_sun. Kepler-538b is a 2.215 +/- (0.040, 0.034) R_earth sub-Neptune with a period of P = 81.73778 +/- 0.00013 d. It is the only known planet in the system. We collected radial velocity (RV) observations with HIRES on Keck I and HARPS-N on the TNG. We characterized stellar activity by a Gaussian process with a quasi-periodic kernel applied to our RV and cross correlation function full width at half maximum (FWHM) observations. By simultaneously modeling Kepler photometry, RV, and FWHM observations, we found a semi-amplitude of K = 1.68 +/- (0.39, 0.38) m s^-1 and a planet mass of M_p = 10.6 +/- (2.5, 2.4) M_earth. Kepler-538b is the smallest planet beyond P = 50 d with an RV mass measurement. The planet likely consists of a significant fraction of ices (dominated by water ice), in addition to rocks/metals, and a small amount of gas. Sophisticated modeling techniques such as those used in this paper, combined with future spectrographs with ultra high-precision and stability will be vital for yielding more mass measurements in this poorly understood exoplanet regime. This in turn will improve our understanding of the relationship between planet composition and insolation flux and how the rocky to gaseous transition depends on planetary equilibrium temperature., Comment: Published in AJ, 17 pages, 8 figures, 3 tables

While the Kepler Mission was designed to look at tens of thousands of faint stars (V > 12), brighter stars that saturated the detector are important because they can be and have been observed very accurately by other instruments. By analyzing the unsaturated scattered-light `halo' around these stars, we have retrieved precise light curves of most of the brightest stars in K2 fields from Campaign~4 onwards. The halo method does not depend on the detailed cause and form of systematics, and we show that it is effective at extracting light curves from both normal and saturated stars. The key methodology is to optimize the weights of a linear combination of pixel time series with respect to an objective function. We test a range of such objective functions, finding that lagged Total Variation, a generalization of Total Variation, performs well on both saturated and unsaturated K2 targets. Applying this to the bright stars across the K2 Campaigns reveals stellar variability ubiquitously, including effects of stellar pulsation, rotation, and binarity. We describe our pipeline and present a catalogue of the 161 bright stars, with classifications of their variability, asteroseismic parameters for red giants with well-measured solar-like oscillations, and remarks on interesting objects. These light curves are publicly available as a High Level Science Product from the Mikulski Archive for Space Telescopes (MAST)., Comment: Accepted ApJS. Paper documents data release of K2 Halo light curves, available from http://archive.stsci.edu/doi/resolve/resolve.html?doi=10.17909/t9-6wj4-eb32 or https://github.com/benjaminpope/k2halo. v2: corrected typo

We present the first data release of the Kepler Smear Campaign, using collateral 'smear' data obtained in the Kepler four-year mission to reconstruct light curves of 102 stars too bright to have been otherwise targeted. We describe the pipeline developed to extract and calibrate these light curves, and show that we attain photometric precision comparable to stars analyzed by the standard pipeline in the nominal Kepler mission. In this paper, aside from publishing the light curves of these stars, we focus on 66 red giants for which we detect solar-like oscillations, characterizing 33 of these in detail with spectroscopic chemical abundances and asteroseismic masses as benchmark stars. We also classify the whole sample, finding nearly all to be variable, with classical pulsations and binary effects. All source code, light curves, TRES spectra, and asteroseismic and stellar parameters are publicly available as a Kepler legacy sample., Comment: 35 pages, accepted ApJS

The open cluster M67 offers the unique opportunity to measure rotation periods for solar-age stars across a range of masses, potentially filling a critical gap in the understanding of angular momentum loss in older main sequence stars. The observation of M67 by NASA K2 Campaign 5 provided light curves with high enough precision to make this task possible, albeit challenging, as the pointing instability, 75d observation window, crowded field, and typically low-amplitude signals mean determining accurate rotation periods on the order of 25 - 30d is inherently difficult. Lingering, non-astrophysical signals with power at >25d found in a set of Campaign 5 A and F stars compounds the problem. To achieve a quantitative understanding of the best-case scenario limits for reliable period detection imposed by these inconveniences, we embarked on a comprehensive set of injection tests, injecting 120,000 sinusoidal signals with periods ranging from 5 to 35d and amplitudes from 0.05% to 3.0% into real Campaign 5 M67 light curves processed using two different pipelines. We attempted to recover the signals using a normalized version of the Lomb-Scargle periodogram and setting a detection threshold. We find that while the reliability of detected periods is high, the completeness (sensitivity) drops rapidly with increasing period and decreasing amplitude, maxing at 15% recovery rate for the solar case (i.e. 25d period, 0.1% amplitude). This study highlights the need for caution in determining M67 rotation periods from Campaign 5 data, but this can be extended to other clusters observed by K2 and, soon, TESS., Comment: 55 pages, 22 figures, published by ApJ

While brown dwarfs show similarities with stars in their early life, their spin evolution is much more akin to that of planets. We have used lightcurves from the K2 mission to measure new rotation periods for 18 young brown dwarfs in the Taurus star-forming region. Our sample spans masses from 0.02 to 0.08 Msol and has been characterised extensively in the past. To search for periods, we utilize three different methods (autocorrelation, periodogram, Gaussian Processes). The median period for brown dwarfs with disks is twice as long as for those without (3.1 vs. 1.6 d), a signature of rotational braking by the disk, albeit with small numbers. With an overall median period of 1.9 d, brown dwarfs in Taurus rotate slower than their counterparts in somewhat older (3-10 Myr) star-forming regions, consistent with spin-up of the latter due to contraction and angular momentum conservation, a clear sign that disk braking overall is inefficient and/or temporary in this mass domain. We confirm the presence of a linear increase of the typical rotation period as a function of mass in the sub-stellar regime. The rotational velocities, when calculated forward to the age of the solar system assuming angular momentum conservation, fit the known spin-mass relation for solar system planets and extra-solar planetary-mass objects. This spin-mass trend holds over six orders of magnitude in mass, including objects from several different formation paths. Our result implies that brown dwarfs by and large retain their primordial angular momentum through the first few Myr of their evolution., Comment: 20 pages, 6 figures, accepted for publication in ApJ

We present an analysis of a slightly eccentric ($e=0.05$), partially eclipsing long-period ($P = 69.73$ d) main sequence binary system (WOCS 12009, Sanders 1247) in the benchmark old open cluster M67. Using Kepler K2 and ground-based photometry along with a large set of new and reanalyzed spectra, we derived highly precise masses ($1.111\pm0.015$ and $0.748\pm0.005 M_\odot$) and radii ($1.071\pm0.008\pm0.003$ and $0.713\pm0.019\pm0.026 R_\odot$, with statistical and systematic error estimates) for the stars. The radius of the secondary star is in agreement with theory. The primary, however, is approximately $15\%$ smaller than reasonable isochrones for the cluster predict. Our best explanation is that the primary star was produced from the merger of two stars, as this can also account for the non-detection of photospheric lithium and its higher temperature relative to other cluster main sequence stars at the same $V$ magnitude. To understand the dynamical characteristics (low measured rotational line broadening of the primary star and the low eccentricity of the current binary orbit), we believe that the most probable (but not the only) explanation is the tidal evolution of a close binary within a primordial triple system (possibly after a period of Kozai-Lidov oscillations), leading to merger approximately 1Gyr ago. This star appears to be a future blue straggler that is being revealed as the cluster ages and the most massive main sequence stars die out., Comment: 33 pages, 11 figures, accepted to AJ, photometry files will be available with the electronic journal article

We set out to study the atmosphere of WASP-80b, a warm inflated gas giant with an equilibrium temperature of $\sim$800~K, using ground-based transmission spectroscopy covering the spectral range from 520~to~910~nm. The observations allow us to probe the existence and abundance of K and Na in WASP-80b's atmosphere, existence of high-altitude clouds, and Rayleigh-scattering in the blue end of the spectrum. We observed two spectroscopic time series of WASP-80b transits with the OSIRIS spectrograph installed in the Gran Telescopio CANARIAS, and use the observations to estimate the planet's transmission spectrum between 520~nm and 910~nm in 20~nm-wide passbands, and around the K~I and Na~I resonance doublets in 6~nm-wide passbands. We model three previously published broadband datasets consisting of 27 light curves jointly prior to the transmission spectroscopy analysis in order to obtain improved prior estimates for the planet's orbital parameters, average radius ratio, and stellar density. We recover a flat transmission spectrum with no evidence of Rayleigh scattering or K~I or Na~I absorption, and obtain an improved system characterisation as a by-product of the broadband- and GTC-dataset modelling. The transmission spectra estimated separately from the two observing runs are consistent with each other, as are the transmission spectra estimated using either a parametric or nonparametric systematics models. The flat transmission spectrum favours an atmosphere model with high-altitude clouds over cloud-free models with stellar or sub-stellar metallicities., Comment: Accepted to A&A 6.9.2017

Initial RV characterisation of the enigmatic planet Kepler-10c suggested a mass of $\sim17$ M$_\oplus$, which was remarkably high for a planet with radius $2.32$ R$_\oplus$; further observations and subsequent analysis hinted at a (possibly much) lower mass, but masses derived using RVs from two different spectrographs (HARPS-N and HIRES) were incompatible at a $3\sigma$-level. We demonstrate here how such mass discrepancies may readily arise from sub-optimal sampling and/or neglecting to model even a single coherent signal (stellar, planetary, or otherwise) that may be present in RVs. We then present a plausible resolution of the mass discrepancy, and ultimately characterise Kepler-10c as having mass $7.37_{-1.19}^{+1.32}$ M$_\oplus$, and mean density $3.14^{+0.63}_{-0.55}$ g cm$^{-3}$., Comment: 7 pages, 4 figures. Accepted for publication in MNRAS Letters

Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasi-periodic --- spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalising over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine-fitting periodogram and an autocorrelation function method. We also demonstrate that it works well on real data, by inferring rotation periods for 275 Kepler stars with previously measured periods. We provide a table of rotation periods for these 1132 Kepler objects of interest and their posterior probability density function samples. Because this method delivers posterior probability density functions, it will enable hierarchical studies involving stellar rotation, particularly those involving population modelling, such as inferring stellar ages, obliquities in exoplanet systems, or characterising star-planet interactions. The code used to implement this method is available online., Comment: Submitted to MNRAS. Replaced 27/06/2017: corrections made to koi_periods.csv

We present the discovery of four low-mass ($M<0.6$ $M_\odot$) eclipsing binary (EB) systems in the sub-Gyr old Praesepe open cluster using Kepler/K2 time-series photometry and Keck/HIRES spectroscopy. We present a new Gaussian process eclipsing binary model, GP-EBOP, as well as a method of simultaneously determining effective temperatures and distances for EBs. Three of the reported systems (AD 3814, AD 2615 and AD 1508) are detached and double-lined, and precise solutions are presented for the first two. We determine masses and radii to 1-3% precision for AD 3814 and to 5-6% for AD 2615. Together with effective temperatures determined to $\sim$50 K precision, we test the PARSEC v1.2 and BHAC15 stellar evolution models. Our EB parameters are more consistent with the PARSEC models, primarily because the BHAC15 temperature scale is hotter than our data over the mid M-dwarf mass range probed. Both ADs 3814 and 2615, which have orbital periods of 6.0 and 11.6 days, are circularized but not synchronized. This suggests that either synchronization proceeds more slowly in fully convective stars than the theory of equilibrium tides predicts or magnetic braking is currently playing a more important role than tidal forces in the spin evolution of these binaries. The fourth system (AD 3116) comprises a brown dwarf transiting a mid M-dwarf, which is the first such system discovered in a sub-Gyr open cluster. Finally, these new discoveries increase the number of characterized EBs in sub-Gyr open clusters by 20% (40%) below $M<1.5$ $M_{\odot}$ ($M<0.6$ $M_{\odot}$)., Comment: 28 pages, 22 figures, 6 tables. Accepted for publication in ApJ

We report the discovery of a repeating photometric signal from a low-mass member of the Praesepe open cluster that we interpret as a Neptune-sized transiting planet. The star is JS 183 (HSHJ 163, EPIC 211916756) with $T_{\rm eff} = 3325\pm100$ K, $M_{*} = 0.44\pm0.04$ $M_{\odot}$, $R_{*} = 0.44\pm0.03$ $R_{\odot}$, and $\log{g_*} = 4.82\pm0.06$. The planet has an orbital period of 10.134588 days and a radius of $R_{P}= 0.32\pm0.02$ $R_J$. Since the star is faint at $V=16.5$ and $J=13.3$, we are unable to obtain a measured radial-velocity orbit, but we can constrain the companion mass to below about 1.7 $M_J$, and thus well below the planetary boundary. JS 183b (since designated as K2-95b) is the second transiting planet found with ${\it K2}$ that resides in a several hundred Myr open cluster; both planets orbit mid-M dwarf stars and are approximately Neptune-sized. With a well-determined stellar density from the planetary transit, and with an independently known metallicity from its cluster membership, JS 183 provides a particularly valuable test of stellar models at the fully convective boundary. We find that JS 183 is the lowest-density transit host known at the fully convective boundary, and that its very low density is consistent with current models of stars just above the fully convective boundary but in tension with the models just below the fully convective boundary., Comment: 30 pages, 11 figures, Accepted for publication in the Astronomical Journal

CoRoT 223992193 is the only known low-mass, pre-main sequence eclipsing binary that shows evidence of a circumbinary disk. The system displays complex photometric and spectroscopic variability over a range of timescales and wavelengths. Using two optical CoRoT runs, and infrared Spitzer 3.6 and 4.5 $\mu$m observations (simultaneous with the second CoRoT run), we model the out-of-eclipse light curves. The large scale structure in both CoRoT light curves is consistent with the constructive and destructive interference of starspot signals at two slightly different periods. Using the stellar $v\sin i$ 's, we infer different rotation periods: the primary is consistent with synchronisation and the secondary is slightly supersynchronous. Comparison of the raw data to the residuals of our spot model in colour-magnitude space indicates additional contributions consistent with variable dust emission and obscuration. We also identify short-duration flux dips preceding secondary eclipse in all three CoRoT and Spitzer bands. We construct a model of the inner regions of the binary and propose that these dips could be caused by partial occultation of the central binary by the accretion stream onto the primary star. Analysis of 15 VLT/FLAMES H$\alpha$ profiles reveal an emission profile associated with each star: the majority is consistent with chromospheric emission but higher velocity emission is also seen, which could be due to prominences or accretion-related phenomena. In addition, simultaneous u and r-band observations from CFHT/MEGACam reveal a short-lived u-band excess consistent with either an accretion hot spot or stellar flare. The photometric and spectroscopic variations are complex but are consistent with the picture of two active stars possibly undergoing non-steady, low-level accretion; the system's very high inclination provides a new view of such variability. [Abridged], Comment: 18 pages, 14 figures, 2 tables. Figures 1 & 2 reproduced from Gillen et al. (2014) for clarity. Accepted for publication in A&A

We present a consistent optimal estimation retrieval analysis of ten hot Jupiter exoplanets, each with transmission spectral data spanning the visible to near-infrared wavelength range. Using the NEMESIS radiative transfer and retrieval tool, we calculate a range of possible atmospheric states for WASP-6b, WASP-12b, WASP-17b, WASP-19b, WASP-31b, WASP-39b, HD 189733b, HD 209458b, HAT-P-1b and HAT-P-12b. We find that the spectra of all ten planets are consistent with the presence of some atmospheric aerosol; WASP-6b, WASP-12b, WASP-17b, WASP-19b, HD 189733b and HAT-P-12b are all fit best by Rayleigh scattering aerosols, whereas WASP-31b, WASP-39b and HD 209458b are better represented by a grey cloud model. HAT-P-1b has solutions that fall into both categories. WASP-6b, HAT-P-12b, HD 189733b and WASP-12b must have aerosol extending to low atmospheric pressures (below 0.1 mbar). In general, planets with equilibrium temperatures between 1300 and 1700 K are best represented by deeper, grey cloud layers, whereas cooler or hotter planets are better fit using high Rayleigh scattering aerosol. We find little evidence for the presence of molecular absorbers other than H$_2$O. Retrieval methods can provide a consistent picture across a range of hot Jupiter atmospheres with existing data, and will be a powerful tool for the interpretation of James Webb Space Telescope observations., Comment: 25 pages, 21 figures. Accepted for publication in ApJ

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 systematics-corrected, 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, WASP-55b (EPIC 212300977b) and Qatar-2b (EPIC 212756297b). We also report the outcome of reconnaissance spectroscopy carried out for all candidates with Kepler magnitude Kp < 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., Comment: Accepted to MNRAS

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 ratios down Rp/Rstar = 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., Comment: 14 pages, 13 figures, accepted for publication in MNRAS. The processed data for campaigns 3 to 6 will be available shortly at https://archive.stsci.edu/prepds/k2sc/ and the K2SC code is available at https://github.com/OxES/k2sc