Search

The aim of our study is to characterise the spin states of the members of the Athor and Zita collisional families and test whether these members have a spin distribution consistent with a common origin from the break up of their respective family parent asteroids. Our method is based on the asteroid family evolution, which indicates that there should be a statistical predominance of retrograde-rotating asteroids on the inward side of family's V-shape, and prograde-rotating asteroids on the outward side. We used photometric data from our campaign and the literature in order to reveal the spin states of the asteroids belonging to these families. We combined dense and sparse photometric data in order to construct lightcurves; we performed the lightcurve inversion method to estimate the sidereal period, spin axis and convex shape of several family members. We obtained 34 new asteroid models for Athor family members and 17 for Zita family members. Along with the literature and revised models, the Athor family contains 60% of retrograde asteroids on the inward side and, 76% of prograde asteroids on the outward side. We also found that the Zita family exhibits 80% of retrograde asteroids on the inward side and an equal amount of prograde and retrograde rotators on the outward side. However, when we applied Kernel density estimation, we also found a clear peak for prograde asteroids on the outward side, as expected from the theory. The spin states of these asteroids validate the existence of both families, with the Athor family exhibiting a stronger signature for the presence of retrograde-rotating and prograde-rotating asteroids on the inner and outer side of the family, respectively. Our work provides an independent confirmation and characterisation of these very old families, whose presence and characteristics offer constraints for theories and models of the Solar System's evolution.

The rotation state of small asteroids is affected in the long term by perturbing torques of gravitational and radiative origin (the YORP effect). Direct observational evidence of the YORP effect is the primary goal of our work. We carried out photometric observations of five near-Earth asteroids: (1862) Apollo, (2100) Ra-Shalom, (85989) 1999 JD6, (138852) 2000 WN10, and (161989) Cacus. Then we applied the light-curve inversion method to all available data to determine the spin state and a convex shape model for each of the five studied asteroids. In the case of (2100) Ra-Shalom, the analysis required that the spin-axis precession due to the solar gravitational torque also be included. We obtained two new detections of the YORP effect: (i) $(2.9 \pm 2.0)\times 10^{-9}\,\mathrm{rad\,d}^{-2}$ for (2100) Ra-Shalom, and (ii) $(5.5\pm 0.7)\times 10^{-8}\,\mathrm{rad\,d}^{-2}$ for (138852) 2000 WN10. The analysis of Ra-Shalom also reveals a precession of the spin axis with a precession constant $\sim 3000''\,\mathrm{yr}^{-1}$. This is the first such detection from Earth-bound photometric data. For the other two asteroids, we improved the accuracy of the previously reported YORP detection: (i) $(4.94 \pm 0.09)\times 10^{-8}\,\mathrm{rad\,d}^{-2}$ for (1862) Apollo, and (ii) $(1.86\pm 0.09)\times 10^{-8}\,\mathrm{rad\,d}^{-2}$ for (161989) Cacus. Despite the recent report of a detected YORP effect for (85989) 1999 JD6, we show that the model without YORP cannot be rejected statistically. Therefore, the detection of the YORP effect for this asteroid requires future observations. The spin-axis precession constant of Ra-Shalom determined from observations matches the theoretically expected value. The total number of asteroids with a YORP detection has increased to 12. In all cases, the rotation frequency increases in time.

As evidenced by recent survey results, majority of asteroids are slow rotators (P>12 h), but lack spin and shape models due to selection bias. This bias is skewing our overall understanding of the spins, shapes, and sizes of asteroids, as well as of their other properties. Also, diameter determinations for large (>60km) and medium-sized asteroids (between 30 and 60 km) often vary by over 30% for multiple reasons. Our long-term project is focused on a few tens of slow rotators with periods of up to 60 hours. We aim to obtain their full light curves and reconstruct their spins and shapes. We also precisely scale the models, typically with an accuracy of a few percent. We used wide sets of dense light curves for spin and shape reconstructions via light-curve inversion. Precisely scaling them with thermal data was not possible here because of poor infrared data: large bodies are too bright for WISE mission. Therefore, we recently launched a campaign among stellar occultation observers, to scale these models and to verify the shape solutions, often allowing us to break the mirror pole ambiguity. The presented scheme resulted in shape models for 16 slow rotators, most of them for the first time. Fitting them to stellar occultations resolved previous inconsistencies in size determinations. For around half of the targets, this fitting also allowed us to identify a clearly preferred pole solution, thus removing the ambiguity inherent to light-curve inversion. We also address the influence of the uncertainty of the shape models on the derived diameters. Overall, our project has already provided reliable models for around 50 slow rotators. Such well-determined and scaled asteroid shapes will, e.g. constitute a solid basis for density determinations when coupled with mass information. Spin and shape models continue to fill the gaps caused by various biases., Comment: Accepted to Astronomy & Astrophysics. 12 pages + appendices

