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

Conditions for Changing-Look AGNs from Accretion Disk-Induced Tidal Disruption Events

Authors: Yihan Wang, Matthew J. Graham, K. E. Saavik Ford, Barry McKernan, Taeho Ryu, Daniel Stern

Abstract: The phenomenon of changing-look (CL) behavior in active galactic nuclei (AGN) is characterized by dramatic changes in luminosity and/or emission line profiles over relatively short periods, ranging from months to years. The origin of CL-AGNs remains a mystery, but one proposed explanation involves the response of the inner AGN disk to tidal disruption events (TDEs) around the supermassive black ho… ▽ More The phenomenon of changing-look (CL) behavior in active galactic nuclei (AGN) is characterized by dramatic changes in luminosity and/or emission line profiles over relatively short periods, ranging from months to years. The origin of CL-AGNs remains a mystery, but one proposed explanation involves the response of the inner AGN disk to tidal disruption events (TDEs) around the supermassive black hole (SMBH). In this Letter, we calculate the predicted frequency of AGN TDEs as a function of SMBH mass and compare the results to the observed CL-AGN distribution. We find that if the fraction of CL-AGNs caused by AGN-TDEs is high, then: (1) most SMBHs in CL-AGN are near maximal spin, with the dimensionless spin parameter $a>0.9$; (2) AGN inner disks have a high surface density ($\geq 10^{7}\, {\rm g\, cm^{-2}}$); (3) typical AGN lifetimes are $\sim 10$-$100$ Myr; and (4) a nuclear star cluster initial mass function (IMF) that scales as $\sim m_*^{-1.6}$ is preferred. Future observations of CL-AGN will help constrain the fraction of CL-AGNs caused by AGN-TDEs, SMBH spins, AGN lifetimes, and the nuclear star cluster IMF. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Comments are welcome

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

The magnetically quiet solar surface dominates HARPS-N solar RVs during low activity

Authors: Ben S. Lakeland, Tim Naylor, Raphaëlle Haywood, Nadège Meunier, Federica Rescigno, Shweta Dalal, Annelies Mortier, Samantha J. Thompson, Andrew Collier Cameron, Xavier Dumusque, Mercedes López-Morales, Francesco Pepe, Ken Rice, Alessandro Sozzetti, Stéphane Udry, Eric Ford, Adriano Ghedina, Marcello Lodi

Abstract: Using images from the Helioseismic and Magnetic Imager aboard the \textit{Solar Dynamics Observatory} (SDO/HMI), we extract the radial-velocity (RV) signal arising from the suppression of convective blue-shift and from bright faculae and dark sunspots transiting the rotating solar disc. We remove these rotationally modulated magnetic-activity contributions from simultaneous radial velocities obser… ▽ More Using images from the Helioseismic and Magnetic Imager aboard the \textit{Solar Dynamics Observatory} (SDO/HMI), we extract the radial-velocity (RV) signal arising from the suppression of convective blue-shift and from bright faculae and dark sunspots transiting the rotating solar disc. We remove these rotationally modulated magnetic-activity contributions from simultaneous radial velocities observed by the HARPS-N solar feed to produce a radial-velocity time series arising from the magnetically quiet solar surface (the 'inactive-region radial velocities'). We find that the level of variability in the inactive-region radial velocities remains constant over the almost 7 year baseline and shows no correlation with well-known activity indicators. With an RMS of roughly 1 m/s, the inactive-region radial-velocity time series dominates the total RV variability budget during the decline of solar cycle 24. Finally, we compare the variability amplitude and timescale of the inactive-region radial velocities with simulations of supergranulation. We find consistency between the inactive-region radial-velocity and simulated time series, indicating that supergranulation is a significant contribution to the overall solar radial velocity variability, and may be the main source of variability towards solar minimum. This work highlights supergranulation as a key barrier to detecting Earth twins. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 12 pages, 11 figures, accepted to MNRAS

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

Tidal Disruption Events from three-body scatterings and eccentricity pumping in the disks of Active Galactic Nuclei

Authors: Chaitanya Prasad, Yihan Wang, Rosalba Perna, K. E. Saavik Ford, Barry McKernan

