Showing 1–50 of 201 results for author: Aigrain, S

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

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

Submitted 8 June, 2024; originally announced June 2024.

Investigating stellar activity through eight years of Sun-as-a-star observations

Authors: Baptiste Klein, Suzanne Aigrain, Michael Cretignier, Khaled Al Moulla, Xavier Dumusque, Oscar Barragán, Haochuan Yu, Annelies Mortier, Federica Rescigno, Andrew Collier Cameron, Mercedes López-Morales, Nadège Meunier, Alessandro Sozzetti, Niamh K. O'Sullivan

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

Submitted 22 May, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: Accepted by MNRAS on the 2024 May 20. 17 pages, 14 figures (plus 8 pages of Appendix, 9 Appendix figures)

Planet Hunters TESS V: a planetary system around a binary star, including a mini-Neptune in the habitable zone

Authors: Nora L. Eisner, Samuel K. Grunblatt, Oscar Barragán, Thea H. Faridani, Chris Lintott, Suzanne Aigrain, Cole Johnston, Ian R. Mason, Keivan G. Stassun, Megan Bedell, Andrew W. Boyle, David R. Ciardi, Catherine A. Clark, Guillaume Hebrard, David W. Hogg, Steve B. Howell, Baptiste Klein, Joe Llama, Joshua N. Winn, Lily L. Zhao, Joseph M. Akana Murphy, Corey Beard, Casey L. Brinkman, Ashley Chontos, Pia Cortes-Zuleta , et al. (39 additional authors not shown)

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

Submitted 29 April, 2024; originally announced April 2024.

Comments: 24 pages, 16 figures, 4 tables

Journal ref: Published in AJ, 2024

TOI-837b is a Young Saturn-sized Exoplanet with a Massive 70 $M_{\oplus}$ Core

Authors: Oscar Barragán, Haochuan Yu, Alix Violet Freckelton, Annabella Meech, Michael Cretignier, Annelies Mortier, Suzanne Aigrain, Baptiste Klein, Niamh K. O'Sullivan, Edward Gillen, Louise Dyregaard Nielsen, Manuel Mallorquín, Norbert Zicher

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

Submitted 24 June, 2024; v1 submitted 21 April, 2024; originally announced April 2024.

Comments: MNRAS published version, MNRAS, 531, 4275

Modelling stochastic and quasi-periodic behaviour in stellar time-series: Gaussian process regression versus power-spectrum fitting

Authors: Niamh K. O'Sullivan, Suzanne Aigrain

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

Submitted 17 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS

REPUBLIC: A variability-preserving systematic-correction algorithm for PLATO's multi-camera light curves

Authors: Oscar Barragán, Suzanne Aigrain, James McCormac

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

Submitted 9 April, 2024; originally announced April 2024.

Comments: 12 pages, 6 figures, accepted for publication in RAS Techniques and Instruments

nuance: Efficient detection of planets transiting active stars

Authors: Lionel Garcia, Daniel Foreman-Mackey, Catriona A. Murray, Suzanne Aigrain, Dax L. Feliz, Francisco J. Pozuelos

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

Submitted 17 May, 2024; v1 submitted 9 February, 2024; originally announced February 2024.

Modelling stellar variability in archival HARPS data: I -- Rotation and activity properties with multi-dimensional Gaussian Processes

Authors: Haochuan Yu, Suzanne Aigrain, Baptiste Klein, Oscar Barragán, Annelies Mortier, Niamh K. O'Sullivan, Michael Cretignier

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

Submitted 10 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS

Stellar surface information from the Ca II H&K lines I. Intensity profiles of the solar activity components

Authors: M. Cretignier, A. G. M. Pietrow, S. Aigrain

Abstract: The detection of Earth-like planets with the radial-velocity method is currently limited by the presence of stellar activity signatures. On rotational timescales, spots and plages (or faculae) are known to introduce different RV signals, but their corrections require better activity proxies. The best-known chromospheric activity proxies in the visible are the Ca II H & K lines, but the physical qu… ▽ More The detection of Earth-like planets with the radial-velocity method is currently limited by the presence of stellar activity signatures. On rotational timescales, spots and plages (or faculae) are known to introduce different RV signals, but their corrections require better activity proxies. The best-known chromospheric activity proxies in the visible are the Ca II H & K lines, but the physical quantities measured by their profiles need to be clarified. We first investigate resolved images of the Sun in order to better understand the spectrum of plages, spots, and the network using the Meudon spectroheliogram. We show that distinct line profiles are produced by plages, spots, and by the network component and we also derived the center-to-limb variations of the three profiles. Some care is required to disentangle their contributions due to their similarities. By combining disk-integrated spectra from the ISS high-resolution spectrograph with SDO direct images of the Sun, we managed to extract a high-resolution emission spectrum of the different components, which tend to confirm the spectra extracted from the Meudon spectroheliogram datacubes. Similar results were obtained with the HARPS-N Sun-as-a-star spectra. We concluded using a three-component model that the temporal variation of the popular S-index contains, on average for the 24th solar cycle: 70 +/- 12% of plage, 26 +/- 12% of network and 4 +/- 4% of spots. This preliminary investigation suggests that a detailed study of the Ca II H & K profiles may provide rich information about the filling factor and distribution of different types of active regions. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 20 pages, 17 figures

HD152843 b & c: the masses and orbital periods of a sub-Neptune and a super-puff Neptune

Authors: B. A. Nicholson, S. Aigrain, N. L. Eisner, M. Cretignier, O. Barragán, L. Kaye, J. Taylor, J. Owen, A. Mortier, L. Affer, W. Boschin, A. Collier Cameron, M. Damasso, L. Di Fabrizio, V. DiTomasso, X. Dumusque, A. Gehdina, A. Harutyunyan, D. W. Latham, M. Lopez-Morales, V. Lorenzi, A. F. Martínez Fiorenzano, E. Molinari, M. Pedani, M. Pinamonti , et al. (2 additional authors not shown)