The Ch-type asteroid (130) Elektra is orbited by three moons, making it the first quadruple system in the main asteroid belt. We aim to characterise the irregular shape of Elektra and construct a complete orbital model of its unique moon system. We applied the All-Data Asteroid Modelling (ADAM) algorithm to 60 light curves of Elektra, including our new measurements, 46 adaptive-optics (AO) images obtained by the VLT/SPHERE and Keck/Nirc2 instruments, and two stellar occultation profiles. For the orbital model, we used an advanced $N$-body integrator, which includes a multipole expansion of the central body (with terms up to the order $\ell = 6$), mutual perturbations, internal tides, as well as the external tide of the Sun acting on the orbits. We fitted the astrometry measured with respect to the central body and also relatively, with respect to the moons themselves. We obtained a revised shape model of Elektra with the volume-equivalent diameter $(201\pm 2)\,{\rm km}$. Out of two pole solutions, $(\lambda, \beta) = (189; -88)\,{\rm deg}$ is preferred, because the other one leads to an incorrect orbital evolution of the moons. We also identified the true orbital period of the third moon S/2014 (130) 2 as $P_2 = (1.642112 \pm 0.000400)\,{\rm d}$, which is in between the other periods, $P_1 \simeq 1.212\,{\rm d}$, $P_3 \simeq 5.300\,{\rm d}$, of S/2014 (130) 1 and S/2003 (130) 1, respectively. The resulting mass of Elektra, $(6.606 \substack{+0.007 \\ -0.013}) \times 10^{18}\,{\rm kg}$, is precisely constrained by all three orbits. Its bulk density is then $(1.536 \pm 0.038)\,{\rm g\,cm}^{-3}$. The expansion with the assumption of homogeneous interior leads to the oblateness $J_2 = -C_{20} \simeq 0.16$. However, the best-fit precession rates indicate a slightly higher value, ${\simeq}\,0.18$., Comment: 19 pages, 23 figures

