Table of contents

Observations from the Kepler and K2 missions have provided the astronomical community with unprecedented amounts of data to search for transiting exoplanets and other astrophysical phenomena. Here, we present K2-288, a low-mass binary system (M2.0 ± 1.0; M3.0 ± 1.0) hosting a small (R_p = 1.9 R_⊕), temperate (T_eq = 226 K) planet observed in K2 Campaign 4. The candidate was first identified by citizen scientists using Exoplanet Explorers hosted on the Zooniverse platform. Follow-up observations and detailed analyses validate the planet and indicate that it likely orbits the secondary star on a 31.39-day period. This orbit places K2-288Bb in or near the habitable zone of its low-mass host star. K2-288Bb resides in a system with a unique architecture, as it orbits at >0.1 au from one component in a moderate separation binary (a_proj ∼ 55 au), and further follow-up may provide insight into its formation and evolution. Additionally, its estimated size straddles the observed gap in the planet radius distribution. Planets of this size occur less frequently and may be in a transient phase of radius evolution. K2-288 is the third transiting planet system identified by the Exoplanet Explorers program and its discovery exemplifies the value of citizen science in the era of Kepler, K2, and the Transiting Exoplanet Survey Satellite.

We report on Cosmic Origin Spectrograph observations of the gamma-ray bright blazar B2 1215+30, collected in 2015 November. These observations allow for the confirmation of the source redshift from the detection of a Lyα emission feature atλ ∼ 1374 Å. The emission feature places the source at a redshift ofz = 0.1305 ± 0.003, confirming the source's ground-based spectral measurement. The gamma-ray emission of the source is discussed in the context of the source distance, required for the accurate reconstruction of the intrinsic gamma-ray emission taking the absorption by the extragalactic background light into account. The source distance is found to be low enough that the previously reported detection of an exceptional flaring event from B2 1215+30 in 2014 cannot be used to investigate opacity-specific spectral and variability characteristics introduced by possible ultra-high-energy cosmic-ray propagation.

We consider the problem of inferring the shape of a transiting object's silhouette from its light curve alone, without assuming a physical model for the object. We model the object as a grid of pixels which transits a star; each pixel has an opacity, ranging from transparent to opaque, which we infer from the light curve. We explore three interesting degeneracies inherent to this problem, in which markedly different transiting shapes can produce identical light curves: (i) the "flip" degeneracy, by which two pixels transiting at the same impact parameter on opposite sides of the star's horizontal midplane generate the same light curve; (ii) the "arc" degeneracy, by which opacity can be redistributed along the semicircular arc of pixels which undergoes ingress or egress at the same time without consequence to the light curve; and (iii) the "stretch" degeneracy, by which a wide shape moving fast can produce the same light curve as a narrow shape moving more slowly. By understanding these degeneracies and adopting some additional assumptions, we are able to numerically recover informative shadow images of transiting objects, and we explore a number of different algorithmic approaches to this problem. We apply our methods to real data, including the TRAPPIST-1c/e/f triple transit and two dips of Boyajian's Star. We provide Python code to calculate the transit light curve of any grid and, conversely, infer the image grid which generates any light curve in the software package accompanying this paper, EightBitTransit (https://github.com/esandford/EightBitTransit).

We present the discovery of three new transiting hot Jupiters by the WASP-South project, WASP-161 b, WASP-163 b, and WASP-170 b. Follow-up radial velocities obtained with the Euler/CORALIE spectrograph and transit light curves obtained with the TRAPPIST-North, TRAPPIST-South, SPECULOOS-South, NITES, and Euler telescopes have enabled us to determine the masses and radii for these transiting exoplanets. WASP-161 b completes an orbit around its V = 11.1 F6V-type host star in 5.406 days, and has a mass M_p = 2.5 ± 0.2M_Jup and radius R_p = 1.14 ± 0.06 R_Jup. WASP-163 b orbits around its host star (spectral type G8V and the magnitude V = 12.5) every 1.609 days, and has a mass of M_P = 1.9 ± 0.2 M_Jup and a radius of R_p = 1.2 ± 0.1 R_Jup. WASP-170 b has a mass of 1.7 ± 0.2 M_Jup and a radius of 1.14 ± 0.09 R_Jup and is on a 2.344 day orbit around a G1V-type star of magnitude V = 12.8. Given their irradiations (∼10⁹ erg s⁻¹ cm⁻²) and masses, the three new planets' sizes are in good agreement with classical models of irradiated giant planets.

We present a near-infrared study of the Seyfert 2 galaxy NGC 6300, based on subarcsecond images and long slit spectroscopy obtained with Flamingos-2 at Gemini South. We have found that the peak of the nuclear continuum emission in the K_s band and the surrounding nuclear disk are 25 pc off-center with respect to the center of symmetry of the larger scale circumnuclear disk, suggesting that this black hole is still not fixed in the galaxy potential well. The molecular gas radial velocity curve yields a central black hole upper mass estimation of ${M}_{\mathrm{SMBH}}^{\mathrm{upper}}=(6.25\pm 2.5)\times {10}^{7}\,{M}_{\odot }$. The Paβ emission line has a strongly asymmetric profile with a blueshifted broad component that we associate with a nuclear ionized gas outflow. We have found in the K_s-band spectra that the slope of the continuum becomes steeper with increasing radii, which can be explained as the presence of large amounts of hot dust not only in the nucleus but also in the circumnuclear region up to r = 27 pc. In fact, the nuclear red excess obtained after subtracting the stellar contribution resembles to that of a blackbody with temperatures around 1200 K. This evidence supports the idea that absorbing material located around the nucleus, but not close enough to be the torus of the unified model, could be responsible for at least part of the nuclear obscuration in this Seyfert 2 nucleus.

We investigate via numerical modeling the effects of two planets locked in resonance, and migrating outward, on the dust distribution of the natal circumstellar disk. We aim to test whether the dust distribution exhibits peculiar features arising from the interplay among the gravitational perturbations of the planets in resonance, the evolution of the gas, and its influence on the dust grain dynamics. We focus on the 3:2 and 2:1 resonance, where the trapping may be caused by the convergent migration of a Jupiter- and Saturn-mass planet, preceding the common gap formation and ensuing outward (or inward) migration. Models show that a common gap also forms in the dust component similarly to what a single, more massive planet would generate and that outward migration leads to a progressive widening of the dust gap and to a decoupling from the gas gap. As the system evolves, a significantly wider gap is observed in the dust distribution, which ceases to overlap with the gas gap in the inner disk regions. At the outer edge of the gas gap, outward migration of the planets produces an overdensity of dust particles, which evolve differently in the 3:2 and 2:1 resonances. For the 3:2 resonance, the dust trap at the gap's outer edge is partly efficient, and a significant fraction of the grains filters through the gap. For the 2:1 resonance, the trap is more efficient, and very few grains cross the gap, while the vast majority accumulate at the outer edge of the gap.