Abstract: We present the characterisation of the two transiting planets around HD 152843 (TOI 2319, TIC 349488688) using an intensive campaign of HARPS-N radial velocities, and two sectors of TESS data. These data reveal a unique and fascinating system: HD 152843 b and c have near equal masses of around 9 M$_\oplus$ but differing radii of $3.05 \pm 0.11$ R$_\oplus$ and $5.94^{+0.18}_{-0.16}$ R$_\oplus$ , re… ▽ More We present the characterisation of the two transiting planets around HD 152843 (TOI 2319, TIC 349488688) using an intensive campaign of HARPS-N radial velocities, and two sectors of TESS data. These data reveal a unique and fascinating system: HD 152843 b and c have near equal masses of around 9 M$_\oplus$ but differing radii of $3.05 \pm 0.11$ R$_\oplus$ and $5.94^{+0.18}_{-0.16}$ R$_\oplus$ , respectively, and orbital periods of $11.62071^{+9.6e-05}_{-0.000106}$ days and $19.502104^{+7.4e-05}_{-8.5e-05}$ days. This indicates that HD 152843 c is in the lowest fifth percentile in density of the known exoplanet population, and has the longest orbital period among these low density planets. Further, HD 152843 c's radius places it in the Saturn valley, the observed lack of planets larger than Neptune, but smaller than Saturn. The orbital periods of these planets indicate they are near a 5:3 mean motion resonance, indicating the possibility of transit timing variations, and hints at the possibility of interaction with a third planet at some point in the evolution of this system. Further, the brightness of the host star and the low density of HD 152843 c make it a key target for atmospheric characterisation. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Submitted for review to the Monthly Notices of the Royal Astronomical Society. 13 pages, 16 figures

PlatoSim: An end-to-end PLATO camera simulator for modelling high-precision space-based photometry

Authors: N. Jannsen, J. De Ridder, D. Seynaeve, S. Regibo, R. Huygen, P. Royer, C. Paproth, D. Grießbach, R. Samadi, D. R. Reese, M. Pertenais, E. Grolleau, R. Heller, S. M. Niemi, J. Cabrera, A. Börner, S. Aigrain, J. McCormac, P. Verhoeve, P. Astier, N. Kutrowski, B. Vandenbussche, A. Tkachenko, C. Aerts

Abstract: PLAnetary Transits and Oscillations of stars (PLATO) is the ESA M3 space mission dedicated to detect and characterise transiting exoplanets including information from the asteroseismic properties of their stellar hosts. The uninterrupted and high-precision photometry provided by space-borne instruments such as PLATO require long preparatory phases. An exhaustive list of tests are paramount to desi… ▽ More PLAnetary Transits and Oscillations of stars (PLATO) is the ESA M3 space mission dedicated to detect and characterise transiting exoplanets including information from the asteroseismic properties of their stellar hosts. The uninterrupted and high-precision photometry provided by space-borne instruments such as PLATO require long preparatory phases. An exhaustive list of tests are paramount to design a mission that meets the performance requirements, and as such, simulations are an indispensable tool in the mission preparation. To accommodate PLATO's need of versatile simulations prior to mission launch - that at the same time describe accurately the innovative but complex multi-telescope design - we here present the end-to-end PLATO simulator specifically developed for the purpose, namely PlatoSim. We show step-by-step the algorithms embedded into the software architecture of PlatoSim that allow the user to simulate photometric time series of CCD images and light curves in accordance to the expected observations of PLATO. In the context of the PLATO payload, a general formalism of modelling, end-to-end, incoming photons from the sky to the final measurement in digital units is discussed. We show the strong predictive power of PlatoSim through its diverse applicability and contribution to numerous working groups within the PLATO Mission Consortium. This involves the on-going mechanical integration and alignment, performance studies of the payload, the pipeline development and assessments of the scientific goals. PlatoSim is a state-of-the-art simulator that is able to produce the expected photometric observations of PLATO to a high level of accuracy. We demonstrate that PlatoSim is a key software tool for the PLATO mission in the preparatory phases until mission launch and prospectively beyond. △ Less

Submitted 8 May, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 30 pages, 22 figures

Journal ref: A&A 681, A18 (2024)

YARARA V2: Reaching sub m/s precision over a decade using PCA on line-by-line RVs

Authors: M. Cretignier, X. Dumusque, S. Aigrain, F. Pepe

Abstract: Context. The detection of Earth-like planets with the radial-velocity (RV) method is extremely challenging today due to the presence of non-Doppler signatures such as stellar activity and instrumental signals that mimic and hide the signals of exoplanets. In a previous paper, we presented the YARARA pipeline, which implements corrections for telluric absorption, stellar activity and instrumental s… ▽ More Context. The detection of Earth-like planets with the radial-velocity (RV) method is extremely challenging today due to the presence of non-Doppler signatures such as stellar activity and instrumental signals that mimic and hide the signals of exoplanets. In a previous paper, we presented the YARARA pipeline, which implements corrections for telluric absorption, stellar activity and instrumental systematics at the spectral level, then extracts line-by-line (LBL) RVs with significantly better precision than standard pipelines. Aims. In this paper, we demonstrate that further gains in RVs precision can be achieved by performing Principal Component Analysis (PCA) decomposition on the LBL RVs. Methods. The mean-insensitive nature of PCA means that it is unaffected by true Doppler shifts, and thus can be used to isolate and correct nuisance signals other than planets. Results. We analysed the data of 20 intensively observed HARPS targets by applying our PCA approach on the LBL RVs obtained by YARARA. The first principal components show similarities across most of the stars and correspond to newly identified instrumental systematics, which we can now correct for. For several targets, this results in an unprecedented RV root-mean-square of around 90 cm/s over the full lifetime of HARPS. We use the corrected RVs to confirm a previously published 120-day signal around 61Vir, and to detect a Super-Earth candidate (K = 60 +/- 6 cm/s, m sin i = 6.6 +/- 0.7 Earth mass) around the G6V star HD20794, which spends part of its 600-day orbit within the habitable zone of the host star. Conclusions. This study highlights the potential of LBL PCA to identify and correct hitherto unknown, long-term instrumental effects and thereby extend the sensitivity of existing and future instruments towards the Earth analogue regime. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 22 pages, 19 figures

Journal ref: A&A 678, A2 (2023)