The satellite Linus orbiting the main-belt asteroid (22) Kalliope exhibited occultation and transit events in late 2021. A photometric campaign was organized and observations were taken by the TRAPPIST-South, SPECULOOS-Artemis, OWL-Net, and BOAO telescopes, with the goal to constrain models of this system. Our dynamical model is complex, with multipoles (up to the order $\ell = 2$), internal tides, and external tides. The model was constrained by astrometry (spanning 2001--2021), occultations, adaptive-optics imaging, calibrated photometry, as well as relative photometry. Our photometric model was substantially improved. A new precise (${<}\,0.1\,{\rm mmag}$) light curve algorithm was implemented, based on polygon intersections, which are computed exactly -- by including partial eclipses and partial visibility of polygons. Moreover, we implemented a `cliptracing' algorithm, based again on polygon intersections, in which partial contributions to individual pixels are computed exactly. Both synthetic light curves and synthetic images are then very smooth. Based on our combined solution, we confirmed the size of Linus, $(28\pm 1)\,{\rm km}$. However, this solution exhibits some tension between the light curves and the PISCO speckle-interferometry dataset. In most solutions, Linus is darker than Kalliope, with the albedos $A_{\rm w} = 0.40$ vs. $0.44$. This is confirmed on deconvolved images. A~detailed revision of astrometric data allowed us to revise also the $J_2 \equiv -C_{20}$ value of Kalliope. Most importantly, a~homogeneous body is excluded. For a differentiated body, two solutions exist: low-oblateness ($C_{20} \simeq -0.12$), with a~spherical iron core, and alternatively, high-oblateness ($C_{20} \simeq -0.22$) with an elongated iron core. These correspond to the low- and high-energy collisions, respectively, studied by means of SPH simulations in our previous work., Comment: Astronomy and Astrophysics, accepted

Cybele asteroids constitute an appealing reservoir of primitive material genetically linked to the outer Solar System, and the physical properties of the largest members can be readily accessed by large telescopes. We took advantage of the bright apparition of (65) Cybele in July and August 2021 to acquire high-angular-resolution images and optical light curves of the asteroid with which we aim to analyse its shape and bulk properties. 7 series of images acquired with VLT/SPHERE were combined with optical light curves to reconstruct the shape of the asteroid using the ADAM, MPCD, and SAGE algorithms. The origin of the shape was investigated by means of N-body simulations. Cybele has a volume-equivalent diameter of 263+/-3km and a bulk density of 1.55+/-0.19g.cm-3. Notably, its shape and rotation state are closely compatible with those of a Maclaurin equilibrium figure. The lack of a collisional family associated with Cybele and the higher bulk density of that body with respect to other large P-type asteroids suggest that it never experienced any large disruptive impact followed by rapid re-accumulation. This would imply that its present-day shape represents the original one. However, numerical integration of the long-term dynamical evolution of a hypothetical family shows that it is dispersed by gravitational perturbations and chaotic diffusion over Gyrs of evolution. The very close match between Cybele and an equilibrium figure opens up the possibility that D>260km small bodies from the outer Solar System all formed at equilibrium. However, we cannot rule out an old impact as the origin of the equilibrium shape. Cybele itself is found to be dynamically unstable, implying that it was recently (<1Ga) placed on its current orbit either through slow diffusion from a relatively stable orbit in the Cybele region or, less likely, from an unstable, JFC orbit in the planet-crossing region., Comment: 19 pages, 14 figures, 4 tables, accepted for publication in A&A

Families of asteroids generated by the collisional fragmentation of a common parent body have been identified using clustering methods of asteroids in their proper orbital element space. An alternative method has been developed in order to identify collisional families from the correlation between the asteroid fragment sizes and their proper semi-major axis distance from the family centre (V-shape). This method has been shown to be effective in the cases of the very diffuse families that formed billions of years ago. We obtained photometric observations of asteroids in order to construct their rotational light curves; we combine them with the literature light curves and sparse-in-time photometry; we input these data in the light curve inversion methods to determine the shape and the spin pole of the asteroids in order to assess whether an object is prograde or retrograde. The ultimate goal is to assess whether we find an excess of retrograde asteroids on the inward side of the V-shape of a 4 Gyr asteroid family identified via the V-shape method. This excess of retrograde rotators is predicted by the theory of asteroid family evolution. We obtained the spin poles for 55 asteroids claimed to belong to a 4 Gyr collisional family of the inner main belt that consists of low-albedo asteroids. After re-evaluating the albedo and spectroscopic information, we found that nine of these asteroids are interlopers in the 4 Gyr family. Of the 46 remaining asteroids, 31 are found to be retrograde and 15 prograde. We also found that these retrograde rotators have a very low probability (1.29%) of being due to random sampling from an underlying uniform distribution of spin poles. Our results constitute corroborating evidence that the asteroids identified as members of a 4 Gyr collisional family have a common origin, thus strengthening their family membership., Comment: 20 pages, 7 figures. A&A (2022)

The orbit of the outer satellite Alexhelios of (216) Kleopatra is already constrained by adaptive-optics astrometry, obtained with the VLT/SPHERE instrument. However, there is also a preceding occultation event in 1980 attributed to this satellite. Hereinafter, we try to link all observations, spanning 1980--2018. We find the nominal orbit exhibits an unexplained shift by $+60^\circ$ in the true longitude. Using both periodogram analysis and an $\ell = 10$ multipole model suitable for the motion of mutually interacting moons about the irregular body, we confirmed that it is not possible to adjust the respective osculating period $P_2$. Instead, we were forced to use a model with tidal dissipation (and increasing orbital periods) to explain the shift. We also analyzed light curves, spanning 1977--2021, and searched for the expected spin deceleration of Kleopatra. According to our best-fit model, the observed period rate is $\dot P_2 = (1.8\pm 0.1)\cdot 10^{-8}\,{\rm d}\,{\rm d}^{-1}$ and the corresponding time lag $\Delta t_2 = 42\,{\rm s}$ of tides, for the assumed value of the Love number $k_2 = 0.3$. It is the first detection of tidal evolution for moons orbiting 100-km asteroids. The corresponding dissipation factor $Q$ is comparable with other terrestrial bodies, albeit at a higher loading frequency $2|\omega-n|$. We also predict a secular evolution of the inner moon, $\dot P_1 = 5.0\cdot 10^{-8}$, as well as a spin deceleration of Kleopatra, $\dot P_0 = 1.9\cdot 10^{-12}$. In alternative models, with moons captured in the 3:2 mean-motion resonance or more massive moons, the respective values of $\Delta t_2$ are a factor of 2--3 lower. Future astrometric observations by direct imaging or occultations should allow to distinguish between these models, which is important for the internal structure and mechanical properties of (216) Kleopatra., Comment: accepted in A&A

Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of 5 g.cm$^{-3}$. Such a high density implies a high metal content and a low porosity which is not easy to reconcile with its peculiar dumbbell shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, we aim to constrain the mass and the shape of Kleopatra with high accuracy, hence its density. Methods. We combined our new VLT/SPHERE observations of Kleopatra recorded in 2017 and 2018 with archival data, as well as lightcurve, occultation, and delay-Doppler images, to derive its 3D shape model using two different algorithms (ADAM, MPCD). Furthermore, an N-body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper. Results. The shape of Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75$\pm$1.40) km and mass (2.97$\pm$0.32) 10$^{18}$ kg imply a bulk density of (3.38$\pm$0.50) g cm$^{-3}$. Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near-equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. Kleopatra's current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding. Conclusions. Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. It deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission., Comment: Accepted for publication in Astronomy and Astrophysics

We present $V$-band photometry of the 20,000 brightest asteroids using data from the All-Sky Automated Survey for Supernovae (ASAS-SN) between 2012 and 2018. We were able to apply the convex inversion method to more than 5,000 asteroids with more than 60 good measurements in order to derive their sidereal rotation periods, spin axis orientations, and shape models. We derive unique spin state and shape solutions for 760 asteroids, including 163 new determinations. This corresponds to a success rate of about 15%, which is significantly higher than the success rate previously achieved using photometry from surveys. We derive the first sidereal rotation periods for additional 69 asteroids. We find good agreement in spin periods and pole orientations for objects with prior solutions. We obtain a statistical sample of asteroid physical properties that is sufficient for the detection of several previously known trends, such as the underrepresentation of slow rotators in current databases, and the anisotropic distribution of spin orientations driven by the nongravitational forces. We also investigate the dependence of spin orientations on the rotation period. Since 2018, ASAS-SN has been observing the sky with higher cadence and deeper limiting magnitude, which will lead to many more new solutions in just a few years., Comment: Accepted for publication in Astronomy and Astrophysics

To interpret adaptive-optics observations of (216) Kleopatra, we need to describe an evolution of multiple moons, orbiting an extremely irregular body and including their mutual interactions. Such orbits are generally non-Keplerian and orbital elements are not constants. Consequently, we use a modified $N$-body integrator, which was significantly extended to include the multipole expansion of the gravitational field up to the order $\ell = 10$. Its convergence was verified against the `brute-force' algorithm. We computed the coefficients $C_{\ell m},S_{\!\ell m}$ for Kleopatra's shape, assuming a~constant bulk density. For solar-system applications, it was also necessary to implement a variable distance and geometry of observations. Our $\chi^2$ metric then accounts for the absolute astrometry, the relative astrometry (2nd moon with respect to 1st), angular velocities, and also silhouettes, constraining the pole orientation. This allowed us to derive the orbital elements of Kleopatra's two moons. Using both archival astrometric data and new VLT/SPHERE observations (ESO LP 199.C-0074), we were able to identify the true periods of the moons, $P_1 = (1.822359\pm0.004156)\,{\rm d}$, $P_2 = (2.745820\pm0.004820)\,{\rm d}$. They orbit very close to the 3:2 mean-motion resonance, but their osculating eccentricities are too small compared to other perturbations (multipole, mutual), so that regular librations of the critical argument are not present. The resulting mass of Kleopatra, $m_1 = (1.49\pm0.16)\cdot10^{-12}\,M_\odot$ or $2.97\cdot10^{18}\,{\rm kg}$, is significantly lower than previously thought. An implication explained in the accompanying paper (Marchis et al.) is that (216) Kleopatra is a critically rotating body., Comment: accepted in A&A

