Search

Results from the TESS mission showed that previous studies strngly underestimated the number of slow rotators, revealing the importance of studying those asteroids. For most slowly rotating asteroids (P > 12), no spin and shape model is available because of observation selection effects. This hampers determination of their thermal parameters and accurate sizes. We continue our campaign in minimising selection effects among main belt asteroids. Our targets are slow rotators with low light-curve amplitudes. The goal is to provide their scaled spin and shape models together with thermal inertia, albedo, and surface roughness to complete the statistics. Rich multi-apparition datasets of dense light curves are supplemented with data from Kepler and TESS. In addition to data in the visible range, we also use thermal data from infrared space observatories (IRAS, Akari and WISE) in a combined optimisation process using the Convex Inversion Thermophysical Model (CITPM). This novel method has so far been applied to only a few targets, and in this work we further validate the method. We present the models of 16 slow rotators. All provide good fits to both thermal and visible data. The obtained sizes are on average accurate at the 5% precision, with diameters in the range from 25 to 145 km. The rotation periods of our targets range from 11 to 59 hours, and the thermal inertia covers a wide range of values, from 2 to <400 SI units, not showing any correlation with the period. With this work we increase the sample of slow rotators with reliable spin and shape models and known thermal inertia by 40%. The thermal inertia values of our sample do not display a previously suggested increasing trend with rotation period, which might be due to their small skin depth., Comment: Accepted to Astronomy & Astrophysics. 10 pages + appendices

Thermal-infrared measurements of asteroids are crucial for deriving the objects' sizes, albedos, and also the thermophysical properties of the surface material. Depending on the available data, a range of simple to complex thermal models are applied to achieve specific science goals. However, testing these models is often a difficult process and the uncertainties of the derived parameters are not easy to estimate. Here, we make an attempt to verify a widely accepted thermophysical model (TPM) against unique thermal infrared (IR), full-disk, and well-calibrated measurements of the Moon. The data were obtained by the High-resolution InfraRed Sounder (HIRS) instruments on board a fleet of Earth weather satellites that serendipitously scan over the Moon. We found 22 Moon intrusions, taken in 19 channels between 3.75 micron and 15.0 micron, and over a wide phase angle range from -73.1 deg to +73.8 deg. The similarity between these Moon data and typical asteroid spectral-IR energy distributions allows us to benchmark the TPM concepts and to point out problematic aspects. The TPM predictions match the HIRS measurements within 5% (10% at the shortest wavelengths below 5 micron when using the Moon's known properties (size, shape, spin, albedo, thermal inertia, roughness) in combination with a newly established wavelength-dependent hemispherical emissivity. In the 5-7.5 micron and in the 9.5 to 11 micron ranges, the global emissivity model deviates considerably from the known lunar sample spectra. Our findings will influence radiometric studies of near-Earth and main-belt asteroids in cases where only short-wavelength data (from e.g., NEOWISE, the warm Spitzer mission, or ground-based M-band measurements) are available. The new, full-disk IR Moon model will also be used for the calibration of IR instrumentation on interplanetary missions (e.g., for Hayabusa-2) and weather satellites., Comment: 21 pages, 9 figures, 7 tables, accepted for publication in Astronomy & Astrophysics in March 2021

We predicted a stellar occultation of the bright star Gaia DR1 4332852996360346368 (UCAC4 385-75921) (m$_{\rm V}$= 14.0 mag) by the centaur 2002 GZ$_{32}$ for 2017 May 20$^{\rm th}$. Our latest shadow path prediction was favourable to a large region in Europe. Observations were arranged in a broad region inside the nominal shadow path. Series of images were obtained with 29 telescopes throughout Europe and from six of them (five in Spain and one in Greece) we detected the occultation. This is the fourth centaur, besides Chariklo, Chiron and Bienor, for which a multi-chord stellar occultation is reported. By means of an elliptical fit to the occultation chords we obtained the limb of 2002 GZ$_{32}$ during the occultation, resulting in an ellipse with axes of 305 $\pm$ 17 km $\times$ 146 $\pm$ 8 km. From this limb, thanks to a rotational light curve obtained shortly after the occultation, we derived the geometric albedo of 2002 GZ$_{32}$ ($p_{\rm V}$ = 0.043 $\pm$ 0.007) and a 3-D ellipsoidal shape with axes 366 km $\times$ 306 km $\times$ 120 km. This shape is not fully consistent with a homogeneous body in hydrostatic equilibrium for the known rotation period of 2002 GZ$_{32}$. The size (albedo) obtained from the occultation is respectively smaller (greater) than that derived from the radiometric technique but compatible within error bars. No rings or debris around 2002 GZ$_{32}$ were detected from the occultation, but narrow and thin rings cannot be discarded., Comment: Accepted for publication in MNRAS (8-Dec.-2020), 15 pages, 9 figures

Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.

