Using an ensemble of N-body simulations, this paper considers the fate of the outer gas giants (Jupiter, Saturn, Uranus, and Neptune) after the Sun leaves the main sequence and completes its stellar evolution. Due to solar mass loss—which is expected to remove roughly half of the star’s mass—the orbits of the giant planets expand. This adiabatic process maintains the orbital period ratios, but the mutual interactions between planets and the width of mean-motion resonances (MMR) increase, leading to the capture of Jupiter and Saturn into a stable 5:2 resonant configuration. The expanded orbits, coupled with the large-amplitude librations of the critical MMR angle, make the system more susceptible to perturbations from stellar flyby interactions. Accordingly, within about 30 Gyr, stellar encounters perturb the planets onto the chaotic subdomain of the 5:2 resonance, triggering a large-scale instability, which culminates in the ejections of all but one planet over the subsequent ∼10 Gyr. After an additional ∼50 Gyr, a close stellar encounter (with a perihelion distance less than ∼200 au) liberates the final planet. Through this sequence of events, the characteristic timescale over which the solar system will be completely dissolved is roughly 100 Gyr. Our analysis thus indicates that the expected dynamical lifetime of the solar system is much longer than the current age of the universe, but is significantly shorter than previous estimates.

Recent analyses have shown that distant orbits within the scattered disk population of the Kuiper Belt exhibit an unexpected clustering in their respective arguments of perihelion. While several hypotheses have been put forward to explain this alignment, to date, a theoretical model that can successfully account for the observations remains elusive. In this work we show that the orbits of distant Kuiper Belt objects (KBOs) cluster not only in argument of perihelion, but also in physical space. We demonstrate that the perihelion positions and orbital planes of the objects are tightly confined and that such a clustering has only a probability of 0.007% to be due to chance, thus requiring a dynamical origin. We find that the observed orbital alignment can be maintained by a distant eccentric planet with mass ≳10 m _⊕ whose orbit lies in approximately the same plane as those of the distant KBOs, but whose perihelion is 180° away from the perihelia of the minor bodies. In addition to accounting for the observed orbital alignment, the existence of such a planet naturally explains the presence of high-perihelion Sedna-like objects, as well as the known collection of high semimajor axis objects with inclinations between 60° and 150° whose origin was previously unclear. Continued analysis of both distant and highly inclined outer solar system objects provides the opportunity for testing our hypothesis as well as further constraining the orbital elements and mass of the distant planet.

The coolest dwarf stars are intrinsically faint at visible wavelengths and exhibit rotationally modulated stellar activity from spots and plages. It is advantageous to observe these stars at near-infrared (NIR) wavelengths (1–2.5 μm) where they emit the bulk of their bolometric luminosity and are most quiescent. In this work, we describe our methodology and results in obtaining precise radial velocity (RV) measurements of low-mass stars using K-band spectra taken with the R ∼ 80,000 iSHELL spectrograph and the NASA Infrared Telescope Facility using a methane isotopologue gas cell in the calibration unit. Our novel analysis pipeline extracts RVs by minimizing the rms of the residuals between the observed spectrum and a forward model. The model accounts for the gas cell, tellurics, blaze function, multiple sources of quasi-sinusoidal fringing, and line spread function of the spectrograph. The stellar template is derived iteratively using the target observations themselves through averaging barycenter-shifted residuals. We have demonstrated 5 m s ⁻¹ precision over one-year timescales for the M4 dwarf Barnard’s Star and K dwarf 61 Cygni A, and 3 m s ⁻¹ over a month for the M2 dwarf GJ 15 A. This work demonstrates the potential for iSHELL to determine dynamical masses for candidate exoplanets discovered with the NASA Transiting Exoplanet Survey Satellite mission, and to search for exoplanets orbiting moderately active and/or young K & M dwarfs.

