Photometric phase curves provide an important window onto exoplanetary atmospheres and potentially even their surfaces. With similar amplitudes to occultations but far longer baselines, they have a higher sensitivity to planetary photons at the expense of a more challenging data reduction in terms of long-term stability. In this work, we introduce a novel non-parametric algorithm dubbed "phasma" to produce clean, robust exoplanet phase curves and apply it to 115 Neptunian and 50 Terran exoplanets observed by Kepler. We stack the signals to further improve signal-to-noise, and measure an average Neptunian albedo of Ag < 0.23 to 95% confidence, indicating a lack of bright clouds consistent with theoretical models. Our Terran sample provides the first constraint on the ensemble albedo of exoplanets which are most likely solid, constraining Ag < 0.42 to 95% confidence. In agreement with our constraint on the greenhouse effect, our work implies that Kepler's solid planets are unlikely to resemble cloudy Venusian analogs, but rather dark Mercurian rocks.

We report the detection of the full orbital phase curve and occultation of the hot-Jupiter WASP-100b using TESS photometry. The phase curve is isolated by suppressing low frequency stellar and instrumental modes using both a non-parametric harmonic notch filter (phasma) and semi-sector long polynomials. This yields a phase curve signal of (73 +/- 9) ppm amplitude, preferred over a null-model by deltaBIC = 25, indicating very strong evidence for an observed effect. We recover the occultation event with a suite of five temporally localized tools, including Gaussian processes and cosine filtering. This allows us to infer an occultation depth of (100 +/- 14) ppm, with an additional +/- 16 ppm systematic error from the differences between methods. We regress a model including atmospheric reflection, emission, ellipsoidal variations and Doppler beaming to the combined phase curve and occultation data. This allows us to infer that WASP-100b has a geometric albedo of A_g = 0.16 (+0.04, -0.03) in the TESS bandpass, with a maximum dayside brightness temperature of (2710 +/- 100) K and a warm nightside temperature of (2380 [+170, -200]) K. Additionally, we find evidence that WASP-100b has a high thermal redistribution efficiency, manifesting as a substantial eastward hotspot offset of (71 [+2, -4]) degrees. These results present the first measurement of a thermal phase shift among the phase curves observed by TESS so far, and challenge the predicted efficiency of heat transport in the atmospheres of ultra-hot Jupiters.

Planetary rotation rate has a significant effect on atmospheric circulation, where the strength of the Coriolis effect in part determines the efficiency of latitudinal heat transport, altering cloud distributions, surface temperatures, and precipitation patterns. In this study we use the ROCKE-3D dynamic-ocean general circulation model to study the effects of slow rotations and increased insolations on the 'fractional habitability' and silicate weathering rate of an Earth-like world. Defining the fractional habitability f_h to be the percentage of a planet's surface which falls in the 0 <= T <= 100 C temperature regime, we find a moderate increase in f_h with a 10% and 20% increase in insolation and a possible maximum in f_h at sidereal day-lengths between 8 and 32 times that of the modern Earth. By tracking precipitation and run-off we further determine that there is a rotational regime centered on a 4-day period in which the silicate weathering rate is maximized and is particularly strongly peaked at higher overall insolations. Because of weathering's integral role in the long-term carbonate-silicate cycle, we suggest that climate stability may be strongly affected by the anticipated rotational evolution of temperate terrestrial-type worlds, and should be considered a major factor in their study. In light of our results we argue that planetary rotation period is an important factor to consider when determining the habitability of terrestrial worlds.

We analyze K2 light curves for 132 low-mass ($1\ \gtrsim\ M_*\ \gtrsim\ 0.1$~${M_{\odot}}$) members of the 600--800~Myr-old Hyades cluster and measure rotation periods ($P_{rot}$) for 116 of these stars. These include 93 stars with no prior $P_{rot}$ measurement; the total number of Hyads with known $P_{rot}$ is now 232. We then combine literature binary data with Gaia DR2 photometry and astrometry to select single star sequences in the Hyades and its roughly coeval Praesepe open cluster, and derive a new reddening value of $A_V = 0.035$$\pm$$0.011$ for Praesepe. Comparing the effective temperature--$P_{rot}$ distributions for the Hyades and Praesepe, we find that solar-type Hyads rotate, on average, 0.4~d slower than their Praesepe counterparts. This $P_{rot}$ difference indicates that the Hyades is slightly older than Praesepe: we apply a new gyrochronology model tuned with Praesepe and the Sun, and find an age difference between the two clusters of 57~Myr. However, this $P_{rot}$ difference decreases and eventually disappears for lower-mass stars. This provides further evidence for stalling in the rotational evolution of these stars, and highlights the need for more detailed analysis of angular-momentum evolution for stars of different masses and ages.

We present the discovery of a pair of transiting giant planets, TOI-216b and c, using four sectors of TESS photometry. TOI-216 is a $0.87 M_{\odot}$ dwarf orbited by two transiters with radii of $8.2 R_{\oplus}$ and $11.3 R_{\oplus}$, and periods of $17.01$d and $34.57$d, respectively. Anti-correlated TTVs are clearly evident indicating that the transiters orbit the same star and interact via a near 2:1 mean motion resonance. By fitting the TTVs with a dynamical model, we infer masses of $26_{-11}^{+24} M_{\oplus}$ and $190_{-80}^{+220} M_{\oplus}$, establishing that the objects are planetary in nature and have likely sub-Kronian and Kronian densities. TOI-216 lies close to the southern ecliptic pole and thus will be observed by TESS throughout the first year, providing an opportunity for continuous dynamical monitoring and considerable refinement of the dynamical masses presented here. TOI-216 closely resembles Kepler-9 in architecture, and we hypothesize that in such systems these Saturn-analogs failed to fully open a gap and thus migrated far deeper into the system before becoming trapped into resonance, which would imply that future detections of new analogs may also have sub-Jupiter masses.

