Source author record

Suzanne Aigrain

Suzanne Aigrain appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2026arXiv

Mass estimates of the young TOI-451 transiting planets: Multidimensional Gaussian Process on stellar spectroscopic and photometric signals

The young TOI-451 planetary system, aged 125 Myr, provides a unique opportunity to test theories of planetary internal structures and atmospheric mass loss through examination of its three transiting planets. We present an exhaustive photometric and spectroscopic follow-up to determine the orbital and physical properties of the system. We perform multidimensional Gaussian Process regression with the code pyaneti on spectroscopic time-series and NGTS/LCO light curves to disentangle the stellar and planetary signal in ESPRESSO radial velocities. We show how contemporaneous photometry serves as an activity indicator to inform RV modelling within a multidimensional Gaussian Processes framework. We argue that this can be exploited when spectroscopic observations are adversely affected by low signal-to-noise and/or poor sampling. We estimate the Doppler semi-amplitudes of Kb = 2.6(+1.1,-1.2) m/s, Kc = 1.2(+1.0,-0.8) m/s and Kd = 2.7 +/- 1.2 m/s. This translates into 2-sigma mass estimates for TOI-451 b and d of Mb = 4.7(+2.1,-2.2) Earth masses and Md = 10.2(+4.6,-4.5) Earth masses, as well as a mass upper limit for TOI-451 c of Mc < 11.5 Earth masses. The derived planetary properties suggest that planets c and d contain significant hydrogen-rich envelopes. The inferred parameters of TOI-451 b are consistent with either a rocky world that still retains a small hydrogen envelope or a water world. These insights make the TOI-451 system an ideal laboratory for future follow-up studies aimed at measuring atmospheric compositions, detecting atmospheric mass-loss signatures, and further exploring planetary formation and evolution processes.

preprint2022arXiv

Applications of a Gaussian Process Framework for Modelling of High-Resolution Exoplanet Spectra

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 $ΔK_\mathrm{p} =1.45 \pm 1.53\,\mathrm{km\,s^{-1}}$ and $Δ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.

preprint2022arXiv

One year of AU Mic with HARPS: I -- measuring the masses of the two transiting planets

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.

preprint2022arXiv

One year of AU Mic with HARPS: II -- stellar activity and star-planet interaction

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-$σ$ 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.

preprint2022arXiv

Planet Hunters TESS IV: A massive, compact hierarchical triple star system TIC 470710327

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.

preprint2022arXiv

Quasi-periodic Gaussian Processes for stellar activity: from physical to kernel parameters

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.

preprint2022arXiv

The EXPRES Stellar Signals Project II. State of the Field in Disentangling Photospheric Velocities

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.

preprint2020arXiv

Mon-735: A new low-mass pre-main sequence eclipsing binary in NGC 2264

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.

preprint2020arXiv

Pleiades or Not? Resolving the Status of the Lithium Rich M Dwarfs HHJ339 and HHJ430

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.

preprint2019arXiv

A robust, template-free approach to precise radial velocity extraction

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.

preprint2016arXiv

A consistent retrieval analysis of 10 Hot Jupiters observed in transmission

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.

preprint2016arXiv

HST hot-Jupiter transmission spectral survey: Clear skies for cool Saturn WASP-39b

We present HST STIS optical transmission spectroscopy of the cool Saturn-mass exoplanet WASP-39b from 0.29-1.025 micron, along with complementary transit observations from Spitzer IRAC at 3.6 and 4.5 micron. The low density and large atmospheric pressure scale height ofWASP-39b make it particularly amenable to atmospheric characterization using this technique. We detect a Rayleigh scattering slope as well as sodium and potassium absorption features; this is the first exoplanet in which both alkali features are clearly detected with the extended wings predicted by cloud-free atmosphere models. The full transmission spectrum is well matched by a clear, H2-dominated atmosphere or one containing a weak contribution from haze, in good agreement with the preliminary reduction of these data presented in Sing et al. (2016). WASP-39b is predicted to have a pressure-temperature profile comparable to that of HD 189733b and WASP-6b, making it one of the coolest transiting gas giants observed in our HST STIS survey. Despite this similarity, WASP-39b appears to be largely cloud-free while the transmission spectra of HD 189733b and WASP-6b both indicate the presence of high altitude clouds or hazes. These observations further emphasize the surprising diversity of cloudy and cloud-free gas giant planets in short-period orbits and the corresponding challenges associated with developing predictive cloud models for these atmospheres.

preprint2016arXiv

K2SC: Flexible systematics correction and detrending of K2 light curves using Gaussian Process regression

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.

preprint2016arXiv

New Pleiades Eclipsing Binaries and a Hyades Transiting System Identified by K2

We present the discovery in Kepler's $K2$ mission observations and our follow-up radial velocity observations from Keck/HIRES for four eclipsing binary (EB) star systems in the young benchmark Pleiades cluster. Based on our modeling results, we announce two new low mass ($M_{tot} < 0.6 M_\odot$) EBs among Pleiades members (HCG 76 and MHO 9) and we report on two previously known Pleiades binaries that are also found to be EB systems (HII 2407 and HD 23642). We measured the masses of the binary HCG 76 to $\lesssim$2.5% precision, and the radii to $\lesssim$4.5% precision, which together with the precise effective temperatures yield an independent Pleiades distance of 132$\pm$5 pc. We discuss another EB towards the Pleiades that is a possible but unlikely Pleiades cluster member (AK II 465). The two new confirmed Pleiades systems extend the mass range of Pleiades EB components to 0.2-2 $M_\odot$. Our initial measurements of the fundamental stellar parameters for the Pleiades EBs are discussed in the context of the current stellar models and the nominal cluster isochrone, finding good agreement with the stellar models of Baraffe et al. (2015) at the nominal Pleiades age of 120 Myr. Finally, in the Hyades, we report a new low mass eclipsing system (vA 50) that was concurrently discovered and studied by Mann et al. (2016). We confirm that the eclipse is likely caused by a Neptune-sized transiting planet, and with the additional radial velocity constraints presented here we improve the constraint on the maximum mass of the planet to be $\lesssim$1.2 $M_\mathrm{Jup}$.

preprint2016arXiv

Telling twins apart: Exo-Earths and Venuses with transit spectroscopy

The planned launch of the James Webb Space Telescope in 2018 will herald a new era of exoplanet spectroscopy. JWST will be the first telescope sensitive enough to potentially characterize terrestrial planets from their transmission spectra. In this work, we explore the possibility that terrestrial planets with Venus-type and Earth-type atmospheres could be distinguished from each other using spectra obtained by JWST. If we find a terrestrial planet close to the liquid water habitable zone of an M5 star within a distance of 10 parsecs, it would be possible to detect atmospheric ozone if present in large enough quantities, which would enable an oxygen-rich atmosphere to be identified. However, the cloudiness of a Venus-type atmosphere would inhibit our ability to draw firm conclusions about the atmospheric composition, making any result ambiguous. Observing small, temperate planets with JWST requires significant investment of resources, with single targets requiring of order 100 transits to achieve sufficient signal to noise. The possibility of detecting a crucial feature such as the ozone signature would need to be carefully weighed against the likelihood of clouds obscuring gas absorption in the spectrum.

