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.

Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V$\sim$7-8) and large angular size of 0.33 arcsec during its apparitions. Iris is believed to belong to the category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system's collisional evolution. We recovered information about the shape and surface topography of Iris from disk-resolved VLT/SPHERE/ZIMPOL images acquired in the frame of our ESO large program. We used the All-Data Asteroid Modeling (ADAM) shape reconstruction algorithm to model the 3D shape of Iris, using optical disk-integrated data and disk-resolved images from SPHERE as inputs. We analyzed the SPHERE images to infer the asteroid's global shape and the morphology of its main craters. We present the volume-equivalent diameter D$_{{\rm eq}}$=214$\pm$5 km, and bulk density $\rho$=2.7$\pm$0.3 g cm$^{-3}$ of Iris. Its shape appears to be consistent with that of an oblate spheroid with a large equatorial excavation. We identified eight putative surface features 20--40 km in diameter detected at several epochs, which we interpret as impact craters. Craters on Iris have depth-to-diameter ratios that are similar to those of analogous 10 km craters on Vesta. The bulk density of Iris is consistent with that of its meteoritic analog, namely LL ordinary chondrites. Considering the absence of a collisional family related to Iris and the number of large craters on its surface, we suggest that its equatorial depression may be the remnant of an ancient (at least 3 Gyr) impact. Iris's shape further opens the possibility that large planetesimals formed as almost perfect oblate spheroids. Finally, we attribute the difference in crater morphology between Iris and Vesta to their different surface gravities.

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.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

Context. CM-like asteroids (Ch and Cgh classes) are a major population within the broader C-complex, encompassing about 10% of the mass of the main asteroid belt. Their internal structure has been predicted to be homogeneous, based on their compositional similarity as inferred from spectroscopy (Vernazza et al., 2016, AJ 152, 154) and numerical modeling of their early thermal evolution (Bland & Travis, 2017, Sci. Adv. 3, e1602514). Aims. Here we aim to test this hypothesis by deriving the density of the CM-like asteroid (41) Daphne from detailed modeling of its shape and the orbit of its small satellite. Methods. We observed Daphne and its satellite within our imaging survey with the Very Large Telescope extreme adaptive-optics SPHERE/ZIMPOL camera (ID 199.C-0074, PI P. Vernazza) and complemented this data set with earlier Keck/NIRC2 and VLT/NACO observations. We analyzed the dynamics of the satellite with our Genoid meta-heuristic algorithm. Combining our high-angular resolution images with optical lightcurves and stellar occultations, we determine the spin period, orientation, and 3-D shape, using our ADAM shape modeling algorithm. Results. The satellite orbits Daphne on an equatorial, quasi-circular, prograde orbit, like the satellites of many other large main-belt asteroids. The shape model of Daphne reveals several large flat areas that could be large impact craters. The mass determined from this orbit combined with the volume computed from the shape model implies a density for Daphne of 1.77+/-0.26 g/cm3 (3 {\sigma}). This density is consistent with a primordial CM-like homogeneous internal structure with some level of macroporosity (~17%). Conclusions. Based on our analysis of the density of Daphne and 75 other Ch/Cgh-type asteroids gathered from the literature, we conclude that the primordial internal structure of the CM parent bodies was homogeneous.

(3200) Phaethon is a compelling object as it has an asteroidal appearance and spectrum, produces a weak dust tail during perihelion at just 0.14 AU, and is the parent body of the Geminid Meteor Shower. A better understanding of the physical properties of Phaethon is needed to understand the nature of its current and previous activity, relationship to potential source populations, and to plan for the upcoming flyby of the DESTINY+ spacecraft of Phaethon in the 2020s. We performed rotationally-resolved spectroscopy of Phaethon at visible and near-infrared wavelengths (0.4-2.5 microns) in 2007 and 2017, respectively, to better understand its surface properties. The visible and near-infrared observations both spanned nearly a full rotation or more and were under similar observing geometries, covering the whole surface with the exception of the north pole. The visible wavelengths show blue slopes with only minor slope variations and no absorption features. The NIR data is minimally varying and concave upwards, from very blue to blue-neutral with increasing wavelength. We fit the short-wavelength tail of Phaethon's thermal emission and retrieve an average visible albedo of pv = 0.08 +/- 0.01, which is lower than previous measurements but plausible in light of the recent larger radar-measured diameter of Phaethon. We retrieve an average infrared beaming parameter of Phaethon of eta = 1.70 +/- 0.05, which is similar to previous results. We discuss the implications of Phaethon's visible and near-infrared spectrum as well as the lower albedo on its origin, source population, and evolutionary history.

