Table of contents

The exoplanet HD 118203 b, orbiting a bright (V = 8.05) host star, was discovered using the radial velocity method by da Silva et al., but was not previously known to transit. Transiting Exoplanet Survey Satellite (TESS) photometry has revealed that this planet transits its host star. Nine planetary transits were observed by TESS, allowing us to measure the radius of the planet to be ${1.136}_{-0.028}^{+0.029}{R}_{J}$, and to calculate the planet mass to be ${2.166}_{-0.079}^{+0.074}{M}_{J}$. The host star is slightly evolved with an effective temperature of ${T}_{\mathrm{eff}}={5683}_{-85}^{+84}$ K and a surface gravity of $\mathrm{log}\,g={3.889}_{0.018}^{0.017}$. With an orbital period of ${6.134985}_{-0.000030}^{+0.000029}$ days and an eccentricity of 0.314 ± 0.017, the planet occupies a transitional regime between circularized hot Jupiters and more dynamically active planets at longer orbital periods. The host star is among the 10 brightest known to have transiting giant planets, providing opportunities for both planetary atmospheric and asteroseismic studies.

We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager (GPI). HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region (∼3 stars per square arcsecond with H < 22), necessary because of the small field of view of the GPI. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (ρ < 2''). We achieved a proper motion precision of ∼0.3 mas yr⁻¹ and constrained the parallax of each star to be ≲1 mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements, we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with excellent astrometric calibration. We find that the predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers.

The discovery of spectral type transition of active galactic nuclei (AGNs), the so-called "changing-look" (CL) phenomenon, challenges the widely accepted AGN paradigm, not only in the orientation-based unified model, but also in the standard disk model. In past decades, only a couple of nearby repeat changing-look active galactic nuclei (CL-AGNs) have been identified. Here we report spectroscopic observations of UGC 3223 over the course of 18 yr, from 2001 onwards. Combining the spectrum taken in 1987 by Stirpe, we have witnessed its type transitions from $1.5\to 2.0\to 1.8$ over 32 yr, and captured a long-lived (at least 10 yr) thorough "turn-off" state with a spectrum typical of a Seyfert 2 galaxy. The long-term thorough turn-off state probably suggests a once-dormant and an awakening central engine in UGC 3223. We argue that the (dis)appearance of the broad Balmer emission lines can be explained by the disk–wind broad-line region model given the evolution of the calculated Eddington ratio of accretion of the supermassive black hole.

We present UBVRI photometry of 8702 stars in a 05 × 05 field in the direction of NGC 7142, taken with the Half Degree Imager at the WIYN 0.9 m telescope, to improve knowledge of this cluster's basic parameters. Our photometry spans the ranges 10.6–20.4 mag in U, 10.6–22.0 mag in B, 10.0–21.8 mag in V, 9.2–20.7 mag in R, and 8.5–19.9 mag in I. Using color–color diagrams that employ all four color combinations that include U, versus B − V, we derive a reddening–metallicity relation for the cluster, with preferred values E(B − V) = 0.338 ± 0.031 mag for the left-edge fiducial of the main sequence and [Fe/H] = 0.0 ± 0.1 dex, where the Hyades cluster has been used as an unreddened reference cluster, the extinction relations of Cardelli have been employed, and the metallicity dependence of the Yonsei-Yale (Y²) isochrones has been assumed. Comparison of our data to the Y² isochrones in multiple color–magnitude diagrams (CMDs) yields distance–metallicity and age–metallicity relations, with preferred values of $m-M=12.65\pm 0.23$ mag and age = ${4.0}_{+1.3}^{-0.7}$ Gyr. Re-evaluation of the parameters of M67 using Stetson's UBVI photometry yields [Fe/H] = −0.02 ± 0.05 dex, E(B − V) = 0.04 ± 0.01 mag, m − M = 9.75 ± 0.03 mag, and age = 3.85 ± 0.17 Gyr; we thus find the metallicity and age of the two clusters to be indistinguishable. A semi-independent analysis adopting the parameters of M67 and shifting the fiducials of the two clusters in six CMDs until they match strongly corroborates the values listed above. The differences between our inferred parallaxes and the Gaia DR2 values are 87 ± 60 μas for NGC 7142 and 48 ± 15 μas for M67, consistent with previous studies.

We investigate a possible correlation between the solid surface density Σ of the minimum-mass extrasolar nebula (MMEN) and the host star mass M_⋆ and metallicity [Fe/H]. Leveraging on the precise host star properties from the California-Kepler Survey (CKS), we found that ${\rm{\Sigma }}={50}_{-20}^{+33}\,{\rm{g}}\,{\mathrm{cm}}^{-2}$ (a/1 au)^{−1.75±0.07} (M_⋆/M_⊙)^1.04±0.22 10^{0.22±0.05[Fe/H]} for Kepler-like systems (1–4R_⊕; a < 1 au). The strong M_⋆ dependence is reminiscent of previous dust continuum results that the solid disk mass scales with M_⋆. The weaker [Fe/H] dependence shows that sub-Neptune planets, unlike giant planets, form readily in lower metallicity environment. The innermost region (a < 0.1 au) of an MMEN maintains a smooth profile despite a steep decline of planet occurrence rate: a result that favors the truncation of disks by corotating magnetospheres with a range of rotation periods, rather than the sublimation of dust. The Σ of Kepler multitransiting systems shows a much stronger correlation with M_⋆ and [Fe/H] than singles. This suggests that the dynamically hot evolution that produced single systems also partially removed the memory of formation in disks. Radial-velocity planets yielded a MMEN very similar to CKS planets; transit-timing-variation planets' postulated convergent migration history is supported by their poorly constrained MMEN. We found that lower mass stars have a higher efficiency of forming/retaining planets: for Sun-like stars, about 20% of the solid mass within ∼1 au are converted/preserved as sub-Neptunes, compared to 70% for late-K to early-M stars. This may be due to the lower binary fraction, lower giant-planet occurrence, or the longer disk lifetime of lower mass stars.

We present exoplanet occurrence rates estimated with approximate Bayesian computation for planets with radii between 0.5 and 16 R_⊕ and orbital periods between 0.78 and 400 days orbiting FGK dwarf stars. We base our results on an independent planet catalog compiled from our search of all ∼200,000 stars observed over the Kepler mission, with precise planetary radii supplemented by Gaia DR2-incorporated stellar radii. We take into account detection and vetting efficiency, planet radius uncertainty, and reliability against transit-like noise signals in the data. By analyzing our FGK occurrence rates as well as those computed after separating F-, G-, and K-type stars, we explore dependencies on stellar effective temperature, planet radius, and orbital period. We reveal new characteristics of the photoevaporation-driven "radius gap" between ∼1.5 and 2 R_⊕, indicating that the bimodal distribution previously revealed for P < 100 days exists only over a much narrower range of orbital periods, above which sub-Neptunes dominate and below which super-Earths dominate. Finally, we provide several estimates of the "eta-Earth" value—the frequency of potentially habitable, rocky planets orbiting Sun-like stars. For planets with sizes 0.75–1.5 R_⊕ orbiting in a conservatively defined habitable zone (0.99–1.70 au) around G-type stars, we place an upper limit (84.1th percentile) of <0.18 planets per star.

