Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.

On 28th January 2018, the large Trans-Neptunian Object (TNO) 2002TC302 occulted a m$_v= $15.3 star with ID 130957813463146112 in the Gaia DR2 stellar catalog. 12 positive occultation chords were obtained from Italy, France, Slovenia and Switzerland. Also, 4 negative detections were obtained near the north and south limbs. This represents the best observed stellar occultation by a TNO other than Pluto, in terms of the number of chords published thus far. From the 12 chords, an accurate elliptical fit to the instantaneous projection of the body, compatible with the near misses, can be obtained. The resulting ellipse has major and minor axes of 543 $\pm$ 18 km and 460 $\pm$ 11 km, respectively, with a position angle of 3 $\pm$ 1 degrees for the minor axis. This information, combined with rotational light curves obtained with the 1.5m telescope at Sierra Nevada Observatory and the 1.23m telescope at Calar Alto observatory, allows us to derive possible 3D shapes and density estimations for the body, based on hydrostatic equilibrium assumptions. The effective area equivalent diameter is $\sim$ 84 km smaller than the radiometrically derived diameter using thermal data from Herschel and Spitzer Space Telescopes. This might indicate the existence of an unresolved satellite of up to $\sim$ 300 km in diameter, to account for all the thermal flux, although the occultation and thermal diameters are compatible within their error bars given the considerable uncertainty of the thermal results. The existence of a potential satellite also appears to be consistent with other ground-based data presented here. From the effective occultation diameter combined with H$_V$ measurements we derive a geometric albedo of 0.147 $\pm$ 0.005, which would be somewhat smaller if 2002TC302 has a satellite. The best occultation light curves do not show any signs of ring features or any signatures of a global atmosphere.

Non-resolved thermal infrared observations enable studies of thermal and physical properties of asteroid surfaces provided the shape and rotational properties of the target are well determined via thermo-physical models. We used calibration-programme Herschel PACS data (70, 100, 160 $\mu$m) and state-of-the-art shape models derived from adaptive-optics observations and/or optical light curves to constrain for the first time the thermal inertia of twelve large main-belt asteroids. We also modelled previously well-characterised targets such as (1) Ceres or (4) Vesta as they constitute important benchmarks. Using the scale as a free parameter, most targets required a re-scaling $\sim$5\% consistent with what would be expected given the absolute calibration error bars. This constitutes a good cross-validation of the scaled shape models, although some targets required larger re-scaling to reproduce the IR data. We obtained low thermal inertias typical of large main belt asteroids studied before, which continues to give support to the notion that these surfaces are covered by fine-grained insulating regolith. Although the wavelengths at which PACS observed are longwards of the emission peak for main-belt asteroids, they proved to be extremely valuable to constrain size and thermal inertia and not too sensitive to surface roughness. Finally, we also propose a graphical approach to help examine how different values of the exponent used for scaling the thermal inertia as a function of heliocentric distance (i.e. temperature) affect our interpretation of the results.

Thanks to the Gaia mission, it will be possible to determine the masses of approximately hundreds of large main belt asteroids with very good precision. We currently have diameter estimates for all of them that can be used to compute their volume and hence their density. However, some of those diameters are still based on simple thermal models, which can occasionally lead to volume uncertainties as high as 20-30%. The aim of this paper is to determine the 3D shape models and compute the volumes for 13 main belt asteroids that were selected from those targets for which Gaia will provide the mass with an accuracy of better than 10%. We used the genetic Shaping Asteroids with Genetic Evolution (SAGE) algorithm to fit disk-integrated, dense photometric lightcurves and obtain detailed asteroid shape models. These models were scaled by fitting them to available stellar occultation and/or thermal infrared observations. We determine the spin and shape models for 13 main belt asteroids using the SAGE algorithm. Occultation fitting enables us to confirm main shape features and the spin state, while thermophysical modeling leads to more precise diameters as well as estimates of thermal inertia values. We calculated the volume of our sample of main-belt asteroids for which the Gaia satellite will provide precise mass determinations. From our volumes, it will then be possible to more accurately compute the bulk density, which is a fundamental physical property needed to understand the formation and evolution processes of small solar system bodies.

With an estimated diameter in the 320 to 350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with $D$ > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with $D$ $\lesssim$ 200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller ($D$ $\lesssim$ 200 km) bodies. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling (ADAM) algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Interamnia's volume-equivalent diameter of 332 $\pm$ 6 km implies a bulk density of $\rho$=1.98 $\pm$ 0.68 gcm$^{-3}$ , which suggests that Interamnia - like Ceres and Hygiea - contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2 $\sigma$ level. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller ($D$ $\lesssim$ 300km) less massive ($m$ $\lesssim$ 3x10$^{19}$ kg) bodies.

We present a new non-convex model of the binary asteroid (809) Lundia. A SAGE (Shaping Asteroids with Genetic Evolution) method using disc-integrated photometry only was used for deriving physical parameters of this binary system. The model of (809) Lundia improves former system's pole solution and gives the ecliptic coordinates of the orbit pole - $\lambda=122^{\circ}$, $\beta=22^{\circ}$, $\sigma=\pm5^{\circ}$ - and the orbital period of $15.41574 \pm 0.00001$ h. For scaling our results we used effective diameter of $D_{eff} = 9.6 \pm 1.1$ km obtained from Spitzer observations. The non-convex shape description of the components permitted a refined calculation of the components' volumes, leading to a density estimation of $2.5\pm0.2$ g/cm$^3$ and macroporosity of 13-23\%. The intermediate-scale features of the model may also offer new clues on the components' origin and evolution.

The need for more accurate asteroid models is perhaps secondary to the need to measure their quality. The uncertainties of models' parameters propagate to quantities like volume or density -- the most important and informative properties of asteroids -- affecting conclusions about their physical nature. Our knowledge on shapes and spins of small solar system bodies comes mostly from visual, disk-integrated photometry. In this work we present a method for asteroid model uncertainty assessment based on visual photometry (lightcurves and sparse-in-time absolute measurements) allowing the determination of realistic volume uncertainty, as well as spin axis orientation, rotational period and local surface features. The sensitivity analysis is conducted by creating clones of the nominal model and accepting the ones that fit the observations within a confidence level. The uncertainties of model parameters are extracted from the extreme values found in the accepted clone population. Creation of such population of clones enables the conversion of a deterministic asteroid model into stochastic one, and can be utilized to create observation predictions with error bars. The method was used to assess the uncertainties of fictitious test models and real targets, i.e. (21) Lutetia, (89) Julia, (243) Ida, (433) Eros and (162173) Ryugu. We conclude that volumes, and subsequently, densities of asteroids derived from lightcurve-based models likely have vastly understated uncertainties, the biggest source of which is the inability to establish the extent of the model along its spin axis.

In this work we present SAGE (Shaping Asteroid models using Genetic Evolution) asteroid modelling algorithm based solely on photometric lightcurve data. It produces non-convex shapes, rotation axes orientati and rotational periods of asteroids. The main concept behind a genetic evolution algorithm is to produce random populations of shapes and spin orientations by mutating a seed shape and iterating the process until it converges to a stable global minimum. To test SAGE we have performed tes on five artificial shapes. We have also modelled (433) Eros and (9) Meti asteroids, as ground truth observations for them exist, allowing us to validate the models. We have compared derived Eros shape with NEAR Shoem model and Metis shape with adaptive optics and stellar occultation observations as with other available Metis models from various inversion methods.

Context. Over the past decades, several interplanetary missions have studied small bodies in situ, leading to major advances in our understanding of their geological and geophysical properties. These missions, however, have had a limited number of targets. Among them, the NASA Dawn mission has characterised in detail the topography and albedo variegation across the surface of asteroid (4) Vesta down to a spatial resolution of approximately 20 m/pixel scale. Aims. Here our aim was to determine how much topographic and albedo information can be retrieved from the ground with VLT/SPHERE in the case of Vesta, having a former space mission (Dawn) providing us with the ground truth that can be used as a benchmark. Methods. We observed Vesta with VLT/SPHERE/ZIMPOL as part of our ESO large programme (ID 199.C-0074) at six different epochs, and deconvolved the collected images with a parametric point spread function (PSF). We then compared our images with synthetic views of Vesta generated from the 3D shape model of the Dawn mission, on which we projected Vesta's albedo information. Results. We show that the deconvolution of the VLT/SPHERE images with a parametric PSF allows the retrieval of the main topographic and albedo features present across the surface of Vesta down to a spatial resolution of approximately 20-30 km. Contour extraction shows an accuracy of approximately 1 pixel (3.6 mas). The present study provides the very first quantitative estimate of the accuracy of ground-based adaptive-optics imaging observations of asteroid surfaces. Conclusions. In the case of Vesta, the upcoming generation of 30-40m telescopes (ELT, TMT, GMT) should in principle be able to resolve all of the main features present across its surface, including the troughs and the north-south crater dichotomy, provided that they operate at the diffraction limit.

