Showing 1–50 of 215 results for author: Ford, E B

Paths to Robust Exoplanet Science Yield Margin for the Habitable Worlds Observatory

Authors: Christopher C. Stark, Bertrand Mennesson, Steve Bryson, Eric B. Ford, Tyler D. Robinson, Ruslan Belikov, Matthew R. Bolcar, Lee D. Feinberg, Olivier Guyon, Natasha Latouf, Avi M. Mandell, Bernard J. Rauscher, Dan Sirbu, Noah W. Tuchow

Abstract: The Habitable Worlds Observatory (HWO) will seek to detect and characterize potentially Earth-like planets around other stars. To ensure that the mission achieves the Astro2020 Decadal's recommended goal of 25 exoEarth candidates (EECs), we must take into account the probabilistic nature of exoplanet detections and provide "science margin" to budget for astrophysical uncertainties with a reasonabl… ▽ More The Habitable Worlds Observatory (HWO) will seek to detect and characterize potentially Earth-like planets around other stars. To ensure that the mission achieves the Astro2020 Decadal's recommended goal of 25 exoEarth candidates (EECs), we must take into account the probabilistic nature of exoplanet detections and provide "science margin" to budget for astrophysical uncertainties with a reasonable level of confidence. In this study, we explore the probabilistic distributions of yields to be expected from a blind exoEarth survey conducted by such a mission. We identify and estimate the impact of all major known sources of astrophysical uncertainty on the exoEarth candidate yield. As expected, eta_Earth uncertainties dominate the uncertainty in EEC yield, but we show that sampling uncertainties inherent to a blind survey are another important source of uncertainty that should be budgeted for during mission design. We adopt the Large UV/Optical/IR Surveyor Design B (LUVOIR-B) as a baseline and modify the telescope diameter to estimate the science margin provided by a larger telescope. We then depart from the LUVOIR-B baseline design and identify six possible design changes that, when compiled, provide large gains in exoEarth candidate yield and more than an order of magnitude reduction in exposure times for the highest priority targets. We conclude that a combination of telescope diameter increase and design improvements could provide robust exoplanet science margins for HWO. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 90 pages, 26 figures, 4 tables, Submitted to JATIS

GRASS II: Simulations of Potential Granulation Noise Mitigation Methods

Authors: Michael L. Palumbo III, Eric B. Ford, Elizabeth B. Gonzalez, Jason T. Wright, Khaled Al Moulla, Rolf Schlichenmaier

Abstract: We present an updated version of GRASS (the GRanulation And Spectrum Simulator, Palumbo et al. 2022) which now uses an expanded library of 22 solar lines to empirically model time-resolved spectral variations arising from solar granulation. We show that our synthesis model accurately reproduces disk-integrated solar line profiles and bisectors, and we quantify the intrinsic granulation-driven radi… ▽ More We present an updated version of GRASS (the GRanulation And Spectrum Simulator, Palumbo et al. 2022) which now uses an expanded library of 22 solar lines to empirically model time-resolved spectral variations arising from solar granulation. We show that our synthesis model accurately reproduces disk-integrated solar line profiles and bisectors, and we quantify the intrinsic granulation-driven radial-velocity (RV) variability for each of the 22 lines studied. We show that summary statistics of bisector shape (e.g., bisector inverse slope) are strongly correlated with the measured anomalous, variability-driven RV at high pixel signal-to-noise ratio (SNR) and spectral resolution. Further, the strength of the correlations vary both line by line and with the summary statistic used. These correlations disappear for individual lines at the typical spectral resolutions and SNRs achieved by current EPRV spectrographs; so we use simulations from GRASS to demonstrate that they can, in principle, be recovered by selectively binning lines that are similarly affected by granulation. In the best-case scenario (high SNR and large number of binned lines), we find that a $\lesssim$30$\%$ reduction in the granulation-induced root mean square (RMS) RV can be achieved, but that the achievable reduction in variability is most strongly limited by the spectral resolution of the observing instrument. Based on our simulations, we predict that existing ultra-high-resolution spectrographs, namely ESPRESSO and PEPSI, should be able to resolve convective variability in other, bright stars. △ Less

Submitted 14 May, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Comments: 20 + 16 pages, 6 + 11 figures, 4 tables; in revision for AJ

Journal ref: AJ 168, 46 (2024)

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.

Comments: 25 pages, 13 figures. Accepted in AJ

Journal ref: AJ 167 243 (2024)

Searching for Giant Exoplanets around M-dwarf Stars (GEMS) I: Survey Motivation

Authors: Shubham Kanodia, Caleb I. Cañas, Suvrath Mahadevan, Eric B. Ford, Ravit Helled, Dana E. Anderson, Alan Boss, William D. Cochran, Megan Delamer, Te Han, Jessica E. Libby-Roberts, Andrea S. J. Lin, Simon Müller, Paul Robertson, Guðmundur Stefánsson, Johanna Teske

Abstract: Recent discoveries of transiting giant exoplanets around M-dwarf stars (GEMS), aided by the all-sky coverage of TESS, are starting to stretch theories of planet formation through the core-accretion scenario. Recent upper limits on their occurrence suggest that they decrease with lower stellar masses, with fewer GEMS around lower-mass stars compared to solar-type. In this paper, we discuss existing… ▽ More Recent discoveries of transiting giant exoplanets around M-dwarf stars (GEMS), aided by the all-sky coverage of TESS, are starting to stretch theories of planet formation through the core-accretion scenario. Recent upper limits on their occurrence suggest that they decrease with lower stellar masses, with fewer GEMS around lower-mass stars compared to solar-type. In this paper, we discuss existing GEMS both through confirmed planets, as well as protoplanetary disk observations, and a combination of tests to reconcile these with theoretical predictions. We then introduce the \textit{Searching for GEMS} survey, where we utilize multi-dimensional nonparameteric statistics to simulate hypothetical survey scenarios to predict the required sample size of transiting GEMS with mass measurements to robustly compare their bulk-density with canonical hot-Jupiters orbiting FGK stars. Our Monte-Carlo simulations predict that a robust comparison requires about 40 transiting GEMS (compared to the existing sample of $\sim$ 15) with 5-$σ$ mass measurements. Furthermore, we discuss the limitations of existing occurrence estimates for GEMS, and provide a brief description of our planned systematic search to improve the occurrence rate estimates for GEMS. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages + references, including 7 figures. Accepted in AAS Journals

Updated Catalog of Kepler Planet Candidates: Focus on Accuracy and Orbital Periods

Authors: Jack J. Lissauer, Jason F. Rowe, Daniel Jontof-Hutter, Daniel C. Fabrycky, Eric B. Ford, Darin Ragozzine, Jason H. Steffen, Kadri M. Nizam

Abstract: We present a new catalog of Kepler planet candidates that prioritizes accuracy of planetary dispositions and properties over uniformity. This catalog contains 4376 transiting planet candidates, including 1791 residing within 709 multi-planet systems, and provides the best parameters available for a large sample of Kepler planet candidates. We also provide a second set of stellar and planetary prop… ▽ More We present a new catalog of Kepler planet candidates that prioritizes accuracy of planetary dispositions and properties over uniformity. This catalog contains 4376 transiting planet candidates, including 1791 residing within 709 multi-planet systems, and provides the best parameters available for a large sample of Kepler planet candidates. We also provide a second set of stellar and planetary properties for transiting candidates that are uniformly-derived for use in occurrence rates studies. Estimates of orbital periods have been improved, but as in previous catalogs, our tabulated values for period uncertainties do not fully account for transit timing variations (TTVs). We show that many planets are likely to have TTVs with long periodicities caused by various processes, including orbital precession, and that such TTVs imply that ephemerides of Kepler planets are not as accurate on multi-decadal timescales as predicted by the small formal errors (typically 1 part in $10^6$ and rarely $ > 10^{-5}$) in the planets' measured mean orbital periods during the Kepler epoch. Analysis of normalized transit durations implies that eccentricities of planets are anti-correlated with the number of companion transiting planets. Our primary catalog lists all known Kepler planet candidates that orbit and transit only one star; for completeness, we also provide an abbreviated listing of the properties of the two dozen non-transiting planets that have been identified around stars that host transiting planets discovered by Kepler. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: 64 pages, 32 figures, 5 tables. Accepted for publication in The Planetary Science Journal. A high-resolution version of Figure 32 can be downloaded from the following link: https://kona.ubishops.ca/architecture/KeplerPlot_large_cbar.png

The Extreme Stellar-Signals Project III. Combining Solar Data from HARPS, HARPS-N, EXPRES, and NEID

Authors: Lily L. Zhao, Xavier Dumusque, Eric B. Ford, Joe Llama, Annelies Mortier, Megan Bedell, Khaled Al Moulla, Chad F. Bender, Cullen H. Blake, John M. Brewer, Andrew Collier Cameron, Rosario Cosentino, Pedro Figueira, Debra A. Fischer, Adriano Ghedina, Manuel Gonzalez, Samuel Halverson, Shubham Kanodia, David W. Latham, Andrea S. J. Lin, Gaspare Lo Curto, Marcello Lodi, Sarah E. Logsdon, Christophe Lovis, Suvrath Mahadevan , et al. (15 additional authors not shown)

Abstract: We present an analysis of Sun-as-a-star observations from four different high-resolution, stabilized spectrographs -- HARPS, HARPS-N, EXPRES, and NEID. With simultaneous observations of the Sun from four different instruments, we are able to gain insight into the radial velocity precision and accuracy delivered by each of these instruments and isolate instrumental systematics that differ from true… ▽ More We present an analysis of Sun-as-a-star observations from four different high-resolution, stabilized spectrographs -- HARPS, HARPS-N, EXPRES, and NEID. With simultaneous observations of the Sun from four different instruments, we are able to gain insight into the radial velocity precision and accuracy delivered by each of these instruments and isolate instrumental systematics that differ from true astrophysical signals. With solar observations, we can completely characterize the expected Doppler shift contributed by orbiting Solar System bodies and remove them. This results in a data set with measured velocity variations that purely trace flows on the solar surface. Direct comparisons of the radial velocities measured by each instrument show remarkable agreement with residual intra-day scatter of only 15-30 cm/s. This shows that current ultra-stabilized instruments have broken through to a new level of measurement precision that reveals stellar variability with high fidelity and detail. We end by discussing how radial velocities from different instruments can be combined to provide powerful leverage for testing techniques to mitigate stellar signals. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 17 pages, 9 figures, accepted for publication

Beyond 2-D Mass-Radius Relationships: A Nonparametric and Probabilistic Framework for Characterizing Planetary Samples in Higher Dimensions

