Issue 
A&A
Volume 541, May 2012



Article Number  A93  
Number of page(s)  17  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201118562  
Published online  04 May 2012 
TNOs are cool: A survey of the transNeptunian region
V. Physical characterization of 18 Plutinos using HerschelPACS observations^{⋆}
^{1}
Deutsches Zentrum für Luft und Raumfahrt (DLR), Institut für
Planetenforschung,
Rutherfordstraße 2,
12489
Berlin,
Germany
email: michael.mommert@dlr.de
^{2} Konkoly Observatory of the Hungarian Academy of Sciences,
1525 Budapest, PO Box 67, Hungary
^{3}
LESIAObservatoire de Paris, CNRS, UPMC Univ Paris 06, Univ. ParisDiderot5 place J.
Janssen, 92195
Meudon Principal Cedex,
France
^{4}
MaxPlanckInstitut für Sonnensystemforschung (MPS),
MaxPlanckStraße
2, 37191
KatlenburgLindau,
Germany
^{5}
The University of Arizona, Tucson, AZ
85721,
USA
^{6}
MaxPlanckInstitut für extraterrestrische Physik (MPE),
Postfach 1312,
Giessenbachstr., 85741
Garching,
Germany
^{7}
Center for Geophysics of the University of Coimbra,
Av. Dr. Dias da
Silva, 3000134
Coimbra,
Portugal
^{8}
Astronomical Observatory of the University of Coimbra,
Almas de Freire,
304004
Coimbra,
Portugal
^{9}
Instituto de Astrofísica de Andalucía (CSIC) C/ Bajo de Huétor 50,
18008
Granada,
Spain
^{10}
SRON Netherlands Institute for Space Research,
Postbus 800,
9700 AV
Groningen, The
Netherlands
^{11}
UNSCNRSObservatoire de la Côte d’Azur, Laboratoire Cassiopée,
BP 4229,
06304
Nice Cedex 04,
France
^{12}
Univ. Paris Diderot, Sorbonne Paris Cité, 4 rue Elsa Morante,
75205
Paris,
France
^{13}
Laboratoire d’Astrophysique de Marseille, CNRS &
Université de Provence, 38 rue
Frédéric JoliotCurie, 13388
Marseille Cedex 13,
France
^{14}
Department of Astronomy, University of Maryland,
College Park, MD
20742,
USA
Received:
1
December
2011
Accepted:
31
January
2012
Context. The Herschel open time key programme TNOs are Cool: A survey of the transNeptunian region aims to derive physical and thermal properties for a set of ~140 Centaurs and transNeptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. One goal of the project is to determine albedo and size distributions for specific classes and the overall population of TNOs.
Aims. We present HerschelPACS photometry of 18 Plutinos and determine sizes and albedos for these objects using thermal modeling. We analyze our results for correlations, draw conclusions on the Plutino size distribution, and compare to earlier results.
Methods. Flux densities are derived from PACS mini scanmaps using specialized data reduction and photometry methods. In order to improve the quality of our results, we combine our PACS data with existing Spitzer MIPS data where possible, and refine existing absolute magnitudes for the targets. The physical characterization of our sample is done using a thermal model. Uncertainties of the physical parameters are derived using customized Monte Carlo methods. The correlation analysis is performed using a bootstrap Spearman rank analysis.
Results. We find the sizes of our Plutinos to range from 150 to 730 km and geometric albedos to vary between 0.04 and 0.28. The average albedo of the sample is 0.08 ± 0.03, which is comparable to the mean albedo of Centaurs, Jupiter family comets and other TNOs. We were able to calibrate the Plutino size scale for the first time and find the cumulative Plutino size distribution to be best fit using a cumulative power law with q = 2 at sizes ranging from 120–400 km and q = 3 at larger sizes. We revise the bulk density of 1999 TC36 and find ρ = 0.64_{0.11}^{+0.15} g cm^{3}. On the basis of a modified Spearman rank analysis technique our Plutino sample appears to be biased with respect to object size but unbiased with respect to albedo. Furthermore, we find biases based on geometrical aspects and color in our sample. There is qualitative evidence that icy Plutinos have higher albedos than the average of the sample.
Key words: Kuiper belt: general / infrared: planetary systems / methods: observational / techniques: photometric
Herschel is an ESA space observatory with science instruments provided by Europeanled Principal Investigator consortia and with important participation from NASA. The Photodetector Array Camera and Spectrometer (PACS) is one of Herschel’s instruments. PACS has been developed by a consortium of institutes led by MPE (Germany) and including UVIE (Austria); KU Leuven, CSL, IMEC (Belgium); CEA, LAM (France); MPIA (Germany); INAFIFSI/OAA/OAP/OAT, LENS, SISSA (Italy); IAC (Spain). This development has been supported by the funding agencies BMVIT (Austria), ESAPRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI/INAF (Italy), and CICYT/MCYT (Spain).
© ESO, 2012
1. Introduction
Since its discovery in 1930, Pluto has been a unique object not only for being the only rocky planetsized object not bound to a planet outside the orbit of Mars, but also for having the most eccentric and inclined orbit of all the planets, which even overlaps the orbit of Neptune. Cohen & Hubbard (1965) were the first to show that, despite that overlap, close approaches between the planets are prevented by the 2:3 mean motion resonance: Pluto’s revolution period equals 3/2 of Neptune’s period, ensuring that conjunctions always occur near Pluto’s aphelion. This leads to a high degree of stability of the orbit. The origin of Pluto’s peculiar orbit (e = 0.25 and i = 17°) was first explained by Malhotra (1993), who showed that encounters of the Jovian planets with residual planetesimals during the late stages of the formation of the solar system could lead to a radial migration of the former. As a result of Neptune’s outward migration, a Plutolike body could have been captured in the 2:3 resonance, excited to its highly inclined and eccentric orbit and transported outward.
In the meantime, the discovery of 1992 QB1 by Jewitt & Luu (1993) showed that Pluto is not the only object beyond the orbit of Neptune. The following years revealed a large population of transNeptunian objects (TNOs), presumably a population of residual planetesimals from the age of the formation of the solar system, as previously proposed by Edgeworth (1949) and Kuiper (1951). The population shows dynamical complexity and recurring orbital characteristics allowing the classification of TNOs in different dynamical groups (cf. Elliot et al. 2005; Gladman et al. 2008). Many of the newly discovered objects show a dynamical behaviour similar to that of Pluto, which led Jewitt & Luu (1996) to dub them Plutinos.
Plutinos reside in Neptune’s 2:3 resonance (Gladman et al. 2008) and are the most numerous resonant population. This observation agrees with dynamical studies by Melita & Brunini (2000), which reveal that the 2:3 resonance is much more stable than other resonances, particularly at low inclinations. While the semimajor axes of Plutinos are strongly concentrated around 39.5 AU, their eccentricities and inclinations vary significantly from Pluto’s (0.03 ≤ e ≤ 0.88 and 0.4° ≤ i ≤ 40.2° with a mean eccentricity and inclination of 0.21 and 12.1°, respectively)^{1}. The origin of the Plutino population is ascribed to the same resonance capture mechanism which is responsible for Pluto’s peculiar orbit (Malhotra 1995). Based on numerical simulations, Duncan et al. (1995) find that the 2:3 mean motion resonance interferes with secular resonances leading to a severe instability for high inclination orbits and a longterm leakage for low inclination orbits of Plutinos, making them a possible source region of Jupiter family comets (JFCs). Morbidelli (1997) found that a slow chaotic diffusion of Plutinos can be provided, which is necessary to explain the observed continuous flux of JFCs. Their results were recently confirmed by di Sisto et al. (2010), who found that Plutinos may also be a secondary source of the Centaur population. A collisional analysis of the 2:3 resonance carried out by de Elía et al. (2008) shows that the Plutino population larger than a few kilometers in diameter is not significantly altered by catastrophic collisions over the age of the solar system. They also pointed out the importance of specifying the number of Plutosized objects among the Plutino population, since the escape frequency of Plutinos strongly depends on this number.
Fig. 1 Our Plutino sample (filled circles) compared to the known Plutino population (open circles, data from MPC, as of 18 Nov. 2011) in e–i phase space. The crosses with solid and dashed error bars refer to the mean and the standard deviation of our sample and the known Plutino population, respectively. This plot demonstrates that our sample is heterogenous and reflects the dynamical spread of the complete Plutino population. (The plot omits 3 known Plutinos, due to their high eccentricities, which are not part of our sample; however, these objects are included in the determination of the mean.) 
Plutino target sample.
PACS observations summary.
Absolute magnitude H, lightcurve information, information on the presence of ices and color indices.
Our sample of 18 Plutinos is selected based on the object classification scheme of Gladman et al. (2008) and represents ~7.5% of the known Plutino population as of 27 October, 2011 (data from MPC). 14 of our 18 sample targets have been classified in Gladman et al. (2008). For the remaining ones dynamical simulations performed by Ch. Ejeta show that their Plutinotype orbits are stable for at least 10^{7} years, which suggests, but does not prove, orbital stability over the lifetime of the solar system. We give the numbers and preliminary designations of the sample objects in Table 1. In order to avoid confusion, we will refer to specific objects in this work only by their preliminary designation or name, where available. In Fig. 1 we plot our dataset in comparison to the sample of all known Plutinos in e − i phase space. We neglect a here, since it is very well confined around 39.5 AU. The plot shows the spread of both samples in e − i space, whereas the means of the two samples show a good agreement. Therefore, we assume that our sample well represents the dynamical variety of the whole Plutino population.
Our Plutino sample is selected solely on the basis of Herschel observability in the sense of observation geometry and detector sensitivity. Herschel observations presented here are conducted over a time span of more than 7 months. In this time span, Herschel is able to observe at all ecliptical longitudes, despite its restrictions on the solar aspect angle, which must not be lower than 60.8° and not be larger than 119.2°. Hence, no bias is introduced based on preferred ecliptical longitudes. Since the target objects of the “TNOs are Cool” project are selected on the basis of optical discoveries, the most important bias in our sample is presumably the optical detection bias, favoring large objects and objects with high albedos. The detectional bias impact on our results will be discussed in Sect. 4.5.
Our target sample includes 3 known multiple systems, which are Pluto, 1999 TC36 and 2003 AZ84 and the expected binary 2002 GN171. We emphasize, that the spatial resolution of PACS is not sufficient to separate any multiple system. Flux densities measured for one object therefore consist of contributions from fluxes of all system components. Diameter and albedo determined with our models refer to respective parameters of an object with the area equivalent diameter of the whole system. We discuss these objects in detail in Sect. 4.2.
The “TNOs are Cool: A Survey of the transNeptunian region” project is a Herschel open time key programme awarded some 370 h of Herschel observing time for investigating about 140 TNOs with known orbits (Müller et al. 2009). The observations include PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) pointsource photometry. The goal is to characterize the diameter and albedo for the individual objects and the full sample using radiometric techniques, in order to probe formation and evolution processes in the solar system. Using thermal modeling we determine sizes and albedos of the 18 Plutino members. The “TNOs are Cool” sample includes data of 7 more Plutinos, which have not been processed by the time of writing this.
2. Observations and data reduction
2.1. Herschel observations
Photometric measurements using the Photometer Array Camera and Spectrometer (PACS, Poglitsch et al. 2010) onboard the Herschel Space Observatory (Pilbratt et al. 2010) have been taken in mini scan map mode covering homogeneously a field of roughly 1′ in diameter. This mode turned out to be best suited for our needs and offers more sensitivity than other observation modes (Müller et al. 2010).
According to the “TNOs are Cool!” open time key programme observation strategy a target is observed at two epochs, separated by a time interval that corresponds to a movement of 30–50′′ of the target, allowing for an optimal background subtraction, to eliminate confusion noise and background sources. At each epoch the target is observed in the “blue” (nominal wavelength of 70 μm) and “green” (100 μm) band twice, using two different scan position angles. “Red” (160 μm) band data are taken in parallel when sampling one of the other bands. This forms a series of 4 measurements in the blue and green bands, and a series of 8 measurements in the red band for a specific target. The maps are taken in the medium scan speed (20′′/s) mode, using a scanleg length of 3′ and 2–4 repetitions. In the case of 2002 VR128 a scanleg length of 2.5′ was used. More details on the observation planning can be found in Vilenius et al. (2012).
Observational circumstances and PACS flux densities are summarized in Table 2. Additional information on the targets is given in Table 3.
2.2. HerschelPACS data processing
The raw PACS measurements were used to produce individual scan maps (level2) using an optimized version of the standard PACS mini scan map pipeline in HIPE (Ott 2010). The individual scan maps of the same epoch and band were mosaicked using the MosaicTask(). In the production of the final maps backgroundmatching and source stacking techniques are applied to correct for the possible relative astrometric uncertainties between the two visits. We created two sorts of backgroundeliminated products from the “pervisit” mosaics: “backgroundsubtracted” (see Stansberry et al. 2008; SantosSanz et al. 2012) and “doubledifferential” maps^{2}. The latter are generated by subtracting the maps of the two visits, yielding a positive and a negative beam of the object on the differential map with all background structures eliminated. A duplicate of this map is then shifted in such a way as to match the positive beam of the original map with the negative one of the duplicate map. In the last step the original and shifted mosaics are subtracted again and averaged, resulting in a doubledifferential mosaic with a positive beam showing the full target flux, and two negative beams on either side showing half the target flux. The final photometry is performed on the central, positive beam. The advantage of this technique is the nearly complete elimination of the sky background that makes it favorable for faint targets. Full details can be found in Kiss et al. (in prep.).
Photometry is performed on the final backgroundsubtracted and/or doubledifferential maps, which are both backgroundeliminated. Flux densities are derived via aperture photometry using either IRAF/Daophot or IDL/Astrolib routines, both producing identical results. We extract the flux at the photocenter position of the target and apply an aperture correction technique (Howell 1989) for each aperture radius based on the encircled energy fraction for the PACS point spread function (Müller et al. 2011a). We construct an aperturecorrected curve of growth from which we derive the optimum synthetic aperture, which usually lies in the “plateau” region of the curve of growth.
Photometric uncertainties are estimated by random implantation of 200 artificial sources in the nearby sky background of the target. The 1σ photometric uncertainty of the target flux is derived as the standard deviation of these artificial source fluxes (SantosSanz et al. 2012). This 1σ limit is also given as an upper limit in the case of nondetections.
The absolute calibration of our data is based on mini scan maps of 5 fiducial stars and 18 large main belt asteroids (Müller et al. 2011a). The absolute calibration uncertainties of the PACS bands are 3% for the 70 and 100 μm bands and 5% for the 160 μm band, respectively. This additional uncertainty adds quadratically to the photometric error and is included in the flux densities given in Table 2. The absolute accuracy of the photometry was checked against the predicted standard star flux densities. In the bright regime (above 30 mJy) the relative accuracy was found to be ~5%, while in the faint regime (below 30 mJy) the accuracy is driven by the uncertainties in the maps and is about or below 2 mJy in all bands, on the individual maps.
A general description of the PACS data reduction steps (including photometry) will be given in Kiss et al. (in prep.).
Spitzer observations summary.
The flux densities given in Table 2 have been color corrected and therefore are monochromatic. The measured flux density is determined by the response function of the PACS band filters convolved with the spectral energy distribution (SED) of the source (Müller et al. 2011b). The SED of TNOs resembles a black body spectrum of a certain temperature; the resulting colorcorrection factors depend weakly on the emission temperature. The black body temperature is approximated by the mean surface temperature, which is given by 2^{0.25}T_{SS}, where T_{SS} is the subsolar temperature, defined in Eq. (3). We determine the color correction factors based on the mean surface temperature in an iterative process during the modeling. The mean color correction factors for all objects in this sample are 0.9810 ± 0.0003, 0.988 ± 0.003 and 1.015 ± 0.005 for the 70, 100 and 160 μm band, respectively: the size of the color corrections are smaller than the uncertainty of the individual flux density measurements (cf. Table 2), which justifies the use of the mean surface temperature as an approximation of the best fit black body temperature.
We have obtained timeresolved lightcurve observations of 2003 AZ84 and 2003 VS2, covering 110% and 106% of the respective lightcurve using 99 and 100 observations, respectively. A detailed analysis of both lightcurves will be subject of an upcoming publication (SantosSanz et al., in prep.). In this work, we make use of the averaged fluxes, which are more precise than single measurements, due to the larger number of observations and the cancelling of lightcurve effects.
2.3. Spitzer MIPS observations
Whenever possible, the Herschel flux densities were combined with existing flux density measurements of the Multiband Imaging Photometer for Spitzer (MIPS, Rieke et al. 2004) onboard the Spitzer Space Telescope (Werner et al. 2004) to improve the results. Useful data were obtained in the MIPS 24 and 70 μm bands, which have effective wavelengths of 23.68 and 71.42 μm. The MIPS 70 μm band overlaps with the Herschel 70 μm channel, providing a consistency check between data from the two observatories. A comparison of the PACS and MIPS 70 μm band flux densities shows that most measurements are consistent within a 1σ range. Significant deviations occur for 2003 AZ84 and Huya, which might have been caused by lightcurve effects or statistical noise.
The reduction of MIPS observations of TNOs has been described in detail in Stansberry et al. (2008) and Brucker et al. (2009), details on the calibration can be found in Gordon et al. (2007) and Engelbracht et al. (2007). We adopt absolute calibration uncertainties of 3% and 6% for the 24 and 70 μm observations (50% larger than the uncertainties derived for observations of stellar calibrators), respectively. The larger uncertainties account for effects from the skysubtraction technique, the faintness of the targets, and uncertainties in the color corrections.
The sky subtraction techniques introduced in Sect. 2.2 are derived from techniques originally developed for the MIPS data reductions. Reprocessed fluxes presented here are based on new reductions of the MIPS data, utilizing updated ephemeris positions for the targets. This allows for more precise masking of the target when generating the image of the sky background, and for more precise placement of the photometric aperture. In most cases the new fluxes are very similar to the previously published values for any given target, but in a few cases significant improvements in the measured flux density and SNR were achieved.
The standard color correction routine (Stansberry et al. 2007) resulting in monochromatic flux densities for the 24 and 70 μm bands requires the measurement of both bands and assumes the temperature of a black body which fits the 24:70 flux density ratio best. However, the two previously unpublished flux densities of 2002 VE95 and 2002 XV93 are for the 24 μm band only, which precludes the application of the standard method. In order to provide approximate color correction to these flux densities, we apply a method similar to the HerschelPACS color correction routine.
2.4. Optical Photometry
In order to derive diameter and albedo estimates, we combine thermal infrared measurements with optical data in the form of the absolute magnitude H, which is the object’s magnitude if it was observed at 1 AU heliocentric distance and 1 AU distance from the observer at a phase angle α = 0. H, albedo p_{V} and diameter d are related via (1)with the magnitude of the Sun V_{⊙} = −26.76 ± 0.02 (Bessell et al. 1998) in the JohnsonCousinsGlass system and δ a numerical constant. The relation returns d in km, if δ equals 1 AU expressed in km. In order to derive reliable albedo estimates, H and its uncertainty have to be well known (Harris & Harris 1997). This is not the case for most small bodies in the solar system, including Plutinos, for which reliable data are usually sparse. We determine H magnitudes from observed magnitudes in literature and observational data from the Minor Planet Center (MPC) (cf. Table 3). We deliberately do not use the H magnitudes provided by the MPC, since there is no uncertainty estimate given for these values and the reliability of these magnitudes is questionable (for instance, Romanishin & Tegler 2005, have determined a systematic uncertainty of MPC H magnitudes of 0.3 mag). Our results include uncertainty assessments and are reasonably close other estimates (for instance Doressoundiram et al. 2007).
H and the geometric albedo are defined at zero phase angle (α). However, existing photometric data were taken at α ≠ 0, which is corrected for by assuming a linear approximation of the phase angle dependence: H_{V} = m_{V}(1,1) − α·β, in which m_{V}(1,1) is the apparent V band magnitude normalized to unity heliocentric and geocentric distance, and β is the linear phase coefficient. Since we want to derive the geometric albedo p_{V}, we use V band photometry when available, otherwise we use R band photometry and convert to V band magnitudes using an averaged ⟨ V − R ⟩ = 0.567 ± 0.118 color index^{3}, and a correction for the intrinsic solar color index of (V − R)_{⊙} = 0.36 (Doressoundiram et al. 2008). ⟨ V − R ⟩ is based on observations of 41 Plutinos and is adopted from Hainaut & Delsanti (2002, data as of July 13, 2010).
If available, we give priority to data drawn from refereed publications, due to the existence of photometric uncertainties and the better calibration. All available apparent magnitudes for a target are converted to m_{V}(1,1) magnitudes and plotted as a function of the phase angle. For targets where the data show a clear trend of m_{V}(1,1) vs. α, we use a weighted linear regression analysis to fit a line to the data, yielding the linear phase coefficient. In the case of high scatter, a canonical β = 0.10 ± 0.04 mag/° is assumed, based on data given by Belskaya et al. (2008) (V and R band data, excluding Pluto). H then represents the average of the values determined from each m_{V}(1,1). MPC photometric data are usually highly scattered due to the coarse photometry, so the fixedβ technique is always applied if only MPC data are available. The results of all computations are given in Table 3, which also gives lightcurve amplitude information and color indices, as far as available. Due to the usage of a twovisit observation strategy, the determined Herschel flux densities are a combination of flux densities taken at two different points in time. Hence, it is not trivial to correct both Herschel measurements for lightcurve effects. However, lightcurve effects are already included in the uncertainty of the absolute magnitude H, assuming that the optical measurements are randomly sampled with respect to the individual lightcurve. In cases where the uncertainty of H is smaller than half the peaktopeak lightcurve amplitude, we account for lightcurve effects by a quadratic addition of the latter value to the uncertainty of H. This was only necessary in three cases (cf. Table 3).
3. Thermal modeling
By combining thermalinfrared and optical data, the physical properties of an object can be estimated using a thermal model. The diskintegrated thermal emission at wavelength λ of a spherical model asteroid is calculated from the surface temperature distribution T(θ,ϕ) as (2)where ϵ is the emissivity, d the object’s diameter, Δ its distance from the observer, α is the solar phase angle and θ and ϕ are surface latitude and longitude measured from the subsolar point, respectively. B(λ,T) is the Planck equation, B(λ,T) = 2hc^{2}/λ^{5} [exp(hc/(λk_{B}T)) − 1] ^{1}, where T is the temperature, h Planck’s constant, c the velocity of light and k_{B} Boltzmann’s constant. The diskintegrated thermal emission is then fitted to the thermalinfrared data by variation of the model surface temperature distribution T. The body’s temperature distribution depends on parameters such as albedo, thermal inertia, surface roughness, observation geometry, spin axis orientation and rotation period. In the modeling process, we take into account the individual observation geometries of each flux density measurement.
One of the first and most simple models is the standard thermal model (STM, cf. Lebofsky & Spencer 1989, and references therein), which assumes the model body to be of spherical shape, nonrotating and/or having zero thermal inertia, with a smooth surface and observed at phase angle α = 0. The surface temperature distribution is described by an instantaneous thermal equilibrium between emitted thermal radiation and absorbed sunlight.
In this work we make use of the nearEarth asteroid thermal model (NEATM, Harris 1998), which is, despite its name, applicable to atmosphereless TNOs (Müller et al. 2010; Stansberry et al. 2008, the latter use a similar approach with their hybrid STM). The main difference of the NEATM compared to the STM is the use of a variable beaming factor η, which adjusts the subsolar temperature (3)where A is the Bond albedo, S_{⊙} the solar constant, σ the StefanBoltzmann constant and r the object’s heliocentric distance. A is related to the geometric albedo, p_{V}, via A ~ A_{V} = qp_{V} (Lebofsky & Spencer 1989). We use a phase integral q = 0.336 p_{V} + 0.479 (Brucker et al. 2009). η accounts for thermal inertia and surface roughness in a firstorder approximation and can serve as a free parameter being derived as part of the fitting process, which usually provides better results than assuming a fixed value of η. Small values of η (η < 1) imply higher surface temperatures compared to that of a Lambertian surface, for instance due to surface roughness or porosity. High values of η (η > 1) lead to a reduction of the model surface temperature, mimicking the effect of thermal inertia. In contrast to the STM, the NEATM accounts for phase angles ≠ 0 (as shown in Eq. (2)). In the case of TNOs, however, this aspect is unimportant, since these objects are always observed at very low phase angles, due to their large heliocentric distances.
Whenever Spitzer data are available we combine them with PACS data. The combination of the data has a large impact on the reliability and accuracy of the model output. The peak emission wavelength in flux density units (Jy) of typical Plutino orbits ranges from 80 μm to 100 μm (compare model SEDs in Fig. 3). Hence, PACS band wavelengths (70 μm, 100 μm and 160 μm) are located in the flat peak plateau and the shallow RayleighJeans tail of the Planck function, respectively, and therefore constrain the model SED and η insufficiently. Adding a Spitzer 24 μm band measurement, which is located in the steeper Wienslope, improves the ability to constrain the model SED and allows for the determination of reliable estimates of η. Therefore, we adopt floatingη fits if additional Spitzer data are available and rely on floatingη fits if only Herschel data are available. However, in some cases of poor data quality or whenever data from different parts of the lightcurve or from very different observation geometries are combined, the floatingη method leads to η values which are too high or too low and therefore unphysical (the range of physically meaningful η values is 0.6 ≤ η ≤ 2.6^{4}). In such cases, we also apply the fixedη method, which makes use of a canonical value of η = 1.20 ± 0.35 derived by Stansberry et al. (2008) from Spitzer observations of a sample of Centaurs and TNOs. The validity of this approach is discussed in Sect. 4.1. We treat nondetections with a 1σ uncertainty as 0 ± σ detections in our models.
Throughout this work we assume the surface emissivity ϵ to be constant at a value of 0.9 for all wavelengths, which is based on laboratory measurements of silicate powder (Hovis & Callahan 1966) and a commonly adopted approximation for small bodies in the solar system.
3.1. Uncertainty assessment
In order to derive uncertainty estimates for diameter and albedo, we make use of a Monte Carlo (MC) simulation in which a sample of 1000 randomized synthetic bodies is created by variation of the observed flux density, the η value, and the absolute magnitude H. The randomized parameters follow a normal distribution centered on the nominal value, within the limits of the respective uncertainties. The uncertainties in flux density and H are taken from Tables 2 and 3, respectively. In the case of floatingη fits, η follows as a result of the modeling process, whereas in the case of fixedη fits, the value of η ± 0.35, determined by Stansberry et al. (2008) is applied.
Uncertainties of a specific parameter are derived from the ensemble of modeled synthetic bodies as the upper and lower values that include 68.2% of the ensemble (centered on the median), respectively. The uncertainties are usually aligned asymmetrically around the median, since most diameter and albedo distributions do not strictly follow a normal distribution. This method was introduced by Mueller et al. (2011). Finally, we combine the uncertainties determined through the MC method with the respective parameter values provided by the best model fit.
In order to improve uncertainty estimates, we make use of the rescaled uncertainties approach, introduced by SantosSanz et al. (2012). The motivation for this technique is the fact, that uncertainty estimates resulting from the MC method might be unrealistically small in cases in which the model fits the observational data only poorly. As a measure for the fit quality we use the reduced χ^{2}, which is a result of the fitting process. It is defined as (4)with F_{obs, i} and F_{model, i} being the observed and modeled flux densities of the different observations i, respectively, and σ_{obs, i} being the respective observational flux density uncertainty. is normalized by the number of degrees of freedom (n.d.f.), which equals the number of datapoints −2 in the case of a floatingη fit and the number of datapoints −1 in the case of a fixedη method, respectively. A poor fit results in .
Fig. 2 Illustration of the rescaled flux density uncertainties technique using synthetic data. a) Filled circles represent the flux densities at the different wavelengths with original flux density uncertainties. The dark line represents the best fit line (); the grey lines are fit lines of 100 out of 1000 synthetical bodies generated for the MC method. The bad quality of the fit is obvious, since the 70 and 100 μm datapoints are hardly fit by any of the single fit lines. Therefore, the reliability of the uncertainties of the output parameters is questionable. b) By rescaling the error bars by a factor of , the broadness of the distribution of single fit lines is increased, which covers now all datapoints. As a result, the output parameter fractional uncertainties are increased and therefore more realistic. 
Modeling results sorted by modeling technique.
The idea is to increase the flux density uncertainties uniformly by a factor S, which leads to an increase in variation of the flux densities of the synthetic bodies. The value of S is determined in such a way as to provide by definition the best possible model fit, i.e. , which is achieved by choosing , with being the minimum resulting from fitting the unaltered observational data. Hence, the rescaled flux density uncertainties are and the resulting , which is obtained by replacing σ_{obs} with in Eq. (4), is by definition unity. Diameter uncertainties are dominated by flux uncertainties and the quality of the fit, whereas albedo uncertainties are dominated by the uncertainty of H. The rescaled uncertainties are used only in the MC method to determine the uncertainties of the physical parameters. In order to obtain the best fit, the original uncertainties are applied.
The method is illustrated in Fig. 2 using synthetic data. Figure 2a shows the original flux densities with error bars of the size of the original uncertainties. The dark line represents the best NEATM fit (floatingη method) of the dataset, whereas each grey line symbolizes the best fit SED of one of 100 of the 1000 randomized synthetic bodies. It is obvious that the distribution of single fit curves misses the 100 μm PACS band completely and hardly fits the 70 μm datapoint. Hence, the fit quality is bad, which results in an underestimation of the model output uncertainties. By rescaling the uncertainties of the flux densities, the broadness of the distribution of single fit lines increases, as depicted in Fig. 2b. As a consequence of the rescaling, all datapoints are covered by the set of lines and the derived uncertainties are significantly increased, leading to much more realistic uncertainty estimates.
Fig. 3 Model fits of the Plutino sample using different modeling approaches and different data samples: the solid and dashed lines are floatingη fits based on combined Herschel and Spitzer and Herschelonly data, respectively; the dashdotted and dashdotdotted lines represent fixedη fits of combined Herschel and Spitzer and Herschelonly data samples, respectively; datapoints at 23.68 μm and 71.42 μm are Spitzer MIPS data (open circles), those at 70 μm, 100 μm and 160 μm are HerschelPACS data (filled circles). The model is fit to each flux measurement using the circumstances appropriate for the epoch of that observation. In order to simplify the figures, we have normalized the Spitzer flux densities to the epoch of the Herschel observations using the ratio of the bestfit model flux densities for the Herschel and Spitzer epoch, and then plotted the measured value times this ratio. Uncertainties are rescaled similarly, preserving the SNR for the Spitzer data. 
Fig. 3 continued. 
Fig. 3 continued. 
4. Results and discussion
An overview of all modeling results is given in Table 5. Where possible, floating and fixedη modeling approaches were applied based on Herschelonly and combined Herschel and Spitzer datasets, in order to test the consistency of the results. The objects’ SEDs determined by the different models are plotted in Fig. 3.
In the following sections we test the validity of the η = 1.2 assumption for the Plutino sample, and discuss the modeling results of the individual objects. After that, the results of the sample as a whole are discussed and used to determine a cumulative size distribution of the Plutino population. Finally, correlations between the different physical, photometric and dynamical parameters and the presence of ices are examined.
4.1. Test of the fixedη assumption
We were able to derive floatingη fits for 12 targets; we adopted the 7 which are based on combined Herschel and Spitzer data. The top section of Table 5 gives the results of all converged floatingη fits. The weighted mean^{5}η determined from the 7 adopted floatingη model fits yields , which is consistent with the adopted value of η = 1.20 ± 0.35 (Stansberry et al. 2008) applied in the fixedη approach. The root mean square fractional residuals in the results of the fixedη method between using η = 1.20 and ⟨ η ⟩ = 1.11 is 5% and 10% for diameter and albedo, respectively. However, the average fractional uncertainties in diameter and albedo of all adopted results are ± 10% and +49/−31%, respectively, which is larger than the residuals. Therefore, we continue to adopt η = 1.20 for the fixedη approach.
4.2. Discussion of individual targets
1996 TP66: Despite the availability of both Herschel and Spitzer data, the floatingη approach was not applicable for 1996 TP66, because of the nondetections in the PACS 70, 100 and 160 μm bands. The best fit fixedη solution is largely constrained by the two 24 μm MIPS flux densities. Therefore, our results ( km, , based on η = 1.20) barely differ from those of Stansberry et al. (2008), who find km and using , solely based on the 24 μm MIPS flux densities. The reason for the nondetections is unclear. However, the fact that no clear detection was possible in the wavelength regime of 70 μm and above using two different instruments (PACS and MIPS) and at different epochs points to a real effect. We based the observation planning on p_{V} = 0.07, which agrees well with the best fit model.
1999 TC36: This target was well detected in all three PACS bands. In combination with MIPS flux densities, a good floatingη fit was possible here, which misses the PACS 160 μm band flux density. This might hint to a more complex surface albedo distribution or a wavelengthdependency of the emissivity. Our adopted diameter and albedo estimates, km and , respectively, agree with earlier estimates by Stansberry et al. (2008): km and . 1999 TC36 is a triple component system (Trujillo & Brown 2002; Jacobson & Margot 2007; Benecchi et al. 2010), consisting of two similarly sized central components A_{1} and A_{2} and a more distant secondary component, B. The nature of this system makes it possible to determine the mean density of the system and, using optical flux differences between the single components, the sizes of the components. We revisit the calculations performed by Stansberry et al. (2006) and Benecchi et al. (2010). Our newly derived estimates of the sizes, based on the flux differences found by Benecchi et al. (2010), are: km, km and km. Using these diameters, we are able to improve the mean system density , which is somewhat higher than the earlier estimate (Benecchi et al. 2010). Assuming material densities of 1.0 ≤ ϱ_{0} ≤ 2.0 gcm^{3}, we find porosities or fractional void space of 36–68%. Our calculations are based on the assumptions of spherical shape of each component and the same albedo for each component.
2000 GN171: Both fixed and floatingη model fits based on the combined Herschel and Spitzer datasets were possible, leading to good fits with similar outcomes. We adopt the floatingη result km, , and , which differs significantly from an earlier estimate by Stansberry et al. (2008), who determined km, and . The large discrepancy stems from an overestimated MIPS 70 μm band flux density in the Stansberry et al. (2008) data, which has been revised for this work. By adding the HerschelPACS data, the SED is better constrained, leading to a more reliable result. 2000 GN171 is either a Roche binary or a Jacobi ellipsoid (Lacerda & Jewitt 2007). With a lightcurve amplitude of 0.61 mag (Sheppard & Jewitt 2002), the size of the companion or the ellipticity of the body is significant. Hence, the real diameter of 2000 GN171, the smallest Plutino in our sample, might be even smaller.
2001 KD77: We modeled this target solely on the basis of PACS data, since existing MIPS flux density measurements suffer from a large positional uncertainty. For this faint target, the fixedη fit mainly relies on the 70 μm band measurement, which has the highest SNR of all bands. Figure 3 clearly shows the mismatch of the 100 μm band measurement, which is scarcely detectable on the sky subtracted maps. Due to its large uncertainty, this datapoint is nearly neglected by the fixedη fit, and therefore scarcely impacts the modeling results. The reason for the low 100 μm flux density is unknown.
2002 VE95, 2002 XV93, 2003 AZ84, 2003 VS2 and Huya: These objects show a good to excellent agreement of both the different model technique fits compared to each other, as well as the model SED compared with the different measured flux densities. For all five targets combined Herschel and Spitzer flux densities are available and we adopt the floatingη modeling results. 2003 AZ84 and 2003 VS2 data include averaged lightcurve data, which strongly constrain the respective SEDs, due to their small uncertainties. Earlier estimates of 2003 VS2 by Stansberry et al. (2008) based on Spitzer data suggested km and using , which is significantly different from our results km, , based on . The discrepancy might stem from the high 70 μm MIPS flux density value, which, together with the 24 μm MIPS data point, suggests a steeper Wien slope of the spectral energy distribution compared to the combined Herschel and Spitzer data set. 2003 AZ84 was observed before by Herschel during the Science Demonstration Phase using the chop/nod technique (Müller et al. 2010), yielding d = 896 ± 55 and p_{V} = 0.065 ± 0.008, using η = 1.31 ± 0.08, which differs from our results km, , based on . Despite the smaller diameter determined using the averaged lightcurve data, 2003 AZ84 is still the largest Plutino in our sample, apart from Pluto. 2003 AZ84 has a moon, which is 5.0 ± 0.3 mag fainter than the primary (Brown & Suer 2007). This large difference in magnitude suggests that the thermal flux of the companion is negligible. The modeling results of Huya vary significantly depending on the model technique and the dataset. A visual inspection of the model fits shows that the floatingη approach based on the combined data set matches the measured SED best and results in and , based on . Earlier estimates by Stansberry et al. (2008) based on Spitzer data suggest and , based on . 2002 XV93 has the lowest albedo in the sample.
2001 QF298, 2002 VR128, 2002 VU130, 2003 UT292, 2004 EW95, 2004 PF115, 2004 UX10 and 2006 HJ123: For these targets we have to rely on Herschelonly data and therefore adopt the results of the fixedη technique. We were able to succesfully apply the floatingη method to 2004 EW95, 2004 PF115, 2004 UX10 and 2006 HJ123 as well. Modeling results in Table 5 and the plots in Fig. 3 of the SEDs show a good agreement between the results of the different techniques, supporting the validity of the fixedη approach and the reliability of the Herschel flux measurements. Differences in the fit quality are induced by the different SNR values of the respective measurements. Both 2002 VU130 and 2006 HJ123 suffer from large albedo uncertainties, which is a result of the large uncertainty of their absolute magnitude H. 2006 HJ123 shows the highest albedo of our sample, excluding Pluto.
Pluto: Paradoxically, Pluto is a clear outlier from the Plutino population, being by far the largest and having the brightest surface. Besides that, Pluto was the first TNO discovered to have an atmosphere (Hubbard et al. 1988; Brosch 1995), has a pronounced optical lightcurve caused by albedo variations, and has a moon half its own size. An additional complication is the effect of N_{2} ice on Pluto’s surface: sublimation and deposition of N_{2} is accompanied by latent heat transport that results in much lower (higher) dayside (nightside) temperatures than would be expected on an airless body (e.g. Spencer et al. 1997). Despite these facts, we attempt to treat Pluto like any other Plutino and compare our modeling results to the known properties of Pluto and Charon, and discuss the reasons for discrepancies that emerge.
The Pluto system has been subject to detailed study via, for instance, direct imaging, stellar occultation observations and thermal radiometry in the past. Mean hemispheric albedo ranges are 0.49 to 0.66 for Pluto and 0.36 to 0.39 for Charon (Buie et al. 1997). Charon’s diameter is well determined to be 1208 km (Gulbis et al. 2006; Sicardy et al. 2006). Pluto’s diameter estimates suffer from uncertainties related to the presence of the atmosphere and range from about 2290–2400 km (Tholen & Buie 1997). For the purpose of comparing with the results of our thermal models, we adopt a diameter of 2322 km (Young et al. 2007). These numbers lead to an equivalent diameter and albedo of the combined system of 2617 km and 0.44–0.73, respectively.
We deliberately refrain from combining our Herschel data set with existing Spitzer measurements (for instance Lellouch et al. 2011), due to the so far unexplained secular fading of the Pluto/Charon system (Stansberry et al. 2009). Instead, we rely on our original Herschel measurements only, and so we adopt the fixedη results: km and (cf. Fig. 3). Our thermal estimate of the effective system diameter deviates by 3σ from the known effective system diameter given above; the albedo derived from our thermal modeling is only barely consistent with the maximum allowable albedo given above. Our attempt to apply a simple thermal model to Pluto shows the susceptibility of such models to violations of the model assumptions. We discuss effects, which may explain the discrepancy of our adopted values with the equivalent values:

