Issue 
A&A
Volume 525, January 2011



Article Number  A145  
Number of page(s)  6  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201015599  
Published online  09 December 2010 
Research Note
Thermophysical properties of 162173 (1999 JU3), a potential flyby and rendezvous target for interplanetary missions
^{1}
MaxPlanckInstitut für extraterrestrische Physik,
Giessenbachstraße, 85748
Garching, Germany
email: tmueller@mpe.mpg.de
^{2}
Astronomical Institute, Faculty of Mathematics and Physics,
Charles University in Prague, V
Holešovičkách 2, 180 00, Praha 8, Czech
Republic
^{3}
Institute of Space and Astronautical Science,
Japan Aerospace Exploration Agency, 311
Yoshinodai, Sagamihara, Kanagawa
2298510,
Japan
^{4}
Department of Physics and Astronomy, The University of Western
Ontario, 1151
Richimond St. London, ON, N6A
3K7, Canada
^{5}
Institute of Astronomy, National Central University,
300 Jhongda RD.,
Jhongli, Taoyuan
32001,
Taiwan
^{6}
Okayama Astrophysical Observatory, National Astronomical
Observatory of Japan, 30375 Honjo,
Kamogata, Asaguchi, Okayama
7190232,
Japan
^{7}
Bisei Spaceguard Center, Japan Spaceguard Association, 17163
Okura, Bisei, Ibara, Okayama
7141411,
Japan
^{8}
Kiso Observatory, Institute of Astronomy, The University of
Tokyo, 1076230 Mitake,
Kiso, Nagano
3970101,
Japan
^{9}
Tokyo Metropolitan Government, 281 Nishishinjyuku, Shinjyuku,
Tokyo
1638001,
Japan
^{10}
Aichi Toho University, 311 Heiwagaoka, Meito, Nagoya, Aichi
4688515,
Japan
^{11}
Jet Propulsion Laboratory, 4800 Oak Grove Drive, MS 183301, Pasadena, CA
91109,
USA
^{12}
Korea Astronomy and Space Science Institute,
611 Hwaamdong, Yueinggu,
Daejeon, Taejon
305348,
Korea
^{13}
Lunar and Planetary Laboratory, University of Arizona,
Tucson,
AZ, 857210092, USA
^{14}
National Astronomical Observatory of Japan,
2211 Osawa,
Mitaka, Tokyo
1818588,
Japan
Received:
17
August
2010
Accepted:
5
November
2010
Context. NearEarth asteroid 162173 (1999 JU3) is a potential flyby and rendezvous target for interplanetary missions because of its easytoreach orbit. The physical and thermal properties of the asteroid are relevant for establishing the scientific mission goals and also important in the context of nearEarth object studies in general.
Aims. Our goal was to derive key physical parameters such as shape, spinvector, size, geometric albedo, and surface properties of 162173 (1999 JU3).
Methods. With three sets of published thermal observations (groundbased Nband, 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 spinvector solutions derived from the available sample of visual lightcurve observations.
Results. The nearEarth asteroid 162173 (1999 JU3) has an effective diameter of 0.87 ± 0.03 km and a geometric albedo of 0.070 ± 0.006. The χ^{2}test reveals a strong preference for a retrograde sense of rotation with a spinaxis orientation of λ_{ecl} = 73°, β_{ecl} = −62° and P_{sid} = 7.63 ± 0.01 h. The most likely thermal inertia ranges between 200 and 600 J m^{2} s^{0.5} K^{1}, about a factor of 2 lower than the value for 25143 Itokawa. This indicates that the surface lies somewhere between a thickdust regolith and a rock/boulder/cmsized, graveldominated surface like that of 25143 Itokawa. Our analysis represents the first time that shape and spinvector information has been derived from a combined data set of visual lightcurves (reflected light) and midinfrared photometry and spectroscopy (thermal emission).
Key words: minor planets, asteroids: individual: 162173 (1999 JU3) / radiation mechanisms: thermal / techniques: photometric / infrared: planetary systems
© ESO, 2010
1. Introduction
Asteroid 162173 (1999 JU3) is currently among the potential targets of future interplanetary exploration missions. The target is relatively easy to reach with stateoftheart mission capabilities, and it offers high scientific potential (Binzel et al. 2004). This nearEarth asteroid belongs to the Cclass objects, which are believed to represent primitive, volatilerich remnants of the early solar system. Various aspects of this small body have been covered in some detail in the recent works by Hasegawa et al. (2008) and by Campins et al. (2009).
Hasegawa et al. (2008) use a spherical shape model for their radiometric analysis, and alternatively use an ellipsoidal shape model, but without knowing the true spinvector orientation. The results (radiometric diameter of 0.92 ± 0.12 km, visual geometric albedo of 0.063$\begin{array}{}\mathrm{+}\mathrm{0.020}\\ 0.015\end{array}$) which indicate a thermal inertia larger than 500 J m^{2} s^{0.5} K^{1}), were based on a set of photometric Subaru and Akari observations and connected to simplified shape and spinaxis assumptions. Campins et al. (2009) have obtained a singleepoch Spitzer infrared spectrograph (IRS) spectrum. They used a spherical shape model, and for the spinpole orientation they used the extreme case of an equatorial retrograde geometry and a prograde solution published by Abe et al. (2008). Their analysis, based on the single IRSspectrum and ignoring the data sets published by Hasegawa et al. (2008), yielded a value for the thermal inertia of 700 ± 200 J m^{2} s^{0.5} K^{1}, a diameter estimate of 0.90 ± 0.14 km, and geometric albedo of 0.07 ± 0.01.
Despite the simplifications in shape and spinvector properties, both sets of published radiometric diameter and albedo values agree within the given uncertainties. Both teams also favoured a relatively high thermal inertia (close to that of 25143 Itokawa). The unknowns of the spinvector orientation cause a large uncertainty in the thermal properties.
Here we reanalyse all available lightcurve observations to derive (on the basis of standard χ^{2} lightcurve inversion techniques) all matching spinvector and shape solutions (Sect. 2). The full possible range for shape, spinaxis orientation, and rotation period was then used as input for a thermophysical χ^{2} analysis of all available thermal observations (Sect. 3) with the goal of deriving radiometric sizes, albedos and thermal inertias. At the same time, we determined the most likely shapesolution, rotation period and spinaxis orientation (Sect. 4).
2. Possible shape and spinvector solutions
A detailed list of the available photometric observations is presented in Table 1. There are about 40 dedicated visual lightcurve data sets spread over more than 270 days. About half of the lightcurves were calibrated; some of the data sets are very noisy. Based on these data, Abe et al. (2008) found a rotation period of 7.6272 ± 0.0072 h and a spin orientation of λ_{ecl} = 331.0°, β_{ecl} = + 20.0°, indicating a prograde rotation. A more recent analysis by the same authors (priv. communication) resulted in a rotation model with slightly different values: λ_{ecl} = 327.3°, β_{ecl} = + 34.7°, P_{sid} = 7.6273922 h. Both solutions were derived using the epoch and amplitude methods described by Magnusson (1986). These methods are reliable for irregularly shaped bodies and for sufficient lightcurve data covering various aspect angles. 1999 JU3 has a comparatively spherical shape and the available lightcurve data were from restricted directions. Nevertheless, both solutions gave a reasonable match to the observed lightcurves, but it turned out that these solutions are not unique and other parameter sets with different rotation periods, spinaxis orientations and shape models could not be excluded (see Fig. 1).
Observation circumstances for the lightcurve measurements.
Fig. 1 Match between observed lightcurves and shape/spinaxis model solution “7_1”. 
Since these values are crucial input parameters for our thermophysical model analysis, we repeated the search for possible shape and spinvector solutions using the lightcurve inversion method developed by Kaasalainen & Torppa (2001). Our goal was to derive a set of the most likely convex shape models that would fit all available lightcurves. Because of the poor quality of some of the data, we were only able to determine the range of the sidereal rotation period to 7.6204–7.6510 h. Different shape models with pole orientations covering almost the entire celestial sphere (without any preference for pro or retrograde solutions) fit the data equally well. By scanning the periodpole parameter space, we derived 77 shape models that were physically acceptable (rotating around the shortest axis) and for which the χ^{2} of the fit was no more than 10% higher than the bestfit χ^{2}. These models corresponded to local minima in the parameter space.
In addition, we included both of the original Abe et al. (2008) solutions and added another five with the pole fixed to the two (λ_{ecl}, β_{ecl})pairs mentioned above and rotation periods in the range given by Abe et al. (2008). Only two of these models pass the χ^{2} + 10% limit on the basis of the visual lightcurves, the other five models have higher χ^{2}values. Six of these shape models are elongated along the “z” axis and therefore unphysical.
For all 84 models we performed a thermophysical model analysis for a wide range of possible parameters (see Table 3).
3. Thermophysical model analysis
3.1. Model and input parameters
The midIR photometric data were already described in Hasegawa et al. (2008). The data set includes 15 Nband Subaru observations and two dedicated Akari observations at 15 and 24 μm. We binned the singleepoch Spitzer IRS data (Campins et al. 2009) into 20 wavelength points (4 for band SL2, 7 for SL1, 4 for LL2, 5 for LL1; see Fig. 5, bottom). The 20 wavelength points were chosen to give about equal weight to the two published data samples in terms of number of observations (17 in Hasegawa et al. (2008) and 20 for the Campins et al. (2009) sample). In this way the derived object properties are better connected to the entire data set, and they do not just match the measurements of one observing epoch. Each observation set also has a mixture of lower and higher quality data: the groundbased Subaru data are of lower quality than the Akari data, while the quality of the IRS spectrum changes with wavelength. This again ensures that the final solutions are not biased towards a single measurement. All observations are listed in Table 2.
Summary of the avaible thermal observations.
All 84 possible spinvector and shape solutions from Sect. 2 have been used in combination with these thermal data.
For our analysis we are using a thermophysical model (TPM) described by Lagerros (1996, 1997, 1998) and Müller & Lagerros (1998). This TPM works with true illumination and observing geometries, accepts irregular shape models and arbitrary spinvector solutions, works with roughness controlled by the rms of the surface slopes, considers heatconduction into the surface as well as multiple scattering of both the solar and the thermally emitted radiation. The model has been tested and validated thoroughly for NEAs (e.g., Müller et al. 2005) and MBAs (e.g., Müller & Lagerros 2002). The TPM input parameters and applied variations are listed in Table 3.
Summary of general TPM input parameters and applied ranges.
3.2. Solving for effective diameter, geometric albedo and thermal inertia
Campins et al. (2009) and Mueller (2007) used a χ^{2} or reduced χ^{2}test to find solutions for the thermal inertia Γ. Here we follow a modified approach to find the most robust solutions with respect to thermal inertia and allowing for the full range in effective diameter and geometric albedo at the same time. The following procedure was executed for all 84 possible shape and spinvector solutions separately:
 (1)