Eos family was created during a catastrophic impact about 1.3 Gyr ago. Rotation states of individual family members contain information about the history of the whole population. We aim to increase the number of asteroid shape models and rotation states within the Eos collision family, as well as to revise previously published shape models from the literature. Such results can be used to constrain theoretical collisional and evolution models of the family, or to estimate other physical parameters by a thermophysical modeling of the thermal infrared data. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as input for the convex inversion method, and derive 3D shape models of asteroids together with their rotation periods and orientations of rotation axes. We present updated shape models for 15 asteroids and new shape model determinations for 16 asteroids. Together with the already published models from the publicly available DAMIT database, we compiled a sample of 56 Eos family members with known shape models that we used in our analysis of physical properties within the family. Rotation states of asteroids smaller than ~20 km are heavily influenced by the YORP effect, whilst the large objects more or less retained their rotation state properties since the family creation. Moreover, we also present a shape model and bulk density of asteroid (423) Diotima, an interloper in the Eos family, based on the disk-resolved data obtained by the Near InfraRed Camera (Nirc2) mounted on the W.M. Keck II telescope.

Context. Many databases on asteroid brightnesses (e.g. ALCDEF, WISE) are potential sources for extensive asteroid shape and spin modelling. Individual lightcurve inversion models require several apparitions and hundreds of data points per target. However, we can analyse the coarse shape and spin distributions over populations of at least thousands of targets even if there are only a few points and one apparition per asteroid. This is done by examining the distribution of the brightness variations observed within the chosen population. Aims. Brightness variation has been proposed as a population-scale rather than individual-target observable in two studies so far. We aim to examine this approach rigorously to establish its theoretical validity, degree of ill-posedness, and practical applicability. Methods. We model the observed brightness variation of a target population by considering its cumulative distribution function (CDF) caused by the joint distribution function of two fundamental shape and spin indicators. These are the shape elongation and the spin latitude of a simple ellipsoidal model. The main advantage of the model is that we can derive analytical basis functions that yield the observed CDF as a function of the shape and spin distribution. The inverse problem can be treated linearly. Even though the inaccuracy of the model is considerable, databases of thousands of targets should yield some information on the distribution. Results. We establish the theoretical soundness and the typical accuracy limits of the approach both analytically and numerically. Using simulations, we derive a practical estimate of the model distribution in the (shape, spin)-plane. We show that databases such as Wide-field Infrared Survey Explorer (WISE) yield coarse but robust estimates of this distribution, and as an example compare various asteroid families with each other.

Asteroid modeling efforts in the last decade resulted in a comprehensive dataset of almost 400 convex shape models and their rotation states. This amount already provided a deep insight into physical properties of main-belt asteroids or large collisional families. We aim to increase the number of asteroid shape models and rotation states. Such results are an important input for various further studies such as analysis of asteroid physical properties in different populations, including smaller collisional families, thermophysical modeling, and scaling shape models by disk-resolved images, or stellar occultation data. This provides, in combination with known masses, bulk density estimates, but constrains also theoretical collisional and evolutional models of the Solar System. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as an input for the convex inversion method, and derive 3D shape models of asteroids, together with their rotation periods and orientations of rotation axes. The key ingredient is the support of more that one hundred observers who submit their optical data to publicly available databases. We present updated shape models for 36 asteroids, for which mass estimates are currently available in the literature or their masses will be most likely determined from their gravitational influence on smaller bodies, which orbital deflection will be observed by the ESA Gaia astrometric mission. This was achieved by using additional optical data from recent apparitions for the shape optimization. Moreover, we also present new shape model determinations for 250 asteroids, including 13 Hungarias and 3 near-Earth asteroids.

In the analysis of thermal infrared data of asteroids by means of thermophysical models (TPMs) it is a common practice to neglect the uncertainty of the shape model and the rotational state, which are taken as an input for the model. Here, we present a novel method of investigating the importance of the shape model and the pole orientation uncertainties in the thermophysical modeling - the varied shape TPM (VS-TPM). Our method uses optical photometric data to generate various shape models that map the uncertainty in the shape and the rotational state. The TPM procedure is then run for all these shape models. We apply the implementation of the classical TPM as well as our VS-TPM to the convex shape models of several asteroids together with their thermal infrared data acquired by the NASA's Wide-field Infrared Survey Explorer (WISE) and compare the results. These show that the uncertainties of the shape model and the pole orientation can be very important (e.g., for the determination of the thermal inertia) and should be considered in the thermophysical analyses. We present thermophysical properties for six asteroids - (624) Hektor, (771) Libera, (1036) Ganymed, (1472) Muonio, (1627) Ivar, and (2606) Odessa.