By replacing the Pluto/Charon system with a single equivalent body, we had to assume a single average albedo for both objects, which is a coarse simplification, in order to apply our models. Furthermore, surface albedo variations of Pluto (Buie et al. 2010) lead to a more complex spectral energy distribution and affect the model results significantly. Our model results would be improved by using a twoterrain model, which was used to model the thermal emission of Makemake (Lim et al. 2010), in which Pluto and Charon would each be represented by one terrain, or applying more sophisticated albedo surface distributions: Lellouch et al. (2011) applied a 3terrain model to model the thermal emission of Pluto and an additional terrain to model Charon’s contribution.

Former thermal infrared measurements of Pluto revealed its thermal lightcurve in different wavelength ranges (Lellouch et al. 2000, 2011), which was not accounted for in our modeling. Given the time lags between the individual observations, this might have introduced a perceptible distortion of the SED, which affects the modeling results.

Millimeter wavelength observations of Pluto (Gurwell & Butler 2005; Gurwell et al. 2010) revealed a lower surface brightness temperature than assumed from thermal equilibrium, which is due to the combination of two effects: (1) the mixing of SEDs at different temperatures, and in particular the increased contribution of the N_{2}dominated areas, whose temperature is ~37 K due to the dominance of sublimation cooling in the thermal budget of N_{2} ice, and (2) the fact that emissivities are expected to decrease at long wavelengths as a result of subsurface sounding in a progressively more transparent medium. Lellouch et al. (2011) find in particular that the emissivity of CH_{4} ice, one of the terrains in their 3terrain Pluto model, shows a steady decrease with increasing wavelength (their Fig. 8). Without additional assumptions, our thermal model is not capable of taking such effects into account.
Fig. 4 Diameter and albedo probability density functions and histograms. The probability density functions describe the probability to find an object of a specific diameter or albedo, when picking a random object. The functions were derived by summing up the probability density functions of each object, being represented by a nonsymmetric normal distribution taking into account the different upper and lower uncertainties of each result. This approach returns a more detailed picture than the histogram. Pluto was excluded from the sample. 
4.3. Sample statistics
Figure 4 shows diameter and albedo probability density functions and histograms based on the results compiled in Table 6. The size distribution ranges from 150 km to 730 km and the albedo distribution from 0.04 to 0.28 with a clear peak around 0.08. The low probability density at high albedos is caused by the large uncertainties of such albedos. Excluding Pluto, the weighted mean^{6} of the albedo, weighted by absolute uncertainties, yields 0.08 ± 0.03, which agrees well with the range of typical geometric albedos of the TNO and Centaur sample of Stansberry et al. (2008). Table 7 shows a comparison of the mean albedos of our Plutino sample with different TNO subpopulation and JFCs. It turns out that the average albedo of the Plutinos agrees well with that of the scattered disk, Centaur, hot classical and JFC population. Through dynamical studies Duncan et al. (1995), Yu & Tremaine (1999) and di Sisto et al. (2010) suggested the Plutino population to be a source of JFCs and/or Centaurs, which is compatible with our results. Furthermore, most of the Plutinos appear to have darker surfaces than the cold classical TNOs and detached objects.
Adopted modeling results.
Average albedos of Centaurs, different TNO subpopulations and JFCs.
4.4. Plutino size distribution
Using the absolute magnitude estimates for all known Plutinos from the MPC and our measured albedos we are able to calibrate the Plutino size scale for the first time and to determine a cumulative size distribution of the known Plutino population. We use two different approaches to determine the cumulative size distribution N(≥ d). Firstly, we assume monochromatic albedo distributions based on the two average albedos determined in Sect. 4.3. This approach enables the direct conversion of H magnitudes into diameter but neglects the measured diversity in albedo. The second approach makes use of a Monte Carlo method to determine the size distribution based on the actually measured albedo distribution. For this purpose, each MPC H magnitude is assigned to a randomly generated albedo, which is used to determine the diameter. Hence, we assume that H is not correlated to the albedo, which is supported by our findings in Sect. 4.5. The distribution of the randomly generated albedos follows the determined probability density function (Fig. 4). Hence, this approach takes into account the measured albedo diversity. Figure 5 shows a comparison of the two different methods. It is clearly visible that the averaged Monte Carlo approach agrees well with the monochromatic p_{V} = 0.08 model, which supports our assumption that this value better represents the average Plutino albedo compared to the value of 0.11.
Fig. 5 Cumulative Plutino size distributions N(≥ d) based on MPC H magnitudes and different assumptions: the solid and dashed line represent monochromatic albedo distributions using p_{V} = 0.08 and p_{V} = 0.11, respectively. The grey line shows the result of an averaged Monte Carlo simulations based on 100 runs and the measured albedo distribution with corresponding 1σ standard deviations. It is obvious that the monochromatic p_{V} = 0.08 size distribution agrees better with the result of the Monte Carlo simulation. Furthermore, we show two cumulative power laws N(≥d) ~ d^{−q} with q = 2 and q = 3 and arbitrary interceptions for comparison. A q = 2 cumulative power law better fits the measured size distribution in the size range 120–450 km, whereas the distribution of larger objects is better fit by a q = 3 cumulative power law. We can give an upper size limit for the kink at the smallsized end around 120 km. The upper abscissa gives the equivalent H magnitude based on a monochromatic albedo distribution of p_{V} = 0.08. 
For comparison, we plot in Fig. 5 the slopes of different cumulative power laws N(≥ d) ~ d^{−q}. We find that the distribution of intermediatesized Plutinos (120–450 km) is welldescribed by q = 2, which is smaller than the presumed value of q ~ 3 (Trujillo et al. 2001; Kenyon & Bromley 2004). The large diameter tail (d ≥ 450 km) seems to be better described by q = 3. However, this region suffers from uncertainties due to small number statistics. A change in the cumulative power law slope at the smallsized tail seems to occur at diameters of 80−120 km, which is larger than the proposed range of 40−80 km (Bernstein et al. 2004; Kenyon et al. 2008; de Elía et al. 2008). However, the number of intrinsically fainter, and therefore smaller Plutinos is likely to be underrepresented, since the MPC sample suffers from an optical discovery bias. A proper debiasing of the size and albedo distributions is planned as part of future work. However, we suppose that a larger number of smallsized Plutinos would shift the kink to smaller diameters. Hence, we are only able to give an upper size limit of the location of the kink, which is 120 km. Furthermore, we note that our size distribution is based on albedo measurements derived from objects of the size range 150−730 km, which might be different from the albedo distribution of smaller Plutinos.
Fig. 6 Plots of some correlations. Upper left: we find no correlation between diameter and albedo; upper right: a clear anticorrelation between eccentricity and diameter is the result of a detection bias; lower left: we find smaller Plutinos to be redder than larger ones, which contradicts what was previously suggested by Peixinho et al. (2004) based on the intrinsical brightness H of a different Plutino sample; lower right: we find that the color of Plutinos is correlated to the heliocentric distance at the time of their discovery, which hints to a color bias. 
4.5. Correlations
The Plutino sample was checked for correlations of the physical parameters, orbital properties and intrinsic optical colors, using a modified Spearman rank correlation analysis technique. The technique is based on a bootstrap method and takes into account the uncertainties (error bars) of the parameters (data points), and computes the confidence interval of the corresponding correlation coefficient. Since we are taking into account data error bars, the results are estimates of the most probable correlation coefficient ⟨ ρ ⟩ accompanied by the probability P (pvalue) of that ⟨ ρ ⟩ value to occur given no relationship in the sample, i.e. to occur by pure chance. The lower the value of P the higher the significance level of the correlation. We present the 68% confidence interval for the correlation value ⟨ ρ ⟩ , i.e. the interval which includes 68% of the bootstrapped values. This confidence interval equals the canonical 1σ interval of a Gaussian distribution. As an additional aid in reading the P values we also compute its Gaussianσ equivalent significance P_{σ}. The details of the technique are described in SantosSanz et al. (2012). Following Efron & Tibshirani (1993), we define a weak correlation by 0.3 <⟨ ρ ⟩< 0.6 and a strong correlation by  ⟨ ρ ⟩> 0.6. P < 0.05 indicates reasonably strong evidence of correlation, P < 0.025 indicates strong evidence of correlation and P < 0.003 describes a clear (3σ) correlation.
It turns out that most parameters are uncorrelated. Only very few parameter pairs show reasonably strong evidence of correlation. Those and a few other interesting results are listed in Table 8 and discussed here. Other parameter pairs lack sufficient evidence for correlation or possess weak correlation values.

