We present a review of the problem of asteroid shape and spin reconstruction from generalised projections; i.e., from lightcurves, disk-resolved images, occultation silhouettes, radar range-Doppler data, and interferometry. The aim of this text is to summarize all important mathematical facts and proofs related to this inverse problem, to describe their implications to observers and modellers, and to provide the reader with all relevant references.

Time-resolved photometry of asteroids can be used for shape and spin reconstruction. If the number of measurements per asteroid is not sufficient to create a model, the whole data set can be used to reconstruct the distribution of shape elongations and pole latitudes in the population. This is done by reconstructing amplitudes of lightcurves that are estimated from dispersion of points observed at (assumed) constant aspect angle. Here, we formulate orbit-averaged approach where the observable is the orbit-averaged dispersion of brightness.

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 Large Synoptic Survey Telescope (LSST) is expected to increase known small solar system object populations by an order of magnitude or more over the next decade, enabling a broad array of transformative solar system science investigations to be performed. In this white paper, we discuss software tools and infrastructure that we anticipate will be needed to conduct these investigations and outline possible approaches for implementing them. Feedback from the community or contributions to future updates of this work are welcome. Our aim is for this white paper to encourage further consideration of the software development needs of the LSST solar system science community, and also to be a call to action for working to meet those needs in advance of the expected start of the survey in late 2022.

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).

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

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.

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.

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.

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.

In the past decade, hundreds of asteroid shape models have been derived using the lightcurve inversion method. At the same time, a new framework of 3-D shape modeling based on the combined analysis of widely different data sources such as optical lightcurves, disk-resolved images, stellar occultation timings, mid-infrared thermal radiometry, optical interferometry, and radar delay-Doppler data, has been developed. This multi-data approach allows the determination of most of the physical and surface properties of asteroids in a single, coherent inversion, with spectacular results. We review the main results of asteroid lightcurve inversion and also recent advances in multi-data modeling. We show that models based on remote sensing data were confirmed by spacecraft encounters with asteroids, and we discuss how the multiplication of highly detailed 3-D models will help to refine our general knowledge of the asteroid population. The physical and surface properties of asteroids, i.e., their spin, 3-D shape, density, thermal inertia, surface roughness, are among the least known of all asteroid properties. Apart for the albedo and diameter, we have access to the whole picture for only a few hundreds of asteroids. These quantities are nevertheless very important to understand as they affect the non-gravitational Yarkovsky effect responsible for meteorite delivery to Earth, or the bulk composition and internal structure of asteroids.

We introduce ADAM, the All-Data Asteroid Modelling algorithm. ADAM is simple and universal since it handles all disk-resolved data types (adaptive optics or other images, interferometry, and range-Doppler radar data) in a uniform manner via the 2D Fourier transform, enabling fast convergence in model optimization. The resolved data can be combined with disk-integrated data (photometry). In the reconstruction process, the difference between each data type is only a few code lines defining the particular generalized projection from 3D onto a 2D image plane. Occultation timings can be included as sparse silhouettes, and thermal infrared data are efficiently handled with an approximate algorithm that is sufficient in practice due to the dominance of the high-contrast (boundary) pixels over the low-contrast (interior) ones. This is of particular importance to the raw ALMA data that can be directly handled by ADAM without having to construct the standard image. We study the reliability of the inversion by using the independent shape supports of function series and control-point surfaces. When other data are lacking, one can carry out fast nonconvex lightcurve-only inversion, but any shape models resulting from it should only be taken as illustrative global-scale ones.

We present the analysis of high angular resolution observations of the triple Asteroid (87) Sylvia collected with three 8-10 m class telescopes (Keck, VLT, Gemini North) and the Hubble Space Telescope. The moons' mutual orbits were derived individually using a purely Keplerian model. We computed the position of Romulus, the outer moon of the system, at the epoch of a recent stellar occultation which was successfully observed at less than 15 km from our predicted position, within the uncertainty of our model. The occultation data revealed that the Moon, with a surface-area equivalent diameter Ds=23.1$\pm$0.7km, is strongly elongated (axes ratio of 2.7$\pm$0.32.7$\pm$0.3), significantly more than single asteroids of similar size in the main-belt. We concluded that its shape is probably affected by the tides from the primary. A new shape model of the primary was calculated combining adaptive-optics observations with this occultation and 40 archived light-curves recorded since 1978. The difference between the J2=0.024-0.009+0.016 derived from the 3-D shape model assuming an homogeneous distribution of mass for the volume equivalent diameter Dv=273$\pm$10km primary and the null J2 implied by the Keplerian orbits suggests a non-homogeneous mass distribution in the asteroid's interior.