The largest asteroids in the Koronis family (sizes $\geq 25$ km) have very peculiar rotation state properties, with the retrograde- and prograde-rotating objects being distinctly different. A recent e-analysis of observations suggests that one of the asteroids formerly thought to be retrograde-rotating, 208~Lacrimosa, in reality exhibits prograde rotation, yet other properties of this object are discrepant with other members this group. We seek to understand whether the new spin solution of Lacrimosa invalidates the previously proposed model of the Koronis large members or simply reveals more possibilities for the long-term evolutionary paths, including some that have not yet been explored. We confirm and substantiate the previously suggested prograde rotation of Lacrimosa. Its spin vector has an ecliptic longitude and latitude of $(\lambda,\beta)=(15^\circ \pm 2^\circ, 67^\circ\pm 2^\circ)$ and a sidereal rotation period $P=14.085734\pm 0.000007$ hr. The thermal and occultation data allow us to calibrate a volume equivalent size of $D=44\pm 2$ km of Lacrimosa. The observations also constrain the shape model relatively well. Assuming uniform density, the dynamical ellipticity is $\Delta=0.35\pm 0.05$. Unlike other large prograde-rotating Koronis members, Lacrimosa spin is not captured in the Slivan state. We propose that Lacrimosa differed from this group in that it had initially slightly larger obliquity and longer rotation period. With those parameters, it jumped over the Slivan state instead of being captured and slowly evolved into the present spin configuration. In the future, it is likely to be captured in the Slivan state corresponding to the proper (instead of forced) mode of the orbital plane precession in the inertial space., Comment: 18 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