On 28th January 2018, the large Trans-Neptunian Object (TNO) 2002TC302 occulted a m$_v= $15.3 star with ID 130957813463146112 in the Gaia DR2 stellar catalog. 12 positive occultation chords were obtained from Italy, France, Slovenia and Switzerland. Also, 4 negative detections were obtained near the north and south limbs. This represents the best observed stellar occultation by a TNO other than Pluto, in terms of the number of chords published thus far. From the 12 chords, an accurate elliptical fit to the instantaneous projection of the body, compatible with the near misses, can be obtained. The resulting ellipse has major and minor axes of 543 $\pm$ 18 km and 460 $\pm$ 11 km, respectively, with a position angle of 3 $\pm$ 1 degrees for the minor axis. This information, combined with rotational light curves obtained with the 1.5m telescope at Sierra Nevada Observatory and the 1.23m telescope at Calar Alto observatory, allows us to derive possible 3D shapes and density estimations for the body, based on hydrostatic equilibrium assumptions. The effective area equivalent diameter is $\sim$ 84 km smaller than the radiometrically derived diameter using thermal data from Herschel and Spitzer Space Telescopes. This might indicate the existence of an unresolved satellite of up to $\sim$ 300 km in diameter, to account for all the thermal flux, although the occultation and thermal diameters are compatible within their error bars given the considerable uncertainty of the thermal results. The existence of a potential satellite also appears to be consistent with other ground-based data presented here. From the effective occultation diameter combined with H$_V$ measurements we derive a geometric albedo of 0.147 $\pm$ 0.005, which would be somewhat smaller if 2002TC302 has a satellite. The best occultation light curves do not show any signs of ring features or any signatures of a global atmosphere., Comment: Accepted for publication in Astronomy & Astrophysics. 14 pages, 9 figures

Non-resolved thermal infrared observations enable studies of thermal and physical properties of asteroid surfaces provided the shape and rotational properties of the target are well determined via thermo-physical models. We used calibration-programme Herschel PACS data (70, 100, 160 $\mu$m) and state-of-the-art shape models derived from adaptive-optics observations and/or optical light curves to constrain for the first time the thermal inertia of twelve large main-belt asteroids. We also modelled previously well-characterised targets such as (1) Ceres or (4) Vesta as they constitute important benchmarks. Using the scale as a free parameter, most targets required a re-scaling $\sim$5\% consistent with what would be expected given the absolute calibration error bars. This constitutes a good cross-validation of the scaled shape models, although some targets required larger re-scaling to reproduce the IR data. We obtained low thermal inertias typical of large main belt asteroids studied before, which continues to give support to the notion that these surfaces are covered by fine-grained insulating regolith. Although the wavelengths at which PACS observed are longwards of the emission peak for main-belt asteroids, they proved to be extremely valuable to constrain size and thermal inertia and not too sensitive to surface roughness. Finally, we also propose a graphical approach to help examine how different values of the exponent used for scaling the thermal inertia as a function of heliocentric distance (i.e. temperature) affect our interpretation of the results., Comment: Accepted for publication in Astronomy & Astrophysics (preprint version)

Here we present the results of visible range light curve observations of ten Centaurs using the Kepler Space Telescope in the framework of the K2 mission. Well defined periodic light curves are obtained in six cases allowing us to derive rotational periods, a notable increase in the number of Centaurs with known rotational properties. The low amplitude light curves of (471931) 2013 PH44 and (250112) 2002 KY14 can be explained either by albedo variegations, binarity or elongated shape. (353222) 2009 YD7 and (514312) 2016 AE193 could be rotating elongated objects, while 2017 CX33 and 2012 VU85 are the most promising binary candidates due to their slow rotations and higher light curve amplitudes. (463368) 2012 VU85 has the longest rotation period, P=56.2h observed among Centaurs. The P>20h rotation periods obtained for the two potential binaries underlines the importance of long, uninterrupted time series photometry of solar system targets that can suitably be performed only from spacecraft, like the Kepler in the K2 mission, and the currently running TESS mission., Comment: 17 pages, 7 figures, 3 tables

Thanks to the Gaia mission, it will be possible to determine the masses of approximately hundreds of large main belt asteroids with very good precision. We currently have diameter estimates for all of them that can be used to compute their volume and hence their density. However, some of those diameters are still based on simple thermal models, which can occasionally lead to volume uncertainties as high as 20-30%. The aim of this paper is to determine the 3D shape models and compute the volumes for 13 main belt asteroids that were selected from those targets for which Gaia will provide the mass with an accuracy of better than 10%. We used the genetic Shaping Asteroids with Genetic Evolution (SAGE) algorithm to fit disk-integrated, dense photometric lightcurves and obtain detailed asteroid shape models. These models were scaled by fitting them to available stellar occultation and/or thermal infrared observations. We determine the spin and shape models for 13 main belt asteroids using the SAGE algorithm. Occultation fitting enables us to confirm main shape features and the spin state, while thermophysical modeling leads to more precise diameters as well as estimates of thermal inertia values. We calculated the volume of our sample of main-belt asteroids for which the Gaia satellite will provide precise mass determinations. From our volumes, it will then be possible to more accurately compute the bulk density, which is a fundamental physical property needed to understand the formation and evolution processes of small solar system bodies., Comment: Accepted for publication in A&A

Context. Earlier work suggests that slowly rotating asteroids should have higher thermal inertias than faster rotators because the heat wave penetrates deeper into the sub-surface. However, thermal inertias have been determined mainly for fast rotators due to selection effects in the available photometry used to obtain shape models required for thermophysical modelling (TPM). Aims. Our aims are to mitigate these selection effects by producing shape models of slow rotators, to scale them and compute their thermal inertia with TPM, and to verify whether thermal inertia increases with the rotation period. Methods. To decrease the bias against slow rotators, we conducted a photometric observing campaign of main-belt asteroids with periods longer than 12 hours, from multiple stations worldwide, adding in some cases data from WISE and Kepler space telescopes. For spin and shape reconstruction we used the lightcurve inversion method, and to derive thermal inertias we applied a thermophysical model to fit available infrared data from IRAS, AKARI, and WISE. Results. We present new models of 11 slow rotators that provide a good fit to the thermal data. In two cases, the TPM analysis showed a clear preference for one of the two possible mirror solutions. We derived the diameters and albedos of our targets in addition to their thermal inertias, which ranged between 3$^{+33}_{-3}$ and 45$^{+60}_{-30}$ Jm$^{-2}$s$^{-1/2}$K$^{-1}$. Conclusions. Together with our previous work, we have analysed 16 slow rotators from our dense survey with sizes between 30 and 150 km. The current sample thermal inertias vary widely, which does not confirm the earlier suggestion that slower rotators have higher thermal inertias., Comment: Accepted for publication in Astronomy & Astrophysics

We study the visible and near-infrared (NIR) spectral properties of different ACO populations and compare them to the independently determined properties of comets. We select our ACOs sample based on published dynamical criteria and present our own observational results obtained using the 10.4m Gran Telescopio Canarias (GTC), the 4.2m William Herschel Telescope (WHT), the 3.56m Telescopio Nazionale Galileo (TNG), and the 2.5m Isaac Newton Telescope (INT), all located at the El Roque de los Muchachos Observatory (La Palma, Spain), and the 3.0m NASA Infrared Telescope Facility (IRTF), located at the Mauna Kea Observatory, in Hawaii. We include in the analysis the spectra of ACOs obtained from the literature. We derive the spectral class and the visible and NIR spectral slopes. We also study the presence of hydrated minerals by studying the 0.7 $\mu$m band and the UV-drop below 0.5 $\mu$m associated with phyllosilicates. We present new observations of 17 ACOs, 11 of them observed in the visible, 2 in the NIR and 4 in the visible and NIR. We also discuss the spectra of 12 ACOs obtained from the literature. All but two ACOs have a primitive-like class spectrum (X or D-type). Almost 100\% of the ACOs in long-period cometary orbits (Damocloids) are D-types. Those in Jupiter family comet orbits (JFC-ACOs) are $\sim$ 60\% D-types and $\sim$ 40\% X-types. The mean spectral slope $S'$ of JFC-ACOs is 9.7 $\pm$ 4.6 \%/1000 \AA \ and for the Damocloids this is 12.2 $\pm$ 2.0 \%/1000 \AA . No evidence of hydration on the surface of ACOs is found from their visible spectra. The slope and spectral class distribution of ACOs is similar to that of comets. The spectral classification, the spectral slope distribution of ACOs, and the lack of spectral features indicative of the presence of hydrated minerals on their surface, strongly suggest that ACOs are likely dormant or extinct comets., Comment: 11 pages, 10 Figures, published in A&A

