The development of the Skipper Charge Coupled Devices (Skipper-CCDs) has been a major technological breakthrough for sensing very weak ionizing particles. The sensor allows to reach the ultimate sensitivity of silicon material as a charge signal sensor by unambiguous determination of the charge signal collected by each cell or pixel, even for single electron-hole pair ionization. Extensive use of the technology was limited by the lack of specific equipment to operate the sensor at the ultimate performance. In this work a simple, single-board Skipper-CCD controller is presented, aimed for the operation of the detector in high sensitivity scientific applications. The article describes the main components and functionality of the Low Threshold Acquisition (LTA) together with experimental results when connected to a Skipper-CCD sensor. Measurements show unprecedented deep sub-electron noise of 0.039 e$^-_{rms}$/pix for 5000 pixel measurements.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

The Coherent Neutrino-Nucleus Interaction Experiment (CONNIE) uses low-noise fully depleted charge-coupled devices (CCDs) with the goal of measuring low-energy recoils from coherent elastic scattering (CE$\nu$NS) of reactor antineutrinos with silicon nuclei and testing nonstandard neutrino interactions (NSI). We report here the first results of the detector array deployed in 2016, considering an active mass 47.6 g (8 CCDs), which is operating at a distance of 30 m from the core of the Angra 2 nuclear reactor, with a thermal power of 3.8 GW. A search for neutrino events is performed by comparing data collected with reactor on (2.1 kg-day) and reactor off (1.6 kg-day). The results show no excess in the reactor-on data, reaching the world record sensitivity down to recoil energies of about 1 keV (0.1 keV electron-equivalent). A 95% confidence level limit for new physics is established at an event rate of 40 times the one expected from the standard model at this energy scale. The results presented here provide a new window to low-energy neutrino physics, allowing one to explore for the first time the energies accessible through the low threshold of CCDs. They will lead to new constrains on NSI from the CE$\nu$NS of antineutrinos from nuclear reactors.

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7].

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

We present direct detection constraints on the absorption of hidden-photon dark matter with particle masses in the range 1.2-30 eV$c^{-2}$ with the DAMIC experiment at SNOLAB. Under the assumption that the local dark matter is entirely constituted of hidden photons, the sensitivity to the kinetic mixing parameter $\kappa$ is competitive with constraints from solar emission, reaching a minimum value of 2.2$\times$$10^{-14}$ at 17 eV$c^{-2}$. These results are the most stringent direct detection constraints on hidden-photon dark matter in the galactic halo with masses 3-12 eV$c^{-2}$ and the first demonstration of direct experimental sensitivity to ionization signals $<$12 eV from dark matter interactions.

We present results of a dark matter search performed with a 0.6 kg day exposure of the DAMIC experiment at the SNOLAB underground laboratory. We measure the energy spectrum of ionization events in the bulk silicon of charge-coupled devices down to a signal of 60 eV electron equivalent. The data are consistent with radiogenic backgrounds, and constraints on the spin-independent WIMP-nucleon elastic-scattering cross section are accordingly placed. A region of parameter space relevant to the potential signal from the CDMS-II Si experiment is excluded using the same target for the first time. This result obtained with a limited exposure demonstrates the potential to explore the low-mass WIMP region (<10 GeV/$c^{2}$) of the upcoming DAMIC100, a 100 g detector currently being installed in SNOLAB.

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

To date, 17 circumbinary planets have been discovered. In this paper, we focus our attention on the stability of the Kepler circumbinary planetary systems with only one planet, i.e. Kepler-16, Kepler-34, Kepler-35, Kepler-38, Kepler-64 and Kepler-413. In addition to their intrinsic interest, the study of such systems is an opportunity to test our understanding of planetary system formation and evolution around binaries. The investigation is done by means of numerical simulations. We perform numerical integrations of the full equations of motion of each system with the aim of checking the stability of the planetary orbit. The investigation of the stability of the above systems consists of three numerical experiments. In the first one we perform a long term (1Gyr) numerical integration of the nominal solution of the six Kepler systems under investigation. In the second experiment, we look for the critical semimajor axis of the six planetary orbits, and finally, in the third experiment, we construct two dimensional stability maps on the eccentricity-pericentre distance plane. Additionally, using numerical integrations of the nominal solutions we checked if this solutions were close to the exact resonance.

We study the effect of a massive planetesimal disk on the dynamical stability of the outer planets assuming, as has been suggested recently, that these were initially locked in a compact and multiresonant configuration as a result of gas-driven migration in a protoplanetary disk. The gravitational interaction among all bodies in our simulations is included self-consistently using the Mercury6.5 code. Several initial multiresonant configurations and planetesimal disk models are considered. Under such conditions a strong dynamical instability, manifested as a rapid giant planet migration and planetesimal disk dispersal, develops on a timescale of less than 40 Myr in most cases. Dynamical disk heating due to the gravitational interactions among planetesimals leads to more frequent interactions between the planetesimals and the ice giants Uranus and Neptune, in comparison to models in which planetesitmal-planetesimal interactions are neglected. On account of the rapid evolution of the multiresonant configurations obtained with fully self-consistent simulations, our results are inconsistent with the dynamical instability origin of the Late Heavy Bombardment as currently considered by the Nice model for the Solar System.

It has been suggested that the ejection to interplanetary space of terrestrial crustal material, accelerated in a large impact, may result in the interchange of biological material between Earth and other Solar System bodies. In this paper, we analyze the fate of debris ejected from Earth by means of numerical simulations of the dynamics of a large collection of test particles. This allows us to determine the probability and conditions for the collision of ejecta with other planets of the Solar System. We also estimate the amount of particles falling-back to Earth as a function of time after being ejected. We find that, in general, the collision rates of Earth ejecta with Venus and the Moon, as well as the fall-back rates, are consistent with results reported in the literature. By considering a larger number of particles than in all previous calculations we have also determined directly the collision probability with Mars and, for the first time, computed collision probabilities with Jupiter. We find that the collision probability with Mars is greater than values determined from collision cross section estimations previously reported.

We study some of the effects that magnetic turbulent fluctuations have on the dynamics of pickup O+ ions in the magnetic polar regions of the Mars ionosheath. In particular we study their effect on the bulk velocity profiles of ions as a function of altitude over the magnetic poles, in order to compare them with recent Mars Express data; that indicate that their average velocity is very low and essentially in the anti-sunward direction. We find that, while magnetic field fluctuations do give rise to deviations from simple ExB-drift gyromotion, even fluctuation amplitudes much greater than those of in situ measurements are {\it not} able to reproduce the vertical velocity profile of O+ ions. We conclude that other physical mechanisms, different from a pure charged particle dynamics, are acting on pickup ions at the Martian terminator. A possibility being a viscous-like interaction between the Solar Wind and the Martian ionosphere at low altitudes.

We study the dynamics of the upsilon-Andromedae planetary system proposed by Curiel et al. (2011). We focus on the study of the 3:1 Mean Motion Resonance between upsilon-Andromedae-d and the recently discovered upsilon-Andromedae-e (hereafter upsilon-And-d and upsilon-And-e). Numerical simulations of the dynamics of the four planet system are conducted. The previously reported apsidal resonance between upsilon-And-c and upsilon-And-d is confirmed. In addition, we find that upsilon-And-d and upsilon-And-e are also in an apsidal resonance condition. Our results further indicate that the upsilon-Andromedae planetary system configuration is in the middle of a stability island in the semimajor axis-eccentricity domain. Additionally, we performed numerical integrations of the planetary configuration over 500 Myr and found it to be stable. We conclude that, within the uncertainties in the value of the orbital parameters, it is likely that upsilon-Andromedae planetary system will remain stable for a long timescale.

We present here results of numerical calculations demonstrating that the very long periodic variability detected in the Cataclysmic Variable FS Aur can be the result of eccentricity modulation of a close binary (CB) orbit induced by the presence of a third body on a circumbinary orbit. A third component with a substellar mass on a circular, relatively close orbit, modulates the mass transfer rate of the binary on much longer time scales than periods within the triple system. We also report here preliminary results of X-ray observations of FS Aur, providing further evidence that it contains magnetic and freely precessing white dwarf. These two findings allow us to incorporate the new and previously stressed hypothesis on the nature of FS Aur into one consistent model.

In this article we explore the aspect of the F ring with respect to the anti-alignment configuration between the ring and Prometheus. We focus our attention on the shape of the F ring's azimuthal channels which were first reported by Porco et al. (2005) and numerically explored by Murray et al. (2005), who found excellent agreement between Cassini's ISS reprojected images and their numerical model via a direct comparison. We find that for anti-alignment the channels are wider and go deeper inside the ring material. From our numerical model we find a new feature, an island in the middle of the channel. This island is made up of the particles that have been perturbed the most by Prometheus and only appears when this satellite is close to apoapsis. In addition, plots of the anti-alignment configuration for different orbital stages of Prometheus are obtained and discussed here.