Asteroid (3200) Phaethon is an active near-Earth asteroid and the parent body of the Geminid Meteor Shower. Because of its small perihelion distance, Phaethon's surface reaches temperatures sufficient to destabilize hydrated materials. We conducted rotationally resolved spectroscopic observations of this asteroid, mostly covering the northern hemisphere and the equatorial region, beyond 2.5-micron to search for evidence of hydration on its surface. Here we show that the observed part of Phaethon does not exhibit the 3-micron hydrated mineral absorption (within 2-sigma). These observations suggest that Phaethon's modern activity is not due to volatile sublimation or devolatilization of phyllosilicates on its surface. It is possible that the observed part of Phaethon was originally hydrated and has since lost volatiles from its surface via dehydration, supporting its connection to the Pallas family, or it was formed from anhydrous material.

Rotation properties (spin-axis direction and rotation period) and coarse shape models of asteroids can be reconstructed from their disk-integrated brightness when measured from various viewing geometries. These physical properties are essential for creating a global picture of structure and dynamical evolution of the main belt. The number of shape and spin models can be increased not only when new data are available, but also by combining independent data sets and inverting them together. Our aim was to derive new asteroid models by processing readily available photometry. We used asteroid photometry compiled in the Lowell Observatory photometry database with photometry from the Gaia Data Release 2. Both data sources are available for about 5400 asteroids. In the framework of the Asteroids@home distributed computing project, we applied the light curve inversion method to each asteroid to find its convex shape model and spin state that fits the observed photometry. Due to the limited number of Gaia DR2 data points and poor photometric accuracy of Lowell data, we were able to derive unique models for only ~1100 asteroids. Nevertheless, 762 of these are new models that significantly enlarge the current database of about 1600 asteroid models. Our results demonstrate the importance of a combined approach to inversion of asteroid photometry. While our models in general agree with those obtained by separate inversion of Lowell and Gaia data, the combined inversion is more robust, model parameters are more constrained, and unique models can be reconstructed in many cases when individual data sets alone are not sufficient.

The recent close approach of the NEA (3200) Phaethon offered a rare opportunity to obtain high-quality observational data. We used the newly obtained optical light curves to improve the spin and shape model of Phaethon and to determine its surface physical properties derived by thermophysical modeling. We also used the available astrometric observations of Phaethon, including those obtained by the Arecibo radar and the Gaia spacecraft, to constrain the secular drift of the orbital semimajor axis. This constraint allowed us to estimate the bulk density by assuming that the drift is dominated by the Yarkovsky effect. We used the convex inversion model to derive the 3D shape model of Phaethon, and a detailed numerical approach for an accurate analysis of the Yarkovsky effect. We obtained a unique solution for Phaethon's pole orientation at $(318,-47)^{\circ}$ ecliptic longitude and latitude (uncertainty of $5^{\circ}$), and confirm the previously reported thermophysical properties ($D=5.1\pm0.2$ km, $\Gamma=600\pm200$ SI). Phaethon has a top-like shape with possible north-south asymmetry. The characteristic size of the regolith grains is 1-2 cm. The orbit analysis reveals a secular drift of the semimajor axis of $-(6.9\pm1.9)\times 10^{-4}$ au Myr$^{-1}$. With the derived volume-equivalent size of 5.1~km, the bulk density $\rho$ is $1.67\pm0.47$ g cm$^{-3}$. If the size is slightly larger $\sim5.7$ km, as suggested by radar data, $\rho$ would decrease to $1.48\pm0.42$ g cm$^{-3}$. We further investigated the suggestion that Phaethon may be in a cluster with asteroids (155140) 2005 UD and (225416) 1999 YC that was formed by rotational fission of a critically spinning parent body. Phaethon's $\rho$ is consistent with typical values for large ($>100$ km) C-complex asteroids and supports its association with asteroid (2) Pallas. These findings render a cometary origin unlikely for Phaethon.

In addition to stellar data, Gaia Data Release 2 (DR2) also contains accurate astrometry and photometry of about 14,000 asteroids covering 22 months of observations. We used Gaia asteroid photometry to reconstruct rotation periods, spin axis directions, and the coarse shapes of a subset of asteroids with enough observations. One of our aims was to test the reliability of the models with respect to the number of data points and to check the consistency of these models with independent data. Another aim was to produce new asteroid models to enlarge the sample of asteroids with known spin and shape. We used the lightcurve inversion method to scan the period and pole parameter space to create final shape models that best reproduce the observed data. To search for the sidereal rotation period, we also used a simpler model of a geometrically scattering triaxial ellipsoid. By processing about 5400 asteroids with at least ten observations in DR2, we derived models for 173 asteroids, 129 of which are new. Models of the remaining asteroids were already known from the inversion of independent data, and we used them for verification and error estimation. We also compared the formally best rotation periods based on Gaia data with those derived from dense lightcurves. We show that a correct rotation period can be determined even when the number of observations $N$ is less than 20, but the rate of false solutions is high. For $N > 30$, the solution of the inverse problem is often successful and the parameters are likely to be correct in most cases. These results are very promising because the final Gaia catalogue should contain photometry for hundreds of thousands of asteroids, typically with several tens of data points per object, which should be sufficient for reliable spin reconstruction.

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.