We do not detect any strong and significant trend between diameter d and albedo p_{V} (cf. Fig. 6). Although we cannot rule out the possibility of some correlation between these two parameters given our sample size, it is very unlikely that a strong correlation between the two parameters might exist. On the other hand, both diameter and albedo seem to be clearly uncorrelated with the η value.

The absolute magnitude H is not correlated with diameter, albedo, or η. This shows that our sample is not biased towards large and high (or low) albedo objects, which supports our approach to randomly assign albedos to H magnitudes in Sect. 4.4. This finding, however, does not exclude a discovery bias on smaller diameters than those examined here or a discovery bias based on geometrical aspects.

We find a significant anticorrelation between eccentricity e and diameter d (cf. Fig. 6), showing that in our sample highly eccentric objects tend to be smaller. This anticorrelation is precisely the opposite of the correlation between perihelion distance q and diameter d, and also the opposite of the correlation between diameter d and the heliocentric distance of each object at the time of its discovery r^{⋆}. We consider that these relationships are merely caused by the fact that objects on more eccentric orbits come closer to the Sun, which significantly improved their detectability. Hence, the observed relation between an object’s diameter and its eccentricity (and perihelion distance) is very likely to be caused by a discovery bias.

The correlation between diameter d and the heliocentric distance of each object at the time of its discovery r^{⋆}, compared to the lack of correlation between albedo p_{V} and r^{⋆}, shows that the likelihood of detection solely depends on size, and is rather independent of the object’s albedo. This follows because the brightness of an object scales as d^{2}p_{V}: the measured range of diameters translates into a brightness change factor of 24, whereas the total effect of the albedo can only account for a factor of 6. We can rule out a bias towards high (and low) albedos as a result of the nature of discovery of our sample targets. This supports the representativeness of the measured albedo distribution for the Plutino population and shows that the nature of the discovery bias is sizedependent.

Most color indices and the spectral slope are anticorrelated with heliocentric distance at the time of an object’s discovery r^{⋆} (Fig. 6). The farther an object is, the less red it seems to be. This points to a probable color bias in our sample, maybe induced by the common use of R band filters and the improved sensitivity in red bands of stateoftheart detectors, which are used in TNO surveys. However, we find a lack of bluer objects at shorter heliocentric distances, which should be detectable, despite their color. This might suggest that such objects do not exist, at least within our sample.
Table 8Results of the modified Spearman rank correlation analysis of the Plutino sample.

We can not confirm a trend found by Hainaut & Delsanti (2002) and Peixinho et al. (2004) of bluer Plutinos being intrinsically fainter using the spectral slope s as a measure of color. The correlation analysis of B − R and H leads to similar results: and P = 0.420. Our sample suggests precisely the opposite, with a trend between spectral slope s and the diameter d showing smaller Plutinos to be redder (Fig. 6). We note however, that Peixinho et al. (2004) find their trend due to a “cluster” of blue Plutinos with H_{R} > 7.5, i.e. H ≳ 8, objects at magnitude/size ranges we do not have in our sample. Hence, both effects might be the result of selection effects.
5. Summary
The analysis of the diameters and albedos of 18 Plutinos using PACS photometry, leads us to the following conclusions:

The diameter and albedo range of our Plutino sample yields 150–730 km and 0.04–0.28, respectively. Excluding Pluto, the weighted mean of the Plutino albedo distribution is 0.08 ± 0.03 and agrees with the average albedos observed for the scattered disk, Centaur and JFC population. This agreement is compatible with the idea that Plutinos are a source population of Centaurs and JFCs.

The floatingη fits yield a weighted mean η for the Plutino sample of , which agrees with the canonical η = 1.20 ± 0.35 (Stansberry et al. 2008) utilized in the fixedη fits.

We refine the size estimates of the components of 1999 TC36: km, km and km; we estimate the bulk density of the multiple system to be .

Pluto is the clear outlier of the Plutino population, being by far the largest object with the brightest surface. We have shown, that canonical simple thermal modeling of Pluto data leads to inadequate results, probably mainly due to its multicomponent nature, thermal lightcurve and atmosphere.

Using our measured albedos we calibrated the Plutino size scale for the first time and determined a cumulative size distribution of the known Plutinos. We find that intermediate sized Plutinos (120 km ≤ D ≤ 400 km) follow a cumulative powerlaw distribution with q = 2, whereas the distribution of larger objects is better described by q = 3. We are able to give an upper size limit for the kink at the smallsized end of the distribution, which is 120 km.

We find no correlations between albedo and diameter, as well as H, and diameter, albedo and η, respectively. This shows that our sample is not biased towards large and high (or low) albedo objects. Furthermore, we find a correlation between diameter and the heliocentric distance at the time of discovery, but not between the albedo and the heliocentric distance at the time of discovery. This leads us to the conclusion that the nature of the discovery bias is mainly sizedependent.

A significant correlation between diameter and eccentricity (and perihelion distance) is very likely to be caused by a detectional bias based on geometrical aspects.

We find hints to color biases in our sample: Plutinos, which have been farther from Sun at the time of their discovery seem to be bluer and smaller sized Plutinos tend to be redder, which contradicts previous finding by Peixinho et al. (2004).

There is qualitative evidence that icy Plutinos have higher albedos than the average of the sample. We are not able to conclude on a correlation between the diameter and the presence of ices.
Statistical data are based on Minor Planet Center data as of 18 Nov. 2011 (http://minorplanetcenter.net/iau/MPCORB.html).
The applied Plutino V − R color agrees within 1σ with our sample mean of 0.51 ± 0.12 (cf. Table 3).
The physical range of η values were probed by using NEATM to model the two extreme cases of (1) a fast rotating object of high thermal inertia with low surface roughness and (2) a slow rotating object of low thermal inertia with high surface roughness, which yield the upper and lower η limit, respectively. Fluxes for both cases were determined using a full thermophysical model (Lagerros 1996, 1997, 1998; Mueller 2007), in which surface roughness is implemented via spherical cratering. The degree of surface roughness is determined via the opening angle (0°–180°) and surface density (0.0–1.0) of the craters.
The weighted mean and its uncertainties were determined using the equation given in Table 3.
Upper and lower uncertainties were calculated corresponding to the calculations applied in Table 3. The calculation of the weighted mean albedo based on fractional uncertainties yields 0.11. However, we prefer to adopt the weighted mean albedo based on absolute errors, since it better agrees with the low albedo of the bulk of the Plutinos and is enforced from our results in Sect. 4.4.
Acknowledgments
We would like to thank Alain Doressoundiram for his useful suggestions as a referee. Furthermore, we thank Chemeda Ejeta (MPS, KatlenburgLindau, Germany) for performing dynamical simulations on some of our targets. M. Mommert acknowledges support through the DFG Special Priority Program 1385: The First 10 Million Years of the solar system – a Planetary Materials Approach. C. Kiss and A. Pal acknowledge the support of the Bolyai Research Fellowship of the Hungarian Academy of Sciences. P. SantosSanz would like to acknowledge financial support by the Centre National de la Recherche Scientifique (CNRS). J. Stansberry acknowledges support for this work provided by NASA through an award issued by JPL/Caltech. R. Duffard acknowledges financial support from the MICINN (contract Ramón y Cajal). Part of this work was supported by the German Deutsches Zentrum für Luft und Raumfahrt (DLR) project numbers 50 OR 0903, 50 OFO 0903 and 50 OR 1108, and the PECS program of the European Space Agency (ESA) and the Hungarian Space Office, PECS98073. J.L. Ortiz acknowledges support from spanish grants AYA200806202C0301, AYA201130106C0201 and 2007FQM2998.
References
 AlvarezCandal, A., Barucci, M. A., Merlin, F., Guilbert, A., & de Bergh, C. 2007, A&A, 475, 369 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barkume, K. M., Brown, M. E., & Schaller, E. L. 2008, AJ, 135, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Barucci, M. A., Doressoundiram, A., Tholen, D., Fulchignoni, M., & Lazzarin, M. 1999, Icarus, 142, 476 [NASA ADS] [CrossRef] [Google Scholar]
 Barucci, M. A., Merlin, F., Dotto, E., Doressoundiram, A., & de Bergh, C. 2006, A&A, 455, 725 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barucci, M. A., AlvarezCandal, A., Merlin, F., et al. 2011, Icarus, 214, 297 [NASA ADS] [CrossRef] [Google Scholar]
 Barucci, M. A., Merlin, F., Perna, D., et al. 2012, A&A, 539, A152 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belskaya, I. N., LevasseurRegourd, A., Shkuratov, Y. G., & Muinonen, K. 2008, Surface Properties of Kuiper Belt Objects and Centaurs from Photometry and Polarimetry, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, A. Morbidelli, & R. Dotson, 115 [Google Scholar]
 Benecchi, S. D., Noll, K. S., Grundy, W. M., et al. 2009, Icarus, 200, 292 [NASA ADS] [CrossRef] [Google Scholar]
 Benecchi, S. D., Noll, K. S., Grundy, W. M., & Levison, H. F. 2010, Icarus, 207, 978 [NASA ADS] [CrossRef] [Google Scholar]
 Bernstein, G. M., Trilling, D. E., Allen, R. L., et al. 2004, AJ, 128, 1364 [NASA ADS] [CrossRef] [Google Scholar]
 Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
 Boehnhardt, H., Tozzi, G. P., Birkle, K., et al. 2001, A&A, 378, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brosch, N. 1995, MNRAS, 276, 571 [NASA ADS] [Google Scholar]
 Brown, M. E. 2000, AJ, 119, 977 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, M. E., & Suer, T.A. 2007, IAU Circ., 8812, 1 [NASA ADS] [Google Scholar]
 Brucker, M. J., Grundy, W. M., Stansberry, J. A., et al. 2009, Icarus, 201, 284 [NASA ADS] [CrossRef] [Google Scholar]
 Buie, M. W., Tholen, D. J., & Wasserman, L. H. 1997, Icarus, 125, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Buie, M. W., Grundy, W. M., Young, E. F., Young, L. A., & Stern, S. A. 2010, ApJ, 139, 1128 [Google Scholar]
 Cohen, C. J., & Hubbard, E. C. 1965, AJ, 70, 10 [NASA ADS] [CrossRef] [Google Scholar]
 CollanderBrown, S. J., Fitzsimmons, A., Fletcher, E., Irwin, M. J., & Williams, I. P. 1999, MNRAS, 308, 588 [NASA ADS] [CrossRef] [Google Scholar]
 Davies, J. K., Green, S., McBride, N., et al. 2000, Icarus, 146, 253 [NASA ADS] [CrossRef] [Google Scholar]
 de Bergh, C., Boehnhardt, H., Barucci, M. A., et al. 2004, A&A, 416, 791 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 de Elía, G. C., Brunini, A., & di Sisto, R. P. 2008, A&A, 490, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Delsanti, A. C., Boehnhardt, H., Barrera, L., et al. 2001, A&A, 380, 347 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 DeMeo, F. E., Fornasier, S., Barucci, M. A., et al. 2009, A&A, 493, 283 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 DeMeo, F. E., Dumas, C., de Bergh, C., et al. 2010, Icarus, 208, 412 [NASA ADS] [CrossRef] [Google Scholar]
 di Sisto, R. P., Brunini, A., & de Elía, G. C. 2010, A&A, 519, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Doressoundiram, A., Barucci, M. A., Romon, J., & Veillet, C. 2001, Icarus, 154, 277 [NASA ADS] [CrossRef] [Google Scholar]
 Doressoundiram, A., Peixinho, N., de Bergh, C., et al. 2002, AJ, 124, 2279 [NASA ADS] [CrossRef] [Google Scholar]
 Doressoundiram, A., Peixinho, N., Doucet, C., et al. 2005, Icarus, 174, 90 [NASA ADS] [CrossRef] [Google Scholar]
 Doressoundiram, A., Peixinho, N., Moullet, A., et al. 2007, AJ, 134, 2186 [NASA ADS] [CrossRef] [Google Scholar]
 Doressoundiram, A., Boehnhardt, H., Tegler, S. C., & Trujillo, C. 2008, Color Properties and Trends of the Transneptunian Objects, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, A. Morbidelli, & R. Dotson, 91 [Google Scholar]
 Dotto, E., Barucci, M. A., Boehnhardt, H., et al. 2003, Icarus, 162, 408 [NASA ADS] [CrossRef] [Google Scholar]
 Duncan, M. J., Levison, H. F., & Budd, S. M. 1995, AJ, 110, 3073 [NASA ADS] [CrossRef] [Google Scholar]
 Edgeworth, K. E. 1949, MNRAS, 109, 600 [NASA ADS] [CrossRef] [Google Scholar]
 Efron, B., & Tibshirani, R. 1993, An introduction to the bootstrap, Monographs on Statistics and Applied Probability (New York: Chapman and Hall/CRC), 57 [Google Scholar]
 Elliot, J. L., Kern, S. D., Clancy, K. B., et al. 2005, ApJ, 129, 1117 [Google Scholar]
 Engelbracht, C. W., Blaylock, M., Su, K. Y. L., et al. 2007, PASP, 119, 994 [NASA ADS] [CrossRef] [Google Scholar]
 Fernández, Y. R., Jewitt, D. C., & Sheppard, S. S. 2005, AJ, 130, 308 [NASA ADS] [CrossRef] [Google Scholar]
 Ferrin, I., Rabinowitz, D., Schaefer, B., et al. 2001, ApJ, 548, L243 [Google Scholar]
 Fornasier, S., Doressoundiram, A., Tozzi, G. P., et al. 2004, A&A, 421, 353 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gladman, B., Marsden, B. G., & Vanlaerhoven, C. 2008, in Nomenclature in the Outer Solar System, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli, 43 [Google Scholar]
 Gordon, K. D., Engelbracht, C. W., Fadda, D., et al. 2007, PASP, 119, 1019 [NASA ADS] [CrossRef] [Google Scholar]
 Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3 [Google Scholar]
 Guilbert, A., AlvarezCandal, A., Merlin, F., et al. 2009, Icarus, 201, 272 [NASA ADS] [CrossRef] [Google Scholar]
 Gulbis, A. A. S., Elliot, J. L., Person, M. J., et al. 2006, Nature, 439, 48 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Gurwell, M. A., & Butler, B. J. 2005, in AAS/Division for Planetary Sciences Meeting Abstracts, BAAS, 37, 743 [Google Scholar]
 Gurwell, M. A., Butler, B. J., & Moullet, A. 2010, in AAS/Division for Planetary Sciences Meeting Abstracts, BAAS, 42, 1014 [Google Scholar]
 Hainaut, O. R., & Delsanti, A. C. 2002, A&A, 389, 641 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Harris, A. W. 1998, Icarus, 131, 291 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Harris, A. W., & Harris, A. W. 1997, Icarus, 126, 450 [NASA ADS] [CrossRef] [Google Scholar]
 Hovis, W. A., & Callahan, W. R. 1966, J. Opt. Soc. Amer., 56, 639 [Google Scholar]
 Howell, S. B. 1989, PASP, 101, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Hubbard, W. B., Hunten, D. M., Dieters, S. W., Hill, K. M., & Watson, R. D. 1988, Nature, 336, 452 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jacobson, S., & Margot, J. L. 2007, in AAS/Division for Planetary Sciences Meeting Abstracts, BAAS, 38, 519 [Google Scholar]
 Jewitt, D., & Luu, J. 1993, Nature, 362, 730 [NASA ADS] [CrossRef] [Google Scholar]
 Jewitt, D., & Luu, J. 1996, in Completing the Inventory of the Solar System, ed. T. Rettig, & J. M. Hahn, ASP Conf. Ser., 107, 255 [Google Scholar]
 Jewitt, D., & Luu, J. 1998, AJ, 115, 1667 [NASA ADS] [CrossRef] [Google Scholar]
 Jewitt, D. C., & Luu, J. X. 2001, AJ, 122, 2099 [NASA ADS] [CrossRef] [Google Scholar]
 Kenyon, S. J., & Bromley, B. C. 2004, AJ, 128, 1916 [NASA ADS] [CrossRef] [Google Scholar]
 Kenyon, S. J., Bromley, B. C., O’Brien, D. P., & Davis, D. R. 2008, in Formation and Collisional Evolution of Kuiper Belt Objects, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, A. Morbidelli, & R. Dotson, 293 [Google Scholar]
 Kuiper, G. P. 1951, Proc. National Academy of Science, 37, 1 [Google Scholar]
 Lacerda, P., & Jewitt, D. C. 2007, ApJ, 133, 1393 [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]
 Lebofsky, L. A., & Spencer, J. R. 1989, in Asteroids II, ed. R. P. Binzel, T. Gehrels, & M. S. Matthews, 128 [Google Scholar]
 Lellouch, E., Laureijs, R., Schmitt, B., et al. 2000, Icarus, 147, 220 [NASA ADS] [CrossRef] [Google Scholar]
 Lellouch, E., Stansberry, J., Emery, J., Grundy, W., & Cruikshank, D. P. 2011, Icarus, 214, 701 [NASA ADS] [CrossRef] [Google Scholar]
 Licandro, J., Oliva, E., & Di Martino, M. 2001, A&A, 373, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lim, T. L., Stansberry, J., Müller, T. G., et al. 2010, A&A, 518, L148 [Google Scholar]
 Malhotra, R. 1993, Nature, 365, 819 [CrossRef] [Google Scholar]
 Malhotra, R. 1995, AJ, 110, 420 [NASA ADS] [CrossRef] [Google Scholar]
 McBride, N., Green, S. F., Davies, J. K., et al. 2003, Icarus, 161, 501 [NASA ADS] [CrossRef] [Google Scholar]
 Melita, M. D., & Brunini, A. 2000, Icarus, 147, 205 [NASA ADS] [CrossRef] [Google Scholar]
 Merlin, F., Barucci, M. A., Dotto, E., de Bergh, C., & Lo Curto, G. 2005, A&A, 444, 977 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Merlin, F., Barucci, M. A., de Bergh, C., et al. 2010, Icarus, 210, 930 [NASA ADS] [CrossRef] [Google Scholar]
 Morbidelli, A. 1997, Icarus, 127, 1 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Mueller, M. 2007, Ph.D. Thesis, FU Berlin, Germany [Google Scholar]
 Mueller, M., Delbo’, M., Hora, J. L., et al. 2011, AJ, 141, 109 [NASA ADS] [CrossRef] [Google Scholar]
 Müller, T. G., Lellouch, E., Böhnhardt, H., et al. 2009, Earth Moon and Planets, 105, 209 [CrossRef] [Google Scholar]
 Müller, T. G., Lellouch, E., Stansberry, J., et al. 2010, A&A, 518, L146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Müller, T. G., Nielbock, M., Balog, Z., Klaas, U., & Vilenius, E. 2011a, PACS Photometer – PointSource Flux Calibration, Document: PICCMETN037, Tech. Rep. [Google Scholar]
 Müller, T. G., Okumura, K., & Klaas, U. 2011b, PACS Photometer Passbands and Colour Correction Factors for Various Source SEDs, Document: PICCMETN038, Tech. Rep. [Google Scholar]
 Ortiz, J. L., Gutiérrez, P. J., Casanova, V., & Sota, A. 2003, A&A, 407, 1149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ortiz, J. L., Gutiérrez, P. J., SantosSanz, P., Casanova, V., & Sota, A. 2006, A&A, 447, 1131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ott, S. 2010, in Astronomical Data Analysis Software and Systems XIX, ed. Y. Mizumoto, K.I. Morita, & M. Ohishi, ASP Conf. Ser., 434, 139 [Google Scholar]
 Owen, T. C., Roush, T. L., Cruikshank, D. P., et al. 1993, Science, 261, 745 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Peixinho, N., Boehnhardt, H., Belskaya, I., et al. 2004, Icarus, 170, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Perna, D., Barucci, M. A., Fornasier, S., et al. 2010, A&A, 510, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
 Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Protopapa, S., AlvarezCandal, A., Barucci, M. A., et al. 2009, A&A, 501, 375 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rabinowitz, D. L., Schaefer, B. E., & Tourtellotte, S. W. 2007, AJ, 133, 26 [NASA ADS] [CrossRef] [Google Scholar]
 Rabinowitz, D. L., Schaefer, B. E., Schaefer, M., & Tourtellotte, S. W. 2008, AJ, 136, 1502 [NASA ADS] [CrossRef] [Google Scholar]
 Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, ApJS, 154, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Romanishin, W., & Tegler, S. C. 2005, Icarus, 179, 523 [NASA ADS] [CrossRef] [Google Scholar]
 Romanishin, W., Tegler, S. C., & Consolmagno, G. J. 2010, AJ, 140, 29 [NASA ADS] [CrossRef] [Google Scholar]
 SantosSanz, P., Ortiz, J. L., Barrera, L., & Boehnhardt, H. 2009, A&A, 494, 693 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 SantosSanz, P., Lellouch, E., Fornasier, S., et al. 2012, A&A, 541, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sheppard, S. S., & Jewitt, D. C. 2002, AJ, 124, 1757 [NASA ADS] [CrossRef] [Google Scholar]
 Sheppard, S. S., & Jewitt, D. C. 2003, Earth Moon and Planets, 92, 207 [Google Scholar]
 Sicardy, B., Bellucci, A., Gendron, E., et al. 2006, Nature, 439, 52 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Spencer, J. R., Stansberry, J. A., Trafton, L. M., et al. 1997, in Volatile Transport, Seasonal Cycles, and Atmospheric Dynamics on Pluto, ed. S. A. Stern, & D. J. Tholen, 435 [Google Scholar]
 Stansberry, J. A., Grundy, W. M., Margot, J. L., et al. 2006, ApJ, 643, 556 [NASA ADS] [CrossRef] [Google Scholar]
 Stansberry, J. A., Gordon, K. D., Bhattacharya, B., et al. 2007, PASP, 119, 1038 [NASA ADS] [CrossRef] [Google Scholar]
 Stansberry, J., Grundy, W., Brown, M., et al. 2008, in Physical Properties of Kuiper Belt and Centaur Objects: Constraints from the Spitzer Space Telescope, ed. M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli, 161 [Google Scholar]
 Stansberry, J. A., Lellouch, E., Grundy, W., et al. 2009, in AAS/Division for Planetary Sciences Meeting Abstracts, 41, 06.02 [Google Scholar]
 Tegler, S. C., & Romanishin, W. 1998, Nature, 392, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Tegler, S. C., Romanishin, W., & Consolmagno, G. J. 2003, ApJ, 599, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Thirouin, A., Ortiz, J. L., Duffard, R., et al. 2010, A&A, 522, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tholen, D. J., & Buie, M. W. 1997, in Bulk Properties of Pluto and Charon, ed. S. A. Stern, & D. J. Tholen, 193 [Google Scholar]
 Trujillo, C. A., & Brown, M. E. 2002, IAU Circ., 7787, 1 [NASA ADS] [Google Scholar]
 Trujillo, C. A., Luu, J. X., Bosh, A. S., & Elliot, J. L. 2001, AJ, 122, 2740 [NASA ADS] [CrossRef] [Google Scholar]
 Vilenius, E., Kiss, C., Mommert, M., et al. 2012, A&A, submitted [Google Scholar]
 Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [Google Scholar]
 Young, E. F., Young, L. A., & Buie, M. 2007, in BAAS, 38, 541 [Google Scholar]
 Yu, Q., & Tremaine, S. 1999, AJ, 118, 1873 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Absolute magnitude H, lightcurve information, information on the presence of ices and color indices.
Results of the modified Spearman rank correlation analysis of the Plutino sample.
All Figures
Fig. 1 Our Plutino sample (filled circles) compared to the known Plutino population (open circles, data from MPC, as of 18 Nov. 2011) in e–i phase space. The crosses with solid and dashed error bars refer to the mean and the standard deviation of our sample and the known Plutino population, respectively. This plot demonstrates that our sample is heterogenous and reflects the dynamical spread of the complete Plutino population. (The plot omits 3 known Plutinos, due to their high eccentricities, which are not part of our sample; however, these objects are included in the determination of the mean.) 

In the text 
Fig. 2 Illustration of the rescaled flux density uncertainties technique using synthetic data. a) Filled circles represent the flux densities at the different wavelengths with original flux density uncertainties. The dark line represents the best fit line (); the grey lines are fit lines of 100 out of 1000 synthetical bodies generated for the MC method. The bad quality of the fit is obvious, since the 70 and 100 μm datapoints are hardly fit by any of the single fit lines. Therefore, the reliability of the uncertainties of the output parameters is questionable. b) By rescaling the error bars by a factor of , the broadness of the distribution of single fit lines is increased, which covers now all datapoints. As a result, the output parameter fractional uncertainties are increased and therefore more realistic. 

In the text 
Fig. 3 Model fits of the Plutino sample using different modeling approaches and different data samples: the solid and dashed lines are floatingη fits based on combined Herschel and Spitzer and Herschelonly data, respectively; the dashdotted and dashdotdotted lines represent fixedη fits of combined Herschel and Spitzer and Herschelonly data samples, respectively; datapoints at 23.68 μm and 71.42 μm are Spitzer MIPS data (open circles), those at 70 μm, 100 μm and 160 μm are HerschelPACS data (filled circles). The model is fit to each flux measurement using the circumstances appropriate for the epoch of that observation. In order to simplify the figures, we have normalized the Spitzer flux densities to the epoch of the Herschel observations using the ratio of the bestfit model flux densities for the Herschel and Spitzer epoch, and then plotted the measured value times this ratio. Uncertainties are rescaled similarly, preserving the SNR for the Spitzer data. 

In the text 
Fig. 3 continued. 

In the text 
Fig. 3 continued. 

In the text 
Fig. 4 Diameter and albedo probability density functions and histograms. The probability density functions describe the probability to find an object of a specific diameter or albedo, when picking a random object. The functions were derived by summing up the probability density functions of each object, being represented by a nonsymmetric normal distribution taking into account the different upper and lower uncertainties of each result. This approach returns a more detailed picture than the histogram. Pluto was excluded from the sample. 

In the text 
Fig. 5 Cumulative Plutino size distributions N(≥ d) based on MPC H magnitudes and different assumptions: the solid and dashed line represent monochromatic albedo distributions using p_{V} = 0.08 and p_{V} = 0.11, respectively. The grey line shows the result of an averaged Monte Carlo simulations based on 100 runs and the measured albedo distribution with corresponding 1σ standard deviations. It is obvious that the monochromatic p_{V} = 0.08 size distribution agrees better with the result of the Monte Carlo simulation. Furthermore, we show two cumulative power laws N(≥d) ~ d^{−q} with q = 2 and q = 3 and arbitrary interceptions for comparison. A q = 2 cumulative power law better fits the measured size distribution in the size range 120–450 km, whereas the distribution of larger objects is better fit by a q = 3 cumulative power law. We can give an upper size limit for the kink at the smallsized end around 120 km. The upper abscissa gives the equivalent H magnitude based on a monochromatic albedo distribution of p_{V} = 0.08. 

In the text 
Fig. 6 Plots of some correlations. Upper left: we find no correlation between diameter and albedo; upper right: a clear anticorrelation between eccentricity and diameter is the result of a detection bias; lower left: we find smaller Plutinos to be redder than larger ones, which contradicts what was previously suggested by Peixinho et al. (2004) based on the intrinsical brightness H of a different Plutino sample; lower right: we find that the color of Plutinos is correlated to the heliocentric distance at the time of their discovery, which hints to a color bias. 

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.