A recent occultation measurement of the dwarf planet Haumea (Ortiz et al. 2017) revealed an elongated shape with the longest axis comparable to Pluto's mean diameter. The chords also indicate a ring around Haumea's equatorial plane, where its largest moon, Hi'iaka, is also located. The Haumea occultation size estimate (equivalent diameter 1595 km) is larger than previous radiometric solutions (in the range between 1150 and 1350 km), which lowers the object's density to about 1.8 gcm-3, a value closer to the densities of other large TNOs. We present unpublished and also reprocessed Herschel and Spitzer MIR and FIR measurements. We compare 100 and 160um thermal lightcurve amplitudes - originating from Haumea itself - with models of the total measured system fluxes (ring, satellite, Haumea) from 24-350um. The combination with results derived from the occultation measurements show that Haumea's crystalline water ice surface must have a thermal inertia of ~5 SIu. We also have indications that the satellites (at least Hi'iaka) must have high geometric albedos >=0.5, otherwise the derived thermal amplitude would be inconsistent with the total measured system fluxes at 24, 70, 100, 160, 250, and 350um. The high albedos imply sizes of about 300 and 150 km for Hi'iaka and Namaka, respectively, indicating unexpectedly high densities >1.0 gcm-3 for TNOs this small, and the assumed collisional formation from Haumea's icy crust. We also estimated the thermal emission of the ring for the time period 1980-2030, showing that the contribution during the Spitzer and Herschel epochs was small, but not negligible. Due to the progressive opening of the ring plane, the ring emission will be increasing in the next decade when JWST is operational. In the MIRI 25.5um band it will also be possible to obtain a very high-quality thermal lightcurve to test the derived Haumea properties., Comment: 13 pages, 5 figures, 9 tables, accepted for publication in Icarus in November 2018; abstract has been shortened with respect to the original one

NASA's OSIRIS-REx and JAXA's Hayabusa2 sample-return missions are currently on their way to encounter primitive near-Earth asteroids (101955) Bennu and (162173) Ryugu, respectively. Spectral and dynamical evidence indicates that these near-Earth asteroids originated in the inner part of the main belt. There are several primitive collisional families in this region, and both these asteroids are most likely to have originated in the Polana-Eulalia family complex. We present the expected spectral characteristics of both targets based on our studies of our primitive collisional families in the inner belt: Polana-Eulalia, Erigone, Sulamitis, and Clarissa. Observations were obtained in the framework of our PRIMitive Asteroids Spectroscopic Survey (PRIMASS). Our results are especially relevant to the planning and interpretation of in-situ images and spectra to be obtained by the two spacecraft during the encounters with their targets., Comment: 22 pages, 11 figures. Accepted for publication in Icarus on May 11, 2018

The AKARI IRC All-sky survey provided more than twenty thousand thermal infrared observations of over five thousand asteroids. Diameters and albedos were obtained by fitting an empirically calibrated version of the standard thermal model to these data. After the publication of the flux catalogue in October 2016, our aim here is to present the AKARI IRC all-sky survey data and discuss valuable scientific applications in the field of small-body physical properties studies. As an example, we update the catalogue of asteroid diameters and albedos based on AKARI using the near-Earth asteroid thermal model (NEATM). We fit the NEATM to derive asteroid diameters and, whenever possible, infrared beaming parameters. We obtained a total of 8097 diameters and albedos for 5170 asteroids, and we fitted the beaming parameter for almost two thousand of them. When it was not possible to fit the beaming parameter, we used a straight line fit to our sample's beaming parameter-versus-phase angle plot to set the default value for each fit individually instead of using a single average value. Our diameters agree with stellar-occultation-based diameters well within the accuracy expected for the model. They also match the previous AKARI-based catalogue at phase angles lower than 50 degrees, but we find a systematic deviation at higher phase angles, at which near-Earth and Mars-crossing asteroids were observed. The AKARI IRC All-sky survey provides observations at different observation geometries, rotational coverages and aspect angles. For example, by comparing in more detail a few asteroids for which dimensions were derived from occultations, we discuss how the multiple observations per object may already provide three-dimensional information about elongated objects even based on an idealised model like the NEATM., Comment: Accepted for publication in Astronomy & Astrophysics