The abrupt brightening of an Energetic Neutral Atom (ENA) blob or cloud has been interpreted as plasma injection in the Kronian magnetosphere (termed ENA injection herein). Morphologically, there appears to be two types of abrupt ENA cloud brightening: (1) the brightening of a large cloud usually seen at distances >10–12 R_s (R_s ∼ 60,268 km) in the midnight or postmidnight region; (2) the brightening of a smaller cloud usually seen at distances <10–12 R_s around 21–03 magnetic local time. Among many radio waves observed at Saturn, type 2 ENA injections correlate best with the 5 kHz narrowband (NB) waves. Using Cassini Ion and Neutral Camera (INCA) and Radio and Plasma Wave Science (RPWS) data, we examine the periodicities of the type 2 ENA injections and the 5 kHz NB emissions as well as their cross-correlations, which have been previously used to measure the lag times or phase differences. Because correlational analysis can only establish linear relationships, we also use mutual information to establish linear and nonlinear relationships. On average, the peak of the 5 kHz NB emission lags those of the type 2 ENA injection by about a few minutes to 2 hr. The injection of hot plasma to the inner magnetosphere can lead to temperature anisotropy, which can generate electrostatic upper hybrid waves, which upon encountering the high-density gradient at the outer edge of the Enceladus density torus, can mode convert to Z mode and then to O mode. The 5 kHz NB waves commonly propagate in the O mode.

Photometric and astrometric calibration of high-contrast images is essential for the characterization of companions at small angular separation from their stellar host. The main challenge to performing accurate relative photometry and astrometry of high-contrast companions with respect to the host star is that the central starlight cannot be directly used as a reference, as it is either blocked by a coronagraphic mask or saturating the detector. Our approach is to add fiducial incoherent faint copies of the host star in the image plane and alternate the pattern of these copies between exposures. Subtracting two frames with different calibration patterns removes measurement bias due to static and slowly varying incoherent speckle halo components, while ensuring that calibration references are inserted on each frame. Each calibration pattern is achieved by high-speed modulation of a pupil-plane deformable mirror to ensure incoherence. We implemented the technique on-sky on the Subaru Coronagraphic Extreme Adaptive Optics instrument with speckles which were of the order of 10³ times fainter than the central host. The achieved relative photometric and astrometric measurement precisions for 10 s exposure were respectively 5% and 20 milliarcsecond. We also demonstrate, over a 540 s measurement span, that residual photometric and astrometric errors are uncorrelated in time, indicating that residual noise averages as the inverse square root of the number of exposures in longer time-series data sets.

We present new H-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager. The disk is detected in total and polarized intensity down to a projected angular separation of 015, or 20 au. On the other hand, the large-scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of ≈100 au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high-fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

Scattered light imaging has revealed nearly a dozen circumstellar disks around young Herbig Ae/Be stars—enabling studies of structures in the upper disk layers as potential signs of ongoing planet formation. We present the first images of the disk around the variable Herbig Ae star PDS 201 (V* V351 Ori) and an analysis of the images and spectral energy distribution through 3D Monte Carlo radiative transfer simulations and forward modeling. The disk is detected in three data sets with the Large Binocular Telescope Interferometer/Large Binocular Telescope mid-infrared camera at the Large Binocular Telescope, including direct observations in the Ks and L' filters, and an L' observation with the 360° vector apodizing phase plate coronagraph. The scattered light disk extends to a very large radius of ∼250 au, which places it among the largest of such disks. Exterior to the disk, we establish detection limits on substellar companions down to ∼5 M_Jup at ≳15 (≳500 au), assuming the Baraffe et al. models. The images show a radial gap extending to ∼04 (∼140 au at a distance of 340 pc) that is also evident in the spectral energy distribution. The large gap is a possible signpost of multiple high-mass giant planets at orbital distances (∼60–100 au) that are unusually massive and widely separated compared to those of planet populations previously inferred from protoplanetary disk substructures.