The available set of spin and shape modelled asteroids is strongly biased against slowly rotating targets and those with low lightcurve amplitudes. As a consequence of these selection effects, the current picture of asteroid spin axis distribution, rotation rates, or radiometric properties, might be affected too. To counteract these selection effects, we are running a photometric campaign of a large sample of main belt asteroids omitted in most previous studies. We determined synodic rotation periods and verified previous determinations. When a dataset for a given target was sufficiently large and varied, we performed spin and shape modelling with two different methods. We used the convex inversion method and the non-convex SAGE algorithm, applied on the same datasets of dense lightcurves. Unlike convex inversion, the SAGE method allows for the existence of valleys and indentations on the shapes based only on lightcurves. We obtained detailed spin and shape models for the first five targets of our sample: (159) Aemilia, (227) Philosophia, (329) Svea, (478) Tergeste, and (487) Venetia. When compared to stellar occultation chords, our models obtained an absolute size scale and major topographic features of the shape models were also confirmed. When applied to thermophysical modelling, they provided a very good fit to the infrared data and allowed their size, albedo, and thermal inertia to be determined. Convex and non-convex shape models provide comparable fits to lightcurves. However, some non-convex models fit notably better to stellar occultation chords and to infrared data in sophisticated thermophysical modelling (TPM). In some cases TPM showed strong preference for one of the spin and shape solutions. Also, we confirmed that slowly rotating asteroids tend to have higher-than-average values of thermal inertia.

The combination of visible and thermal data from the ground and astrophysics space missions is key to improving the scientific understanding of near-Earth, main-belt, trojans, centaurs, and trans-Neptunian objects. To get full information on a small sample of selected bodies we combine different methods and techniques: lightcurve inversion, stellar occultations, thermophysical modeling, radiometric methods, radar ranging and adaptive optics imaging. The SBNAF project will derive size, spin and shape, thermal inertia, surface roughness, and in some cases bulk densities and even internal structure and composition, for objects out to the most distant regions in the Solar System. The applications to objects with ground-truth information allows us to advance the techniques beyond the current state-of-the-art and to assess the limitations of each method. We present results from our project's first phase: the analysis of combined Herschel-KeplerK2 data and Herschel-occultation data for TNOs; synergy studies on large MBAs from combined high-quality visual and thermal data; establishment of well-known asteroids as celestial calibrators for far-infrared, sub-millimetre, and millimetre projects; first results on near-Earth asteroids properties from combined lightcurve, radar and thermal measurements, as well as the Hayabusa-2 mission target characterisation. We also introduce public web-services and tools for studies of small bodies in general.

Among asteroids there exist ambiguities in their rotation period determinations. They are due to incomplete coverage of the rotation, noise and/or aliases resulting from gaps between separate lightcurves. To help to remove such uncertainties, basic characteristic of the lightcurves resulting from constraints imposed by the asteroid shapes and geometries of observations should be identified. We simulated light variations of asteroids which shapes were modelled as Gaussian random spheres, with random orientations of spin vectors and phase angles changed every $5^\circ$ from $0^\circ$ to $65^\circ$. This produced 1.4 mln lightcurves. For each simulated lightcurve Fourier analysis has been made and the harmonic of the highest amplitude was recorded. From the statistical point of view, all lightcurves observed at phase angles $\alpha < 30^\circ$, with peak-to-peak amplitudes $A>0.2$ mag are bimodal. Second most frequently dominating harmonic is the first one, with the 3rd harmonic following right after. For 1% of lightcurves with amplitudes $A < 0.1$ mag and phase angles $\alpha < 40^\circ$ 4th harmonic dominates.

Despite the large amount of high quality data generated in recent space encounters with asteroids, the majority of our knowledge about these objects comes from ground based observations. Asteroids travelling in orbits that are potentially hazardous for the Earth form an especially interesting group to be studied. In order to predict their orbital evolution, it is necessary to investigate their physical properties. This paper briefly describes the data requirements and different techniques used to solve the lightcurve inversion problem. Although photometry is the most abundant type of observational data, models of asteroids can be obtained using various data types and techniques. We describe the potential of radar imaging and stellar occultation timings to be combined with disk-integrated photometry in order to reveal information about physical properties of asteroids.

We present a new non-convex model of the 90 Antiope binary asteroid, derived with a modified version of the SAGE (Shaping Asteroids with Genetic Evolution) method using disk-integrated photometry only. A new variant of the SAGE algorithm capable of deriving models of binary systems is described. The model of 90 Antiope confirms the system's pole solution ($\lambda=199^{\circ}$, $\beta=38^{\circ}$, $\sigma=\pm5^{\circ}$) and the orbital period ($16.505046 \pm 0.000005$ h). A comparison between the stellar occultation chords obtained during the 2011 occultation and the projected shape solution has been used to scale the model. The resulting scaled model allowed us to obtain the equivalent radii ($R_{1}=40.4\pm0.9$ km and $R_{2}=40.2\pm0.9$ km) and the distance between the two system components ($176\pm4$ km), leading to a total system mass of ($9.14\pm0.62$)$\cdot10^{17}$ kg. The non-convex shape description of the components permitted a refined calculation of the components' volumes, leading to a density estimation of $1.67\pm0.23$ g cm$^{-3}$. The intermediate-scale features of the model may also offer new clues on the components' origin and evolution.