Context. The study of non-principal axis (NPA) rotators can provide important clues to the evolution of the spin state of asteroids. However, so far, very few studies have focused on NPA-rotating main-belt asteroids (MBAs). One of MBAs that are known to be in an excited rotation state is asteroid (5247) Krylov. Aims. By using disk-integrated photometric data, we construct a physical model of (5247) Krylov including shape and spin state. Methods. We apply the light curve convex inversion method employing optical light curves obtained by using ground-based telescopes in three apparitions during 2006, 2016, and 2017, along with infrared light curves obtained by the Wide-field Infrared Survey Explorer (WISE) satellite in 2010. Results. Asteroid (5247)~Krylov is spinning in a short axis mode (SAM) characterized by rotation and precession periods of 368.7 hr and 67.27 hr, respectively. The angular momentum vector orientation of Krylov is found to be $\lambda_{L} = 298^\circ$ and $\beta_{L} = -58^\circ$. The ratio of the rotational kinetic energy to the basic spin state energy $E/E_{0} \simeq 1.02$ shows that the (5247) Krylov is about 2% excited state compared to the Principal Axis (PA) rotation state. The shape of (5247) Krylov can be approximated by an elongated prolate ellipsoid with a ratio of moments of inertia of $I_{a}:I_{b}:I_{c}=0.36:0.96:1$. This is the first physical model of NPA rotator among MBAs. The physical processes that led to the current NPA rotation cannot be unambiguously reconstructed.

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.

The nature and physical properties of asteroids, in particular those orbiting in the near-Earth space, are of scientific interest and practical importance. Exoplanet surveys can be excellent resources to detect asteroids, both already known and new objects. This is due their similar observing requirements: large fields of view, long sequences, and short cadence. If the targeted fields are not located far from the ecliptic, many asteroids will cross occasionally the field of view. We present two complementary methodologies to identify asteroids serendipitously observed in large-area astronomical surveys. One methodology focuses on detecting already known asteroids using the Virtual Observatory tool SkyBoT, which predicts their positions and motions in the sky at a specific epoch. The other methodology applies the ssos pipeline, which is able to identify known and new asteroids based on their apparent motion. The application of these methods to the 6.4 deg 2 of the sky covered by the Wide-Field CAMera Transit Survey in the J-band is described. We identified 15 661 positions of 1 821 different asteroids. Of them, 182 are potential new discoveries. A publicly accessible online, Virtual Observatory compliant catalogue was created. We obtained the shapes and periods for five of our asteroids from their light-curves built with additional photometry taken from external archives. We demonstrated that our methodologies are robust and reliable approaches to find, at zero cost of observing time, asteroids observed by chance in astronomical surveys. Our future goal is to apply them to other surveys with adequate temporal coverage.

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.

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.

Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models. Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $\beta$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data). Methods. We use the LEADER algorithm to compare the $b/a$ and $\beta$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $\beta$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method. Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h.

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.

The JAXA Hayabusa-2 mission was approved in 2010 and launched on December 3, 2014. The spacecraft will arrive at the near-Earth asteroid 162173 Ryugu in 2018 where it will perform a survey, land and obtain surface material, then depart in Dec 2019 and return to Earth in Dec 2020. We observed Ryugu with the Herschel Space Observatory in Apr 2012 at far-IR thermal wavelengths, supported by several ground-based observations to obtain optical lightcurves. We reanalysed previously published Subaru-COMICS and AKARI-IRC observations and merged them with a Spitzer-IRS data set. In addition, we used a large set of Spitzer-IRAC observations obtained in the period Jan to May, 2013. The data set includes two complete rotational lightcurves and a series of ten "point-and-shoot" observations. The almost spherical shape of the target together with the insufficient lightcurve quality forced us to combine radiometric and lightcurve inversion techniques in different ways to find the object's key physical and thermal parameters. We find that the solution which best matches our data sets leads to this C class asteroid having a retrograde rotation with a spin-axis orientation of (lambda = 310-340 deg; beta = -40+/-15 deg) in ecliptic coordinates, an effective diameter (of an equal-volume sphere) of 850 to 880 m, a geometric albedo of 0.044 to 0.050 and a thermal inertia in the range 150 to 300 Jm-2s-0.5K-1. Based on estimated thermal conductivities of the top-layer surface in the range 0.1 to 0.6 WK-1m-1, we calculated that the grain sizes are approximately equal to between 1 and 10 mm. The finely constrained values for this asteroid serve as a `design reference model', which is currently used for various planning, operational and modelling purposes by the Hayabusa2 team.