We report the detection of the first circumbinary planet (CBP) found by Transiting Exoplanet Survey Satellite (TESS). The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30 minute cadence and in sectors 4 through 12 at 2 minute cadence. It consists of two stars with masses of 1.1 M_⊙ and 0.3 M_⊙ on a slightly eccentric (0.16), 14.6 day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ∼6.9 R_⊕ and was observed to make three transits across the primary star of roughly equal depths (∼0.2%) but different durations—a common signature of transiting CBPs. Its orbit is nearly circular (e ≈ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ∼1°. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for CBPs and provides further understanding of the formation and evolution of planets orbiting close binary stars.

The most massive and complex globular clusters in the Galaxy are thought to have originated as the nuclear cores of now tidally disrupted dwarf galaxies, but the connection between globular clusters and dwarf galaxies is tenuous with the M54/Sagittarius system representing the only unambiguous link. The globular cluster Omega Centauri (ω Cen) is more massive and chemically diverse than M 54, and is thought to have been the nuclear star cluster of either the Sequoia or Gaia-Enceladus galaxy. Local Group dwarf galaxies with masses equivalent to these systems often host significant populations of very metal-poor stars ([Fe/H] < −2.5), and one might expect to find such objects in ω Cen. Using high-resolution spectra from Magellan-M2FS, we detected 11 stars in a targeted sample of 395 that have [Fe/H] ranging from −2.30 to −2.52. These are the most metal-poor stars discovered in the cluster, and are five times more metal-poor than ω Cen's dominant population. However, these stars are not so metal-poor as to be unambiguously linked to a dwarf galaxy origin. The cluster's metal-poor tail appears to contain two populations near [Fe/H] ∼ −2.1 and −2.4, which are very centrally concentrated but do not exhibit any peculiar kinematic signatures. Several possible origins for these stars are discussed.

We report the discovery of two transiting exoplanets from the WASP survey, WASP-150b and WASP-176b. WASP-150b is an eccentric (e = 0.38) hot Jupiter on a 5.6 day orbit around a V = 12.03, F8 main-sequence host. The host star has a mass and radius of 1.4 ${M}_{\odot }$ and 1.7 ${R}_{\odot }$ respectively. WASP-150b has a mass and radius of 8.5 ${M}_{{\rm{J}}}$ and 1.1 R_J, leading to a large planetary bulk density of 6.4 ρ_J. WASP-150b is found to be ∼3 Gyr old, well below its circularization timescale, supporting the eccentric nature of the planet. WASP-176b is a hot Jupiter planet on a 3.9 day orbit around a V = 12.01, F9 sub-giant host. The host star has a mass and radius of 1.3 M_⊙ and 1.9 R_⊙. WASP-176b has a mass and radius of 0.86 M_J and 1.5 R_J, respectively, leading to a planetary bulk density of 0.23 ρ_J.

We report the analysis of the microlensing event OGLE-2018-BLG-0677. A small feature in the light curve of the event leads to the discovery that the lens is a star–planet system. Although there are two degenerate solutions that could not be distinguished for this event, both lead to a similar planet-host mass ratio. We perform a Bayesian analysis based on a Galactic model to obtain the properties of the system and find that the planet corresponds to a super-Earth/sub-Neptune with a mass of ${M}_{\mathrm{planet}}={3.96}_{-2.66}^{+5.88}{M}_{\oplus }$. The host star has a mass of ${M}_{\mathrm{host}}={0.12}_{-0.08}^{+0.14}\,{M}_{\odot }$. The projected separation for the inner and outer solutions are ${0.63}_{-0.17}^{+0.20}$ au and ${0.72}_{-0.19}^{+0.23}$ au respectively. At ${\rm{\Delta }}{\chi }^{2}={\chi }^{2}(1{\rm{L}}1{\rm{S}})-{\chi }^{2}(2{\rm{L}}1{\rm{S}})=46$, this is by far the lowest Δχ² for any securely detected microlensing planet to date, a feature that is closely connected to the fact that it is detected primarily via a "dip" rather than a "bump."

We present parallaxes, proper motions, and J-band photometry for 348 L and T dwarfs measured using the wide-field near-infrared camera WFCAM on the United Kingdom Infrared Telescope. This is the largest single batch of infrared parallaxes for brown dwarfs to date. Our parallaxes have a median uncertainty of 3.5 mas, similar to most previous ground-based infrared parallax surveys. Our target list was designed to complete a volume-limited parallax sample of L0–T8 dwarfs out to 25 pc spanning declinations −30° to +60° (68% of the sky). We report the first parallaxes for 165 objects, and we improve on previous measurements for another 53 objects. Our targets include 104 objects (mostly early-L dwarfs) having Gaia DR2 parallax measurements with which our parallaxes are consistent. We include an extensive comparison of previous literature parallaxes for L and T dwarfs with both our results and Gaia DR2 measurements, identifying systematic offsets for some previous surveys. Our parallaxes confirm that 14 objects previously identified as wide common proper motion companions to main-sequence stars have distances consistent with companionship. We also report new ${J}_{\mathrm{MKO}}$ photometry for our targets, including the first measurements for 193 of our targets and improvements over previously published ${J}_{\mathrm{MKO}}$ photometry for another 60 targets. Altogether, our parallaxes will enable the first population studies using a volume-limited sample spanning spectral types L0–T8 defined entirely by parallaxes.

Accurate photometric redshift (photo-z) estimates are essential to the cosmological science goals of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this work, we use simulated photometry for mock galaxy catalogs to explore how LSST photo-z estimates can be improved by the addition of near-infrared (NIR) and/or ultraviolet (UV) photometry from the Euclid, Wide-Field InfrarRed Survey Telescope (WFIRST), and/or Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR) space telescopes. Generally, we find that deeper optical photometry can reduce the standard deviation of the photo-z estimates more than adding NIR or UV filters but that additional filters are the only way to significantly lower the fraction of galaxies with catastrophically under- or overestimated photo-z. For Euclid, we find that the addition of JH 5σ photometric detections can reduce the standard deviation for galaxies with z > 1 (z > 0.3) by ∼20% (∼10%), and the fraction of outliers by ∼40% (∼25%). For WFIRST, we show how the addition of deep YJHK photometry could reduce the standard deviation by ≳50% at z > 1.5 and drastically reduce the fraction of outliers to just ∼2% overall. For CASTOR, we find that the addition of its UV- and u-band photometry could reduce the standard deviation by ∼30% and the fraction of outliers by ∼50% for galaxies with z < 0.5. We also evaluate the photo-z results within sky areas that overlap with both the NIR and UV surveys and when spectroscopic training sets built from the surveys' small-area deep fields are used.

We monitored the z = 0.158 quasar 3C 273 between 2015 and 2019 in the optical (BVrz) and near-infrared (JHK) with the aim to perform dust reverberation mapping. Accounting for host galaxy and accretion disk contributions, we obtained pure dust light curves in JHK. Cross correlations between the V band and the dust light curves yield an average rest-frame delay for the hot dust of τ_cent ∼ 410 days. This is a factor of two shorter than that expected from the the dust ring radius R_x ∼ 900 lt-day reported from interferometric studies. The dust covering factor (CF) is about 8%, much smaller than that predicted from the half covering angle of 45° found for active galactic nuclei (AGNs). We analyze the asymmetric shape of the correlation functions and explore whether an inclined biconical bowl-shaped dust torus geometry could bring these findings (τ_cent, R_x and CF) into a consistent picture. The hot varying dust emission originates from the edge of the bowl rim with a small covering angle 40° < θ < 45°, and we see only the near side of the biconus. Such a dust gloriole with R_x = 900 ± 200 lt-day and an inclination 12° matches the data remarkably well. Comparing the results of 3C 273 with literature for less luminous AGN, we find a lag–luminosity relation τ ∝ L^α with α = 0.33–0.40, flatter than the widely adopted relation with α ∼ 0.5. We address several explanations for the new lag–luminosity relation.

An intriguing question in the context of dynamics arises: could a moon possess a moon itself? Such a configuration does not exist in the solar system, although this may be possible in theory. Kollmeier & Raymond determined the critical size of a satellite necessary to host a long-lived subsatellite, or submoon. However, the orbital constraints for these submoons to exist are still undetermined. Domingos et al. indicated that moons are stable out to a fraction of the host planet's Hill radius R_H,p, which in turn depend on the eccentricity of its host's orbit. Motivated by this, we simulate systems of exomoons and submoons for 10⁵ planetary orbits, while considering many initial orbital phases to obtain the critical semimajor axis in terms of R_H,p or the host satellite's Hill radius R_H,sat, respectively. We find that, assuming circular coplanar orbits, the stability limit for an exomoon is 0.40 R_H,p and for a submoon is 0.33 R_H,sat. Additionally, we discuss the observational feasibility of detecting these subsatellites through photometric, radial velocity, or direct imaging observations using the Neptune-sized exomoon candidate Kepler 1625b-I and identify how stability can shape the identification of future candidates.

We present the discovery of a planet on a very wide orbit in the microlensing event OGLE-2012-BLG-0838. The signal of the planet is well separated from the main peak of the event and the planet–star projected separation is found to be twice the Einstein ring radius, which corresponds to a projected separation of ≈4 au. Similar planets around low-mass stars are very hard to find using any technique other than microlensing. We discuss microlensing model fitting in detail and discuss the prospects for measuring the mass and distance of the lens system directly.