Abstract: Tidal Disruption Events (TDEs) are routinely observed in quiescent galaxies, as stars from the nuclear star cluster are scattered into the loss cone of the central supermassive black hole (SMBH). TDEs are also expected to occur in Active Galactic Nuclei (AGN), due to scattering or orbital eccentricity pumping of stars embedded in the innermost regions of the AGN accretion disk. Encounters with emb… ▽ More Tidal Disruption Events (TDEs) are routinely observed in quiescent galaxies, as stars from the nuclear star cluster are scattered into the loss cone of the central supermassive black hole (SMBH). TDEs are also expected to occur in Active Galactic Nuclei (AGN), due to scattering or orbital eccentricity pumping of stars embedded in the innermost regions of the AGN accretion disk. Encounters with embedded stellar-mass black holes (BH) can result in AGN $μ$TDEs. AGN TDEs and $μ$TDEs could therefore account for a fraction of observed AGN variability. Here, by performing scattering experiments with the few-body code {\tt SpaceHub}, we compute the probability of AGN TDEs and $μ$TDEs as a result of 3-body interactions between stars and binary BHs. We find that AGN TDEs are more probable during the early life of the AGNs, when rates are $\sim (6\times 10^{-5}-5 \times 10^{-2}) (f_\bullet/0.01)$ $~\rm{AGN}^{-1}$~yr$^{-1}$ (where $f_\bullet$ is the ratio between the number density of BHs and stars), generally higher than in quiescent galactic nuclei. By contrast, $μ$TDEs should occur throughout the AGN lifetime at a rate of $\sim (1\times 10^{-4} - 4\times 10^{-2} (f_\bullet/0.01)$ $~\rm{AGN}^{-1}$~yr$^{-1}$. Detection and characterization of AGN TDEs and $μ$ AGN TDEs with future surveys using {\em Rubin} and {\em Roman} will help constrain the populations of stars and compact objects embedded in AGN disks, a key input for the LVK AGN channel. △ Less

Submitted 23 May, 2024; v1 submitted 26 September, 2023; originally announced October 2023.

Comments: Final version accepted by MNRAS

Constraining the LIGO/Virgo AGN channel with black hole spins

Authors: B. McKernan, K. E. S. Ford

Abstract: Merging black holes (BH) are expected to produce remnants with large dimensionless spin parameters ($a_{\rm spin} \sim 0.7$). However, gravitational wave (GW) observations with LIGO/Virgo suggest that merging BH are consistent with modestly positive but not high spin ($a_{\rm spin} \sim 0.2$), causing tension with models suggesting that high mass mergers are produced by hierarchical merger channel… ▽ More Merging black holes (BH) are expected to produce remnants with large dimensionless spin parameters ($a_{\rm spin} \sim 0.7$). However, gravitational wave (GW) observations with LIGO/Virgo suggest that merging BH are consistent with modestly positive but not high spin ($a_{\rm spin} \sim 0.2$), causing tension with models suggesting that high mass mergers are produced by hierarchical merger channels. Some BH also show evidence for strong in-plane spin components. Here we point out that \emph{spin down} of BH due to eccentric prograde post-merger orbits within the gas of an active galactic nucleus (AGN) disk can yield BH with masses in the upper mass gap, but only modestly positive $a_{\rm spin}$, and thus observations of BH with low spin \emph{do not} rule out hierarchical models. We also point out that the fraction of BBH mergers with significant in-plane spin components is a strong test of interactions between disk binary black holes (BBH) and nuclear spheroid orbiters. Spin magnitude and spin tilt constraints from LIGO/Virgo observations of BBH are an excellent test of dynamics of black holes in AGN disks, disk properties and the nuclear clusters interacting with AGN. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 7 pages, 3 figures, MNRAS submitted

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

NuSTAR Observations of Candidate Subparsec Binary Supermassive Black Holes

Authors: M. Lynne Saade, Murray Brightman, Daniel Stern, Thomas Connor, S. G. Djorgovski, Daniel J. D'Orazio, K. E. S. Ford, Matthew J. Graham, Zoltan Haiman, Hyunsung D. Jun, Elias Kammoun, Ralph P. Kraft, Barry McKernan, Alexei Vikhlinin, Dominic J. Walton

Abstract: We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously o… ▽ More We present analysis of NuSTAR X-ray observations of three AGN that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low energy X-ray properties and the larger AGN population. However some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs. △ Less

Submitted 27 March, 2024; v1 submitted 12 April, 2023; originally announced April 2023.

Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Authors: Joshua Fagin, Ji Won Park, Henry Best, James Hung-Hsu Chan, K. E Saavik Ford, Matthew J. Graham, V. Ashley Villar, Shirley Ho, Matthew O'Dowd

Abstract: Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described… ▽ More Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling. △ Less

Submitted 16 April, 2024; v1 submitted 9 April, 2023; originally announced April 2023.

Comments: 22 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 965, Number 2, April 2024

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.

An extreme test case for planet formation: a close-in Neptune orbiting an ultracool star

Authors: Gudmundur Stefansson, Suvrath Mahadevan, Yamila Miguel, Paul Robertson, Megan Delamer, Shubham Kanodia, Caleb Cañas, Joshua Winn, Joe Ninan, Ryan Terrien, Rae Holcomb, Eric Ford, Brianna Zawadzki, Brendan P. Bowler, Chad Bender, William Cochran, Scott Diddams, Michael Endl, Connor Fredrick, Samuel Halverson, Fred Hearty, Gary J. Hill, Andrea Lin, Andrew Metcalf, Andrew Monson , et al. (5 additional authors not shown)

Abstract: In current theories of planet formation, close-orbiting planets as massive as Neptune are expected to be very rare around low-mass stars. We report the discovery of a Neptune-mass planet orbiting the `ultracool' star LHS 3154, which is nine times less massive than the Sun. The planet's orbital period is 3.7 days and its minimum mass is 13.2 Earth masses, giving it the largest known planet-to-star… ▽ More In current theories of planet formation, close-orbiting planets as massive as Neptune are expected to be very rare around low-mass stars. We report the discovery of a Neptune-mass planet orbiting the `ultracool' star LHS 3154, which is nine times less massive than the Sun. The planet's orbital period is 3.7 days and its minimum mass is 13.2 Earth masses, giving it the largest known planet-to-star mass ratio among short-period planets ($<$\,100 days) orbiting ultracool stars. Both the core accretion and gravitational instability theories for planet formation struggle to account for this system. In the core-accretion scenario, in particular, the dust mass of the protoplanetary disk would need to be an order of magnitude higher than typically seen in protoplanetary disk observations of ultracool stars. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: Original Manuscript as submitted to Science on Oct 17, 2022. In review

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- IV. Aperture Masking Interferometry

Authors: Anand Sivaramakrishnan, Peter Tuthill, James P. Lloyd, Alexandra Z. Greenbaum, Deepashri Thatte, Rachel A. Cooper, Thomas Vandal, Jens Kammerer, Joel Sanchez-Bermudez, Benjamin J. S. Pope, Dori Blakely, Loïc Albert, Neil J. Cook, Doug Johnstone, André R. Martel, Kevin Volk, Anthony Soulain, Étienne Artigau, David Lafrenière, Chris J. Willott, Sébastien Parmentier, K. E. Saavik Ford, Barry McKernan, M. Begoña Vila, Neil Rowlands , et al. (14 additional authors not shown)

Abstract: The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, pres… ▽ More The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, present some sample science cases, explain its operational observing strategies, indicate how AMI proposals can be developed with data simulations, and how AMI data can be analyzed. We also present key results from commissioning AMI. Since the allied Kernel Phase Imaging (KPI) technique benefits from AMI operational strategies, we also cover NIRISS KPI methods and analysis techniques, including a new user-friendly KPI pipeline. The NIRISS KPI bright limit is $\simeq 8$ W2 magnitudes. AMI (and KPI) achieve an inner working angle of $\sim 70$ mas that is well inside the $\sim 400$ mas NIRCam inner working angle for its circular occulter coronagraphs at comparable wavelengths. △ Less

Submitted 7 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 30 pages, 10 figures

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

A light in the dark: searching for electromagnetic counterparts to black hole-black hole mergers in LIGO/Virgo O3 with the Zwicky Transient Facility

Authors: Matthew J. Graham, Barry McKernan, K. E. Saavik Ford, Daniel Stern, S. G. Djorgovski, Michael Coughlin, Kevin B. Burdge, Eric C. Bellm, George Helou, Ashish A. Mahabal, Frank J. Masci, Josiah Purdum, Philippe Rosnet, Ben Rusholme

Abstract: The accretion disks of active galactic nuclei (AGN) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high density environment, mergers within AGN are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGN with unusual flaring activity observed by… ▽ More The accretion disks of active galactic nuclei (AGN) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high density environment, mergers within AGN are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGN with unusual flaring activity observed by the Zwicky Transient Facility (ZTF), we describe a search for candidate EM counterparts to binary black hole (BBH) mergers detected by LIGO/Virgo in O3. After removing probable false positives, we find nine candidate counterparts to BBH mergers mergers during O3 (seven in O3a, two in O3b) with a $p$-value of 0.019. Based on ZTF sky coverage, AGN geometry, and merger geometry, we expect $\approx 3(N_{\rm BBH}/83)(f_{\rm AGN}/0.5)$ potentially detectable EM counterparts from O3, where $N_{\rm BBH}$ is the total number of observed BBH mergers and $f_{\rm AGN}$ is the fraction originating in AGN. Further modeling of breakout and flaring phenomena in AGN disks is required to reduce our false positive rate. Two of the events are also associated with mergers with total masses $> 100M_\odot$, which is the expected rate for O3 if hierarchical (large mass) mergers occur in the AGN channel. Candidate EM counterparts in future GW observing runs can be better constrained by coverage of the Southern sky as well as spectral monitoring of unusual AGN flaring events in LIGO/Virgo alert volumes. A future set of reliable AGN EM counterparts to BBH mergers will yield an independent means of measuring cosmic expansion ($H_0$) as a function of redshift. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 24 pages, 6 figures, submitted to ApJ

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/)