Transit Timing Variations in the three-planet system: TOI-270

Authors: Laurel Kaye, Shreyas Vissapragada, Maximilian N. Gunther, Suzanne Aigrain, Thomas Mikal-Evans, Eric L. N. Jensen, Hannu Parviainen, Francisco J. Pozuelos, Lyu Abe, Jack S. Acton, Abdelkrim Agabi, Douglas R. Alves, David R. Anderson, David J. Armstrong, Khalid Barkaoui, Oscar Barragan, Bjorn Benneke, Patricia T. Bo yd, Rafael Brahm, Ivan Bruni, Edward M. Bryant, Matthew R. Burleigh, Sarah L. Casewell, David Ciardi, Ryan Cloutier , et al. (47 additional authors not shown)

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

Submitted 21 August, 2023; originally announced August 2023.

Comments: 22 pages, 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 510, Issue 4, pp.5464-5485 (2022)

Revisiting K2-233 spectroscopic time-series with multidimensional Gaussian Processes

Authors: Oscar Barragán, Edward Gillen, Suzanne Aigrain, Annabella Meech, Baptiste Klein, Louise Dyregaard Nielsen, Haochuan Yu, Niamh K. O'Sullivan, Belinda A. Nicholson, Jorge Lillo-Box

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

Submitted 13 April, 2023; originally announced April 2023.

Comments: Accepted for publication in MNRAS

Discovering planets with PLATO: Comparison of algorithms for stellar activity filtering

Authors: G. Canocchi, L. Malavolta, I. Pagano, O. Barragán, G. Piotto, S. Aigrain, S. Desidera, S. Grziwa, J. Cabrera, H. Rauer

Abstract: Context. To date, stellar activity is one of the main limitations in detecting small exoplanets via transit photometry. Since this activity is enhanced in young stars, traditional filtering algorithms may severely under-perform in detecting such exoplanets. Aims.This paper aims to compare the relative performances of four algorithms developed by independent research groups specifically for the f… ▽ More Context. To date, stellar activity is one of the main limitations in detecting small exoplanets via transit photometry. Since this activity is enhanced in young stars, traditional filtering algorithms may severely under-perform in detecting such exoplanets. Aims.This paper aims to compare the relative performances of four algorithms developed by independent research groups specifically for the filtering of activity in the light curves (LCs) of young active stars, prior to the search for planetary transit signals: Notch and LOCoR(N&L), Young Stars Detrending(YSD), K2 Systematics Correction(K2SC) and VARLET. We include in the comparison also the two best-performing algorithms implemented in Wotan, namely the Tukey's biweight and the Huber Spline. Methods. We performed a series of injection-retrieval tests of planetary transits of different types, from Jupiter down to Earth-sized planets, moving both on circular and eccentric orbits. The tests were carried out over 100 simulated LCs of both quiet and active solar-like stars that will be observed by the ESA space telescope PLATO. Results. We found that N&L is the best choice in many cases, since it misses the lowest number of transits. However, it under-performs if the planetary orbital period closely matches the stellar rotation period, especially in the case of small planets for which the biweight and VARLET algorithms work better. For LCs with a large number of data, the combined results of YSD and Huber Spline yield the highest recovery percentage. Filtering algorithms allow us to get a very precise estimate of the orbital period and the mid-transit time of the detected planets, while the planet-to-star radius is under-estimated most of the time, especially in the case of grazing transits or eccentric orbits. A refined filtering taking into account the presence of the planet is compulsory for a proper planetary characterization. △ Less

Submitted 28 February, 2023; v1 submitted 4 February, 2023; originally announced February 2023.

Comments: 25 pages, accepted for publication by A&A

Journal ref: A&A 672, A144 (2023)

Gaussian Process regression for astronomical time-series

Authors: Suzanne Aigrain, Daniel Foreman-Mackey

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

Submitted 10 November, 2022; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: Updated in response to community feedback; source code: https://github.com/dfm/araa-gps

Direct discovery of the inner exoplanet in the HD206893 system. Evidence for deuterium burning in a planetary-mass companion

Authors: S. Hinkley, S. Lacour, G. -D. Marleau, A. M. Lagrange, J. J. Wang, J. Kammerer, A. Cumming, M. Nowak, L. Rodet, T. Stolker, W. -O. Balmer, S. Ray, M. Bonnefoy, P. Mollière, C. Lazzoni, G. Kennedy, C. Mordasini, R. Abuter, S. Aigrain, A. Amorim, R. Asensio-Torres, C. Babusiaux, M. Benisty, J. -P. Berger, H. Beust , et al. (89 additional authors not shown)

Abstract: Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoi… ▽ More Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $μ$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au. △ Less

Submitted 3 April, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted to A&A

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

Authors: Belinda A. Nicholson, Suzanne Aigrain

Abstract: 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 interpr… ▽ More 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. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 18 pages, 24 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

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

Authors: Annabella Meech, Suzanne Aigrain, Matteo Brogi, Jayne Birkby

Abstract: 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 di… ▽ More 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. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 14 pages, 7 figures, 3 tables. Accepted for publication in MNRAS

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

Authors: Baptiste Klein, Norbert Zicher, Robert D. Kavanagh, Louise D. Nielsen, Suzanne Aigrain, Aline A. Vidotto, Oscar Barragán, Antoine Strugarek, Belinda Nicholson, Jean-françois Donati, Jérôme Bouvier

Abstract: 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… ▽ More 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. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS (15/03/2022); 19 pages including 6 pages of appendices; 14 figures including 5 in appendices

Report number: MN-22-0224-MJ.R1

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

Authors: Norbert Zicher, Oscar Barragán, Baptiste Klein, Suzanne Aigrain, James E. Owen, Davide Gandolfi, Anne-Marie Lagrange, Luisa Maria Serrano, Laurel Kaye, Louise Dyregaard Nielsen, Vinesh Maguire Rajpaul, Antoine Grandjean, Elisa Goffo, Belinda Nicholson

Abstract: 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 c… ▽ More 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. △ Less

Submitted 3 March, 2022; originally announced March 2022.

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

