The exploration of the binary Pluto-Charon and its small satellites during the New Horizons flyby in 2015 revealed not only widespread geologic and compositional diversity across Pluto, but surprising complexity, a wide range of surface unit ages, evidence for widespread activity stretching across billion of years to the near-present, as well as numerous atmospheric puzzles, and strong atmospheric coupling with its surface. New Horizons also found an unexpected diversity of landforms on its binary companion, Charon. Pluto's four small satellites yielded surprises as well, including their unexpected rapid and high obliquity rotation states, high albedos, and diverse densities. Here we briefly review the findings made by New Horizons and the case for a follow up mission to investigate the Pluto system in more detail. As the next step in the exploration of this spectacular planet-satellite system, we recommend an orbiter to study it in considerably more detail, with new types of instrumentation, and to observe its changes with time. We further call for the in-depth study of Pluto orbiter missions as a precursor to the 2023 Planetary Science Decadal Survey.

Tidal heating is the major source of heat in the outer Solar System. Because of its strong tidal interaction with Jupiter and the other Galilean Satellites, Io is incredibly volcanically active. We use the directly measured volcanic activity level of Io's volcanoes as a proxy for surface heat flow and compare to tidal heating model predictions. Volcanic activity is a better proxy for heat flow than simply the locations of volcanic constructs. We determine the volcanic activity level using three data sets: Galileo PPR, Galileo NIMS, and New Horizons LEISA. We also present a systematic reanalysis of the Galileo NIMS observations to determine the 3.5 micron brightness of 51 active volcanoes. We find that potential differences in volcanic style between high and low latitudes make high latitude observations unreliable in distinguishing between tidal heating models. Observations of Io's polar areas, such as those by JUNO, are necessary to unambiguously understand Io's heat flow. However, all three of the data sets examined show a relative dearth of volcanic brightness near 180 W (anti-Jovian point) and the equator and the only data set with good observations of the sub-Jovian point (LEISA) also shows a lack of volcanic brightness in that region. These observations are more consistent with the mantle-heating model than the asthenospheric-heating model. Furthermore, all three of the data sets are consistent with four-fold symmetry in longitude and peak heat flow at mid-latitudes, which best matches with the combined heating case of Tackley et al. (2001).

We investigate the relationship between the turnaround radius ($R_t$) and the virial mass ($M_v$) of cosmic structures in the context of $\Lambda$CDM model and in an $f(R)$ model of modified gravity -- namely, the Hu-Sawicki model. The $R_t$ is the distance from the center of the cosmic structure to the shell that is detaching from the Hubble flow at a given time, while the $M_v$ is defined, for this work, as the mass enclosed within the volume where the density is $200$ times the background density. We consider that gravitationally bound astrophysical systems follow a Navarro-Frenk-White density profile, while beyond the virial radius ($R_v$) the profile is approximated by the 2-halo term of the matter correlation function. By combining them together with the information drawn from solving the spherical collapse for the structures, we are able to connect two observables: the $R_t$ and the $M_v$. We show that, in $\Lambda$CDM, the turnaround mass ($M_t$) at $z=0$ is related to the $M_v$ of that same structure by $M_t \simeq 3.07 \, M_v$, while in terms of the radii we have that $R_t \simeq 3.7 \, R_v$ (for $M_v$ of $10^{13} \, h^{-1} \, M_\odot$). In the $f(R)$ model, we have $M_t \simeq 3.43 \, M_v$ and $R_t \simeq 4.1 \, R_v$, for $|f_{R0}|=10^{-6}$ and the same mass scale. Therefore, the difference between $\Lambda$CDM and $f(R)$ in terms of these observable relations is of order $\sim 10-20\%$ even for a relatively mild strength of the modification of gravity ($|f_{R0}|=10^{-6}$). For the $R_t$ itself we find a difference of $\sim 9\%$ between the weakly modification in gravity considered in this work ($|f_{R0}|=10^{-6}$) and $\Lambda$CDM for a mass of $10^{13} \, h^{-1} \, M_\odot$. Once observations allow precisions of this order or better in measurements $R_t$, as well as the $M_v$, these quantities will become powerful tests of modified gravity.

We investigate the turnaround radius in the spherical collapse model, both in General Relativity and in modified gravity, in particular $f(R)$ scenarios. The phases of spherical collapse are marked by the density contrast in the instant of turnaround $\delta_t$, and by the linear density contrast in the moment of collapse, $\delta_c$. We find that the effective mass of the extra scalar degree of freedom which arises in modified gravity models has an impact on $\delta_t$ of up to $\sim10\%$, and that $\delta_c$ can increase by $\sim1.0\%$. We also compute the turnaround radius, $R_t$, which in modified gravity models can increase by up to $\sim 6\%$ at $z \simeq 0$.