High-cadence observations of the Galactic bulge by the microlensing surveys led to the discovery of a handful of extremely short-timescale microlensing events that can be attributed to free-floating or wide-orbit planets. Here, we report the discovery of another strong free-floating planet candidate, which was found from the analysis of the gravitational microlensing event OGLE-2019-BLG-0551. The light curve of the event is characterized by a very short duration (≲3 days) and a very small amplitude (≲0.1 mag). From modeling of the light curve, we find that the Einstein timescale, ${t}_{{\rm{E}}}=0.381\pm 0.017$ day, is much shorter, and the angular Einstein radius, ${\theta }_{{\rm{E}}}=4.35\pm 0.34$μas, is much smaller than those of typical lensing events produced by stellar-mass lenses (${t}_{{\rm{E}}}\sim 20$ days, ${\theta }_{{\rm{E}}}\sim 0.3$ mas), indicating that the lens is very likely to be a planetary-mass object. We conduct an extensive search for possible signatures of a companion star in the light curve of the event, finding no significant evidence for the putative host star. For the first time, we also demonstrate that the angular Einstein radius of the lens does not depend on blending in the low-magnification events with strong finite source effects.

We present L'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wave front sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semimajor axis of ${20}_{-4}^{+3}$ au and PDS 70 c to have a semimajor axis of ${34}_{-6}^{+12}$ au (95% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2–3 R_Jup. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 M_Jup and a mean mass accretion rate between 3 × 10⁻⁷ and 8 × 10⁻⁷M_Jup/yr. For PDS 70 c, we computed a mass between 1 and 3 M_Jup and mean mass accretion rate between 1 × 10⁻⁷ and 5 × 10⁻⁷M_Jup/yr. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

In this work, we investigate the size, thermal inertia, surface roughness, and geometric albedo of 10 Vesta family asteroids using the Advanced Thermophysical Model, based on the thermal-infrared data acquired by mainly NASA's Wide-field Infrared Survey Explorer. Here, we show that the average thermal inertia and geometric albedo of the investigated Vesta family members are 42 J m⁻² s^−1/2 K⁻¹ and 0.314, respectively, where the derived effective diameters are less than 10 km. Moreover, the family members have a relatively low roughness fraction on their surfaces. The similarity in thermal inertia and geometric albedo among the V-type Vesta family members may reveal their close connection in origin and evolution. As the fragments of the cratering event of Vesta, the family members may have undergone a similar evolutionary process, thereby leading to very close thermal properties. Finally, we estimate their regolith grain sizes with different volume filling factors.

Architecture and parameters of two wide, nearby hierarchical systems containing five solar-type stars each, κ Tuc and ξ Sco, are studied. Using Gaia astrometry and photometry, masses are determined from visual orbits and isochrones, and effective temperatures from spectra or colors. Both systems are ∼2 Gyr old. Their spatial motion corresponds to a young disk but does not match any known kinematic group. Internal proper motions relative to the center of mass and radial velocities show that wide ∼8 kau outer pairs are bound. No correlation between orbit orientations in the inner subsystems is observed. All masses except one are confined to the narrow range from 0.8 to 1.5 solar. Strongly correlated masses and wide orbits can be explained if those systems formed by fragmentation in relative isolation and their components accreted gas from a common source, as expected in a hierarchical collapse. Young moving groups could be formed in similar environments, and many of them contain high-order hierarchies.

The orbits of eight systems with low-mass components (HIP 14524, HIP 16025, HIP 28671, HIP 46199, HIP 47791, HIP 60444, HIP 61100, and HIP 73085) are presented. Speckle interferometric data were obtained at the 6 m Big Telescope Alt-azimuth Special Astrophysical Observatory of the Russian Academy of Sciences (BTA SAO RAS) from 2007 to 2019. New data, together with measures already in the literature, made it possible to improve upon previous orbital solutions in six cases and to construct orbits for the first time in the two remaining cases (HIP 14524 and HIP 60444). Mass sums are obtained using both Hipparcos and Gaia parallaxes, and a comparison with previously published values is made. Using the Worley & Heintz criteria, the classification of the orbits constructed is carried out.

We report the discovery of HATS-71b, a transiting gas giant planet on a $P=3.7955$ day orbit around a $G=15.35$ mag M3 dwarf star. HATS-71 is the coolest M dwarf star known to host a hot Jupiter. The loss of light during transits is 4.7%, more than in any other confirmed transiting planet system. The planet was identified as a candidate by the ground-based HATSouth transit survey. It was confirmed using ground-based photometry, spectroscopy, and imaging, as well as space-based photometry from the NASA Transiting Exoplanet Survey Satellite mission (TIC 234523599). Combining all of these data, and utilizing Gaia DR2, we find that the planet has a radius of $1.024\pm 0.018$${R}_{{\rm{J}}}$ and mass of $0.37\,\pm 0.24$${M}_{{\rm{J}}}$ (95% confidence upper limit of $\lt 0.80$${M}_{{\rm{J}}}$), while the star has a mass of $0.4861\pm 0.0060$${M}_{\odot }$ and a radius of $0.4783\pm 0.0060$${R}_{\odot }$.

We present the details of the Bayesian analysis of the planetary microlensing event MOA-2016-BLG-227, whose excess flux is likely due to a source/lens companion or an unrelated ambient star, as well as of the assumed prior distributions. Furthermore, we apply this method to four reported planetary events, MOA-2008-BLG-310, MOA-2011-BLG-293, OGLE-2012-BLG-0527, and OGLE-2012-BLG-0950, where adaptive optics observations have detected excess flux at the source star positions. For events with small angular Einstein radii, our lens mass estimates are more uncertain than those of previous analyses, which assumed that the excess was due to the lens. Our predictions for MOA-2008-BLG-310 and OGLE-2012-BLG-0950 are consistent with recent results on these events obtained via Keck and Hubble Space Telescope observations when the source star is resolvable from the lens star. For events with small angular Einstein radii, we find that it is generally difficult to conclude whether the excess flux comes from the host star. Therefore, it is necessary to identify the lens star by measuring its proper motion relative to the source star to determine whether the excess flux comes from the lens star. Even without such measurements, our method can be used to statistically test the dependence of the planet-hosting probability on the stellar mass.