An exoplanet's optical phase curve constrains the thermal emission and albedo of the planet's surface and/or atmosphere, as well as potentially constraining the mass via gravitational influences on the host star. Recently, Jansen & Kipping (2018) demonstrated that exoplanets with precisely constrained orbital periods, as is typical of transiting planets, enable one to exploit a Fourier-argument to non-parametrically separate optical phase curves from long-term photometric time series in the presence of realistic noise structures. This algorithm, dubbed phasma, was applied to 477 confirmed exoplanets in the original paper. In this research note, we extend the sample to all non-false positive Kepler Objects of Interest with periods less than ten days. Our sample of 1998 detrended phase curves are made publicly available to help the community search for various effects. A Python port of the phasma code is also made available.

We present two state-of-the-art models of the solar system, one corresponding to the present day and one to the Archean Eon 3.5 billion years ago. Each model contains spatial and spectral information for the star, the planets, and the interplanetary dust, extending to 50 AU from the sun and covering the wavelength range 0.3 to 2.5 micron. In addition, we created a spectral image cube representative of the astronomical backgrounds that will be seen behind deep observations of extrasolar planetary systems, including galaxies and Milky Way stars. These models are intended as inputs to high-fidelity simulations of direct observations of exoplanetary systems using telescopes equipped with high-contrast capability. They will help improve the realism of observation and instrument parameters that are required inputs to statistical observatory yield calculations, as well as guide development of post-processing algorithms for telescopes capable of directly imaging Earth-like planets.

Direct imaging of extrasolar planets with future space-based coronagraphic telescopes may provide a means of detecting companion moons at wavelengths where the moon outshines the planet. We propose a detection strategy based on the positional variation of the center of light with wavelength, "spectroastrometry." This new application of this technique could be used to detect an exomoon, to determine the exomoon's orbit and the mass of the host exoplanet, and to disentangle of the spectra of the planet and moon. We consider two model systems, for which we discuss the requirements for detection of exomoons around nearby stars. We simulate the characterization of an Earth-Moon analog system with spectroastrometry, showing that the orbit, the planet mass, and the spectra of both bodies can be recovered. To enable the detection and characterization of exomoons we recommend that coronagraphic telescopes should extend in wavelength coverage to 3 micron, and should be designed with spectroastrometric requirements in mind.

We describe the design, manufacture, and performance of bare-fiber integral field units (IFUs) for the SDSS-IV survey MaNGA (Mapping Nearby Galaxies at APO) on the the Sloan 2.5 m telescope at Apache Point Observatory (APO). MaNGA is a luminosity-selected integral-field spectroscopic survey of 10,000 local galaxies covering 360-1030 nm at R ~ 2200. The IFUs have hexagonal dense packing of fibers with packing regularity of 3 um (RMS), and throughput of 96+/-0.5% from 350 nm to 1 um in the lab. Their sizes range from 19 to 127 fibers (3-7 hexagonal layers) using Polymicro FBP 120:132:150 um core:clad:buffer fibers to reach a fill fraction of 56%. High throughput (and low focal-ratio degradation) is achieved by maintaining the fiber cladding and buffer intact, ensuring excellent surface polish, and applying a multi-layer AR coating of the input and output surfaces. In operations on-sky, the IFUs show only an additional 2.3% FRD-related variability in throughput despite repeated mechanical stressing during plate plugging (however other losses are present). The IFUs achieve on-sky throughput 5% above the single-fiber feeds used in SDSS-III/BOSS, attributable to equivalent performance compared to single fibers and additional gains from the AR coating. The manufacturing process is geared toward mass-production of high-multiplex systems. The low-stress process involves a precision ferrule with hexagonal inner shape designed to lead inserted fibers to settle in a dense hexagonal pattern. The ferrule inner diameter is tapered at progressively shallower angles toward its tip and the final 2 mm are straight and only a few um larger than necessary to hold the desired number of fibers. This process scales to accommodate other fiber sizes and to IFUs with substantially larger fiber count. (Abridged)

We present the largest sample of flares ever compiled for a single M dwarf, the active M4 star GJ 1243. Over 6100 individual flare events, with energies ranging from $10^{29}$ to $10^{33}$ erg, are found in 11 months of 1-minute cadence data from Kepler. This sample is unique for its completeness and dynamic range. We have developed automated tools for finding flares in short-cadence Kepler light curves, and performed extensive validation and classification of the sample by eye. From this pristine sample of flares we generate a median flare template. This template shows that two exponential cooling phases are present during the white-light flare decay, providing fundamental constraints for models of flare physics. The template is also used as a basis function to decompose complex multi-peaked flares, allowing us to study the energy distribution of these events. Only a small number of flare events are not well fit by our template. We find that complex, multi-peaked flares occur in over 80% of flares with a duration of 50 minutes or greater. The underlying distribution of flare durations for events 10 minutes and longer appears to follow a broken power law. Our results support the idea that sympathetic flaring may be responsible for some complex flare events.