We describe test observations made with a customized 640 × 512 pixel indium gallium arsenide (InGaAs) prototype astronomical camera on the 2.5 m DuPont telescope. This is the first test of InGaAs as a cost-effective alternative to HgCdTe for research-grade astronomical observations. The camera exhibits a background of 113 e^– s^–1/pixel (dark + thermal) at operating sensor temperature T = −40°C, maintained by a simple thermoelectric cooler. The optomechanical structure floats at ambient temperature with no cold stop, unlike most IR instruments which must be cooled to mitigate thermal backgrounds. Measurements of the night sky using a reimager scaled at 04/pixel show that the sky flux in Y is comparable to the dark current. At J the sky exceeds dark current by a factor of four, dominating the noise budget. The read noise (∼43 e⁻) falls below sky + dark noise for exposures of t > 7 s in Y and 3.5 s in J. We observed several representative science targets, including two supernovae, a z = 6.3 quasar, two local galaxies monitored for IR transients, and a galaxy cluster at z = 0.87. We observe a partial transit of the hot Jupiter HATS-34b, demonstrating the photometric stability to detect a 1.2% transit. A tiling of larger-format sensors would produce an IR survey instrument with significant cost savings relative to HgCdTe-based cameras, if one is willing to forego the K band. Such a camera would be sensitive for a week or more to isotropic emission from r-process kilonova ejecta similar to that observed in GW170817, over the full 190 Mpc horizon of Advanced LIGO.

We conducted a variable star search on the metal-rich Galactic globular cluster NGC 6652 using archival Gemini-S/Gemini Multiobject Spectrograph data. We report the discovery of nine new variable stars in the NGC 6652 field, of which we classify six as eclipsing binaries and one as an SX Phoenicis star, leaving two variables without classification. Using proper motions from Gaia DR2 and Hubble Space Telescope, albeit with some uncertainties, we find that the cluster, the field, and the background Sagittarius stream each have three of these variables. We also reassess the membership of known variables based on the Gaia proper motions, confirming the existence of one RR Lyrae star in the cluster.

Narrow-line Seyfert 1 galaxies (NLS1s) are a class of active galactic nuclei that are known to be one of the few sources of γ-rays, which originate in a relativistic beamed jet. Because of their relatively large distance, a poorly investigated aspect of these jetted NLS1s is their environment, and in particular, their host galaxy. In this work, we present the results of a morphological analysis of the host galaxy of the jetted NLS1 IRAS 20181-2244 observed with the 6.5 m Baade Telescope of the Las Campanas Observatory. The GALFIT analysis run on the Ks image, along with additional spectroscopic observations performed with the Nordic Optical Telescope, clearly revealed the presence of an interacting system of two galaxies. The data suggest that this NLS1 is hosted by a late-type galaxy, although the result is not conclusive. This analysis, along with other results in the literature, might suggest that two populations of jetted NLS1 exist. Further morphological studies are needed to confirm or disprove this hypothesis.

We report on the spread of [Fe/H] values in the massive Large Magellanic Cloud cluster NGC 1978, recently confirmed to harbor multiple populations of nearly the same age. We used accurate Strömgren photometry of carefully selected cluster red giant branch stars along with a high-dispersion spectroscopy-based calibration of the metallicity-sensitive index m₁. Once we accounted for the photometry quality, assessed from extensive artificial star tests to trace the photometric uncertainties as a function of the position of the cluster's center as well as the stellar brightness, and those from the metallicity calibration, we found that NGC 1978 exhibits a small metallicity spread of 0.035 dex (±0.019–0.023), depending on whether stars with individual σ[Fe/H] ≤ 0.15 dex or those located in the cluster's outer areas are considered. Such a spread in [Fe/H] is consistent with a cluster formation model with self-enrichment, if mass loss higher than 90% due to stellar evolutionary and galactic tidal effects is assumed. Nevertheless, scenarios in which the apparent [Fe/H] variation reflects CN abundance anomalies or less extreme mass-loss models with environmentally dependent self-enrichment should not be ruled out.

Using integral field data we extract the optical spectra of shocked interstellar clouds in Kepler's supernova remnant located in the inner regions of our Galaxy, as well as in the Large Magellanic Cloud, the Small Magellanic Cloud, NGC 6822, and IC 1613. Using self-consistent shock modeling, we make a new determination of the chemical composition of the interstellar medium in N, O, Ne, S, Cl, and Ar in these galaxies and obtain accurate estimates of the fraction of refractory grains destroyed in the shock. By comparing our derived abundances with those obtained in recent works using observations of B-stars, F supergiant stars, and H ii regions, we provide a new calibration for abundance scaling in the range of $7.9\lesssim 12+\mathrm{log}\,{\rm{O}}/{\rm{H}}\lesssim 9.1$.

We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial velocity (RV) measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest and most massive known to host a hot Jupiter. Based on the 27 days of TESS photometry and RV data from the CHIRON, HARPS, and Tillinghast Reflector Echelle Spectrograph, the planet has a mass of ${1.017}_{-0.068}^{+0.070}\,{M}_{{\rm{J}}}$ and radius of ${1.545}_{-0.060}^{+0.052}\,{R}_{{\rm{J}}}$, making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known.

We use radial velocity (RV) observations to search for long-period gas giant companions in systems hosting inner super-Earth (1–4 R_⊕, 1–10 M_⊕) planets to constrain formation and migration scenarios for this population. We consistently refit published RV data sets for 65 stars and find nine systems with statistically significant trends indicating the presence of an outer companion. We combine these RV data with AO images to constrain the masses and semi-major axes of these companions. We quantify our sensitivity to the presence of long-period companions by fitting the sample with a power-law distribution and find an occurrence rate of 39% ± 7% for companions 0.5–20 M_Jup and 1–20 au. Half of our systems were discovered by the transit method, and half were discovered by the RV method. While differences in the RV baselines and number of data points between the two samples lead to different sensitivities to distant companions, we find that occurrence rates of gas giant companions in each sample are consistent at the 0.5σ level. We compare the frequency of Jupiter analogs in these systems to the equivalent rate from field star surveys and find that Jupiter analogs are more common around stars hosting super-Earths. We conclude that the presence of outer gas giants does not suppress the formation of inner super-Earths, and that these two populations of planets instead appear to be correlated. We also find that the stellar metallicities of systems with gas giant companions are higher than those without companions, in agreement with the well-established metallicity correlation from RV surveys of field stars.

We report the discovery of two mid-infrared nebulae in the northern hemisphere with the Wide-field Infrared Survey Explorer and the results of optical spectroscopy of their central stars, BD+60° 2668 (composed of two components, separated from each other by ≈3 arcsec) and ALS 19653, with the Calar Alto 3.5 m telescope and the Southern African Large Telescope (SALT), respectively. We classify the components of BD+60° 2668 as stars of spectral types B0.5 II and B1.5 III. ALS 19653 is indicated in the SIMBAD database as a planetary nebula, while our observations show that it is a massive B0.5 Ib star, possibly in a binary system. Using the stellar atmosphere code fastwind, we derived fundamental parameters of the three stars as well as their surface element abundances, implying that all of them are either on the main sequence or only recently left it. This provides further evidence that massive stars can produce circumstellar nebulae while they are still relatively unevolved. We also report the detection of optical counterparts to the mid-infrared nebulae and a second, more extended optical nebula around ALS 19653, and present the results of SALT spectroscopy of both nebulae associated with this star. The possible origin of the nebulae is discussed.