The diffuse ultraviolet background radiation has been mapped over most of the sky with 2' resolution using data from the Galaxy Evolution Explorer survey. We utilize this map to study the correlation between the UV background and clusters of galaxies discovered via the Sunyaev–Zeldovich effect in the Planck survey. We use only high Galactic latitude ($| b| \gt 60^\circ$) galaxy clusters to avoid contamination by Galactic foregrounds, and we only analyze clusters with a measured redshift. This leaves us with a sample of 142 clusters over the redshift range of 0.02 ≤ z ≤ 0.72, which we further subdivide into four redshift bins. In analyzing our stacked samples binned by redshift, we find evidence for a central excess of UV background light compared to local backgrounds for clusters with z < 0.3. We then stacked these z < 0.3 clusters to find a statistically significant excess of 12 ± 2.3 photon cm⁻² s⁻¹ sr⁻¹ Å⁻¹  over the median of ∼380 photon cm⁻² s⁻¹ sr⁻¹ Å⁻¹  measured around random blank fields. We measure the stacked radial profile of these clusters, and find that the excess UV radiation decays to the level of the background at a radius of ∼1 Mpc, roughly consistent with the maximum radial extent of the clusters. Analysis of possible physical processes contributing to the excess UV brightness indicates that non-thermal emission from relativistic electrons in the intracluster medium and faint, unresolved UV emission from cluster member galaxies and intracluster light are likely the dominant contributors.

The radial velocities of 2768 fundamental-mode RR Lyrae stars (RRLs) toward the southern Galactic bulge are presented, spanning the southern bulge from −8° < l < + 8° and −3° < b  < −6°. Distances derived from the pulsation properties of the RRLs are combined with Gaia proper motions to give constraints on the orbital motions of 1389 RRLs. The majority (∼75%) of the bulge RRLs have orbits consistent with these stars being permanently bound to <3.5 kpc from the Galactic Center, similar to the bar. However, unlike the bulge giants, the RRLs exhibit slower rotation and a higher velocity dispersion. The higher velocity dispersion arises almost exclusively from halo interlopers passing through the inner Galaxy. We present 82 stars with space velocities ≳500 km s⁻¹ and find that the majority of these high-velocity stars are halo interlopers; it is unclear if a subsample of these stars with similar space velocities has a common origin. Once the 25% of the sample represented by halo interlopers is cleaned, we can clearly discern two populations of bulge RRLs in the inner Galaxy. The first population of RRLs is not as tightly bound to the Galaxy (but is still confined to the inner ∼3.5 kpc) and is both spatially and kinematically consistent with the barred bulge. The second population is more centrally concentrated and does not trace the bar. One possible interpretation is that this population was born prior to bar formation, as their spatial location, kinematics, and pulsation properties suggest, possibly from an accretion event at high redshift.

Stellar activity remains a limiting factor in measuring precise planet parameters from radial velocity spectroscopy, not least in the search for Earth-mass planets orbiting in the habitable zones of Sun-like stars. One approach to mitigate stellar activity is to use combined analyses of both radial velocity and time-series photometry. We present an analysis of simultaneous disk-integrated photometry and radial velocity data of the Sun in order to determine the useful limits of a combined analysis. We find that simple periodogram or autocorrelation analysis of solar photometry give the correct rotation period <50% of the time. We therefore use a Gaussian process to investigate the time variability of solar photometry and to directly compare simultaneous photometry with radial velocity data. We find that the hyperparameter posteriors are relatively stable over 70 yr of solar photometry and the amplitude tracks the solar cycle. We observe good agreement between the hyperparameter posteriors for the simultaneous photometry and radial velocity data. Our primary conclusion is a recommendation to include an additional prior in Gaussian process fits to constrain the evolutionary timescale to be greater than the recurrence timescale (i.e., the rotation period) to recover more physically plausible and useful results. Our results indicate that such simultaneous monitoring may be a useful tool in enhancing the precision of radial velocity surveys.

We present a catalog of high-precision proper motions in the Orion Nebula Cluster (ONC), based on Treasury Program observations with the Hubble Space Telescope's (HST) ACS/WFC camera. Our catalog contains 2454 objects in the magnitude range of 14.2 < m_F775W < 24.7, thus probing the stellar masses of the ONC from ∼0.4 M_☉ down to ∼0.02 M_☉ over an area of ∼550 arcmin². We provide a number of internal velocity dispersion estimates for the ONC that indicate a weak dependence on stellar location and mass. There is good agreement with the published velocity dispersion estimates, although nearly all of them (including ours at ${\sigma }_{v,x}=0.94$ and ${\sigma }_{v,y}=1.25$ mas yr⁻¹) might be biased by the overlapping young stellar populations of Orion A. We identified four new ONC candidate runaways based on HST and the Gaia DR 2 data, all with masses less than ∼1 M_☉. The total census of known candidate runaway sources is 10—one of the largest samples ever found in any Milky Way open star cluster. Surprisingly, none of them have tangential velocities exceeding 20 km s⁻¹. If most of them indeed originated in the ONC, it may compel the re-examination of dynamical processes in very young star clusters. It appears that the mass function of the ONC is not significantly affected by the lost runaways.

We present an analysis of K2 light curves (LCs) from Campaigns 4 and 13 for members of the young (∼3 Myr) Taurus association, in addition to an older (∼30 Myr) population of stars that is largely in the foreground of the Taurus molecular clouds. Out of 156 of the highest-confidence Taurus members, we find that 81% are periodic. Our sample of young foreground stars is biased and incomplete, but nearly all stars (37/38) are periodic. The overall distribution of rotation rates as a function of color (a proxy for mass) is similar to that found in other clusters: the slowest rotators are among the early M spectral types, with faster rotation toward both earlier FGK and later M types. The relationship between period and color/mass exhibited by older clusters such as the Pleiades is already in place by Taurus age. The foreground population has very few stars but is consistent with the USco and Pleiades period distributions. As found in other young clusters, stars with disks rotate on average slower, and few with disks are found rotating faster than ∼2 days. The overall amplitude of the LCs decreases with age, and higher-mass stars have generally lower amplitudes than lower-mass stars. Stars with disks have on average larger amplitudes than stars without disks, though the physical mechanisms driving the variability and the resulting LC morphologies are also different between these two classes.

During the New Horizons spacecraft's encounter with Pluto, the Alice ultraviolet spectrograph conducted a series of observations that detected emissions from both the interplanetary medium (IPM) and Pluto. In the direction of Pluto, the IPM was found to be 133.4 ± 0.6 R at Lyα, 0.24 ± 0.02 R at Lyβ, and <0.10 R at He i 584 Å. We analyzed 3900 s of data obtained shortly before closest approach to Pluto and detect airglow emissions from H i, N i, N ii, N₂, and CO above the disk of Pluto. We find Pluto's brightness at Lyα to be 29.3 ± 1.9 R, in good agreement with preencounter estimates. The detection of the N ii multiplet at 1085 Å marks the first direct detection of ions in Pluto's atmosphere. We do not detect any emissions from noble gases and place a 3σ upper limit of 0.14 R on the brightness of the Ar i 1048 Å line. We compare preencounter model predictions and predictions from our own airglow model, based on atmospheric profiles derived from the solar occultation observed by New Horizons, to the observed brightness of Pluto's airglow. Although completely opaque at Lyα, Pluto's atmosphere is optically thin at wavelengths longer than 1425 Å. Consequently, a significant amount of solar far-UV light reaches the surface, where it can participate in space weathering processes. From the brightness of sunlight reflected from Pluto, we find the surface has a reflectance factor (I/F) of 17% between 1400 and 1850 Å. We also report the first detection of a C₃ hydrocarbon molecule, methylacetylene, in absorption, at a column density of ∼5 × 10¹⁵ cm⁻², corresponding to a column-integrated mixing ratio of 1.6 × 10⁻⁶.

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in X-ray and ultraviolet (XUV) irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on an ultraviolet continuum flux measurement (centered around ∼1900 Å) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photodissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 × 10⁻⁴L_*, while the full flux including non-flaring periods is 6.1 × 10⁻⁵L_*, when L_* is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass loss and tidal evolution and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio $Q^{\prime} =3Q/2{k}_{2}$ for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator VPLanet to study ocean retention for these planets using our derived system parameters.