The available set of spin and shape modelled asteroids is strongly biased against slowly rotating targets and those with low lightcurve amplitudes. As a consequence of these selection effects, the current picture of asteroid spin axis distribution, rotation rates, or radiometric properties, might be affected too. To counteract these selection effects, we are running a photometric campaign of a large sample of main belt asteroids omitted in most previous studies. We determined synodic rotation periods and verified previous determinations. When a dataset for a given target was sufficiently large and varied, we performed spin and shape modelling with two different methods. We used the convex inversion method and the non-convex SAGE algorithm, applied on the same datasets of dense lightcurves. Unlike convex inversion, the SAGE method allows for the existence of valleys and indentations on the shapes based only on lightcurves. We obtained detailed spin and shape models for the first five targets of our sample: (159) Aemilia, (227) Philosophia, (329) Svea, (478) Tergeste, and (487) Venetia. When compared to stellar occultation chords, our models obtained an absolute size scale and major topographic features of the shape models were also confirmed. When applied to thermophysical modelling, they provided a very good fit to the infrared data and allowed their size, albedo, and thermal inertia to be determined. Convex and non-convex shape models provide comparable fits to lightcurves. However, some non-convex models fit notably better to stellar occultation chords and to infrared data in sophisticated thermophysical modelling (TPM). In some cases TPM showed strong preference for one of the spin and shape solutions. Also, we confirmed that slowly rotating asteroids tend to have higher-than-average values of thermal inertia., Comment: Accepted in astronomy astrophys. 39 pages, 49 figures

The Cybele and Hilda dynamical groups delimit the outer edge of the asteroid belt. Their compositional distribution is a key element to constrain evolutionary models of the Solar System. In this paper, we present a compositional analysis of these populations using spectroscopic observations, SDSS and NEOWISE data. As part of the PRIMASS (Primitive Asteroids Spectroscopic Survey), we acquired visible spectra of 18 objects in Hilda or Cybele groups with the Goodman High Throughput Spectrometer at the 4.1m SOAR telescope and 20 near-IR spectra of Hilda objects with Near Infrared Camera Spectrograph at the 3.56m TNG. The sample is enlarged with spectra taken from the literature in order to increase our statistical analysis. The spectra were inspected for aqueous alteration bands and other spectral features that can be linked to compositional constraints. The analysis shows a continuous distribution of compositions from the main-belt to the Cybele, Hilda and Trojan regions. We also identify a population in the Trojans group not present in Hilda or Cybele objects., Comment: 26 pages, 10 figures

The combination of visible and thermal data from the ground and astrophysics space missions is key to improving the scientific understanding of near-Earth, main-belt, trojans, centaurs, and trans-Neptunian objects. To get full information on a small sample of selected bodies we combine different methods and techniques: lightcurve inversion, stellar occultations, thermophysical modeling, radiometric methods, radar ranging and adaptive optics imaging. The SBNAF project will derive size, spin and shape, thermal inertia, surface roughness, and in some cases bulk densities and even internal structure and composition, for objects out to the most distant regions in the Solar System. The applications to objects with ground-truth information allows us to advance the techniques beyond the current state-of-the-art and to assess the limitations of each method. We present results from our project's first phase: the analysis of combined Herschel-KeplerK2 data and Herschel-occultation data for TNOs; synergy studies on large MBAs from combined high-quality visual and thermal data; establishment of well-known asteroids as celestial calibrators for far-infrared, sub-millimetre, and millimetre projects; first results on near-Earth asteroids properties from combined lightcurve, radar and thermal measurements, as well as the Hayabusa-2 mission target characterisation. We also introduce public web-services and tools for studies of small bodies in general., Comment: Accepted for publication in Advances in Space Research, 43 pages, 5 figures