Authors: Shubham Kanodia, Matthias Y. He, Eric B. Ford, Sujit K. Ghosh, Angie Wolfgang

Abstract: Fundamental to our understanding of planetary bulk compositions is the relationship between their masses and radii, two properties that are often not simultaneously known for most exoplanets. However, while many previous studies have modeled the two-dimensional relationship between planetary mass and radii, this approach largely ignores the dependencies on other properties that may have influenced… ▽ More Fundamental to our understanding of planetary bulk compositions is the relationship between their masses and radii, two properties that are often not simultaneously known for most exoplanets. However, while many previous studies have modeled the two-dimensional relationship between planetary mass and radii, this approach largely ignores the dependencies on other properties that may have influenced the formation and evolution of the planets. In this work, we extend the existing nonparametric and probabilistic framework of \texttt{MRExo} to jointly model distributions beyond two dimensions. Our updated framework can now simultaneously model up to four observables, while also incorporating asymmetric measurement uncertainties and upper limits in the data. We showcase the potential of this multi-dimensional approach to three science cases: (i) a 4-dimensional joint fit to planetary mass, radius, insolation, and stellar mass, hinting of changes in planetary bulk density across insolation and stellar mass; (ii) a 3-dimensional fit to the California Kepler Survey sample showing how the planet radius valley evolves across different stellar masses; and (iii) a 2-dimensional fit to a sample of Class-II protoplanetary disks in Lupus while incorporating the upper-limits in dust mass measurements. In addition, we employ bootstrap and Monte-Carlo sampling to quantify the impact of the finite sample size as well as measurement uncertainties on the predicted quantities. We update our existing open-source user-friendly \texttt{MRExo} \texttt{Python} package with these changes, which allows users to apply this highly flexible framework to a variety of datasets beyond what we have shown here. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: Accepted in ApJ. Updated MRExo package and sample scripts available here: https://github.com/shbhuk/mrexo/tree/v1.0dev. Package will be released on PyPI (pip) along with full documentation upon publication in ApJ

Statistical methods for exoplanet detection with radial velocities

Authors: Nathan C. Hara, Eric B. Ford

Abstract: Exoplanets can be detected with various observational techniques. Among them, radial velocity (RV) has the key advantages of revealing the architecture of planetary systems and measuring planetary mass and orbital eccentricities. RV observations are poised to play a key role in the detection and characterization of Earth twins. However, the detection of such small planets is not yet possible due t… ▽ More Exoplanets can be detected with various observational techniques. Among them, radial velocity (RV) has the key advantages of revealing the architecture of planetary systems and measuring planetary mass and orbital eccentricities. RV observations are poised to play a key role in the detection and characterization of Earth twins. However, the detection of such small planets is not yet possible due to very complex, temporally correlated instrumental and astrophysical stochastic signals. Furthermore, exploring the large parameter space of RV models exhaustively and efficiently presents difficulties. In this review, we frame RV data analysis as a problem of detection and parameter estimation in unevenly sampled, multivariate time series. The objective of this review is two-fold: to introduce the motivation, methodological challenges, and numerical challenges of RV data analysis to nonspecialists, and to unify the existing advanced approaches in order to identify areas for improvement. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: Published in Annual Reviews of Statistics and Its Application, contains the main text and supplemental appendix

A High-Eccentricity Warm Jupiter Orbiting TOI-4127

Authors: Arvind F. Gupta, Jonathan M. Jackson, Guillaume Hebrard, Andrea S. Lin, Keivan G. Stassun, Jiayin Dong, Steven Villanueva, Diana Dragomir, Suvrath Mahadevan, Jason T. Wright, Jose Manuel Almenara, Cullen H. Blake, Isabelle Boisse, Pia Cortes-Zuleta, Paul A. Dalba, Rodrigo F. Diaz, Eric B. Ford, Thierry Forveille, Robert Gagliano, Samuel P. Halverson, Neda Heidari, Shubham Kanodia, Flavien Kiefer, David W. Latham, Michael W. McElwain , et al. (14 additional authors not shown)

Abstract: We report the discovery of TOI-4127 b, a transiting, Jupiter-sized exoplanet on a long-period ($P = 56.39879^{+0.00010}_{-0.00010}$ d), high-eccentricity orbit around a late F-type dwarf star. This warm Jupiter was first detected and identified as a promising candidate from a search for single-transit signals in TESS Sector 20 data, and later characterized as a planet following two subsequent tran… ▽ More We report the discovery of TOI-4127 b, a transiting, Jupiter-sized exoplanet on a long-period ($P = 56.39879^{+0.00010}_{-0.00010}$ d), high-eccentricity orbit around a late F-type dwarf star. This warm Jupiter was first detected and identified as a promising candidate from a search for single-transit signals in TESS Sector 20 data, and later characterized as a planet following two subsequent transits (TESS Sectors 26 and 53) and follow-up ground-based RV observations with the NEID and SOPHIE spectrographs. We jointly fit the transit and RV data to constrain the physical ($R_p = 1.096^{+0.039}_{-0.032} R_J$, $M_p = 2.30^{+0.11}_{-0.11} M_J$) and orbital parameters of the exoplanet. Given its high orbital eccentricity ($e=0.7471^{+0.0078}_{-0.0086}$), TOI-4127 b is a compelling candidate for studies of warm Jupiter populations and of hot Jupiter formation pathways. We show that the present periastron separation of TOI-4127 b is too large for high-eccentricity tidal migration to circularize its orbit, and that TOI-4127 b is unlikely to be a hot Jupiter progenitor unless it is undergoing angular momentum exchange with an undetected outer companion. Although we find no evidence for an external companion, the available observational data are insufficient to rule out the presence of a perturber that can excite eccentricity oscillations and facilitate tidal migration. △ Less

Submitted 25 March, 2023; originally announced March 2023.

Detection of p-mode Oscillations in HD 35833 with NEID and TESS

Authors: Arvind F. Gupta, Jacob K. Luhn, Jason T. Wright, Suvrath Mahadevan, Eric B. Ford, Gudmundur Stefansson, Chad F. Bender, Cullen H. Blake, Samuel Halverson, Fred R. Hearty, Shubham Kanodia, Sarah E. Logsdon, Michael W. McElwain, Joe P. Ninan, Paul Robertson, Arpita Roy, Christian Schwab, Ryan C. Terrien