We present a geomorphologic map of Titan's polar terrains. The map was generated from a combination of Cassini Synthetic Aperture Radar (SAR) and Imaging Science Subsystem imaging products, as well as altimetry, SARTopo and radargrammetry topographic datasets. In combining imagery with topographic data, our geomorphologic map reveals a stratigraphic sequence from which we infer process interactions between units. In mapping both polar regions with the same geomorphologic units, we conclude that processes that formed the terrains of the north polar region also acted to form the landscape we observe at the south. Uniform, SAR-dark plains are interpreted as sedimentary deposits, and are bounded by moderately dissected uplands. These plains contain the highest density of filled and empty lake depressions, and canyons. These units unconformably overlay a basement rock that outcrops as mountains and SAR-bright dissected terrains at various elevations across both poles. All these units are then superposed by surficial units that slope towards the seas, suggestive of subsequent overland transport of sediment. From estimates of the depths of the embedded empty depressions and canyons that drain into the seas, the SAR-dark plains must be >600 m thick in places, though the thickness may vary across the poles. At the lowest elevations of each polar region, there are large seas, which are currently liquid methane/ethane filled at the north and empty at the south. The large plains deposits and the surrounding hillslopes may represent remnant landforms that are a result of previously vast polar oceans, where larger liquid bodies may have allowed for a sustained accumulation of soluble and insoluble sediments, potentially forming layered sedimentary deposits. Coupled with vertical crustal movements, the resulting layers would be of varying solubilities and erosional resistances.

We develop the theory of radiative association of an atom and a diatomic molecule within a close-coupling framework. We apply it to the formation of H$_3^-$ after the low energy collision (below 0.5 eV) of H$_2$ with H$^-$. Using recently obtained potential energy and permanent dipole moment surfaces of H$_3^-$, we calculate the lowest rovibrational levels of the H$_3^-$ electronic ground state, and the cross section for the formation of H$_3^-$ by radiative association between H$^-$ and ortho- and para-H$_2$. We discuss the possibility for the H$_3^-$ ion to be formed and observed in the cold and dense interstellar medium in an environment with a high ionization rate. Such an observation would be a probe for the presence of H$^-$ in the interstellar medium.

We discuss two-dimensional mapping of the near-infrared emission-line intensity distributions and kinematics of the narrow-line region (NLR) of NGC 4151, obtained with the Gemini Near-Infrared Integral Field Spectrograph, with a projected spatial resolution of ~8pc. The ionized gas intensity distribution follows the projected bi-cone morphology observed in previous optical narrow-band images, and its kinematics reveal outflows along the bi-cone. We propose a kinematic model in which the gas in the NLR moves at a velocity of ~600 km s^-1 up to ~100 pc from the nucleus. A completely distinct morphology and kinematic structure is observed for the molecular gas, which avoids the region of the bi-cone and has velocities close to systemic, and is consistent with an origin in the galaxy plane. The molecular gas thus traces the AGN feeding, while the ionized gas traces its feedback.

We present a detailed investigation of the incidence of circumnuclear dust structure in a large, well-matched sample of early-type galaxies with and without Active Galactic Nuclei (AGN). All 34 early-type AGN hosts in our sample have circumnuclear dust, while dust is only observed in 26% (nine) of a pair-matched sample of 34 early-type, inactive galaxies. This result demonstrates a strong correlation between the presence of circumnuclear dust and accretion onto the central, supermassive black hole in elliptical and lenticular galaxies. This correlation is not present at later Hubble types, where a sample of 31 active and 31 inactive galaxies all contain circumnuclear dust. These archival, HST observations reveal a wide range of mostly chaotic dust morphologies. Current estimates suggest the dust settling or destruction time is on order of 10^8 years and therefore the presence of dust in ~ 50% of early-type galaxies requires frequent replenishment and similarly frequent fueling of their central, supermassive black holes. The observed dust could be internally-produced (via stellar winds) or externally-accreted, although there are observational challenges for both of these scenarios. Our analysis also reveals that approximately a third of the early-type galaxies without circumnuclear dust have nuclear stellar disks. These nuclear stellar disks may provide a preferred kinematic axis to externally-accreted material and this material may in turn form new stars in these disks. The observed incidence of nuclear stellar disks and circumnuclear dust suggests that episodic replenishment of nuclear stellar disks occurs and is approximately concurrent with the fueling of the central AGN.