From Data to Software to Science with the Rubin Observatory LSST

Authors: Katelyn Breivik, Andrew J. Connolly, K. E. Saavik Ford, Mario Jurić, Rachel Mandelbaum, Adam A. Miller, Dara Norman, Knut Olsen, William O'Mullane, Adrian Price-Whelan, Timothy Sacco, J. L. Sokoloski, Ashley Villar, Viviana Acquaviva, Tomas Ahumada, Yusra AlSayyad, Catarina S. Alves, Igor Andreoni, Timo Anguita, Henry J. Best, Federica B. Bianco, Rosaria Bonito, Andrew Bradshaw, Colin J. Burke, Andresa Rodrigues de Campos , et al. (75 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the po… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

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

Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc

Authors: Syeda S. Nasim, Gaia Fabj, Freddy Caban, Amy Secunda, K. E. Saavik Ford, Barry McKernan, Jillian M. Bellovary, Nathan W. C. Leigh, Wladimir Lyra

Abstract: Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population w… ▽ More Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population which drives sBH-sBH mergers detectable in gravitational waves using LIGO-Virgo-KAGRA (LVK) or sBH-SMBH mergers detectable with LISA (Laser Interferometer Space Antenna). Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive 'immortal' stars. Here, we investigate the behaviour of polar and retrograde orbiters $(i \geq 90^{\circ})$ interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation $(i<90^{\circ})$ during capture, and decrease $a$ dramatically towards the SMBH. For sBH, we find a critical angle $i_{\rm ret} \sim 113^{\circ}$, below which retrograde sBH decay towards embedded prograde orbits $(i \to 0^{\circ})$, while for $i_{\rm o}>i_{\rm ret}$ sBH decay towards embedded retrograde orbits $(i \to 180^{\circ})$. sBH near polar orbits $(i \sim 90^{\circ})$ and stars on nearly embedded retrograde orbits $(i \sim 180^{\circ})$ show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture-time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution. △ Less

