Showing papers by "Suzanne Aigrain published in 2018"

Inferring probabilistic stellar rotation periods using Gaussian processes

Abstract: 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.

K2 photometry and HERMES spectroscopy of the blue supergiant ρ Leo: rotational wind modulation and low-frequency waves★

Abstract: We present an 80-d long uninterrupted high-cadence K2 light curve of the B1Iab supergiant rho Leo (HD 91316), deduced with the method of halo photometry This light curve reveals a dominant frequency of $f_{\rmrot}=00373$d$^{-1}$ and its harmonics This dominant frequency corresponds with a rotation period of 268d and is subject to amplitude and phase modulation The K2 photometry additionally reveals multiperiodic low-frequency variability ($<15 $d$^{-1}$) and is in full agreement with low-cadence high-resolution spectroscopy assembled during 1800 days The spectroscopy reveals rotational modulation by a dynamic aspherical wind with an amplitude of about 20km s$^{-1}$ in the H$\alpha$ line, as well as photospheric velocity variations of a few km s$^{-1}$ at frequencies in the range 02 to 06 d$^{-1}$ in the SiIII 4567\AA\ line Given the large macroturbulence needed to explain the spectral line broadening of the star, we interpret the detected photospheric velocity as due to travelling super-inertial low-degree large-scale gravity waves with dominant tangential amplitudes and discuss why $\rho$~Leo is an excellent target to study how the observed photospheric variability propagates into the wind

A universal spin mass relation for brown dwarfs and planets

Abstract: While brown dwarfs show similarities to stars early in their lives, their spin evolutions are much more akin to those of planets. We have used light curves 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 Moand has been characterized 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 versus 1.6 days), a signature of rotational braking by the disk, albeit with small numbers. With an overall median period of 1.9 days, 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 substellar 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.

Planets, candidates, and binaries from the CoRoT/Exoplanet programme. The CoRoT transit catalogue

Abstract: The CoRoT space mission observed 163 665 stars over 26 stellar fields in the faint star channel. The exoplanet teams detected a total of 4123 transit-like features in the 177 454 light curves. We present the complete re-analysis of all these detections carried out with the same softwares so that to ensure their homogeneous analysis. Although the vetting process involves some human evaluation, it also involves a simple binary flag system over basic tests: Detection significance, presence of a secondary, difference between odd and even depths, colour dependence, V-shape transit, and duration of the transit. We also gathered the information from the large accompanying ground-based programme carried out on the planet candidates and checked how useful the flag system could have been at the vetting stage of the candidates. From the initial list of transit-like features, we identified and separated 824 false alarms of various kind, 2269 eclipsing binaries among which 616 are contact binaries and 1653 are detached ones, 37 planets and brown dwarfs, and 557 planet candidates. We provide the catalogue of all these transit-like features, including false alarms. For the planet candidates, the catalogue gives not only their transit parameters but also the products of their light curve modelling: Reduced radius, reduced semi-major axis, and impact parameter, together with a summary of the outcome of follow-up observations when carried out and their current status. For the detached eclipsing binaries, the catalogue provides, in addition to their transit parameters, a simple visual classification. Among the planet candidates whose nature remains unresolved, we estimate that eight (within an error of three) planets are still to be identified. After correcting for geometric and sensitivity biases, we derived planet and brown dwarf occurrences and confirm disagreements with Kepler estimates, as previously reported by other authors from the analysis of the first runs: Small-size planets with orbital period less than ten days are underabundant by a factor of three in the CoRoT fields whereas giant planets are overabundant by a factor of two. These preliminary results would however deserve further investigations using the recently released CoRoT light curves that are corrected of the various instrumental effects and a homogeneous analysis of the stellar populations observed by the two missions.

Planets, candidates, and binaries from the CoRoT/Exoplanet programme: the CoRoT transit catalogue