preprint2016arXiv

Transiting exoplanet candidates from K2 Campaigns 5 and 6

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.

preprint2015arXiv

A continuum from clear to cloudy hot-Jupiter exoplanets without primordial water depletion

Thousands of transiting exoplanets have been discovered, but spectral analysis of their atmospheres has so far been dominated by a small number of exoplanets and data spanning relatively narrow wavelength ranges (such as 1.1 to 1.7 μm). Recent studies show that some hot-Jupiter exoplanets have much weaker water absorption features in their near-infrared spectra than predicted. The low amplitude of water signatures could be explained by very low water abundances, which may be a sign that water was depleted in the protoplanetary disk at the planet's formation location, but it is unclear whether this level of depletion can actually occur. Alternatively, these weak signals could be the result of obscuration by clouds or hazes, as found in some optical spectra. Here we report results from a comparative study of ten hot Jupiters covering the wavelength range 0.3-5 micrometres, which allows us to resolve both the optical scattering and infrared molecular absorption spectroscopically. Our results reveal a diverse group of hot Jupiters that exhibit a continuum from clear to cloudy atmospheres. We find that the difference between the planetary radius measured at optical and infrared wavelengths is an effective metric for distinguishing different atmosphere types. The difference correlates with the spectral strength of water, so that strong water absorption lines are seen in clear-atmosphere planets and the weakest features are associated with clouds and hazes. This result strongly suggests that primordial water depletion during formation is unlikely and that clouds and hazes are the cause of weaker spectral signatures.

preprint2015arXiv

A Gaussian process framework for modelling stellar activity signals in radial velocity data

To date, the radial velocity (RV) method has been one of the most productive techniques for detecting and confirming extrasolar planetary candidates. Unfortunately, stellar activity can induce RV variations which can drown out or even mimic planetary signals - and it is notoriously difficult to model and thus mitigate the effects of these activity-induced nuisance signals. This is expected to be a major obstacle to using next-generation spectrographs to detect lower mass planets, planets with longer periods, and planets around more active stars. Enter Gaussian processes (GPs) which, we note, have a number of attractive features that make them very well suited to disentangling stellar activity signals from planetary signals. We present here a GP framework we developed to model RV time series jointly with ancillary activity indicators (e.g. bisector velocity spans, line widths, chromospheric activity indices), allowing the activity component of RV time series to be constrained and disentangled from e.g. planetary components. We discuss the mathematical details of our GP framework, and present results illustrating its encouraging performance on both synthetic and real RV datasets, including the publicly-available Alpha Centauri B dataset.

preprint2015arXiv

A uniform analysis of HD209458b Spitzer/IRAC lightcurves with Gaussian process models

We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3.6 micron and 4.5 micron channels and the ramp effect in the 5.8 micron and 8.0 micron channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 micron channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of $1484\pm 18$ K.

preprint2015arXiv

Calibrating Gyrochronology using Kepler Asteroseismic targets

Among the available methods for dating stars, gyrochronology is a powerful one because it requires knowledge of only the star's mass and rotation period. Gyrochronology relations have previously been calibrated using young clusters, with the Sun providing the only age dependence, and are therefore poorly calibrated at late ages. We used rotation period measurements of 310 Kepler stars with asteroseismic ages, 50 stars from the Hyades and Coma Berenices clusters and 6 field stars (including the Sun) with precise age measurements to calibrate the gyrochronology relation, whilst fully accounting for measurement uncertainties in all observable quantities. We calibrated a relation of the form $P=A^n\times(B-V-c)^b$, where $P$ is rotation period in days, $A$ is age in Myr, $B$ and $V$ are magnitudes and $a$, $b$ and $n$ are the free parameters of our model. We found $a = 0.40^{+0.3}_{-0.05}$, $b = 0.31^{+0.05}_{-0.02}$ and $n = 0.55^{+0.02}_{-0.09}$. Markov Chain Monte Carlo methods were used to explore the posterior probability distribution functions of the gyrochronology parameters and we carefully checked the effects of leaving out parts of our sample, leading us to find that no single relation beween rotation period, colour and age can adequately describe all the subsets of our data. The Kepler asteroseismic stars, cluster stars and local field stars cannot all be described by the same gyrochronology relation. The Kepler asteroseismic stars may be subject to observational biases, however the clusters show unexpected deviations from the predicted behaviour, providing concerns for the overall reliability of gyrochronology as a dating method.

preprint2015arXiv

CSI 2264: Probing the inner disks of AA Tau-like systems in NGC 2264

The classical T Tauri star AA Tau presented photometric variability attributed to an inner disk warp, caused by the interaction between the inner disk and an inclined magnetosphere. Previous studies of NGC 2264 have shown that similar photometric behavior is common among CTTS. The goal of this work is to investigate the main causes of the observed photometric variability of CTTS in NGC 2264 that present AA Tau-like light curves, and verify if an inner disk warp could be responsible for their variability. We investigate veiling variability in their spectra and u-r color variations and estimate parameters of the inner disk warp using an occultation model proposed for AA Tau. We compare infrared and optical light curves to analyze the dust responsible for the occultations. AA Tau-like variability is transient on a timescale of a few years. We ascribe it to stable accretion regimes and aperiodic variability to unstable accretion regimes and show that a transition, and even coexistence, between the two is common. We find evidence of hot spots associated with occultations, indicating that the occulting structures could be located at the base of accretion columns. We find average values of warp maximum height of 0.23 times its radial location, consistent with AA Tau, with variations of on average 11% between rotation cycles. We show that extinction laws in the inner disk indicate the presence of grains larger than interstellar grains. The inner disk warp scenario is consistent with observations for all but one periodic star in our sample. AA Tau-like systems comprise 14% of CTTS observed in NGC 2264, though this increases to 35% among systems of mass 0.7M_sun<M<2.0M_sun. Assuming random inclinations, we estimate that nearly all systems in this mass range likely possess an inner disk warp, possibly because of a change in magnetic field configurations among stars of lower mass.

preprint2015arXiv

Exoplanet Transmission Spectroscopy using KMOS