For each value of Γ in a wide range (see Table 3) we calculate the radiometric diameter and albedo solution via the TPM for each individually observed thermal flux (37 individual diameter and albedo solutions). Diameter and albedo are linked by the absolute magnitude H_{V} which was kept constant (the rotational amplitude is only about 0.1 mag).
 (2)
We calculate the weighted mean radiometric diameter and albedo solution for each given Γ ($\mathit{x\u0305}\mathrm{=}\frac{\mathrm{\Sigma}{\mathit{x}}_{\mathit{i}}\mathit{/}{\mathit{\sigma}}_{\mathit{i}}^{\mathrm{2}}}{\mathrm{\Sigma}\mathrm{1}\mathit{/}{\mathit{\sigma}}_{\mathit{i}}^{\mathrm{2}}}$, with diameter/albedo errors σ_{i} connected to the observational errors).
 (3)
For each individual observation we predict TPM fluxes based on the given Γ and the corresponding weighted mean radiometric diameter and albedo from step (2).
 (4)
The most robust solutions occur when the observations and the TPM predictions agree best (taking the uncertainties of the measurements into account in a weighted mean sense, see step 2). This can be expressed as $\frac{\mathrm{1}}{\mathit{N}}{\mathrm{\Sigma}}_{\mathit{i}\mathrm{=}\mathrm{1}}^{\mathit{N}}\mathrm{\left(}\mathrm{\right(}{\mathrm{obs}}_{\mathit{i}}\mathrm{}{\mathrm{mod}}_{\mathit{i}}\mathrm{)}\mathit{/}{\mathit{\sigma}}_{\mathit{i}}{\mathrm{)}}^{\mathrm{2}}$, a modified reduced χ^{2} method. The most likely thermal inertia is found at the smallest χ^{2} values; the connected effective diameter and geometric albedo values are the ones calculated at step (2).
In a first round of executing this procedure we kept the surface roughness constant at values which were specified as “default values” for large, regolithcovered asteroids (Müller et al. 1999). The corresponding roughness parameter values are ρ = 0.7, the rms of the surface slopes, and f = 0.6, the fraction of surface covered by craters. The results are shown in Fig. 2.
Fig. 2 TPM χ^{2}optimization process to find robust solutions for diameter, albedo and thermal inertia simultaneously. Each line represents the reduced χ^{2} values for an individual shape/spinvector solution as a function of thermal inertia. The surface roughness was kept constant using the “default values” of ρ = 0.7 and f = 0.6 as proposed by Müller et al. (1999). Model 7_1 solutions marked with squares on the solid line correspond to the two cases with thermal inertias of 20 and 1000 J m^{2} s^{0.5} K^{1} shown in Fig. 3 (left and middle). 
The best shape and spinvector solutions (with lowest values for the reduced χ^{2} and clear minima in Fig. 2) were then the starting point for further tests:

(i)
Are these solutions robust against subsets of the thermal data?
 (ii)
How does the surface roughness influence the solutions?
 (iii)
Do the solutions explain the thermal behaviour over the observed phase angle range (from ~20° to ~55°)?
 (iv)
Is the TPM match of equal quality at all observed wavelengths?
 (v)
Are there large discrepancies at certain rotational phases?
4. Results and discussion
4.1. Solution for shape and spinvector
The shape and spinvector solutions which produce the lowest χ^{2}values in Fig. 2 are listed in Table 4. The model IDs represent a full shapemodel, each with more than 2000 surface elements and more than 1000 vertices. The Julian date at zero rotational phase γ_{0} is in all cases T_{0} = 2454289.0.
All three solutions are retrograde solutions and, in fact, the eight best solutions in Fig. 2 are retrograde solutions. The best prograde solutions in the χ^{2} picture are models with ID 5_2 and ID 13_5, both have χ^{2}minima at 2.3 (almost a factor of 4 higher than the best retrograde solution 7_1) and would require an extremely high thermal inertia (>1000 J m^{2} s^{0.5} K^{1}) to match the observations. Our most likely solution can be summarized as: λ_{ecl} = 73°, β_{ecl} = −62°, P_{sid} = 7.63 ± 0.01 h. Figure 3 (left, middle) shows the modelID 7_1, as seen from Spitzer during the IRSobservations and for the two thermal inertias marked in Fig. 2 with boxes. The match between model “7_1” with observed lightcurves is shown in Fig. 1.
We also analysed the thermal data set against the Abe et al. (2008) spinpole solutions discussed above. The corresponding χ^{2}minima are more than a factor five worse than our two best models above (see dashed lines in Fig. 2). Some of these solutions are compatiple with the inertia range given by Campins et al. (2009) and even produced an excellent match to the IRSspectrum. Nevertheless, these solutions can be excluded with very high confidence:

(i)
the match to the rest of the data set (Akari and groundbaseddata) is very poor (reflected in the high χ^{2}values);
 (ii)
the corresponding shape models are unphysical with elongations along the spin axis. Such rotational states would not be stable. The best of these models (in terms of χ^{2}minima) is shown in Fig. 3 on the right side (with the rotation along the zaxis with the object’s largest extension).
4.2. Thermal inertia and surface roughness
The thermal inertia is clearly a key parameter when modelling the midIR data for NEAs; it strongly influences the shape of the spectral energy distribution (SED). This can be seen in the Spitzer IRS spectrum, especially in the Wienpart of the spectrum. But the thermal inertia also drives the thermal behaviour as a function of phase angle (e.g., Müller 2002), the thermal phase curve. The temperature of the unilluminated fraction of the surface changes strongly with thermal inertia. In Fig. 3 the unilluminated fraction has just rotated out of the solar insolation. In one case (left) we assumed a low thermal inertia of ~20 J m^{2} s^{0.5} K^{1}, which was found to be typical for large MBAs (Müller & Lagerros 2002), while in the second case (middle) the thermal inertia is 1000 J m^{2} s^{0.5} K^{1}, close to the value found by Müller et al. (2005) for 25143 Itokawa.
Fig. 3 Left and middle: TPM implementation of the shape model 7_1 at the time of the Spitzer IRS observations on 2008May02 02:01 UT, as seen from Spitzer in asteroidcoordinates, i.e., the zaxis goes along the rotation axis. Left: low thermal inertia. Middle: high thermal inertia. Both solutions are marked with squares in Fig. 2. Right: shape model fixed on the Abe et al. (2008) spinvector and tuned to match the Campins et al. (2009) findings. The rotation along the largest object extension (zaxis) is unphysical. The surface temperatures are given in Kelvin. 
The importance of the thermal inertia in the modelling is also visible in the χ^{2}solutions (Figs. 2 and 4): the χ^{2}values change significantly when going through the whole grid of physically meaningful thermal inertias.
But Fig. 4 demonstrates that thermal inertia and surface roughness are not easy to disentangle, at least on the basis of this data set. Both parameters influence the shortwavelength SEDpart where the hottest surface temperatures dominate the SEDshape. One way of solving for both parameters would be to obtain data with a larger phase angle coverage and larger wavelength coverage. The roughness plays a much bigger role at small phase angles (beaming effect), but thermal data close to opposition are not available for 162173 (1999 JU3). The thermal inertia is more important at large phase angles and at longer wavelengths where the diskaveraged temperature dominates the SED shape. In general, the larger the coverage in phase angle and wavelength, the more accurate is the determination of diameter, albedo, thermal inertia and roughness.
Figure 4 also shows that similarly low χ^{2}values can be reached independent of the surface roughness. This demonstrates that roughness effects are still important for the interpretation of thermal data (via the thermal phase curves), even at these relatively large phase angles (Müller 2002). But the data set does not allow us to constrain the surface roughness which broadens the possible range of thermal inertias as can be seen in Fig. 4. Based on our three best shape/spinvector solutions and considering the uncertainties due to roughness, we conclude that the thermal inertia is in the range 200–600 J m^{2} s^{0.5} K^{1}. This range is in agreement with the lower limit of 150 J m^{2} s^{0.5} K^{1} given by Campins et al. (2009) but lower than their best fit value of 700 ± 200 J m^{2} s^{0.5} K^{1}. Our value is about a factor of two lower than the one found for 25143 Itokawa (Müller et al. 2005). We expect that the surface of 162173 (1999 JU3) is therefore different in the sense that there might be fewer rocks and boulders and that the surface includes millimetre sized particles (as opposed to the cmsized gravel on Itokawa). We also find a rigorous lower limit to the thermal inertia of about 100 J m^{2} s^{0.5} K^{1}, similar to Campins et al. (2009). This limit is not compatible with a thick dusty regolith covering the entire surface which would result in a very low thermal conductivity and thermal inertia values well below 100 J m^{2} s^{0.5} K^{1}, which would in principle be possible for an object of that size and the relativly slow rotation rate. For comparison, the Moon’s thick regolith gives a value below 40 J m^{2} s^{0.5} K^{1} (Keihm 1984, calculated for T = 300 K).
4.3. Radiometric diameter and albedo solution
Our best fit to all observations, as represented in Fig. 2 by the solid line, resulted in a radiometric effective diameter of 0.87 ± 0.02 km and 0.070 ± 0.003 for the visual geometric albedo. Both values are within the error bars of the solutions found by Hasegawa et al. (2008) and by Campins et al. (2009), but now with much smaller errors. The uncertainties are based on the best χ^{2}values for modelIDs “7_1” and “8_2” and the full variation in roughness (as shown in Fig. 4). The radiometric effective diameter is connected to the most likely shape model and spinvector solution from above and corresponds to the size of a spherical object with equal volume. Using the determined possible thermal inertia range of 200–600 J m^{2} s^{0.5} K^{1} and uncertainties in roughness and H_{V}magnitude lead to a possible value for the effective diameter of 0.87 ± 0.03 km. The derived geometric albedo provides the best solution between all thermal observations over the full wavelength and phase angle range and the absolute H_{V}magnitude. The H_{V}magnitude was given with only 0.02 mag error. Using a more realistic H_{V}error of ± 0.1 mag leads to a final geometric albedo value of 0.070 ± 0.006. The small uncertainties reflect the importance of multiepoch, multiwavelength and large phase angle coverage for thermophysical studies of small bodies. In a similar study by Müller et al. (2005), also based on a large thermal data set and a shape model from lightcurve inversion techniques, the derived effective diameter agreed within 2% of the true, insitu diameter (Müller et al. 2005; Fujiwara et al. 2006). The quoted uncertainties above are formal errors from the χ^{2} optimization, including the possible range in thermal inertia, roughness and H_{V}.
Fig. 4 TPM χ^{2}optimization process for the shape and spinaxis solutions with the lowest χ^{2}values (models 7_1, 8_2, 14_8). The solid lines are the “default roughness” values derived for MBAs. 
The remaining discrepancy between the Spitzer and the Akariflux at 15 μm (seen in Fig. 5 bottom) might be related to an additional error introduced by the flux scaling done by Campins et al. (2009) to match the shortwavelength part of the spectrum (<14 μm) to the longwavelength part of the spectrum (>14 μm). The mismatch is caused mainly by the placement of the object within the IRS slit and reflected in the specified 10% systematic absolute calibration uncertainty.
Although the Fig. 4 solutions and the model match in Fig. 5 look very convincing, there are some uncertainties remaining: The spinvector and shape solution from lightcurve inversion techniques is not very robust; more lightcurve observations are needed to improve the quality. The χ^{2}test works best if the thermal observations cover a wide range of wavelengths, phase angles (before and after opposition) and rotational phases. But all thermal observations have been taken at preopposition (positive) phase angles (leading the Sun) where the unilluminated part of the surface visible to the observer is either warm (prograde rotation) or cold (retrograde rotation); observations after opposition are not yet available. A combined before and after opposition data set would constrain the sense of rotation and the thermal properties much better, observations close to opposition would determine the surface roughness a bit better, hence constrain the thermal inertia further.
It is also important to note here that higher thermal inertias (>600 J m^{2} s^{0.5} K^{1}, Fig. 3 middle) would make a slightly better fit to the Spitzer IRS spectrum (improvement mainly at the shortest wavelengths in Fig. 5 bottom), but would cause a significant dependency in the diameter and albedo solutions with phase angle. The measurements taken at around 20° phase angle (Subaru) would then have fluxes that are about 40% higher than the corresponding model predictions (i.e., values >1.4 in Fig. 5, top). The Akari fluxes would still be ~20% higher than the model predictions. Our observational data set covering about 30° in phase angle constrains the possible thermal inertia to values below 700 J m^{2} s^{0.5} K^{1}.
Fig. 5 All observations divided by the corresponding TPM predictions based on our optimized radiometric solution. Top: as a function of phase angle. Bottom: as a function of wavelength. The full set of Spitzer IRS data are shown with little circles (Campins et al. 2009), the triangles are the rebinned data, the Akari IRC data are represented by squares; the SubaruCOMICS observations by diamond symbols. 
5. Conclusions
The radiometric analysis provides the following results: (i) a strong indication of a retrograde sense of rotation; (ii) a spinvector with λ_{ecl} = 73°, β_{ecl} = −62°, P_{sid} = 7.63 ± 0.01 h and γ_{0} = 0 at T_{0} = 2454289.0, derived for the first time based on a combined analysis of visual lightcurve data and thermal observations; (iii) a shape model (here labelled with 7_1) as shown in Fig. 3 (left & middle); (iv) a thermal inertia in the range 200 to 600 J m^{2} s^{0.5} K^{1}; (v) a radiometric effective diameter (of an equal volume sphere) of D_{eff} = 0.87 ± 0.03 km; (vi) a radiometric geometric albedo of p_{V} = 0.070 ± 0.006; (vii) a lower thermal inertia than for Itokawa, suggesting the presence of smaller particles, <cmsized, in the regolith, though likely not fine dust; (viii) very good agreement in the radiometric solutions between the Spitzer, the Akari and the Subaru observations; (ix) an excellent match of the flux changes with phase angle (the phase angle range covered here is from ~20° to ~55°).
The example of 162173 (1999 JU3) shows that a combination of visual lightcurves (reflected sunlight) and mid/farIR photometry or photospectroscopy (thermal emission) can improve the quality of shape and spinvector solutions significantly.
Acknowledgments
J. D. received grants from the Czech Science Foundation (GACR P209/10/0537) and the Research Program MSM0021620860 of the Ministry of education. S. H. was supported by the Space Plasma Laboratory, ISAS, JAXA. We are also grateful to Professor N. Kawai and the gammaray bursts project members for furnishing their optical camera at Ishigakijima Astronomical Observatory. Development of the optical CCD camera at Ishigakijima Astronomical Observatory was supported by the Ministry of Education, Science, Sports and Culture, GrantinAid for Creative Scientific Research. T.K. thanks to the JSPS Research Fellowships for Research Abroad for their financial support. This work was supported in part by the NASA Planetary Astronomy Program and performed in part at the Jet Propulsion Laboratory. We thank the Steward Observatory of the University of Arizona for its allocation of telescope time. We also thank the referee Dr. J. Emery for very helpful comments.
References
 Abe, M., Kawakami, K., Hasegawa, S., et al. 2008, COSPAR Scientific Assembly, B04006108 [Google Scholar]
 Binzel, R. P., Perozzi, E., Rivkin, A. S., Rossi, A., et al. 2004, Meteorit. Planet. Sci., 39, 351 [NASA ADS] [CrossRef] [Google Scholar]
 Campins, H., Emery, J. P., Kelley, M., et al. 2009, A&A, 503, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fujiwara, A., Kawaguchi, J., Uesugi, K., et al. 2006, Science, 312, 1330 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Hasegawa, S., Müller, T. G., Kawakami, K., et al. 2008, PASJ, 60, 399 [Google Scholar]
 Kaasalainen, M., & Torppa, J. 2001, Icarus, 153, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Keihm, S. J. 1984, Icarus, 60, 568 [NASA ADS] [CrossRef] [Google Scholar]
 Lagerros, J. S. V. 1996, A&A, 310, 1011 [NASA ADS] [Google Scholar]
 Lagerros, J. S. V. 1997, A&A, 325, 1226 [NASA ADS] [Google Scholar]
 Lagerros, J. S. V. 1998, A&A, 332, 1123 [NASA ADS] [Google Scholar]
 Magnusson, P. 1986, Icarus, 68, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Müller, T. G., & Lagerros, J. S. V. 1998, A&A, 338, 340 [NASA ADS] [Google Scholar]
 Müller, T. G., Lagerros, J. S. V., Burgdorf, M., et al. 1999, in The Universe as Seen by ISO, ed. P. Cox, & M. F. Kessler, ESA SP427, 141 [Google Scholar]
 Müller, T. G., & Lagerros, J. S. V. 2002, A&A, 381, 324 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Müller, T. G. 2002, M&PS, 37, 1919 [Google Scholar]
 Müller, T. G., Sekiguchi, T., Kaasalainen, M., Abe, M., & Hasegawa, S. 2005, A&A, 443, 347 (M05) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mueller, M. 2007, DLR, Thesis FU Berlin (http://www.diss.fuberlin.de/diss/receive/FUDISS_thesis_000000002596) [Google Scholar]
All Tables
All Figures
Fig. 1 Match between observed lightcurves and shape/spinaxis model solution “7_1”. 

In the text 
Fig. 2 TPM χ^{2}optimization process to find robust solutions for diameter, albedo and thermal inertia simultaneously. Each line represents the reduced χ^{2} values for an individual shape/spinvector solution as a function of thermal inertia. The surface roughness was kept constant using the “default values” of ρ = 0.7 and f = 0.6 as proposed by Müller et al. (1999). Model 7_1 solutions marked with squares on the solid line correspond to the two cases with thermal inertias of 20 and 1000 J m^{2} s^{0.5} K^{1} shown in Fig. 3 (left and middle). 

In the text 
Fig. 3 Left and middle: TPM implementation of the shape model 7_1 at the time of the Spitzer IRS observations on 2008May02 02:01 UT, as seen from Spitzer in asteroidcoordinates, i.e., the zaxis goes along the rotation axis. Left: low thermal inertia. Middle: high thermal inertia. Both solutions are marked with squares in Fig. 2. Right: shape model fixed on the Abe et al. (2008) spinvector and tuned to match the Campins et al. (2009) findings. The rotation along the largest object extension (zaxis) is unphysical. The surface temperatures are given in Kelvin. 

In the text 
Fig. 4 TPM χ^{2}optimization process for the shape and spinaxis solutions with the lowest χ^{2}values (models 7_1, 8_2, 14_8). The solid lines are the “default roughness” values derived for MBAs. 

In the text 
Fig. 5 All observations divided by the corresponding TPM predictions based on our optimized radiometric solution. Top: as a function of phase angle. Bottom: as a function of wavelength. The full set of Spitzer IRS data are shown with little circles (Campins et al. 2009), the triangles are the rebinned data, the Akari IRC data are represented by squares; the SubaruCOMICS observations by diamond symbols. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.