Abstract: We provide the catalogue of all transit-like features, including false alarms, detected by the CoRoT exoplanet teams in the 177 454 light curves of the mission. All these detections have been re-analysed with the same softwares so that to ensure their homogeneous analysis. Although the vetting process involves some human evaluation, it also involves a simple binary flag system over basic tests: detection significance, presence of a secondary, difference between odd and even depths, colour dependence, V-shape transit, and duration of the transit. We also gathered the information from the large accompanying ground-based programme carried out on the planet candidates and checked how useful the flag system could have been at the vetting stage of the candidates. In total, we identified and separated 824 false alarms of various kind, 2269 eclipsing binaries among which 616 are contact binaries and 1653 are detached ones, 37 planets and brown dwarfs, and 557 planet candidates. For the planet candidates, the catalogue gives not only their transit parameters but also the products of their light curve modelling, together with a summary of the outcome of follow-up observations when carried out and their current status. Among the planet candidates whose nature remains unresolved, we estimate that 8 +/- 3 planets are still to be identified. We derived planet and brown dwarf occurrences and confirm disagreements with Kepler estimates: small-size planets with orbital period less than ten days are underabundant by a factor of three in the CoRoT fields whereas giant planets are overabundant by a factor of two. These preliminary results would however deserve further investigations using the recently released CoRoT light curves that are corrected of the various instrumental effects and a homogeneous analysis of the stellar populations observed by the two missions.

The K2 M67 Study: Establishing the Limits of Stellar Rotation Period Measurements in M67 with K2 Campaign 5 Data

Abstract: The open cluster M67 offers a 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, 75 day observation window, crowded field, and typically low-amplitude signals mean that determining accurate rotation periods on the order of 25-30 days is inherently difficult. Lingering, non-astrophysical signals with power at ≥25 days 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 35 days 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 a 15% recovery rate for the solar case (i.e., 25 day 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).

The GTC exoplanet transit spectroscopy survey. VIII. Flat transmission spectrum for the warm gas giant WASP-80b

Abstract: Aims. We set out to study the atmosphere of WASP-80b, a warm inflated gas giant with an equilibrium temperature of ~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.Methods. We observed two spectroscopic time series of WASP-80b transits with the OSIRIS spectrograph installed in the Gran Telescopio Canarias (GTC), 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 jointly model three previously published broadband datasets consisting of 27 light curves, prior to a transmission spectroscopy analysis in order to obtain improved estimates of the planet’s orbital parameters, average radius ratio, and stellar density. The parameter posteriors from the broadband analysis are used to set informative priors on the transmission spectroscopy analysis. The final transmission spectroscopy analyses are carried out jointly for the two nights using a divide-by-white approach to remove the common-mode systematics, and Gaussian processes to model the residual wavelength-dependent systematics.Results. 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 model. The flat transmission spectrum favours an atmosphere model with high-altitude clouds over cloud-free models with stellar or sub-stellar metallicities.Conclusions. Our results disagree with the recently published discovery of strong K I absorption in WASP-80b’s atmosphere based on ground-based transmission spectroscopy with FORS2 at VLT.

The K2 M67 Study: A Curiously Young Star in an Eclipsing Binary in an Old Open Cluster

Abstract: We present an analysis of a slightly eccentric (e = 0.05), partially eclipsing, long-period (P = 69.73 days) 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 ± 0.015 and 0.748 ± 0.005 Mo) and radii (1.071 ± 0.008 ± 0.003 and 0.713 ± 0.019 ± 0.026 Ro, 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 nondetection 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 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 1 Gyr 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.

A universal spin-mass relation for brown dwarfs and planets

Abstract: 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.

The K2 M67 Study: Establishing the Limits of Stellar Rotation Period Measurements in M67 with K2 Campaign 5 Data

Abstract: 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.

VizieR Online Data Catalog: CoRoT transit catalogue (Deleuil+, 2018)

Abstract: We provide the full catalogue of all the transit-like events identified in the CoRoT light curves during the mission lifetime by the CoRoT collaboration. It includes planet candidates that were identified after the release of each run, their status once follow-up observations were completed, and also binaries. The later were separated in two groups: detached eclipsing binaries and contact binaries. All these events, whose detection was performed by different teams in a first stage, have been re-analysed with a unified procedure and using the same softwares so that to derive, in a homogeneous and consistent way, the complete set of transit parameters, using fixed limb-darkening coefficients because of the uncertainties on the spectral classification. In addition to this modelling, we also carried out a vetting process of both candidates and binaries. It is based on a simple binary flag system over basic tests: detection significance, presence of a secondary, odd or even depth differences, colour dependence, Vshape transit, and duration of the transit. Among the 4123 transitlike signals that were analysed, it allowed us to identify periodic signals that are instrumental false positives and false detections. In total, the presented catalogue comes up with 594 events initially flagged as planet candidates which include published planets and brown dwarfs, 1653 EBs, 616 contact binaries, and 824 false alarms. (6 data files).

The K2 M67 Study: A Curiously Young Star in an Eclipsing Binary in an Old Open Cluster

Abstract: 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.