We present spectral and photometric observations of 10 Type Ia supernovae (SNe Ia) in the redshift range 0.16 ≤ z ≤ 0.62. The luminosity distances of these objects are determined by methods that employ relations between SN Ia luminosity and light curve shape. Combined with previous data from our High- z Supernova Search Team and recent results by Riess et al., this expanded set of 16 high-redshift supernovae and a set of 34 nearby supernovae are used to place constraints on the following cosmological parameters: the Hubble constant ( H ₀), the mass density (Ω _M ), the cosmological constant (i.e., the vacuum energy density, Ω _Λ), the deceleration parameter ( q ₀), and the dynamical age of the universe ( t ₀). The distances of the high-redshift SNe Ia are, on average, 10%–15% farther than expected in a low mass density (Ω _M = 0.2) universe without a cosmological constant. Different light curve fitting methods, SN Ia subsamples, and prior constraints unanimously favor eternally expanding models with positive cosmological constant (i.e., Ω _Λ > 0) and a current acceleration of the expansion (i.e., q ₀ < 0). With no prior constraint on mass density other than Ω _M ≥ 0, the spectroscopically confirmed SNe Ia are statistically consistent with q ₀ < 0 at the 2.8 σ and 3.9 σ confidence levels, and with Ω _Λ > 0 at the 3.0 σ and 4.0 σ confidence levels, for two different fitting methods, respectively. Fixing a "minimal" mass density, Ω _M = 0.2, results in the weakest detection, Ω _Λ > 0 at the 3.0 σ confidence level from one of the two methods. For a flat universe prior (Ω _M + Ω _Λ = 1), the spectroscopically confirmed SNe Ia require Ω _Λ > 0 at 7 σ and 9 σ formal statistical significance for the two different fitting methods. A universe closed by ordinary matter (i.e., Ω _M = 1) is formally ruled out at the 7 σ to 8 σ confidence level for the two different fitting methods. We estimate the dynamical age of the universe to be 14.2 ± 1.7 Gyr including systematic uncertainties in the current Cepheid distance scale. We estimate the likely effect of several sources of systematic error, including progenitor and metallicity evolution, extinction, sample selection bias, local perturbations in the expansion rate, gravitational lensing, and sample contamination. Presently, none of these effects appear to reconcile the data with Ω _Λ = 0 and q ₀ ≥ 0.

The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

We model the settlement of the Galaxy by space-faring civilizations in order to address issues related to the Fermi Paradox. We are motivated to explore the problem in a way that avoids assumptions about the agency (i.e., questions of intent and motivation) of any exo-civilization seeking to settle other planetary systems. We begin by considering the speed of an advancing settlement front to determine if the Galaxy can become inhabited with space-faring civilizations on timescales shorter than its age. Our models for the front speed include the directed settlement of nearby settleable systems through the launching of probes with a finite velocity and range. We also include the effect of stellar motions on the long-term behavior of the settlement front which adds a diffusive component to its advance. As part of our model we also consider that only a fraction, f, of planets will have conditions amenable to settlement by the space-faring civilization. The results of these models demonstrate that the Milky Way can be readily filled-in with settled stellar systems under conservative assumptions about interstellar spacecraft velocities and launch rates. We then move on to consider the question of the Galactic steady state achieved in terms of the fraction X of settled planets. We do this by considering the effect of finite settlement civilization lifetimes on the steady states. We find a range of parameters for which 0 <  X < 1, i.e., the Galaxy supports a population of interstellar space-faring civilizations even though some settleable systems are uninhabited. In addition we find that statistical fluctuations can produce local overabundances of settleable worlds. These generate long-lived clusters of settled systems immersed in large regions that remain unsettled. Both results point to ways in which Earth might remain unvisited in the midst of an inhabited galaxy. Finally we consider how our results can be combined with the finite horizon for evidence of previous settlements in Earth’s geologic record. Using our steady-state model we constrain the probabilities for an Earth visit by a settling civilization before a given time horizon. These results break the link between Hart’s famous “Fact A” (no interstellar visitors on Earth now) and the conclusion that humans must, therefore, be the only technological civilization in the Galaxy. Explicitly, our solutions admit situations where our current circumstances are consistent with an otherwise settled, steady-state galaxy.