KMOS (K-Band Multi Object Spectrograph) is a novel integral field spectrograph installed in the VLT's ANTU unit. The instrument offers an ability to observe 24 2.8"$\times$2.8" sub-fields positionable within a 7.2' patrol field, each sub-field producing a spectrum with a 14$\times$14-pixel spatial resolution. The main science drivers for KMOS are the study of galaxies, star formation, and molecular clouds, but its ability to simultaneously measure spectra of multiple stars makes KMOS an interesting instrument for exoplanet atmosphere characterization via transmission spectroscopy. We set to test whether transmission spectroscopy is practical with KMOS, and what are the conditions required to achieve the photometric precision needed, based on observations of a partial transit of WASP-19b, and full transits of GJ 1214b and HD 209458b. Our analysis uses the simultaneously observed comparison stars to reduce the effects from instrumental and atmospheric sources, and Gaussian processes to model the residual systematics. We show that KMOS can, in theory, deliver the photometric precision required for transmission spectroscopy. However, this is shown to require a) pre-imaging to ensure accurate centering and b) a very stable night with optimal observing conditions (seeing $\sim$0.8"). Combining these two factors with the need to observe several transits, each with a sufficient out-of-transit baseline (and with the fact that similar or better precision can be reached with telescopes and instruments with smaller pressure,) we conclude that transmission spectroscopy is not the optimal science case to take advantage of the abilities offered by KMOS and VLT.

preprint2015arXiv

From Cosmic Birth to Living Earths: The Future of UVOIR Space Astronomy

For the first time in history, humans have reached the point where it is possible to construct a revolutionary space-based observatory that has the capability to find dozens of Earth-like worlds, and possibly some with signs of life. This same telescope, designed as a long-lived facility, would also produce transformational scientific advances in every area of astronomy and astrophysics from black hole physics to galaxy formation, from star and planet formation to the origins of the Solar System. The Association of Universities for Research in Astronomy (AURA) commissioned a study on a next-generation UVOIR space observatory with the highest possible scientific impact in the era following JWST. This community-based study focuses on the future space-based options for UV and optical astronomy that significantly advance our understanding of the origin and evolution of the cosmos and the life within it. The committee concludes that a space telescope equipped with a 12-meter class primary mirror can find and characterize dozens of Earth-like planets and make fundamental advances across nearly all fields of astrophysics. The concept is called the High Definition Space Telescope (HDST). The telescope would be located at the Sun-Earth L2 point and would cover a spectral range that, at a minimum, runs from 0.1 to 2 microns. Unlike JWST, HDST will not need to operate at cryogenic temperatures. HDST can be made to be serviceable on orbit but does not require servicing to complete its primary scientific objectives. We present the scientific and technical requirements for HDST and show that it could allow us to determine whether or not life is common outside the Solar System. We do not propose a specific design for such a telescope, but show that designing, building and funding such a facility is feasible beginning in the next decade - if the necessary strategic investments in technology begin now.

preprint2015arXiv

Ghost in the time series: no planet for Alpha Cen B

We re-analyse the publicly available radial velocity (RV) measurements for Alpha Cen B, a star hosting an Earth-mass planet candidate, Alpha Cen Bb, with 3.24 day orbital period. We demonstrate that the 3.24 d signal observed in the Alpha Cen B data almost certainly arises from the window function (time sampling) of the original data. We show that when stellar activity signals are removed from the RV variations, other significant peaks in the power spectrum of the window function are coincidentally suppressed, leaving behind a spurious yet apparently-significant 'ghost' of a signal that was present in the window function's power spectrum to begin with. Even when fitting synthetic data with time sampling identical to the original data, but devoid of any genuine periodicities close to that of the planet candidate, the original model used to infer the presence of Alpha Cen Bb leads to identical conclusions: viz., the 3$σ$ detection of a half-a-metre-per-second signal with 3.236 day period. Our analysis underscores the difficulty of detecting weak planetary signals in RV data, and the importance of understanding in detail how every component of an RV data set, including its time sampling, influences final statistical inference.

preprint2015arXiv

HII 2407: A Low-Mass Eclipsing Binary Revealed by K2 Observations of the Pleiades

The star HII 2407 is a member of the relatively young Pleiades star cluster and was previously discovered to be a single-lined spectroscopic binary. It is newly identified here within $Kepler$/$K2$ photometric time series data as an eclipsing binary system. Mutual fitting of the radial velocity and photometric data leads to an orbital solution and constraints on fundamental stellar parameters. While the primary has arrived on the main sequence, the secondary is still pre-main-sequence and we compare our results for the $M/M_\odot$ and $R/R_\odot$ values with stellar evolutionary models. We also demonstrate that the system is likely to be tidally synchronized. Follow-up infrared spectroscopy is likely to reveal the lines of the secondary, allowing for dynamically measured masses and elevating the system to benchmark eclipsing binary status.

preprint2015arXiv

LDTk: Limb Darkening Toolkit

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.

preprint2015arXiv

Photometry of Very Bright Stars with Kepler and K2 Smear Data

High-precision time series photometry with the Kepler satellite has been crucial to our understanding both of exoplanets, and via asteroseismology, of stellar physics. After the failure of two reaction wheels, the Kepler satellite has been repurposed as Kepler-2 (K2), observing fields close to the ecliptic plane. As these fields contain many more bright stars than the original Kepler field, K2 provides an unprecedented opportunity to study nearby objects amenable to detailed follow-up with ground-based instruments. Due to bandwidth constraints, only a small fraction of pixels can be downloaded, with the result that most bright stars which saturate the detector are not observed. We show that engineering data acquired for photometric calibration, consisting of collateral `smear' measurements, can be used to reconstruct light curves for bright targets not otherwise observable with Kepler/K2. Here we present some examples from Kepler Quarter 6 and K2 Campaign 3, including the delta Scuti variables HD 178875 and 70 Aqr, and the red giant HR 8500 displaying solar-like oscillations. We compare aperture and smear photometry where possible, and also study targets not previously observed. These encouraging results suggest this new method can be applied to most Kepler and K2 fields.

preprint2015arXiv

The GTC exoplanet transit spectroscopy survey II: An overly-large Rayleigh-like feature for exoplanet TrES-3b

We set to search for Rayleigh scattering and K and Na absorption signatures from the atmosphere of TrES-3b using ground-based transmission spectroscopy covering the wavelength range from 530 to 950 nm as observed with OSIRIS@GTC. Our analysis is based on a Bayesian approach where the light curves covering a set of given passbands are fitted jointly with PHOENIX-calculated stellar limb darkening profiles. The analysis is carried out assuming both white and red -- temporally correlated -- noise, with two approaches (Gaussian processes and divide-by-white) to account for the red noise. An initial analysis reveals a transmission spectrum that shows a strong Rayleigh-like increase in extinction towards the blue end of the spectrum, and enhanced extinction around the K I resonance doublet near 767 nm. However, the signal amplitudes are significantly larger than expected from theoretical considerations. A detailed analysis reveals that the K I-like feature is entirely due to variability in the telluric O$_2$ absorption, but the Rayleigh-like feature remains unexplained.

preprint2015arXiv

Transit spectroscopy with JWST: Systematics, starspots and stitching

The James Webb Space Telescope (JWST) is predicted to make great advances in the field of exoplanet atmospheres. Its 25 m2 mirror means that it can reach unprecedented levels of precision in observations of transit spectra, and can thus characterise the atmospheres of planets orbiting stars several hundred pc away. Its coverage of the infrared spectral region between 0.6 and 28 μm allows the abundances of key molecules to be probed during the transit of a planet in front of the host star, and when the same planet is eclipsed constraints can be placed on its temperature structure. In this work, we explore the possibility of using low-spectral-resolution observations by JWST/NIRSpec and JWST/MIRI-LRS together to optimise wavelength coverage and break degeneracies in the atmospheric retrieval problem for a range of exoplanets from hot Jupiters to super Earths. This approach involves stitching together non-simultaneous observations in different wavelength regions, rendering it necessary to consider the effect of time-varying instrumental and astrophysical systematics. We present the results of a series of retrieval feasibility tests examining the effects of instrument systematics and star spots on the recoverability of the true atmospheric state, and demonstrate that correcting for these systematics is key for successful exoplanet science with JWST.

preprint2014arXiv

Clouds on the hot Jupiter HD189733b: constraints from the reflection spectrum

The hot Jupiter HD 189733b is probably the best studied of the known extrasolar planets, with published transit and eclipse spectra covering the near UV to mid-IR range. Recent work on the transmission spectrum has shown clear evidence for the presence of clouds in its atmosphere, which significantly increases the model atmosphere parameter space that must be explored in order to fully characterise this planet. In this work, we apply the NEMESIS atmospheric retrieval code to the recently published HST/STIS reflection spectrum, and also to the dayside thermal emission spectrum in the light of new Spitzer/IRAC measurements, as well as our own re-analysis of the HST/NICMOS data. We first use the STIS data to place some constraints on the nature of the cloud on HD189733b, and explore solution degeneracy between different cloud properties and the abundance of Na in the atmosphere; as already noted in previous work, absorption due to Na plays a significant role in determining the shape of the reflection spectrum. We then perform a new retrieval of the temperature profile and abundances of H2O, CO2, CO and CH4 from the dayside thermal emission spectrum. Finally, we investigate the effect of including cloud in the model on this retrieval process. We find that the current quality of data does not warrant the extra complexity introduced by including cloud in the model; however, future data are likely to be of sufficient resolution and signal-to-noise that a more complete model, including scattering particles, will be required.

preprint2014arXiv

CSI 2264: Characterizing Accretion-Burst Dominated Light Curves for Young Stars in NGC 2264

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

preprint2014arXiv

CSI 2264: Simultaneous optical and infrared light curves of young disk-bearing stars in NGC 2264 with CoRoT and Spitzer-- evidence for multiple origins of variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$α$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

preprint2014arXiv

Exoplanet atmospheres with EChO: spectral retrievals using EChOSim

We demonstrate the effectiveness of the Exoplanet Characterisation Observatory mission concept for constraining the atmospheric properties of hot and warm gas giants and super Earths. Synthetic primary and secondary transit spectra for a range of planets are passed through EChOSim (Waldmann & Pascale 2014) to obtain the expected level of noise for different observational scenarios; these are then used as inputs for the NEMESIS atmospheric retrieval code and the retrieved atmospheric properties (temperature structure, composition and cloud properties) compared with the known input values, following the method of Barstow et al. (2013a). To correctly retrieve the temperature structure and composition of the atmosphere to within 2 σ, we find that we require: a single transit or eclipse of a hot Jupiter orbiting a sun-like (G2) star at 35 pc to constrain the terminator and dayside atmospheres; 20 transits or eclipses of a warm Jupiter orbiting a similar star; 10 transits/eclipses of a hot Neptune orbiting an M dwarf at 6 pc; and 30 transits or eclipses of a GJ1214b-like planet.

preprint2014arXiv

Exploring the Diversity of Jupiter-Class Planets (Discussion Meeting Contribution)

Royal Society Discussion Meeting (2013) `Characterizing exoplanets'. Of the 900+ confirmed exoplanets discovered since 1995 for which we have constraints on their mass (i.e., not including Kepler candidates), 75% have masses larger than Saturn (0.3MJ), 53% are more massive than Jupiter, and 67% are within 1 AU of their host stars. And yet the term `hot Jupiter' fails to account for the incredible diversity of this class of object, which exists on a continuum of giant planets from the cool jovians of our own solar system to the highly-irradiated, tidally-locked hot roasters. We review theoretical expectations for the temperatures, molecular composition and cloud properties of Jupiter-class objects under a variety of different conditions. We discuss the classification schemes for these Jupiter-class planets proposed to date, including the implications for our own Solar System giant planets and the pitfalls associated with classification at this early stage of exoplanetary spectroscopy. We discuss the range of planetary types described by previous authors, accounting for: (i) thermochemical equilibrium expectations for cloud condensation and favoured chemical stability fields; (ii) the metallicity and formation mechanism for these giant planets; (iii) the importance of optical absorbers for energy partitioning and the generation of a temperature inversion; (iv) the favoured photochemical pathways and expectations for minor species (e.g., saturated hydrocarbons and nitriles); (v) the unexpected presence of molecules due to vertical mixing of species above their quench levels; and (vi) methods for energy and material redistribution throughout the atmosphere (e.g., away from the highly irradiated daysides of close-in giants). Finally, we will discuss the benefits and flaws of retrieval techniques for establishing a family of atmospheric solutions that reproduce the available data.

preprint2014arXiv

Precise time-series photometry for the Kepler-2.0 mission

The recently approved NASA K2 mission has the potential to multiply by an order of magnitude the number of short-period transiting planets found by Kepler around bright and low-mass stars, and to revolutionise our understanding of stellar variability in open clusters. However, the data processing is made more challenging by the reduced pointing accuracy of the satellite, which has only two functioning reaction wheels. We present a new method to extract precise light curves from K2 data, combining list-driven, soft-edged aperture photometry with a star-by-star correction of systematic effects associated with the drift in the roll-angle of the satellite about its boresight. The systematics are modelled simultaneously with the stars' intrinsic variability using a semi-parametric Gaussian process model. We test this method on a week of data collected during an engineering test in January 2014, perform checks to verify that our method does not alter intrinsic variability signals, and compute the precision as a function of magnitude on long-cadence (30-min) and planetary transit (2.5-hour) timescales. In both cases, we reach photometric precisions close to the precision reached during the nominal Kepler mission for stars fainter than 12th magnitude, and between 40 and 80 parts per million for brighter stars. These results confirm the bright prospects for planet detection and characterisation, asteroseismology and stellar variability studies with K2. Finally, we perform a basic transit search on the light curves, detecting 2 bona fide transit-like events, 7 detached eclipsing binaries and 13 classical variables.

preprint2014arXiv

Rotation Periods of 34,030 Kepler Main-Sequence Stars: The Full Autocorrelation Sample

We analyzed 3 years of data from the Kepler space mission to derive rotation periods of main-sequence stars below 6500 K. Our automated autocorrelation-based method detected rotation periods between 0.2 and 70 days for 34,030 (25.6%) of the 133,030 main-sequence Kepler targets (excluding known eclipsing binaries and Kepler Objects of Interest), making this the largest sample of stellar rotation periods to date. In this paper we consider the detailed features of the now well-populated period-temperature distribution and demonstrate that the period bimodality, first seen by McQuillan, Aigrain & Mazeh (2013) in the M-dwarf sample, persists to higher masses, becoming less visible above 0.6 M_sun. We show that these results are globally consistent with the existing ground-based rotation-period data and find that the upper envelope of the period distribution is broadly consistent with a gyrochronological age of 4.5 Gyrs, based on the isochrones of Barnes (2007), Mamajek & Hillenbrand (2008) and Meibom et al. (2009). We also performed a detailed comparison of our results to those of Reinhold et al. (2013) and Nielsen et al. (2013), who have measured rotation periods of field stars observed by Kepler. We examined the amplitude of periodic variability for the stars with detected rotation periods, and found a typical range between ~950 ppm (5th percentile) and ~22,700 ppm (95th percentile), with a median of ~5,600 ppm. We found typically higher amplitudes for shorter periods and lower effective temperatures, with an excess of low-amplitude stars above ~5400 K.

preprint2014arXiv

The K2 Mission: Characterization and Early results

The K2 mission will make use of the Kepler spacecraft and its assets to expand upon Kepler's groundbreaking discoveries in the fields of exoplanets and astrophysics through new and exciting observations. K2 will use an innovative way of operating the spacecraft to observe target fields along the ecliptic for the next 2-3 years. Early science commissioning observations have shown an estimated photometric precision near 400 ppm in a single 30 minute observation, and a 6-hour photometric precision of 80 ppm (both at V=12). The K2 mission offers long-term, simultaneous optical observation of thousands of objects at a precision far better than is achievable from ground-based telescopes. Ecliptic fields will be observed for approximately 75-days enabling a unique exoplanet survey which fills the gaps in duration and sensitivity between the Kepler and TESS missions, and offers pre-launch exoplanet target identification for JWST transit spectroscopy. Astrophysics observations with K2 will include studies of young open clusters, bright stars, galaxies, supernovae, and asteroseismology.

preprint2014arXiv

The Transit Spectra of Earth and Jupiter

In recent years, a number of observations have been made of the transits of 'Hot Jupiters', such as HD 189733b, which have been modelled to derive atmospheric structure and composition. As measurement techniques improve, the transit spectra of 'Super-Earths' such as GJ 1214b are becoming better constrained, allowing model atmospheres to be fitted for this class of planet also. While it is not yet possible to constrain the atmospheric states of small planets such as the Earth or cold planets like Jupiter, this may become practical in the coming decades and if so, it is of interest to determine what we might infer from such measurements. Here we have constructed atmospheric models of the Solar System planets from 0.4 - 15.5 microns that are consistent with ground-based and satellite observations and from these calculate the primary transit and secondary eclipse spectra (with respect to the Sun and typical M-dwarfs) that would be observed by a 'remote observer', many light years away. From these spectra we test what current retrieval models might infer about their atmospheres and compare these with the 'ground truths' in order to assess: a) the inherent uncertainties in transit spectra observations; b) the relative merits of primary transit and secondary eclipse spectra; and c) the advantages of directly imaged spectra. We find that secondary eclipses would not give sufficient information, but that primary transits give much better determination. We find that a single transit of Jupiter in front of the Sun could potentially be used to determine temperature and stratospheric composition, but for the Earth the mean atmospheric composition could only be determined if it were orbiting an M-dwarf. For both planets we note that direct imaging with sufficient nulling of the light from the parent star provides the best method of determining the atmospheric properties of such planets.