Observations of active asteroid P/2017 S5 when near perihelion reveal the ejection of large (10²–10⁴μm) particles at 0.2–2 m s⁻¹ speeds, with estimated mass-loss rates of a few kg s⁻¹. The protracted nature of the mass loss (continuous over ∼150 days) is compatible with a sublimation origin, meaning that this object is likely an ice-bearing main-belt comet. Equilibrium sublimation of exposed water ice covering as little as 0.1 km² can match the data. Observations a year after perihelion show the object in an inactive state from which we deduce a nucleus effective radius ${450}_{-60}^{+100}\,{\rm{m}}$ (albedo 0.06 ± 0.02 assumed). The gravitational escape speed from a body of this size is just ∼0.3 m s⁻¹, comparable to the inferred ejection speed of the dust. Time-series photometry provides tentative evidence for rapid rotation (lightcurve period 1.4 hr) that may also play a role in the loss of mass and which, if real, is a likely consequence of spin-up by sublimation torques. P/2017 S5 shares both physical and orbital similarities with the split active asteroid pair P/2016 J1-A and J1-B, and all three objects are likely members of the ∼7 Myr old, collisionally produced, Theobalda family.

We report the discovery of 10 transiting extrasolar planets by the HATSouth survey. The planets range in mass from the super-Neptune HATS-62b, with ${M}_{p}\lt 0.179$${M}_{{\rm{J}}}$, to the super-Jupiter HATS-66b, with ${M}_{p}=5.33$${M}_{{\rm{J}}}$, and in size from the Saturn HATS-69b, with ${R}_{p}=0.94$${R}_{{\rm{J}}}$, to the inflated Jupiter HATS-67b, with ${R}_{p}=1.69$${R}_{{\rm{J}}}$. The planets have orbital periods between $1.6092$ days (HATS-67b) and $7.8180$ days (HATS-61b). The hosts are dwarf stars with masses ranging from $0.89$${M}_{\odot }$ (HATS-69) to $1.56$${M}_{\odot }$ (HATS-64) and have apparent magnitudes between $V=12.276\pm 0.020$ mag (HATS-68) and $V=14.095\pm 0.030$ mag (HATS-66). The super-Neptune HATS-62b is the least massive planet discovered to date with a radius larger than Jupiter. Based largely on the Gaia DR2 distances and broadband photometry, we identify three systems (HATS-62, HATS-64, and HATS-65) as having possible unresolved binary star companions. We discuss in detail our methods for incorporating the Gaia DR2 observations into our modeling of the system parameters and into our blend analysis procedures.

We present 248 speckle observations of 43 binary and 19 trinary star systems chosen to make progress in two main areas of investigation: the fundamental properties of metal-poor stars and star formation mechanisms. The observations were taken at the Gemini North and South telescopes during the period 2015 July to 2018 April, mainly with the Differential Speckle Survey Instrument but also with a few early results from the new 'Alopeke speckle camera at Gemini North. We find that the astrometry and photometry of these observations as a whole are consistent with previous work at Gemini. We present five new visual orbits for systems important in understanding metal-poor stars, three of which have orbital periods of less than 4 yr, and we indicate the degree to which these and future observations can impact our knowledge of stellar properties and star formation. In particular, we find a decrease in mass at fixed spectral type for metal-poor stars versus their solar-metallicity analogs that is consistent with predictions that are made from current stellar models.

We present high spatial resolution LBTI/NOMIC 9–12 μm images of VY CMa and its massive outflow feature, the Southwest (SW) Clump. Combined with high-resolution imaging from the Hubble Space Telescope (0.4–1 μm) and LBT/LMIRCam (1–5μm), we isolate the spectral energy distribution (SED) of the clump from the star itself. Using radiative-transfer code DUSTY, we model both the scattered light from VY CMa and the thermal emission from the dust in the clump to estimate the optical depth, mass, and temperature of the SW Clump. The SW Clump is optically thick at 8.9 μm with a brightness temperature of ∼200 K. With a dust chemistry of equal parts silicates and metallic iron, as well as assumptions on grain size distribution, we estimate a dust mass of 5.4 × 10⁻⁵M_⊙. For a gas-to-dust ratio of 100, this implies a total mass of 5.4 × 10⁻³M_⊙. Compared to the typical mass-loss rate of VY CMa, the SW Clump represents an extreme, localized mass-loss event from ≲300 yr ago.

We present high-resolution optical transmission spectroscopy of two sub-Saturn mass transiting exoplanets, HAT-P-12b and WASP-69b. With relatively low densities and high atmospheric scale heights, these planets are particularly well-suited to characterization through transit spectroscopy, and serve as ideal candidates for extending previously tested methods to lower planetary masses. Using a single transit for each planet, we take advantage of the Doppler cross-correlation technique to search for sodium, potassium, and water absorption features. Our analysis reveals a likely (3.2σ) detection of sodium absorption features in the atmosphere of HAT-P-12b, and enables us to place constraints on the presence of alkaline and molecular species in the atmospheres of both planets. With our results, we highlight the efficacy of ground-based campaigns for characterizing exoplanetary atmospheres and pave the way for future analyses of low-mass planets.

A relatively massive and moderately eccentric disk of trans-Neptunian objects (TNOs) can effectively counteract apse precession induced by the outer planets, and in the process shepherd highly eccentric members of its population into nearly stationary configurations that are antialigned with the disk itself. We were sufficiently intrigued by this remarkable feature to embark on an extensive exploration of the full spatial dynamics sustained by the combined action of giant planets and a massive trans-Neptunian debris disk. In the process, we identified ranges of disk mass, eccentricity, and precession rate that allow apse-clustered populations that faithfully reproduce key orbital properties of the much-discussed TNO population. The shepherding disk hypothesis is, to be sure, complementary to any potential ninth member of the solar system pantheon, and could obviate the need for it altogether. We discuss its essential ingredients in the context of solar system formation and evolution, and argue for their naturalness in view of the growing body of observational and theoretical knowledge about self-gravitating disks around massive bodies, extra-solar debris disks included.

We present an analytic description of tides raised on a star by a small orbiting body. In particular, we highlight the disproportionate effect of eccentricity and thus the scope for using these tides to detect and characterize the orbits of exoplanets and brown dwarfs. The tidal distortions of the star produced by an eccentric orbit are, in comparison to a circular orbit, much richer in detail and potentially visible from any viewing angle. The magnitude of these variations is much larger than that in a circular orbit of the same semimajor axis. These variations are visible in both photometric and spectroscopic data and dominate other regular sources of phase variability (e.g., reflection and Doppler beaming) over a particularly interesting portion of parameter space. These tidal signatures will be a useful tool for planet detection on their own and, used in concert with other methods, provide powerful constraints on planetary and stellar properties.

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