We present 114 trigonometric parallaxes for 107 nearby white dwarf (WD) systems from both the Cerro Tololo Inter-American Observatory Parallax Investigation (CTIOPI) and the U. S. Naval Observatory Flagstaff Station (NOFS) parallax programs. Of these, 76 parallaxes for 69 systems were measured by the CTIOPI program and 38 parallaxes for as many systems were measured by the NOFS program. A total of 50 systems are confirmed to be within the 25-pc horizon of interest. Coupled with a spectroscopic confirmation of a common proper-motion companion to a Hipparcos star within 25 pc as well as confirmation parallax determinations for two WD systems included in the recently released Tycho Gaia Astrometric Solution catalog, we add 53 new systems to the 25-pc WD sample—a 42% increase. Our sample presented here includes four strong candidate halo systems, a new metal-rich DAZ WD, a confirmation of a recently discovered nearby short-period ( P = 2.85 hr) double degenerate, a WD with a new astrometric perturbation (long period, unconstrained with our data), and a new triple system where the WD companion main-sequence star has an astrometric perturbation ( P ∼ 1.6 year).

Numerous research topics rely on an improved cosmic distance scale (e.g., cosmology, gravitational waves) and the NASA/IPAC Extragalactic Database of Distances (NED-D) supports those efforts by tabulating multiple redshift-independent distances for 12,000 galaxies (e.g., Large Magellanic Cloud (LMC) zero-point). Six methods for securing a mean estimate distance (MED) from the data are presented (e.g., indicator and Decision Tree). All six MEDs yield surprisingly consistent distances for the cases examined, including for the key benchmark LMC and M106 galaxies. The results underscore the utility of the NED-D MEDs in bolstering the cosmic distance scale and facilitating the identification of systematic trends.