Context: Large all-sky surveys provide us with a lot of photometric data that are sparse in time (typically few measurements per night) and can be potentially used for the determination of shapes and rotational states of asteroids. The method which is generally used to derive these parameters is the lightcurve inversion. However, for most asteroids their sparse data are not accurate enough to derive a unique model and the lightcurve inversion method is thus not very efficient. Aims: To fully utilize photometry sparse in time, we developed a new simplified model and applied it on the data from the Lowell photometric database. Our aim was to derive spin axis orientations and shape elongations of asteroids and to reconstruct distributions of these parameters for selected subpopulations to find if there are some differences. Methods: We model asteroids as geometrically scattering triaxial ellipsoids. Observed values of mean brightness and the dispersion of brightness are compared with computed values obtained from the parameters of the model -- ecliptical longitude $\lambda$ and latitude $\beta$ of the pole and the ratios $a/b$, $b/c$ of axes of the ellipsoid. These parameters are optimized to get the best agreement with the observation. Results: We found that the distribution of $\lambda$ for main-belt asteroids is not uniform (in agreement with findings of Bowell et al., 2014, M&PS, 49, 95) and is dependent on the inclination of orbit. Surprisingly, the non-uniformity of $\lambda$ distribution is larger for asteroids residing on low-inclination orbits. We also studied distributions of $a/b$ for several groups of asteroids and found that small asteroids ($D<25\,$km) are on average more elongated than large ones.

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.

Infrared radiation emitted from an asteroid surface causes a torque that can significantly affect rotational state of the asteroid. The influence of small topographic features on this phenomenon, called the YORP effect, seems to be of utmost importance. In this work, we show that a lateral heat diffusion in boulders of suitable sizes leads to an emergence of a local YORP effect which magnitude is comparable to the YORP effect due to the global shape. We solve a three-dimensional heat diffusion equation in a boulder and its surroundings by the finite element method, using the FreeFem++ code. The contribution to the total torque is inferred from the computed temperature distribution. Our general approach allows us to compute the torque induced by a realistic irregular boulder. For an idealized boulder, our result is consistent with an existing one-dimensional model. We also estimated (and extrapolated) a size distribution of boulders on (25143) Itokawa from close-up images of its surface. We realized that topographic features on Itokawa can potentially induce a torque corresponding to a rotational acceleration of the order 10^-7 rad day^-2 and can therefore explain the observed phase shift in light curves.

Using our photometric observations taken between 1996 and 2013 and other published data, we derived properties of the binary near-Earth asteroid (175706) 1996 FG3 including new measurements constraining evolution of the mutual orbit with potential consequences for the entire binary asteroid population. We also refined previously determined values of parameters of both components, making 1996 FG3 one of the most well understood binary asteroid systems. We determined the orbital vector with a substantially greater accuracy than before and we also placed constraints on a stability of the orbit. Specifically, the ecliptic longitude and latitude of the orbital pole are 266{\deg} and -83{\deg}, respectively, with the mean radius of the uncertainty area of 4{\deg}, and the orbital period is 16.1508 +/- 0.0002 h (all quoted uncertainties correspond to 3sigma). We looked for a quadratic drift of the mean anomaly of the satellite and obtained a value of 0.04 +/- 0.20 deg/yr^2, i.e., consistent with zero. The drift is substantially lower than predicted by the pure binary YORP (BYORP) theory of McMahon and Scheeres (McMahon, J., Scheeres, D. [2010]. Icarus 209, 494-509) and it is consistent with the theory of an equilibrium between BYORP and tidal torques for synchronous binary asteroids as proposed by Jacobson and Scheeres (Jacobson, S.A., Scheeres, D. [2011]. ApJ Letters, 736, L19). Based on the assumption of equilibrium, we derived a ratio of the quality factor and tidal Love number of Q/k = 2.4 x 10^5 uncertain by a factor of five. We also derived a product of the rigidity and quality factor of mu Q = 1.3 x 10^7 Pa using the theory that assumes an elastic response of the asteroid material to the tidal forces. This very low value indicates that the primary of 1996 FG3 is a 'rubble pile', and it also calls for a re-thinking of the tidal energy dissipation in close asteroid binary systems.

From a set of adaptive optics (AO) observations collected with the W.M. Keck telescope between August and September 2009, we derived the orbital parameters of the most recently discovered satellites of the large C-type asteroid (93) Minerva. The satellites of Minerva, which are approximately 3 and 4 km in diameter, orbit very close to the primary $\sim$5 & $\sim$8 $\times$ Rp and $\sim$1% & $\sim$2% $\times$ RHill) in a circular manner, sharing common characteristics with most of the triple asteroid systems in the main-belt. Combining these AO observations with lightcurve data collected since 1980 and two stellar occultations in 2010 & 2011, we removed the ambiguity of the pole solution of Minerva's primary and showed that it has an almost regular shape with an equivalent diameter Deq = 154 $\pm$ 6 km in agreement with IRAS observations. The surprisingly high bulk density of 1.75 $\pm$ 0.30 g/cm$\^3$ for this C-type asteroid, suggests that this taxonomic class is composed of asteroids with different compositions, For instance, Minerva could be made of the same material as dry CR, CO, and CV meteorites. We discuss possible scenarios on the origin of the system and conclude that future observations may shine light on the nature and composition of this fifth known triple main-belt asteroid.

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.