Dynamical models of Solar System evolution have suggested that P-/D-type volatile-rich asteroids formed in the outer Solar System and may be genetically related to the Jupiter Trojans, the comets and small KBOs. Indeed, their spectral properties resemble that of anhydrous cometary dust. High-angular-resolution images of P-type asteroid (87) Sylvia with VLT/SPHERE were used to reconstruct its 3D shape, and to study the dynamics of its two satellites. We also model Sylvia's thermal evolution. The shape of Sylvia appears flattened and elongated. We derive a volume-equivalent diameter of 271 +/- 5 km, and a low density of 1378 +/- 45 kg.m-3. The two satellites orbit Sylvia on circular, equatorial orbits. The oblateness of Sylvia should imply a detectable nodal precession which contrasts with the fully-Keplerian dynamics of the satellites. This reveals an inhomogeneous internal structure, suggesting that Sylvia is differentiated. Sylvia's low density and differentiated interior can be explained by partial melting and mass redistribution through water percolation. The outer shell would be composed of material similar to interplanetary dust particles (IDPs) and the core similar to aqueously altered IDPs or carbonaceous chondrite meteorites such as the Tagish Lake meteorite. Numerical simulations of the thermal evolution of Sylvia show that for a body of such size, partial melting was unavoidable due to the decay of long-lived radionuclides. In addition, we show that bodies as small as 130-150 km in diameter should have followed a similar thermal evolution, while smaller objects, such as comets and the KBO Arrokoth, must have remained pristine, in agreement with in situ observations of these bodies. NASA Lucy mission target (617) Patroclus (diameter~140 km) may, however, be differentiated., Comment: Accepted for publication in A&A