preprint2014arXiv

VISTA Variables in the Vía Láctea (VVV): Halfway Status and Results

The VISTA Variables in the Vía Láctea (VVV) survey is one of six public ESO surveys, and is now in its 4th year of observing. Although far from being complete, the VVV survey has already delivered many results, some directly connected to the intended science goals (detection of variables stars, microlensing events, new star clusters), others concerning more exotic objects, e.g. novae. Now, at the end of the fourth observing period, and comprising roughly 50% of the proposed observations, the actual status of the survey, as well some of the results based on the VVV data, are presented.

preprint2013arXiv

Constraining the atmosphere of GJ 1214b using an optimal estimation technique

We explore cloudy, extended H2-He atmosphere scenarios for the warm super-Earth GJ 1214b using an optimal estimation retrieval technique. This planet, orbiting an M4.5 star only 13 pc from the Earth, is of particular interest because it lies between the Earth and Neptune in size and may be a member of a new class of planet that is neither terrestrial nor gas giant. Its relatively flat transmission spectrum has so far made atmospheric characterisation difficult. The NEMESIS algorithm (Irwin et al. 2008) is used to explore the degenerate model parameter space for a cloudy, H2-He-dominated atmosphere scenario. Optimal estimation is a data-led approach that allows solutions beyond the range permitted by ab initio equilibrium model atmosphere calculations, and as such prevents restriction from prior expectations. We show that optimal estimation retrieval is a powerful tool for this kind of study, and present an exploration of the degenerate atmospheric scenarios for GJ 1214b. Whilst we find a family of solutions that provide a very good fit to the data, the quality and coverage of these data are insufficient for us to more precisely determine the abundances of cloud and trace gases given an H2-He atmosphere, and we also cannot rule out the possibility of a high molecular weight atmosphere. Future ground- and space-based observations will provide the opportunity to confirm or rule out an extended H2-He atmosphere, but more precise constraints will be limited by intrinsic degeneracies in the retrieval problem, such as variations in cloud top pressure and temperature.