Submitted 17 May, 2023; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures, 1 table

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, July 2023, Pages 5393-5401

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

Effects of an Immortal Stellar Population in AGN Disks

Authors: Adam S. Jermyn, Alexander J. Dittmann, B. McKernan, K. E. S. Ford, Matteo Cantiello

Abstract: Stars are likely embedded in the gas disks of Active Galactic Nuclei (AGN). Theoretical models predict that in the inner regions of the disk these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial… ▽ More Stars are likely embedded in the gas disks of Active Galactic Nuclei (AGN). Theoretical models predict that in the inner regions of the disk these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial AGN disk. We find numbers of order $10^{2-4}$, confined to the inner $r_{\rm cap} \sim 3000 r_s \sim 0.03\rm pc$. These stars can profoundly alter the chemistry of AGN disks, enriching them in helium and depleting them in hydrogen, both by order-unity amounts. We further consider mergers between these stars and other disk objects, suggesting that star-star mergers result in rapid mass loss from the remnant to restore an equilibrium mass, while star-compact object mergers may result in exotic outcomes and even host binary black hole mergers within themselves. Finally, we examine how these stars react as the disk dissipates towards the end of its life, and find that they may return mass to the disk fast enough to extend its lifetime by a factor of several and/or may drive powerful outflows from the disk. Post-AGN, these stars rapidly lose mass and form a population of stellar mass black holes around $10M_{\odot}$. Due to the complex and uncertain interactions between embedded stars and the disk, their plausible ubiquity, and their order unity impact on disk structure and evolution, they must be included in realistic disk models. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 12 pages, 2 figures. Accepted in ApJ

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.

GW190521 as a black-hole merger coincident with the ZTF19abanrhr flare

Authors: Juan Calderón Bustillo, Samson H. W. Leong, Koustav Chandra, Barry McKernan, K. E. S. Ford

