The physical properties of cometary nuclei observed today relate to their complex history and help to constrain their formation and evolution. In this article, we review some of the main physical properties of cometary nuclei and focus in particular on the thermal, mechanical, structural and dielectric properties, emphasizing the progress made during the Rosetta mission. Comets have a low density of 480 $\pm$ 220 kg m-3 and a low permittivity of 1.9 - 2.0, consistent with a high porosity of 70 - 80 %, are weak with a very low global tensile strength $<$100 Pa, and have a low bulk thermal inertia of 0 - 60 J K-1 m-2 s-1/2 that allowed them to preserve highly volatiles species (e.g. CO, CO2, CH4, N2) into their interior since their formation. As revealed by 67P/Churyumov-Gerasimenko, the above physical properties vary across the nucleus, spatially at its surface but also with depth. The broad picture is that the bulk of the nucleus consists of a weakly bonded, rather homogeneous material that preserved primordial properties under a thin shell of processed material, and possibly covered by a granular material; this cover might in places reach a thickness of several meters. The properties of the top layer (the first meter) are not representative of that of the bulk nucleus. More globally, strong nucleus heterogeneities at a scale of a few meters are ruled out on 67P small lobe.

Aims. We use four observational data sets, mainly from the Rosetta mission, to constrain the activity pattern of the nucleus of comet 67P/Churyumov-Gerasimenko. Methods. We develop a numerical model that computes the production rate and non-gravitational acceleration of the nucleus of comet 67P as a function of time, taking into account its complex shape with a shape model reconstructed from OSIRIS imagery. We use this model to fit three observational data sets: the trajectory data from flight dynamics; the rotation state, as reconstructed from OSIRIS imagery; and the water production measurements from ROSINA, of 67P. The two key parameters of our model, adjusted to fit the three data sets all together, are the activity pattern and the momentum transfer efficiency (i.e., the so-called "$\eta$ parameter" of the non-gravitational forces). Results. We find an activity pattern able to successfully reproduce the three data sets simultaneously. The fitted activity pattern exhibits two main features: a higher effective active fraction in two southern super-regions ($\sim 10$~\%) outside perihelion compared to the northern ones ($< 4$~\%), and a drastic rise of the effective active fraction of the southern regions ($\sim 25-35$~\%) around perihelion. We interpret the time-varying southern effective active fraction by cyclic formation and removal of a dust mantle in these regions. Our analysis supports moderate values of the momentum transfer coefficient $\eta$ in the range $0.6-0.7$; values $\eta\leq0.5$ or $\eta\geq0.8$ degrade significantly the fit to the three data sets. Our conclusions reinforce the idea that seasonal effects linked to the orientation of the spin axis play a key role in the formation and evolution of dust mantles, and in turn largely control the temporal variations of the gas flux.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

We simulate the stresses induced by temperature changes in a putative hard layer near the surface of comet 67P/Churyumov--Gerasimenko with a thermo-viscoelastic model. Such a layer could be formed by the recondensation or sintering of water ice (and dust grains), as suggested by laboratory experiments and computer simulations, and would explain the high compressive strength encountered by experiments on board the Philae lander. Changes in temperature from seasonal insolation variation penetrate into the comet's surface to depths controlled by the thermal inertia, causing the material to expand and contract. Modelling this with a Maxwellian viscoelastic response on a spherical nucleus, we show that a hard, icy layer with similar properties to Martian permafrost will experience high stresses: up to tens of MPa, which exceed its material strength (a few MPa), down to depths of centimetres to a metre. The stress distribution with latitude is confirmed qualitatively when taking into account the comet's complex shape but neglecting thermal inertia. Stress is found to be comparable to the material strength everywhere for sufficient thermal inertia ($\gtrsim50$ J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) and ice content ($\gtrsim 45\%$ at the equator). In this case, stresses penetrate to a typical depth of $\sim0.25$ m, consistent with the detection of metre-scale thermal contraction crack polygons all over the comet. Thermal fracturing may be an important erosion process on cometary surfaces which breaks down material and weakens cliffs.

We use a laboratory facility to study the sputtering properties of centimeter-thick porous water ice subjected to the bombardment of ions and electrons to better understand the formation of exospheres of the icy moons of Jupiter. Our ice samples are as similar as possible to the expected moon surfaces but surface charging of the samples during ion irradiation may distort the experimental results. We therefore monitor the time scales for charging and dis- charging of the samples when subjected to a beam of ions. These experiments allow us to derive an electric conductivity of deep porous ice layers. The results imply that electron irradiation and sputtering play a non-negligible role for certain plasma conditions at the icy moons of Jupiter. The observed ion sputtering yields from our ice samples are similar to previous experiments where compact ice films were sputtered off a micro-balance.

With hundreds of exoplanets detected, it is necessary to revisit giant planets accretion models to explain their mass distribution. In particular, formation of sub-jovian planets remains unclear, given the short timescale for the runaway accretion of massive atmospheres. However, gas needs to pass through a circum-planetary disc. If the latter has a low viscosity (as expected if planets form in "dead zones"), it might act as a bottleneck for gas accretion. We investigate what the minimum accretion rate is for a planet under the limit assumption that the circum-planetary disc is totally inviscid, and the transport of angular momentum occurs solely because of the gravitational perturbations from the star. To estimate the accretion rate, we present a steady-state model of an inviscid circum-planetary disc, with vertical gas inflow and external torque from the star. Hydrodynamical simulations of a circum-planetary disc were conducted in 2D, in a planetocentric frame, with the star as an external perturber in order to measure the torque exerted by the star on the disc. The disc shows a two-armed spiral wave caused by stellar tides, propagating all the way in from the outer edge of the disc towards the planet. The stellar torque is small and corresponds to a doubling time for a Jupiter mass planet of the order of 5 Myrs. Given the limit assumptions, this is clearly a lower bound of the real accretion rate. This result shows that gas accretion onto a giant planet can be regulated by circum-planetary discs. This suggests that the diversity of masses of extra-solar planets may be the result of different viscosities in these discs.