The most distant Kuiper Belt objects (KBOs) appear to be clustered in longitude of perihelion and in orbital pole position. To date, the only two suggestions for the cause of these apparent clusterings have been either the effects of observational bias or the existence of a distant giant planet in an eccentric inclined orbit known as Planet Nine. To determine if observational bias can be the cause of these apparent clusterings, we develop a rigorous method of quantifying the observational biases in the observations of longitude of perihelion and orbital pole position. From this now more complete understanding of the biases, we calculate that the probability that these distant KBOs would be clustered as strongly as observed in both longitude of perihelion and in orbital pole position is only 0.2%. While explanations other than Planet Nine may someday be found, the statistical significance of this clustering is now difficult to discount.

M subdwarfs are low-metallicity M dwarfs that typically inhabit the halo population of the Galaxy. Metallicity controls the opacity of stellar atmospheres; in metal-poor stars, hydrostatic equilibrium is reached at a smaller radius, leading to smaller radii for a given effective temperature. We compile a sample of 88 stars that span spectral classes K7 to M6 and include stars with metallicity classes from solar-metallicity dwarf stars to the lowest metallicity ultra subdwarfs to test how metallicity changes the stellar radius. We fit models to Palomar Double Spectrograph (DBSP) optical spectra to derive effective temperatures (T_eff) and we measure bolometric luminosities (L_bol) by combining broad wavelength-coverage photometry with Gaia parallaxes. Radii are then computed by combining the T_eff and L_bol using the Stefan–Boltzman law. We find that for a given temperature, ultra subdwarfs can be as much as five times smaller than their solar-metallicity counterparts. We present color-radius and color-surface brightness relations that extend down to [Fe/H] of −2.0 dex, in order to aid the radius determination of M subdwarfs, which will be especially important for the WFIRST exoplanetary microlensing survey.