Authors: Nora L. Eisner, Cole Johnston, Silvia Toonen, Abigail J. Frost, Soetkin Janssens, Chris J. Lintott, Suzanne Aigrain, Hugues Sana, Michael Abdul-Masih, Karla Z. Arellano-Córdova, Paul G. Beck, Emma Bordier, Emily Canon, Ana Escorza, Mattias Fabry, Lars Hermansson, Steve Howell, Grant Miller, Shreeya Sheyte, Safaa Alhassan, Elisabeth M. L. Baeten, Frank Barnet, Stewart. J. Bean, Mikael Bernau, David M. Bundy , et al. (15 additional authors not shown)

Abstract: 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, w… ▽ More 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. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS (11 pages, 8 figures)

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

Authors: Lily L. Zhao, Debra A. Fischer, Eric B. Ford, Alex Wise, Michaël Cretignier, Suzanne Aigrain, Oscar Barragan, Megan Bedell, Lars A. Buchhave, João D. Camacho, Heather M. Cegla, Jessi Cisewski-Kehe, Andrew Collier Cameron, Zoe L. de Beurs, Sally Dodson-Robinson, Xavier Dumusque, João P. Faria, Christian Gilbertson, Charlotte Haley, Justin Harrell, David W. Hogg, Parker Holzer, Ancy Anna John, Baptiste Klein, Marina Lafarga , et al. (18 additional authors not shown)

Abstract: 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-consist… ▽ More 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. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 33 pages (+12 pages of Appendix), 10 figures, 8 tables, accepted for publication in AJ

The young HD 73583 (TOI-560) planetary system: Two 10-M$_\oplus$ mini-Neptunes transiting a 500-Myr-old, bright, and active K dwarf

Authors: O. Barragán, D. J. Armstrong, D. Gandolfi, I. Carleo, A. A. Vidotto, C. Villarreal D'Angelo, A. Oklopčić, H. Isaacson, D. Oddo, K. Collins, M. Fridlund, S. G. Sousa, C. M. Persson, C. Hellier, S. Howell, A. Howard, S. Redfield, N. Eisner, I. Y. Georgieva, D. Dragomir, D. Bayliss, L. D. Nielsen, B. Klein, S. Aigrain, M. Zhang , et al. (82 additional authors not shown)

Abstract: We present the discovery and characterisation of two transiting planets observed by \textit{TESS} in the light curves of the young and bright (V=9.67) star HD73583 (TOI-560). We perform an intensive spectroscopic and photometric space- and ground-based follow-up in order to confirm and characterise the system. We found that HD73583 is a young ($\sim 500$~Myr) active star with a rotational period o… ▽ More We present the discovery and characterisation of two transiting planets observed by \textit{TESS} in the light curves of the young and bright (V=9.67) star HD73583 (TOI-560). We perform an intensive spectroscopic and photometric space- and ground-based follow-up in order to confirm and characterise the system. We found that HD73583 is a young ($\sim 500$~Myr) active star with a rotational period of $12.08 \pm 0.11 $\,d, and a mass and radius of $ 0.73 \pm 0.02 M_\odot$ and $0.65 \pm 0.02 R_\odot$, respectively. HD73583 b ($P_b=6.3980420 _{ - 0.0000062 }^{+0.0000067}$ d) has a mass and radius of $10.2 _{-3.1}^{+3.4} M_\oplus$ and$2.79 \pm 0.10 R_\oplus$, respectively, that gives a density of $2.58 _{-0.81}^{ 0.95} {\rm g\,cm^{-3}}$. HD73583 c ($P_c= 18.87974 _{-0.00074 }^{+0.00086}$) has a mass and radius of $9.7_{-1.7} ^ {+1.8} M_\oplus$ and $2.39_{-0.09}^{+0.10} R_\oplus$, respectively, this translates to a density of $3.88 _{-0.80}^{+0.91} {\rm g\,cm^{-3}}$. Both planets are consistent with worlds made of a solid core surrounded by a volatile envelope. Because of their youth and host star brightness, they both are excellent candidates to perform transmission spectroscopy studies. We expect ongoing atmospheric mass-loss for both planets caused by stellar irradiation. We estimate that the detection of evaporating signatures on H and He would be challenging, but doable with present and future instruments. △ Less

Submitted 7 March, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS

pyaneti II: A multidimensional Gaussian process approach to analysing spectroscopic time-series

Authors: Oscar Barragán, Suzanne Aigrain, Vinesh M. Rajpaul, Norbert Zicher

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

Submitted 22 November, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: This version looks closer to the published one: MNRAS, 509, 866

A HARPS-N mass for the elusive Kepler-37d: a case study in disentangling stellar activity and planetary signals

Authors: V. M. Rajpaul, L. A. Buchhave, G. Lacedelli, K. Rice, A. Mortier, L. Malavolta, S. Aigrain, L. Borsato, A. W. Mayo, D. Charbonneau, M. Damasso, X. Dumusque, A. Ghedina, D. W. Latham, M. López-Morales, A. Magazzù, G. Micela, E. Molinari, F. Pepe, G. Piotto, E. Poretti, S. Rowther, A. Sozzetti, S. Udry, C. A. Watson

Abstract: To date, only 18 exoplanets with radial velocity (RV) semi-amplitudes $<2$ m/s have had their masses directly constrained. The biggest obstacle to RV detection of such exoplanets is variability intrinsic to stars themselves, e.g. nuisance signals arising from surface magnetic activity such as rotating spots and plages, which can drown out or even mimic planetary RV signals. We use Kepler-37 - know… ▽ More To date, only 18 exoplanets with radial velocity (RV) semi-amplitudes $<2$ m/s have had their masses directly constrained. The biggest obstacle to RV detection of such exoplanets is variability intrinsic to stars themselves, e.g. nuisance signals arising from surface magnetic activity such as rotating spots and plages, which can drown out or even mimic planetary RV signals. We use Kepler-37 - known to host three transiting planets, one of which, Kepler-37d, should be on the cusp of RV detectability with modern spectrographs - as a case study in disentangling planetary and stellar activity signals. We show how two different statistical techniques - one seeking to identify activity signals in stellar spectra, and another to model activity signals in extracted RVs and activity indicators - can enable detection of the hitherto elusive Kepler-37d. Moreover, we show that these two approaches can be complementary, and in combination, facilitate a definitive detection and precise characterisation of Kepler-37d. Its RV semi-amplitude of $1.22\pm0.31$ m/s (mass $5.4\pm1.4$ $M_\oplus$) is formally consistent with TOI-178b's $1.05^{+0.25}_{-0.30}$ m/s, the latter being the smallest detected RV signal of any transiting planet to date, though dynamical simulations suggest Kepler-37d's mass may be on the lower end of our $1σ$ credible interval. Its consequent density is consistent with either a water-world or that of a gaseous envelope ($\sim0.4\%$ by mass) surrounding a rocky core. Based on RV modelling and a re-analysis of Kepler-37 TTVs, we also argue that the putative (non-transiting) planet Kepler-37e should probably be stripped of its 'confirmed' status. △ Less