Current amount of ~500 asteroid models derived from the disk-integrated photometry by the lightcurve inversion method allows us to study not only the spin-vector properties of the whole population of MBAs, but also of several individual collisional families. We create a data set of 152 asteroids that were identified by the HCM method as members of ten collisional families, among them are 31 newly derived unique models and 24 new models with well-constrained pole-ecliptic latitudes of the spin axes. The remaining models are adopted from the DAMIT database or the literature. We revise the preliminary family membership identification by the HCM method according to several additional criteria - taxonomic type, color, albedo, maximum Yarkovsky semi-major axis drift and the consistency with the size-frequency distribution of each family, and consequently we remove interlopers. We then present the spin-vector distributions for eight asteroidal families. We use a combined orbital- and spin-evolution model to explain the observed spin-vector properties of objects among collisional families. In general, we observe for studied families similar trends in the (a_p, \beta) space: (i) larger asteroids are situated in the proximity of the center of the family; (ii) asteroids with \beta>0{\deg} are usually found to the right from the family center; (iii) on the other hand, asteroids with \beta<0{\deg} to the left from the center; (iv) majority of asteroids have large pole-ecliptic latitudes (|\beta|\gtrsim 30{\deg}); and finally (v) some families have a statistically significant excess of asteroids with \beta>0{\deg} or \beta<0{\deg}. Our numerical simulation of the long-term evolution of a collisional family is capable of reproducing well the observed spin-vector properties. Using this simulation, we also independently constrain the age of families Flora (1.0\pm0.5 Gyr) and Koronis (2.5-4 Gyr).

We select 50 main-belt asteroids with a diameter between 20 and 400 km for which we have (i) shape models derived by the lightcurve inversion method (LI) and (ii) resolved observations of good quality collected with the Keck II adaptive optics (AO) system in the near-infrared. We derive the size of these asteroids by minimizing the difference between the contours from deconvolved AO images and the projected silhouettes calculated from the shape model at the time of the AO observations. We compute the volume-equivalent diameters for 48 of these asteroids. For 15 of them, we remove the ambiguity of the pole orientation typical for shape models derived by the LI. We have found that our equivalent diameters are smaller by 3%, 7%, and 2% compared with the effective diameters derived from mid-IR photometric observations provided by IRAS, WISE and AKARI. For 40 asteroids with previously determined mass estimates, we compute their bulk densities and discuss the mass-density dependence with respect to taxonomic types.

We investigate the potential of the sparse data produced by the Catalina Sky Survey astrometric project (CSS for short) in asteroid shape and rotational state determination by the lightcurve inversion method. We show that although the photometric quality of the CSS data, compared to the dense data, is significantly worse, it is in principle possible that these data are for some asteroids with high lightcurve amplitudes sufficient for a unique shape determination. CSS data are available for $\sim$180 asteroids for which shape models were previously derived from different photometric data sets. For 13 asteroids from this sample, we derive their unique shape models based only on CSS data, compare the two independent shape models together and discuss the reliability of models derived from only CSS data. We also use CSS data to determine shape models for asteroids with already known rotational period values, derive 12 unique models and compare previously published periods with periods determined from the full 3D modeling by the lightcurve inversion method. Finally, we test different shape resolutions used in the lightcurve inversion method in order to find reliable asteroid models.

The larger number of models of asteroid shapes and their rotational states derived by the lightcurve inversion give us better insight into both the nature of individual objects and the whole asteroid population. With a larger statistical sample we can study the physical properties of asteroid populations, such as main-belt asteroids or individual asteroid families, in more detail. Shape models can also be used in combination with other types of observational data (IR, adaptive optics images, stellar occultations), e.g., to determine sizes and thermal properties. We use all available photometric data of asteroids to derive their physical models by the lightcurve inversion method and compare the observed pole latitude distributions of all asteroids with known convex shape models with the simulated pole latitude distributions. We used classical dense photometric lightcurves from several sources and sparse-in-time photometry from the U.S. Naval Observatory in Flagstaff, Catalina Sky Survey, and La Palma surveys (IAU codes 689, 703, 950) in the lightcurve inversion method to determine asteroid convex models and their rotational states. We also extended a simple dynamical model for the spin evolution of asteroids used in our previous paper. We present 119 new asteroid models derived from combined dense and sparse-in-time photometry. We discuss the reliability of asteroid shape models derived only from Catalina Sky Survey data (IAU code 703) and present 20 such models. By using different values for a scaling parameter cYORP (corresponds to the magnitude of the YORP momentum) in the dynamical model for the spin evolution and by comparing synthetics and observed pole-latitude distributions, we were able to constrain the typical values of the cYORP parameter as between 0.05 and 0.6.