preprint2013arXiv

CoRoT 223992193: A new, low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk

We present the discovery of CoRoT 223992193, a double-lined, detached eclipsing binary, comprising two pre-main sequence M dwarfs, discovered by the CoRoT space mission during a 23-day observation of the 3 Myr old NGC 2264 star-forming region. Using multi-epoch optical and near-IR follow-up spectroscopy with FLAMES on the Very Large Telescope and ISIS on the William Herschel Telescope we obtain a full orbital solution and derive the fundamental parameters of both stars by modelling the light curve and radial velocity data. The orbit is circular and has a period of $3.8745745 \pm 0.0000014$ days. The masses and radii of the two stars are $0.67 \pm 0.01$ and $0.495 \pm 0.007$ $M_{\odot}$ and $1.30 \pm 0.04$ and $1.11 ~^{+0.04}_{-0.05}$ $R_{\odot}$, respectively. This system is a useful test of evolutionary models of young low-mass stars, as it lies in a region of parameter space where observational constraints are scarce; comparison with these models indicates an apparent age of $\sim$3.5-6 Myr. The systemic velocity is within $1σ$ of the cluster value which, along with the presence of lithium absorption, strongly indicates cluster membership. The CoRoT light curve also contains large-amplitude, rapidly evolving out-of-eclipse variations, which are difficult to explain using starspots alone. The system's spectral energy distribution reveals a mid-infrared excess, which we model as thermal emission from a small amount of dust located in the inner cavity of a circumbinary disk. In turn, this opens up the possibility that some of the out-of-eclipse variability could be due to occultations of the central stars by material located at the inner edge or in the central cavity of the circumbinary disk.