Abstract: We report the results of observations of p-mode oscillations in the G0 subgiant star HD 35833 in both radial velocities and photometry with NEID and TESS, respectively. We achieve separate, robust detections of the oscillation signal with both instruments (radial velocity amplitude $A_{\rm RV}=1.11\pm0.09$ m s$^{-1}$, photometric amplitude $A_{\rm phot}=6.42\pm0.60$ ppm, frequency of maximum power… ▽ More We report the results of observations of p-mode oscillations in the G0 subgiant star HD 35833 in both radial velocities and photometry with NEID and TESS, respectively. We achieve separate, robust detections of the oscillation signal with both instruments (radial velocity amplitude $A_{\rm RV}=1.11\pm0.09$ m s$^{-1}$, photometric amplitude $A_{\rm phot}=6.42\pm0.60$ ppm, frequency of maximum power $ν_{\rm max} = 595.71\pm17.28$ $μ$Hz, and mode spacing $Δν= 36.65\pm0.96$ $μ$Hz) as well as a non-detection in a TESS sector concurrent with the NEID observations. These data shed light on our ability to mitigate the correlated noise impact of oscillations with radial velocities alone, and on the robustness of commonly used asteroseismic scaling relations. The NEID data are used to validate models for the attenuation of oscillation signals for exposure times $t<ν_{\rm max}^{-1}$, and we compare our results to predictions from theoretical scaling relations and find that the observed amplitudes are weaker than expected by $>4σ$, hinting at gaps in the underlying physical models. △ Less

Submitted 2 October, 2022; originally announced October 2022.

Comments: 19 Pages, 14 Figures, Appendix

Debiasing the Minimum-Mass Extrasolar Nebula: On the Diversity of Solid Disk Profiles

Authors: Matthias Y. He, Eric B. Ford

Abstract: A foundational idea in the theory of in situ planet formation is the "minimum mass extrasolar nebula" (MMEN), a surface density profile ($Σ$) of disk solids that is necessary to form the planets in their present locations. While most previous studies have fit a single power-law to all exoplanets in an observed ensemble, it is unclear whether most exoplanetary systems form from a universal disk tem… ▽ More A foundational idea in the theory of in situ planet formation is the "minimum mass extrasolar nebula" (MMEN), a surface density profile ($Σ$) of disk solids that is necessary to form the planets in their present locations. While most previous studies have fit a single power-law to all exoplanets in an observed ensemble, it is unclear whether most exoplanetary systems form from a universal disk template. We use an advanced statistical model for the underlying architectures of multi-planet systems to reconstruct the MMEN. The simulated physical and Kepler-observed catalogs allows us to directly assess the role of detection biases, and in particular the effect of non-transiting or otherwise undetected planets, in altering the inferred MMEN. We find that fitting a power-law of the form $Σ= Σ_0^* (a/a_0)^β$ to each multi-planet system results in a broad distribution of disk profiles; $Σ_0^* = 336_{-291}^{+727}$ g/cm$^2$ and $β= -1.98_{-1.52}^{+1.55}$ encompass the 16th-84th percentiles of the marginal distributions in an underlying population, where $Σ_0^*$ is the normalization at $a_0 = 0.3$ AU. Around half of inner planet-forming disks have minimum solid masses of $\gtrsim 40 M_\oplus$ within 1 AU. While transit observations do not tend to bias the median $β$, they can lead to both significantly over- and under-estimated $Σ_0^*$ and thus broaden the inferred distribution of disk masses. Nevertheless, detection biases cannot account for the full variance in the observed disk profiles; there is no universal MMEN if all planets formed in situ. The great diversity of solid disk profiles suggests that a substantial fraction ($\gtrsim 23\%$) of planetary systems experienced a history of migration. △ Less

Submitted 7 October, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted to AJ. 14 pages, 6 figures, 1 table. Accompanying code is available via SysSimPyMMEN, a pip-installable Python package (see https://syssimpymmen.readthedocs.io/en/latest/)

GJ 3929: High Precision Photometric and Doppler Characterization of an Exo-Venus and its Hot, Mini-Neptune-mass Companion

Authors: Corey Beard, Paul Robertson, Shubham Kanodia, Jack Lubin, Caleb I. Cañas, Arvind F. Gupta, Rae Holcomb, Sinclaire Jones, Jessica E. Libby-Roberts, Andrea S. J. Lin, Suvrath Mahadevan, Guðmundur Stefánsson, Chad F. Bender, Cullen H. Blake, William D. Cochran, Michael Endl, Mark Everett, Eric B. Ford, Connor Fredrick, Samuel Halverson, Leslie Hebb, Dan Li, Sarah E. Logsdon, Jacob Luhn, Michael W. McElwain , et al. (9 additional authors not shown)

Abstract: We detail the follow up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b, (TOI-2013b) and its non-transiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star's Venus-zone (P$_{b}$ = 2.616272 $\pm$ 0.000005 days; S$_{b}$ = 17.3$^{+0.8}_{-0.7}$ S$_{\oplus}$) orbiting a nearby M dwarf. GJ 3929c is most likely a non-transiting sub-Nept… ▽ More We detail the follow up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b, (TOI-2013b) and its non-transiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star's Venus-zone (P$_{b}$ = 2.616272 $\pm$ 0.000005 days; S$_{b}$ = 17.3$^{+0.8}_{-0.7}$ S$_{\oplus}$) orbiting a nearby M dwarf. GJ 3929c is most likely a non-transiting sub-Neptune. Using the new, ultra-precise NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak National Observatory, we are able to modify the mass constraints of planet b reported in previous works and consequently improve the significance of the mass measurement to almost 4$σ$ confidence (M$_{b}$ = 1.75 $\pm$ 0.45 M$_{\oplus}$). We further adjust the orbital period of planet c from its alias at 14.30 $\pm$ 0.03 days to the likely true period of 15.04 $\pm$ 0.03 days, and adjust its minimum mass to m$\sin i$ = 5.71 $\pm$ 0.92 M$_{\oplus}$. Using the diffuser-assisted ARCTIC imager on the ARC 3.5 m telescope at Apache Point Observatory, in addition to publicly available TESS and LCOGT photometry, we are able to constrain the radius of planet b to R$_{p}$ = 1.09 $\pm$ 0.04 R$_{\oplus}$. GJ 3929b is a top candidate for transmission spectroscopy in its size regime (TSM = 14 $\pm$ 4), and future atmospheric studies of GJ 3929b stand to shed light on the nature of small planets orbiting M dwarfs. △ Less

Submitted 30 July, 2022; v1 submitted 21 July, 2022; originally announced July 2022.

Comments: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2204.09063

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.

Comments: 30 pages, 9 figures, accepted in AJ

Modeling Stellar Oscillations and Granulation in Radial Velocity Time Series: A Fourier-based Method

Authors: Zhao Guo, Eric B. Ford, Dennis Stello, Jacob K. Luhn, Suvrath Mahadevan, Arvind F. Gupta, Jie Yu

Abstract: Tens of thousands of solar-like oscillating stars have been observed by space missions. Their photometric variability in the Fourier domain can be parameterized by a sum of two super-Lorentizian functions for granulation and a Gaussian-shaped power excess for oscillation. The photometric granulation/oscillation parameters scale with stellar parameters and they can also make predictions for corresp… ▽ More Tens of thousands of solar-like oscillating stars have been observed by space missions. Their photometric variability in the Fourier domain can be parameterized by a sum of two super-Lorentizian functions for granulation and a Gaussian-shaped power excess for oscillation. The photometric granulation/oscillation parameters scale with stellar parameters and they can also make predictions for corresponding parameters in radial velocity measurements. Based on scaling relations, we simulate realistic radial velocity time series and examine how the root-mean-square scatter of radial velocity measurements varies with stellar parameters and different observation strategies such as the length of integration time and gaps in the time series. Using stars with extensive spectroscopic observations from the spectrographs (SONG and HARPS), we measure the granulation amplitude and timescale from the power spectrum of the radial velocity time series. We compare these measurements with literature values based on Kepler photometry. We find that the granulation amplitude in radial velocity can be well predicted from the photometry and scaling relations. Both granulation timescales in radial velocity agree with those predicted from photometry for giants and sub-giants. However, for main-sequence stars, only one granulation timescale in radial velocity is in agreement with the photometric-based values, while the other timescale generally lies at lower frequencies compared to the result of photometry. In conclusion, we show the photometric scaling relations from Kepler photometry and the scaling relationship to Doppler observations can be very useful for predicting the photometric and radial velocity stellar variabilities due to stellar granulation and oscillation. △ Less

Submitted 12 February, 2022; originally announced February 2022.

Comments: Submitted to ApJ, comments welcome

Migration traps as the root cause of the Kepler dichotomy

Authors: Brianna Zawadzki, Daniel Carrera, Eric B. Ford

Abstract: It is often assumed that the "Kepler dichotomy" -- the apparent excess of planetary systems with a single detected transiting planet in the Kepler catalog -- reflects an intrinsic bimodality in the mutual inclinations of planetary orbits. After conducting 600 simulations of planet formation followed by simulated Kepler observations, we instead propose that the apparent dichotomy reflects a diverge… ▽ More It is often assumed that the "Kepler dichotomy" -- the apparent excess of planetary systems with a single detected transiting planet in the Kepler catalog -- reflects an intrinsic bimodality in the mutual inclinations of planetary orbits. After conducting 600 simulations of planet formation followed by simulated Kepler observations, we instead propose that the apparent dichotomy reflects a divergence in the amount of migration and the separation of planetary semimajor axes into distinct "clusters". We find that our simulated high-mass systems migrate rapidly, bringing more planets into orbital periods of less than 200 days. The outer planets are often caught in a migration trap -- a range of planet masses and locations in which a dominant co-rotation torque prevents inward migration -- which splits the system into two clusters. If clusters are sufficiently separated, the inner cluster remains dynamically cold, leading to low mutual inclinations and a higher probability of detecting multiple transiting planets. Conversely, our simulated low-mass systems typically bring fewer planets inside 200 days, forming a single cluster that quickly becomes dynamically unstable, leading to collisions and high mutual inclinations. We propose an alternative explanation for the apparent Kepler dichotomy in which migration traps during formation lead to fewer planets inside the Kepler detection window, and where mutual inclinations play only a secondary role. If our scenario is correct, then Kepler's STIPs (Systems with Tightly-packed Inner Planets) are a sample of planets that escaped capture by co-rotation traps, and their sizes may be a valuable probe into the structure of protoplanetary discs. △ Less

Submitted 17 August, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: 21 pages, 10 figures. Accepted for publication in The Astrophysical Journal

NEID Rossiter-McLaughlin Measurement of TOI-1268b: A Young Warm Saturn Aligned with Its Cool Host Star

Authors: Jiayin Dong, Chelsea X. Huang, George Zhou, Rebekah I. Dawson, Gudmundur K. Stefánsson, Chad F. Bender, Cullen H. Blake, Eric B. Ford, Samuel Halverson, Shubham Kanodia, Suvrath Mahadevan, Michael W. McElwain, Joe P. Ninan, Paul Robertson, Arpita Roy, Christian Schwab, Daniel J. Stevens, Ryan C. Terrien, Andrew Vanderburg, Adam L. Kraus, Stephanie Douglas, Elisabeth Newton, Rayna Rampalli, Daniel M. Krolikowski, Karen A. Collins , et al. (34 additional authors not shown)

Abstract: Close-in gas giants present a surprising range of stellar obliquity, the angle between a planet's orbital axis and its host star's spin axis. It is unclear whether the obliquities reflect the planets' dynamical history (e.g., aligned for in situ formation or disk migration versus misaligned for high-eccentricity tidal migration) or whether other mechanisms (e.g., primordial misalignment or planet-… ▽ More Close-in gas giants present a surprising range of stellar obliquity, the angle between a planet's orbital axis and its host star's spin axis. It is unclear whether the obliquities reflect the planets' dynamical history (e.g., aligned for in situ formation or disk migration versus misaligned for high-eccentricity tidal migration) or whether other mechanisms (e.g., primordial misalignment or planet-star interactions) are more important in sculpting the obliquity distribution. Here we present the stellar obliquity measurement of TOI-1268 (TIC-142394656, $V_{\rm mag} {\sim} 10.9$), a young K-type dwarf hosting an 8.2-day period, Saturn-sized planet. TOI-1268's lithium abundance and rotation period suggest the system age between the ages of Pleiades cluster (${\sim}120$ Myr) and Praesepe cluster (${\sim}670$ Myr). Using the newly commissioned NEID spectrograph, we constrain the stellar obliquity of TOI-1268 via the Rossiter-McLaughlin (RM) effect from both radial velocity (RV) and Doppler Tomography (DT) signals. The 3$σ$ upper bounds of the projected stellar obliquity $|λ|$ from both models are below 60$^\circ$. The large host star separation ($a/R_\star {\sim} 17$), combined with the system's young age, makes it unlikely that the planet has realigned its host star. The stellar obliquity measurement of TOI-1268 probes the architecture of a young gas giant beyond the reach of tidal realignment ($a/R_\star {\gtrsim} 10$) and reveals an aligned or slightly misaligned system. △ Less

Submitted 30 January, 2022; originally announced January 2022.

Comments: 14 pages, 6 figures, 1 table, accepted for publication in ApJL; see independent work by Subjak et al. for RV follow-up of TOI-1268

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

FIESTA II. Disentangling stellar and instrumental variability from exoplanetary Doppler shifts in Fourier domain

Authors: Jinglin Zhao, Eric B. Ford, Chris G. Tinney

Abstract: The radial velocity (RV) detection of exoplanets is challenged by stellar spectroscopic variability that can mimic the presence of planets and by instrumental instability that can further obscure the detection. Both stellar and instrumental changes can distort the spectral line profiles and be misinterpreted as apparent RV shifts. We present an improved FourIEr \textit{phase} SpecTrum Analysis (… ▽ More The radial velocity (RV) detection of exoplanets is challenged by stellar spectroscopic variability that can mimic the presence of planets and by instrumental instability that can further obscure the detection. Both stellar and instrumental changes can distort the spectral line profiles and be misinterpreted as apparent RV shifts. We present an improved FourIEr \textit{phase} SpecTrum Analysis (FIESTA a.k.a. $\mathitφ$ESTA) to disentangle apparent velocity shifts due to a line deformation from a true Doppler shift. $\mathitφ$ESTA projects stellar spectrum's cross correlation function (CCF) onto a truncated set of Fourier basis functions. Using the amplitude and phase information from each Fourier mode, we can trace the line variability at different CCF width scales to robustly identify and mitigate multiple sources of RV contamination. For example, in our study of the 3 years of HARPS-N solar data, $\mathitφ$ESTA reveals the solar rotational effect, the long-term trend due to solar magnetic cycle, instrumental instability and apparent solar rotation rate changes. Applying a multiple linear regression model on $\mathitφ$ESTA metrics, we reduce the weighted rms noise from 1.89 m/s to 0.98 m/s. In addition, we observe a $\sim$3 days lag in the $\mathitφ$ESTA metrics, similar to the findings from previous studies on BIS and FWHM. △ Less

Submitted 18 July, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Observing the Sun as a star: Design and early results from the NEID solar feed

Authors: Andrea S. J. Lin, Andrew Monson, Suvrath Mahadevan, Joe P. Ninan, Samuel Halverson, Colin Nitroy, Chad F. Bender, Sarah E. Logsdon, Shubham Kanodia, Ryan C. Terrien, Arpita Roy, Jacob K. Luhn, Arvind F. Gupta, Eric B. Ford, Fred Hearty, Russ R. Laher, Emily Hunting, William R. McBride, Noah Isaac Salazar Rivera, Jayadev Rajagopal, Marsha J. Wolf, Paul Robertson, Jason T. Wright, Cullen H. Blake, Caleb I. Canas , et al. (5 additional authors not shown)

Abstract: Efforts with extreme-precision radial velocity (EPRV) instruments to detect small-amplitude planets are largely limited, on many timescales, by the effects of stellar variability and instrumental systematics. One avenue for investigating these effects is the use of small solar telescopes which direct disk-integrated sunlight to these EPRV instruments, observing the Sun at high cadence over months… ▽ More Efforts with extreme-precision radial velocity (EPRV) instruments to detect small-amplitude planets are largely limited, on many timescales, by the effects of stellar variability and instrumental systematics. One avenue for investigating these effects is the use of small solar telescopes which direct disk-integrated sunlight to these EPRV instruments, observing the Sun at high cadence over months or years. We have designed and built a solar feed system to carry out "Sun-as-a-star" observations with NEID, a very high precision Doppler spectrometer recently commissioned at the WIYN 3.5m Telescope at Kitt Peak National Observatory. The NEID solar feed has been taking observations nearly every day since December 2020; data is publicly available at the NASA Exoplanet Science Institute (NExScI) NEID Solar Archive: \url{https://neid.ipac.caltech.edu/search_solar.php}. In this paper, we present the design of the NEID solar feed and explanations behind our design intent. We also present early radial velocity (RV) results which demonstrate NEID's RV stability on the Sun over 4 months of commissioning: 0.66~m/s RMS under good sky conditions and improving to 0.41~m/s RMS under best conditions. △ Less

Submitted 15 February, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: Accepted to AJ

A hot Mars-sized exoplanet transiting an M dwarf

Authors: Caleb I. Cañas, Suvrath Mahadevan, William D. Cochran, Chad F. Bender, Eric D. Feigelson, C. E. Harman, Ravi Kumar Kopparapu, Gabriel A. Caceres, Scott A. Diddams, Michael Endl, Eric B. Ford, Samuel Halverson, Fred Hearty, Sinclaire Jones, Shubham Kanodia, Andrea S. J. Lin, Andrew J. Metcalf, Andrew Monson, Joe P. Ninan, Lawrence W. Ramsey, Paul Robertson, Arpita Roy, Christian Schwab, Guðmundur Stefánsson

Abstract: We validate the planetary nature of an ultra-short period planet orbiting the M dwarf KOI-4777. We use a combination of space-based photometry from Kepler, high-precision, near-infrared Doppler spectroscopy from the Habitable-zone Planet Finder, and adaptive optics imaging to characterize this system. KOI-4777.01 is a Mars-sized exoplanet ($\mathrm{R}_{p}=0.51 \pm 0.03R_{\oplus}$) orbiting the hos… ▽ More We validate the planetary nature of an ultra-short period planet orbiting the M dwarf KOI-4777. We use a combination of space-based photometry from Kepler, high-precision, near-infrared Doppler spectroscopy from the Habitable-zone Planet Finder, and adaptive optics imaging to characterize this system. KOI-4777.01 is a Mars-sized exoplanet ($\mathrm{R}_{p}=0.51 \pm 0.03R_{\oplus}$) orbiting the host star every 0.412-days ($\sim9.9$-hours). This is the smallest validated ultra-short period planet known and we see no evidence for additional massive companions using our HPF RVs. We constrain the upper $3σ$ mass to $M_{p}<0.34~\mathrm{M_\oplus}$ by assuming the planet is less dense than iron. Obtaining a mass measurement for KOI-4777.01 is beyond current instrumental capabilities. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 25 pages, 8 figures, accepted for publication in AJ

Journal ref: AJ, 163, 3 (2022)

The Warm Neptune GJ 3470b has a Polar Orbit

Authors: Gudmundur Stefansson, Suvrath Mahadevan, Cristobal Petrovich, Joshua N. Winn, Shubham Kanodia, Sarah C. Millholland, Marissa Maney, Caleb I. Cañas, John Wisniewski, Paul Robertson, Joe P. Ninan, Eric B. Ford, Chad F. Bender, Cullen H. Blake, Heather Cegla, William D. Cochran, Scott A. Diddams, Jiayin Dong, Michael Endl, Connor Fredrick, Samuel Halverson, Fred Hearty, Leslie Hebb, Teruyuki Hirano, Andrea S. J. Lin , et al. (12 additional authors not shown)

Abstract: The warm Neptune GJ 3470b transits a nearby ($d=29$pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5m Telescope at Kitt Peak Observatory, we model the classical Rossiter-Mclaughlin effect yielding a sky-projected obliquity of $λ=98_{-12}^{+15\:\circ}$ and a… ▽ More The warm Neptune GJ 3470b transits a nearby ($d=29$pc) bright slowly rotating M1.5-dwarf star. Using spectroscopic observations during two transits with the newly commissioned NEID spectrometer on the WIYN 3.5m Telescope at Kitt Peak Observatory, we model the classical Rossiter-Mclaughlin effect yielding a sky-projected obliquity of $λ=98_{-12}^{+15\:\circ}$ and a $v \sin i = 0.85_{-0.33}^{+0.27}$km/s. Leveraging information about the rotation period and size of the host star, our analysis yields a true obliquity of $ψ=95_{-8}^{+9\:\circ}$, revealing that GJ 3470b is on a polar orbit. Using radial velocities from HIRES, HARPS and the Habitable-zone Planet Finder, we show that the data are compatible with a long-term RV slope of $\dotγ = -0.0022 \pm 0.0011$m/s/day over a baseline of 12.9 years. If the RV slope is due to acceleration from another companion in the system, we show that such a companion is capable of explaining the polar and mildly eccentric orbit of GJ 3470b using two different secular excitation models. The existence of an outer companion can be further constrained with additional RV observations, Gaia astrometry, and future high-contrast imaging observations. Lastly, we show that tidal heating from GJ 3470b's mild eccentricity has most likely inflated the radius of GJ 3470b by a factor of $\sim$1.5-1.7, which could help account for its evaporating atmosphere. △ Less

Submitted 1 May, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 21 pages, 7 figures. Accepted in ApJL

$\texttt{GRASS}$: Distinguishing Planet-induced Doppler Signatures from Granulation with a Synthetic Spectra Generator

Authors: Michael L. Palumbo III, Eric B. Ford, Jason T. Wright, Suvrath Mahadevan, Alexander W. Wise, Johannes Löhner-Böttcher

Abstract: Owing to recent advances in radial-velocity instrumentation and observation techniques, the detection of Earth-mass planets around Sun-like stars may soon be primarily limited by intrinsic stellar variability. Several processes contribute to this variability, including starspots, pulsations, and granulation. Although many previous studies have focused on techniques to mitigate signals from pulsati… ▽ More Owing to recent advances in radial-velocity instrumentation and observation techniques, the detection of Earth-mass planets around Sun-like stars may soon be primarily limited by intrinsic stellar variability. Several processes contribute to this variability, including starspots, pulsations, and granulation. Although many previous studies have focused on techniques to mitigate signals from pulsations and other types of magnetic activity, granulation noise has to date only been partially addressed by empirically-motivated observation strategies and magnetohydrodynamic simulations. To address this deficit, we present the GRanulation And Spectrum Simulator ($\texttt{GRASS}$), a new tool designed to create time-series synthetic spectra with granulation-driven variability from spatially- and temporally-resolved observations of solar absorption lines. In this work, we present $\texttt{GRASS}$, detail its methodology, and validate its model against disk-integrated solar observations. As a first-of-its-kind empirical model for spectral variability due to granulation in a star with perfectly known center-of-mass radial-velocity behavior, $\texttt{GRASS}$ is an important tool for testing new methods of disentangling granular line-shape changes from true Doppler shifts. △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: 18 pages, 7 figures, 2 tables, accepted for publication in The Astronomical Journal

Extreme Precision Radial Velocity Working Group Final Report

Authors: Jonathan Crass, B. Scott Gaudi, Stephanie Leifer, Charles Beichman, Chad Bender, Gary Blackwood, Jennifer A. Burt, John L. Callas, Heather M. Cegla, Scott A. Diddams, Xavier Dumusque, Jason D. Eastman, Eric B. Ford, Benjamin Fulton, Rose Gibson, Samuel Halverson, Raphaëlle D. Haywood, Fred Hearty, Andrew W. Howard, David W. Latham, Johannes Löhner-Böttcher, Eric E. Mamajek, Annelies Mortier, Patrick Newman, Peter Plavchan , et al. (11 additional authors not shown)

Abstract: Precise mass measurements of exoplanets discovered by the direct imaging or transit technique are required to determine planet bulk properties and potential habitability. Furthermore, it is generally acknowledged that, for the foreseeable future, the Extreme Precision Radial Velocity (EPRV) measurement technique is the only method potentially capable of detecting and measuring the masses and orbit… ▽ More Precise mass measurements of exoplanets discovered by the direct imaging or transit technique are required to determine planet bulk properties and potential habitability. Furthermore, it is generally acknowledged that, for the foreseeable future, the Extreme Precision Radial Velocity (EPRV) measurement technique is the only method potentially capable of detecting and measuring the masses and orbits of habitable-zone Earths orbiting nearby F, G, and K spectral-type stars from the ground. In particular, EPRV measurements with a precision of better than approximately 10 cm/s (with a few cm/s stability over many years) are required. Unfortunately, for nearly a decade, PRV instruments and surveys have been unable to routinely reach RV accuracies of less than roughly 1 m/s. Making EPRV science and technology development a critical component of both NASA and NSF program plans is crucial for reaching the goal of detecting potentially habitable Earthlike planets and supporting potential future exoplanet direct imaging missions such as the Habitable Exoplanet Observatory (HabEx) or the Large Ultraviolet Optical Infrared Surveyor (LUVOIR). In recognition of these facts, the 2018 National Academy of Sciences (NAS) Exoplanet Science Strategy (ESS) report recommended the development of EPRV measurements as a critical step toward the detection and characterization of habitable, Earth-analog planets. In response to the NAS-ESS recommendation, NASA and NSF commissioned the EPRV Working Group to recommend a ground-based program architecture and implementation plan to achieve the goal intended by the NAS. This report documents the activities, findings, and recommendations of the EPRV Working Group. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: Full report: 103 pages. Executive summary: 7 pages. More information about the NASA-NSF Exoplanet Observational Research (NN-EXPLORE) program, including the NASA-NSF Extreme Precision Radial Velocity Initiative, can be found here: https://exoplanets.nasa.gov/exep/NNExplore/

Evidence for a Non-Dichotomous Solution to the Kepler Dichotomy: Mutual Inclinations of Kepler Planetary Systems from Transit Duration Variations

Authors: Sarah C. Millholland, Matthias Y. He, Eric B. Ford, Darin Ragozzine, Daniel Fabrycky, Joshua N. Winn

Abstract: Early analyses of exoplanet statistics from the Kepler Mission revealed that a model population of multiple-planet systems with low mutual inclinations (${\sim1^{\circ}-2^{\circ}}$) adequately describes the multiple-transiting systems but underpredicts the number of single-transiting systems. This so-called "Kepler dichotomy" signals the existence of a sub-population of multi-planet systems posses… ▽ More Early analyses of exoplanet statistics from the Kepler Mission revealed that a model population of multiple-planet systems with low mutual inclinations (${\sim1^{\circ}-2^{\circ}}$) adequately describes the multiple-transiting systems but underpredicts the number of single-transiting systems. This so-called "Kepler dichotomy" signals the existence of a sub-population of multi-planet systems possessing larger mutual inclinations. However, the details of these inclinations remain uncertain. In this work, we derive constraints on the intrinsic mutual inclination distribution by statistically exploiting Transit Duration Variations (TDVs) of the Kepler planet population. When planetary orbits are mutually inclined, planet-planet interactions cause orbital precession, which can lead to detectable long-term changes in transit durations. These TDV signals are inclination-sensitive and have been detected for roughly two dozen Kepler planets. We compare the properties of the Kepler observed TDV detections to TDV detections of simulated planetary systems constructed from two population models with differing assumptions about the mutual inclination distribution. We find strong evidence for a continuous distribution of relatively low mutual inclinations that is well-characterized by a power law relationship between the median mutual inclination ($\tildeμ_{i,n}$) and the intrinsic multiplicity ($n$): $\tildeμ_{i,n} = \tildeμ_{i,5}(n/5)^α$, where $\tildeμ_{i,5} = 1.10^{+0.15}_{-0.11}$ and $α= -1.73^{+0.09}_{-0.08}$. These results suggest that late-stage planet assembly and possibly stellar oblateness are the dominant physical origins for the excitation of Kepler planet mutual inclinations. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: Accepted to AJ, 19 pages, 10 figures (excluding appendix)

Friends and Foes: Conditional Occurrence Rates of Exoplanet Companions and their Impact on Radial Velocity Follow-up Surveys

Authors: Matthias Y. He, Eric B. Ford, Darin Ragozzine

Abstract: Population studies of Kepler's multi-planet systems have revealed a surprising degree of structure in their underlying architectures. Information from a detected transiting planet can be combined with a population model to make predictions about the presence and properties of additional planets in the system. Using a statistical model for the distribution of planetary systems (He et al. 2020; arXi… ▽ More Population studies of Kepler's multi-planet systems have revealed a surprising degree of structure in their underlying architectures. Information from a detected transiting planet can be combined with a population model to make predictions about the presence and properties of additional planets in the system. Using a statistical model for the distribution of planetary systems (He et al. 2020; arXiv:2007.14473), we compute the conditional occurrence of planets as a function of the period and radius of Kepler-detectable planets. About half ($0.52 \pm 0.03$) of the time, the detected planet is not the planet with the largest semi-amplitude $K$ in the system, so efforts to measure the mass of the transiting planet with radial velocity (RV) follow-up will have to contend with additional planetary signals in the data. We simulate RV observations to show that assuming a single-planet model to measure the $K$ of the transiting planet often requires significantly more observations than in the ideal case with no additional planets, due to the systematic errors from unseen planet companions. Our results show that planets around 10-day periods with $K$ close to the single-measurement RV precision ($σ_{1,\rm obs}$) typically require $\sim 100$ observations to measure their $K$ to within 20% error. For a next generation RV instrument achieving $σ_{1,\rm obs} = 10$ cm/s, about $\sim 200$ ($600$) observations are needed to measure the $K$ of a transiting Venus in a Kepler-like system to better than 20% (10%) error, which is $\sim 2.3$ times as many as what would be necessary for a Venus without any planetary companions. △ Less

Submitted 6 October, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: Accepted to AJ. 28 pages, 13 figures

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.

Comments: 20 pages, 12 figures, accepted for publication in The Astronomical Journal

A Harsh Test of Far-Field Scrambling with the Habitable Zone Planet Finder and the Hobby Eberly Telescope

Authors: Shubham Kanodia, Samuel Halverson, Joe P. Ninan, Suvrath Mahadevan, Gudmundur Stefansson, Arpita Roy, Lawrence W. Ramsey, Chad F. Bender, Steven Janowiecki, William D. Cochran, Scott A. Diddams, Niv Drory, Michael Endl, Eric B. Ford, Fred Hearty, Andrew J. Metcalf, Andrew Monson, Paul Robertson, Christian Schwab, Ryan C. Terrien, Jason T. Wright

Abstract: The Habitable zone Planet Finder (HPF) is a fiber fed precise radial velocity spectrograph at the 10 m Hobby Eberly Telescope (HET). Due to its fixed altitude design, the HET pupil changes appreciably across a track, leading to significant changes of the fiber far-field illumination. HPF's fiber scrambler is designed to suppress the impact of these illumination changes on the radial velocities --… ▽ More The Habitable zone Planet Finder (HPF) is a fiber fed precise radial velocity spectrograph at the 10 m Hobby Eberly Telescope (HET). Due to its fixed altitude design, the HET pupil changes appreciably across a track, leading to significant changes of the fiber far-field illumination. HPF's fiber scrambler is designed to suppress the impact of these illumination changes on the radial velocities -- but the residual impact on the radial velocity measurements has yet to be probed on sky. We use GJ 411, a bright early type (M2) M dwarf to probe the effects of far-field input trends due to these pupil variations on HPF radial velocities (RVs). These large changes ($\sim$ 2x) in pupil area and centroid present a harsh test of HPF's far-field scrambling. Our results show that the RVs are effectively decoupled from these extreme far-field input changes due to pupil centroid offsets, attesting to the effectiveness of the scrambler design. This experiment allows us to test the impact of these changes with large pupil variation on-sky, something we would not easily be able to do at a conventional optical telescope. While the pupil and illumination changes expected at these other telescopes are small, scaling from our results enables us to estimate and bound these effects, and show that they are controllable even for the new and next generation of RV instruments in their quest to beat down instrumental noise sources towards the goal of a few cm/s. △ Less

Submitted 5 May, 2021; v1 submitted 8 March, 2021; originally announced March 2021.

Comments: 16 pages. Published in The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 912, Number 1, 2021

Target Prioritization and Observing Strategies for the NEID Earth Twin Survey

Authors: Arvind F. Gupta, Jason T. Wright, Paul Robertson, Samuel Halverson, Jacob Luhn, Arpita Roy, Suvrath Mahadevan, Eric B . Ford, Chad F. Bender, Cullen H. Blake, Fred Hearty, Shubham Kanodia, Sarah E. Logsdon, Michael W. McElwain, Andrew Monson, Joe P. Ninan, Christian Schwab, Gudmundur Stefansson, Ryan C. Terrien

Abstract: NEID is a high-resolution optical spectrograph on the WIYN 3.5-m telescope at Kitt Peak National Observatory and will soon join the new generation of extreme precision radial velocity instruments in operation around the world. We plan to use the instrument to conduct the NEID Earth Twin Survey (NETS) over the course of the next 5 years, collecting hundreds of observations of some of the nearest an… ▽ More NEID is a high-resolution optical spectrograph on the WIYN 3.5-m telescope at Kitt Peak National Observatory and will soon join the new generation of extreme precision radial velocity instruments in operation around the world. We plan to use the instrument to conduct the NEID Earth Twin Survey (NETS) over the course of the next 5 years, collecting hundreds of observations of some of the nearest and brightest stars in an effort to probe the regime of Earth-mass exoplanets. Even if we take advantage of the extreme instrumental precision conferred by NEID, it will remain difficult to disentangle the weak (~10 cm s$^{-1}$) signals induced by such low-mass, long-period exoplanets from stellar noise for all but the quietest host stars. In this work, we present a set of quantitative selection metrics which we use to identify an initial NETS target list consisting of stars conducive to the detection of exoplanets in the regime of interest. We also outline a set of observing strategies with which we aim to mitigate uncertainty contributions from intrinsic stellar variability and other sources of noise. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: 19 pages, 10 figures, 1 table

Following up the Kepler field: Masses of Targets for transit timing and atmospheric characterization

Authors: Daniel Jontof-Hutter, Angie Wolfgang, Eric B. Ford, Jack J. Lissauer, Daniel C. Fabrycky, Jason F. Rowe

Abstract: We identify a set of planetary systems observed by Kepler that merit transit timing variation (TTV) analysis given the orbital periods of transiting planets, the uncertainties for their transit times and the number of transits observed during the Kepler mission. We confirm the planetary nature of 4 KOIs within multicandidate systems. We forward model each of the planetary systems identified to det… ▽ More We identify a set of planetary systems observed by Kepler that merit transit timing variation (TTV) analysis given the orbital periods of transiting planets, the uncertainties for their transit times and the number of transits observed during the Kepler mission. We confirm the planetary nature of 4 KOIs within multicandidate systems. We forward model each of the planetary systems identified to determine which systems are likely to yield mass constraints that may be significantly improved upon with follow-up transit observations. We find projected TTVs diverge by more than 90 minutes after 6000 days in 27 systems, including 22 planets with orbital periods exceeding 25 days. Such targets would benefit the most from additional transit timing data. TTV follow-up could push exoplanet characterization to lower masses, at greater orbital periods and at cooler equilibrium temperatures than is currently possible from the Kepler dataset alone. Combining TTVs and recently revised stellar parameters, we characterize an ensemble of homogeneously selected planets and identify planets in the Kepler field with large enough estimated transmission annuli for atmospheric characterization with JWST. △ Less

Submitted 27 May, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 70 pages, 33 figures

Journal ref: AJ,161,246J (2021)

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

The Occurrence of Rocky Habitable Zone Planets Around Solar-Like Stars from Kepler Data

Authors: Steve Bryson, Michelle Kunimoto, Ravi K. Kopparapu, Jeffrey L. Coughlin, William J. Borucki, David Koch, Victor Silva Aguirre, Christopher Allen, Geert Barentsen, Natalie. M. Batalha, Travis Berger, Alan Boss, Lars A. Buchhave, Christopher J. Burke, Douglas A. Caldwell, Jennifer R. Campbell, Joseph Catanzarite, Hema Chandrasekharan, William J. Chaplin, Jessie L. Christiansen, Jorgen Christensen-Dalsgaard, David R. Ciardi, Bruce D. Clarke, William D. Cochran, Jessie L. Dotson , et al. (57 additional authors not shown)

Abstract: We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orb… ▽ More We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orbiting stars with effective temperatures between 4800 K and 6300 K. We find that $η_\oplus$ for the conservative HZ is between $0.37^{+0.48}_{-0.21}$ (errors reflect 68\% credible intervals) and $0.60^{+0.90}_{-0.36}$ planets per star, while the optimistic HZ occurrence is between $0.58^{+0.73}_{-0.33}$ and $0.88^{+1.28}_{-0.51}$ planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates using both a Poisson likelihood Bayesian analysis and Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with $95\%$ confidence that, on average, the nearest HZ planet around G and K dwarfs is about 6 pc away, and there are about 4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun. △ Less

Submitted 3 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: To appear in The Astronomical Journal

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.

Comments: 19 pages, 4 figures, published in ApJ

Journal ref: ApJ 905 155 (2020)

Architectures of Exoplanetary Systems. III: Eccentricity and Mutual Inclination Distributions of AMD-stable Planetary Systems

Authors: Matthias Y. He, Eric B. Ford, Darin Ragozzine, Daniel Carrera

Abstract: The angular momentum deficit (AMD) of a planetary system is a measure of its orbital excitation and a predictor of long-term stability. We adopt the AMD-stability criteria to constrain the orbital architectures for exoplanetary systems. Previously, He, Ford, & Ragozzine (2019) (arXiv:1907.07773v2) showed through forward modelling (SysSim) that the observed multiplicity distribution can be well rep… ▽ More The angular momentum deficit (AMD) of a planetary system is a measure of its orbital excitation and a predictor of long-term stability. We adopt the AMD-stability criteria to constrain the orbital architectures for exoplanetary systems. Previously, He, Ford, & Ragozzine (2019) (arXiv:1907.07773v2) showed through forward modelling (SysSim) that the observed multiplicity distribution can be well reproduced by two populations consisting of a low and a high mutual inclination component. Here, we show that a broad distribution of mutual inclinations arising from systems at the AMD-stability limit can also match the observed Kepler population. We show that distributing a planetary system's maximum AMD amongst its planets results in a multiplicity-dependent distribution of eccentricities and mutual inclinations. Systems with intrinsically more planets have lower median eccentricities and mutual inclinations, and this trend is well described by power-law functions of the intrinsic planet multiplicity ($n$): $\tildeμ_{e,n} \propto n^{-1.74_{-0.07}^{+0.11}}$ and $\tildeμ_{i,n} \propto n^{-1.73_{-0.08}^{+0.09}}$, where $\tildeμ_{e,n}$ and $\tildeμ_{i,n}$ are the medians of the eccentricity and inclination distributions. We also find that intrinsic single planets have higher eccentricities ($σ_{e,1} \sim 0.25$) than multi-planet systems, and that the trends with multiplicity appear in the observed distributions of period-normalized transit duration ratios. We show that the observed preferences for planet size orderings and uniform spacings are more extreme than what can be produced by the detection biases of the Kepler mission alone. Finally, we find that for systems with detected transiting planets between 5 and 10 days, there is another planet with a greater radial velocity signal $\simeq~53\%$ of the time. △ Less

Submitted 9 December, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: Published in AJ; 41 pages, 15 figures, plus supplemental online materials with 8 figures

Journal ref: AJ, 160, 276 (2020)

The Habitable-zone Planet Finder Reveals A High Mass and a Low Obliquity for the Young Neptune K2-25b

Authors: Gudmundur Stefansson, Suvrath Mahadevan, Marissa Maney, Joe P. Ninan, Paul Robertson, Jayadev Rajagopal, Flynn Haase, Lori Allen, Eric B. Ford, Joshua Winn, Angie Wolfgang, Rebekah I. Dawson, John Wisniewski, Chad F. Bender, Caleb Cañas, William Cochran, Scott A. Diddams, Connor Fredrick, Samuel Halverson, Fred Hearty, Leslie Hebb, Shubham Kanodia, Eric Levi, Andrew J. Metcalf, Andrew Monson , et al. (5 additional authors not shown)

Abstract: Using radial-velocity data from the Habitable-zone Planet Finder, we have measured the mass of the Neptune-sized planet K2-25b, as well as the obliquity of its M4.5-dwarf host star in the 600-800MYr Hyades cluster. This is one of the youngest planetary systems for which both of these quantities have been measured, and one of the very few M dwarfs with a measured obliquity. Based on a joint analysi… ▽ More Using radial-velocity data from the Habitable-zone Planet Finder, we have measured the mass of the Neptune-sized planet K2-25b, as well as the obliquity of its M4.5-dwarf host star in the 600-800MYr Hyades cluster. This is one of the youngest planetary systems for which both of these quantities have been measured, and one of the very few M dwarfs with a measured obliquity. Based on a joint analysis of the radial velocity data, time-series photometry from the K2 mission, and new transit light curves obtained with diffuser-assisted photometry, the planet's radius and mass are $3.44\pm 0.12 \mathrm{R_\oplus}$ and $24.5_{-5.2}^{+5.7} \mathrm{M_\oplus}$. These properties are compatible with a rocky core enshrouded by a thin hydrogen-helium atmosphere (5% by mass). We measure an orbital eccentricity of $e=0.43 \pm 0.05$. The sky-projected stellar obliquity is $λ=3 \pm 16^{\circ}$, compatible with spin-orbit alignment, in contrast to other "hot Neptunes" that have been studied around older stars. △ Less

Submitted 24 July, 2020; originally announced July 2020.

Comments: Accepted for publication in AJ, 31 pages, 14 figures

A warm Jupiter transiting an M dwarf: A TESS single transit event confirmed with the Habitable-zone Planet Finder

Authors: Caleb I. Cañas, Gudmundur Stefansson, Shubham Kanodia, Suvrath Mahadevan, William D. Cochran, Michael Endl, Paul Robertson, Chad F. Bender, Joe P. Ninan, Corey Beard, Jack Lubin, Arvind F. Gupta, Mark E. Everett, Andrew Monson, Robert F. Wilson, Hannah M. Lewis, Mary Brewer, Steven R. Majewski, Leslie Hebb, Rebekah I. Dawson, Scott A. Diddams, Eric B. Ford, Connor Fredrick, Samuel Halverson, Fred Hearty , et al. (8 additional authors not shown)

Abstract: We confirm the planetary nature of a warm Jupiter transiting the early M dwarf TOI-1899, using a combination of available TESS photometry; high-precision, near-infrared spectroscopy with the Habitable-zone Planet Finder; and speckle and adaptive optics imaging. The data reveal a transiting companion on an $\sim29$-day orbit with a mass and radius of $0.66\pm0.07\ \mathrm{M_{J}}$ and… ▽ More We confirm the planetary nature of a warm Jupiter transiting the early M dwarf TOI-1899, using a combination of available TESS photometry; high-precision, near-infrared spectroscopy with the Habitable-zone Planet Finder; and speckle and adaptive optics imaging. The data reveal a transiting companion on an $\sim29$-day orbit with a mass and radius of $0.66\pm0.07\ \mathrm{M_{J}}$ and $1.15_{-0.05}^{+0.04}\ \mathrm{R_{J}}$, respectively. The star TOI-1899 is the lowest-mass star known to host a transiting warm Jupiter, and we discuss the follow-up opportunities afforded by a warm ($\mathrm{T_{eq}}\sim362$ K) gas giant orbiting an M0 star. Our observations reveal that TOI-1899.01 is a puffy warm Jupiter, and we suggest additional transit observations to both refine the orbit and constrain the true dilution observed in TESS. △ Less

Submitted 3 September, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: 24 pages, 5 figures, 3 tables, published in AJ

TOI-1728b: The Habitable-zone Planet Finder confirms a warm super Neptune orbiting an M dwarf host

Authors: Shubham Kanodia, Caleb I. Canas, Gudmundur Stefansson, Joe P. Ninan, Leslie Hebb, Andrea S. J. Lin, Helen Baran, Marissa Maney, Ryan C. Terrien, 7 Suvrath Mahadevan, William D. Cochran, Michael Endl, Jiayin Dong, Chad F. Bender, Scott A. Diddams, Eric B. Ford, Connor Fredrick, Samuel Halverson, Fred Hearty, Andrew J. Metcalf, Andrew Monson, Lawrence W. Ramsey, Paul Robertson, Arpita Roy, Christian Schwab , et al. (1 additional authors not shown)

Abstract: We confirm the planetary nature of TOI-1728b using a combination of ground-based photometry, near-infrared Doppler velocimetry and spectroscopy with the Habitable-zone Planet Finder.TOI-1728 is an old, inactive M0 star with \teff{} $= 3980^{+31}_{-32}$ K, which hosts a transiting super Neptune at an orbital period of $\sim$ 3.49 days. Joint fitting of the radial velocities and TESS and ground-base… ▽ More We confirm the planetary nature of TOI-1728b using a combination of ground-based photometry, near-infrared Doppler velocimetry and spectroscopy with the Habitable-zone Planet Finder.TOI-1728 is an old, inactive M0 star with \teff{} $= 3980^{+31}_{-32}$ K, which hosts a transiting super Neptune at an orbital period of $\sim$ 3.49 days. Joint fitting of the radial velocities and TESS and ground-based transits yields a planetary radius of $5.05_{-0.17}^{+0.16}$ R$_{\oplus}$, mass $26.78_{-5.13}^{+5.43}$ M$_{\oplus}$ and eccentricity $0.057_{-0.039}^{+0.054}$. We estimate the stellar properties, and perform a search for He 10830 Åabsorption during the transit of this planet and claim a null detection with an upper limit of 1.1$\%$ with 90\% confidence. A deeper level of He 10830 Å~ absorption has been detected in the planet atmosphere of GJ 3470b, a comparable gaseous planet. TOI-1728b is the largest super Neptune -- the intermediate subclass of planets between Neptune and the more massive gas-giant planets -- discovered around an M dwarf. With its relatively large mass and radius, TOI-1728 represents a valuable datapoint in the M-dwarf exoplanet mass-radius diagram, bridging the gap between the lighter Neptune-sized planets and the heavier Jovian planets known to orbit M-dwarfs. With a low bulk density of $1.14_{-0.24}^{+0.26}$ g/cm$^3$, and orbiting a bright host star (J $\sim 9.6$, V $\sim 12.4$), TOI-1728b is also a promising candidate for transmission spectroscopy both from the ground and from space, which can be used to constrain planet formation and evolutionary models. △ Less

Submitted 25 June, 2020; originally announced June 2020.

Comments: 21 pages, 12 figures, 4 tables: Accepted for publication

Persistent starspot signals on M dwarfs: multi-wavelength Doppler observations with the Habitable-zone Planet Finder and Keck/HIRES

Authors: Paul Robertson, Gudmundur Stefansson, Suvrath Mahadevan, Michael Endl, William D. Cochran, Corey Beard, Chad F. Bender, Scott A. Diddams, Nicholas Duong, Eric B. Ford, Connor Fredrick, Samuel Halverson, Fred Hearty, Rae Holcomb, Lydia Juan, Shubham Kanodia, Jack Lubin, Andrew J. Metcalf, Andrew Monson, Joe P. Ninan, Jonathan Palafoutas, Lawrence W. Ramsey, Arpita Roy, Christian Schwab, Ryan C. Terrien , et al. (1 additional authors not shown)

Abstract: Young, rapidly-rotating M dwarfs exhibit prominent starspots, which create quasiperiodic signals in their photometric and Doppler spectroscopic measurements. The periodic Doppler signals can mimic radial velocity (RV) changes expected from orbiting exoplanets. Exoplanets can be distinguished from activity-induced false positives by the chromaticity and long-term incoherence of starspot signals, bu… ▽ More Young, rapidly-rotating M dwarfs exhibit prominent starspots, which create quasiperiodic signals in their photometric and Doppler spectroscopic measurements. The periodic Doppler signals can mimic radial velocity (RV) changes expected from orbiting exoplanets. Exoplanets can be distinguished from activity-induced false positives by the chromaticity and long-term incoherence of starspot signals, but these qualities are poorly constrained for fully-convective M stars. Coherent photometric starspot signals on M dwarfs may persist for hundreds of rotations, and the wavelength dependence of starspot RV signals may not be consistent between stars due to differences in their magnetic fields and active regions. We obtained precise multi-wavelength RVs of four rapidly-rotating M dwarfs (AD Leo, G 227-22, GJ 1245B, GJ 3959) using the near-infrared (NIR) Habitable-zone Planet Finder, and the optical Keck/HIRES spectrometer. Our RVs are complemented by photometry from Kepler, TESS, and the Las Cumbres Observatory (LCO) network of telescopes. We found that all four stars exhibit large spot-induced Doppler signals at their rotation periods, and investigated the longevity and optical-to-NIR chromaticity for these signals. The phase curves remain coherent much longer than is typical for Sunlike stars. Their chromaticity varies, and one star (GJ 3959) exhibits optical and NIR RV modulation consistent in both phase and amplitude. In general, though, we find that the NIR amplitudes are lower than their optical counterparts. We conclude that starspot modulation for rapidly-rotating M stars frequently remains coherent for hundreds of stellar rotations, and gives rise to Doppler signals that, due to this coherence, may be mistaken for exoplanets. △ Less

Submitted 19 May, 2020; originally announced May 2020.

Comments: Accepted for publication in the Astrophysical Journal

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.

Comments: 3 pages, 1 figure, published in RNAAS

Journal ref: RNAAS volume 4, number 4, April 2020, page 59

Architectures of Exoplanetary Systems. II: An Increase in Inner Planetary System Occurrence Toward Later Spectral Types for Kepler's FGK Dwarfs

Authors: Matthias Y. He, Eric B. Ford, Darin Ragozzine

Abstract: The Kepler mission observed thousands of transiting exoplanet candidates around hundreds of thousands of FGK dwarf stars. He, Ford, & Ragozzine (2019) (arXiv:1907.07773) applied forward modeling to infer the distribution of intrinsic architectures of planetary systems, developed a clustered Poisson point process model for exoplanetary systems (SysSim) to reproduce the marginal distributions of the… ▽ More The Kepler mission observed thousands of transiting exoplanet candidates around hundreds of thousands of FGK dwarf stars. He, Ford, & Ragozzine (2019) (arXiv:1907.07773) applied forward modeling to infer the distribution of intrinsic architectures of planetary systems, developed a clustered Poisson point process model for exoplanetary systems (SysSim) to reproduce the marginal distributions of the observed Kepler population, and they showed that orbital periods and planet radii are clustered within a given planetary system. Here, we extend the clustered model to explore correlations between planetary systems and their host star properties. We split the sample of Kepler FGK dwarfs into two halves and model the fraction of stars with planets (between 0.5-10 $R_\oplus$ and 3--300 days), $f_{\rm swpa}$, as a linear function of the Gaia DR2 color. We confirm previous findings that the occurrence of these planetary systems rises significantly toward later type (redder) stars. The fraction of stars with planets increases from $f_{\rm swpa} = 0.32_{-0.11}^{+0.12}$ for F2V dwarfs to $f_{\rm swpa} = 0.96_{-0.19}^{+0.04}$ for mid K-dwarfs. About half ($f_{\rm swpa} = 0.57_{-0.10}^{+0.14}$) of all solar-type (G2V) dwarfs harbor a planetary system between 3 and 300 days. This simple model can closely match the observed multiplicity distributions of both the bluer and redder halves in our sample, suggesting that the architectures of planetary systems around stars of different spectral types may be similar aside from a shift in the overall fraction of planet hosting stars. △ Less

Submitted 9 December, 2020; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: Published in AJ; 28 pages, 9 figures, plus supplemental online materials with 5 figures

Journal ref: AJ, 161, 16 (2021)

Occurrence Rates of Planets Orbiting M Stars: Applying ABC to Kepler DR25, Gaia DR2, and 2MASS Data

Authors: Danley C. Hsu, Eric B. Ford, Ryan Terrien

Abstract: We present robust planet occurrence rates for Kepler planet candidates around M stars for planet radii $R_p = 0.5-4~\textrm{R}_\oplus$ and orbital periods $P = 0.5-256$ days using the approximate Bayesian computation (ABC) technique. This work incorporates the final Kepler DR25 planet candidate catalog and data products and augments them with updated stellar properties using Gaia DR2 and 2MASS PSC… ▽ More We present robust planet occurrence rates for Kepler planet candidates around M stars for planet radii $R_p = 0.5-4~\textrm{R}_\oplus$ and orbital periods $P = 0.5-256$ days using the approximate Bayesian computation (ABC) technique. This work incorporates the final Kepler DR25 planet candidate catalog and data products and augments them with updated stellar properties using Gaia DR2 and 2MASS PSC. We apply a set of selection criteria to select a sample of 1,746 Kepler M dwarf targets that host 89 associated planet candidates. These early M dwarfs and late K dwarfs were selected from cross-referenced targets using several photometric quality flags from Gaia DR2 and color-magnitude cuts using 2MASS magnitudes. We estimate a habitable zone occurrence rate of $f_{\textrm{M,HZ}} = 0.33^{+0.10}_{-0.12}$ for planets with $0.75-1.5$ R$_\oplus$ size. We caution that occurrence rate estimates for Kepler M stars are sensitive to the choice of prior due to the small sample of target stars and planet candidates. For example, we find an occurrence rate of $4.2^{+0.6}_{-0.6}$ or $8.4^{+1.2}_{-1.1}$ planets per M dwarf (integrating over $R_p = 0.5-4~\textrm{R}_\oplus$ and $P = 0.5-256$ days) for our two choices of prior. These occurrence rates are greater than those for FGK dwarf target when compared at the same range of orbital periods, but similar to occurrence rates when computed as a function of equivalent stellar insolation. Combining our result with recent studies of exoplanet architectures indicates that most, and potentially all, early-M dwarfs harbor planetary systems. △ Less

Submitted 2 August, 2020; v1 submitted 6 February, 2020; originally announced February 2020.

Comments: 15 pages, 4 figures; submitted to MNRAS

Superabundance of Exoplanet Sub-Neptunes Explained by Fugacity Crisis

Authors: Edwin S. Kite, Bruce Fegley Jr., Laura Schaefer, Eric B. Ford

Abstract: Transiting planets with radii 2-3 $R_\bigoplus$ are much more numerous than larger planets. We propose that this drop-off is so abrupt because at $R$ $\sim$ 3 $R_\bigoplus$, base-of-atmosphere pressure is high enough for the atmosphere to readily dissolve into magma, and this sequestration acts as a strong brake on further growth. The viability of this idea is demonstrated using a simple model. Ou… ▽ More Transiting planets with radii 2-3 $R_\bigoplus$ are much more numerous than larger planets. We propose that this drop-off is so abrupt because at $R$ $\sim$ 3 $R_\bigoplus$, base-of-atmosphere pressure is high enough for the atmosphere to readily dissolve into magma, and this sequestration acts as a strong brake on further growth. The viability of this idea is demonstrated using a simple model. Our results support extensive magma-atmosphere equilibration on sub-Neptunes, with numerous implications for sub-Neptune formation and atmospheric chemistry. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: Accepted by Astrophysical Journal Letters

A sub-Neptune sized planet transiting the M2.5-dwarf G 9-40: Validation with the Habitable-zone Planet Finder

Authors: Gudmundur Stefansson, Caleb Cañas, John Wisniewski, Paul Robertson, Suvrath Mahadevan, Marissa Maney, Shubham Kanodia, Corey Beard, Chad F. Bender, Peter Brunt, J. Christopher Clemens, William Cochran, Scott A. Diddams, Michael Endl, Eric B. Ford, Connor Fredrick, Samuel Halverson, Fred Hearty, Leslie Hebb, Joseph Huehnerhoff, Jeff Jennings, Kyle Kaplan, Eric Levi, Emily Lubar, Andrew J. Metcalf , et al. (10 additional authors not shown)

Abstract: We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distan… ▽ More We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distance of $d=27.9\mathrm{pc}$, G 9-40b is the second closest transiting planet discovered by K2 to date. The planet's large transit depth ($\sim$3500ppm), combined with the proximity and brightness of the host star at NIR wavelengths (J=10, K=9.2) makes G 9-40b one of the most favorable sub-Neptune-sized planet orbiting an M-dwarf for transmission spectroscopy with JWST, ARIEL, and the upcoming Extremely Large Telescopes. The star is relatively inactive with a rotation period of $\sim$29 days determined from the K2 photometry. To estimate spectroscopic stellar parameters, we describe our implementation of an empirical spectral matching algorithm using the high-resolution NIR HPF spectra. Using this algorithm, we obtain an effective temperature of $T_{\mathrm{eff}}=3404\pm73$K, and metallicity of $\mathrm{[Fe/H]}=-0.08\pm0.13$. Our RVs, when coupled with the orbital parameters derived from the transit photometry, exclude planet masses above $11.7 M_\oplus$ with 99.7% confidence assuming a circular orbit. From its radius, we predict a mass of $M=5.0^{+3.8}_{-1.9} M_\oplus$ and an RV semi-amplitude of $K=4.1^{+3.1}_{-1.6}\mathrm{m\:s^{-1}}$, making its mass measurable with current RV facilities. We urge further RV follow-up observations to precisely measure its mass, to enable precise transmission spectroscopic measurements in the future. △ Less

Submitted 30 November, 2019; originally announced December 2019.

Comments: Accepted for publication in AJ, 22 pages, 15 figures

Evidence for He I 10830 Å~ absorption during the transit of a warm Neptune around the M-dwarf GJ 3470 with the Habitable-zone Planet Finder

Authors: Joe P. Ninan, Gudmundur Stefansson, Suvrath Mahadevan, Chad Bender, Paul Robertson, Lawrence Ramsey, Ryan Terrien, Jason Wright, Scott A. Diddams, Shubham Kanodia, William Cochran, Michael Endl, Eric B. Ford, Connor Fredrick, Samuel Halverson, Fred Hearty, Jeff Jennings, Kyle Kaplan, Emily Lubar, Andrew J. Metcalf, Andrew Monson, Colin Nitroy, Arpita Roy, Christian Schwab

Abstract: Understanding the dynamics and kinematics of out-flowing atmospheres of hot and warm exoplanets is crucial to understanding the origins and evolutionary history of the exoplanets near the evaporation desert. Recently, ground based measurements of the meta-stable Helium atom's resonant absorption at 10830 Å~has become a powerful probe of the base environment which is driving the outflow of exoplane… ▽ More Understanding the dynamics and kinematics of out-flowing atmospheres of hot and warm exoplanets is crucial to understanding the origins and evolutionary history of the exoplanets near the evaporation desert. Recently, ground based measurements of the meta-stable Helium atom's resonant absorption at 10830 Å~has become a powerful probe of the base environment which is driving the outflow of exoplanet atmospheres. We report evidence for the He I 10830 Å~in absorption (equivalent width $\sim$ $0.012 \pm 0.002$ Å) in the exosphere of a warm Neptune orbiting the M-dwarf GJ 3470, during three transits using the Habitable Zone Planet Finder (HPF) near infrared spectrograph. This marks the first reported evidence for He I 10830 Å\, atmospheric absorption for a planet orbiting an M-dwarf. Our detected absorption is broad and its blueshifted wing extends to -36 km/sec, the largest reported in the literature to date. We modelled the state of Helium atoms in the exosphere of GJ3470b based on assumptions on the UV and X-ray flux of GJ 3470, and found our measurement of flux-weighted column density of meta-stable state Helium $(N_{He^2_3S} = 2.4 \times 10^{10} \mathrm{cm^{-2}})$, derived from our transit observations, to be consistent with model, within its uncertainties. The methodology developed here will be useful to study and constrain the atmospheric outflow models of other exoplanets like GJ 3470b which are near the edge of the evaporation desert. △ Less

Submitted 30 March, 2020; v1 submitted 4 October, 2019; originally announced October 2019.

Comments: Accepted in ApJ

Realizing the potential of astrostatistics and astroinformatics

Authors: Gwendolyn Eadie, Thomas J. Loredo, Ashish A. Mahabal, Aneta Siemiginowska, Eric Feigelson, Eric B. Ford, S. G. Djorgovski, Matthew Graham, Zeljko Ivezic, Kirk Borne, Jessi Cisewski-Kehe, J. E. G. Peek, Chad Schafer, Padma A. Yanamandra-Fisher, C. Alex Young

Abstract: This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy. This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: 14 pages, 1 figure; submitted to the Decadal Survey on Astronomy and Astrophysics (Astro2020) on 10 July 2019; see https://sites.nationalacademies.org/DEPS/Astro2020/DEPS_192906

An Automated Method to Detect Transiting Circumbinary Planets

Authors: Diana Windemuth, Eric Agol, Josh Carter, Eric B. Ford, Nader Haghighipour, Jerome A. Orosz, William F. Welsh

Abstract: To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we descri… ▽ More To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we describe an efficient and automated technique for detecting circumbinary planets that transit their binary hosts in Kepler light curves. Our method accounts for large transit timing and duration variations (TTVs and TDVs), induced by binary reflex motion, in two ways: 1) We directly correct for large-scale TTVs and TDVs in the light curves by using Keplerian models to approximate binary and CBP orbits; and 2) We allow additional aperiodicities on the corrected light curves by employing the Quasi-periodic Automated Transit Search algorithm (QATS). We demonstrate that our method dramatically improves detection significance using simulated data and two previously identified CBP systems, Kepler-35 and Kepler-64. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 14 pages, 8 figures. Accepted to MNRAS

Sensitivity Analyses of Exoplanet Occurrence Rates from Kepler and Gaia

Authors: Megan Shabram, Natalie Batalha, Susan E. Thompson, Danley C. Hsu, Eric B. Ford, Jessie L. Christiansen, Daniel Huber, Travis Berger, Joseph Catanzarite, Benjamin E. Nelson, Steve Bryson, Ruslan Belikov, Chris Burke, Doug Caldwell

Abstract: We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the genera… ▽ More We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the general features of the outcome of the planet formation and evolution around GK stars, and provides an infrastructure to use the $Kepler$ results to constrain analytic planet distribution models. We report an increased mean and variance in the marginal posterior distributions for the number of planets per $GK$ star when including planet radius measurement uncertainties. We estimate the number of planets-per-$GK$ star between 0.75 and 2.5 $R_{\oplus}$ and 50 to 300 day orbital periods to have a $68\%$ credible interval of $0.49$ to $0.77$ and a posterior mean of $0.63$. This posterior has a smaller mean and a larger variance than the occurrence rate calculated in this work and in Burke et al. (2015) for the same parameter space using the $Q1-Q16$ (previous $Kepler$ planet candidate and stellar catalog). We attribute the smaller mean to many of the instrumental false positives at longer orbital periods being removed from the $DR25$ catalog. We find that the accuracy and precision of our hierarchical Bayesian model posterior distributions are less sensitive to the total number of planets in the sample, and more so on the characteristics of the catalog completeness and reliability and the span of the planet parameter space. △ Less

Submitted 26 September, 2020; v1 submitted 1 August, 2019; originally announced August 2019.

Comments: 16 pages, 3 figures, published in the Astrophysical Journal

Journal ref: Astrophys. J., 160, 1 (2020)

Scheduling Discovery in the 2020s

Authors: Eric C. Bellm, Eric B. Ford, Aaron Tohuvavohu, Michael W. Coughlin, Brett Morris, Bryan Miller, Jennifer Sobeck, Reed Riddle, Chuanfei Dong, Peter Yoachim

Abstract: The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis. The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: 10 pages; Astro2020 APC White Paper

Architectures of Exoplanetary Systems. I: A Clustered Forward Model for Exoplanetary Systems around Kepler's FGK Stars

Authors: Matthias Y. He, Eric B. Ford, Darin Ragozzine

Abstract: Observations of exoplanetary systems provide clues about the intrinsic distribution of planetary systems, their architectures, and how they formed. We develop a forward modelling framework for generating populations of planetary systems and "observed" catalogues by simulating the Kepler detection pipeline (SysSim). We compare our simulated catalogues to the Kepler DR25 catalogue of planet candidat… ▽ More Observations of exoplanetary systems provide clues about the intrinsic distribution of planetary systems, their architectures, and how they formed. We develop a forward modelling framework for generating populations of planetary systems and "observed" catalogues by simulating the Kepler detection pipeline (SysSim). We compare our simulated catalogues to the Kepler DR25 catalogue of planet candidates, updated to include revised stellar radii from Gaia DR2. We constrain our model based on the observed 1D marginal distributions of orbital periods, period ratios, transit depths, transit depth ratios, transit durations, transit duration ratios, and transit multiplicities. Models assuming planets with independent periods and sizes do not adequately account for the properties of the multiplanet systems. Instead, a clustered point process model for exoplanet periods and sizes provides a significantly better description of the Kepler population, particularly the observed multiplicity and period ratio distributions. We find that $0.56^{+0.18}_{-0.15}$ of FGK stars have at least one planet larger than $0.5 R_\oplus$ between 3 and 300 d. Most of these planetary systems ($\sim 98\%$) consist of one or two clusters with a median of three planets per cluster. We find that the Kepler dichotomy is evidence for a population of highly inclined planetary systems and is unlikely to be solely due to a population of intrinsically single planet systems. We provide a large ensemble of simulated physical and observed catalogues of planetary systems from our models, as well as publicly available code for generating similar catalogues given user-defined parameters. △ Less

Submitted 14 November, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: Published in MNRAS, 31 pages, 9 figures, plus supplemental online materials with 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2019, 490, 4575-4605

Discovery of a Third Transiting Planet in the Kepler-47 Circumbinary System

Authors: Jerome A. Orosz, William F. Welsh, Nader Haghighipour, Billy Quarles, Donald R. Short, Sean M. Mills, Suman Satyal, Guillermo Torres, Eric Agol, Daniel C. Fabrycky, Daniel Jontof-Hutter, Gur Windmiller, Tobias W. A. Müller, Tobias C. Hinse, William D. Cochran, Michael Endl, Eric B. Ford, Tsevi Mazeh, Jack J. Lissauer

Abstract: Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the… ▽ More Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1σ$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $ρ_{\rm middle} < 0.68$ g cm$^{-3}$ and $ρ_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1σ$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 68 pages, 30 figures

Journal ref: AJ, 157, 174 (2019)

Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics

Authors: A. Siemiginowska, G. Eadie, I. Czekala, E. Feigelson, E. B. Ford, V. Kashyap, M. Kuhn, T. Loredo, M. Ntampaka, A. Stevens, A. Avelino, K. Borne, T. Budavari, B. Burkhart, J. Cisewski-Kehe, F. Civano, I. Chilingarian, D. A. van Dyk, G. Fabbiano, D. P. Finkbeiner, D. Foreman-Mackey, P. Freeman, A. Fruscione, A. A. Goodman, M. Graham , et al. (27 additional authors not shown)

Abstract: Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New met… ▽ More Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers