Aims. The goal of this work is to characterize the ensemble thermal properties of the Centaurs / trans-Neptunian population. Methods. Thermal flux measurements obtained with Herschel/PACS and Spitzer/MIPS provide size, albedo, and beaming factors for 85 objects (13 of which are presented here for the first time) by means of standard radiometric techniques. The measured beaming factors are influenced by the combination of surface roughness and thermal inertia effects. They are interpreted within a thermophysical model to constrain, in a statistical sense, the thermal inertia in the population and to study its dependence on several parameters. We use in particular a Monte-Carlo modeling approach to the data whereby synthetic datasets of beaming factors are created using random distributions of spin orientation and surface roughness. Results. Beaming factors η range from values <1 to ∼2.5, but high η values (>2) are lacking at low heliocentric distances (r h < 30 AU). Beaming factors lower than 1 occur frequently (39% of the objects), indicating that surface roughness effects are important. We determine a mean thermal inertia for Centaurs/ TNO of Γ = (2.5 ± 0.5) J m −2 s −1/2 K −1 , with evidence of a trend toward decreasing Γ with increasing heliocentric (by a factor ∼2.5 from 8-25 AU to 41-53 AU). These thermal inertias are 2-3 orders of magnitude lower than expected for compact ices, and generally lower than on Saturn's satellites or in the Pluto/Charon system. Most high-albedo objects are found to have unusually low thermal inertias. Our results suggest highly porous surfaces, in which the heat transfer is affected by radiative conductivity within pores and increases with depth in the subsurface.

Context. Physical characterization of trans-Neptunian objects, a primitive population of the outer solar system, may provide constraints on their formation and evolution. Aims. The goal of this work is to characterize a set of 15 scattered disk (SDOs) and detached objects, in terms of their size, albedo, and thermal properties. Methods. Thermal flux measurements obtained with the Herschel-PACS instrument at 70, 100 and 160 μm, and whenever applicable, with Spitzer-MIPS at 24 and 70 μm, are modeled with radiometric techniques, in order to derive the objects' individual size, albedo and when possible beaming factor. Error bars are obtained from a Monte-Carlo approach. We look for correlations between these and other physical and orbital parameters. Results. Diameters obtained for our sample range from 100 to 2400 km, and the geometric albedos (in V band) vary from 3.8% to 84.5%. The unweighted mean V geometric albedo for the whole sample is 11.2% (excluding Eris); 6.9% for the SDOs, and 17.0% for the detached objects (excluding Eris). We obtain new bulk densities for three binary systems: Ceto/Phorcys, Typhon/Echidna and Eris/Dysnomia. Apart from correlations clearly due to observational bias, we find significant correlations between albedo and diameter (more reflective objects being bigger), and between albedo, diameter and perihelion distance (brighter and bigger objects having larger perihelia). We discuss possible explanations for these correlations.

Since the discovery of the Kuiper belt, broadband surface colors were thoroughly studied as a first approximation to the object reflectivity spectra. Visible colors (BVRI) have proven to be a reasonable proxy for real spectra, which are rather linear in this range. In contrast, near-IR colors (JHK bands) could be misleading when absorption features of ices are present in the spectra. Although the physical and chemical information provided by colors are rather limited, broadband photometry remains the best tool for establishing the bulk surface properties of Kuiper belt objects (KBOs) and Centaurs. In this work, we explore for the first time general, recurrent effects in the study of visible colors that could affect the interpretation of the scientific results: i) how a correlation could be missed or weakened as a result of the data error bars; ii) the "risk" of missing an existing trend because of low sampling, and the possibility of making quantified predictions on the sample size needed to detect a trend at a given significance level -assuming the sample is unbiased; iii) the use of partial correlations to distinguish the mutual effect of two or more (physical) parameters; and iv) the sensitivity of the "reddening line" tool to the central wavelength of the filters used. To illustrate and apply these new tools, we have compiled the visible colors and orbital parameters of about 370 objects available in the literature − assumed, by default, as unbiased samplesand carried out a traditional analysis per dynamical family. Our results show in particular how a) data error bars impose a limit on the detectable correlations regardless of sample size and that therefore, once that limit is achieved, it is important to diminish the error bars, but it is pointless to enlarge the sampling with the same or larger errors; b) almost all dynamical families still require larger samplings to ensure the detection of correlations stronger than ±0.5, that is, correlations that may explain ∼25% or more of the color variability; c) the correlation strength between (V − R) vs. (R − I) is systematically lower than the one between (B − V) vs. (V − R) and is not related with error-bar differences between these colors; d) it is statistically equivalent to use any of the different flavors of orbital excitation or collisional velocity parameters regarding the famous color-inclination correlation among classical KBOs − which no longer appears to be a strong correlation -whereas the inclination and Tisserand parameter relative to Neptune cannot be separated from one another; and e) classical KBOs are the only dynamical family that shows neither (B − V) vs. (V − R) nor (V − R) vs. (R − I) correlations. It therefore is the family with the most unpredictable visible surface reflectivities.

No abstract