Information about the spin state of asteroids is important for our understanding of the dynamical processes affecting them. However, spin properties of asteroids are known for only a small fraction of the whole population. To enlarge the sample of asteroids with a known rotation state and basic shape properties, we combined sparse-in-time photometry from the Lowell Observatory Database with flux measurements from NASA's WISE satellite. We applied the light curve inversion method to the combined data. The thermal infrared data from WISE were treated as reflected light because the shapes of thermal and visual light curves are similar enough for our purposes. While sparse data cover a wide range of geometries over many years, WISE data typically cover an interval of tens of hours, which is comparable to the typical rotation period of asteroids. The search for best-fitting models was done in the framework of the Asteroids@home distributed computing project. By processing the data for almost 75,000 asteroids, we derived unique shape models for about 900 of them. Some of them were already available in the DAMIT database and served us as a consistency check of our approach. In total, we derived new models for 662 asteroids, which significantly increased the total number of asteroids for which their rotation state and shape are known. For another 789 asteroids, we were able to determine their sidereal rotation period and estimate the ecliptic latitude of the spin axis direction. We studied the distribution of spins in the asteroid population. We revealed a significant discrepancy between the number of prograde and retrograde rotators for asteroids smaller than about 10 km. Combining optical photometry with thermal infrared light curves is an efficient approach to obtaining new physical models of asteroids.

By means of a varied-shape thermophysical model (VS-TPM) of Hanus et al. (2015) that takes into account asteroid shape and pole uncertainties, we analyze the thermal IR data acquired by the NASA's WISE satellite of about 300 asteroids with derived convex shape models. We utilize publicly available convex shape models and rotation states as input for the TPM. For more than one hundred asteroids, the TPM gives us an acceptable fit to the thermal IR data allowing us to report their size, thermal inertia, surface roughness or visible geometric albedo. This work more than doubles the number of asteroids with determined thermophysical properties. In the remaining cases, the shape model and pole orientation uncertainties, specific rotation or thermophysical properties, poor thermal IR data or their coverage prevent the determination of reliable thermophysical properties. Finally, we present the main results of the statistical study of derived thermophysical parameters within the whole population of main-belt asteroids and within few asteroid families. Our sizes are, in average, consistent with the radiometric sizes reported by Mainzer et al. (2016). The thermal inertia increases with decreasing size, but a large range of thermal inertia values is observed within the similar size ranges between D~10-100 km. We derived unexpectedly low thermal inertias (<20 SI) for several asteroids with sizes 10<D<50 km, indicating a very fine and mature regolith on their surface. The thermal inertia values seem to be consistent within several collisional families. The fast rotators with P<4 h tend to have slightly larger thermal inertia values, so probably are not covered by a fine regolith. This could be explained, for example, by the loss of the fine regolith due to the centrifugal force, or by the ineffectiveness of the regolith production (e.g., by the thermal cracking mechanism of Delbo' et al. 2014).

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.

We present shape models and volume estimates of twenty asteroids based on relative photometry and adaptive optics images. We discuss error estimation and the effects of myopic deconvolution on shape solutions. For further analysis of the information capacities of data sources, we also present and discuss ambiguity and uniqueness results for the reconstruction of nonconvex shapes from photometry.

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.

We present physical models, i.e. convex shapes, directions of the rotation axis, and sidereal rotation periods, of 18 asteroids out of which 10 are new models and 8 are refined models based on much larger data sets than in previous work. The models were reconstructed by the lightcurve inversion method from archived publicly available lightcurves and our new observations with BlueEye600 robotic observatory. One of the new results is the shape model of asteroid (1663)~van~den~Bos with the rotation period of 749\,hr, which makes it the slowest rotator with known shape. We describe our strategy for target selection that aims at fast production of new models using the enormous potential of already available photometry stored in public databases. We also briefly describe the control software and scheduler of the robotic observatory and we discuss the importance of building a database of asteroid models for studying asteroid physical properties in collisional families.

Optical lightcurves can be used for the shape and spin reconstruction of asteroids. Due to unknown albedo, these models are scale-free. When thermal infrared data are available, they can be used for scaling the shape models and for deriving thermophysical properties of the surface by applying a thermophysical model. We introduce a new method of simultaneous inversion of optical and thermal infrared data that allows the size of an asteroid to be derived along with its shape and spin state. The method optimizes all relevant parameters (shape and its size, spin state, light-scattering properties, thermal inertia, surface roughness) by gradient-based optimization. The thermal emission is computed by solving the 1-D heat diffusion equation. Calibrated optical photometry and thermal fluxes at different wavelengths are needed as input data. We demonstrate the reliability and test the accuracy of the method on selected targets with different amount and quality of data. Our results in general agree with those obtained by independent methods. Combining optical and thermal data into one inversion method opens a new possibility for processing photometry from large optical sky surveys with the data from WISE. It also provides more realistic estimates of errors of thermophysical parameters.

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.