preprint2013arXiv

Measuring the rotation period distribution of field M-dwarfs with Kepler

We have analysed 10 months of public data from the Kepler space mission to measure rotation periods of main-sequence stars with masses between 0.3 and 0.55 M_sun. To derive the rotational period we introduce the autocorrelation function and show that it is robust against phase and amplitude modulation and residual instrumental systematics. Of the 2483 stars examined, we detected rotation periods in 1570 (63.2%), representing an increase of a factor ~ 30 in the number of rotation period determination for field M-dwarfs. The periods range from 0.37-69.7 days, with amplitudes ranging from 1.0-140.8 mmags. The rotation period distribution is clearly bimodal, with peaks at ~ 19 and ~ 33 days, hinting at two distinct waves of star formation, a hypothesis that is supported by the fact that slower rotators tend to have larger proper motions. The two peaks of the rotation period distribution form two distinct sequences in period-temperature space, with the period decreasing with increasing temperature, reminiscent of the Vaughan-Preston gap. The period-mass distribution of our sample shows no evidence of a transition at the fully convective boundary. On the other hand, the slope of the upper envelope of the period-mass relation changes sign around 0.55 M_sun, below which period rises with decreasing mass.

preprint2013arXiv

Radial Velocity Variations of Photometrically Quiet, Chromospherically Inactive Kepler Stars: A Link Between RV Jitter and Photometric Flicker

We compare stellar photometric variability, as measured from Kepler light curves by Basri et al. (2011), with measurements of radial velocity (RV) root-mean-square (RMS) variations of all California Planet Search overlap stars. We newly derive rotation periods from the Kepler light curves for all of the stars in our study sample. The RV variations reported herein range from less than 4 m/s to 135 m/s, yet the stars all have amplitudes of photometric variability less than 3 mmag, reflecting the preference of the RV program for chromospherically "quiet" stars. Despite the small size of our sample, we find with high statistical significance that the RV RMS manifests strongly in the Fourier power spectrum of the light curve: stars that are noisier in RV have a greater number of frequency components in the light curve. We also find that spot models of the observed light curves systematically underpredict the observed RV variations by factors of ~2--1000, likely because the low level photometric variations in our sample are driven by processes not included in simple spot models. The stars best fit by these models tend to have simpler light curves, dominated by a single relatively high amplitude component of variability. Finally, we demonstrate that the RV RMS behavior of our sample can be explained in the context of the photometric variability evolutionary diagram introduced by Bastien et al. (2013). We use this diagram to derive the surface gravities of the stars in our sample, revealing many of them to have moved off the main-sequence. More generally, we find that the stars with the largest RV RMS are those that have evolved onto the "flicker floor" sequence in that diagram, characterized by relatively low amplitude but highly complex photometric variations which grow as the stars evolve to become subgiants.

preprint2013arXiv

Stellar Rotation Periods of the Kepler Objects of Interest: A Dearth of Close-in Planets around Fast Rotators

We present a large sample of stellar rotation periods for Kepler Objects of Interest (KOIs), based on three years of public Kepler data. These were measured by detecting periodic photometric modulation caused by star spots, using an algorithm based on the autocorrelation function (ACF) of the light curve, developed recently by McQuillan, Aigrain & Mazeh (2013). Of the 1919 main-sequence exoplanet hosts analyzed, robust rotation periods were detected for 737. Comparing the detected stellar periods to the orbital periods of the innermost planet in each system reveals a notable lack of close-in planets around rapid rotators. It appears that only slowly spinning stars, with rotation periods longer than 5-10 days, host planets on orbits shorter than 3 days, although the mechanism(s) that lead(s) to this is not clear.

preprint2013arXiv

The Deep Blue Color of HD189733b: Albedo Measurements with Hubble Space Telescope/Space Telescope Imaging Spectrograph at Visible Wavelengths

We present a secondary eclipse observation for the hot Jupiter HD189733b across the wavelength range 290-570nm made using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. We measure geometric albedos of Ag = 0.40 \pm 0.12 across 290-450nm and Ag < 0.12 across 450-570nm at 1-sigma confidence. The albedo decrease toward longer wavelengths is also apparent when using six wavelength bins over the same wavelength range. This can be interpreted as evidence for optically thick reflective clouds on the dayside hemisphere with sodium absorption suppressing the scattered light signal beyond ~450nm. Our best-fit albedo values imply that HD189733b would appear a deep blue color at visible wavelengths.

preprint2012arXiv

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

preprint2012arXiv

Detection of Neptune-size planetary candidates with CoRoT data. Comparison with the planet occurrence rate derived from Kepler

[Abridged] Context. The CoRoT space mission has been searching for transiting planets since the end of December 2006. Aims. We aim to investigate the capability of CoRoT to detect small-size transiting planets in short-period orbits, and to compare the number of CoRoT planets with 2 \leq R_p \leq 4 Rearth with the occurrence rate of small-size planets provided by the distribution of Kepler planetary candidates (Howard et al. 2012). Methods. We performed a test that simulates transits of super-Earths and Neptunes in real CoRoT light curves and searches for them blindly by using the LAM transit detection pipeline. Results. The CoRoT detection rate of planets with radius between 2 and 4 Rearth and orbital period P \leq 20 days is 59% (31%) around stars brighter than r'=14.0 (15.5). By properly taking the CoRoT detection rate for Neptune-size planets and the transit probability into account, we found that according to the Kepler planet occurrence rate, CoRoT should have discovered 12 \pm 2 Neptunes orbiting G and K dwarfs with P \leq 17 days in six observational runs. This estimate must be compared with the validated Neptune CoRoT-24b and five CoRoT planetary candidates in the considered range of planetary radii. We thus found a disagreement with expectations from Kepler at 3 σor 5 σ, assuming a blend fraction of 0% (six Neptunes) and 100% (one Neptune) for these candidates. Conclusions. This underabundance of CoRoT Neptunes with respect to Kepler may be due to several reasons. Regardless of the origin of the disagreement, which needs to be investigated in more detail, the noticeable deficiency of CoRoT Neptunes at short orbital periods seems to indirectly support the general trend found in Kepler data, i.e. that the frequency of small-size planets increases with increasing orbital periods and decreasing planet radii.