Asteroids with satellites are natural laboratories to constrain the formation and evolution of our solar system. The binary Trojan asteroid (624) Hektor is the only known Trojan asteroid to possess a small satellite. Based on W.M. Keck adaptive optics observations, we found a unique and stable orbital solution, which is uncommon in comparison to the orbits of other large multiple asteroid systems studied so far. From lightcurve observations recorded since 1957, we showed that because the large Req=125-km primary may be made of two joint lobes, the moon could be ejecta of the low-velocity encounter, which formed the system. The inferred density of Hektor's system is comparable to the L5 Trojan doublet (617) Patroclus but due to their difference in physical properties and in reflectance spectra, both captured Trojan asteroids could have a different composition and origin.

With the adaptive optics (AO) system on the 10 m Keck-II telescope, we acquired a high quality set of 84 images at 14 epochs of asteroid (52) Europa on 2005 January 20. The epochs covered its rotation period and, by following its changing shape and orientation on the plane of sky, we obtained its triaxial ellipsoid dimensions and spin pole location. An independent determination from images at three epochs obtained in 2007 is in good agreement with these results. By combining these two data sets, along with a single epoch data set obtained in 2003, we have derived a global fit for (52) Europa of diameters (379x330x249) +/- (16x8x10) km, yielding a volume-equivalent spherical-diameter of 315 +/- 7 km, and a rotational pole within 7 deg of [RA; Dec] = [257,+12] in an Equatorial J2000 reference frame (ECJ2000: 255,+35). Using the average of all mass determinations available forEuropa, we derive a density of 1.5 +/- 0.4, typical of C-type asteroids. Comparing our images with the shape model of Michalowski et al. (A&A 416, 2004), derived from optical lightcurves, illustrates excellent agreement, although several edge features visible in the images are not rendered by the model. We therefore derived a complete 3-D description of Europa's shape using the KOALA algorithm by combining our imaging epochs with 4 stellar occultations and 49 lightcurves. We use this 3-D shape model to assess these departures from ellipsoidal shape. Flat facets (possible giant craters) appear to be less distinct on (52) Europa than on other C-types that have been imaged in detail. We show that fewer giant craters, or smaller craters, is consistent with its expected impact history. Overall, asteroid (52) Europa is still well modeled as a smooth triaxial ellipsoid with dimensions constrained by observations obtained over several apparitions.

The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time d\omega /dt. Values of d\omega/ dt derived from observations were compared with the values predicted by theory. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger d\omega / dt = (1.4 +/- 0.6) x 10^{-8} rad/d (3\sigma error), which corresponds to a decrease in the rotation period by 4.2 ms/yr. The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |d\omega / dt| < 8 x 10^{-9} rad/d for Cerberus, and |d\omega / dt| < 3 x 10^{-8} rad/d for Ra-Shalom.

We present a comparison of our results from ground-based observations of asteroid (21) Lutetia with imaging data acquired during the flyby of the asteroid by the ESA Rosetta mission. This flyby provided a unique opportunity to evaluate and calibrate our method of determination of size, 3-D shape, and spin of an asteroid from ground-based observations. We present our 3-D shape-modeling technique KOALA which is based on multi-dataset inversion. We compare the results we obtained with KOALA, prior to the flyby, on asteroid (21) Lutetia with the high-spatial resolution images of the asteroid taken with the OSIRIS camera on-board the ESA Rosetta spacecraft, during its encounter with Lutetia. The spin axis determined with KOALA was found to be accurate to within two degrees, while the KOALA diameter determinations were within 2% of the Rosetta-derived values. The 3-D shape of the KOALA model is also confirmed by the spectacular visual agreement between both 3-D shape models (KOALA pre- and OSIRIS post-flyby). We found a typical deviation of only 2 km at local scales between the profiles from KOALA predictions and OSIRIS images, resulting in a volume uncertainty provided by KOALA better than 10%. Radiometric techniques for the interpretation of thermal infrared data also benefit greatly from the KOALA shape model: the absolute size and geometric albedo can be derived with high accuracy, and thermal properties, for example the thermal inertia, can be determined unambiguously. We consider this to be a validation of the KOALA method. Because space exploration will remain limited to only a few objects, KOALA stands as a powerful technique to study a much larger set of small bodies using Earth-based observations.

