We present a point-source detection algorithm that employs the second order Spherical Mexican Hat wavelet filter (SMHW2), and use it on C-BASS northern intensity data to produce a catalogue of point-sources. This catalogue allows us to cross-check the C-BASS flux-density scale against existing source surveys, and provides the basis for a source mask which will be used in subsequent C-BASS and cosmic microwave background (CMB) analyses. The SMHW2 allows us to filter the entire sky at once, avoiding complications from edge effects arising when filtering small sky patches. The algorithm is validated against a set of Monte Carlo simulations, consisting of diffuse emission, instrumental noise, and various point-source populations. The simulated source populations are successfully recovered. The SMHW2 detection algorithm is used to produce a $4.76\,\mathrm{GHz}$ northern sky source catalogue in total intensity, containing 1784 sources and covering declinations $\delta\geq-10^{\circ}$. The C-BASS catalogue is matched with the Green Bank 6\,cm (GB6) and Parkes-MIT-NRAO (PMN) catalogues over their areas of common sky coverage. From this we estimate the $90$ per cent completeness level to be approximately $610\,\mathrm{mJy}$, with a corresponding reliability of $98$ per cent, when masking the brightest $30$ per cent of the diffuse emission in the C-BASS northern sky map. We find the C-BASS and GB6 flux-density scales to be consistent with one another to within approximately $4$ per cent.

We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $\nu_1$ and $\nu_3$ vibrational bands at 2.7 $\mu$m, from the $\nu_2$ vibrational band at 6.1 $\mu$m, and from pure rotational transitions near 10-13 $\mu$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission.

Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range $10$-$60\,$GHz and a new window into the properties of sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the $\approx10^{\circ}$ diameter $\lambda$ Orionis ring by combining intensity data from the QUIJOTE experiment at $11$, $13$, $17$ and $19\,$GHz and the C-Band All Sky Survey (C-BASS) at $4.76\,$GHz, together with 19 ancillary datasets between $1.42$ and $3000\,$GHz. Maps of physical parameters at $1^{\circ}$ resolution are produced through Markov Chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a log-normal distribution. AME is detected in excess of $20\,\sigma$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx35\,$GHz near the free-free region to $\approx21\,$GHz in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalised by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role.

The cosmic microwave background $B$-mode signal is potentially weaker than the diffuse Galactic foregrounds over most of the sky at any frequency. A common method of separating the CMB from these foregrounds is via pixel-based parametric-model fitting. There are not currently enough all-sky maps to fit anything more than the most simple models of the sky. By simulating the emission in seven representative pixels, we demonstrate that the inclusion of a 5 GHz data point allows for more complex models of low-frequency foregrounds to be fitted than at present. It is shown that the inclusion of the CBASS data will significantly reduce the uncertainties in a number of key parameters in the modelling of both the galactic foregrounds and the CMB. The extra data allow estimates of the synchrotron spectral index to be constrained much more strongly than is presently possible, with corresponding improvements in the accuracy of the recovery of the CMB amplitude. However, we show that to place good limits on models of the synchrotron spectral curvature will require additional low-frequency data.

The chemical composition of high-mass protostars reflects the physical evolution associated with different stages of star formation. In this study, the molecular inventory of the forming high-mass star AFGL 4176 is studied in detail at high angular resolution (~0.35 arcsec) using ALMA. This high resolution makes it possible to separate the emission associated with the inner hot envelope and disk around the forming star from that of its cool outer envelope. In addition, the high sensitivity of ALMA makes it possible to identify weak and optically thin lines and allows for many isotopologues to be detected, providing a more complete and accurate inventory of the source. A total of 23 different molecular species and their isotopologues are detected in the spectrum towards AFGL 4176. The most abundant species is methanol (CH3OH), remaining species are present at levels between 0.003 % and 15 % with respect to CH3OH. Hints that N-bearing species peak slightly closer to the location of the peak continuum emission than the O-bearing species are seen. AFGL 4176 comprises a rich chemical inventory including many complex species present on disk-scales. On average, the derived column density ratios with respect to methanol of O-bearing species are higher than those derived for N-bearing species by a factor of three. This may indicate that AFGL 4176 is a relatively young source since nitrogen chemistry generally takes longer to evolve in the gas-phase. Taking methanol as a reference, the composition of AFGL 4176 more closely resembles that of the low-mass protostar IRAS 16293-2422B than that of high-mass star-forming regions located near the Galactic centre. This similarity hints that the chemical composition of complex species is already set in the cold cloud stage and implies that AFGL 4176 is a young source whose chemical composition has not yet been strongly processed by the central protostar.

The TRAPPIST-1 planetary system consists of seven planets within 0.05 au of each other, five of which are in a multi-resonant chain. {These resonances suggest the system formed via planet migration; subsequent tidal evolution has damped away most of the initial eccentricities. We used dynamical N-body simulations to estimate how long it takes for the multi-resonant configuration that arises during planet formation to break. From there we use secular theory to pose limits on the tidal parameters of planets b and c. We calibrate our results against multi-layered interior models constructed to fit the masses and radii of the planets, from which the tidal parameters are computed independently.} The dynamical simulations show that the planets typically go unstable 30 Myr after their formation. {Assuming synchronous rotation throughout} we compute $\frac{k_2}{Q} \gtrsim 2\times 10^{-4}$ for planet b and $\frac{k_2}{Q} \gtrsim 10^{-3}$ for planet c. Interior models yield $(0.075-0.37) \times 10^{-4}$ for TRAPPIST-1 b and $(0.4-2)\times 10^{-4}$ for TRAPPIST-1 c. The agreement between the {dynamical and interior} models is not too strong, but is still useful to constrain the dynamical history of the system. We suggest that this two-pronged approach could be of further use in other multi-resonant systems if the planet's orbital and interior parameters are sufficiently well known.

The satellite of (225088) 2007 OR10 was discovered on archival Hubble Space Telescope images and along with new observations with the WFC3 camera in late 2017 we have been able to determine the orbit. The orbit's notable eccentricity, e$\approx$0.3, may be a consequence of an intrinsically eccentric orbit and slow tidal evolution, but may also be caused by the Kozai mechanism. Dynamical considerations also suggest that the moon is small, D$_{eff}$ $<$ 100 km. Based on the newly determined system mass of 1.75x10$^{21}$ kg, 2007 OR10 is the fifth most massive dwarf planet after Eris, Pluto, Haumea and Makemake. The newly determined orbit has also been considered as an additional option in our radiometric analysis, provided that the moon orbits in the equatorial plane of the primary. Assuming a spherical shape for the primary this approach provides a size of 1230$\pm$50 km, with a slight dependence on the satellite orbit orientation and primary rotation rate chosen, and a bulk density of 1.75$\pm$0.07 g cm$^{-3}$ for the primary. A previous size estimate that assumed an equator-on configuration (1535$^{+75}_{-225}$ km) would provide a density of 0.92$^{+0.46}_{-0.14}$ g cm$^{-3}$, unexpectedly low for a 1000 km-sized dwarf planet.

The C-Band All-Sky Survey C-BASS is a high-sensitivity all-sky radio survey at an angular resolution of 45 arcmin and a frequency of 4.7 GHz. We present a total intensity 4.7 GHz map of the North Celestial Pole (NCP) region of sky, above declination +80 deg, which is limited by source confusion at a level of ~0.6 mK rms. We apply the template-fitting (cross-correlation) technique to WMAP and Planck data, using the C-BASS map as the synchrotron template, to investigate the contribution of diffuse foreground emission at frequencies ~20-40 GHz. We quantify the anomalous microwave emission (AME) that is correlated with far-infrared dust emission. The AME amplitude does not change significantly (<10%) when using the higher frequency C-BASS 4.7 GHz template instead of the traditional Haslam 408 MHz map as a tracer of synchrotron radiation. We measure template coefficients of $9.93\pm0.35$ and $9.52\pm0.34$ K per unit $\tau_{353}$ when using the Haslam and C-BASS synchrotron templates, respectively. The AME contributes $55\pm2\,\mu$K rms at 22.8 GHz and accounts for ~60% of the total foreground emission. Our results suggest that a harder (flatter spectrum) component of synchrotron emission is not dominant at frequencies >5 GHz; the best-fitting synchrotron temperature spectral index is $\beta=-2.91\pm0.04$ from 4.7 to 22.8 GHz and $\beta=-2.85\pm0.14$ from 22.8 to 44.1 GHz. Free-free emission is weak, contributing ~$7\,\mu$K rms (~7%) at 22.8 GHz. The best explanation for the AME is still electric dipole emission from small spinning dust grains.