preprint2012arXiv

On the potential of the EChO mission to characterise gas giant atmospheres

Space telescopes such as EChO (Exoplanet Characterisation Observatory) and JWST (James Webb Space Telescope) will be important for the future study of extrasolar planet atmospheres. Both of these missions are capable of performing high sensitivity spectroscopic measurements at moderate resolutions in the visible and infrared, which will allow the characterisation of atmospheric properties using primary and secondary transit spectroscopy. We use the NEMESIS radiative transfer and retrieval tool (Irwin et al. 2008, Lee et al. 2012) to explore the potential of the proposed EChO mission to solve the retrieval problem for a range of H2-He planets orbiting different stars. We find that EChO should be capable of retrieving temperature structure to ~200 K precision and detecting H2O, CO2 and CH4 from a single eclipse measurement for a hot Jupiter orbiting a Sun-like star and a hot Neptune orbiting an M star, also providing upper limits on CO and NH3. We provide a table of retrieval precisions for these quantities in each test case. We expect around 30 Jupiter-sized planets to be observable by EChO; hot Neptunes orbiting M dwarfs are rarer, but we anticipate observations of at least one similar planet.

preprint2012arXiv

Statistics of Stellar Variability from Kepler - I: Revisiting Quarter 1 with an Astrophysically Robust Systematics Correction

We investigate the variability properties of main sequence stars in the first month of Kepler data, using a new astrophysically robust systematics correction, and find that 60% of stars are more variable then the active Sun. We define low and high variability samples, with a cut corresponding to twice the variability level of the active Sun, and compare the properties of the stars belonging to each sample. We show tentative evidence that the more active stars have lower proper motions and may be located closer to the galactic plane. We also investigate the frequency content of the variability, finding clear evidence for periodic or quasi-periodic behaviour in 16% of stars, and showing that there exist significant differences in the nature of variability between spectral types. Of the periodic objects, most A and F stars have short periods (< 2 days) and highly sinusoidal variability, suggestive of pulsations, whilst G, K and M stars tend to have longer periods (> 5 days, with a trend towards longer periods at later spectral types) and show a mixture of periodic and stochastic variability, indicative of activity. Finally, we use auto-regressive models to characterise the stochastic component of the variability, and show that its typical amplitude and time-scale both increase towards later spectral types, which we interpret as a corresponding increase in the characteristic size and life-time of active regions.

preprint2011arXiv

Spitzer Infrared Observations and Independent Validation of the Transiting Super-Earth CoRoT-7b

The detection and characterization of the first transiting super-Earth, CoRoT-7 b, has required an unprecedented effort in terms of telescope time and analysis. Although the star does display a radial velocity signal at the period of the planet, this has been difficult to disentangle from the intrinsic stellar variability, and pinning down the velocity amplitude has been very challenging. As a result, the precise value of the mass of the planet - and even the extent to which it can be considered to be confirmed - have been debated in the recent literature, with six mass measurements published so far based on the same spectroscopic observations, ranging from about 2 to 8 Earth masses. Here we report on an independent validation of the planet discovery, using one of the fundamental properties of a transit signal: its achromaticity. We observed four transits of CoRoT-7 b with Spitzer, in order to determine whether the depth of the transit signal in the near-infrared is consistent with that observed in the CoRoT bandpass, as expected for a planet. We detected the transit and found an average depth of 0.426 {\pm} 0.115 mmag at 4.5 μm, which is in good agreement with the depth of 0.350 {\pm} 0.011 mmag found by CoRoT. These observations place important constraints on the kinds of astrophysical false positives that could mimic the signal. Combining this with additional constraints reported earlier, we performed an exhaustive exploration of possible blends scenarios for CoRoT-7 b using the BLENDER technique. We are able to rule out the vast majority of false positives, and the remaining ones are found to be much less likely than a true transiting planet. We thus validate CoRoT-7 b as a bona-fide planet with a very high degree of confidence, independently of any radial-velocity information. Our Spitzer observations have additionally allowed us to significantly improve the ephemeris of the planet.

preprint2010arXiv

A new look at NICMOS transmission spectroscopy of HD189733, GJ-436 and XO-1: no conclusive evidence for molecular features

We present a re-analysis of archival HST/NICMOS transmission spectroscopy of three exoplanet systems; HD 189733, GJ-436 and XO-1. Detections of several molecules, including H20, CH4 and CO2, have been claimed for HD 189733 and XO-1, but similarly sized features are attributed to systematic noise for GJ-436. The data consist of time-series grism spectra covering a planetary transit. After extracting light curves in independent wavelength channels, we use a linear decorrelation technique account for instrumental systematics (which is becoming standard in the field), and measure the planet-to-star radius ratio as a function of wavelength. For HD 189733, the uncertainties in the transmission spectrum are significantly larger than those previously reported. We also find the transmission spectrum is considerably altered when using different out-of-transit orbits to remove the systematics, when some parameters are left out of the decorrelation procedure, or when we perform the decorrelation with quadratic functions rather than linear functions. Given that there is no physical reason to believe the baseline flux should be modelled as a linear function of any particular set of parameters, we interpret this as evidence that the linear decorrelation technique is not a robust method to remove systematic effects from the light curves for each wavelength channel. For XO-1, the parameters measured to decorrelate the light curves would require extrapolation to the in-transit orbit to remove the systematics, and we cannot reproduce the previously reported results. We conclude that the resulting NICMOS transmission spectra are too dependent on the method used to remove systematics to be considered robust detections of molecular species in planetary atmospheres, although the presence of these molecules is not ruled out.

preprint2010arXiv

A new look at NICMOS transmission spectroscopy: no conclusive evidence for molecular features

We present a re-analysis of archival HST/NICMOS transmission spectroscopy of the exoplanet system, HD 189733, from which detections of several molecules have been claimed. As expected, we can replicate the transmission spectrum previously published when we use an identical model for the systematic effects, although the uncertainties are larger as we use a residual permutation algorithm in an effort to account for instrumental systematics. We also find that the transmission spectrum is considerably altered when slightly changing the instrument model, and conclude that the NICMOS transmission spectrum is too dependent on the method used to remove systematics to be considered a robust detection of molecular species, given that there is no physical reason to believe that the baseline flux should be modelled as a linear function of any chosen set of parameters.