Four planetary nebulae (PNe) are considered to be probable or possible members of Galactic globular clusters (GCs). These are Ps 1 = K648 in M15, GJJC 1 = IRAS 18333−2357 in M22, JaFu 1 in Palomar 6, and JaFu 2 in NGC 6441. In addition to lying close to the host GCs on the sky, these PNe have radial velocities that are consistent, within the errors and stellar velocity dispersions, with cluster membership. The remaining membership criterion is whether the proper motions (PMs) of the central stars are in agreement with those of the host clusters. We have carried out the PM test for all four PNe. Two of the central stars—those of Ps 1 and GJJC 1—have PMs listed in the recent Gaia Data Release 2 (DR2). We updated the PM of the Ps 1 central star to a more precise value using archival Hubble Space Telescope (HST) frames. Both PMs are statistically consistent with cluster membership. For the other two PNe, we used archival HST images to derive the PMs of their nuclei. For JaFu 2, there are HST images at several epochs, and the measured PM of the nucleus is in excellent agreement with that of the host cluster. For JaFu 1 the available archival HST images are less optimal and the results are less conclusive; the measured PM for the central star is marginally consistent with cluster membership, but additional astrometric observations are desirable for a more robust membership test.

Planets and satellites orbiting a binary system exist in the solar system and extrasolar planetary systems. Their orbits can be significantly different from Keplerian orbits, if they are close to the binary and the secondary-to-primary mass ratio is high. A proper description of a circumbinary orbit is in terms of the free eccentricity e_free at the epicyclic frequency κ₀, forced eccentricity e_forced at the mean motion n₀, and oscillations at higher frequencies forced by the non-axisymmetric components of the binary's potential. We show that accurate numerical values for the amplitudes and frequencies of these terms can be extracted from numerical orbit integrations by applying fast Fourier transformation (FFT) to the cylindrical distance between the circumbinary object and the center of mass of the binary as a function of time. We apply this method to three Kepler circumbinary planets and the satellites of Pluto-Charon. For the satellite Styx of Pluto-Charon, the FFT results for κ₀ and e_free differ significantly from the first-order analytic value and the value reported by Showalter & Hamilton, respectively. We show that the deviation in κ₀ is likely due to the effect of the 3:1 mean-motion resonance and discuss the implications of the lower value for e_free.

Infrared observations of metastable 2³S helium absorption with ground- and space-based spectroscopy are rapidly maturing, as this species is a unique probe of exoplanet atmospheres. Specifically, the transit depth in the triplet feature (with vacuum wavelengths near 1083.3 nm) can be used to constrain the temperature and mass-loss rate of an exoplanet's upper atmosphere. Here, we present a new photometric technique to measure metastable 2³S helium absorption using an ultranarrowband filter (FWHM 0.635 nm) coupled to a beam-shaping diffuser installed in the Wide-field Infrared Camera on the 200 inch Hale Telescope at Palomar Observatory. We use telluric OH lines and a helium arc lamp to characterize refractive effects through the filter and to confirm our understanding of the filter transmission profile. We benchmark our new technique by observing a transit of WASP-69b and detect an excess absorption of 0.498% ± 0.045% (11.1σ), consistent with previous measurements after considering our bandpass. We then use this method to study the inflated gas giant WASP-52b and place a 95th percentile upper limit on excess absorption in our helium bandpass of 0.47%. Using an atmospheric escape model, we constrain the mass-loss rate for WASP-69b to be ${5.25}_{-0.46}^{+0.65}\times {10}^{-4}\,{M}_{{\rm{J}}}\,{\mathrm{Gyr}}^{-1}$ (${3.32}_{-0.56}^{+0.67}\times {10}^{-3}\,{M}_{{\rm{J}}}\,{\mathrm{Gyr}}^{-1}$) at 7000 K (12,000 K). Additionally, we set an upper limit on the mass-loss rate of WASP-52b at these temperatures of $2.1\times {10}^{-4}\,{M}_{{\rm{J}}}\,{\mathrm{Gyr}}^{-1}$ ($2.1\times {10}^{-3}\,{M}_{{\rm{J}}}\,{\mathrm{Gyr}}^{-1}$). These results show that ultranarrowband photometry can reliably quantify absorption in the metastable helium feature.

Exoplanet catalogs produced by surveys suffer from a lack of completeness (not every planet is detected) and less than perfect reliability (not every planet in the catalog is a true planet), particularly near the survey's detection limit. Exoplanet occurrence rate studies based on such a catalog must be corrected for completeness and reliability. The final Kepler data release, DR25, features a uniformly vetted planet candidate catalog and data products that facilitate corrections. We present a new probabilistic approach to the characterization of Kepler completeness and reliability, making full use of the Kepler DR25 products. We illustrate the impact of completeness and reliability corrections with a Poisson-likelihood occurrence rate method, using a recent stellar properties catalog that incorporates Gaia stellar radii and essentially uniform treatment of the stellar population. Correcting for reliability has a significant impact: the exoplanet occurrence rate for orbital period and radius within 20% of Earth's around GK dwarf stars, corrected for reliability, is ${0.015}_{-0.007}^{+0.011}$, whereas not correcting results in ${0.034}_{-0.012}^{+0.018}$—correcting for reliability reduces this occurrence rate by more than a factor of two. We further show that using Gaia-based versus DR25 stellar properties impacts the same occurrence rate by a factor of two. We critically examine the the DR25 catalog and the assumptions behind our occurrence rate method. We propose several ways in which confidence in both the Kepler catalog and occurrence rate calculations can be improved. This work provides an example of how the community can use the DR25 completeness and reliability products.