Abstract: We present an analysis that reconciles the gravitational-wave signal GW190521 observed by the Advanced LIGO and Advanced Virgo detectors with the electromagnetic flare ZTF19abanrhr observed by the Zwicky Transient Facility. We analyze GW190521 under a mass-ratio prior uniform in $Q \in [1,4]$ and using the state-of-the-art waveform model for black-hole mergers \texttt{NRSur7dq4}. We find a $90\%$… ▽ More We present an analysis that reconciles the gravitational-wave signal GW190521 observed by the Advanced LIGO and Advanced Virgo detectors with the electromagnetic flare ZTF19abanrhr observed by the Zwicky Transient Facility. We analyze GW190521 under a mass-ratio prior uniform in $Q \in [1,4]$ and using the state-of-the-art waveform model for black-hole mergers \texttt{NRSur7dq4}. We find a $90\%$ credible region for the black-hole masses extending far outside what originally reported by \cite{GW190521D}, where our maximum likelihood masses reside. We find a $15\%$ probability that both black holes avoid the pair-instability supernova gap. We infer a three-dimensional sky-location highly consistent with ZTF19abanrhr, obtaining an odds-ratio ${\cal{O}}_{C/R}=72:1$ that strongly favors the hypothesis of a true coincidence over a random one. Combining this event with the neutron-star merger GW170817, we estimate a Hubble constant H$_0=72.1^{+10.6}_{-6.4}\mathrm{km\,s^{-1}\,Mpc^{-1}}$ at the $68\%$ credible level. △ Less

Submitted 23 December, 2021; originally announced December 2021.

Comments: 14 Pages, 11 Figures

Report number: LIGO-P2100467

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)

Science with the Ultraviolet Explorer (UVEX)

Authors: S. R. Kulkarni, Fiona A. Harrison, Brian W. Grefenstette, Hannah P. Earnshaw, Igor Andreoni, Danielle A. Berg, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Jessie L. Christiansen, Michael W. Coughlin, Alexander Wuollet Criswell, Behnam Darvish, Kaustav K. Das, Kishalay De, Luc Dessart, Don Dixon, Bas Dorsman, Kareem El-Badry, Christopher Evans, K. E. Saavik Ford, Christoffer Fremling, Boris T. Gansicke, Suvi Gezari, Y. Goetberg , et al. (31 additional authors not shown)

Abstract: UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky… ▽ More UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program. △ Less

Submitted 17 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 69 pages, 43 figures

Measuring the properties of active galactic nuclei disks with gravitational waves

Authors: Avi Vajpeyi, Eric Thrane, Rory Smith, Barry McKernan, K. E. Saavik Ford

Abstract: Active galactic nuclei (AGN) are promising environments for the assembly of merging binary black hole (BBH) systems. Interest in AGNs as nurseries for merging BBH is rising following the detection of gravitational waves from a BBH system from the purported pair-instability mass gap, most notably, GW190521. Active galactic nuclei have also been invoked to explain the formation of the high-mass-rati… ▽ More Active galactic nuclei (AGN) are promising environments for the assembly of merging binary black hole (BBH) systems. Interest in AGNs as nurseries for merging BBH is rising following the detection of gravitational waves from a BBH system from the purported pair-instability mass gap, most notably, GW190521. Active galactic nuclei have also been invoked to explain the formation of the high-mass-ratio system, GW190814. We draw on simulations of BBH systems in AGN to propose a phenomenological model for the distribution of black hole spins of merging binaries in AGN disks. The model incorporates distinct features that make the AGN channel potentially distinguishable from other channels, such as assembly in the field and in globular clusters. The model parameters can be mapped heuristically to the age and density of AGN disks. We estimate the extent to which different populations of mergers in AGNs can be distinguished. If most merging black holes are assembled in AGNs, future gravitational-wave observations may provide insights into the dynamics of AGN disks. △ Less

Submitted 2 March, 2023; v1 submitted 6 November, 2021; originally announced November 2021.

Comments: 10 pages, 4 figures

Journal ref: The Astrophysical Journal, Volume 931, Issue 2, id.82, 8 pp., June 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

Starfall: A heavy rain of stars in 'turning on' AGN

Authors: B. McKernan, K. E. S. Ford, M. Cantiello, M. J. Graham, A. S. Jermyn, N. W. C. Leigh, T. Ryu, D. Stern

