Decadal-long radial velocity surveys have recently started to discover analogs to the most influential planet of our solar system, Jupiter. Detecting and characterizing these worlds is expected to shape our understanding of our uniqueness in the cosmos. Despite the great successes of recent transit surveys, Jupiter analogs represent a terra incognita, owing to the strong intrinsic bias of this method against long orbital periods. We here report on the first validated transiting Jupiter analog, Kepler-167e (KOI-490.02), discovered using Kepler archival photometry orbiting the K4-dwarf KIC-3239945. With a radius of (0.91 ± 0.02) R J , a low orbital eccentricity (0.06 +0.10 −0.04 ) and an equilibrium temperature of (131 ± 3) K, Kepler-167e bears many of the basic hallmarks of Jupiter. Kepler-167e is accompanied by three Super-Earths on compact orbits, which we also validate, leaving a large cavity of transiting worlds around the habitable-zone. With two transits and continuous photometric coverage, we are able to uniquely and precisely measure the orbital period of this post snow-line planet (1071.2323 ± 0.0006 d), paving the way for follow-up of this K = 11.8 mag target.

The analysis of Proxima Centauri's radial velocities recently led Anglada-Escudé et al. (2016) to claim the presence of a low mass planet orbiting the Sun's nearest star once every 11.2 days. Although the a-priori probability that Proxima b transits its parent star is just 1.5%, the potential impact of such a discovery would be considerable. Independent of recent radial velocity efforts, we observed Proxima Centauri for 12.5 days in 2014 and 31 days in 2015 with the MOST space telescope. We report here that we cannot make a compelling case that Proxima b transits in our precise photometric time series. Imposing an informative prior on the period and phase, we do detect a candidate signal with the expected depth. However, perturbing the phase prior across 100 evenly spaced intervals reveals one strong false-positive and one weaker instance. We estimate a false-positive rate of at least a few percent and a much higher false-negative rate of 20-40%, likely caused by the very high flare rate of Proxima Centauri. Comparing our candidate signal to HATSouth ground-based photometry reveals that the signal is somewhat, but not conclusively, disfavored (1-2 σ) leading us to argue that the signal is most likely spurious. We expect that infrared photometric follow-up could more conclusively test the existence of this candidate signal, owing to the suppression of flare activity and the impressive infrared brightness of the parent star.

Observational biases distort our view of nature, such that the patterns we see within a surveyed population of interest are often unrepresentative of the truth we seek. Transiting planets currently represent the most informative data set on the ensemble properties of exoplanets within 1 au of their star. However, the transit method is inherently biased due to both geometric and detection-driven effects. In this work, we derive the overall observational biases affecting the most basic transit parameters from first principles. By assuming a trapezoidal transit and using conditional probability, we infer the expected distribution of these terms both as a joint distribution and in a marginalized form. These general analytic results provide a baseline against which to compare trends predicted by mission-tailored injection/recovery simulations and offer a simple way to correct for observational bias. Our results explain why the observed population of transiting planets displays a non-uniform impact parameter distribution, with a bias towards near-equatorial geometries. We also find that the geometric bias towards observed planets transiting near periastron is attenuated by the longer durations which occur near apoastron. Finally, we predict that the observational bias with respect to ratio-of-radii is super-quadratic, scaling as (R P /R ) 5/2 , driven by an enhanced geometric transit probability and modestly longer durations.

The true multiplicity distribution of transiting planet systems is obscured by strong observational biases, leading low-multiplicity systems to be overrepresented in the observed sample. Using the Kepler FGK planet hosts, we employ approximate Bayesian computation to infer the multiplicity distribution by comparing simulated catalogs to the observed one. After comparing a total of ten different multiplicity distributions, half of which were two-population models, to the observed data, we find that a singlepopulation model following a Zipfian distribution is able to explain the Kepler data as well as any of the dichotomous models we test. Our work provides another example of a way to explain the observed Kepler multiplicities without invoking a dichotomous planet population. Using our preferred Zipfian model, we estimate that an additional 2393 +904 −717 planets likely reside in the 1537 FGK Kepler systems studied in this work, which would increase the planet count by a factor of 2.22 +0.46 −0.36 . Of these hidden worlds, 663 +158 −151 are expected to reside in ostensibly single-transiting-planet systems, meaning that an additional planet(s) is expected for approximately 1-in-2 such Kepler systems.

Future liquid-argon DarkSide-20k and Argo detectors, designed for direct dark matter search, will be sensitive also to core-collapse supernova neutrinos, via coherent elastic neutrino-nucleus scattering. This interaction channel is flavor-insensitive with a high-cross section, enabling for a high-statistics neutrino detection with target masses of ∼50 t and ∼360 t for DarkSide-20k and Argo respectively. Thanks to the low-energy threshold of ∼0.5 keV nr achievable by exploiting the ionization channel, DarkSide-20k and Argo have the potential to discover supernova bursts throughout our galaxy and up to the Small Magellanic Cloud, respectively, assuming a 11-M progenitor star. We report also on the sensitivity to the neutronization burst, whose electron neutrino flux is suppressed by oscillations when detected via charged current and elastic scattering. Finally, the accuracies in the reconstruction of the average and total neutrino energy in the different phases of the supernova burst, as well as its time profile, are also discussed, taking into account the expected background and the detector response.