Context. New estimates of the masses and radii of the seven planets orbiting the ultracool M-dwarf TRAPPIST-1 star permit improved modelling of their compositions, heating by tidal dissipation, and removal of tidal heat by solid-state convection. Aims. Here, we compute the heat flux due to insolation and tidal heating for the inner four planets. Methods. We apply a Maxwell viscoelastic rheology to compute the tidal response of the planets using the volume-weighted average of the viscosities and rigidities of the metal, rock, high-pressure ice and liquid water/ice I layers. Results. We show that TRAPPIST-1d and e can avoid entering a runaway greenhouse state. Planet e is the most likely to support a habitable environment, with Earth-like surface temperatures and possibly liquid water oceans. Planet d also avoids a runaway greenhouse, if its surface reflectance is at least as high as that of the Earth. Planets b and c, closer to the star, have heat fluxes high enough to trigger a runaway greenhouse and support volcanism on the surfaces of their rock layers, rendering them too warm for life. Planets f, g, and h are too far from the star to experience significant tidal heating, and likely have solid ice surfaces with possible subsurface liquid water oceans.

Ammonia, if present in the ice shells of icy satellites, could lower the temperature for the onset of melting to 176 K and create a large temperature range where partial melt is thermally stable. The evolution of regions of ammonia-rich partial melt could strongly influence the geological and thermal evolution of icy bodies. For melt to be extracted from partially molten regions, the surrounding solid matrix must deform and compact. Whether ammonia-rich melts sink to the subsurface ocean or become frozen into the ice shell depends on the compaction rate and thermal evolution. Here we construct a model for the compaction and thermal evolution of a partially molten, ammonia-rich ice shell in a one-dimensional geometry. We model the thickening of an initially thin ice shell above an ocean with $10\%$ ammonia. We find that ammonia-rich melts can freeze into the upper $5$ to $10$ kilometers of the ice shell, when ice shell thickening is rapid compared to the compaction rate. The trapping of near-surface volatiles suggests that, upon reheating of the ice shell, eutectic melting events are possible. However, as the ice shell thickening rate decreases, ammonia-rich melt is efficiently excluded from the ice shell and the bulk of the ice shell is pure water ice. We apply our results to the thermal evolution of Neptune's moon Triton. As Triton's ice shell thickens, the gradual increase of ammonia concentration in Triton's subsurface ocean helps to prevent freezing and increases the predicted final ocean thickness by up to $50$ km.

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $\nu_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO$_2$ gas in the massive YSO MonR2 IRS3. This yields a SO$_2$/$\mathrm{H}$ abundance lower limit of $5.6\pm0.5\times10^{-7}$, or $>\!4\%$ of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of $b<3.20\;\mathrm{km\;s}^{-1}$. The small line widths and high temperature ($T_\mathrm{ex}=234\pm15\;\mathrm{K}$) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov--Gerasimenko), which is followed by SO$_2$ formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO$_2$ from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

We have performed a 5-8 $\mu$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $\mu$m using the iSHELL instrument on the Infrared Telescope Facility (IRTF), and the full N band window from 8-13 $\mu$m using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the IRTF. Here we present the first detection of ro-vibrational transitions of CS in this source. The absorption lines are centred on average around -10 kms$^{-1}$ and the line widths of CS compare well with the hot component of $^{13}$CO (around 10 kms$^{-1}$). Temperatures for CS, hot $^{13}$CO and $^{12}$CO v=1-2 agree well and are around 700 K. We derive a CS abundance of 8$\times$10$^{-3}$ and 2$\times$10$^{-6}$ with respect to CO and H$_2$ respectively. This enhanced CS abundance with respect to the surrounding cloud (1$\times$10$^{-8}$) may reflect sublimation of H$_2$S ice followed by gas-phase reactions to form CS. Transitions are in LTE and we derive a density of $>$10$^7$ cm$^{-3}$, which corresponds to an absorbing region of $<$0.04$''$. EXES observations of CS are likely to probe deeply into the hot core, to the base of the outflow. Submillimeter and infrared observations trace different components of the hot core as revealed by the difference in systemic velocities, line widths and temperatures, as well as the CS abundance.

Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016.

Fatigue can cause materials that undergo cyclic loading to experience brittle failure at much lower stresses than under monotonic loading. We propose that the lithospheres of icy satellites could become fatigued and thus weakened by cyclical tidal stresses. To test this hypothesis, we performed a series of laboratory experiments to measure the fatigue of water ice at temperatures of $198$ K and $233$ K and at a loading frequency of $1$ Hz. We find that ice is \textit{not} susceptible to fatigue at our experimental conditions and that the brittle failure stress does not decrease with increasing number of loading cycles. Even though fatigue was not observed at our experimental conditions, colder temperatures, lower loading frequencies, and impurities in the ice shells of icy satellites may increase the likelihood of fatigue crack growth. We also explore other mechanisms that may explain the weak behavior of the lithospheres of some icy satellites.

