Table of contents

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.

For years, scientists have used data from NASA's Kepler Space Telescope to look for and discover thousands of transiting exoplanets. In its extended K2 mission, Kepler observed stars in various regions of the sky all across the ecliptic plane, and therefore in different galactic environments. Astronomers want to learn how the populations of exoplanets are different in these different environments. However, this requires an automatic and unbiased way to identify exoplanets in these regions and rule out false-positive signals that mimic transiting planet signals. We present a method for classifying these exoplanet signals using deep learning, a class of machine learning algorithms that have become popular in fields ranging from medical science to linguistics. We modified a neural network previously used to identify exoplanets in the Kepler field to be able to identify exoplanets in different K2 campaigns that exist in a range of galactic environments. We train a convolutional neural network, called AstroNet-K2, to predict whether a given possible exoplanet signal is really caused by an exoplanet or a false positive. AstroNet-K2 is highly successful at classifying exoplanets and false positives, with accuracy of 98% on our test set. It is especially efficient at identifying and culling false positives, but for now, it still needs human supervision to create a complete and reliable planet candidate sample. We use AstroNet-K2 to identify and validate two previously unknown exoplanets. Our method is a step toward automatically identifying new exoplanets in K2 data and learning how exoplanet populations depend on their galactic birthplace.

We present new methodological features and physical ingredients included in the one-dimensional radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY andPHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature–pressure (T–P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, and alkali atoms (and ionized hydrogen) playing an important role in the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce nonmonotonic T–P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80°, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ∼10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative–convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations.

We visually analyzed the transit timing variation (TTV) data of 5930 Kepler Objects of Interest (KOIs) homogeneously. Using data from Rowe et al. and Holczer et al., we investigated TTVs for nearly all KOIs in Kepler's Data Release 24 catalog. Using TTV plots, periodograms, and phase-folded quadratic plus sinusoid fits, we visually rated each KOI's TTV data in five categories. Our ratings emphasize the hundreds of planets with TTVs that are weaker than the ∼200 that have been studied in detail. Our findings are consistent with statistical methods for identifying strong TTVs, though we found some additional systems worth investigation. Between about 3–50 days and 1.3–6 Earth radii, the frequency of strong TTVs increases with period and radius. As expected, strong TTVs are very common when period ratios are near a resonance, but there is not a one-to-one correspondence. The observed planet-by-planet frequency of strong TTVs is only somewhat lower in systems with one or two known planets (7% ± 1%) than in systems with three or more known planets (11% ± 2%). We attribute TTVs to known planets in multitransiting systems but find ∼30 cases where the perturbing planet is unknown. Our conclusions are valuable as an ensemble for learning about planetary system architectures and individually as stepping stones toward more-detailed mass–radius constraints. We also discuss Data Release 25 TTVs, investigate ∼100 KOIs with transit duration and/or depth variations, and estimate that the Transiting Exoplanet Survey Satellite will likely find only ∼10 planets with strong TTVs.

We measure the rotation periods of 19 stars in the Kepler transiting planetary systems, P_rot,astero from asteroseismology and P_rot,phot from the photometric variation of their light curves. Two stars exhibit two clear peaks in the Lomb–Scargle periodogram, neither of which agrees with the seismic rotation period. Other four systems do not show any clear peak, whose stellar rotation period is impossible to estimate reliably from the photometric variation; their stellar equators may be significantly inclined with respect to the planetary orbital plane. For the remaining 13 systems, P_rot,astero and ${P}_{\mathrm{rot},\mathrm{phot}}$ agree within 30%. Interestingly, 3 out of the 13 systems are in the spin–orbit resonant state in which ${P}_{\mathrm{orb},{\rm{b}}}/{P}_{\mathrm{rot},\mathrm{astero}}\approx 1$ with P_orb,b being the orbital period of the innermost planet of each system. The corresponding chance probability is (0.2–4.7)% based on the photometric rotation period data for 464 Kepler transiting planetary systems. While further analysis of stars with reliable rotation periods is required to examine the statistical significance, the spin–orbit resonance between the star and planets, if confirmed, has important implications for the star–planet tidal interaction, in addition to the origin of the spin–orbit (mis-)alignment of transiting planetary systems.

Following our previous detection of ubiquitous ${{\rm{H}}}_{2}{\rm{O}}$ and ${{\rm{O}}}_{2}$ absorption against the far-ultraviolet continuum of stars located near the nucleus of Comet 67P/Churyumov–Gerasimenko, we present a serendipitously observed stellar occultation that occurred on 2015 September 13, approximately one month after the comet's perihelion passage. The occultation appears in two consecutive 10-minute spectral images obtained by Alice, Rosetta's ultraviolet (700–2100 Å) spectrograph, both of which show ${{\rm{H}}}_{2}{\rm{O}}$ absorption with column density >10^17.5 cm⁻² and significant ${{\rm{O}}}_{2}$ absorption (${{\rm{O}}}_{2}$/${{\rm{H}}}_{2}{\rm{O}}$ ≈ 5%–10%). Because the projected distance from the star to the nucleus changes between exposures, our ability to study the ${{\rm{H}}}_{2}{\rm{O}}$ column density profile near the nucleus (impact parameters <1 km) is unmatched by our previous observations. We find that the ${{\rm{H}}}_{2}{\rm{O}}$ and ${{\rm{O}}}_{2}$ column densities decrease with increasing impact parameter, in accordance with expectations, but the ${{\rm{O}}}_{2}$ column decreases ∼3 times more quickly than ${{\rm{H}}}_{2}{\rm{O}}$. When combined with previously published results from stellar appulses, we conclude that the ${{\rm{O}}}_{2}$ and ${{\rm{H}}}_{2}{\rm{O}}$ column densities are highly correlated, and ${{\rm{O}}}_{2}$/${{\rm{H}}}_{2}{\rm{O}}$ decreases with the increasing ${{\rm{H}}}_{2}{\rm{O}}$ column.

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

N-body simulations have shown that a bar in a galaxy can be triggered by two processes: (1) by its own instabilities in the disk, or (2) by interactions with other galaxies. Both mechanisms have been widely studied. However, the literature has not shown measurements of the critical limits of the disk stability parameters (DSPs). We show measurements of those parameters through the whole evolution in isolated disk models and find that the initial rotation configuration of those models stays in the stable or unstable regime from the initial to the final evolution. Then we perturbed the isolated models to study the evolution of DSPs under perturbation. We find that the critical limits of DSPs are not much affected in barred models, but when the bar is triggered by a perturbation, the disk falls into the unstable regimen. We show in our models that a bar triggered by a light perturbation grows in two phases: first, the bar appears as a slow rotator, and then it evolves to be a fast rotator; second, when the perturbation is far from the target galaxy, the bar evolves from fast to slow rotator. When the bar is triggered by a heavy perturbation, it appears as a fast rotator and evolves to be a slow rotator, similar to classical bar models.

We present a stacking analysis of 2.61 Ms of archival Chandra observations of stellar wind bow shocks. We place an upper limit on the X-ray luminosity of infrared-detected bow shocks of <2 × 10²⁹ erg s⁻¹, a more stringent constraint than has been found in previous archival studies and dedicated observing campaigns of nearby bow shocks. We compare the X-ray luminosities and L_X/L_bol ratios of bow shock driving stars to those of other OB stars within the Chandra field of view. Driving stars are, on average, of later spectral type than the field-of-view OB stars, and we do not observe any unambiguously high L_X/L_bol ratios indicative of magnetic stars in our sample. We additionally assess the feasibility of detecting X-rays from stellar wind bow shocks with the proposed Lynx X-ray Observatory. If the X-ray flux originating from the bow shocks is just below our Chandra detection limit, the nearest bow shock in our sample (at ∼0.4 kpc with an absorbing column of ∼10²¹ cm⁻²) should be observable with Lynx in exposure times on the order of ∼100 ks.

Snowden & Young suggested that the reason why there are GK subgiants is because they are members of binaries, which would bring them above the main sequence in an Hertzsprung–Russell (HR) diagram. They studied a sample of 30 G0-K1 IV stars and were disappointed to find only two to be spectroscopic binaries. With more accurate radial velocities I found seven binaries in their samples of subgiants and control stars; orbital elements are given for those seven. Using Hipparcos parallaxes and SIMBAD data, I found that nearly all of the G0-K1 IV stars fall on the evolutionary tracks by Garardi et al. for Population I stars with masses of 0.9–1.9 M_⊙ and ages of up to 10¹⁰ yr, which are normal parameters for nearby field stars. Therefore there is no problem regarding the existence of GK subgiants.

We present a visible-light full orbital phase curve of the transiting planet WASP-18b measured by the TESS mission. The phase curve includes the transit, secondary eclipse, and sinusoidal modulations across the orbital phase shaped by the planet's atmospheric characteristics and the star–planet gravitational interaction. We measure the beaming (Doppler boosting) and tidal ellipsoidal distortion phase modulations and show that the amplitudes of both agree with theoretical expectations. We find that the light from the planet's dayside hemisphere occulted during secondary eclipse, with a relative brightness of ${341}_{-18}^{+17}$ ppm, is dominated by thermal emission, leading to an upper limit on the geometric albedo in the TESS band of 0.048 ($2\sigma$). We also detect the phase modulation due to the planet's atmosphere longitudinal brightness distribution. We find that its maximum is well aligned with the substellar point to within 29 ($2\sigma$). We do not detect light from the planet's nightside hemisphere, with an upper limit of 43 ppm ($2\sigma$), which is 13% of the dayside brightness. The low albedo, lack of atmospheric phase shift, and inefficient heat distribution from the day to night hemispheres that we deduce from our analysis are consistent with theoretical expectations and similar findings for other strongly irradiated gas giant planets. This work demonstrates the potential of TESS data for studying the full orbital phase curves of transiting systems. Finally, we complement our study by looking for transit timing variations (TTVs) in the TESS data combined with previously published transit times, although we do not find a statistically significant TTV signal.

The altitude distribution of meteors detected by a radar is sensitive to the instrument's response function and can thus provide insight into the physical processes involved in radar measurements. This, in turn, can be used to determine the rate of ablation and ionization of the meteoroids and ultimately the input flux on Earth. In this work, we model the radar meteor head echo altitude distribution for three High Power and Large Aperture radar systems, by considering meteoroid populations from the main cometary family sources. In this simulation, we first use the results of a dynamical model of small meteoroids impacting Earth's upper atmosphere to model the incoming mass, velocity, and entry angular distributions. We then combine these with the Chemical Ablation Model and establish the meteoroid ionization rates as a function of mass, velocity, and entry angle in order to determine the altitude at which these radars should detect the produced meteors and the portion of produced meteors from each population that are detected by these radars. We explore different sizes of head plasma as well as the possible effects on radar scattering of the head echo aspect sensitivity. We find that the modeled altitude distributions are generally in good agreement with measurements, particularly for ultra-high-frequency radars. In addition, our results indicate that the number of particles from Jupiter Family Comets (JFCs) required to fit the observations is lower than predicted by astronomical models. It is not clear yet if this discrepancy is due to the overprediction of JFC meteoroids by dynamical models or due to unaccounted physical processes in the treatment of ablation, ionization, and detections of meteoroids as they pass through Earth's atmosphere.

Over 100 rocky planets orbiting Sun-like stars in very short orbital periods (≲1 day) have been discovered by the Kepler mission. The origin of these planets, known as ultra-short-period (USP) planets, remains elusive. Here, we propose that most of these planets, originally at periods of ∼5–10 days, reach their current orbits via high-eccentricity migration. In a scaled-down version of the dynamics that may have been experienced by their high-mass analogs, the hot Jupiters, these planets reach high eccentricities via chaotic secular interactions with their companion planets and then undergo orbital circularization due to dissipation from tides raised on the planet. This proposal is motivated by the following observations: planetary systems observed by Kepler often contain several super-Earths with non-negligible eccentricities and inclinations, possibly extending beyond ∼au distances; by contrast, only a small fraction of USP planets have known transiting companions, which are generally not closely spaced, and we argue that most of them should have companions with periods ≳10 days. The proposed scenario naturally explains the observation that most USP planets have significantly more distant transiting companions compared to their counterparts at slightly longer periods (1–3 days). Our model predicts that USP planets should have: (i) spin–orbit angles, and inclinations relative to outer planets, in the range of ∼10–50°; (ii) several outer planetary companions extending beyond ∼1 au distances. Both of these predictions may be tested by TESS and its follow-up observations.

We develop an evolutionary model of the long-period comet (LPC) population, starting from their birthplace in a massive trans-Neptunian disk that was dispersed by migrating giant planets. Most comets that remain bound to the solar system are stored in the Oort cloud. Galactic tides and passing stars make some of these bodies evolve into observable comets in the inner solar system. Our approach models each step in a full-fledged numerical framework. Subsequent analysis consists of applying plausible fading models and computing the original orbits to compare with observations. Our results match the observed semimajor axis distribution of LPCs when Whipple's power-law fading scheme with an exponent $\kappa ={0.6}_{-0.2}^{+0.1}$ is adopted. The cumulative perihelion (q) distribution is well fit by a linear increase plus a weak quadratic term. Beyond q = 15 au, however, the population increases steeply, and the isotropy of LPC orbital planes breaks. We find tentative evidence from the perihelion distribution of LPCs that the returning comets are depleted in supervolatiles and become active due to water ice sublimation for q ≤ 3 au. Using an independent calibration of the population of the initial disk, our predicted LPC flux is smaller than observations suggest by a factor of ≃2. Current data only characterize comets from the outer Oort cloud (semimajor axes ≳10⁴ au). A true boost in understanding the Oort cloud's structure should result from future surveys when they detect LPCs with perihelia beyond 15 au. Our results provide observational predictions of what can be expected from these new data.

Refraction by the atmosphere causes the positions of sources to depend on the airmass through which an observation was taken. This shift is dependent on the underlying spectral energy of the source and the filter or bandpass through which it is observed. Wavelength-dependent refraction within a single passband is often referred to as differential chromatic refraction (DCR). With a new generation of astronomical surveys undertaking repeated observations of the same part of the sky over a range of different airmasses and parallactic angles, DCR should be a detectable and measurable astrometric signal. In this paper we introduce a novel procedure that takes this astrometric signal and uses it to infer the underlying spectral energy distribution of a source; we solve for multiple latent images at specific wavelengths via a generalized deconvolution procedure built on robust statistics. We demonstrate the utility of such an approach for estimating a partially deconvolved image, at higher spectral resolution than the input images, for surveys such as the Large Synoptic Survey Telescope.

We present results from the analysis of WIYN pODI imaging of 23 ultracompact high-velocity clouds (UCHVCs), which were identified in the ALFALFA H i survey as possible dwarf galaxies in or near the Local Group. To search for a resolved stellar population associated with the H i gas in these objects, we carried out a series of steps designed to identify stellar overdensities in our optical images. We identify five objects that are likely stellar counterparts to the UCHVCs, at distances of ∼350 kpc to ∼1.6 Mpc. Two of the counterparts were already described in Janesh et al.; the estimated distance and detection significance for one of them changed in the final analysis of the full pODI data set. At their estimated distances, the detected objects have H i masses from 2 × 10⁴ to 3 × 10⁶M_⊙, M_V from −1.4 to −7.1, and stellar masses from 4 × 10² to 4 × 10⁵M_⊙. None of the objects shows evidence of a young stellar population. Their properties would make the UCHVCs some of the most extreme objects in and around the Local Group, comparable to ultrafaint dwarf galaxies in their stellar populations, but with significant gas content. Such objects probe the extreme end of the galaxy mass function and provide a test bed for theories regarding the baryonic feedback processes that impact star formation and galaxy evolution in this low-mass regime.

We report 18 new primary minima timing observations of the short-period eclipsing binary system NY Virginis. We combined these minima with previously published primary minima to update circumbinary exoplanet models in this system based on O − C timing variations. We performed a nonlinear least-squares minimization search using a quadratic ephemeris and either one or two exoplanets. The only model with an acceptable fit includes a period derivative $\dot{P}=2.83\times {10}^{-12}$ and two planets in eccentric orbits e = 0.15, 0.15 with minimum masses of 2.7 and 5.5 Jovian masses. Analysis of the orbit stability shows that this solution is stable for at least 10⁸ years, but a small increase in eccentricity (e ≥ 0.20) for either planet renders the orbits unstable in less than 10⁶ years. A number of model parameters are significantly degenerate, so additional observations are required to determine planetary parameters with high statistical confidence.

We present new observations of the multiplanet system HIP 41378, a bright star (V = 8.9, K_s = 7.7) with five known transiting planets. Previous K2 observations showed multiple transits of two Neptune-sized bodies and single transits of three larger planets (R_P = 0.33R_J, 0.47R_J, 0.88R_J). K2 recently observed the system again in Campaign 18 (C18). We observe one new transit each of two of the larger planets d/f, giving maximal orbital periods of 1114/1084 days, as well as integer divisions of these values down to a lower limit of about 50 days. We use all available photometry to determine the eccentricity distributions of HIP 41378 d & f, finding that periods ≲300 days require non-zero eccentricity. We check for overlapping orbits of planets d and f to constrain their mutual periods, finding that short periods (P < 300 days) for planet f are disfavored. We also observe transits of planets b and c with Spitzer/Infrared Array Camera (IRAC), which we combine with the K2 observations to search for transit timing variations (TTVs). We find a linear ephemeris for planet b, but see a significant TTV signal for planet c. The ability to recover the two smaller planets with Spitzer shows that this fascinating system will continue to be detectable with Spitzer, CHEOPS, TESS, and other observatories, allowing us to precisely determine the periods of d and f, characterize the TTVs of planet c, recover the transits of planet e, and further enhance our view of this remarkable dynamical laboratory.

We observed comet 96P/Machholz 1 on a total of nine nights before and after perihelion during its 2017/2018 apparition. Both its unusually small perihelion distance and the observed fragmentation during multiple apparitions make 96P an object of great interest. Our observations show no evidence of a detectable dust coma, implying that we are observing a bare nucleus at distances ranging from 2.3 to 3.8 au. Based on this assumption, we calculated its color and found average values of g'–r' = 0.50 ± 0.04, r'–i' = 0.17 ± 0.03, and i'–z' = 0.06 ± 0.04. These are notably more blue than those of the nuclei of other Jupiter-family and long-period comets. Furthermore, assuming a bare nucleus, we found an equivalent nuclear radius of 3.4 ± 0.2 km with an axial ratio of at least 1.6 ± 0.1. The lightcurve clearly displays one large peak, one broad flat peak, and two distinct troughs, with a clear asymmetry that suggests that the shape of the nucleus deviates from that of a simple triaxial ellipsoid. This asymmetry in the lightcurve allowed us to constrain the nuclear rotation period to 4.10 ± 0.03 hr and 4.096 ± 0.002 hr before and after perihelion, respectively. Within the uncertainties, 96P's rotation period does not appear to have changed throughout the apparition, and we conclude a maximum possible change in rotation period of 130 s. The observed properties were compared to those of comet 322P and interstellar object 1I/'Oumuamua in an attempt to study the effects of close perihelion passages on cometary surfaces and their internal structure and the potential interstellar origin of 96P.

To discover Earth analogs around other stars, next generation spectrographs must measure radial velocity with 10 cm s⁻¹ precision. Since even microtellurics can induce RV errors of up to 50 cm s⁻¹, achieving 10 cm s⁻¹ precision requires precise modeling of telluric absorption features. The standard approaches to telluric modeling are (a) observing a standard star and (b) using a radiative transfer code. Observing standard stars, however, takes valuable observing time away from science targets. Radiative transfer codes, meanwhile, may omit microtelluric features, which are an important contributor to the RV error budget at 10 cm s⁻¹. To address these issues, we present a telluric model of the self-calibrating, empirical, light-weight linear regression telluric (SELENITE) model for high-resolution spectra. The model exploits two simple observations: (a) water tellurics grow proportionally to precipitable water vapor and therefore proportionally to each other and (b) non-water tellurics grow proportionally to airmass. Water tellurics can be identified by looking for pixels whose growth correlates with a known calibration water telluric and modeled by regression against it, and likewise non-water tellurics with airmass. The model does not require line data, water vapor measurements, or additional observations (beyond one-time calibration observations), achieves fits with a ${\chi }_{\mathrm{red}}^{2}$ of 1.17 on B stars and 2.95 on K dwarfs, and leaves residuals of 1% (B stars) and 1.1% (K dwarfs) of continuum. Fitting takes seconds on laptop PCs; SELENITE is light-weight enough to guide observing runs.

Space-based direct imaging missions (HabEx, LUVOIR) would observe reflected light from exoplanets in the habitable zones of Sun-like stars. The ultimate—but not sole—goal of these concept missions is to characterize such planets. Knowing an exoplanet's orbit would help twofold: (i) its semimajor axis informs whether the planet might harbor surface liquid water, making it a priority target; and (ii) predicting the planet's future location would tell us where and when to look. The science yields of HabEx and LUVOIR depend on the number, cadence, and precision of observations required to establish a planet's orbit. We produce mock observations using realistic distributions for the six Keplerian orbital parameters, experimenting with both beta and uniform eccentricity distributions, and accounting for imperfect astrometry (σ = 3.5 mas) and obscuration due to the inner working angle of a high-contrast imaging system (inner working angle = 31 mas). Using Markov chain Monte Carlo methods, we fit the orbital parameters, and retrieve their average precisions and accuracies as functions of cadence, number of epochs, and distance to the target. Given the time at which it was acquired, each image provides two data: the x and y position of the planet with respect to its star. Parameter retrieval based on one or two images is formally underconstrained, yet the semimajor axis posterior can be obtained semi-analytically. For a planet at 1 au around a star at a distance of 10 pc, three epochs constrain the semimajor axis to within ≲5%, if each image is taken at least 90 days apart.

Rotation and orbital eccentricity both strongly influence planetary climate. Eccentricities can often be measured for exoplanets, but rotation rates are currently difficult or impossible to constrain. Here we examine how the combined effects of rotation and eccentricity on observed emission from ocean-rich terrestrial planets can be used to infer their rotation rates in circumstances where their eccentricities are known. We employ an Earth climate model with no land and a slab ocean, and consider two eccentricities (e = 0.3 and 0.6) and two rotation rates: a fast Earth-like period of 24 hr, and a slower pseudo-synchronous period that generalizes spin synchronization for eccentric orbits. We adopt bandpasses of the Mid-Infrared Instrument on the James Webb Space Telescope as a template for future photometry. At e = 0.3 the rotation rates can be distinguished if the planet transits near periastron, because slow rotation produces a strong day–night contrast and thus an emission minimum during periastron. However, light curves behave similarly if the planet is eclipsed near periastron, as well as for either viewing geometry at e = 0.6. Rotation rates can nevertheless be distinguished using ratios of emission in different bands, one in the water vapor window with another in a region of strong water absorption. These ratios vary over an orbit by ≲0.1 dex for Earth-like rotation, but by 0.3–0.5 dex for pseudo-synchronous rotation because of large day–night contrast in upper-tropospheric water. For planets with condensible atmospheric constituents in eccentric orbits, rotation regimes might thus be distinguished with infrared observations for a range of viewing geometries.

As part of our multi-observatory, multifilter campaign, we present r–i color observations of 82 near-Earth objects (NEOs) obtained with the reionization and transients infrared camera (RATIR) instrument on the 1.5 m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid-response observations of small (≲850 m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of ∼30–850 m is 1.1, which is in accordance with our previous results. We also find that 10% of all NEOs in our sample are neither C- nor S-type asteroids

We report the discovery of TOI-172 b from the Transiting Exoplanet Survey Satellite (TESS) mission, a massive hot Jupiter transiting a slightly evolved G star with a 9.48-day orbital period. This is the first planet to be confirmed from analysis of only the TESS full frame images, because the host star was not chosen as a two-minute cadence target. From a global analysis of the TESS photometry and follow-up observations carried out by the TESS Follow-up Observing Program Working Group, TOI-172 (TIC 29857954) is a slightly evolved star with an effective temperature of T_eff = 5645 ± 50 K, a mass of M_⋆ = ${1.128}_{-0.061}^{+0.065}$M_⊙, radius of R_⋆ = ${1.777}_{-0.044}^{+0.047}$R_⊙, a surface gravity of log g_⋆ = ${3.993}_{-0.028}^{+0.027}$, and an age of ${7.4}_{-1.5}^{+1.6}\ \mathrm{Gyr}$. Its planetary companion (TOI-172 b) has a radius of R_P = ${0.965}_{-0.029}^{+0.032}$R_J, a mass of M_P = ${5.42}_{-0.20}^{+0.22}$M_J, and is on an eccentric orbit ($e={0.3806}_{-0.0090}^{+0.0093}$). TOI-172 b is one of the few known massive giant planets on a highly eccentric short-period orbit. Future study of the atmosphere of this planet and its system architecture offer opportunities to understand the formation and evolution of similar systems.

The discovery of thousands of planetary systems by Kepler has demonstrated that planets are ubiquitous. However, a major challenge has been the confirmation of Kepler planet candidates, many of which still await confirmation. One of the most enigmatic examples is KOI 4.01, Kepler's first discovered planet candidate detection (as KOI 1.01, 2.01, and 3.01 were known prior to launch). Here we present the confirmation and characterization of KOI 4.01 (now Kepler-1658), using a combination of asteroseismology and radial velocities. Kepler-1658 is a massive, evolved subgiant (M_⋆ = 1.45 ± 0.06 M_⊙, R_⋆ = 2.89 ± 0.12 R_⊙) hosting a massive (${M}_{{\rm{p}}}$ = 5.88 ± 0.47 ${M}_{{\rm{J}}}$, ${R}_{{\rm{p}}}$ = 1.07 ± 0.05 ${R}_{{\rm{J}}}$) hot Jupiter that orbits every 3.85 days. Kepler-1658 joins a small population of evolved hosts with short-period (≤100 days) planets and is now the closest known planet in terms of orbital period to an evolved star. Because of its uniqueness and short orbital period, Kepler-1658 is a new benchmark system for testing tidal dissipation and hot Jupiter formation theories. Using all four years of the Kepler data, we constrain the orbital decay rate to be $\dot{P}$ ≤ −0.42 s yr⁻¹, corresponding to a strong observational limit of ${Q}_{\star }^{{\prime} }$ ≥ 4.826 × ${10}^{3}$ for the tidal quality factor in evolved stars. With an effective temperature of ${T}_{\mathrm{eff}}$ ∼ 6200 K, Kepler-1658 sits close to the spin–orbit misalignment boundary at ∼6250 K, making it a prime target for follow-up observations to better constrain its obliquity and to provide insight into theories for hot Jupiter formation and migration.

We present thermal infrared observations of the active asteroid (and Geminid meteoroid stream parent) 3200 Phaethon using the Very Large Telescope. The images, at 10.7 μm wavelength, were taken with Phaethon at its closest approach to Earth (separation of 0.07 au) in 2017 December, at a linear resolution of about 14 km. We probe the Hill sphere (of radius ∼66 km) for trapped dust and macroscopic bodies, finding neither, and we set limits to the presence of unbound dust. The derived limits to the optical depth of dust near Phaethon depend somewhat on the assumed geometry, but are of an order of 10⁻⁵. The upper limit to the rate of loss of mass in dust is ≲14 kg s⁻¹. This is ∼50 times smaller than the rate needed to sustain the Geminid meteoroid stream in steady state. The observations thus show that the production of the Geminids does not proceed in a steady state.

A procedure is presented to improve on measurement of total H i fluxes for extended sources in the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey of neutral hydrogen sources in the nearby universe. A number of tests of the procedure are detailed, and we verify that we recover all of the flux measured with much larger telescope beams. Total fluxes are reported for all sources (1) exceeding 10 Jy km s⁻¹ in the α.100 catalog, or (2) with Uppsala General Catalog diameters 3.0 arcmin or more, or (3) ALFALFA pipeline isophotal ellipse area more than 3.0 times the Arecibo beam. Total fluxes are also provided for a number of confused pairs and small groups including one or more of those high-flux sources. These data should be of use in baryonic Tully–Fisher studies and other applications where the measurement of the total reservoir of neutral atomic gas is important.

We present a low-frequency imaging study of an extended sample of X-shaped radio sources using the Giant Metrewave Radio Telescope (GMRT) at two frequencies (610 and 240 MHz). The sources were drawn from a Very Large Array Radio Images of the Sky at Twenty-Centimeters (FIRST) selected sample and extend an initial GMRT study at the same frequencies of 12 X-shaped radio galaxies predominantly from the 3CR catalog (Lal & Rao 2007). Both the intensity maps and spectral index maps of the 16 newly observed sources are presented. With the combined sample of 28 X-shaped radio sources, we found no systematic differences in the spectral properties of the higher surface brightness, active lobes versus the lower surface brightness, or off-axis emission. The properties of the combined sample are discussed, including the possible role of a twin active galactic nuclei model in the formation of such objects.

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from Apache Point Observatory Galactic Evolution Experiment-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double-lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of ∼0.05–1.5 M_⊙, we find 399 binaries, including both radial velocity (RV) variables and SB2s. The mass ratio distribution of SB2s is consistent with being uniform for q < 0.95 with an excess of twins for q > 0.95. The period distribution is consistent with what has been observed in close binaries (<10 au) in the evolved populations. Three systems are found to have q ∼ 0.1, with a companion located within the brown dwarf desert. There are no strong trends in the multiplicity fraction as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at Σ_* ∼ 30 stars/pc⁻² and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of ∼2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal-mass systems that is less effective in binaries with lower q.

We present scope (simulated CCD observations for photometric experimentation), a Python package to create a forward model of telescope detectors and simulate stellar targets with motion relative to the CCD. The primary application of this package is the simulation of the Kepler Space Telescope detector to predict and characterize increased instrumental noise in the spacecraft's final campaigns of observation. As the fuel powering the spacecraft's stabilizing thrusters ran out and thruster fires began to sputter and fail, stellar point-spread functions experienced more extreme and less predictable motion relative to regions of varied sensitivity on the spacecraft detector, generating more noise in transiting exoplanet light curves. Using our simulations, we demonstrate that current de-trending techniques effectively capture and remove systematics caused by sensitivity variation for spacecraft motion as high as about 10 times that typically experienced by K2. The scope package is open source and has been generalized to allow custom detector and stellar target parameters. Future applications include simulating observations made by the Transiting Exoplanet Survey Satellite and ground-based observations with synthetic atmospheric interference as testbeds for noise removal techniques.

Characterizing the dependence of the orbital architectures and formation environments on the eccentricity distribution of planets is vital for understanding planet formation. In this work, we perform statistical eccentricity studies of transiting exoplanets using transit durations measured via Kepler combined with precise and accurate stellar radii from the California-Kepler Survey and Gaia. Compared to previous works that characterized the eccentricity distribution from transit durations, our analysis benefits from both high-precision stellar radii (∼3%) and a large sample of ∼1000 planets. We observe that systems with only a single observed transiting planet have a higher mean eccentricity ($\bar{e}\sim 0.21$) than systems with multiple transiting planets ($\bar{e}\sim 0.05$), in agreement with previous studies. We confirm the preference for high- and low-eccentricity subpopulations among the single transiting systems. Finally, we show suggestive new evidence that high-e planets in the Kepler sample are preferentially found around high-metallicity ([Fe/H] > 0) stars. We conclude by discussing the implications on planetary formation theories.

The asteroid 2014 JO25, considered to be potentially hazardous by the Minor Planet Center, was spectroscopically followed during its close-Earth encounter on 2017 April 19 and 20. The spectra of the asteroid were taken with the low-resolution spectrograph (LISA), mounted on the 1.2 m telescope at the Mount Abu Infrared Observatory, India. Coming from a region close to the Hungaria population of asteroids, this asteroid follows a comet-like orbit with a relatively high inclination and large eccentricity. Hence, we carried out optical spectroscopic observations of the asteroid to look for comet-like molecular emissions or outbursts. However, the asteroid showed a featureless spectrum, devoid of any comet-like features. The light curve of the asteroid was analyzed using V-band magnitudes derived from the spectra and the most likely solution for the rotation of the asteroid was obtained. The absolute magnitude H and the slope parameter G were determined for the asteroid in the V filter band using the IAU accepted standard two-parameter H–G model. A peculiar, rarely found result from these observations is its phase bluing trend. The relative B–V color index seems to decrease with increasing phase angle, which indicates a phase bluing trend. Such trends have seldom been reported in the literature. However, phase reddening in asteroids is very common. The asymmetry parameter g and the single-scattering albedo w were estimated for the asteroid by fitting the Hapke phase function to the observed data. The asteroid shows a relatively large value for the single-scattering albedo and a highly back-scattering surface.

We report trigonometric parallax and proper motion measurements of 6.7 GHz CH₃OH and 22 GHz H₂O masers in eight high-mass star-forming regions (HMSFRs) based on Very Long Baseline Array (VLBA) observations as part of the Bar and Spiral Structure Legacy (BeSSeL) Survey. The distances of these HMSFRs combined with their Galactic coordinates, radial velocities, and proper motions, allow us to assign them to a segment of the Perseus arm with ℓ ≲ 70°. These HMSFRs are clustered in Galactic longitude from ≈30° to ≈50° neighboring a dearth of such sources between longitudes ≈50° to ≈90°.

Cometary impacts pose a long-term hazard to life on Earth. Impact mitigation techniques have been studied extensively, but they tend to focus on asteroid diversion. Typical asteroid interdiction schemes involve spacecraft physically intercepting the target, a task feasible only for targets identified decades in advance and in a narrow range of orbits—criteria unlikely to be satisfied by a threatening comet. Comets, however, are naturally perturbed from purely gravitational trajectories via solar heating of their surfaces, which activates sublimation-driven jets. Artificial heating of a comet, such as by a laser, may supplement natural heating by the Sun to purposefully manipulate its path and thereby avoid an impact. Deflection effectiveness depends on the comet's heating response, which varies dramatically depending on factors including nucleus size, orbit, and dynamical history. These factors are incorporated into a numerical orbital model to assess the effectiveness and feasibility of using high-powered laser arrays in Earth orbit and on the ground for comet deflection. Simulation results suggest that a diffraction-limited 500 m orbital or terrestrial laser array operating at 10 GW for 1% of each day over 1 yr is sufficient to fully avert the impact of a typical 500 m diameter comet with primary nongravitational parameter A₁ = 2 × 10⁻⁸ au day⁻². Strategies to avoid comet fragmentation during deflection are also discussed.

The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100–200 Myr, when violent impacts occur between large planetesimals of sizes up to protoplanets. Such events are marked by the production of large amounts of debris, as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-sized bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system suggest two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation.

We present the first complete multiband observations of a binary asteroid mutual event. We obtained high-cadence, high-signal-to-noise photometry of the UT 2018 April 9 inferior shadowing event in the Jupiter Trojan binary system Patroclus–Menoetius in four Sloan bands—g', r', i', and z'. We use an eclipse light-curve model to fit for a precise mideclipse time and estimate the minimum separation of the two eclipsing components during the event. Our best-fit mideclipse time of ${2458217.80943}_{-0.00050}^{+0.00057}$ is 19 minutes later than the prediction of Grundy et al. The minimum separation between the center of Menoetius's shadow and the center of Patroclus is 72.5 ± 0.7 km—slightly larger than the predicted 69.5 km. Using the derived light curves, we find no evidence for significant albedo variations or large-scale topographic features on the Earth-facing hemisphere and limb of Patroclus. We also apply the technique of eclipse mapping to place an upper bound of ∼0.15 mag on wide-scale surface color variability across Patroclus.

We describe a new method to derive clean, iodine-free spectra directly from observations acquired using high-resolution echelle spectrographs equipped with iodine cells. The main motivation to obtain iodine-free spectra is to use portions of the spectrum that are superimposed with the dense forest of iodine absorption lines, in order to retrieve lines that can be used to monitor the magnetic activity of the star, helping to validate candidate planets. In short, we provide a straightforward methodology to clean the spectra using the forward model used to derive radial velocities, the line spread function information plus the stellar spectrum without iodine to reconstruct and subtract the iodine spectrum from the observations. We show our results using observations of the star τ Ceti acquired with the Planet Finder Spectrograph (PFS), High Resolution Echelle Spectrometer (HIRES), and University College London Echelle Spectrograph (UCLES), reaching an iodine-free spectrum correction at the ∼1% rms level. We additionally discuss the limitations and further applications of the method.

We present here a reanalysis of the Spitzer Space Telescope phase curves of the hot Jupiter WASP43 b, using the wavelet pixel-independent component analysis, a blind signal-source separation method. The data analyzed were recorded with the Infrared Array Camera and consisted of two visits at 3.6 μm, and one visit at 4.5 μm, each visit covering one transit and two eclipse events. To test the robustness of our technique we repeated the analysis on smaller portions of the phase curves, and by employing different instrument ramp models. Our reanalysis presents significant updates of the planetary parameters compared to those reported in the original phase curve study of WASP43 b. In particular, we found (1) higher nightside temperatures, (2) smaller hotspot offsets, (3) a greater consistency (∼1σ) between the two 3.6 μm visits, and (4) a greater similarity with the predictions of the atmospheric circulation models. Our parameter results are consistent within 1σ with those reported by a recent reanalysis of the same data sets. For each visit we studied the variation of the retrieved transit parameters as a function of various sets of stellar limb-darkening coefficients, finding significant degeneracy between the limb-darkening models and the analysis output. Furthermore, we performed the analysis of the single transit and eclipse events, and we examined the differences between these results with the ones obtained with the whole phase curve. Finally we provide a formula useful to optimize the trade-off between precision and duration of observations of transiting exoplanets.

Accurate estimations of atmospheric properties of exoplanets from transmission spectra require the understanding of degeneracies between model parameters and observations that can resolve them. We conduct a systematic investigation of such degeneracies using a combination of detailed atmospheric retrievals and a range of model assumptions, focusing on H₂-rich atmospheres. As a case study, we consider the well-studied hot Jupiter HD 209458 b. We perform extensive retrievals with models ranging from simple isothermal and isobaric atmospheres to those with full pressure–temperature profiles, inhomogeneous cloud/haze coverage, multiple-molecular species, and data in the optical–infrared wavelengths. Our study reveals four key insights. First, we find that a combination of models with minimal assumptions and broadband transmission spectra with current facilities allows precise estimates of chemical abundances. In particular, high-precision optical and infrared spectra, along with models including variable cloud coverage and prominent opacity sources, with Na and K being important in the optical, provide joint constraints on cloud/haze properties and chemical abundances. Second, we show that the degeneracy between planetary radius and its reference pressure is well characterized and has little effect on abundance estimates, contrary to previous claims using semi-analytic models. Third, collision-induced absorption due to H₂–H₂ and H₂–He interactions plays a critical role in correctly estimating atmospheric abundances. Finally, our results highlight the inadequacy of simplified semi-analytic models with isobaric assumptions for reliable retrievals of transmission spectra. Transmission spectra obtained with current facilities such as the Hubble Space Telescope and Very Large Telescope can provide strong constraints on atmospheric abundances of exoplanets.

VZ Lib is a southern triple-lined system. By analyzing all available times of light minima, orbital period changes are revisited in detail. We discovered that the observed–calculated (O − C) curve shows a long-term period decrease at a rate of ${dP}/{dt}=-2.25\times {10}^{-7}\,{\mathrm{days\; yr}}^{-1}$, revealing a mass transfer from the more massive component to the less massive one. A cyclic variation covering more than three cycles was discovered, which was analyzed for the light-travel-time effect via the presence of the tertiary companion. The cyclic variation has a short period of 2.96 (±0.04) yr and a small semiamplitude of 0.0039 (±0.0004) days. The mass of the third body was determined to be ${M}_{3}\sin {i}_{3}=0.52(\pm 0.07)$M_⊙ and an orbital semimajor axis shorter than $1.93(\pm 0.31)$ au was obtained. Orbital properties of this close-in companion should provide valuable information on the formation of close binaries and stellar dynamical interaction. New complete ${{BVR}}_{c}{I}_{c}$ light curves of VZ Lib were obtained and modeled with the Wilson–Devinney code. The light curves show a small but significant O'Connell effect that was not detectable in 1981 and 2007 but in 2004, so we derived a new photometric solution with assuming spot and a third light in the system. It is found that the light-curve subtype changed from A-type to W-type, which was possibly caused by a dark spot on the massive component. Our photometric solutions are in agreement with the spectroscopic results given by previous authors.

Of the solar system's four terrestrial planets, the origin of Mercury is perhaps the most mysterious. Modern numerical simulations designed to model the dynamics of terrestrial planet formation systematically fail to replicate Mercury, which possesses just 5% of the mass of Earth and the highest orbital eccentricity and inclination among the planets. However, Mercury's large iron-rich core and low volatile inventory stand out among the inner planets, and seem to imply a violent collisional origin. Because most algorithms used for simulating terrestrial accretion do not consider the effects of collisional fragmentation, it has been difficult to test these collisional hypotheses within the larger context of planet formation. Here, we analyze a large suite of terrestrial accretion models that account for the fragmentation of colliding bodies. We find that planets with core mass fractions boosted as a result of repeated hit-and-run collisions are produced in 90% of our simulations. While many of these planets are similar to Mercury in mass, they rarely lie on Mercury-like orbits. Furthermore, we perform an additional batch of simulations designed to specifically test the single giant impact origin scenario. We find less than a 1% probability of simultaneously replicating the Mercury–Venus dynamical spacing and the terrestrial system's degree of orbital excitation after such an event. While dynamical models have made great strides in understanding Mars' low mass, their inability to form accurate Mercury analogs remains a glaring problem.

The signature of wind patterns caused by the interplay of rotation and energy redistribution in hot Jupiters is detectable at high spectral resolution, yet no direct comparison has been attempted between predictions from general circulation models (GCMs) and observed high-resolution spectra. We present the first such comparison on near-infrared transmission spectra of the hot Jupiter HD 189733b. Exploring 12 rotation rates and two chemical regimes, we have created model spectra from 3D GCMs and cross-correlated them with the observed spectra. Comparing our models against those of HD 189733b, we obtain three key results: (1) we confirm CO and H₂O in the planet's atmosphere at a detection significance of 8.2σ; (2) we recover the signature of day-to-night winds with speeds of several km s⁻¹ at pressures of several millibars; and (3) we constrain the rotation period of the planet to between 1.2 and 4.69 days (synchronous rotation (2.2 days) remains consistent with existing observations). Our results do not suffer from the shortcomings of 1D models as cross-correlation templates—these models mainly tend to overconstrain the slower rotation rates and show evidence for anomalous blueshifts. Our 3D models instead match the observed line-of-sight velocity of this planet by self-consistently including the effects of high-altitude day-to-night winds. Overall, we find a high degree of consistency between observations of HD 189733b and our GCM-based spectra, implying that the physics and chemistry are adequately described in current 3D forward models for the purpose of interpreting observations at high spectral resolution.

(349) Dembowska is a big R-type main-belt asteroid. A new and more precise determination of its mass, M, is presented. The high precision is achieved, in particular, using our generally applied criterion of selecting model parameters. This criterion takes into account the sensitivity to M of the available observables, i.e., the position of another much smaller body involved in a close encounter with (349) Dembowska. Using this criterion, as well as observations from the Minor Planet Center (MPC) and the ephemeris DE431, M is redetermined to be (4.1 ± 0.4) × 10¹⁸ kg, which is a weighted mean of compatible masses independently determined with five encounters. From the most recently determined effective diameter, we get a bulk density of ρ = 2.08 ± 0.24 g cm⁻³ for (349) Dembowska. A macroporosity of 36 is derived from a grain density of 3.25 g cm⁻³.

The Transiting Exoplanet Survey Satellite (TESS) is conducting a two-year wide-field survey searching for transiting exoplanets around nearby bright stars that will be ideal for follow-up characterization. To facilitate studies of planet compositions and atmospheric properties, accurate and precise planetary radii need to be derived from the transit light curves. Since 40%–50% of exoplanet host stars are in multiple star systems, however, the observed transit depth may be diluted by the flux of a companion star, causing the radius of the planet to be underestimated. High angular resolution imaging can detect companion stars that are not resolved in the TESS Input Catalog, or by seeing-limited photometry, to validate exoplanet candidates and derive accurate planetary radii. We examine the population of stellar companions that will be detectable around TESS planet candidate host stars, and those that will remain undetected, by applying the detection limits of speckle imaging to the simulated host star populations of Sullivan et al. and Barclay et al. By detecting companions with contrasts of Δm ≲ 7–9 and separations of ∼002–12, speckle imaging can detect companion stars as faint as early M stars around A–F stars and stars as faint as mid-M around G–M stars, as well as up to 99% of the expected binary star distribution for systems located within a few hundred parsecs.

We measure near-ultraviolet (NUV) aperture magnitudes from Galaxy Evolution Explorer images for 258 ultra-diffuse galaxy (UDG) candidates drawn from the initial Systematically Measuring Ultra-Diffuse Galaxies (SMUDGes) survey of ∼300 square degrees surrounding, and including, the Coma galaxy cluster. For the vast majority, 242 of them, we present flux upper limits due either to a lack of significant flux in the aperture or confusion with other objects projected within the aperture. These limits often place interesting constraints on the UDG candidates, indicating that they are non-star-forming or quiescent. In particular, we identify field, quiescent UDG candidates, which are a challenge for formation models and are, therefore, compelling prospects for spectroscopic follow-up and distance determinations. We present far-ultraviolet (FUV) and NUV magnitudes for 16 detected UDG candidates and compare those galaxies to the local population of galaxies on color–magnitude and specific star formation rate diagrams. The NUV-detected UDG candidates form mostly an extension toward lower stellar masses of the star-forming galaxy sequence, and none of these lie within regions of high local galaxy density. UDG candidates span a range of properties, although almost all are consistent with being quiescent, low surface brightness galaxies, regardless of environment.

In the near-future, atmospheric characterization of Earth-like planets in the habitable zone will become possible via reflectance spectroscopy with future telescopes such as the proposed LUVOIR and HabEx missions. While previous studies have considered the effect of clouds on the reflectance spectra of Earth-like planets, the molecular detectability considering a wide range of cloud properties has not been previously explored in detail. In this study, we explore the effect of cloud altitude and coverage on the reflectance spectra of Earth-like planets at different geological epochs and examine the detectability of ${{\rm{O}}}_{2},{{\rm{H}}}_{2}{\rm{O}}$, and CH₄ with test parameters for the future mission concept, LUVOIR, using a coronagraph noise simulator previously designed for WFIRST-AFTA. Considering an Earth-like planet located at 5 pc away, we have found that for the proposed LUVOIR telescope, the detection of the O₂A-band feature (0.76 μm) will take approximately 100, 30, and 10 hr for the majority of the cloud parameter space modeled for the atmospheres with 10%, 50%, and 100% of modern Earth O₂ abundances, respectively. In particular, for the case of ≳50% of modern Earth O₂ abundance, the feature will be detectable with an integration time ≲10 hr as long as there are lower-altitude (≲8 km) clouds with a global coverage of ≳20%. For the 1% of the modern Earth O₂ abundance case, however, it will take more than 100 hr for all the cloud parameters we modeled.