
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 socalled YORP effect can change
asteroid spin directions and their rotation periods. We analyzed lightcurves of
four selected nearEarth 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, RaShalom, 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 RaShalom 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}$. RaShalom 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 socalled 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 lightcurve 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 diskintegrated optical data
(i.e., classical denseintime photometry obtained from public databases and
through a large collaboration network as well as sparseintime 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
diskresolved data obtained by the Near InfraRed Camera (Nirc2) mounted on the
W.M. Keck II telescope.

Context. The highangularresolution capability of the newgeneration
groundbased adaptiveoptics camera SPHERE at ESO VLT allows us to assess, for
the very first time, the cratering record of mediumsized (D~100200 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
adaptiveoptics images and optical light curves, we model the spin and
threedimensional (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 subpixel residuals and the light
curves to 0.02 mag. The rotation period (7.274 47 h), spin (343 deg,+47 deg),
and volumeequivalent diameter (193 +/ 6km) are consistent with previous
determinations and thermophysical modeling. Hebe's inferred density is 3.48 +/
0.64 g.cm3 , in agreement with an intact interior based on its Hchondrite
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 closeby Stype Hchondritelike 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
highprecision 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.

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 populationscale rather
than individualtarget observable in two studies so far. We aim to examine this
approach rigorously to establish its theoretical validity, degree of
illposedness, 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 Widefield Infrared Survey Explorer
(WISE) yield coarse but robust estimates of this distribution, and as an
example compare various asteroid families with each other.

The nearEarth 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 diskintegrated optical data to derive a reliable
convex shape model of Phaethon. By interpreting the available space and
groundbased 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 longterm stability of Phaethon's orbit and spin axis with a numerical
orbital and rotationstate 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 nonpreferential heating.

We use the lightcurve inversion method to derive new shape models and spin
states of asteroids from the sparseintime photometry compiled in the Lowell
Photometric Database. To speed up the timeconsuming 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
periodsearch 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
bestfit 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 convexmodel 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 falsepositive 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 nonuniform
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 mainbelt
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 diskresolved 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 diskintegrated optical data (i.e.,
classical denseintime photometry obtained from public databases and through a
large collaboration network as well as sparseintime 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 nearEarth
asteroids.

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 midinfrared with the ESO VLTI suggested
that (234) Barbara might have a bilobated 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, spinvector coordinates, and 3D 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.

Current amount of ~500 asteroid models derived from the diskintegrated
photometry by the lightcurve inversion method allows us to study not only the
spinvector 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 wellconstrained
poleecliptic 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 semimajor axis
drift and the consistency with the sizefrequency distribution of each family,
and consequently we remove interlopers. We then present the spinvector
distributions for eight asteroidal families. We use a combined orbital and
spinevolution model to explain the observed spinvector 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 poleecliptic 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 longterm
evolution of a collisional family is capable of reproducing well the observed
spinvector properties. Using this simulation, we also independently constrain
the age of families Flora (1.0\pm0.5 Gyr) and Koronis (2.54 Gyr).

We select 50 mainbelt 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 nearinfrared. 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 volumeequivalent 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 midIR photometric observations provided by IRAS, WISE
and AKARI. For 40 asteroids with previously determined mass estimates, we
compute their bulk densities and discuss the massdensity 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 mainbelt 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
sparseintime 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 sparseintime 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 polelatitude distributions, we were able to constrain the typical
values of the cYORP parameter as between 0.05 and 0.6.