We present measurements of the near-infrared brightness of Io's hot spots derived from 2-5 micron imaging with adaptive optics on the Keck and Gemini N telescopes. The data were obtained on 271 nights between August 2013 and the end of 2018, and include nearly 1000 detections of over 75 unique hot spots. The 100 observations obtained between 2013 and 2015 have been previously published in de Kleer and de Pater (2016a); the observations since the start of 2016 are presented here for the first time, and the analysis is updated to include the full five-year dataset. These data provide insight into the global properties of Io's volcanism. Several new hot spots and bright eruptions have been detected, and the preference for bright eruptions to occur on Io's trailing hemisphere noted in the 2013-2015 data (de Kleer and de Pater 2016a) is strengthened by the larger dataset and remains unexplained. The program overlapped in time with Sprint-A/EXCEED and Juno observations of the jovian system, and correlations with transient phenomena seen in other components of the system have the potential to inform our understanding of the impact of Io's volcanism on Jupiter and its neutral/plasma environment.

The Euclid telescope, due for launch in 2021, will perform an imaging and slitless spectroscopy survey over half the sky, to map baryon wiggles and weak lensing. During the survey Euclid is expected to resolve 100,000 strong gravitational lens systems. This is ideal to find rare lens configurations, provided they can be identified reliably and on a reasonable timescale. For this reason we have developed a Convolutional Neural Network (CNN) that can be used to identify images containing lensing systems. CNNs have already been used for image and digit classification as well as being used in astronomy for star-galaxy classification. Here our CNN is trained and tested on Euclid-like and KiDS-like simulations from the Euclid Strong Lensing Group, successfully classifying 77% of lenses, with an area under the ROC curve of up to 0.96. Our CNN also attempts to classify the lenses in COSMOS HST F814W-band images. After convolution to the Euclid resolution, we find we can recover most systems that are identifiable by eye. The Python code is available on Github.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The masses, atmospheric makeups, spin-orbit alignments, and system architectures of extrasolar planets can be best studied when the planets orbit bright stars. We report the discovery of three bodies orbiting HD 106315, a bright (V = 8.97 mag) F5 dwarf targeted by our K2 survey for transiting exoplanets. Two small, transiting planets have radii of 2.23 (+0.30/-0.25) R_Earth and 3.95 (+0.42/-0.39) R_Earth and orbital periods of 9.55 d and 21.06 d, respectively. A radial velocity (RV) trend of 0.3 +/- 0.1 m/s/d indicates the likely presence of a third body orbiting HD 106315 with period >160 d and mass >45 M_Earth. Transits of this object would have depths of >0.1% and are definitively ruled out. Though the star has v sin i = 13.2 km/s, it exhibits short-timescale RV variability of just 6.4 m/s, and so is a good target for RV measurements of the mass and density of the inner two planets and the outer object's orbit and mass. Furthermore, the combination of RV noise and moderate v sin i makes HD 106315 a valuable laboratory for studying the spin-orbit alignment of small planets through the Rossiter-McLaughlin effect. Space-based atmospheric characterization of the two transiting planets via transit and eclipse spectroscopy should also be feasible. This discovery demonstrates again the power of K2 to find compelling exoplanets worthy of future study.

We are conducting an experiment to search for WIMPs, or weakly-interacting massive particles, in the galactic halo using terrestrial detectors. This generic class of hypothetical particles, whose properties are similar to those predicted by extensions of the standard model of particle physics, could comprise the cold component of non-baryonic dark matter. We describe our experiment, which is based on cooled germanium and silicon detectors in a shielded low-background cryostat. The detectors achieve a high degree of background rejection through the simultaneous measurement of the energy in phonons and ionization. Using exposures on the order of one kilogram-day from initial runs of our experiment, we have achieved (preliminary) upper limits on the WIMP-nucleon cross section that are comparable to much longer runs of other experiments.

We propose and investigate a new algorithm for quantifying the topological properties of cosmological density fluctuations. We first motivate this algorithm by drawing a formal distinction between two definitions of relevant topological characteristics, based on concepts, on the one hand, from differential topology and, on the other, from integral geometry. The former approach leads one to concentrate on properties of the contour surfaces which, in turn, leads to the algorithms CONTOUR2D and CONTOUR3D familiar to cosmologists. The other approach, which we adopt here, actually leads to much simpler algorithms in both two and three dimensions. (The 2D algorithm has already been introduced to the astronomical literature.) We discuss the 3D case in some detail and compare results obtained with it to analogous results using the CONTOUR3D algorithm.