Eos family was created during a catastrophic impact about 1.3 Gyr ago. Rotation states of individual family members contain information about the history of the whole population. We aim to increase the number of asteroid shape models and rotation states within the Eos collision family, as well as to revise previously published shape models from the literature. Such results can be used to constrain theoretical collisional and evolution models of the family, or to estimate other physical parameters by a thermophysical modeling of the thermal infrared data. We use all available disk-integrated optical data (i.e., classical dense-in-time photometry obtained from public databases and through a large collaboration network as well as sparse-in-time individual measurements from a few sky surveys) as input for the convex inversion method, and derive 3D shape models of asteroids together with their rotation periods and orientations of rotation axes. We present updated shape models for 15 asteroids and new shape model determinations for 16 asteroids. Together with the already published models from the publicly available DAMIT database, we compiled a sample of 56 Eos family members with known shape models that we used in our analysis of physical properties within the family. Rotation states of asteroids smaller than ~20 km are heavily influenced by the YORP effect, whilst the large objects more or less retained their rotation state properties since the family creation. Moreover, we also present a shape model and bulk density of asteroid (423) Diotima, an interloper in the Eos family, based on the disk-resolved data obtained by the Near InfraRed Camera (Nirc2) mounted on the W.M. Keck II telescope., Comment: Accepted for publication in ICARUS Special Issue - Asteroids: Origin, Evolution & Characterization

Context. Many databases on asteroid brightnesses (e.g. ALCDEF, WISE) are potential sources for extensive asteroid shape and spin modelling. Individual lightcurve inversion models require several apparitions and hundreds of data points per target. However, we can analyse the coarse shape and spin distributions over populations of at least thousands of targets even if there are only a few points and one apparition per asteroid. This is done by examining the distribution of the brightness variations observed within the chosen population. Aims. Brightness variation has been proposed as a population-scale rather than individual-target observable in two studies so far. We aim to examine this approach rigorously to establish its theoretical validity, degree of ill-posedness, and practical applicability. Methods. We model the observed brightness variation of a target population by considering its cumulative distribution function (CDF) caused by the joint distribution function of two fundamental shape and spin indicators. These are the shape elongation and the spin latitude of a simple ellipsoidal model. The main advantage of the model is that we can derive analytical basis functions that yield the observed CDF as a function of the shape and spin distribution. The inverse problem can be treated linearly. Even though the inaccuracy of the model is considerable, databases of thousands of targets should yield some information on the distribution. Results. We establish the theoretical soundness and the typical accuracy limits of the approach both analytically and numerically. Using simulations, we derive a practical estimate of the model distribution in the (shape, spin)-plane. We show that databases such as Wide-field Infrared Survey Explorer (WISE) yield coarse but robust estimates of this distribution, and as an example compare various asteroid families with each other., Comment: 16 pages, 21 figures, manuscript accepted in Astronomy & Astrophysics, to be published in section 10. Planets and planetary systems

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., Comment: Astronomy and Astrophysics. In press

The potentially hazardous asteroid (PHA) (99942) Apophis is one of the most remarkable near-Earth asteroids (NEA) in terms of impact hazard. A good determination of its surface thermal inertia is very important in order to evaluate the Yarkovsky effect on its orbital evolution. We present thermal infrared observations obtained on January 29, 2013, with CanariCam mid-infrared camera/spectrograph attached to the Gran Telescopio CANARIAS (GTC, Roque de los Muchachos Observatory, La Palma, Spain) using the Si2-8.7, Si6-12.5, and Q1-17.65 filters with the aim of deriving Apophis' diameter ($D$), geometric albedo ($p_V$), and thermal inertia ($\Gamma$). We performed a detailed thermophysical model analysis of the GTC data combined with previously published thermal data obtained using Herschel Space Observatory PACS instrument at 70, 100, and 160 $\mu$m.The thermophysical model fit of the data favors low surface roughness solutions (within a range of roughness slope angles $rms$ between 0.1 and 0.5), and constrains the effective diameter, visible geometric albedo, and thermal inertia of Apophis to be $D_{eff} =$~380 -- 393 m, $p_V = $~0.24--0.33 (assuming absolute magnitude $H = 19.09 \pm 0.19$) and $\Gamma =$~50 -- 500 Jm$^{-2}$ s$^{-0.5}$ K$^{-1}$, respectively., Comment: 4 pages, 3 figures