Disk-integrated photometric data of asteroids do not contain accurate information on shape details or size scale. Additional data such as disk-resolved images or stellar occultation measurements further constrain asteroid shapes and allow size estimates. We aim to use all available disk-resolved images of about forty asteroids obtained by the Near-InfraRed Camera (Nirc2) mounted on the W.M. Keck II telescope together with the disk-integrated photometry and stellar occultation measurements to determine their volumes. We can then use the volume, in combination with the known mass, to derive the bulk density. We download and process all asteroid disk-resolved images obtained by the Nirc2 that are available in the Keck Observatory Archive (KOA). We combine optical disk-integrated data and stellar occultation profiles with the disk-resolved images and use the All-Data Asteroid Modeling (ADAM) algorithm for the shape and size modeling. Our approach provides constraints on the expected uncertainty in the volume and size as well. We present shape models and volume for 41 asteroids. For 35 asteroids, the knowledge of their mass estimates from the literature allowed us to derive their bulk densities. We clearly see a trend of lower bulk densities for primitive objects (C-complex) than for S-complex asteroids. The range of densities in the X-complex is large, suggesting various compositions. Moreover, we identified a few objects with rather peculiar bulk densities, which is likely a hint of their poor mass estimates. Asteroid masses determined from the Gaia astrometric observations should further refine most of the density estimates.

Asteroid (130) Elektra belongs to one of the six known triple asteroids in the main belt, so its mass has been reliably determined. We aim to use all available disk-resolved images of (130) Elektra obtained by the SPHERE instrument at VLT and by the Nirc2 of the Keck telescope together with the disk-integrated photometry to determine its shape model and its size. The volume can be then used in combination with the known mass to derive the bulk density of the primary. We apply the All-Data Asteroid Modeling (ADAM) algorithm to the optical disk-integrated data, 2 disk-resolved images obtained by the SPHERE instrument and 13 disk-resolved images from the Nirc2 of the Keck telescope, and derive the shape model and size of Elektra. We present the shape model, volume-equivalent diameter (199$\pm$7 km) and bulk density (1.60$\pm$0.13 g cm$^{-3}$) of the C-type asteroid Elektra.

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.

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.

We present the project Asteroids@home that uses distributed computing to solve the time-consuming inverse problem of shape reconstruction of asteroids. The project uses the Berkeley Open Infrastructure for Network Computing (BOINC) framework to distribute, collect, and validate small computational units that are solved independently at individual computers of volunteers connected to the project. Shapes, rotational periods, and orientations of the spin axes of asteroids are reconstructed from their disk-integrated photometry by the lightcurve inversion method.

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.

Asteroid sizes can be directly measured by observing occultations of stars by asteroids. When there are enough observations across the path of the shadow, the asteroid's projected silhouette can be reconstructed. Asteroid shape models derived from photometry by the lightcurve inversion method enable us to predict the orientation of an asteroid for the time of occultation. By scaling the shape model to fit the occultation chords, we can determine the asteroid size with a relative accuracy of typically ~ 10%. We combine shape and spin state models of 44 asteroids (14 of them are new or updated models) with the available occultation data to derive asteroid effective diameters. In many cases, occultations allow us to reject one of two possible pole solutions that were derived from photometry. We show that by combining results obtained from lightcurve inversion with occultation timings, we can obtain unique physical models of asteroids.

Tens of thousands of sparse-in-time lightcurves from astrometric projects are publicly available. We investigate these data and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population. We use sparse photometry from selected observatories from the AstDyS database, either alone or in combination with dense lightcurves, to determine new asteroid models by the lightcurve inversion method. We present 80 new asteroid models derived from combined data sets where sparse photometry is taken from the AstDyS database and dense lightcurves are from the Uppsala Asteroid Photometric Catalogue (UAPC) and from several individual observers. For 18 asteroids, we present updated shape solutions based on new photometric data. The addition of new models increases the total number of models derived by the lightcurve inversion method to ~200. We also present a simple statistical analysis of physical properties of asteroids where we look for possible correlations between various physical parameters with an emphasis on the spin vector. We present the observed and de-biased distributions of ecliptic latitudes with respect to different size ranges of asteroids as well as a simple theoretical model of the latitude distribution and then compare its predictions with the observed distributions. From this analysis we find that the latitude distribution of small asteroids (D < 30 km) is clustered towards ecliptic poles and can be explained by the YORP thermal effect while the latitude distribution of larger asteroids (D > 60 km) exhibits an evident excess of prograde rotators, probably of primordial origin.