With seven planets, the TRAPPIST-1 system has the largest number of exoplanets discovered in a single system so far. The system is of astrobiological interest, because three of its planets orbit in the habitable zone of the ultracool M dwarf. Assuming the planets are composed of non-compressible iron, rock, and H$_2$O, we determine possible interior structures for each planet. To determine how much tidal heat may be dissipated within each planet, we construct a tidal heat generation model using a single uniform viscosity and rigidity for each planet based on the planet's composition. With the exception of TRAPPIST-1c, all seven of the planets have densities low enough to indicate the presence of significant H$_2$O in some form. Planets b and c experience enough heating from planetary tides to maintain magma oceans in their rock mantles; planet c may have eruptions of silicate magma on its surface, which may be detectable with next-generation instrumentation. Tidal heat fluxes on planets d, e, and f are lower, but are still twenty times higher than Earth's mean heat flow. Planets d and e are the most likely to be habitable. Planet d avoids the runaway greenhouse state if its albedo is $\gtrsim$ 0.3. Determining the planet's masses within $\sim0.1$ to 0.5 Earth masses would confirm or rule out the presence of H$_2$O and/or iron in each planet, and permit detailed models of heat production and transport in each planet. Understanding the geodynamics of ice-rich planets f, g, and h requires more sophisticated modeling that can self-consistently balance heat production and transport in both rock and ice layers.

The Giant Impact is currently accepted as the leading theory for the formation of Earth's Moon. Successful scenarios for lunar origin should be able to explain the chemical composition of the Moon (volatile content and stable isotope ratios), the Moon's initial thermal state, and the system's bulk physical and dynamical properties. Hydrocode simulations of the formation of the Moon have long been able to match the bulk properties, but recent, more detailed work on the evolution of the protolunar disk has yielded great insight into the origin of the Moon's chemistry, and its early thermal history. Here, I show that the community has constructed the elements of an end-to-end theory for lunar origin that matches the overwhelming majority of observational constraints. In spite of the great progress made in recent years, new samples of the Moon, clarification of processes in the impact-generated disk, and a broader exploration of impact parameter space could yield even more insights into this fundamental and uniquely challenging geophysical problem.

The formation of satellites is thought to be a natural by-product of planet formation in our Solar System, and thus, moons of extrasolar planets (exomoons) may be abundant in extrasolar planetary systems, as well. Exomoons have yet to be discovered. However, moons larger than 0.1 Earth masses can be detected and characterized using current transit techniques. Here, we show that collisions between rocky planets with masses between a quarter to ten Earth masses can create impact-generated debris disks that could accrete into moons. Collisions between like-sized objects, at oblique impact angles, and velocities near escape speed create disks massive enough to form satellites that are dynamically stable against planetary tides. Impacts of this type onto a superearth between 2 to 7 Earth masses can launch into orbit enough mass to create a satellite large enough to be detected in Kepler transit data. Impact velocity is a crucial controlling factor on disk mass, which has been overlooked in all prior studies of moon formation via planetary collision.

Satellite formation is a natural by-product of planet formation. With the discovery of nu- merous extrasolar planets, it is likely that moons of extrasolar planets (exomoons) will soon be discovered. Some of the most promising techniques can yield both the mass and radius of the moon. Here, I review recent ideas about the formation of moons in our Solar System, and discuss the prospects of extrapolating these theories to predict the sizes of moons that may be discovered around extrasolar planets. It seems likely that planet-planet collisions could create satellites around rocky or icy planets which are large enough to be detected by currently available techniques. Detectable exomoons around gas giants may be able to form by co-accretion or capture, but determining the upper limit on likely moon masses at gas giant planets requires more detailed, modern simulations of these processes.

The New Horizons spacecraft has found evidence for geologic activity on the surface of Pluto, including extensional tectonic deformation of its water ice bedrock (see Moore et al., 2016). One mechanism that could drive extensional tectonic activity is global surface expansion due to the partial freezing of an ocean. We use updated physical properties for Pluto and simulate its thermal evolution to understand the survival of a possible subsurface ocean. For thermal conductivities of rock less than 3 W m$^{-1}$ K$^{-1}$, an ocean forms and at least partially freezes, leading to recent extensional stresses in the ice shell. In scenarios where the ocean freezes and the ice shell is thicker than $260$ km, ice II forms and causes global volume contraction. Since there is no evidence for recent compressional tectonic features, we argue that ice II has not formed and that Pluto's ocean has likely survived to present day.