Main belt asteroids (6070) Rheinland and (54827) 2001NQ8 belong to a small population of couples of bodies which reside on very similar heliocentric orbits. Vokrouhlicky & Nesvorny (2008, AJ 136, 280) promoted a term "asteroid pairs", pointing out their common origin within the past tens to hundreds of ky. Previous attempts to reconstruct the initial configuration of Rheinland and 2001NQ8 at the time of their separation have led to the prediction that Rheinland's rotation should be retrograde. Here we report extensive photometric observations of this asteroid and use the lightcurve inversion technique to directly determine its rotation state and shape. We confirm the retrograde sense of rotation of Rheinland, with obliquity value constrained to be >= 140 deg. The ecliptic longitude of the pole position is not well constrained as yet. The asymmetric behavior of Rheinland's lightcurve reflects a sharp, near-planar edge in our convex shape representation of this asteroid. Our calibrated observations in the red filter also allow us to determine $H_R = 13.68\pm 0.05$ and $G = 0.31\pm 0.05$ values of the H-G system. With the characteristic color index $V-R = 0.49\pm 0.05$ for the S-type asteroids, we thus obtain $H = 14.17\pm 0.07$ for the absolute magnitude of (6070) Rheinland. This a significantly larger value than previously obtained from analysis of the astrometric survey observations. We next use the obliquity constraint for Rheinland to eliminate some degree of uncertainty in the past propagation of its orbit. This is because the sign of the past secular change of its semimajor axis due to the Yarkovsky effect is now constrained. Determination of the rotation state of the secondary component, asteroid (54827) 2001NQ8, is the key element in further constraining the age of the pair and its formation process.

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.

Near-Earth asteroid 162173 (1999 JU3) is a potential flyby and rendezvous target for interplanetary missions because of its easy to reach orbit. The physical and thermal properties of the asteroid are relevant for establishing the scientific mission goals and also important in the context of near-Earth object studies in general. Our goal was to derive key physical parameters such as shape, spin-vector, size, geometric albedo, and surface properties of 162173 (1999 JU3). With three sets of published thermal observations (ground-based N-band, Akari IRC, Spitzer IRS), we applied a thermophysical model to derive the radiometric properties of the asteroid. The calculations were performed for the full range of possible shape and spin-vector solutions derived from the available sample of visual lightcurve observations. The near-Earth asteroid 162173 (1999 JU3) has an effective diameter of 0.87 +/- 0.03 km and a geometric albedo of 0.070 +/- 0.006. The chi2-test reveals a strong preference for a retrograde sense of rotation with a spin-axis orientation of lambda_ecl = 73 deg, beta_ecl = -62 deg and P_sid = 7.63 +/- 0.01 h. The most likely thermal inertia ranges between 200 and 600 Jm-2s-0.5K-1, about a factor of 2 lower than the value for 25143 Itokawa. This indicates that the surface lies somewhere between a thick-dust regolith and a rock/boulder/cm-sized, gravel-dominated surface like that of 25143 Itokawa. Our analysis represents the first time that shape and spin-vector information has been derived from a combined data set of visual lightcurves (reflected light) and mid-infrared photometry and spectroscopy (thermal emission).

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present the results of our Thousand Asteroid Light Curve Survey (TALCS) conducted with the Canada-France-Hawaii Telescope in September 2006. Our untargeted survey detected 828 Main Belt asteroids to a limiting magnitude of g'~22.5 corresponding to a diameter range of 0.4 km <= D <= 10 km. Of these, 278 objects had photometry of sufficient quality to perform rotation period fits. We debiased the observations and light curve fitting process to determine the true distribution of rotation periods and light curve amplitudes of Main Belt asteroids. We confirm a previously reported excess in the fraction of fast rotators but find a much larger excess of slow rotating asteroids (~15% of our sample). A few percent of objects in the TALCS size range have large light curve amplitudes of ~1 mag. Fits to the debiased distribution of light curve amplitudes indicate that the distribution of triaxial ellipsoid asteroid shapes is proportional to the square of the axis-ratio, (b/a)^2, and may be bi-modal. Finally, we find six objects with rotation periods that may be less than 2 hours with diameters between 400 m and 1.5 km, well above the break-up limit for a gravitationally-bound aggregate. Our debiased data indicate that this population represents <4% of the Main Belt in the 1-10 km size range.

We report on the results of a six-month photometric study of the main-belt binary C-type asteroid 121 Hermione, performed during its 2007 opposition. We took advantage of the rare observational opportunity afforded by one of the annual equinoxes of Hermione occurring close to its opposition in June 2007. The equinox provides an edge-on aspect for an Earth-based observer, which is well suited to a thorough study of Hermione's physical characteristics. The catalog of observations carried out with small telescopes is presented in this work, together with new adaptive optics (AO) imaging obtained between 2005 and 2008 with the Yepun 8-m VLT telescope and the 10-m Keck telescope. The most striking result is confirmation that Hermione is a bifurcated and elongated body, as suggested by Marchis et al., (2005). A new effective diameter of 187 +/- 6 km was calculated from the combination of AO, photometric and thermal observations. The new diameter is some 10% smaller than the hitherto accepted radiometric diameter based on IRAS data. The reason for the discrepancy is that IRAS viewed the system almost pole-on. New thermal observations with the Spitzer Space Telescope agree with the diameter derived from AO and lightcurve observations. On the basis of the new AO astrometric observations of the small 32-km diameter satellite we have refined the orbit solution and derived a new value of the bulk density of Hermione of 1.4 +0.5/-0.2 g cm-3. We infer a macroscopic porosity of ~33 +5/-20%.