-
The death of Vulcan: NEID reveals the planet candidate orbiting HD 26965 is stellar activity
Authors:
Abigail Burrows,
Samuel Halverson,
Jared C. Siegel,
Christian Gilbertson,
Jacob Luhn,
Jennifer Burt,
Chad F. Bender,
Arpita Roy,
Ryan C. Terrien,
Selma Vangstein,
Suvrath Mahadevan,
Jason T. Wright,
Paul Robertson,
Eric B. Ford,
Guðmundur Stefánsson,
Joe P. Ninan,
Cullen H. Blake,
Michael W. McElwain,
Christian Schwab,
Jinglin Zhao
Abstract:
We revisit the long-studied radial velocity (RV) target HD26965 using recent observations from the NASA-NSF 'NEID' precision Doppler facility. Leveraging a suite of classical activity indicators, combined with line-by-line RV analyses, we demonstrate that the claimed 45-day signal previously identified as a planet candidate is most likely an activity-induced signal. Correlating the bulk (spectrall…
▽ More
We revisit the long-studied radial velocity (RV) target HD26965 using recent observations from the NASA-NSF 'NEID' precision Doppler facility. Leveraging a suite of classical activity indicators, combined with line-by-line RV analyses, we demonstrate that the claimed 45-day signal previously identified as a planet candidate is most likely an activity-induced signal. Correlating the bulk (spectrally-averaged) RV with canonical line activity indicators confirms a multi-day 'lag' between the observed activity indicator time series and the measured RV. When accounting for this lag, we show that much of the observed RV signal can be removed by a linear detrending of the data. Investigating activity at the line-by-line level, we find a depth-dependent correlation between individual line RVs and the bulk RVs, further indicative of periodic suppression of convective blueshift causing the observed RV variability, rather than an orbiting planet. We conclude that the combined evidence of the activity correlations and depth dependence is consistent with a radial velocity signature dominated by a rotationally-modulated activity signal at a period of $\sim$42 days. We hypothesize that this activity signature is due to a combination of spots and convective blueshift suppression. The tools applied in our analysis are broadly applicable to other stars, and could help paint a more comprehensive picture of the manifestations of stellar activity in future Doppler RV surveys.
△ Less
Submitted 26 April, 2024;
originally announced April 2024.
-
Developing a Drift Rate Distribution for Technosignature Searches of Exoplanets
Authors:
Megan G. Li,
Sofia Z. Sheikh,
Christian Gilbertson,
Matthias Y. He,
Howard Isaacson,
Steve Croft,
Evan L. Sneed
Abstract:
A stable-frequency transmitter with relative radial acceleration to a receiver will show a change in received frequency over time, known as a "drift rate''. For a transmission from an exoplanet, we must account for multiple components of drift rate: the exoplanet's orbit and rotation, the Earth's orbit and rotation, and other contributions. Understanding the drift rate distribution produced by exo…
▽ More
A stable-frequency transmitter with relative radial acceleration to a receiver will show a change in received frequency over time, known as a "drift rate''. For a transmission from an exoplanet, we must account for multiple components of drift rate: the exoplanet's orbit and rotation, the Earth's orbit and rotation, and other contributions. Understanding the drift rate distribution produced by exoplanets relative to Earth, can a) help us constrain the range of drift rates to check in a Search for Extraterrestrial Intelligence (SETI) project to detect radio technosignatures and b) help us decide validity of signals-of-interest, as we can compare drifting signals with expected drift rates from the target star. In this paper, we modeled the drift rate distribution for $\sim$5300 confirmed exoplanets, using parameters from the NASA Exoplanet Archive (NEA). We find that confirmed exoplanets have drift rates such that 99\% of them fall within the $\pm$53 nHz range. This implies a distribution-informed maximum drift rate $\sim$4 times lower than previous work. To mitigate the observational biases inherent in the NEA, we also simulated an exoplanet population built to reduce these biases. The results suggest that, for a Kepler-like target star without known exoplanets, $\pm$0.44 nHz would be sufficient to account for 99\% of signals. This reduction in recommended maximum drift rate is partially due to inclination effects and bias towards short orbital periods in the NEA. These narrowed drift rate maxima will increase the efficiency of searches and save significant computational effort in future radio technosignature searches.
△ Less
Submitted 2 November, 2023;
originally announced November 2023.
-
A Green Bank Telescope search for narrowband technosignatures between 1.1-1.9 GHz during 12 Kepler planetary transits
Authors:
Sofia Z. Sheikh,
Shubham Kanodia,
Emily Lubar,
William P. Bowman,
Caleb I. Cañas,
Christian Gilbertson,
Mariah G. MacDonald,
Jason Wright,
David MacMahon,
Steve Croft,
Danny Price,
Andrew Siemion,
Jamie Drew,
S. Pete Worden,
Elizabeth Trenholm
Abstract:
A growing avenue for determining the prevalence of life beyond Earth is to search for "technosignatures" from extraterrestrial intelligences/agents. Technosignatures require significant energy to be visible across interstellar space and thus intentional signals might be concentrated in frequency, in time, or in space, to be found in mutually obvious places. Therefore, it could be advantageous to s…
▽ More
A growing avenue for determining the prevalence of life beyond Earth is to search for "technosignatures" from extraterrestrial intelligences/agents. Technosignatures require significant energy to be visible across interstellar space and thus intentional signals might be concentrated in frequency, in time, or in space, to be found in mutually obvious places. Therefore, it could be advantageous to search for technosignatures in parts of parameter space that are mutually-derivable to an observer on Earth and a distant transmitter. In this work, we used the L-band (1.1-1.9 GHz) receiver on the Robert C. Byrd Green Bank Telescope (GBT) to perform the first technosignature search pre-synchronized with exoplanet transits, covering 12 Kepler systems. We used the Breakthrough Listen turboSETI pipeline to flag narrowband hits ($\sim$3 Hz) using a maximum drift rate of $\pm$614.4 Hz/s and a signal-to-noise threshold of 5 - the pipeline returned $\sim 3.4 \times 10^5$ apparently-localized features. Visual inspection by a team of citizen scientists ruled out 99.6% of them. Further analysis found 2 signals-of-interest that warrant follow-up, but no technosignatures. If the signals-of-interest are not re-detected in future work, it will imply that the 12 targets in the search are not producing transit-aligned signals from 1.1-1.9 GHz with transmitter powers $>$60 times that of the former Arecibo radar. This search debuts a range of innovative technosignature techniques: citizen science vetting of potential signals-of-interest, a sensitivity-aware search out to extremely high drift rates, a more flexible method of analyzing on-off cadences, and an extremely low signal-to-noise threshold.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Impact of Correlated Noise on the Mass Precision of Earth-analog Planets in Radial Velocity Surveys
Authors:
Jacob K. Luhn,
Eric B. Ford,
Zhao Guo,
Christian Gilbertson,
Patrick Newman,
Peter Plavchan,
Jennifer A. Burt,
Johanna Teske,
Arvind F. Gupta
Abstract:
Characterizing the masses and orbits of near-Earth-mass planets is crucial for interpreting observations from future direct imaging missions (e.g., HabEx, LUVOIR). Therefore, the Exoplanet Science Strategy report (National Academies of Sciences, Engineering, and Medicine 2018) recommended further research so future extremely precise radial velocity surveys could contribute to the discovery and/or…
▽ More
Characterizing the masses and orbits of near-Earth-mass planets is crucial for interpreting observations from future direct imaging missions (e.g., HabEx, LUVOIR). Therefore, the Exoplanet Science Strategy report (National Academies of Sciences, Engineering, and Medicine 2018) recommended further research so future extremely precise radial velocity surveys could contribute to the discovery and/or characterization of near-Earth-mass planets in the habitable zones of nearby stars prior to the launch of these future imaging missions. Newman et al. (2021) simulated such 10-year surveys under various telescope architectures, demonstrating they can precisely measure the masses of potentially habitable Earth-mass planets in the absence of stellar variability. Here, we investigate the effect of stellar variability on the signal-to-noise ratio (SNR) of the planet mass measurements in these simulations. We find that correlated noise due to active regions has the largest effect on the observed mass SNR, reducing the SNR by a factor of $\sim$5.5 relative to the no-variability scenario -- granulation reduces by a factor of $\sim$3, while p-mode oscillations has little impact on the proposed survey strategies. We show that in the presence of correlated noise, 5-cm s$^{-1}$ instrumental precision offers little improvement over 10-cm s$^{-1}$ precision, highlighting the need to mitigate astrophysical variability. With our noise models, extending the survey to 15 years doubles the number of Earth-analogs with mass SNR $>$ 10, and reaching this threshold for any Earth-analog orbiting a star $>$ 0.76 M$_{\odot}$ in a 10-year survey would require an increase in number of observations per star from that in Newman et al. (2021).
△ Less
Submitted 26 April, 2022;
originally announced April 2022.
-
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.
-
A Stellar Activity F-statistic for Exoplanet Surveys (SAFE)
Authors:
Parker H. Holzer,
Jessi Cisewski-Kehe,
Lily Zhao,
Eric B. Ford,
Christian Gilbertson,
Debra A. Fischer
Abstract:
In the search for planets orbiting distant stars the presence of stellar activity in the atmospheres of observed stars can obscure the radial velocity signal used to detect such planets. Furthermore, this stellar activity contamination is set by the star itself and cannot simply be avoided with better instrumentation. Various stellar activity indicators have been developed that may correlate with…
▽ More
In the search for planets orbiting distant stars the presence of stellar activity in the atmospheres of observed stars can obscure the radial velocity signal used to detect such planets. Furthermore, this stellar activity contamination is set by the star itself and cannot simply be avoided with better instrumentation. Various stellar activity indicators have been developed that may correlate with this contamination. We introduce a new stellar activity indicator called the Stellar Activity F-statistic for Exoplanet surveys (SAFE) that has higher statistical power (i.e., probability of detecting a true stellar activity signal) than many traditional stellar activity indicators in a simulation study of an active region on a Sun-like star with moderate to high signal-to-noise. Also through simulation, the SAFE is demonstrated to be associated with the projected area on the visible side of the star covered by active regions. We also demonstrate that the SAFE detects statistically significant stellar activity in most of the spectra for HD 22049, a star known to have high stellar variability. Additionally, the SAFE is calculated for recent observations of the three low-variability stars HD 34411, HD 10700, and HD 3651, the latter of which is known to have a planetary companion. As expected, the SAFE for these three only occasionally detects activity. Furthermore, initial exploration appears to indicate that the SAFE may be useful for disentangling stellar activity signals from planet-induced Doppler shifts.
△ Less
Submitted 10 April, 2021;
originally announced April 2021.
-
Stability Constrained Characterization of Multiplanet Systems
Authors:
Daniel Tamayo,
Christian Gilbertson,
Daniel Foreman-Mackey
Abstract:
Many discovered multiplanet systems are tightly packed. This implies that wide parameter ranges in masses and orbital elements can be dynamically unstable and ruled out. We present a case study of Kepler-23, a compact three-planet system where constraints from stability, transit timing variations (TTVs), and transit durations can be directly compared. We find that in this tightly packed system, st…
▽ More
Many discovered multiplanet systems are tightly packed. This implies that wide parameter ranges in masses and orbital elements can be dynamically unstable and ruled out. We present a case study of Kepler-23, a compact three-planet system where constraints from stability, transit timing variations (TTVs), and transit durations can be directly compared. We find that in this tightly packed system, stability can place upper limits on the masses and orbital eccentricities of the bodies that are comparable to or tighter than current state of the art methods. Specifically, stability places 68% upper limits on the orbital eccentricities of 0.09, 0.04, and 0.05 for planets $b$, $c$ and $d$, respectively. These constraints correspond to radial velocity signals $\lesssim 20$ cm/s, are significantly tighter to those from transit durations, and comparable to those from TTVs. Stability also yields 68% upper limits on the masses of planets $b$, $c$ and $d$ of 2.2, 16.1, and 5.8 $M_\oplus$, respectively, which were competitive with TTV constraints for the inner and outer planets. Performing this stability constrained characterization is computationally expensive with N-body integrations. We show that SPOCK, the Stability of Planetary Orbital Configurations Klassifier, is able to faithfully approximate the N-body results over 4000 times faster. We argue that such stability constrained characterization of compact systems is a challenging "needle-in-a-haystack" problem (requiring removal of 2500 unstable configurations for every stable one for our adopted priors) and we offer several practical recommendations for such stability analyses.
△ Less
Submitted 1 December, 2020; v1 submitted 24 September, 2020;
originally announced September 2020.
-
Toward Extremely Precise Radial Velocities: II. A Tool For Using Multivariate Gaussian Processes to Model Stellar Activity
Authors:
Christian Gilbertson,
Eric B. Ford,
David E. Jones,
David C. Stenning
Abstract:
The radial velocity method is one of the most successful techniques for the discovery and characterization of exoplanets. Modern spectrographs promise measurement precision of ~0.2-0.5 m/s for an ideal target star. However, the intrinsic variability of stellar spectra can mimic and obscure true planet signals at these levels. Rajpaul et al. (2015) and Jones et al. (2017) proposed applying a physic…
▽ More
The radial velocity method is one of the most successful techniques for the discovery and characterization of exoplanets. Modern spectrographs promise measurement precision of ~0.2-0.5 m/s for an ideal target star. However, the intrinsic variability of stellar spectra can mimic and obscure true planet signals at these levels. Rajpaul et al. (2015) and Jones et al. (2017) proposed applying a physically motivated, multivariate Gaussian process (GP) to jointly model the apparent Doppler shift and multiple indicators of stellar activity as a function of time, so as to separate the planetary signal from various forms of stellar variability. These methods are promising, but performing the necessary calculations can be computationally intensive and algebraically tedious. In this work, we present a flexible and computationally efficient software package, GPLinearOdeMaker.jl, for modeling multivariate time series using a linear combination of univariate GPs and their derivatives. The package allows users to easily and efficiently apply various multivariate GP models and different covariance kernel functions. We demonstrate GPLinearOdeMaker.jl by applying the Jones et al. (2017) model to fit measurements of the apparent Doppler shift and activity indicators derived from simulated active solar spectra time series affected by many evolving star spots. We show how GPLinearOdeMaker.jl makes it easy to explore the effect of different choices for the GP kernel. We find that local kernels could significantly increase the sensitivity and precision of Doppler planet searches relative to the widely used quasi-periodic kernel.
△ Less
Submitted 18 January, 2021; v1 submitted 2 September, 2020;
originally announced September 2020.
-
Predicting the long-term stability of compact multiplanet systems
Authors:
Daniel Tamayo,
Miles Cranmer,
Samuel Hadden,
Hanno Rein,
Peter Battaglia,
Alysa Obertas,
Philip J. Armitage,
Shirley Ho,
David Spergel,
Christian Gilbertson,
Naireen Hussain,
Ari Silburt,
Daniel Jontof-Hutter,
Kristen Menou
Abstract:
We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations…
▽ More
We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.
△ Less
Submitted 14 July, 2020; v1 submitted 13 July, 2020;
originally announced July 2020.
-
Towards Extremely Precise Radial Velocities: I. Simulated Solar Spectra for Testing Exoplanet Detection Algorithms
Authors:
Christian Gilbertson,
Eric B. Ford,
Xavier Dumusque
Abstract:
Recent and upcoming stabilized spectrographs are pushing the frontier for Doppler spectroscopy to detect and characterize low-mass planets. Specifications for these instruments are so impressive that intrinsic stellar variability is expected to limit their Doppler precision for most target stars (Fischer et al. 2016). To realize their full potential, astronomers must develop new strategies for dis…
▽ More
Recent and upcoming stabilized spectrographs are pushing the frontier for Doppler spectroscopy to detect and characterize low-mass planets. Specifications for these instruments are so impressive that intrinsic stellar variability is expected to limit their Doppler precision for most target stars (Fischer et al. 2016). To realize their full potential, astronomers must develop new strategies for distinguishing true Doppler shifts from intrinsic stellar variability. Stellar variability due to star spots, faculae and other rotationally-linked variability are particularly concerning, as the stellar rotation period is often included in the range of potential planet orbital periods. To robustly detect and accurately characterize low-mass planets via Doppler planet surveys, the exoplanet community must develop statistical models capable of jointly modeling planetary perturbations and intrinsic stellar variability. Towards this effort, this note presents simulations of extremely high resolution, solar-like spectra created with SOAP 2.0 (arXiv:1409.3594) that includes multiple evolving star spots. We anticipate this data set will contribute to future studies developing, testing, and comparing statistical methods for measuring physical radial velocities amid contamination by stellar variability.
△ Less
Submitted 30 April, 2020;
originally announced May 2020.
-
Improving Exoplanet Detection Power: Multivariate Gaussian Process Models for Stellar Activity
Authors:
David E. Jones,
David C. Stenning,
Eric B. Ford,
Robert L. Wolpert,
Thomas J. Loredo,
Christian Gilbertson,
Xavier Dumusque
Abstract:
The radial velocity method is one of the most successful techniques for detecting exoplanets. It works by detecting the velocity of a host star induced by the gravitational effect of an orbiting planet, specifically the velocity along our line of sight, which is called the radial velocity of the star. Low-mass planets typically cause their host star to move with radial velocities of 1 m/s or less.…
▽ More
The radial velocity method is one of the most successful techniques for detecting exoplanets. It works by detecting the velocity of a host star induced by the gravitational effect of an orbiting planet, specifically the velocity along our line of sight, which is called the radial velocity of the star. Low-mass planets typically cause their host star to move with radial velocities of 1 m/s or less. By analyzing a time series of stellar spectra from a host star, modern astronomical instruments can in theory detect such planets. However, in practice, intrinsic stellar variability (e.g., star spots, convective motion, pulsations) affects the spectra and often mimics a radial velocity signal. This signal contamination makes it difficult to reliably detect low-mass planets. A principled approach to recovering planet radial velocity signals in the presence of stellar activity was proposed by Rajpaul et al. (2015). It uses a multivariate Gaussian process model to jointly capture time series of the apparent radial velocity and multiple indicators of stellar activity. We build on this work in two ways: (i) we propose using dimension reduction techniques to construct new high-information stellar activity indicators; and (ii) we extend the Rajpaul et al. (2015) model to a larger class of models and use a power-based model comparison procedure to select the best model. Despite significant interest in exoplanets, previous efforts have not performed large-scale stellar activity model selection or attempted to evaluate models based on planet detection power. In the case of main sequence G2V stars, we find that our method substantially improves planet detection power compared to previous state-of-the-art approaches.
△ Less
Submitted 25 August, 2020; v1 submitted 3 November, 2017;
originally announced November 2017.
-
A Study of the Gamma-Ray Burst Fundamental Plane
Authors:
Maria Dainotti,
Christian Gilbertson,
Sergey Postnikov,
Shigehiro Nagataki,
Richard Willingale
Abstract:
A class of long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak. We extended the original analysis with X-ray data from July 2014 to July 2016 achieving a total…
▽ More
A class of long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak. We extended the original analysis with X-ray data from July 2014 to July 2016 achieving a total sample of 183 Swift GRBs with afterglow plateaus and known redshifts. We added the most recent GRBs to the previous 'gold sample' (now including 45 GRBs) and obtained a relation plane with intrinsic scatter compatible within one sigma with the previous result. We compared several GRB categories, such as short with extended emission, X-ray Flashes, GRBs associated with SNe, long-duration GRBs, and the gold sample, composed only by GRBs with light curves with good data coverage and relatively flat plateaus and evaluated their relation planes. We found that they are not statistically different from the fundamental plane derived from the gold sample and that the fundamental plane still has the smallest scatter. We compared the jet opening angles tabulated in literature with the angles derived using the Eiso-Egamma relation of the method in Pescalli et al. (2015) and calculated the relation plane for a sample of long GRBs accounting for the different jet opening angles. We observed that this correction does not significantly reduce the scatter. In an extended analysis, we found that the fundamental plane is independent from several prompt and afterglow parameters.
△ Less
Submitted 31 October, 2016; v1 submitted 28 October, 2016;
originally announced October 2016.