Kuiper Belt objects with absolute magnitude less than 3 (radius $\gtrsim$500 km), the dwarf planets, have a range of different ice/rock ratios, and are more rock-rich than their smaller counterparts. Many of these objects have moons, which suggests that collisions may have played a role in modifying their compositions. We show that the dwarf planets fall into two categories when analysed by their mean densities and satellite-to-primary size ratio. Systems with large moons, such as Pluto/Charon and Orcus/Vanth, can form in low-velocity grazing collisions in which both bodies retain their compositions. We propose that these systems retain a primordial composition, with a density of about 1.8 g/cm$^3$. Triton, thought to be a captured KBO, could have lost enough ice during its early orbital evolution to explain its rock-enrichment relative to the primordial material. Systems with small moons, Eris, Haumea, and Quaoar, formed from a different type of collision in which icy material, perhaps a few tens of percent of the total colliding mass, is lost. The fragments would not remain in physical or dynamical proximity to the parent body. The ice loss process has not yet been demonstrated numerically, which could be due to the paucity of KBO origin simulations, or missing physical processes in the impact models. If our hypothesis is correct, we predict that large KBOs with small moons should be denser than the primordial material, and that the mean density of Orcus should be close to the primordial value.

Ridge-and-trough terrain is a common landform on outer Solar System icy satellites. Examples include Ganymede's grooved terrain, Europa's gray bands, Miranda's coronae, and several terrains on Enceladus. The conditions associated with the formation of each of these terrains are similar: heat flows of order tens to a hundred milliwatts per meter squared, and deformation rates of order $10^{-16}$ to $10^{-12}$ s$^{-1}$. Our prior work shows that the conditions associated with the formation of these terrains on Ganymede and the south pole of Enceladus are consistent with vigorous solid-state ice convection in a shell with a weak surface. We show that sluggish lid convection, an intermediate regime between the isoviscous and stagnant lid regimes, can create the heat flow and deformation rates appropriate for ridge and trough formation on a number of satellites, regardless of the ice shell thickness. For convection to deform their surfaces, the ice shells must have yield stresses similar in magnitude to the daily tidal stresses. Tidal and convective stresses deform the surface, and the spatial pattern of tidal cracking controls the locations of ridge-and-trough terrain.

We present a Radio Recombination Line (RRL) survey of the Galactic Plane from the HI Parkes All-sky Survey and associated Zone of Avoidance survey, which mapped the region l=196degr -- 0degr --52degr and |b| < 5degr at 1.4 GHz and 14.4 arcmin resolution. We combine three RRLs, H168$\alpha$, H167$\alpha$, and H166$\alpha$ to derive fully sampled maps of the diffuse ionized emission along the inner Galactic plane. The velocity information, at a resolution of 20 km/s, allows us to study the spatial distribution of the ionized gas and compare it with that of the molecular gas, as traced by CO. The longitude-velocity diagram shows that the RRL emission is mostly associated with CO gas from the molecular ring and is concentrated within the inner 30degr of longitude. A map of the free-free emission in this region of the Galaxy is derived from the line-integrated RRL emission, assuming an electron temperature gradient with Galactocentric radius of $496\pm100$ K/kpc. Based on the thermal continuum map we extracted a catalogue of 317 compact (<15 arcmin) sources, with flux densities, sizes and velocities. We report the first RRL observations of the southern ionized lobe in the Galactic centre. The line profiles and velocities suggest that this degree-scale structure is in rotation. We also present new evidence of diffuse ionized gas in the 3-kpc arm. Helium and carbon RRLs are detected in this survey. The He line is mostly observed towards HII regions, whereas the C line is also detected further away from the source of ionization. These data represent the first observations of diffuse C RRLs in the Galactic plane at a frequency of 1.4 GHz.

Iapetus, the outermost regular satellite of Saturn, has a drastic albedo dichotomy and an equatorial circumferential ridge that reaches heights of 20 km and widths of 70 km. This moon is thought to have formed concurrently with Saturn, and so would have experienced an intense bombardment after its formation. The ridge, which has been inferred to be one of the most ancient features on Iapetus' surface, could reasonably be expected to have been eroded by impacts; however, it has retained long continuous sections and a nearly pristine triangular shape with ridge slopes reaching $\sim$ 40$^{\circ}$. We use these observations, along with crater counts on Iapetus' surface, to constrain the total bombardment mass experienced by the satellite since its formation. The ridge morphology and the crater population recorded on Iapetus indicate it received less than 20% of the bombardment predicted by the classic Nice model for early Solar System evolution. Under the recently proposed scenarios of planetsimal-driven migration of the young outer planets including more realistic disk conditions, our results would imply a planetesimal disk mass of $M_{D}\sim12-34M_{\oplus}$.