The Euphrosyne asteroid family occupies a unique zone in orbital element space around 3.15 au and may be an important source of the low-albedo near-Earth objects. The parent body of this family may have been one of the planetesimals that delivered water and organic materials onto the growing terrestrial planets. We aim to characterize the compositional properties as well as the dynamical properties of the family. We performed a systematic study to characterize the physical properties of the Euphrosyne family members via low-resolution spectroscopy using the IRTF telescope. In addition, we performed smoothed-particle hydrodynamics (SPH) simulations and N-body simulations to investigate the collisional origin, determine a realistic velocity field, study the orbital evolution, and constrain the age of the Euphrosyne family. Our spectroscopy survey shows that the family members exhibit a tight taxonomic distribution, suggesting a homogeneous composition of the parent body. Our SPH simulations are consistent with the Euphrosyne family having formed via a reaccumulation process instead of a cratering event. Finally, our N-body simulations indicate that the age of the family is 280 Myr +180/-80 Myr, which is younger than a previous estimate., Comment: 10 pages, 13 figures, accepted to be published in A&A

Asteroid (31) Euphrosyne is one of the biggest objects in the asteroid main belt and the Euphrosyne family occupies a highly inclined region in the outer main belt and contains a remarkably large number of members, which is interpreted as an outcome of a disruptive cratering event. The goals of this adaptive-optics imaging study were threefold: to characterize the shape of Euphrosyne, to constrain its density, and to search for the large craters that may be associated with the family formation event. We obtained disk-resolved images of Euphrosyne using SPHERE/ZIMPOL at ESO's 8.2-m VLT as part of our large program (ID: 199.C-0074, PI: Vernazza). We reconstructed its 3D-shape using the adam shape modeling algorithm based on the SPHERE images and the available lightcurves of this asteroid. We analyzed the dynamics of the satellite with the genoid meta-heuristic algorithm. Finally, we studied the shape of Euphrosyne using hydrostatic equilibrium models. Our SPHERE observations show that Euphrosyne has a nearly spherical shape with the sphericity index of 0.9888 and its surface lacks large impact craters. Euphrosyne's diameter is 268+/-6 km, making it one of the top 10 largest main belt asteroids. We detected a satellite of Euphrosyne -- S/2019 (31) 1-- that is about 4 km across, on an circular orbit. The mass determined from the orbit of the satellite together with the volume computed from the shape model imply a density of 1665+/-242 kg/m^3, suggesting that Euphrosyne probably contain a large fraction of water ice in its interior. We find that the spherical shape of Euphrosyne is a result of the reaccumulation process following the impact, as in the case of (10) Hygiea. However, our shape analysis reveals that, contrary to Hygiea, the axis ratios of Euphrosyne significantly differ from the ones suggested by fluid hydrostatic equilibrium following reaccumulation., Comment: 12 pages, 6 figures, accepted to be published in A&A

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: Submitted to Icarus on 2019 January 7. The Electronic Supplementary Information to this paper is at http://www.asu.cas.cz/~asteroid/astpairscomplex_si.pdf

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.

Asteroid (16) Psyche is the target of the NASA Psyche mission. It is considered one of the few main-belt bodies that could be an exposed proto-planetary metallic core and that would thus be related to iron meteorites. Such an association is however challenged by both its near- and mid-infrared spectral properties and the reported estimates of its density. Here, we aim to refine the density of (16) Psyche to set further constraints on its bulk composition and determine its potential meteoritic analog. We observed (16) Psyche with ESO VLT/SPHERE/ZIMPOL as part of our large program (ID 199.C-0074). We used the high angular resolution of these observations to refine Psyche's three-dimensional (3D) shape model and subsequently its density when combined with the most recent mass estimates. In addition, we searched for potential companions around the asteroid. We derived a bulk density of 3.99\,$\pm$\,0.26\,g$\cdot$cm$^{-3}$ for Psyche. While such density is incompatible at the 3-sigma level with any iron meteorites ($\sim$7.8\,g$\cdot$cm$^{-3}$), it appears fully consistent with that of stony-iron meteorites such as mesosiderites (density $\sim$4.25\,$\cdot$cm$^{-3}$). In addition, we found no satellite in our images and set an upper limit on the diameter of any non-detected satellite of 1460\,$\pm$\,200}\,m at 150\,km from Psyche (0.2\%\,$\times$\,R$_{Hill}$, the Hill radius) and 800\,$\pm$\,200\,m at 2,000\,km (3\%\,$\times$\,$R_{Hill}$). Considering that the visible and near-infrared spectral properties of mesosiderites are similar to those of Psyche, there is merit to a long-published initial hypothesis that Psyche could be a plausible candidate parent body for mesosiderites., Comment: 16 pages