Abstract: As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and… ▽ More As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and onto the SMBH. A large AGN loss cone ($θ_{\rm AGN,lc}$) can result from binary scatterings in the inner disk and yield tidal disruption events (TDEs). Signatures of starfall include optical/UV flares that rise in luminosity over time, particularly in the inner disk. If the SMBH mass is $M_{\rm SMBH} \ge 10^{8}M_{\odot}$, flares truncate abruptly and the star is swallowed. If $M_{\rm SMBH}<10^{8}M_{\odot}$, and if the infalling orbit lies within $θ_{\rm AGN,lc}$, the flare is followed by a TDE which can be prograde or retrograde relative to the AGN inner disk. Retrograde AGN TDEs are over-luminous and short-lived as in-plane ejecta collide with the inner disk and a lower AGN state follows. Prograde AGN TDEs add angular momentum to inner disk gas and so start off looking like regular TDEs but are followed by an AGN high state. Searches for such flare signatures test models of AGN `turn on', SMBH mass, as well as disk properties and the embedded population. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: 10 pages, 1 figure, MNRAS submitted

Binary Black Hole Merger Rates in AGN Disks versus Nuclear Star Clusters: Loud beats Quiet

Authors: K. E. Saavik Ford, Barry McKernan

Abstract: Galactic nuclei are promising sites for stellar origin black hole (BH) mergers, as part of merger hierarchies in deep potential wells. We show that binary black hole (BBH) merger rates in active galactic nuclei (AGN) should always exceed merger rates in quiescent galactic nuclei (nuclear star clusters, NSCs) around supermassive BHs (SMBHs) without accretion disks. This is primarily due to average… ▽ More Galactic nuclei are promising sites for stellar origin black hole (BH) mergers, as part of merger hierarchies in deep potential wells. We show that binary black hole (BBH) merger rates in active galactic nuclei (AGN) should always exceed merger rates in quiescent galactic nuclei (nuclear star clusters, NSCs) around supermassive BHs (SMBHs) without accretion disks. This is primarily due to average binary lifetimes in AGN that are significantly shorter than in NSCs. The lifetime difference comes from rapid hardening of BBHs in AGN, such that their semi-major axes are smaller than the hard-soft boundary of their parent NSC; this contrasts with the large average lifetime to merger for BBHs in NSCs around SMBHs, due to binary ionization mechanisms. Secondarily, merger rates in AGNs are enhanced by gas-driven binary formation mechanisms. Formation of new BHs in AGN disks are a minor contributor to the rate differences. With the gravitational wave detection of several BBHs with at least one progenitor in the upper mass gap, and signatures of dynamical formation channels in the $χ_{\rm eff}$ distribution, we argue that AGN could contribute $\sim 25\%-80\%$ of the LIGO-Virgo measured rate of $\sim 24 \rm{Gpc}^{-3} \rm{yr}^{-1}$. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 8 pages, 1 figure, submitted to MNRAS

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/

LIGO--Virgo correlations between mass ratio and effective inspiral spin: testing the active galactic nuclei channel

Authors: B. McKernan, K. E. S. Ford, T. Callister, W. M. Farr, R. O'Shaughnessy, R. Smith, E. Thrane, A. Vajpeyi

Abstract: Observations by LIGO--Virgo of binary black hole mergers suggest a possible anti-correlation between black hole mass ratio ($q=m_{2}/m_{1}$) and the effective inspiral spin parameter $χ_{\rm eff}$, the mass-weighted spin projection onto the binary orbital angular momentum (Callister et al. 2021). We show that such an anti-correlation can naturally occur for binary black holes assembled in active g… ▽ More Observations by LIGO--Virgo of binary black hole mergers suggest a possible anti-correlation between black hole mass ratio ($q=m_{2}/m_{1}$) and the effective inspiral spin parameter $χ_{\rm eff}$, the mass-weighted spin projection onto the binary orbital angular momentum (Callister et al. 2021). We show that such an anti-correlation can naturally occur for binary black holes assembled in active galactic nuclei (AGN) due to spherical and planar symmetry-breaking effects. We describe a phenomenological model in which: 1) heavier black holes live in the AGN disk and tend to spin up into alignment with the disk; 2) lighter black holes with random spin orientations live in the nuclear spheroid; 3) the AGN disk is dense enough to rapidly capture a fraction of the spheroid component. but small in radial extent to limit the number of bulk disk mergers; 4) migration within the disk is non-uniform, likely disrupted by feedback from migrators or disk turbulence; 5) dynamical encounters in the disk are common and preferentially disrupt binaries that are retrograde around their center of mass, particularly at stalling orbits, or traps. This model may explain trends in LIGO--Virgo data while offering falsifiable predictions. Comparisons of predictions in ($q,χ_{\rm eff}$) parameter space for the different channels may allow us to distinguish their fractional contributions to the observed merger rates. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 7 pages, 4 Figures, MNRAS (submitted)

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)

Stellar Activity Manifesting at a One Year Alias Explains Barnard b as a False Positive

Authors: Jack Lubin, Paul Robertson, Gudmundur Stefansson, Joe Ninan, Suvrath Mahadevan, Michael Endl, Eric Ford, Jason T. Wright, Corey Beard, Chad Bender, William D. Cochran, Scott A. Diddams, Connor Fredrick, Samuel Halverson, Shubham Kanodia, Andrew J. Metcalf, Lawrence Ramsey, Arpita Roy, Christian Schwab, Ryan Terrien

Abstract: Barnard's star is among the most studied stars given its proximity to the Sun. It is often considered $the$ Radial Velocity (RV) standard for fully convective stars due to its RV stability and equatorial declination. Recently, an $M \sin i = 3.3 M_{\oplus}$ super-Earth planet candidate with a 233 day orbital period was announced by Ribas et al. (2018). New observations from the near-infrared Habit… ▽ More Barnard's star is among the most studied stars given its proximity to the Sun. It is often considered $the$ Radial Velocity (RV) standard for fully convective stars due to its RV stability and equatorial declination. Recently, an $M \sin i = 3.3 M_{\oplus}$ super-Earth planet candidate with a 233 day orbital period was announced by Ribas et al. (2018). New observations from the near-infrared Habitable-zone Planet Finder (HPF) Doppler spectrometer do not show this planetary signal. We ran a suite of experiments on both the original data and a combined original + HPF data set. These experiments include model comparisons, periodogram analyses, and sampling sensitivity, all of which show the signal at the proposed period of 233 days is transitory in nature. The power in the signal is largely contained within 211 RVs that were taken within a 1000 day span of observing. Our preferred model of the system is one which features stellar activity without a planet. We propose that the candidate planetary signal is an alias of the 145 day rotation period. This result highlights the challenge of analyzing long-term, quasi-periodic activity signals over multi-year and multi-instrument observing campaigns. △ Less

Submitted 14 May, 2021; originally announced May 2021.

Comments: 25 pages, 8 figures. Accepted for publication in The Astronomical Journal

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

Predicting the spectrum of UGC 2885, Rubin's Galaxy with machine learning

Authors: Benne W. Holwerda, John F. Wu, William C. Keel, Jason Young, Ren Mullins, Joannah Hinz, K. E. Saavik Ford, Pauline Barmby, Rupali Chandar, Jeremy Bailin, Josh Peek, Tim Pickering, Torsten Böker

Abstract: Wu & Peek (2020) predict SDSS-quality spectra based on Pan-STARRS broad-band \textit{grizy} images using machine learning (ML). In this letter, we test their prediction for a unique object, UGC 2885 ("Rubin's galaxy"), the largest and most massive, isolated disk galaxy in the local Universe ($D<100$ Mpc). After obtaining the ML predicted spectrum, we compare it to all existing spectroscopic inform… ▽ More Wu & Peek (2020) predict SDSS-quality spectra based on Pan-STARRS broad-band \textit{grizy} images using machine learning (ML). In this letter, we test their prediction for a unique object, UGC 2885 ("Rubin's galaxy"), the largest and most massive, isolated disk galaxy in the local Universe ($D<100$ Mpc). After obtaining the ML predicted spectrum, we compare it to all existing spectroscopic information that is comparable to an SDSS spectrum of the central region: two archival spectra, one extracted from the VIRUS-P observations of this galaxy, and a new, targeted MMT/Binospec observation. Agreement is qualitatively good, though the ML prediction prefers line ratios slightly more towards those of an active galactic nucleus (AGN), compared to archival and VIRUS-P observed values. The MMT/Binospec nuclear spectrum unequivocally shows strong emission lines except H$β$, the ratios of which are consistent with AGN activity. The ML approach to galaxy spectra may be a viable way to identify AGN supplementing NIR colors. How such a massive disk galaxy ($M^* = 10^{11}$ M$_\odot$), which uncharacteristically shows no sign of interaction or mergers, manages to fuel its central AGN remains to be investigated. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 9 pages, 5 figures, submitted to ApJL

The missing link in gravitational-wave astronomy: A summary of discoveries waiting in the decihertz range

Authors: Manuel Arca Sedda, Christopher P L Berry, Karan Jani, Pau Amaro-Seoane, Pierre Auclair, Jonathon Baird, Tessa Baker, Emanuele Berti, Katelyn Breivik, Chiara Caprini, Xian Chen, Daniela Doneva, Jose M Ezquiaga, K E Saavik Ford, Michael L Katz, Shimon Kolkowitz, Barry McKernan, Guido Mueller, Germano Nardini, Igor Pikovski, Surjeet Rajendran, Alberto Sesana, Lijing Shao, Nicola Tamanini, Niels Warburton , et al. (3 additional authors not shown)

Abstract: Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of m… ▽ More Since 2015 the gravitational-wave observations of LIGO and Virgo have transformed our understanding of compact-object binaries. In the years to come, ground-based gravitational-wave observatories such as LIGO, Virgo, and their successors will increase in sensitivity, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will provide gravitational-wave observations of massive black holes binaries. Between the $\sim 10$-$10^3~\mathrm{Hz}$ band of ground-based observatories and the $\sim10^{-4}$-$10^{-1}~\mathrm{Hz}$ band of LISA lies the uncharted decihertz gravitational-wave band. We propose a Decihertz Observatory to study this frequency range, and to complement observations made by other detectors. Decihertz observatories are well suited to observation of intermediate-mass ($\sim10^2$-$10^4 M_\odot$) black holes; they will be able to detect stellar-mass binaries days to years before they merge, providing early warning of nearby binary neutron star mergers and measurements of the eccentricity of binary black holes, and they will enable new tests of general relativity and the Standard Model of particle physics. Here we summarise how a Decihertz Observatory could provide unique insights into how black holes form and evolve across cosmic time, improve prospects for both multimessenger astronomy and multiband gravitational-wave astronomy, and enable new probes of gravity, particle physics and cosmology. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: 13 pages, 1 figure. Published in Experimental Astronomy. Summarising white paper arXiv:1908.11375

Journal ref: Experimental Astronomy; 51(3):1427-1440; 2021

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

Rapid Formation of Super-Earths Around Low-Mass Stars

Authors: Brianna Zawadzki, Daniel Carrera, Eric Ford

Abstract: NASA's TESS mission is expected to discover hundreds of M dwarf planets. However, few studies focus on how planets form around low-mass stars. We aim to better characterize the formation process of M dwarf planets to fill this gap and aid in the interpretation of TESS results. We use ten sets of N-body planet formation simulations which vary in whether a gas disc is present, initial range of embry… ▽ More NASA's TESS mission is expected to discover hundreds of M dwarf planets. However, few studies focus on how planets form around low-mass stars. We aim to better characterize the formation process of M dwarf planets to fill this gap and aid in the interpretation of TESS results. We use ten sets of N-body planet formation simulations which vary in whether a gas disc is present, initial range of embryo semi-major axes, and initial solid surface density profile. Each simulation begins with 147 equal-mass embryos around a 0.2 solar mass star and runs for 100 Myr. We find that planets form rapidly, with most collisions occurring within the first 1 Myr. The presence of a gas disc reduces the final number of planets relative to a gas-free environment and causes planets to migrate inward. We find that roughly a quarter of planetary systems experience their final giant impact inside the gas disc, suggesting that some super-Earths may be able to reaccrete an extended gaseous envelope after their final giant impact, though these may be affected by additional process such as photoevaporation. In addition, we find that the final distribution of planets does not retain a memory of the slope of the initial surface density profile, regardless of whether or not a gas disc is present. Thus, our results suggest that present-day observations are unlikely to provide sufficient information to accurately reverse-engineer the initial distribution of solids. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 17 pages, 8 figures, accepted for publication in MNRAS

Prospects of gravitational-waves detections from common-envelope evolution with LISA

Authors: M. Renzo, T. Callister, K. Chatziioannou, L. A. C. van Son, C. M. F. Mingarelli, M. Cantiello, K. E. S. Ford, B. McKernan, G. Ashton

Abstract: Understanding common envelope (CE) evolution is an outstanding problem in binary evolution. Although the CE phase is not driven by gravitational-wave (GW) emission, the in-spiraling binary emits GWs that passively trace the CE dynamics. Detecting this GW signal would provide direct insight into the gas-driven physics. Even a non-detection might offer invaluable constraints. We investigate the pros… ▽ More Understanding common envelope (CE) evolution is an outstanding problem in binary evolution. Although the CE phase is not driven by gravitational-wave (GW) emission, the in-spiraling binary emits GWs that passively trace the CE dynamics. Detecting this GW signal would provide direct insight into the gas-driven physics. Even a non-detection might offer invaluable constraints. We investigate the prospects of detection of a Galactic CE by LISA. While the dynamical phase of the CE is likely sufficiently loud for detection, it is short and thus rare. We focus instead on the self-regulated phase that proceeds on a thermal timescale. Based on population synthesis calculations and the (unknown) signal duration in the LISA band, we expect $\sim 0.1-100$ sources in the Galaxy during the mission duration. We map the GW observable parameter space of frequency $f_\mathrm{GW}$ and its derivative $\dot f_\mathrm{GW}$ remaining agnostic on the specifics of the inspiral, and find that signals with $\mathrm{SNR}>10$ are possible if the CE stalls at separations such that $f_\mathrm{GW}\gtrsim2\times10^{-3}\,\mathrm{Hz}$. We investigate the possibility of misidentifying the signal with other known sources. If the second derivative $\ddot f_\mathrm{GW}$ can also be measured, the signal can be distinguished from other sources using a GW braking-index. Alternatively, coupling LISA with electromagnetic observations of peculiar red giant stars and/or infrared and optical transients might allow for the disentangling of a Galactic CE from other Galactic and extra-galactic GW sources. △ Less

Submitted 1 July, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: Accepted by ApJ -- author's version. Code available at: https://github.com/tcallister/LISA-and-CE-Evolution