We study the distributions of effective diameter ($D$), beaming parameter ($\eta$), and visible geometric albedo ($p_V$) of asteroids in cometry orbits (ACOs) populations, derived from NASA's Wide-field Infrared Explorer (WISE) observations, and compare these with the same, independently determined properties of the comets. The near-Earth asteroid thermal model (NEATM) is used to compute the $D$, $p_V$ and $\eta$. We obtained $D$ and $p_V$ for 49 ACOs in Jupiter family cometary orbits (JF-ACOs) and 16 ACOs in Halley-type orbits (Damocloids). We also obtained $\eta$ for 45 of them. All but three JF-ACOs (95% of the sample) present a low albedo compatible with a cometary origin. The $p_V$ and $\eta$ distributions of both ACO populations are very similar. For the entire sample of ACOs, the mean geometric albedo is $\bar{p_V} = 0.05 \pm 0.02$, ($\bar{p_V} = 0.05 \pm 0.01$ and $\bar{p_V} =0.05 \pm 0.02$ for JF-ACOs and Damocloids, respectively) compatible with a narrow albedo distribution similar to that of the Jupiter family comets (JFCs), with a $\bar{p_V} \sim 0.04$. The $\bar{\eta} =1.0 \pm 0.2$. We find no correlations between $D$, $p_V$ , or $\eta$. We compare the cumulative size distribution (CSD) of ACOs, Centaurs, and JFCs. Although the Centaur sample contains larger objects, the linear parts in their log-log plot of the CSDs presents a similar cumulative exponent ($\beta = 1.85 \pm 0.30$ and $1.76 \pm 0.35$, respectively). The CSD for Damocloids presents a much shallower exponent $\beta = 0.89 \pm 0.17$. The CSD for JF-ACOs is shallower and shifted towards larger diameters with respect to the CSD of active JFCs, which suggests that the mantling process has a size dependency whereby large comets tend to reach an inactive stage faster than small ones. Finally, the population of JF-ACOs is comparable in number that of JFCs, although there are more tens-km JF-ACOs than JFCs., Comment: 10 pages, 10 figures, A&A in press