Submitted 30 July, 2021; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS. 22 pages, 10 figures, 13 tables

Planet Hunters TESS III: two transiting planets around the bright G dwarf HD 152843

Authors: Nora L. Eisner, Belinda A. Nicholson, Oscar Barragán, Suzanne Aigrain, Chris Lintott, Laurel Kaye, Baptiste Klein, Grant Miller, Jake Taylor, Norbert Zicher, Lars A. Buchhave, Douglas A. Caldwell, Jonti Horner, Joe Llama, Annelies Mortier, Vinesh M. Rajpaul, Keivan Stassun, Avi Sporer, Andrew Tkachenko, Jon M. Jenkins, David W. Latham, George R. Ricker, Sara Seager, Joshua N. Winn, Safaa Alhassan , et al. (15 additional authors not shown)

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

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS (15 pages, 10 figures, 3 tables)

TESS re-observes the young multi-planet system TOI-451: refined ephemeris and activity evolution

Authors: Oscar Barragán, Suzanne Aigrain, Edward Gillen, Fernando Gutiérrez-Canales

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

Submitted 17 March, 2021; originally announced March 2021.

Comments: To appear as an RNAAS, 5 pages, 1 Figure

A self-lensing binary massive black hole interpretation of quasi-periodic eruptions

Authors: Adam Ingram, Sara Motta, Suzanne Aigrain, Aris Karastergiou

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

Submitted 26 February, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS

Planet Hunters TESS II: Findings from the first two years of TESS

Authors: Nora L. Eisner, Oscar Barragán, Chris Lintott, Suzanne Aigrain, Belinda Nicholson, Tabetha S. Boyajian, Steve B. Howell, Cole Johnston, Ben Lakeland, Grant Miller, Adam McMaster, Hannu Parviainen, Emily J. Safron, Megan E. Schwamb, Laura Trouille, Sophia Vaughan, Norbert Zicher, Campbell Allen, Sarah Allen, Mark Bouslog, Cliff Johnson, Molly N. Simon, Zach Wolfenbarger, Elisabeth M. L. Baeten, David M. Bundy , et al. (1 additional authors not shown)

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

Submitted 27 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS (22 pages, 12 figures, 3 tables)

Separating planetary reflex Doppler shifts from stellar variability in the wavelength domain

Authors: A. Collier Cameron, E. B. Ford, S. Shahaf, S. Aigrain, X. Dumusque, R. D. Haywood, A. Mortier, D. F. Phillips, L. Buchhave, M. Cecconi, H. Cegla, R. Cosentino, M. Cretignier, A. Ghedina, M. Gonzalez, D. W. Latham, M. Lodi, M. Lopez-Morales, G. Micela, E. Molinari, F. Pepe, G. Piotto, E. Poretti, D. Queloz, J. San Juan , et al. (6 additional authors not shown)

Abstract: Stellar magnetic activity produces time-varying distortions in the photospheric line profiles of solar-type stars. These lead to systematic errors in high-precision radial-velocity measurements, which limit efforts to discover and measure the masses of low-mass exoplanets with orbital periods of more than a few tens of days. We present a new data-driven method for separating Doppler shifts of dyna… ▽ More Stellar magnetic activity produces time-varying distortions in the photospheric line profiles of solar-type stars. These lead to systematic errors in high-precision radial-velocity measurements, which limit efforts to discover and measure the masses of low-mass exoplanets with orbital periods of more than a few tens of days. We present a new data-driven method for separating Doppler shifts of dynamical origin from apparent velocity variations arising from variability-induced changes in the stellar spectrum. We show that the autocorrelation function (ACF) of the cross-correlation function used to measure radial velocities is effectively invariant to translation. By projecting the radial velocities on to a subspace labelled by the observation identifiers and spanned by the amplitude coefficients of the ACF's principal components, we can isolate and subtract velocity perturbations caused by stellar magnetic activity. We test the method on a 5-year time sequence of 853 daily 15-minute observations of the solar spectrum from the HARPS-N instrument and solar-telescope feed on the 3.58-m Telescopio Nazionale Galileo. After removal of the activity signals, the heliocentric solar velocity residuals are found to be Gaussian and nearly uncorrelated. We inject synthetic low-mass planet signals with amplitude $K=40$ cm s$^{-1}$ into the solar observations at a wide range of orbital periods. Projection into the orthogonal complement of the ACF subspace isolates these signals effectively from solar activity signals. Their semi-amplitudes are recovered with a precision of $\sim~6.6$ cm s$^{-1}$, opening the door to Doppler detection and characterization of terrestrial-mass planets around well-observed, bright main-sequence stars across a wide range of orbital periods. △ Less

Submitted 6 May, 2021; v1 submitted 30 October, 2020; originally announced November 2020.

Comments: 19 pages, 15 figures, accepted 2 May 2021 for publication in MNRAS

How Does Thermal Scattering Shape the Infrared Spectra of Cloudy Exoplanets? A Theoretical Framework and Consequences for Atmospheric Retrievals in the JWST era

Authors: Jake Taylor, Vivien Parmentier, Michael R. Line, Elspeth K. H. Lee, Patrick G. J. Irwin, Suzanne Aigrain

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

Submitted 1 July, 2021; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: Accepted to MNRAS

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