An accurate and precise Kepler Stellar Properties Catalog is essential for the interpretation of the Kepler exoplanet survey results. Previous Kepler Stellar Properties Catalogs have focused on reporting the best-available parameters for each star, but this has required combining data from a variety of heterogeneous sources. We present the Gaia–Kepler Stellar Properties Catalog, a set of stellar properties of 186,301 Kepler stars, homogeneously derived from isochrones and broadband photometry, Gaia Data Release 2 parallaxes, and spectroscopic metallicities, where available. Our photometric effective temperatures, derived from $g\,\mathrm{to}\,{K}_{s}$ colors, are calibrated on stars with interferometric angular diameters. Median catalog uncertainties are 112 K for ${T}_{\mathrm{eff}}$, 0.05 dex for $\mathrm{log}\,g$, 4% for ${R}_{\star }$, 7% for ${M}_{\star }$, 13% for ${\rho }_{\star }$, 10% for ${L}_{\star }$, and 56% for stellar age. These precise constraints on stellar properties for this sample of stars will allow unprecedented investigations into trends in stellar and exoplanet properties as a function of stellar mass and age. In addition, our homogeneous parameter determinations will permit more accurate calculations of planet occurrence and trends with stellar properties.

We propose several descriptive measures to characterize the arrangements of planetary masses, periods, and mutual inclinations within exoplanetary systems. These measures are based on complexity theory and capture the global, system-level trends of each architecture. Our approach considers all planets in a system simultaneously, facilitating both intrasystem and intersystem analysis. We find that based on these measures, Kepler's high-multiplicity (N ≥ 3) systems can be explained if most systems belong to a single intrinsic population, with a subset of high-multiplicity systems (∼20%) hosting additional, undetected planets intermediate in period between the known planets. We confirm prior findings that planets within a system tend to be roughly the same size and approximately coplanar. We find that forward modeling has not yet reproduced the high degree of spacing similarity (in log-period) actually seen in the Kepler data. Although our classification scheme was developed using compact Kepler multis as a test sample, our methods can be immediately applied to any other population of exoplanetary systems. We apply this classification scheme to (1) quantify the similarity between systems, (2) resolve observational biases from physical trends, and (3) identify which systems to search for additional planets and where to look for these planets.

We have performed a survey for new members of the Ophiuchus cloud complex using high-precision astrometry from the second data release of Gaia, proper motions measured with multi-epoch images from the Spitzer Space Telescope, and color–magnitude diagrams constructed with photometry from various sources. Through spectroscopy of candidates selected with those data, we have identified 155 new young stars. Based on available measurements of kinematics, we classify 102, 47, and 6 of those stars as members of Ophiuchus, Upper Sco, and other populations in Sco–Cen, respectively. We have also assessed the membership of all other stars in the vicinity of Ophiuchus that have spectroscopic evidence of youth from previous studies, arriving at a catalog of 373 adopted members of the cloud complex. For those adopted members, we have compiled mid-infrared photometry from Spitzer and the Wide-field Infrared Survey Explorer and have used mid-infrared colors to identify and classify circumstellar disks. We find that 210 of the members show evidence of disks, including 48 disks that are in advanced stages of evolution. Finally, we have estimated the relative median ages of the populations near the Ophiuchus clouds and the surrounding Upper Sco association using absolute K-band magnitudes (M_K) based on Gaia parallaxes. If we adopt an age 10 Myr for Upper Sco, then the relative values of M_K imply median ages of ∼2 Myr for L1689 and embedded stars in L1688, 3–4 Myr for low-extinction stars near L1688, and ∼6 Myr for the group containing ρ Oph.

The treatment of systematic noise is a significant aspect of transit exoplanet data processing due to the signal strength of systematic noise relative to a transit signal. Typically, the standard approach to transit detection is to estimate and remove systematic noise independently of and prior to a transit detection test. If a transit signal is present in a light curve, the process of systematic noise removal may distort the transit signal by overfitting and thereby reduce detection efficiency. We present a Bayesian framework for joint detection of transit signals and systematic noise characterization and describe the implementation of these detectors as optimal Neyman–Pearson likelihood ratio tests. The joint detectors reduce to closed form as matched filters under the assumption of a Gaussian Bayesian prior for the systematic noise. The performance of the exploratory detectors was evaluated in injection tests and show ∼2% improvement in overall detection efficiency relative to the standard approach. We find that joint detection efficiency is specifically improved for short-period, low transit-depth exoplanet transits, providing evidence in support of the hypothesis that joint detection may indeed help to mitigate overfitting. In addition, an initial feasibility test to detect known exoplanets in Kepler data using the joint detectors produced encouraging preliminary results.

We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in ray-tracing dark matter (DM)-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum (PS) and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counterintuitively, we find that using thin lens planes (< 60 h⁻¹ Mpc on a 240 h⁻¹ Mpc simulation box) suppresses the PS over a broad range of scales beyond what would be acceptable for a survey comparable to the Large Synoptic Survey Telescope (LSST). A mass resolution of 7.2 × 10¹¹h⁻¹M_⊙ per DM particle (or 256³ particles in a (240 h⁻¹ Mpc)³ box) is sufficient to extract information using the PS and non-Gaussian statistics from weak lensing data at angular scales down to 1' with LSST-like levels of shape noise.

The longitudes of perihelia and orbital poles of the solar system's dozen or so most remote detected objects are clustered in a manner inconsistent with that of a random sample of uniformly distributed orbits. While small number statistics and observational biases may explain these features, the statistical significance of the clustering has led to the recent development of the "Planet Nine hypothesis." In the proposed scenario, orbits in the distant solar system are shepherded via secular perturbations from an undetected massive planet on an eccentric orbit. However, the precession of perihelia and nodes in the outer Kuiper Belt and inner Oort cloud are also affected by the giant planets, passing stars, and the galactic tide. We perform a large suite of numerical simulations designed to study the orbital alignment of extreme trans-Neptunian objects (ETNOs) and inner Oort cloud objects (IOCOs). In our various integrations that include Planet Nine, we consistently find that ≳60% of ETNOs and IOCOs that are detectable after 4 Gyr are also anti-aligned in perihelia with the distant massive perturber. However, when we randomly select 17 objects from this sample of remaining orbits, there is significant scatter in the degree of longitude of perihelion and orbital pole clustering that might be observed. Furthermore, we argue that, in the absence of Planet Nine, 17 randomly drawn orbits should still exhibit some clustering even if the underlying distribution is uniform. Thus, we find that still more ETNO and IOCO detections are required to confidently infer the presence of Planet Nine.

High-contrast direct imaging of exoplanets can provide many important observables, including measurements of the orbit, spectra that probe the lower layers of the atmosphere, and phase variations of the planet, but cannot directly measure planet radius or mass. Our future understanding of directly imaged exoplanets will therefore rely on extrapolated models of planetary atmospheres and bulk composition, which need robust calibration. We estimate the population of extrasolar planets that could serve as calibrators for these models. Critically, this population of "standard planets" must be accessible to both direct imaging and the transit method, allowing for radius measurement. We show that the search volume of a direct imaging mission eventually overcomes the transit probability falloff with semimajor axis, so that as long as cold planets are not exceedingly rare, the population of transiting planets and directly imageable planets overlaps. Using current extrapolations of Kepler occurrence rates, we estimate that ∼8 standard planets could be characterized shortward of 800 nm with an ambitious future direct imaging mission like LUVOIR-A and several dozen could be detected at the V band. We show the design space that would expand the sample size and discuss the extent to which ground- and space-based surveys could detect this small but crucial population of planets.