The so-called "early activity" of comet 67P/Churyumov-Gerasimenko has been observed to originate mostly in parts of the concave region or "neck" between its two lobes. Since activity is driven by the sublimation of volatiles, this is a puzzling result because this area is less exposed to the Sun and is therefore expected to be cooler on average (Sierks et al. 2015). We used a thermophysical model that takes into account thermal inertia, global self-heating, and shadowing, to compute surface temperatures of the comet. We found that, for every rotation in the August--December 2014 period, some parts of the neck region undergo the fastest temperature variations of the comet's surface precisely because they are shadowed by their surrounding terrains. Our work suggests that these fast temperature changes are correlated to the early activity of the comet, and we put forward the hypothesis that erosion related to thermal cracking is operating at a high rate on the neck region due to these rapid temperature variations. This may explain why the neck contains some ice --as opposed to most other parts of the surface (Capaccioni et al. 2015)-- and why it is the main source of the comet's early activity (Sierks et al. 2015). In a broader context, these results indicate that thermal cracking can operate faster on atmosphereless bodies with significant concavities than implied by currently available estimates (Delbo' et al. 2014)., Comment: 11 pages, 4 figures

We present THERMAP, a mid-infrared (8-16 {\mu}m) spectro-imager for space missions to small bodies in the inner solar system, developed in the framework of the MarcoPolo-R asteroid sample return mission. THERMAP is very well suited to characterize the surface thermal environment of a NEO and to map its surface composition. The instrument has two channels, one for imaging and one for spectroscopy: it is both a thermal camera with full 2D imaging capabilities and a slit spectrometer. THERMAP takes advantage of the recent technological developments of uncooled microbolometers detectors, sensitive in the mid-infrared spectral range. THERMAP can acquire thermal images (8-18 {\mu}m) of the surface and perform absolute temperature measurements with a precision better than 3.5 K above 200 K. THERMAP can acquire mid-infrared spectra (8-16 {\mu}m) of the surface with a spectral resolution {\Delta}{\lambda} of 0.3 {\mu}m. For surface temperatures above 350 K, spectra have a signal-to-noise ratio >60 in the spectral range 9-13 {\mu}m where most emission features occur., Comment: 21 pages, 12 figures, 3 tables; Experimental Astronomy: in press

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

Aims. To derive the thermal inertia of 2008 EV$_5$, the baseline target for the Marco Polo-R mission proposal, and infer information about the size of the particles on its surface. Methods. Values of thermal inertia are obtained by fitting an asteroid thermophysical model to NASA's Wide-field Infrared Survey Explorer (WISE) infrared data. From the constrained thermal inertia and a model of heat conductivity that accounts for different values of the packing fraction (a measure of the degree of compaction of the regolith particles), grain size is derived. Results. We obtain an effective diameter $D = 370 \pm 6\,\mathrm{m}$, geometric visible albedo $p_V = 0.13 \pm 0.05$ (assuming $H=20.0 \pm 0.4$), and thermal inertia $\Gamma = 450 \pm 60$ J/m2/s(1/2)/K at the 1-$\sigma$ level of significance for its retrograde spin pole solution. The regolith particles radius is $r = 6.6^{+1.3}_{-1.3}$ mm for low degrees of compaction, and $r = 12.5^{+2.7}_{-2.6}$ mm for the highest packing densities., Comment: 16 pages, 8 figures; accepted for publication in Astronomy & Astrophysics

Near-Earth binary asteroid (175706) 1996 FG3 is the current backup target of the ESA MarcoPolo-R mission, selected for the study phase of ESA M3 missions. It is a primitive (C-type) asteroid that shows significant variation in its visible and near-infrared spectra. Here we present new spectra of 1996 FG3 and we compare our new data with other published spectra, analysing the variation in the spectral slope. The asteroid will not be observable again over the next three years at least. We obtained the spectra using DOLORES and NICS instruments at the Telescopio Nazionale Galileo (TNG), a 3.6m telescope located at El Roque de los Muchachos Observatory in La Palma, Spain. To compare with other published spectra of the asteroid, we computed the spectral slope S', and studied any plausible correlation of this quantity with the phase angle (alpha). In the case of visible spectra, we find a variation in spectral slope of Delta S' = 0.15 +- 0.10 %/10^3 A/degree for 3 < alpha < 18 degrees, in good agreement with the values found in the literature for the phase reddening effect. In the case of the near-infrared, we find a variation in the slope of Delta S' = 0.04 +- 0.08 %/10^3 A/degree for 6 < alpha < 51 degrees. Our computed variation in S' agrees with the only two values found in the literature for the phase reddening in the near-infrared. The variation in the spectral slope of asteroid 1996 FG3 shows a trend with the phase angle at the time of the observations, both in the visible and the near-infrared. It is worth noting that, to fully explain this spectral variability we should take into account other factors, like the position of the secondary component of the binary asteroid 1999 FG3 with respect to the primary, or the spin axis orientation at the time of the observations. More data are necessary for an analysis of this kind., Comment: 4 pages, 3 figures, Accepted in A&A 25 June 2013

Most of the objects in the trans-Neptunian belt (TNb) and related populations move in prograde orbits with low eccentricity and inclination. However, the list of icy minor bodies moving in orbits with an inclination above 40 deg. has increased in recent years. The origin of these bodies, and in particular of those objects in retrograde orbits, is not well determined, and different scenarios are considered. In this paper, we present new observational and dynamical data of two objects in retrograde orbits, 2008 YB3 and 2005 VD. We find that the surface of these extreme objects is depleted of ices and does not contain the 'ultra-red' matter typical of some Centaurs. Despite small differences, these objects share common colors and spectral characteristics with the Trojans, comet nuclei, and the group of grey Centaurs. All of these populations are supposed to be covered by a mantle of dust responsible for their reddish- to neutral-color. To investigate if the surface properties and dynamical evolution of these bodies are related, we integrate their orbits for 10^(8) years to the past. We find a remarkable difference in their dynamical evolutions: 2005 VD' s evolution is dominated by a Kozai resonance with planet Jupiter while that of 2008 YB3 is dominated by close encounters with planets Jupiter and Saturn. Our models suggest that the immediate site of provenance of 2005 VD is the in the Oort cloud, whereas for 2008 YB3 it is in the trans-Neptunian region. Additionally, the study of their residence time shows that 2005 VD has spent a larger lapse of time moving in orbits in the region of the giant planets than 2008 YB3. Together with the small differences in color between these two objects, with 2005 VD being more neutral than 2008 YB3, this fact suggests that the surface of 2005 VD has suffered a higher degree of processing, probably related to cometary activity episodes., Comment: 9 pages, 5 figures, Accepted by Astronomy and Astrophysics: 28 November 2012