The rotation states of small asteroids are affected by a net torque arising from an anisotropic sunlight reflection and thermal radiation from the asteroids' surfaces. On long timescales, this so-called YORP effect can change asteroid spin directions and their rotation periods. We analyzed lightcurves of four selected near-Earth asteroids with the aim of detecting secular changes in their rotation rates that are caused by YORP. We use the lightcurve inversion method to model the observed lightcurves and include the change in the rotation rate $\mathrm{d} \omega / \mathrm{d} t$ as a free parameter of optimization. We collected more than 70 new lightcurves. For asteroids Toro and Cacus, we used thermal infrared data from the WISE spacecraft and estimated their size and thermal inertia. We also used the currently available optical and radar astrometry of Toro, Ra-Shalom, and Cacus to infer the Yarkovsky effect. We detected a YORP acceleration of $\mathrm{d}\omega / \mathrm{d} t = (1.9 \pm 0.3) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ for asteroid Cacus. For Toro, we have a tentative ($2\sigma$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}\omega / \mathrm{d} t = 3.0 \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$. For asteroid Eger, we confirmed the previously published YORP detection with more data and updated the YORP value to $(1.1 \pm 0.5) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. We also updated the shape model of asteroid Ra-Shalom and put an upper limit for the change of the rotation rate to $|\mathrm{d}\omega / \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3\sigma$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values.

Context. The so-called Barbarian asteroids share peculiar, but common polarimetric properties, probably related to both their shape and composition. They are named after (234) Barbara, the first on which such properties were identified. As has been suggested, large scale topographic features could play a role in the polarimetric response, if the shapes of Barbarians are particularly irregular and present a variety of scattering/incidence angles. This idea is supported by the shape of (234) Barbara, that appears to be deeply excavated by wide concave areas revealed by photometry and stellar occultations. Aims. With these motivations, we started an observation campaign to characterise the shape and rotation properties of Small Main- Belt Asteroid Spectroscopic Survey (SMASS) type L and Ld asteroids. As many of them show long rotation periods, we activated a worldwide network of observers to obtain a dense temporal coverage. Methods. We used light-curve inversion technique in order to determine the sidereal rotation periods of 15 asteroids and the con- vergence to a stable shape and pole coordinates for 8 of them. By using available data from occultations, we are able to scale some shapes to an absolute size. We also study the rotation periods of our sample looking for confirmation of the suspected abundance of asteroids with long rotation periods. Results. Our results show that the shape models of our sample do not seem to have peculiar properties with respect to asteroids with similar size, while an excess of slow rotators is most probably confirmed.

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., Comment: Accepted for publication in ICARUS Special Issue - Asteroids: Origin, Evolution & Characterization

Context. The high-angular-resolution capability of the new-generation ground-based adaptive-optics camera SPHERE at ESO VLT allows us to assess, for the very first time, the cratering record of medium-sized (D~100-200 km) asteroids from the ground, opening the prospect of a new era of investigation of the asteroid belt's collisional history. Aims. We investigate here the collisional history of asteroid (6) Hebe and challenge the idea that Hebe may be the parent body of ordinary H chondrites, the most common type of meteorites found on Earth (~34% of the falls). Methods. We observed Hebe with SPHERE as part of the science verification of the instrument. Combined with earlier adaptive-optics images and optical light curves, we model the spin and three-dimensional (3D) shape of Hebe and check the consistency of the derived model against available stellar occultations and thermal measurements. Results. Our 3D shape model fits the images with sub-pixel residuals and the light curves to 0.02 mag. The rotation period (7.274 47 h), spin (343 deg,+47 deg), and volume-equivalent diameter (193 +/- 6km) are consistent with previous determinations and thermophysical modeling. Hebe's inferred density is 3.48 +/- 0.64 g.cm-3 , in agreement with an intact interior based on its H-chondrite composition. Using the 3D shape model to derive the volume of the largest depression (likely impact crater), it appears that the latter is significantly smaller than the total volume of close-by S-type H-chondrite-like asteroid families. Conclusions. Our results imply that (6) Hebe is not the most likely source of H chondrites. Over the coming years, our team will collect similar high-precision shape measurements with VLT/SPHERE for ~40 asteroids covering the main compositional classes, thus providing an unprecedented dataset to investigate the origin and collisional evolution of the asteroid belt., Comment: 11 pages, 13 figures, accepted for publication in A&A

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., Comment: 16 pages, 21 figures, manuscript accepted in Astronomy & Astrophysics, to be published in section 10. Planets and planetary systems