preprint2010arXiv

Re-assessing the radial-velocity evidence for planets around CoRoT-7

CoRoT-7 is an 11th magnitude K-star whose light curve shows transits with depth of 0.3 mmag and a period of 0.854 d, superimposed on variability at the 1% level, due to the modulation of evolving active regions with the star's 23 d rotation period. In this paper, we revisit the published HARPS radial velocity measurements of the object, which were previously used to estimate the companion mass, but have been the subject of ongoing debate. We build a realistic model of the star's activity during the HARPS observations, by fitting simultaneously the line width and the line bisector, and use it to evaluate the contribution of activity to the RV variations. The data show clear evidence of errors above the level of the formal uncertainties, which are accounted for either by activity, nor by any plausible planet model, and which increase rapidly with decreasing signal-to-noise of the spectra. We cite evidence of similar systematics in mid-SNR spectra of other targets obtained with HARPS and other high-precision RV spectrographs, and discuss possible sources. Allowing for these, we re-evaluate the semi-amplitude of the CoRoT-7b signal, finding Kb=1.6 +-1.3 m/s, a tentative detection with a much reduced significance (1.2-sigma) compared to previous estimates. We also argue that the combined presence of activity and additional errors preclude a meaningful search for additional low-mass companions, despite previous claims to the contrary. Our analysis points to a lower density for CoRoT-7b, the 1-sigma mass range spanning 1-4 MEarth, allowing for a wide range of bulk compositions. In particular, an ice-rich composition is compatible with the RV constraints. This study highlights the importance of a realistic treatment of both activity and uncertainties, particularly in the medium signal-to-noise ratio regime, which applies to most small planet candidates from CoRoT and Kepler.

preprint2010arXiv

The SARS algorithm: detrending CoRoT light curves with Sysrem using simultaneous external parameters

Surveys for exoplanetary transits are usually limited not by photon noise but rather by the amount of red noise in their data. In particular, although the CoRoT spacebased survey data are being carefully scrutinized, significant new sources of systematic noises are still being discovered. Recently, a magnitude-dependant systematic effect was discovered in the CoRoT data by Mazeh & Guterman et al. and a phenomenological correction was proposed. Here we tie the observed effect a particular type of effect, and in the process generalize the popular Sysrem algorithm to include external parameters in a simultaneous solution with the unknown effects. We show that a post-processing scheme based on this algorithm performs well and indeed allows for the detection of new transit-like signals that were not previously detected.

preprint2010arXiv

The thermal emission of the exoplanets WASP-1b and WASP-2b

We present a comparative study of the thermal emission of the transiting exoplanets WASP-1b and WASP-2b using the Spitzer Space Telescope. The two planets have very similar masses but suffer different levels of irradiation and are predicted to fall either side of a sharp transition between planets with and without hot stratospheres. WASP-1b is one of the most highly irradiated planets studied to date. We measure planet/star contrast ratios in all four of the IRAC bands for both planets (3.6-8.0um), and our results indicate the presence of a strong temperature inversion in the atmosphere of WASP-1b, particularly apparent at 8um, and no inversion in WASP-2b. In both cases the measured eclipse depths favor models in which incident energy is not redistributed efficiently from the day side to the night side of the planet. We fit the Spitzer light curves simultaneously with the best available radial velocity curves and transit photometry in order to provide updated measurements of system parameters. We do not find significant eccentricity in the orbit of either planet, suggesting that the inflated radius of WASP-1b is unlikely to be the result of tidal heating. Finally, by plotting ratios of secondary eclipse depths at 8um and 4.5um against irradiation for all available planets, we find evidence for a sharp transition in the emission spectra of hot Jupiters at an irradiation level of 2 x 10^9 erg/s/cm^2. We suggest this transition may be due to the presence of TiO in the upper atmospheres of the most strongly irradiated hot Jupiters.

preprint2009arXiv

ASTEP South: An Antarctic Search for Transiting ExoPlanets around the celestial South pole

ASTEP South is the first phase of the ASTEP project (Antarctic Search for Transiting ExoPlanets). The instrument is a fixed 10 cm refractor with a 4kx4k CCD camera in a thermalized box, pointing continuously a 3.88 degree x 3.88 degree field of view centered on the celestial South pole. ASTEP South became fully functional in June 2008 and obtained 1592 hours of data during the 2008 Antarctic winter. The data are of good quality but the analysis has to account for changes in the point spread function due to rapid ground seeing variations and instrumental effects. The pointing direction is stable within 10 arcseconds on a daily timescale and drifts by only 34 arcseconds in 50 days. A truly continuous photometry of bright stars is possible in June (the noon sky background peaks at a magnitude R=15 arcsec-2 on June 22), but becomes challenging in July (the noon sky background magnitude is R=12.5 arcsec?2 on July 20). The weather conditions are estimated from the number of stars detected in the field. For the 2008 winter, the statistics are between 56.3 % and 68.4 % of excellent weather, 17.9 % to 30 % of veiled weather and 13.7 % of bad weather. Using these results in a probabilistic analysis of transit detection, we show that the detection efficiency of transiting exoplanets in one given field is improved at Dome C compared to a temperate site such as La Silla. For example we estimate that a year-long campaign of 10 cm refractor could reach an efficiency of 69 % at Dome C versus 45 % at La Silla for detecting 2-day period giant planets around target stars from magnitude 10 to 15. This shows the high potential of Dome C for photometry and future planet discoveries. [Short abstract]

preprint2008arXiv

The Monitor project: the search for transits in the open cluster NGC 2362

We present the results of a systematic search for transiting planets in a ~5 Myr open cluster, NGC 2362. We observed ~1200 candidate cluster members, of which ~475 are believed to be genuine cluster members, for a total of ~100 hours. We identify 15 light curves with reductions in flux that pass all our detection criteria, and 6 of the candidates have occultation depths compatible with a planetary companion. The variability in these six light curves would require very large planets to reproduce the observed transit depth. If we assume that none of our candidates are in fact planets then we can place upper limits on the fraction of stars with hot Jupiters (HJs) in NGC 2362. We obtain 99% confidence upper limits of 0.22 and 0.70 on the fraction of stars with HJs (f_p) for 1-3 and 3-10 day orbits, respectively, assuming all HJs have a planetary radius of 1.5R_Jup. These upper limits represent observational constraints on the number of stars with HJs at an age <~10 Myr, when the vast majority of stars are thought to have lost their protoplanetary discs. Finally, we extend our results to the entire Monitor Project, a survey searching young, open clusters for planetary transits, and find that the survey as currently designed should be capable of placing upper limits on f_p near the observed values of f_p in the solar neighbourhood.