We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon, or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet–planet/planet–moon occultation light curves, as well as thermal (rotational) phase curves. In this paper, we also introduce starry, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in starry is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. starry also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map. (Please see https://github.com/rodluger/starry, https://rodluger.github.io/starry/, and https://doi.org/10.5281/zenodo.1312286).

The distant long-period comet C/2017 K2 (PANSTARRS) has been outside the planetary region of the solar system for ∼3 Myr, negating the possibility that heat retained from the previous perihelion could be responsible for its activity. This inbound comet is also too cold for water ice to sublimate and too cold for amorphous water ice, if present, to crystallize. C/2017 K2 thus presents an ideal target in which to investigate the mechanisms responsible for activity in distant comets. We have used the Hubble Space Telescope to study the comet in the pre-perihelion heliocentric distance range 13.8 ≤ r_H ≤ 15.9 au. In this range, the coma maintains a logarithmic surface brightness gradient m = −1.010 ± 0.004, consistent with mass loss proceeding in steady state. The absence of a radiation pressure swept tail indicates that the effective particle size is large (radius ≳0.1 mm) and the mass-loss rate is ∼200 kg s⁻¹, remarkable for a comet still beyond the orbit of Saturn. Extrapolation of the photometry indicates that activity began in 2012.1 ± 0.5, at r_H = 25.9 ± 0.9 au, where the isothermal blackbody temperature is only T_BB = 55 K. This large distance and low temperature suggest that cometary activity is driven by the sublimation of a super-volatile ice (e.g., CO), presumably preserved by K2's long-term residence in the Oort cloud. The mass-loss rate can be sustained by CO sublimation from an area ≲2 km², if located near the hot subsolar point on the nucleus. However, while the drag force from sublimated CO is sufficient to lift millimeter-sized particles against the gravity of the cometary nucleus, it is 10²–10³ times too small to eject these particles against interparticle cohesion. Our observations thus require either a new understanding of the physics of interparticle cohesion or the introduction of another mechanism to drive distant cometary mass loss. We suggest thermal fracture and electrostatic supercharging in this context.

We have identified a set of optical emission-line features 700'' (12 kpc) to the southwest of the nucleus of Centaurus A, roughly opposite to the radio jet and well-known optical emission filaments associated with the northern radio structure. This location is roughly along the axis of the southwestern radio lobes, although there is no coherent jet at this radius. We use integral-field optical spectroscopy to trace the ratios of strong emission lines, showing changes in excitation across the region, and significant local reddening. The emission regions are spatially associated with far-infrared emission peaks in one of two cold dust clouds identified using Herschel and Spitzer data, and there may be a mismatch between the low temperature of the dust and the expected heating effect of young stars. The strong emission lines have ratios consistent with photoionization in normal H ii regions, requiring only modest numbers of OB stars; these stars and their cooler accompanying populations must be largely obscured along our line of sight. These data fit with a picture of fairly ordinary formation of clusters in a large giant molecular cloud, or network of such clouds. The location, projected near the radio-source axis and within the radius where a starburst wind has been inferred on the other side of the galaxy, raises the question of whether this star-forming episode was enhanced or indeed triggered by an outflow from the central parts of Centaurus A. However, the level of star formation is normal for the associated cold-gas mass and column density, and optical emission-line ratios and line widths limit the role of shocks on the gas, so any interaction with an outflow, associated either with the radio source or star formation in the gas-rich disk of Centaurus A, can at most have compressed the gas weakly. We speculate that the presence of similar star-forming regions on both sides of the galaxy, contrasted with the difference in the character of the emission-line clouds, reflects the presence of a collimated radio jet to the northeast and perhaps anisotropic escape of ionizing radiation from the AGN as well. In this view, the star formation on the southwestern side of Cen A could be enhanced indirectly via compression by a broad outflow (whether originated by a starburst or AGN), distinct from the radio jet and lobes.

Modern studies of the early solar system routinely invoke the possibility of an orbital instability among the giant planets triggered by gravitational interactions between the planets and a massive exterior disk of planetesimals. Previous works have suggested that this instability can be substantially delayed (about hundreds of megayears) after the formation of the giant planets. Bodies in the disk are typically treated in a semi-active manner, wherein their gravitational force on the planets is included, but interactions between the planetesimals are ignored. We perform N-body numerical simulations using GENGA, which makes use of GPUs to allow for the inclusion of all gravitational interactions between bodies. Although our simulated Kuiper Belt particles are more massive than the probable masses of real primordial Kuiper Belt objects, our simulations indicate that the self-stirring of the primordial Kuiper Belt is very important to the dynamics of the giant planet instability. We find that interactions between planetesimals dynamically heat the disk and typically prevent the outer solar system's instability from being delayed by more than a few tens of megayears after giant planet formation. Longer delays occur in a small fraction of systems that have at least 3.5 au gaps between the planets and planetesimal disk. Our final planetary configurations match the solar system at a rate consistent with other previous works in most regards. Pre-instability heating of the disk typically yields final Jovian eccentricities comparable to the modern solar system's value, which has been a difficult constraint to match in past works.

We present an optical transmission spectrum of the atmosphere of WASP-4b obtained through observations of four transits with Magellan/IMACS, as part of the Arizona-CfA-Católica-Carnegie Exoplanet Spectroscopy Survey (ACCESS). Using a Bayesian approach to atmospheric retrieval, we find no evidence for scattering or absorption features in our transit spectrum. Our models include a component to model the transit light source effect (spectral contamination from unocculted spots on the stellar photosphere), which we show can have a marked impact on the observed transmission spectrum for reasonable spot-covering fractions (<5%); this is the first such analysis for WASP-4b. We are also able to fit for the size and temperature contrast of spots observed during the second and third transits, finding evidence for both small, cool and large, warm spot-like features on the photosphere. Finally, we compare our results to those published by Huitson et al. using Gemini/GMOS and May et al. using IMACS, and we find that our data are in agreement.

We present a spectrally and temporally resolved detection of the optical Mg i triplet at 7.8σ in the extended atmosphere of the ultra-hot Jupiter KELT-9 b, adding to the list of detected metal species in the hottest gas giant currently known. Constraints are placed on the density and radial extent of the excited hydrogen envelope using simultaneous observations of Hα and Hβ under the assumption of a spherically symmetric atmosphere. We find that planetary rotational broadening of ${v}_{\mathrm{rot}}={8.2}_{-0.7}^{+0.6}$ km s⁻¹ is necessary to reproduce the Balmer line transmission profile shapes, where the model including rotation is strongly preferred over the non-rotating model using a Bayesian information criterion comparison. The time series of both metal line and hydrogen absorption show remarkable structure, suggesting that the atmosphere observed during this transit is dynamic rather than static. We detect a relative emission feature near the end of the transit which exhibits a P-Cygni-like shape, evidence of material moving at ≈50–100 km s⁻¹ away from the planet. We hypothesize that the in-transit variability and subsequent P-Cygni-like profiles are due to a flaring event that caused the atmosphere to expand, resulting in unbound material being accelerated to high speeds by stellar radiation pressure. Further spectroscopic transit observations will help establish the frequency of such events.

Sodium-rich stars are often found in globular clusters giants. However, some sodium-rich stars have been found among field metal-poor stars. These stars are considered as evaporated from globular clusters. Identified such kind of stars among the field stars in the Galaxy may provide insights of which mechanism was responsible for the ejection from a globular cluster and may reveal some chemical peculiarity. Therefore, we started a search, using high-resolution spectroscopy, among metal-poor stars from several sources of the literature to find a sodium-rich star. Here we present the results for the temperature, gravity, metallicity, and sodium abundances for the stars of our sample. For many of them we determined the temperature, gravity, metallicity, and sodium abundances for the first time. As a result of our search we found one star, CD-23°16310, which has a [Na/Fe] ratio of +1.09. We also show that CD-23°16310 is not a carbon-enhanced metal-poor (CEMP) star since carbon is not enriched but is nitrogen-rich. We did not detect any variation of the radial velocity that would support the hypothesis of mass transfer. Thus, the high sodium and nitrogen abundance could be due to a strong internal mixing process, suggesting that CD-23°16310 is an early asymptotic giant branch star.

Direct-imaging exoplanet surveys have discovered a class of 5–20 ${M}_{\mathrm{Jup}}$ substellar companions at separations >100 au from their host stars, which present a challenge to planet and star formation models. Detailed analysis of the orbital architecture of these systems can provide constraints on possible formation mechanisms, including the possibility that they were dynamically ejected onto a wide orbit. We present astrometry for the wide planetary-mass companion GSC 6214-210 b (240 au; ≈14 ${M}_{\mathrm{Jup}}$) obtained using NIRC2 with adaptive optics at the Keck telescope over 10 years. Our measurements achieved astrometric uncertainties of ≈1 mas per epoch. We determined a relative motion of 1.12 ± 0.15 mas yr⁻¹ (0.61 ± 0.09 km s⁻¹), the first detection of orbital motion for this companion. We compute the minimum periastron for the companion due to our measured velocity vector and derive constraints on the orbital parameters through our modified implementation of the Orbits for the Impatient rejection sampling algorithm. We find that close periastron orbits, which could indicate that the companion was dynamically scattered, are present in our posterior but have low likelihoods. For all orbits in our posterior, we assess the detectability of close-in companions that could have scattered GSC 6214-210 b from a closer orbit, and find that most potential scatterers would have been detected in previous imaging. We conclude that formation at small orbital separation and subsequent dynamical scattering through interaction with another potential close-in object is an unlikely formation pathway for this companion. We also update stellar and substellar properties for the system using the new parallax from Gaia DR2.

We report the discovery of a low-mass-ratio planet (q = 1.3 × 10⁻⁴), i.e., 2.5 times higher than the Neptune/Sun ratio. The planetary system was discovered from the analysis of the KMT-2017-BLG-0165 microlensing event, which has an obvious short-term deviation from the underlying light curve produced by the host of the planet. Although the fit improvement with the microlens parallax effect is relatively low, one component of the parallax vector is strongly constrained from the light curve, making it possible to narrow down the uncertainties of the lens physical properties. A Bayesian analysis yields that the planet has a super-Neptune mass $({M}_{2}={34}_{-12}^{+15}\,{M}_{\oplus })$ orbiting a Sun-like star $({M}_{1}={0.76}_{-0.27}^{+0.34}\,{M}_{\odot })$ located at 4.5 kpc. The blended light is consistent with these host properties. The projected planet-host separation is ${a}_{\perp }={3.45}_{-0.95}^{+0.98}\,\mathrm{au}$, implying that the planet is located outside the snow line of the host, i.e., a_sl ∼ 2.1 au. KMT-2017-BLG-0165Lb is the sixteenth microlensing planet with mass ratio q < 3 × 10⁻⁴. Using the fifteen of these planets with unambiguous mass-ratio measurements, we apply a likelihood analysis to investigate the form of the mass-ratio function in this regime. If we adopt a broken power law for the form of this function, then the break is at q_br ≃ 0.55 × 10⁻⁴, which is much lower than previously estimated. Moreover, the change of the power-law slope, ζ > 3.3, is quite severe. Alternatively, the distribution is also suggestive of a pileup of planets at Neptune-like mass ratios, below which there is a dramatic drop in frequency.

Photometric and spectral observations of the W UMa binaries NSVS 4161544 and 1SWASP J034501.24+493659.9 with periods of around 8 hr are presented. The simultaneous light-curve and radial-velocity-curve solutions revealed that the two targets have over-contact configurations, and their components undergo partial eclipses. The derived parameters of NSVS 4161544 are: mass ratio q = 3.377, orbital inclination i = 65.0⁰, temperatures T₁ = 6016 K and T₂ = 5726 K, masses M₁ = 0.43 ${M}_{\odot }$ and M₂ = 1.47 ${M}_{\odot },$ radii R₁ = 0.74 ${R}_{\odot }$ and R₂ = 1.28 ${R}_{\odot },$ luminosities L₁ = 0.64 ${L}_{\odot }$ and L₂ = 1.51 ${L}_{\odot }$. The derived parameters of 1SWASP J034501.24+493659.9 are: mass ratio q = 2.378, orbital inclination i = 60.3⁰, temperatures T₁ = 6514 K and T₂ = 6494 K, masses M₁ = 0.275 ${M}_{\odot }$ and M₂ = 0.654 ${M}_{\odot },$ radii R₁ = 0.693 ${R}_{\odot }$ and R₂ = 1.012 ${R}_{\odot },$ luminosities L₁ = 0.775 ${L}_{\odot }$ and L₂ = 1.632 ${L}_{\odot }$. Hence, the masses of the 1SWASP J034501.24+493659.9 components are quite small for the rest global parameters, and their sum of 0.93 ${M}_{\odot }$ is slightly below the lower mass limit of 1.0–1.2 M_⊙ for the contact binaries. The calculated distance of NSVS 4161544 almost coincides with that determined by GAIA. The calculated distance of 1SWASP J034501.24+493659.9 is considerably smaller than the GAIA value that is most likely due to overestimation of the interstellar extinction to this star from the Galactic disk.

We present the discovery of Qatar-7b—a very hot and inflated giant gas planet orbiting close to its parent star. The host star is a relatively massive main-sequence F-star with mass and radius $\,{M}_{\star }=1.41\pm 0.03\,{M}_{\odot }$ and $\,{R}_{\star }=1.56\pm 0.02\,{R}_{\odot }$, respectively, at a distance d = 726 ± 26 pc, and an estimated age ∼1 Gyr. With its orbital period of P = 2.032 days, the planet is located less than five stellar radii from its host star and is heated to a high temperature T_eq ≈ 2100 K. From a global solution to the available photometric and radial velocity observations, we calculate the mass and radius of the planet to be $\,{M}_{{\rm{P}}}$ = 1.88 ± 0.25 $\,{M}_{{\rm{J}}}$ and $\,{R}_{{\rm{P}}}$ = 1.70 ± 0.03 $\,{R}_{{\rm{J}}}$, respectively. The planet radius and equilibrium temperature put Qatar-7b in the top 6% of the hottest and largest known exoplanets. With its large radius and high temperature, Qatar-7b is a valuable addition to the short list of targets that offer the best opportunity for studying their atmospheres through transmission spectroscopy.

OSIRIS is a near-infrared (1.0–2.4 μm) integral field spectrograph operating behind the adaptive optics system at Keck Observatory and one of the first lenslet-based integral field spectrographs. Since its commissioning in 2005, it has been a productive instrument, producing nearly half the laser guide star adaptive optics papers on Keck. The complexity of its raw data format necessitated a custom data reduction pipeline (DRP) delivered with the instrument in order to iteratively assign flux in overlapping spectra to the proper spatial and spectral locations in a data cube. Other than bug fixes and updates required for hardware upgrades, the bulk of the DRP has not been updated since initial instrument commissioning. We report on the first major comprehensive characterization of the DRP using on-sky and calibration data. We also detail improvements to the DRP, including characterization of the flux assignment algorithm, exploration of spatial rippling in the reduced data cubes, and improvements to several calibration files, including the rectification matrix, bad-pixel mask, and wavelength solution. We present lessons learned from over a decade of OSIRIS data reduction that are relevant to the next generation of integral field spectrograph hardware and data reduction software design.

We present Hubble Space Telescope (HST) observations of the low surface brightness (SB) galaxy Coma P. This system was first discovered in the Arecibo Legacy Fast ALFA H i survey and was cataloged as an (almost) dark galaxy because it did not exhibit any obvious optical counterpart in the available survey data (e.g., Sloan Digital Sky Survey). Subsequent WIYN pODI imaging revealed an ultra-low SB stellar component located at the center of the H i detection. We use the HST images to produce a deep color–magnitude diagram of the resolved stellar population present in Coma P. We clearly detect a red stellar sequence that we interpret to be a red giant branch and use it to infer a tip of the red giant branch distance of ${5.50}_{-0.53}^{+0.28}$ Mpc. The new distance is substantially lower than earlier estimates and shows that Coma P is an extreme dwarf galaxy. Our derived stellar mass is only 4.3 × 10⁵M_⊙, meaning that Coma P has an extreme H i-to-stellar mass ratio of 81. We present a detailed analysis of the galaxy environment within which Coma P resides. We hypothesize that Coma P formed within a local void and has spent most of its lifetime in a low-density environment. Over time, the gravitational attraction of the galaxies located in the void wall has moved it to the edge, where it had a recent "fly-by" interaction with M64. We investigate the possibility that Coma P is at a farther distance and conclude that the available data are best fit by a distance of 5.5 Mpc.

We exploit high-quality photometry from the EVEREST pipeline to evaluate false-positive exoplanet candidates from the K2 mission. We compare the practical capabilities of EVEREST's pixel-level decorrelation scheme to the data analysis pipelines widely used at the time of these planet candidates' discovery. Removing stellar variability from the EVEREST-corrected light curves, we search for potential secondary eclipses. For each object exhibiting a secondary eclipse, we compare the implied brightness temperature of the planet candidate to its calculated equilibrium temperature. We thereby identify objects whose brightness temperature is too high to be consistent with a planet. We identify seven systems previously flagged as planetary candidates in preliminary vetting pipelines, and use EVEREST to instead identify six of them as eclipsing binaries. We also project the importance of optimal photometric vetting for TESS data. We find that the majority of blended eclipsing binaries could be identified using TESS photometry, and a systematic study of that kind could in principle also yield valuable information on the mass ratio distribution in stellar eclipsing binaries.

Binary and multiple stars have long provided an effective empirical method of testing stellar formation and evolution theories. In particular, the existence of wide binary systems (separations >20,000 au) is particularly challenging to binary formation models as their physical separations are beyond the typical size of a collapsing cloud core (∼5000–10,000 au). We mined the recently published Gaia-DR2 catalog to identify bright comoving systems in the five-dimensional space (sky position, parallax, and proper motion). We identified 3741 comoving binary and multiple stellar candidate systems, out of which 575 have compatible radial velocities for all the members of the system. The candidate systems have separations between ∼400 and 500,000 au. We used the analysis tools of the Virtual Observatory to characterize the comoving system members and to assess their reliability. The comparison with previous comoving systems catalogs obtained from TGAS showed that these catalogs contain a large number of false systems. In addition, we were not able to confirm the ultra-wide binary population presented in these catalogs. The robustness of our methodology is demonstrated by the identification of well known comoving star clusters and by the low contamination rate for comoving binary systems with projected physical separations <50,000 au. These last constitute a reliable sample for further studies. The catalog is available online at the Spanish Virtual Observatory portal (http://svo2.cab.inta-csic.es/vocats/v2/comovingGaiaDR2/).

Using a large suite of n-body simulations, we explore the discovery space for new satellites in the Pluto–Charon system. For the adopted masses and orbits of the known satellites, there are few stable prograde or polar orbits with semimajor axes $a\lesssim 1.1\,{a}_{H}$, where a_H is the semimajor axis of the outermost moon Hydra. Small moons with radii $r$ ≲ 2 km and a ≲ 1.1 a_H are ejected on timescales ranging from several years to more than 100 Myr. Orbits with a ≳ 1.1 a_H are stable on timescales exceeding 150–300 Myr. Near-infrared (IR) and mid-IR imaging with several instruments on James Webb Space Telescope and ground-based occultation campaigns with 2–3 m class telescopes can detect 1–2 km satellites outside the orbit of Hydra. Searches for these moons enable new constraints on the masses of the known satellites and on theories for circumbinary satellite formation.

Atmospheric nitrogen may be a necessary ingredient for the habitability of a planet as its presence helps to prevent water loss from a planet. The present-day nitrogen isotopic ratio, ¹⁵N/¹⁴N, in the Earth's atmosphere is a combination of the primitive Earth's ratio and the ratio that might have been delivered in comets and asteroids. Asteroids have a nitrogen isotopic ratio that is close to the Earth's. This indicates either a similar formation environment to the Earth or that the main source of nitrogen was delivery by asteroids. However, according to geological records, the Earth's atmosphere could have been enriched in ¹⁵N during the Archean era. Comets have a higher ¹⁵N/¹⁴N ratio than the current atmosphere of the Earth, and we find that about 5% ∼ 10% of nitrogen in the atmosphere of the Earth may have been delivered by comets to explain the current atmosphere of the Earth or the enriched ¹⁵N atmosphere of the Earth. We model the evolution of the radii of the snow lines of molecular nitrogen and ammonia in a protoplanetary disk and find that both have radii that put them farther from the Sun than the main asteroid belt. With an analytic secular resonance model and N-body simulations we find that the ν₈ apsidal precession secular resonance with Neptune, which is located in the Kuiper Belt, is a likely origin for the nitrogen-delivering comets that impact the Earth.

We report a multi-objective campaign of targeted 21 cm H i line observations of sources selected from the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey and galaxies identified by their morphological and photometric properties in the Sloan Digital Sky Survey. The aims of this program have been (1) to confirm the reality of some ALFALFA sources whose enigmatic nature suggest additional multiwavelength observations; (2) to probe the low signal-to-noise ratio (S/N) regime, below the ALFALFA reliability limit; and (3) to explore the feasibility of using optical morphology, color, and surface brightness to identify gas-rich objects in the region of the Pisces–Perseus Supercluster (PPS) whose H i fluxes are below the ALFALFA sensitivity limit at that distance. As expected, the reliability of ALFALFA detections depends strongly on the S/N of the H i line signal and its coincidence with a probable stellar counterpart identified by its optical properties, suggestive of ongoing star formation. The identification of low-mass, star-forming populations enables targeted H i line observations to detect galaxies with H i line fluxes below the ALFALFA sensitivity limits in fixed local volumes (D < 100 Mpc). The method explored here serves as the basis for extending the sample of gas-bearing objects as part of the ongoing Arecibo Pisces–Perseus Supercluster Survey (APPSS).

During the past five years, 6, 7, and 26 transit observations were carried out for the HAT-P-9b, HAT-P-32b, and HAT-P-36b systems, respectively, through the Transiting Exoplanet Monitoring Project network. Combined with the published photometric data and radial-velocity measurements, our new photometry allows us to revisit the system parameters and search for additional close-in planetary companions in these hot Jupiter systems. We measure an updated R_P/R_* = 0.1260 ± 0.0011 for HAT-P-36 system in the R band, which is 4.5σ larger than the published i-band radius ratio of 0.1186 ± 0.0012. We also perform a transit timing variation (TTV) analysis for each system. Because no significant TTVs were found, we place an upper mass limit on an additional planet for each system.

The presence of dust grains profoundly affects the diffusion of the magnetic field in molecular clouds. When the electrons and ions are well coupled to the magnetic field and charged grains are only indirectly coupled, emergent Hall diffusion may dominate over all the other non-ideal magnetohydrodynamic (MHD) effects in a partially ionized dusty cloud. The low-frequency, long (∼0.01–1 pc) wavelength dispersive MHD waves will propagate in such a medium with the polarization of the waves determined by the dust charge density or the dust size distribution. In the presence of shear flows, these waves may become Kelvin–Helmholtz unstable with the dust charge density or the grain size distribution operating as a switch to the instability. When Hall diffusion time is long (compared to the time over which waves are sheared), the growth rate of the instability in the presence of sub-Alfvénic flow increases with the charge number $| Z|$ on the grain, while it is quenched in the presence of Alfvénic or super-Alfvénic flows. However, when Hall diffusion is fast, the growth rate of the instability depends on the dust charge only indirectly.

We present a semianalytic estimate of the expected yield of single-transit planets from the Transiting Exoplanet Survey Satellite (TESS). We use the TESS Candidate Target List-6 (CTL-6) as an input catalog of over four million sources. We predict that from the 200,000 stars selected to be observed with the high-cadence postage stamps (PSs) with the highest CTL-6 priority, there will be 241 single-transit events caused by planets detectable at a signal-to-noise ratio of S/N ≥ 7.3. We find a lower limit of an additional 977 events caused by single-transit planets in the full frame images (FFIs); this is a lower limit because the CTL-6 is incomplete below a TESS magnitude of $T\gt 12$. Of the single-transit events from the PSs and FFIs, 1091/1218 will have transit depths deeper than 0.1% and will thus be amenable for photometric follow-up from the ground, and 1195/1218 will have radial velocity signals greater than 1 m s⁻¹. We estimate that the periods of 146 single transits will be constrained to better than 10% using the TESS photometry assuming circular orbits. We find that the number of planets detected by TESS in the PSs with periods $P\gt 25$ days can be doubled by including single-transiting planets, while the number of planets with $P\gt 250$ days can be increased by an order of magnitude. We predict 79 habitable zone planets in the TESS light curves from single transits, with 18 orbiting FGK stars.

We present results of our large-scale, optical, multi-epoch photometric survey across ∼180 square degrees in the Orion OB1 association, complemented with extensive follow-up spectroscopy. Our focus is mapping and characterizing the off-cloud, low-mass, pre-main-sequence (PMS) populations. We report 2062 K- and M-type confirmed T Tauri members; 59% are located in the OB1a subassociation, 27% in the OB1b subassociation, and the remaining 14% in the A and B molecular clouds. We characterize two new clusterings of T Tauri stars, the HD 35762 and HR 1833 groups, both located in OB1a not far from the 25 Ori cluster. We also identify two stellar overdensities in OB1b, containing 231 PMS stars, and find that the OB1b region is composed of two populations at different distances, possibly due to the OB1a subassociation overlapping with the front of OB1b. A ∼2 deg wide halo of young stars surrounds the Orion Nebula Cluster, corresponding in part to the low-mass populations of NGC 1977 and NGC 1980. We use the strength of Hα in emission, combined with the IR excess and optical variability, to define a new type of T Tauri star, the C/W class, stars we propose may be nearing the end of their accretion phase, in an evolutionary state between classical and weak-lined T Tauri stars. The evolution of the ensemble-wide equivalent width of Li iλ6707 indicates a Li depletion timescale of ∼8.5 Myr. Disk accretion declines with an e-folding timescale of ∼2 Myr, consistent with previous studies.

1I/'Oumuamua is the first detected interstellar interloper. We test the hypothesis that it is representative of a background population of exo-Oort cloud objects ejected under the effect of post-main sequence mass loss and stellar encounters. We do this by comparing the cumulative number density of interstellar objects inferred from the detection of 1I/'Oumuamua to that expected from these two clearing processes. We consider the 0.08–8 M_⊙ mass range, take into account the dependencies with stellar mass, Galactocentric distance, and evolutionary state, and consider a wide range of size distributions for the ejected objects. Our conclusion is that 1I/'Oumuamua is likely not representative of this background population, even though there are large uncertainties in the masses and size distributions of the exo-Oort Clouds. We discuss whether the number density of free-floating, planetary-mass objects derived from gravitational microlensing surveys could be used as a discriminating measurement regarding 1I/'Oumuamua's origin (given their potential common origin). We conclude that this is challenged by the mass limitation of the surveys and the resulting uncertainty of the mass distribution of the free floaters. The detection of interlopers may be one of the few observational constraints of the small end of this population, with the caveat that, as we conclude here and in Moro-Martín (2018), in the case of 1I/'Oumuamua, it might not be appropriate to assume this object is representative of an isotropic background population, which makes the derivation of a number density very challenging.

We present the apsidal motion and light-curve analyses of 21 eccentric eclipsing binaries located in the Small Magellanic Cloud. Most of these systems have never been studied before, hence their orbital and physical properties as well as the apsidal motion parameters are given here for the first time. All the systems are of early spectral type, having orbital periods up to 4 days. The apsidal motion periods were derived to be from 7.2 to 200 yr (OGLE-SMC-ECL-2194 having the shortest apsidal period among known main-sequence systems). The orbital eccentricities are usually rather mild (median of about 0.06), the maximum eccentricity being 0.33. For the period analysis using O − C diagrams of eclipse timings, in total 951 minima were derived from survey photometry as well as our new data. Moreover, six systems show some additional variation in their O − C diagrams, which should indicate the presence of hidden additional components in them. According to our analysis these third-body variations have periods from 6.9 to 22 yr.

It has been suggested that centaurs may lose their red surfaces and become bluer due to the onset of cometary activity, but the way in which cometary outbursts affect the surface composition and albedo of active centaurs is poorly understood. We obtained consistent visual-near-infrared (VNIR) reflectance spectra of the sporadically active centaur 174P/Echeclus during a period of inactivity in 2014 and six weeks after its outburst in 2016 to see if activity had observably changed the surface properties of the nucleus. We observed no change in the surface reflectance properties of Echeclus following the outburst compared to before, indicating that, in this case, any surface changes due to cometary activity were not sufficiently large to be observable from Earth. Our spectra and post-outburst imaging have revealed, however, that the remaining dust coma is not only blue compared to Echeclus, but also bluer than solar, with a spectral gradient of −7.7 ± 0.6% per 0.1 μm measured through the 0.61–0.88 μm wavelength range that appears to continue up to λ ∼ 1.3 μm before becoming neutral. We conclude that the blue visual color of the dust is likely not a scattering effect, and instead may be indicative of the dust's carbon-rich composition. Deposition of such blue, carbon-rich, comatic dust onto a red active centaur may be a mechanism by which its surface color could be neutralized.

Rotational modulations are observed on brown dwarfs and directly imaged exoplanets, but the underlying mechanism is not well understood. Here we analyze Jupiter's rotational light curves at 12 wavelengths from the ultraviolet (UV) to the mid-infrared (mid-IR). The peak-to-peak amplitudes of Jupiter's light curves range from subpercent to 4% at most wavelengths, but the amplitude exceeds 20% at 5 μm, a wavelength sensing Jupiter's deep troposphere. Jupiter's rotational modulations are primarily caused by discrete patterns in the cloudless belts instead of the cloudy zones. The light-curve amplitude is controlled by the sizes and brightness contrasts of the Great Red Spot (GRS), expansions of the North Equatorial Belt (NEB), patchy clouds in the North Temperate Belt (NTB), and a train of hot spots in the NEB. In reflection, the contrast is controlled by upper tropospheric and stratospheric hazes, clouds, and chromophores in the clouds. In thermal emission, the small rotational variability is caused by the spatial distribution of temperature and opacities of gas and aerosols; the large variation is caused by the NH₃ cloud holes and thin-thick clouds. The methane-band light curves exhibit opposite-shape behavior compared with the UV and visible wavelengths, caused by a wavelength-dependent brightness change of the GRS. Light-curve evolution is induced by periodic events in the belts and longitudinal drifting of the GRS and patchy clouds in the NTB. This study suggests several interesting mechanisms related to distributions of temperature, gas, hazes, and clouds for understanding the observed rotational modulations on brown dwarfs and exoplanets.

Juno is the first polar orbiter around Jupiter. Juno possesses a suite of instruments designed to measure the electron and ion populations in the Jupiter magnetosphere, leading to the powerful Jovian aurorae. The Ultraviolet Spectrograph onboard Juno (Juno-UVS) is a photon-counting imaging spectrograph (68–210 nm), designed to observe and characterize Jupiter's far-ultraviolet aurorae. The instrument borrows heavily from previous Alice and UVS instruments led by Southwest Research Institute (New Horizons and Rosetta Alices, LRO-LAMP), with several major improvements. The pointing flexibility offered by the UVS scan mirror combined with Juno's spin allows UVS access to half of the sky at any given moment. This paper describes how we leverage this extensive database to track the evolution of Juno-UVS calibration with time throughout the mission. UVS observes 7.2° × 360°-long swaths of the sky for each rotation of the spacecraft (nominally 2 rpm). This paper describes how the very substantial amount of stellar spectra has been used to monitor the health of the instrument over the mission. As of PJ14 (2018 July 16), more than 8700 spectra of O, A, and B stars have been extracted in the V-magnitude range of ∼0–7, and more than 99% of the sky was mapped. Selected stars among this list were used to calibrate the UVS bandpass, using observations from the International Ultraviolet Explorer and the Hubble Space Telescope. The retrieved effective area of the instrument is 0.30 ± 0.03 cm² at 125 nm, 0.15 ± 0.02 cm² at 140 nm, and 0.05 ± 0.01 cm² at 160 nm.