Authors: John Stauffer, Davide Barrado, Trevor David, Luisa Rebull, Lynne Hillenbrand, Eric Mamajek, Rebecca Oppenheimer, Suzanne Aigrain, Herve Bouy, Jorge Lillo-Box

Abstract: 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 see… ▽ More 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. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: Accepted by AJ

The spectral impact of magnetic activity on disk-integrated HARPS-N solar observations: exploring new activity indicators

Authors: A. P. G. Thompson, C. A. Watson, R. D. Haywood, J. C. Costes, E. de Mooij, A. Collier Cameron, X. Dumusque, D. F. Phillips, S. H. Saar, A. Mortier, T. W. Milbourne, S. Aigrain, H. M. Cegla, D. Charbonneau, R. Cosentino, A. Ghedina, D. W. Latham, M. López-Morales, G. Micela, E. Molinari, E. Poretti, A. Sozzetti, S. Thompson, R. Walsworth

Abstract: Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the… ▽ More Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the HARPS-N solar telescope. We divide high-activity spectral echelle orders by low-activity master templates (as defined using both log R'HK and images from the Solar Dynamics Observatory, SDO), creating "relative spectra". With resolved images of the surface of the Sun (via SDO), the faculae and spot filling factors can be calculated, giving a measure of activity independent of, and in addition to, log R'HK. We find pseudo-emission (and pseudo-absorption) features in the relative spectra that are similar to those reported in our previous work on alpha Cen B. In alpha Cen B, the features are shown to correlate better to changes in faculae filling factor than spot filling factor. In this work we more confidently identify changes in faculae coverage of the visible hemisphere of the Sun as the source of features produced in the relative spectra. Finally, we produce trailed spectra to observe the RV component of the features, which show that the features move in a redward direction as one would expect when tracking active regions rotating on the surface of a star. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 13 pages, 10 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society

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

Authors: Edward Gillen, Lynne A. Hillenbrand, John Stauffer, Suzanne Aigrain, Luisa Rebull, Ann Marie Cody

Abstract: 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… ▽ More 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. △ Less

Submitted 9 April, 2020; originally announced April 2020.

Comments: 18 pages, 14 figures, 5 tables, accepted for publication in MNRAS

Understanding and Mitigating Biases when Studying Inhomogeneous Emission Spectra with JWST

Authors: Jake Taylor, Vivien Parmentier, Patrick G. J. Irwin, Suzanne Aigrain, Elspeth K. H. Lee, Joshua Krissansen-Totton

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

Submitted 9 February, 2020; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: Accepted for publication by MNRAS

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

Authors: Vinesh M. Rajpaul, Suzanne Aigrain, Lars A. Buchhave

Abstract: 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 improv… ▽ More 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. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: 27 pages, 16 figures. Accepted for publication in MNRAS

Planet Hunters TESS I: TOI 813, a subgiant hosting a transiting Saturn-sized planet on an 84-day orbit

Authors: N. L. Eisner, O. Barragán, S. Aigrain, C. Lintott, G. Miller, N. Zicher, T. S. Boyajian, C. Briceño, E. M. Bryant, J. L. Christiansen, A. D. Feinstein, L. M. Flor-Torres, M. Fridlund, D. Gandolfi, J. Gilbert, N. Guerrero, J. M. Jenkins, K. Jones, M. H. Kristiansen, A. Vanderburg, N. Law, A. R. López-Sánchez, A. W. Mann, E. J. Safron, M. E. Schwamb , et al. (25 additional authors not shown)

Abstract: We report on the discovery and validation of TOI 813b (TIC 55525572 b), a transiting exoplanet identified by citizen scientists in data from NASA's Transiting Exoplanet Survey Satellite (TESS) and the first planet discovered by the Planet Hunters TESS project. The host star is a bright (V = 10.3 mag) subgiant ($R_\star=1.94\,R_\odot$, $M_\star=1.32\,M_\odot$). It was observed almost continuously b… ▽ More We report on the discovery and validation of TOI 813b (TIC 55525572 b), a transiting exoplanet identified by citizen scientists in data from NASA's Transiting Exoplanet Survey Satellite (TESS) and the first planet discovered by the Planet Hunters TESS project. The host star is a bright (V = 10.3 mag) subgiant ($R_\star=1.94\,R_\odot$, $M_\star=1.32\,M_\odot$). It was observed almost continuously by TESS during its first year of operations, during which time four individual transit events were detected. The candidate passed all the standard light curve-based vetting checks, and ground-based follow-up spectroscopy and speckle imaging enabled us to place an upper limit of $2 M_{Jup}$ (99 % confidence) on the mass of the companion, and to statistically validate its planetary nature. Detailed modelling of the transits yields a period of $83.8911_{ - 0.0031 } ^ { + 0.0027 }$ days, a planet radius of $6.71 \pm 0.38$ $R_{\oplus}$, and a semi major axis of $0.423_{ - 0.037 } ^ { + 0.031 }$ AU. The planet's orbital period combined with the evolved nature of the host star places this object in a relatively under-explored region of parameter space. We estimate that TOI-813b induces a reflex motion in its host star with a semi-amplitude of $\sim6$ ms$^{-1}$, making this system a promising target to measure the mass of a relatively long-period transiting planet. △ Less

Submitted 13 January, 2020; v1 submitted 19 September, 2019; originally announced September 2019.

Comments: Accepted for publication in MNRAS (16 pages, 10 figures, 3 tables)

Radial velocity confirmation of K2-100b: a young, highly irradiated, and low density transiting hot Neptune

Authors: O. Barragán, S. Aigrain, D. Kubyshkina, D. Gandolfi, J. Livingston, M. C. V. Fridlund, L. Fossati, J. Korth, H. Parviainen, L. Malavolta, E. Palle, H. J. Deeg, G. Nowak, V. M. Rajpaul, N. Zicher, G. Antoniciello, N. Narita, S. Albrecht, L. R. Bedin, J. Cabrera, W. D. Cochran, J. de Leon, Ph. Eigmüller, A. Fukui, V. Granata , et al. (15 additional authors not shown)

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

Submitted 11 September, 2019; originally announced September 2019.

Comments: Accepted for publications in MNRAS

2.5-D retrieval of atmospheric properties from exoplanet phase curves: Application to WASP-43b observations

Authors: Patrick G. J. Irwin, Vivien Parmentier, Jake Taylor, Jo Barstow, Suzanne Aigrain, Elspeth K. H. Lee, Ryan Garland

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

Submitted 9 December, 2019; v1 submitted 7 September, 2019; originally announced September 2019.

Comments: Submitted to M.N.R.A.S. Revised Submission 9th December 2019

An 11 Earth-Mass, Long-Period Sub-Neptune Orbiting a Sun-like Star

Authors: Andrew W. Mayo, Vinesh M. Rajpaul, Lars A. Buchhave, Courtney D. Dressing, Annelies Mortier, Li Zeng, Charles D. Fortenbach, Suzanne Aigrain, Aldo S. Bonomo, Andrew Collier Cameron, David Charbonneau, Adrien Coffinet, Rosario Cosentino, Mario Damasso, Xavier Dumusque, A. F. Martinez Fiorenzano, Raphaëlle D. Haywood, David W. Latham, Mercedes López-Morales, Luca Malavolta, Giusi Micela, Emilio Molinari, Logan Pearce, Francesco Pepe, David Phillips , et al. (5 additional authors not shown)

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

Submitted 11 December, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: Published in AJ, 17 pages, 8 figures, 3 tables

Journal ref: 2019 AJ 158 4

The K2 Bright Star Survey I: Methodology and Data Release

Authors: Benjamin J. S. Pope, Timothy R. White, Will M. Farr, Jie Yu, Michael Greklek-McKeon, Daniel Huber, Conny Aerts, Suzanne Aigrain, Timothy R. Bedding, Tabetha Boyajian, Orlagh L. Creevey, David W. Hogg

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

Submitted 26 August, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Comments: Accepted ApJS. Paper documents data release of K2 Halo light curves, available from http://archive.stsci.edu/doi/resolve/resolve.html?doi=10.17909/t9-6wj4-eb32 or https://github.com/benjaminpope/k2halo. v2: corrected typo

The Search for Living Worlds and the Connection to Our Cosmic Origins

Authors: M. A. Barstow, S. Aigrain, J. Barstow, M. Barthelemy, B. Biller, A. Bonanos, L. Buchhave, S. Casewell, C. Charbonnel, S. Charlot, R. Davies, N. Devaney, C. Evans, M. Ferrari, L. Fossatti, B. Gaensicke, M. Garcia, A. Gomez de Castro, T. Henning, C. Lintott, C. Knigge, C. Neiner, L. Rossi, C. Snodgrass, D. Stam , et al. (2 additional authors not shown)

Abstract: One of the most exciting scientific challenges is to detect Earth-like planets in the habitable zones of other stars in the galaxy and search for evidence of life. The ability to observe and characterise dozens of potentially Earth-like planets now lies within the realm of possibility due to rapid advances in key space and imaging technologies. The associated challenge of directly imaging very f… ▽ More One of the most exciting scientific challenges is to detect Earth-like planets in the habitable zones of other stars in the galaxy and search for evidence of life. The ability to observe and characterise dozens of potentially Earth-like planets now lies within the realm of possibility due to rapid advances in key space and imaging technologies. The associated challenge of directly imaging very faint planets in orbit around nearby very bright stars is now well understood, with the key instrumentation also being perfected and developed. Such advances will allow us to develop large transformative telescopes, covering a broad UV-optical-IR spectral range, which can carry out the detailed research programmes designed to answer the questions we wish to answer: Carry out high contrast imaging surveys of nearby stars to search for planets within their habitable zones. Characterise the planets detected to determine masses and radii from photometric measurements. Through spectroscopic studies of their atmospheres and surfaces, search for habitability indicators and for signs of an environment that has been modified by the presence of life. Active studies of potential missions have been underway for a number of years. The latest of these is the Large UV Optical IR space telescope (LUVOIR), one of four flagship mission studies commissioned by NASA in support of the 2020 US Decadal Survey. LUVOIR, if selected, will be of interest to a wide scientific community and will be the only telescope capable of searching for and characterizing a sufficient number of exoEarths to provide a meaningful answer to the question - Are we alone?. This paper is a submission to the European Space Agency Voyage 2050 call for white papers outlining the case for an ESA contribution to a Large UVOIR telescope. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: Submission to the European Space Agency Voyage 2050 call for white papers

Constraining the properties of HD 206893 B. A combination of radial velocity, direct imaging, and astrometry data

Authors: A. Grandjean, A. -M. Lagrange, H. Beust, L. Rodet, J. Milli, P. Rubini, C. Babusiaux, N. Meunier, P. Delorme, S. Aigrain, N. Zicher, M. Bonnefoy, B. A. Biller, J. -L. Baudino, M. Bonavita, A. Boccaletti, A. Cheetham, J. H. Girard, J. Hagelberg, M. Janson, J. Lannier, C. Lazzoni, R. Ligi, A. -L. Maire, D. Mesa , et al. (3 additional authors not shown)

Abstract: High contrast imaging enables the determination of orbital parameters for substellar companions (planets, brown dwarfs) from the observed relative astrometry and the estimation of model and age-dependent masses from their observed magnitudes or spectra. Combining astrometric positions with radial velocity gives direct constraints on the orbit and on the dynamical masses of companions. A brown dwar… ▽ More High contrast imaging enables the determination of orbital parameters for substellar companions (planets, brown dwarfs) from the observed relative astrometry and the estimation of model and age-dependent masses from their observed magnitudes or spectra. Combining astrometric positions with radial velocity gives direct constraints on the orbit and on the dynamical masses of companions. A brown dwarf was discovered with the VLT/SPHERE instrument in 2017, which orbits at $\sim$ 11 au around HD 206893. Its mass was estimated between 12 and 50 $M_{Jup}$ from evolutionary models and its photometry. However, given the significant uncertainty on the age of the system and the peculiar spectrophotometric properties of the companion, this mass is not well constrained. We aim at constraining the orbit and dynamical mass of HD 206893 B. We combined radial velocity data obtained with HARPS spectra and astrometric data obtained with the high contrast imaging VLT/SPHERE and VLT/NaCo instruments, with a time baseline less than three years. We then combined those data with astrometry data obtained by Hipparcos and Gaia with a time baseline of 24 years. We used a MCMC approach to estimate the orbital parameters and dynamical mass of the brown dwarf from those data. We infer a period between 21 and 33° and an inclination in the range 20-41° from pole-on from HD 206893 B relative astrometry. The RV data show a significant RV drift over 1.6 yrs. We show that HD 206893 B cannot be the source of this observed RV drift as it would lead to a dynamical mass inconsistent with its photometry and spectra and with Hipparcos and Gaia data. An additional inner (semimajor axis in the range 1.4-2.6 au) and massive ($\sim$ 15 $M_{Jup}$) companion is needed to explain the RV drift, which is compatible with the available astrometric data of the star, as well as with the VLT/SPHERE and VLT/NaCo nondetection. △ Less

Submitted 24 July, 2019; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: 14 pages, 16 figures

The Kepler Smear Campaign: Light curves for 102 Very Bright Stars

Authors: Benjamin J. S. Pope, Guy R. Davies, Keith Hawkins, Timothy R. White, Amalie Stokholm, Allyson Bieryla, David W. Latham, Madeline Lucey, Conny Aerts, Suzanne Aigrain, Victoria Antoci, Timothy R. Bedding, Dominic M. Bowman, Douglas A. Caldwell, Ashley Chontos, Gilbert A. Esquerdo, Daniel Huber, Paula Jofre, Simon J. Murphy, Timothy van Reeth, Victor Silva Aguirre, Jie Yu

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

Submitted 23 May, 2019; originally announced May 2019.

Comments: 35 pages, accepted ApJS

Stellar activity and rotation of the planet host Kepler-17 from long-term space-borne photometry

Authors: A. F. Lanza, Y. Netto, A. S. Bonomo, H. Parviainen, A. Valio, S. Aigrain

Abstract: The study of young Sun-like stars is of fundamental importance to understand the magnetic activity and rotational evolution of the Sun. Space-borne photometry by the Kepler telescope provides unprecedented datasets to investigate these phenomena in Sun-like stars. We present a new analysis of the entire Kepler photometric time series of the moderately young Sun-like star Kepler-17 that is accompan… ▽ More The study of young Sun-like stars is of fundamental importance to understand the magnetic activity and rotational evolution of the Sun. Space-borne photometry by the Kepler telescope provides unprecedented datasets to investigate these phenomena in Sun-like stars. We present a new analysis of the entire Kepler photometric time series of the moderately young Sun-like star Kepler-17 that is accompanied by a transiting hot Jupiter. We applied a maximum-entropy spot model to the long-cadence out-of-transit photometry of the target to derive maps of the starspot filling factor versus the longitude and the time. These maps are compared to the spots occulted during transits to validate our reconstruction and derive information on the latitudes of the starspots. We find two main active longitudes on the photosphere of Kepler-17, one of which has a lifetime of at least $\sim 1400$ days, although with a varying level of activity. The latitudinal differential rotation is of solar type, that is, with the equator rotating faster than the poles. We estimate a minimum relative amplitude $ΔΩ/ Ω$ between $\sim 0.08 \pm 0.05$ and $0.14 \pm 0.05$, our determination being affected by the finite lifetime of individual starspots and depending on the adopted spot model parameters. We find marginal evidence of a short-term intermittent activity cycle of $\sim 48$ days and an indication of a longer cycle of $400-600$ days characterized by an equatorward migration of the mean latitude of the spots as in the Sun. The rotation of Kepler-17 is likely to be significantly affected by the tides raised by its massive close-by planet. We confirm the reliability of maximum-entropy spot models to map starspots in young active stars and characterize the activity and differential rotation of this young Sun-like planetary host. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted by Astronomy and Astrophysics - 22 pages, 29 figure, 1 table, 2 appendixes

Journal ref: A&A 626, A38 (2019)

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

Authors: Rebecca Esselstein, Suzanne Aigrain, Andrew Vanderburg, Jeffrey C. Smith, Soren Meibom, Jennifer Van Saders, Robert Mathieu

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 instab… ▽ More 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. △ Less

Submitted 6 June, 2018; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: 55 pages, 22 figures, published by ApJ

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

Authors: M. Deleuil, S. Aigrain, C. Moutou, J. Cabrera, F. Bouchy, H. J. Deeg, J. -M. Almenara, G. Hébrard, A. Santerne, R. Alonso, A. S. Bonomo, P. Bordé, Sz. Csizmadia, A. Erikson, M. Fridlund, D. Gandolfi, E. Guenther, T. Guillot, P. Guterman, S. Grziwa, A. Hatzes, A. Léger, T. Mazeh, A. Ofir, M. Ollivier , et al. (8 additional authors not shown)

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: de… ▽ More 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. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Comments: Full tables will be provided online at CDS

Journal ref: A&A 619, A97 (2018)

A universal spin-mass relation for brown dwarfs and planets

Authors: Aleks Scholz, Keavin Moore, Ray Jayawardhana, Suzanne Aigrain, Dawn Peterson, Beate Stelzer

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 t… ▽ More 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. △ Less

Submitted 24 April, 2018; v1 submitted 19 April, 2018; originally announced April 2018.

Comments: 20 pages, 6 figures, accepted for publication in ApJ

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

Authors: Eric L. Sandquist, Robert D. Mathieu, Samuel N. Quinn, Maxwell L. Pollack, David W. Latham, Timothy M. Brown, Rebecca Esselstein, Suzanne Aigrain, Hannu Parviainen, Andrew Vanderburg, Dennis Stello, Garrett Somers, Marc H. Pinsonneault, Jamie Tayar, Jerome A. Orosz, Luigi R. Bedin, Mattia Libralato, Luca Malavolta, Domenico Nardiello

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 (… ▽ More 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. △ Less

Submitted 16 February, 2018; originally announced February 2018.

Comments: 33 pages, 11 figures, accepted to AJ, photometry files will be available with the electronic journal article