Saturn's moon Iapetus has an equatorial ridge system, which may be as high as 20 km, that may have formed by endogenic forces, such as tectonic and convective forces, or exogenic processes such as debris infall. We use high-resolution topographic data to conduct a topographic analysis of the ridge, which suggests a predominantly triangular morphology, with some ridge face slopes reaching 40 degrees, allowing for an exogenic formation mechanism.

The Pluto-Charon system, likely formed from an impact, has reached the endpoint of its tidal evolution. During its evolution into the dual-synchronous state, the equilibrium tidal figures of Pluto and Charon would have also evolved as angular momentum was transferred from Pluto's spin to Charon's orbit. The rate of tidal evolution is controlled by Pluto's interior physical and thermal state. We examine three interior models for Pluto: an undifferentiated rock/ice mixture, differentiated with ice above rock, and differentiated with an ocean. For the undifferentiated case without an ocean, the Pluto-Charon binary does not evolve to its current state unless its internal temperature $T_i>200$ K, which would likely lead to strong tidal heating, melting, and differentiation. Without an ocean, Pluto's interior temperature must be higher than 240 K for Charon to evolve on a time scale less than the age of the solar system. Further tidal heating would likely create an ocean. If New Horizons finds evidence of ancient tidally-driven tectonic activity on either body, the most likely explanation is that Pluto had an internal ocean during Charon's orbital evolution.

We report about an ongoing photometric and spectroscopic monitoring survey of about 250 O- and 540 B-type stars in the southern Milky Way with the aim to determine the fraction of close binary systems as a function of mass and to determine the physical parameters of the individual components in the multiple systems. Preliminary results suggest that the multiplicity rate drops from 80% for the highest masses to 20% for stars of 3 solar masses. Our analysis indicates that the binary systems often contain close pairs with components of similar mass. This coincidence cannot originate from a random tidal capture in a dense cluster but is likely due to a particular formation process for high-mass stars. The large percentage of multiple systems requires a new photometric calibration for the absolute magnitudes of O-type stars.

We present a scenario for building the equatorial ridge and de-spinning Iapetus through an impact-generated disk and satellite. This impact puts debris into orbit, forming a ring inside the Roche limit and a satellite outside. This satellite rapidly pushes the ring material down to the surface of Iapetus, and then itself tidally evolves outward, thereby helping to de-spin Iapetus. This scenario can de-spin Iapetus an order of magnitude faster than when tides due to Saturn act alone, almost independently of its interior geophysical evolution. Eventually, the satellite is stripped from its orbit by Saturn. The range of satellite and impactor masses required is compatible with the estimated impact history of Iapetus.

This document covers five broad topics from theory of small bodies: planetesimal formation, cosmochemistry, thermal evolution, collisions and dynamics. Each of these topics is described in a separate section, where we prioritize selected main issues over completeness. The text points out the principal unresolved problems in each area, and suggests ways how progress can be made. This includes support for new code development, observations and experimental work that can be used to constrain theory, new research directions, and studies of cross-over regimes where the system's behavior is determined by several competing processes. The suggested development areas are placed in the context of NASA space exploration.

Asteroid pairs sharing similar heliocentric orbits were found recently. Backward integrations of their orbits indicated that they separated gently with low relative velocities, but did not provide additional insight into their formation mechanism. A previously hypothesized rotational fission process4 may explain their formation - critical predictions are that the mass ratios are less than about 0.2 and, as the mass ratio approaches this upper limit, the spin period of the larger body becomes long. Here we report photometric observations of a sample of asteroid pairs revealing that primaries of pairs with mass ratios much less than 0.2 rotate rapidly, near their critical fission frequency. As the mass ratio approaches 0.2, the primary period grows long. This occurs as the total energy of the system approaches zero requiring the asteroid pair to extract an increasing fraction of energy from the primary's spin in order to escape. We do not find asteroid pairs with mass ratios larger than 0.2. Rotationally fissioned systems beyond this limit have insufficient energy to disrupt. We conclude that asteroid pairs are formed by the rotational fission of a parent asteroid into a proto-binary system which subsequently disrupts under its own internal system dynamics soon after formation.