We present new kinematic data for the Ophiuchus stellar stream. Spectra have been taken of member candidates at the MMT telescope using Hectospec, Hectochelle, and Binospec, which provide more than 1800 new velocities. Combined with proper-motion measurements of stars in the field by the Gaia—DR2 catalog, we have derived stream membership probabilities, resulting in the detection of more than 200 likely members. These data show the stream extends to more than three times the length shown in the discovery data. A spur to the main stream is also detected. The high-resolution spectra allow us to resolve the stellar velocity dispersion, found to be 1.6 ± 0.3 km s⁻¹.

Analyzing Spitzer and Herschel archival measurements we identified a debris disk around the young K7/M0 star CP−72 2713. The system belongs to the 24 Myr old β Pic moving group. Our new 1.33 mm continuum observation, obtained with the Atacama Large Millimeter/submillimeter Array 7 m array, revealed an extended dust disk with a peak radius of 140 au, probably tracing the location of the planetesimal belt in the system. The disk is outstandingly large compared to known spatially resolved debris disks and protoplanetary disks around stars of comparable masses. The dynamical excitation of the belt at this radius is found to be reconcilable with planetary stirring, while self-stirring by large planetesimals embedded in the belt can work only if these bodies form very rapidly, e.g., via pebble concentration. By analyzing the spectral energy distribution, we derived a characteristic dust temperature of 43 K and a fractional luminosity of 1.1 × 10⁻³. The latter value is prominently high; we know of only four other similarly dust-rich Kuiper Belt analogs within 40 pc of the Sun.

Atmospheric retrievals are now a standard tool to analyze observations of exoplanet atmospheres. This data-driven approach quantitatively compares atmospheric models to observations in order to estimate atmospheric properties and their uncertainties. In this paper, we introduce a new retrieval package, the PHOENIX Exoplanet Retrieval Algorithm (PETRA). PETRA places the PHOENIX atmosphere model in a retrieval framework, allowing us to combine the strengths of a well-tested and widely-used atmosphere model with the advantages of retrieval algorithms. We validate PETRA by retrieving on simulated data for which the true atmospheric state is known. We also show that PETRA can successfully reproduce results from previously published retrievals of WASP-43b and HD 209458b. For the WASP-43b results, we show the effect that different line lists and line profile treatments have on the retrieved atmospheric properties. Lastly, we describe a novel technique for retrieving the temperature structure and e⁻ density in ultrahot Jupiters using H⁻ opacity, allowing us to probe atmospheres devoid of most molecular features with the James Webb Space Telescope.

We present the spectroscopic orbits of 11 nearby, mid-to-late M dwarf binary systems in a variety of configurations: 2 single-lined binaries (SB1s), 7 double-lined binaries (SB2s), 1 double-lined triple (ST2), and 1 triple-lined triple (ST3). Eight of these orbits are the first published for these systems, while five are newly identified multiples. We obtained multi-epoch, high-resolution spectra with the TRES instrument on the 1.5 m Tillinghast Reflector at the Fred Lawrence Whipple Observatory located on Mt. Hopkins in AZ. Using the TiO molecular bands at 7065−7165 Å, we calculated radial velocities for these systems, from which we derived their orbits. We find LHS 1817 to have in a 7 hr period a companion that is likely a white dwarf, due to the ellipsoidal modulation we see in our MEarth-North light-curve data. We find G 123-45 and LTT 11586 to host companions with minimum masses of 41 M_Jup and 44 M_Jup with orbital periods of 35 and 15 days, respectively. We find 2MA 0930+0227 to have a rapidly rotating stellar companion in a 917 day orbital period. GJ 268, GJ 1029, LP 734-34, GJ 1182, G 258-17, and LTT 7077are SB2s with stellar companions with orbital periods of 10, 96, 34, 154, 5, and 84 days; LP 655-43 is an ST3 with one companion in an 18 day orbital period and an outer component in a longer undetermined period. In addition, we present radial velocities for both components of L 870-44AB and for the outer components of LTT 11586 and LP 655-43.

At wavelengths between 0.15 and 0.19 μm, the far-ultraviolet spectrum of Jupiter is dominated by the scattered solar spectrum, attenuated by molecular absorptions primarily by acetylene and ethane, and to a lesser extent ammonia and phosphine. We describe the development of our radiative transfer code that enables the retrieval of abundances of these molecular species from ultraviolet reflectance spectra. As a proof-of-concept we present an analysis of Cassini Ultraviolet Imaging Spectrograph (UVIS) observations of the disk of Jupiter during the 2000/2001 flyby. The ultraviolet-retrieved acetylene abundances in the upper stratosphere are lower than those predicted by models based solely on infrared thermal emission from the mid-stratosphere observed by the Composite Infrared Spectrometer (CIRS), requiring an adjustment to the vertical profiles above 1 mbar. We produce a vertical acetylene abundance profile that is compatible with both CIRS and UVIS, with reduced abundances at pressures <1 mbar: the 0.1 mbar abundances are 1.21 ± 0.07 ppm for acetylene and 20.8 ± 5.1 ppm for ethane. Finally, we perform a sensitivity study for the JUICE ultraviolet spectrograph, which has extended wavelength coverage out to 0.21 μm, enabling the retrieval of ammonia and phosphine abundances, in addition to acetylene and ethane.

To continue our study on chromospheric activity and detection for possible prominence events of the very active RS Canum Venaticorum star SZ Piscium (SZ Psc), long-term high-resolution spectroscopic observations were obtained during several observing runs from 2014 to 2018. Based on the spectral subtraction technique, the chromospheric emission of the Ca ii IRT (λ8662, λ8542, and λ8498), H_α, Na i D₁, D₂ doublet, H_β, and Ca ii H & K lines is mainly associated with the K1 IV primary star of the SZ Psc system, in good agreement with the previous studies, and the F8 V secondary star also shows some chromospheric emission, implying its active chromosphere. Moreover, an optical flare characterized by the He i D₃ line emission together with stronger emission in the other indicators was detected. Furthermore, two chromospheric active longitudes around the two quadratures of the system were identified for most of the time, and the chromospheric activity shows significant changes during a few orbital cycles. The chromospheric activity level seems to show a long-term variation during our observations. There were some excess absorption features in the subtracted H_α line and the other activity indicators, which would be caused by prominence-like materials associated with the K1 IV primary star of the system. Prominence materials could absorb the chromospheric emission and continuum from the K1 IV primary star and even the F8 V secondary one.