The DESI Legacy Imaging Surveys ( http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing–Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg ² of the extragalactic sky visible from the northern hemisphere in three optical bands ( g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

ImageJ is a graphical user interface (GUI) driven, public domain, Java-based, software package for general image processing traditionally used mainly in life sciences fields. The image processing capabilities of ImageJ are useful and extendable to other scientific fields. Here we present AstroImageJ (AIJ), which provides an astronomy specific image display environment and tools for astronomy specific image calibration and data reduction. Although AIJ maintains the general purpose image processing capabilities of ImageJ, AIJ is streamlined for time-series differential photometry, light curve detrending and fitting, and light curve plotting, especially for applications requiring ultra-precise light curves (e.g., exoplanet transits). AIJ reads and writes standard Flexible Image Transport System (FITS) files, as well as other common image formats, provides FITS header viewing and editing, and is World Coordinate System aware, including an automated interface to the astrometry.net web portal for plate solving images. AIJ provides research grade image calibration and analysis tools with a GUI driven approach, and easily installed cross-platform compatibility. It enables new users, even at the level of undergraduate student, high school student, or amateur astronomer, to quickly start processing, modeling, and plotting astronomical image data with one tightly integrated software package.

The DESI Legacy Imaging Surveys ( http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing–Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg ² of the extragalactic sky visible from the northern hemisphere in three optical bands ( g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

We describe the catalogs assembled and the algorithms used to populate the revised TESS Input Catalog (TIC), based on the incorporation of the Gaia second data release. We also describe a revised ranking system for prioritizing stars for 2 minute cadence observations, and we assemble a revised Candidate Target List (CTL) using that ranking. The TIC is available on the Mikulski Archive for Space Telescopes server, and an enhanced CTL is available through the Filtergraph data visualization portal system at  http://filtergraph.vanderbilt.edu/tess_ctl.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey. This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semimajor axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M _* > 1.5 M _⊙ more likely to host planets with masses between 2 and 13 M _Jup and semimajor axes of 3–100 au at 99.92% confidence. We fit a double power-law model in planet mass ( m) and semimajor axis ( a) for planet populations around high-mass stars ( M _* > 1.5 M _⊙) of the form , finding α = −2.4 ± 0.8 and β = −2.0 ± 0.5, and an integrated occurrence rate of % between 5–13 M _Jup and 10–100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with % of stars hosting a brown dwarf companion between 13–80 M _Jup and 10–100 au. Brown dwarfs also appear to be distributed differently in mass and semimajor axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semimajor axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the radial velocity method, our results are consistent with a peak in occurrence of giant planets between ∼1 and 10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

We determine the orbital eccentricities of individual small Kepler planets, through a combination of asteroseismology and transit light-curve analysis. We are able to constrain the eccentricities of 51 systems with a single transiting planet, which supplement our previous measurements of 66 planets in multi-planet systems. Through a Bayesian hierarchical analysis, we find evidence that systems with only one detected transiting planet have a different eccentricity distribution than systems with multiple detected transiting planets. The eccentricity distribution of the single-transiting systems is well described by the positive half of a zero-mean Gaussian distribution with a dispersion σ _e  = 0.32 ± 0.06, while the multiple-transit systems are consistent with . A mixture model suggests a fraction of of single-transiting systems have a moderate eccentricity, represented by a Rayleigh distribution that peaks at . This finding may reflect differences in the formation pathways of systems with different numbers of transiting planets. We investigate the possibility that eccentricities are self-excited in closely packed planetary systems, as well as the influence of long-period giant companion planets. We find that both mechanisms can qualitatively explain the observations. We do not find any evidence for a correlation between eccentricity and stellar metallicity, as has been seen for giant planets. Neither do we find any evidence that orbital eccentricity is linked to the detection of a companion star. Along with this paper, we make available all of the parameters and uncertainties in the eccentricity distributions, as well as the properties of individual systems, for use in future studies.

We characterize the occurrence rate of planets, ranging in size from 0.5 to 16 R _⊕, orbiting FGK stars with orbital periods from 0.5 to 500 days. Our analysis is based on results from the “DR25” catalog of planet candidates produced by NASA’s Kepler mission and stellar radii from Gaia “DR2.” We incorporate additional Kepler data products to accurately characterize the efficiency of planets being recognized as “threshold crossing events” by Kepler’s Transiting Planet Search pipeline and labeled as planet candidates by the robovetter. Using a hierarchical Bayesian model, we derive planet occurrence rates for a wide range of planet sizes and orbital periods. For planets with sizes 0.75–1.5 R _⊕ and orbital periods of 237–500 days, we find a rate of planets per FGK star of <0.27 (84.13th percentile). While the true rate of such planets could be lower by a factor of ∼2 (primarily due to potential contamination of planet candidates by false alarms), the upper limits on the occurrence rate of such planets are robust to ∼10%. We recommend that mission concepts aiming to characterize potentially rocky planets in or near the habitable zone of Sun-like stars prepare compelling science programs that would be robust for a true rate in the range f _R,P  = 0.03–0.40 for 0.75–1.5 R _⊕ planets with orbital periods in 237–500 days, or a differential rate of 0.06–0.76.

Analysis of new precision radial velocity (RV) measurements from the Lick Automated Planet Finder and Keck HIRES has yielded the discovery of three new exoplanet candidates orbiting the nearby stars HD 190007 and HD 216520. We also report new velocities from the APF and the Planet Finder Spectrograph and updated orbital fits for the known exoplanet host stars GJ 686 and HD 180617. Of the newly discovered planets, HD 190007 b has a period of P = 11.72 days, an RV semiamplitude of K = 5.64 ± 0.55 m s ⁻¹, a minimum mass of M _pl = 16.46 ± 1.66 M _⊕, and orbits the slightly metal-rich, active K4V star HD 190007. For HD 216520 b, we find P = 35.45 days, K = 2.28 ± 0.20 m s ⁻¹, and M _pl = 10.26 ± 0.99 M _⊕, while for HD 216520 c, P = 154.43 days, K = 1.29 ± 0.22 m s ⁻¹, and M _pl = 9.44 ± 1.63 M _⊕. Both planets orbit the slightly metal-poor, inactive K0V star HD 216520. Our updated best-fit models for HD 180617 b and GJ 686 b are in good agreement with the published results. For HD 180617 b, we obtain P = 105.91 days and M _pl = 12.214 ± 1.05 M _⊕. For GJ 686 b, we find P = 15.53 days and M _pl = 6.624 ± 0.432 M _⊕. Using an injection-recovery exercise, we find that HD 190007 b and HD 216520 b are unlikely to have additional planets with masses and orbital periods within a factor of 2, in marked contrast to ∼85% of planets in this mass and period range discovered by Kepler.

The legacy imaging surveys for the Dark Energy Spectroscopic Instrument project provides multiple-color photometric data, which are about 2 mag deeper than those from the SDSS. In this study, we redetermine the fundamental properties for an old halo globular cluster of Palomar 5 based on these new imaging data, including structure parameters, stellar population parameters, and luminosity and mass functions. These characteristics, together with its tidal tails, are key for dynamical studies of the cluster and constraining the mass model of the Milky Way. By fitting the King model to the radial surface density profile of Palomar 5, we derive the core radius of , tidal radius of , and concentration parameter of c = 0.78 ± 0.04. We apply a Bayesian analysis method to derive the stellar population properties and get an age of 11.508 ± 0.027 Gyr, metallicity of [Fe/H]  =   −1.798 ± 0.014, reddening of E( B −  V) = 0.0552 ± 0.0005, and distance modulus of . The main-sequence luminosity and mass functions for both the cluster center and tidal tails are investigated. The luminosity and mass functions at different distances from the cluster center suggest that there is obvious spatial mass segregation. Many faint low-mass stars have been evaporated at the cluster center, and the tidal tails are enhanced by low-mass stars. Both the concentration and relaxation times suggest that Palomar 5 is a totally relaxed system.

The Voyager 1 (V1) and Voyager 2 (V2) spacecraft were launched in 1977 on a mission to explore the outer planets and reach the heliopause, the boundary between the hot solar plasma and the relatively cool interstellar plasma. V1 reached the heliopause on 2012 August 25, at 121.6 au, and V2 reached the heliopause on 2018 November 5, at 119.0 au. One of their remarkable discoveries was the detection of shocks propagating into the interstellar plasma from energetic solar events. These shocks are typically preceded by electron plasma oscillations excited by electron beams streaming along interstellar magnetic field lines ahead of the shocks. The frequencies of the plasma oscillations have now provided radial electron density profiles in the outer heliosphere and in the interstellar medium to radial distances of more than 145 au. The oscillations are typically preceded by bursts of high-energy ∼5–100 MeV electrons. These electron bursts are interpreted as being due to the reflection (and acceleration) of cosmic-ray electrons by the shock at the time the shock first contacts the magnetic field line that passes through the spacecraft. Relative timing between the cosmic rays reflected by the shock and the onset of the plasma oscillations allow us, for the first time, to estimate the energy, ∼20–100 eV, of the electron beams responsible for the plasma oscillations. These observations are combined into a self-consistent model called the foreshock model that describes the interaction of shocks of solar origin with the interstellar plasma.

High-resolution spectroscopy has opened the way for new, detailed study of exoplanet atmospheres. There is evidence that this technique can be sensitive to the complex, three-dimensional (3D) atmospheric structure of these planets. In this work, we perform cross-correlation analysis of high-resolution ( R ∼ 100,000) CRIRES/VLT emission spectra of the hot Jupiter HD 209458b. We generate template emission spectra from a 3D atmospheric circulation model of the planet, accounting for temperature structure and atmospheric motions—winds and planetary rotation—missed by spectra calculated from one-dimensional models. In this first-of-its-kind analysis, we find that using template spectra generated from a 3D model produces a more significant detection (6.9 σ) of the planet’s signal than any of the hundreds of one-dimensional models we tested (maximum of 5.1 σ). We recover the planet’s thermal emission, its orbital motion, and the presence of CO in its atmosphere at high significance. Additionally, we analyzed the relative influences of 3D temperature and chemical structures in this improved detection, including the contributions from CO and H ₂O, as well as the role of atmospheric Doppler signatures from winds and rotation. This work shows that the hot Jupiter’s 3D atmospheric structure has a first-order influence on its emission spectra at high resolution and motivates the use of multidimensional atmospheric models in high-resolution spectral analysis.

Planets around young stars trace the early evolution of planetary systems. We report the discovery and validation of two planetary systems with ages ≲300 Myr from observations by the Transiting Exoplanet Survey Satellite (TESS). The Myr old G star TOI-251 hosts a mini-Neptune with a day period. The Myr old K star TOI-942 hosts a system of inflated Neptune-sized planets, with TOI-942b orbiting in a period of days with a radius of and TOI-942c orbiting in a period of days with a radius of . Though we cannot place either host star into a known stellar association or cluster, we can estimate their ages via their photometric and spectroscopic properties. Both stars exhibit significant photometric variability due to spot modulation, with measured rotation periods of ∼3.5 days. These stars also exhibit significant chromospheric activity, with age estimates from the chromospheric calcium emission lines and X-ray fluxes matching that estimated from gyrochronology. Both stars also exhibit significant lithium absorption, similar in equivalent width to well-characterized young cluster members. TESS has the potential to deliver a population of young planet-bearing field stars, contributing significantly to tracing the properties of planets as a function of their age.

State-of-the-art atomic and optical clocks have the great potential to precisely test fundamental physical assumptions and enhance our understanding of nature. Their widespread applications require us to rigorously deduce the relativistic frequency shift in the framework of general relativity. One interesting question for clocks is tidal field effect which contains variously periodic variations. By introducing tidally deformed Earth, the general algorithms for clock and frequency comparisons are derived. We investigate the effects of external gravitational bodies and Earth’s tidal deformation on the ground- and space-based clocks. The orbital elements of satellite and locations of laboratory are introduced for corresponding clocks, in which the secular and long-period terms and short-period terms are subsequently presented. These effects called tidal clock effects produce the non-negligible contributions in the modern clock experiments and can be directly evaluated from our parameterized formulas. In addition, we also demonstrate for tidal clock effects the position and distance dependences, as well as periodicity. These tidal effects and dependences can provide valuable information for the clock comparison experiments.

The Kepler DR25 planet candidate catalog was produced using an automated method of planet candidate identification based on various tests. These tests were tuned to obtain a reasonable but arbitrary balance between catalog completeness and reliability. We produce new catalogs with differing balances of completeness and reliability by varying these tests, and study the impact of these alternative catalogs on occurrence rates. We find that if there is no correction for reliability, different catalogs give statistically inconsistent occurrence rates, while if we correct for both completeness and reliability, we get statistically consistent occurrence rates. This is a strong indication that correction for completeness and reliability is critical for the accurate computation of occurrence rates. Additionally, we find that this result is the same whether using Bayesian Poisson-likelihood Markov Chain Monte Carlo or Approximate Bayesian Computation methods. We also examine the use of a Robovetter disposition score cut as an alternative to reliability correction, and find that while a score cut does increase the reliability of the catalog, it is not as accurate as performing a full reliability correction. We get the same result when performing a reliability correction with and without a score cut. Therefore removing low-score planets removes data without providing any advantage, and should be avoided when possible. We make our alternative catalogs publicly available, and propose that these should be used as a test of occurrence rate methods, with the requirement that a method should provide statistically consistent occurrence rates for all these catalogs.

Understanding the atmospheres of exoplanets is a milestone to decipher their formation history and potential habitability. High-contrast imaging and spectroscopy of exoplanets is the major pathway toward the goal. Direct imaging of an exoplanet requires high spatial resolution. Interferometry has proven to be an effective way of improving spatial resolution. However, means of combining interferometry, high-contrast imaging, and high-resolution spectroscopy have been rarely explored. To fill in the gap, we present the dual-aperture fiber nuller (FN) for current-generation 8–10 m telescopes, which provides the necessary spatial and spectral resolution to (1) conduct follow-up spectroscopy of known exoplanets and (2) detect planets in debris-disk systems. The concept of feeding an FN to a high-resolution spectrograph can also be used for future space and ground-based missions. We present a case study of using the dual-aperture FN to search for biosignatures in rocky planets around M stars for a future space interferometry mission. Moreover, we discuss how an FN can be equipped on future extremely large telescopes by using the Giant Magellan Telescope as an example.

We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6 cm and at 7 mm at resolutions of 0.″15 (1650 au) and 0.″04 (440 au), respectively. The 3.6 cm data reveal new morphological detail in the ultracompact H ii region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6 cm and 7 mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2 mm. Three of the smallest ultracompact H ii regions (sources A, B2, and G2) all show broad kinematic linewidths, with Δ V _FWHM ≳ 40 km s ⁻¹. A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 yr time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7 ± 1.8 km s ⁻¹ to 55.6 ± 2.7 km s ⁻¹, an increase of +23.9 ± 3.4 km s ⁻¹. The G2 source was previously reported to have shown a 3.6 cm continuum flux density decrease of 40% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes.

Numerous research topics rely on an improved cosmic distance scale (e.g., cosmology, gravitational waves) and the NASA/IPAC Extragalactic Database of Distances (NED-D) supports those efforts by tabulating multiple redshift-independent distances for 12,000 galaxies (e.g., Large Magellanic Cloud (LMC) zero-point). Six methods for securing a mean estimate distance (MED) from the data are presented (e.g., indicator and Decision Tree). All six MEDs yield surprisingly consistent distances for the cases examined, including for the key benchmark LMC and M106 galaxies. The results underscore the utility of the NED-D MEDs in bolstering the cosmic distance scale and facilitating the identification of systematic trends.

Stratospheric Observatory for Infrared Astronomy High-resolution Airborne Wideband Camera Plus polarimetry at 154 μm is reported for the face-on galaxy M51 and the edge-on galaxy NGC 891. For M51, the polarization vectors generally follow the spiral pattern defined by the molecular gas distribution, the far-infrared (FIR) intensity contours, and other tracers of star formation. The fractional polarization is much lower in the FIR-bright central regions than in the outer regions, and we rule out loss of grain alignment and variations in magnetic field strength as causes. When compared with existing synchrotron observations, which sample different regions with different weighting, we find the net position angles are strongly correlated, the fractional polarizations are moderately correlated, but the polarized intensities are uncorrelated. We argue that the low fractional polarization in the central regions must be due to significant numbers of highly turbulent segments across the beam and along lines of sight in the beam in the central 3 kpc of M51. For NGC 891, the FIR polarization vectors within an intensity contour of 1500 are oriented very close to the plane of the galaxy. The FIR polarimetry is probably sampling the magnetic field geometry in NGC 891 much deeper into the disk than is possible with NIR polarimetry and radio synchrotron measurements. In some locations in NGC 891, the FIR polarization is very low, suggesting we are preferentially viewing the magnetic field mostly along the line of sight, down the length of embedded spiral arms. There is tentative evidence for a vertical field in the polarized emission off the plane of the disk.

The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for 27 contiguous days in a year, making it difficult to measure rotation periods with traditional methods. This is especially problematic for field M dwarfs, which are ideal candidates for exoplanet searches, but which tend to have periods in excess of the 27 day observing baseline. We present a new tool, Astraea, for predicting long rotation periods from short-duration light curves combined with stellar parameters from Gaia DR2. Using Astraea, we can predict the rotation periods from Kepler 4 yr light curves with 13% uncertainty overall (and a 9% uncertainty for periods >30 days). By training on 27 day Kepler light-curve segments, Astraea can predict rotation periods up to 150 days with 9% uncertainty (5% for periods >30 days). After training this tool on these 27 day Kepler light-curve segments, we applied Astraea to real TESS data. For the 195 stars that were observed by both Kepler and TESS, we were able to predict the rotation periods with 55% uncertainty despite the wild differences in systematics.

We present the discoveries of KELT-25 b (TIC 65412605, TOI-626.01) and KELT-26 b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A stars. The transit signals were initially detected by the KELT survey and subsequently confirmed by Transiting Exoplanet Survey Satellite (TESS) photometry. KELT-25 b is on a 4.40 day orbit around the V = 9.66 star CD-24 5016 ( K, M _⋆ =  M _⊙), while KELT-26 b is on a 3.34 day orbit around the V = 9.95 star HD 134004 ( = K, M _⋆ = M _⊙), which is likely an Am star. We have confirmed the substellar nature of both companions through detailed characterization of each system using ground-based and TESS photometry, radial velocity measurements, Doppler tomography, and high-resolution imaging. For KELT-25, we determine a companion radius of R _P = R _J and a 3 σ upper limit on the companion’s mass of ∼64 M _J. For KELT-26 b, we infer a planetary mass and radius of M _P = and R _P = R _J. From Doppler tomographic observations, we find KELT-26 b to reside in a highly misaligned orbit. This conclusion is weakly corroborated by a subtle asymmetry in the transit light curve from the TESS data. KELT-25 b appears to be in a well-aligned, prograde orbit, and the system is likely a member of the cluster Theia 449.