The near-Earth asteroid (3200) Phaethon is an intriguing object: its perihelion is at only 0.14 au and is associated with the Geminid meteor stream. We aim to use all available disk-integrated optical data to derive a reliable convex shape model of Phaethon. By interpreting the available space- and ground-based thermal infrared data and Spitzer spectra using a thermophysical model, we also aim to further constrain its size, thermal inertia, and visible geometric albedo. We applied the convex inversion method to the new optical data obtained by six instruments and to previous observations. The convex shape model was then used as input for the thermophysical modeling. We also studied the long-term stability of Phaethon's orbit and spin axis with a numerical orbital and rotation-state integrator. We present a new convex shape model and rotational state of Phaethon: a sidereal rotation period of 3.603958(2) h and ecliptic coordinates of the preferred pole orientation of (319$^{\circ}$, $-$39$^{\circ}$) with a 5$^{\circ}$ uncertainty. Moreover, we derive its size ($D$=5.1$\pm$0.2 km), thermal inertia ($\Gamma$=600$\pm$200 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$), geometric visible albedo ($p_{\mathrm{V}}$=0.122$\pm$0.008), and estimate the macroscopic surface roughness. We also find that the Sun illumination at the perihelion passage during the past several thousand years is not connected to a specific area on the surface, which implies non-preferential heating., Comment: Astronomy and Astrophysics. In press

We use the lightcurve inversion method to derive new shape models and spin states of asteroids from the sparse-in-time photometry compiled in the Lowell Photometric Database. To speed up the time-consuming process of scanning the period parameter space through the use of convex shape models, we use the distributed computing project Asteroids@home, running on the Berkeley Open Infrastructure for Network Computing (BOINC) platform. This way, the period-search interval is divided into hundreds of smaller intervals. These intervals are scanned separately by different volunteers and then joined together. We also use an alternative, faster, approach when searching the best-fit period by using a model of triaxial ellipsoid. By this, we can independently confirm periods found with convex models and also find rotation periods for some of those asteroids for which the convex-model approach gives too many solutions. From the analysis of Lowell photometric data of the first 100,000 numbered asteroids, we derived 328 new models. This almost doubles the number of available models. We tested the reliability of our results by comparing models that were derived from purely Lowell data with those based on dense lightcurves, and we found that the rate of false-positive solutions is very low. We also present updated plots of the distribution of spin obliquities and pole ecliptic longitudes that confirm previous findings about a non-uniform distribution of spin axes. However, the models reconstructed from noisy sparse data are heavily biased towards more elongated bodies with high lightcurve amplitudes.

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., Comment: Accepted for publication in A&A

Asteroid (234) Barbara is the prototype of a category of asteroids that has been shown to be extremely rich in refractory inclusions, the oldest material ever found in the Solar System. It exhibits several peculiar features, most notably its polarimetric behavior. In recent years other objects sharing the same property (collectively known as "Barbarians") have been discovered. Interferometric observations in the mid-infrared with the ESO VLTI suggested that (234) Barbara might have a bi-lobated shape or even a large companion satellite. We use a large set of 57 optical lightcurves acquired between 1979 and 2014, together with the timings of two stellar occultations in 2009, to determine the rotation period, spin-vector coordinates, and 3-D shape of (234) Barbara, using two different shape reconstruction algorithms. By using the lightcurves combined to the results obtained from stellar occultations, we are able to show that the shape of (234) Barbara exhibits large concave areas. Possible links of the shape to the polarimetric properties and the object evolution are discussed. We also show that VLTI data can be modeled without the presence of a satellite., Comment: 10 pages, 6 figures

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)., Comment: Accepted for publication in A&A (September 16, 2013)

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., Comment: Accepted for publication in Icarus, July 16, 2013

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., Comment: Published in Planetary and Space Science, Volume 73, September 2012

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., Comment: Accepted for publication in A&A, January 15, 2013