In the Solar System, minor bodies and dust deliver various materials to planetary surfaces. Several exoplanetary systems are known to host inner and outer belts, analogues of the main asteroid belt and the Kuiper belt. We study the possibility that exominor bodies and exodust deliver volatiles and refractories to the exoplanets in the system HR8799 by performing N-body simulations. The model consists of the host star, four giant planets (HR8799 e, d, c, and b), 650000 test particles representing the inner belt, and 1450000 test particles representing the outer belt. Moreover we modelled dust populations that originate from both belts. Within a million years, the two belts evolve towards the expected dynamical structure (also derived in other works), where mean-motion resonances with the planets carve the analogues of Kirkwood gaps. We find that, after this point, the planets suffer impacts by objects from the inner and outer belt at rates that are essentially constant with time, while dust populations do not contribute significantly to the delivery process. We convert the impact rates to volatile and refractory delivery rates using our best estimates of the total mass contained in the belts and their volatile and refractory content. Over their lifetime, the four giant planets receive between $10^{-4}$ and $10^{-3}M_\bigoplus$ of material from both belts. The total amount of delivered volatiles and refractories, ${5\times10^{-3}\textrm{M}_\bigoplus}$, is small compared to the total mass of the planets, $11\times10^{3}\textrm{M}_\bigoplus$. However, if the planets were formed to be volatile-rich, their exogenous enrichment in refractory material may well be significant and observable, for example with JWST-MIRI. If terrestrial planets exist within the snow line of the system, volatile delivery would be an important astrobiological mechanism and may be observable as atmospheric trace gases.

Assuming our Solar System as typical, exomoons may outnumber exoplanets. If their habitability fraction is similar, they would thus constitute the largest portion of habitable real estate in the Universe. Icy moons in our Solar System, such as Europa and Enceladus, have already been shown to possess liquid water, a prerequisite for life on Earth. We intend to investigate under what circumstances small, icy moons may sustain subsurface oceans and thus be "subsurface habitable". We pay specific attention to tidal heating. We made use of a phenomenological approach to tidal heating. We computed the orbit averaged flux from both stellar and planetary (both thermal and reflected stellar) illumination. We then calculated subsurface temperatures depending on illumination and thermal conduction to the surface through the ice shell and an insulating layer of regolith. We adopted a conduction only model, ignoring volcanism and ice shell convection as an outlet for internal heat. In doing so, we determined at which depth, if any, ice melts and a subsurface ocean forms. We find an analytical expression between the moon's physical and orbital characteristics and the melting depth. Since this expression directly relates icy moon observables to the melting depth, it allows us to swiftly put an upper limit on the melting depth for any given moon. We reproduce the existence of Enceladus' subsurface ocean; we also find that the two largest moons of Uranus (Titania & Oberon) could well sustain them. Our model predicts that Rhea does not have liquid water. Habitable exomoon environments may be found across an exoplanetary system, largely irrespective of the distance to the host star. Small, icy subsurface habitable moons may exist anywhere beyond the snow line. This may, in future observations, expand the search area for extraterrestrial habitable environments beyond the circumstellar habitable zone.

Binospec is a high throughput, 370 to 1000 nm, imaging spectrograph that addresses two adjacent 8' by 15' fields of view. Binospec was commissioned in late 2017 at the f/5 focus of the 6.5m MMT and is now available to all MMT observers. Aperture masks cut from stainless steel with a laser cutter are used to define the entrance apertures that range from 15' long slits to hundreds of 2" slitlets. System throughputs, including the MMT's mirrors and the f/5 wide-field corrector peak at ~30%. Three reflection gratings, duplicated for the two beams, provide resolutions ($\lambda$/$\Delta \lambda$) between 1300 and $>$5000 with a 1" wide slit. Two through-the-mask guiders are used for target acquisition, mask alignment, guiding, and precision offsets. A full-time Shack-Hartmann wave front sensor allows continuous adjustment of primary mirror support forces, telescope collimation and focus. Active flexure control maintains spectrograph alignment and focus under varying gravity and thermal conditions.

A group of trans-Neptunian objects (TNO) are dynamically related to the dwarf planet 136108 Haumea. Ten of them show strong indications of water ice on their surfaces, are assumed to have resulted from a collision, and are accepted as the only known TNO collisional family. Nineteen other dynamically similar objects lack water ice absorptions and are hypothesized to be dynamical interlopers. We have made observations to determine sizes and geometric albedos of six of the accepted Haumea family members and one dynamical interloper. Ten other dynamical interlopers have been measured by previous works. We compare the individual and statistical properties of the family members and interlopers, examining the size and albedo distributions of both groups. We also examine implications for the total mass of the family and their ejection velocities. We use far-infrared space-based telescopes to observe the target TNOs near their thermal peak and combine these data with optical magnitudes to derive sizes and albedos using radiometric techniques. We determine the power-law slope of ejection velocity as a function of effective diameter. The detected Haumea family members have a diversity of geometric albedos $\sim$ 0.3-0.8, which are higher than geometric albedos of dynamically similar objects without water ice. The median geometric albedo for accepted family members is $p_V=0.48_{-0.18}^{+0.28}$, compared to 0.08$_{-0.05}^{+0.07}$ for the dynamical interlopers. In the size range $D=175-300$ km, the slope of the cumulative size distribution is $q$=3.2$_{-0.4}^{+0.7}$ for accepted family members, steeper than the $q$=2.0$\pm$0.6 slope for the dynamical interlopers with D$< $500 km. The total mass of Haumea's moons and family members is 2.4% of Haumea's mass. The ejection velocities required to emplace them on their current orbits show a dependence on diameter, with a power-law slope of 0.21-0.50.

1I/`Oumuamua is the first confirmed interstellar body in our Solar System. Here we report on observations of `Oumuamua made with the Spitzer Space Telescope on 2017 November 21--22 (UT). We integrated for 30.2~hours at 4.5 micron (IRAC channel 2). We did not detect the object and place an upper limit on the flux of 0.3 uJy (3sigma). This implies an effective spherical diameter less than [98, 140, 440] meters and albedo greater than [0.2, 0.1, 0.01] under the assumption of low, middle, or high thermal beaming parameter eta, respectively. With an aspect ratio for `Oumuamua of 6:1, these results correspond to dimensions of [240:40, 341:57, 1080:180] meters, respectively. We place upper limits on the amount of dust, CO, and CO2 coming from this object that are lower than previous results; we are unable to constrain the production of other gas species. Both our size and outgassing limits are important because `Oumuamua's trajectory shows non-gravitational accelerations that are sensitive to size and mass and presumably caused by gas emission. We suggest that `Oumuamua may have experienced low-level post-perihelion volatile emission that produced a fresh, bright, icy mantle. This model is consistent with the expected eta value and implied high albedo value for this solution, but, given our strict limits on CO and CO2, requires another gas species --- probably H2O --- to explain the observed non-gravitational acceleration. Our results extend the mystery of `Oumuamua's origin and evolution.

We present lightcurves and derive periods and amplitudes for a subset of 38 near earth objects (NEOs) observed at 4.5 microns with the IRAC camera on the the Spitzer Space Telescope, many of them having no previously reported rotation periods. This subset was chosen from about 1800 IRAC NEO observations as having obvious periodicity and significant amplitude. For objects where the period observed did not sample the full rotational period, we derived lower limits to these parameters based on sinusoidal fits. Lightcurve durations ranged from 42 to 544 minutes, with derived periods from 16 to 400 minutes. We discuss the effects of lightcurve variations on the thermal modeling used to derive diameters and albedos from Spitzer photometry. We find that both diameters and albedos derived from the lightcurve maxima and minima agree with our previously published results, even for extreme objects, showing the conservative nature of the thermal model uncertainties. We also evaluate the NEO rotation rates, sizes, and their cohesive strengths.

Given rapid photodissociation and photodegradation, the recently discovered organics in the Martian subsurface and atmosphere were probably delivered in geologically recent times. Possible parent bodies are C-type asteroids, comets, and interplanetary dust particles (IDPs). The dust infall rate was estimated, using different methods, to be between $0.71$ and $2.96 \times 10^6$ kg/yr (Nesvorny et al., 2011, Borin et al., 2017, Crismani et al., 2017); assuming a carbon content of 10% (Flynn, 1996), this implies an IDP carbon flux of $0.07 - 0.3 \times 10^6$ kg/yr. We calculate for the first time the carbon flux from impacts of asteroids and comets. To this end, we perform dynamical simulations of impact rates on Mars. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We separately add comets and asteroids to the simulations as massless test particles, based on their current orbital elements, yielding Mars impact rates of $4.34\times10^{-3}$ comets/Myr and 3.3 asteroids/Myr. We estimate the global carbon flux on Mars from cometary impacts to be $\sim 0.013 \times 10^{6}$~kg/yr within an order of magnitude, while asteroids deliver $\sim 0.05 \times 10^6$~kg/yr. These values correspond to $\sim 4-19 \%$ and $\sim 17-71 \%$, respectively, of the IDP-borne carbon flux estimated by Nesvorny et al. 2011, Borin et al. 2017 and Crismani et al. 2017. Unlike the spatially homogeneous IDP infall, impact ejecta are distributed locally, concentrated around the impact site. We find organics from asteroids and comets to dominate over IDP-borne organics at distances up to 150~km from the crater center. Our results may be important for the interpretation of in situ detections of organics on Mars.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.

{Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness.

The aim of the Asteroid Science Intersections with In-Space Mine Engineering (ASIME) 2016 conference on September 21-22, 2016 in Luxembourg City was to provide an environment for the detailed discussion of the specific properties of asteroids, with the engineering needs of space missions that utilize asteroids. The ASIME 2016 Conference produced a layered record of discussions from the asteroid scientists and the asteroid miners to understand each other's key concerns and to address key scientific questions from the asteroid mining companies: Planetary Resources, Deep Space Industries and TransAstra. These Questions were the focus of the two day conference, were addressed by scientists inside and outside of the ASIME Conference and are the focus of this White Paper. The Questions from the asteroid mining companies have been sorted into the three asteroid science themes: 1) survey, 2) surface and 3) subsurface and 4) Other. The answers to those Questions have been provided by the scientists with their conference presentations or edited directly into an early open-access collaborative Google document (August 2016-October 2016), or inserted by A. Graps using additional reference materials. During the ASIME 2016 last two-hours, the scientists turned the Questions from the Asteroid Miners around by presenting their own key concerns: Questions from the Asteroid Scientists . These answers in this White Paper will point to the Science Knowledge Gaps (SKGs) for advancing the asteroid in-space resource utilisation domain.

We compute the entanglement and R\'enyi entropy growth after a global quench in various dimensions in free scalar field theory. We study two types of quenches: a boundary state quench and a global mass quench. Both of these quenches are investigated for a strip geometry in 1, 2, and 3 spatial dimensions, and for a spherical geometry in 2 and 3 spatial dimensions. We compare the numerical results for massless free scalars in these geometries with the predictions of the analytical quasiparticle model based on EPR pairs, and find excellent agreement in the limit of large region sizes. At subleading order in the region size, we observe an anomalous logarithmic growth of entanglement coming from the zero mode of the scalar.

The JAXA Hayabusa-2 mission was approved in 2010 and launched on December 3, 2014. The spacecraft will arrive at the near-Earth asteroid 162173 Ryugu in 2018 where it will perform a survey, land and obtain surface material, then depart in Dec 2019 and return to Earth in Dec 2020. We observed Ryugu with the Herschel Space Observatory in Apr 2012 at far-IR thermal wavelengths, supported by several ground-based observations to obtain optical lightcurves. We reanalysed previously published Subaru-COMICS and AKARI-IRC observations and merged them with a Spitzer-IRS data set. In addition, we used a large set of Spitzer-IRAC observations obtained in the period Jan to May, 2013. The data set includes two complete rotational lightcurves and a series of ten "point-and-shoot" observations. The almost spherical shape of the target together with the insufficient lightcurve quality forced us to combine radiometric and lightcurve inversion techniques in different ways to find the object's key physical and thermal parameters. We find that the solution which best matches our data sets leads to this C class asteroid having a retrograde rotation with a spin-axis orientation of (lambda = 310-340 deg; beta = -40+/-15 deg) in ecliptic coordinates, an effective diameter (of an equal-volume sphere) of 850 to 880 m, a geometric albedo of 0.044 to 0.050 and a thermal inertia in the range 150 to 300 Jm-2s-0.5K-1. Based on estimated thermal conductivities of the top-layer surface in the range 0.1 to 0.6 WK-1m-1, we calculated that the grain sizes are approximately equal to between 1 and 10 mm. The finely constrained values for this asteroid serve as a `design reference model', which is currently used for various planning, operational and modelling purposes by the Hayabusa2 team.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Near Earth Objects (NEOs) are small Solar System bodies whose orbits bring them close to the Earth's orbit. We are carrying out a Warm Spitzer Cycle 11 Exploration Science program entitled NEOSurvey --- a fast and efficient flux-limited survey of 597 known NEOs in which we derive diameter and albedo for each target. The vast majority of our targets are too faint to be observed by NEOWISE, though a small sample has been or will be observed by both observatories, which allows for a cross-check of our mutual results. Our primary goal is to create a large and uniform catalog of NEO properties. We present here the first results from this new program: fluxes and derived diameters and albedos for 80 NEOs, together with a description of the overall program and approach, including several updates to our thermal model. The largest source of error in our diameter and albedo solutions, which derive from our single band thermal emission measurements, is uncertainty in eta, the beaming parameter used in our thermal modeling; for albedos, improvements in Solar System absolute magnitudes would also help significantly. All data and derived diameters and albedos from this entire program are being posted on a publicly accessible webpage at nearearthobjects.nau.edu .

Photometric lightcurves were also obtained for 14 of them during the Spitzer observations to provide the context of the observations and reliable estimates of their absolute magnitudes. The extracted mid-IR spectra were analyzed using a modified standard thermal model (STM) and a thermophysical model (TPM) that takes into account the shape and geometry of the large primary at the time of the Spitzer observation. We derived a reliable estimate of the size, albedo, and beaming factor for each of these asteroids, representing three main taxonomic groups: C, S, and X. For large (volume-equivalent system diameter Deq $\lt$ 130 km) binary asteroids, the TPM analysis indicates a low thermal inertia ($\Gamma$ < $\sim$100 J s-1/2K-1m-2) and their emissivity spectra display strong mineral features, implying that they are covered with a thick layer of thermally insulating regolith. The smaller (surface-equivalent system diameter Deff $\lt$17 km) asteroids also show some emission lines of minerals, but they are significantly weaker, consistent with regoliths with coarser grains, than those of the large binary asteroids. The average bulk densities of these multiple asteroids vary from 0.7-1.7 g/cm3 (P-, C- type) to $\sim$2 g/cm3 (S-type). The highest density is estimated for the M-type (22) Kalliope (3.2 $\pm$ 0.9 g/cm3). The spectral energy distributions (SED) and emissivity spectra, made available as a supplement document, could help to constrain the surface compositions of these asteroids.

Studies of the physical properties of Trans-Neptunian Objects (TNOs) are a powerful probe into the processes of planetesimal formation and solar system evolution. JWST will provide unique new capabilities for such studies. Here we outline where the capabilities of JWST open new avenues of investigation, potential valuable observations and surveys, and conclude with a discussion of community actions that may serve to enhance the eventual science return of JWSTs TNO observations.

The James Webb Space Telescope (JWST) has the potential to enhance our understanding of near-Earth objects (NEOs). We present results of investigations into the observability of NEOs given the nominal observing requirements of JWST on elongation (85-135 degrees) and non-sidereal rates ($<$30mas/s). We find that approximately 75% of NEOs can be observed in a given year. However, observers will need to wait for appropriate observing windows. We find that JWST can easily execute photometric observations of meter-sized NEOs which will enhance our understanding of the small NEO population.

The field of asteroid thermophysical modeling has experienced an extraordinary growth in the last ten years, as new thermal infrared data became available for hundreds of thousands of asteroids. The infrared emission of asteroids depends on the body's size, shape, albedo, thermal inertia, roughness and rotational properties. These parameters can therefore be derived by thermophysical modeling of infrared data. Thermophysical modeling led to asteroid size estimates that were confirmed at the few-percent level by later spacecraft visits. We discuss how instrumentation advances now allow mid-infrared interferometric observations as well as high-accuracy spectro-photometry, posing their own set of thermal-modeling challenges.We present major breakthroughs achieved in studies of the thermal inertia, a sensitive indicator for the nature of asteroids soils, allowing us, for instance, to determine the grain size of asteroidal regoliths. Thermal inertia also governs non-gravitational effects on asteroid orbits, requiring thermophysical modeling for precise asteroid dynamical studies. The radiative heating of asteroids, meteoroids, and comets from the Sun also governs the thermal stress in surface material; only recently has it been recognized as a significant weathering process. Asteroid space missions with thermal infrared instruments are currently undergoing study at all major space agencies. This will require a high level of sophistication of thermophysical models in order to analyze high-quality spacecraft data.

We perform a search for dormant comets, asteroidal objects of cometary origin, in the near-Earth asteroid (NEA) population based on dynamical and physical considerations. Our study is based on albedos derived within the ExploreNEOs program and is extended by adding data from NEOWISE and the Akari asteroid catalog. We use a statistical approach to identify asteroids on orbits that resemble those of short-period near-Earth comets using the Tisserand parameter with respect to Jupiter, the aphelion distance, and the minimum orbital intersection distance with respect to Jupiter. From the sample of NEAs on comet-like orbits, we select those with a geometric albedo $p_V \leq 0.064$ as dormant comet candidates, and find that only $\sim$50% of NEAs on comet-like orbits also have comet-like albedos. We identify a total of 23 NEAs from our sample that are likely to be dormant short-period near-Earth comets and, based on a de-biasing procedure applied to the cryogenic NEOWISE survey, estimate both magnitude-limited and size-limited fractions of the NEA population that are dormant short-period comets. We find that 0.3-3.3% of the NEA population with $H \leq 21$, and $9^{+2}_{-5}$% of the population with diameters $d \geq 1$ km, are dormant short-period near-Earth comets.

We report on observations of near-Earth asteroid 2011 MD with the Spitzer Space Telescope. We have spent 19.9 h of observing time with channel 2 (4.5 {\mu}m) of the Infrared Array Camera and detected the target within the 2{\sigma} positional uncertainty ellipse. Using an asteroid thermophysical model and a model of nongravitational forces acting upon the object we constrain the physical properties of 2011 MD, based on the measured flux density and available astrometry data. We estimate 2011 MD to be 6 (+4/-2) m in diameter with a geometric albedo of 0.3 (+0.4/-0.2) (uncertainties are 1{\sigma}). We find the asteroid's most probable bulk density to be 1.1 (+0.7/-0.5) g cm^{-3}, which implies a total mass of (50-350) t and a macroporosity of >=65%, assuming a material bulk density typical of non-primitive meteorite materials. A high degree of macroporosity suggests 2011 MD to be a rubble-pile asteroid, the rotation of which is more likely to be retrograde than prograde.

We report on Spitzer Space Telescope IRAC observations of near-Earth object (NEO) 2009 BD that were carried out in support of the NASA Asteroid Robotic Retrieval Mission (ARRM) concept. We did not detect 2009 BD in 25 hrs of integration at 4.5 micron. Based on an upper-limit flux density determination from our data, we present a probabilistic derivation of the physical properties of this object. The analysis is based on the combination of a thermophysical model with an orbital model accounting for the non-gravitational forces acting upon the body. We find two physically possible solutions. The first solution shows 2009 BD as a 2.9+/-0.3 m diameter rocky body (rho = 2.9+/-0.5 g cm-3) with an extremely high albedo of 0.85(+0.20/-0.10) that is covered with regolith-like material, causing it to exhibit a low thermal inertia (Gamma = 30(+20/-10) SI units). The second solution suggests 2009 BD to be a 4+/-1 m diameter asteroid with pV = 0.45(+0.35/-0.15) that consists of a collection of individual bare rock slabs (Gamma = 2000+/-1000 SI units, rho = 1.7(+0.7/-0.4) g cm-3). We are unable to rule out either solution based on physical reasoning. 2009 BD is the smallest asteroid for which physical properties have been constrained, in this case using an indirect method and based on a detection limit, providing unique information on the physical properties of objects in the size range smaller than 10 m.

The classical Kuiper belt contains objects both from a low-inclination, presumably primordial, distribution and from a high-inclination dynamically excited population. Based on a sample of classical TNOs with observations at thermal wavelengths we determine radiometric sizes, geometric albedos and thermal beaming factors as well as study sample properties of dynamically hot and cold classicals. Observations near the thermal peak of TNOs using infra-red space telescopes are combined with optical magnitudes using the radiometric technique with near-Earth asteroid thermal model (NEATM). We have determined three-band flux densities from Herschel/PACS observations at 70.0, 100.0 and 160.0 $\mu$m and Spitzer/MIPS at 23.68 and 71.42 $\mu$m when available. We have analysed 18 classical TNOs with previously unpublished data and re-analysed previously published targets with updated data reduction to determine their sizes and geometric albedos as well as beaming factors when data quality allows. We have combined these samples with classical TNOs with radiometric results in the literature for the analysis of sample properties of a total of 44 objects. We find a median geometric albedo for cold classical TNOs of 0.14 and for dynamically hot classical TNOs, excluding the Haumea family and dwarf planets, 0.085. We have determined the bulk densities of Borasisi-Pabu (2.1 g/cm^3), Varda-Ilmare (1.25 g/cm^3) and 2001 QC298 (1.14 g/cm^3) as well as updated previous density estimates of four targets. We have determined the slope parameter of the debiased cumulative size distribution of dynamically hot classical TNOs as q=2.3 +- 0.1 in the diameter range 100<D<500 km. For dynamically cold classical TNOs we determine q=5.1 +- 1.1 in the diameter range 160<D<280 km as the cold classical TNOs have a smaller maximum size.

Centaurs are the transitional population between trans-Neptunian objects (TNOs) and Jupiter-family comets. For this reason it is possible to access the smaller ones, which is more difficult to do with the TNO population. The goal of this work is to characterize a set of 16 Centaurs in terms of their size, albedo, and thermal properties. We study the correlations, for a more extended sample obtained from the literature, of diameter, albedo, orbital parameters, and spectral slopes. We performed three-band photometric observations using Herschel-PACS and used a consistent method for the data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70, 100, and 160 $\mu$m. Additionally, we used Spitzer-MIPS flux densities at 24 and 70 $\mu$m when available. We also included in our Centaur sample scattered disk objects (SDOs), a dynamical family of TNOS, using results previously published by our team, and some Centaurs observed only with the Spitzer/MIPS instrument. We have determined new radiometric sizes and albedos of 16 Centaurs. The first conclusion is that the albedos of Centaur objects are not correlated with their orbital parameters. Similarly, there is no correlation between diameter and orbital parameters. Most of the objects in our sample are dark (pv < 7%) and most of them are small (D < 120km). However, there is no correlation between albedo and diameter, in particular for the group of the small objects as albedo values are homogeneously distributed between 4 to 16%. The correlation with the color of the objects showed that red objects are all small (mean diameter 65 km), while the gray ones span a wide range of sizes (mean diameter 120 km). Moreover, the gray objects tend to be darker, with a mean albedo of 5.6%, compared with a mean of 8.5% (ranging from 5 to 15%) for the red objects.

The near-Earth object (NEO) population, which mainly consists of fragments from collisions between asteroids in the main asteroid belt, is thought to include contributions from short-period comets as well. One of the most promising NEO candidates for a cometary origin is near-Earth asteroid (3552) Don Quixote, which has never been reported to show activity. Here we present the discovery of cometary activity in Don Quixote based on thermal-infrared observations made with the Spitzer Space Telescope in its 3.6 and 4.5 {\mu}m bands. Our observations clearly show the presence of a coma and a tail in the 4.5 {\mu}m but not in the 3.6 {\mu}m band, which is consistent with molecular band emission from CO2. Thermal modeling of the combined photometric data on Don Quixote reveals a diameter of 18.4 (-0.4/+0.3) km and an albedo of 0.03 (-0.01/+0.02), which confirms Don Quixote to be the third-largest known NEO. We derive an upper limit on the dust production rate of 1.9 kg s^-1 and derive a CO2 gas production rate of (1.1+-0.1)10^26 molecules s^-1. Spitzer IRS spectroscopic observations indicate the presence of fine-grained silicates, perhaps pyroxene rich, on the surface of Don Quixote. Our discovery suggests that CO2 can be present in near-Earth space over a long time. The presence of CO2 might also explain that Don Quixote's cometary nature remained hidden for nearly three decades.

Near-infrared spectroscopy of Near-Earth Objects (NEOs) connects diagnostic spectral features to specific surface mineralogies. The combination of spectroscopy with albedos and diameters derived from thermal infrared observations can increase the scientific return beyond that of the individual datasets. To that end, we have completed a spectroscopic observing campaign to complement the ExploreNEOs Warm Spitzer program that obtained albedos and diameters of nearly 600 NEOs (Trilling et al. 2010). Here we present the results of observations using the low-resolution prism mode (~0.7-2.5 microns) of the SpeX instrument on the NASA Infrared Telescope Facility (IRTF). We also include near-infrared observations of ExploreNEOs targets from the MIT-UH-IRTF Joint Campaign for Spectral Reconnaissance. Our dataset includes near-infrared spectra of 187 ExploreNEOs targets (125 observations of 92 objects from our survey and 213 observations of 154 objects from the MIT survey). We identify a taxonomic class for each spectrum and use band parameter analysis to investigate the mineralogies for the S-, Q-, and V-complex objects. Our analysis suggests that for spectra that contain near-infrared data but lack the visible wavelength region, the Bus-DeMeo system misidentifies some S-types as Q-types. We find no correlation between spectral band parameters and ExploreNEOs albedos and diameters. We find slightly negative Band Area Ratio (BAR) correlations with phase angle for Eros and Ivar, but a positive BAR correlation with phase angle for Ganymed. We find evidence for spectral phase reddening for Eros, Ganymed, and Ivar. We identify the likely ordinary chondrite type analog for a subset of our sample. Our resulting proportions of H, L, and LL ordinary chondrites differ from those calculated for meteorite falls and in previous studies of ordinary chondrite-like NEOs.

We observed nine primary transits of the super-Earth exoplanet GJ 1214b in several optical photometric bands from March to August 2012, with the goal of constraining the short-wavelength slope of the spectrum of GJ 1214b. Our observations were conducted on the Kuiper 1.55 m telescope in Arizona and the STELLA-I robotic 1.2 m telescope in Tenerife, Spain. From the derived light curves we extracted transit depths in R (0.65 {\mu}m), V (0.55 {\mu}m), and g' (0.475 {\mu}m) bands. Most previous observations of this exoplanet suggest a flat spectrum varying little with wavelength from the near-infrared to the optical, corresponding to a low-scale-height, high-molecular-weight atmosphere. However, a handful of observations around Ks band (~2.15 {\mu}m) and g-band (~0.46 {\mu}m) are inconsistent with this scenario and suggest a variation on a hydrogen- or water-dominated atmosphere that also contains a haze layer of small particles. In particular, the g-band observations of de Mooij et al. (2012), consistent with Rayleigh scattering, limit the potential atmosphere compositions of GJ 1214b due to the increasing slope at optical wavelengths (Howe & Burrows 2012). We find that our results overlap within errors the short-wavelength observations of de Mooij et al. (2012), but are also consistent with a spectral slope of zero in GJ 1214b in the optical wavelength region. Our observations thus allow for a larger suite of possible atmosphere compositions, including those with a high-molecular-weight and/or hazes.

The subject of this work is the physical characterization of asteroids, focusing on the thermal inertia of near-Earth asteroids (NEAs). Thermal inertia governs the Yarkovsky effect, a non-gravitational force which significantly alters the orbits of asteroids up to \sim 20 km in diameter. Yet, very little has previously been known about the thermal inertia of small asteroids including NEAs. Observational and theoretical work is reported. The thermal emission of asteroids has been observed in the mid-infrared (5-35 {\mu}m) wavelength range using the Spitzer Space Telescope and the 3.0m IRTF. A detailed thermophysical model (TPM) has been developed and extensively tested; this is the first detailed TPM shown to be applicable to NEA data. Our main result is the determination of the thermal inertia of 5 NEAs, increasing the total number of NEAs with measured thermal inertia to 6. For two of our targets, previously available estimates are refined. Our results allow the first determination of the typical thermal inertia of NEAs, which is around 300 J s^{-1/2} K^{-1} m^{-2}, larger than the typical thermal inertia of large main-belt asteroids (MBAs) by more than an order of magnitude. In particular, thermal inertia appears to increase with decreasing asteroid diameter. Our results have been used by colleagues to estimate the size dependence of the Yarkovsky effect, thus explaining the apparent difference in the size-frequency distribution of NEAs and similarly sized MBAs. Thermal inertia is a very sensitive indicator for the presence or absence of particulate material on the surface. Our results indicate that even sub-km asteroids are covered with coarse regolith.

Trans-Neptunian objects (TNO) represent the leftovers of the formation of the Solar System. Their physical properties provide constraints to the models of formation and evolution of the various dynamical classes of objects in the outer Solar System. Based on a sample of 19 classical TNOs we determine radiometric sizes, geometric albedos and beaming parameters. Our sample is composed of both dynamically hot and cold classicals. We study the correlations of diameter and albedo of these two subsamples with each other and with orbital parameters, spectral slopes and colors. We have done three-band photometric observations with Herschel/PACS and we use a consistent method for data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70.0, 100.0 and 160.0 \mu m. Additionally, we use Spitzer/MIPS flux densities at 23.68 and 71.42 \mu m when available, and we present new Spitzer flux densities of eight targets. We derive diameters and albedos with the near-Earth asteroid thermal model (NEATM). As auxiliary data we use reexamined absolute visual magnitudes from the literature and data bases, part of which have been obtained by ground based programs in support of our Herschel key program. We have determined for the first time radiometric sizes and albedos of eight classical TNOs, and refined previous size and albedo estimates or limits of 11 other classicals. The new size estimates of 2002 MS4 and 120347 Salacia indicate that they are among the 10 largest TNOs known. Our new results confirm the recent findings that there are very diverse albedos among the classical TNOs and that cold classicals possess a high average albedo (0.17 +/- 0.04). Diameters of classical TNOs strongly correlate with orbital inclination in our sample. We also determine the bulk densities of six binary TNOs.

The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interesting spacecraft target: a binary asteroid with a low-DeltaV heliocentric orbit. The orbit of its satellite has provided valuable information about its mass density while its albedo and colors suggest it is primitive or part of the C-complex taxonomic grouping. We extend the physical characterization of this object with new observations of its emission at mid-Infrared (IR) wavelengths and with near-IR reflection spectroscopy. We derive an area-equivalent system diameter of 1.90 \pm 0.28 km (corresponding to approximate component diameters of 1.83 km and 0.51 km, respectively) and a geometric albedo of 0.039 \pm 0.012. 1996 FG3 was previously classified as a C-type asteroid, though the combined 0.4--2.5 micron spectrum with thermal correction indicates classification as B-type; both are consistent with the low measured albedo. Dynamical studies show that 1996 FG3 has most probably originated in the inner main asteroid belt. Recent work has suggested the inner Main Belt (142) Polana family as the possible origin of another low-DeltaV B-type NEA, (101955) 1999 RQ36. A similar origin for 1996 FG3 would require delivery by the overlapping Jupiter 7:2 and Mars 5:9 mean motion resonances rather than the nu-6 resonance, and we find this to be a low probability, but possible, origin.

We present Herschel PACS 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. Flux densities are derived from PACS mini scan-maps 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. 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 \pm 0.03, which is comparable to the mean albedo of Centaurs, Jupiter Family comets and other Trans-Neptunian Objects. 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 a density of 0.64 (+0.15/-0.11) 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.

Physical characterization of Trans-Neptunian objects, a primitive population of the outer solar system, may provide constraints on their formation and evolution. 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. Thermal flux measurements obtained with the Herschel-PACS instrument at 70, 100 and 160 \mu m, and whenever applicable, with Spitzer-MIPS at 24 and 70 \mu 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. 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.

Space missions to NEOs are being planned at all major space agencies, and recently a manned mission to an NEO was announced as a NASA goal. Efforts to find and select suitable targets (plus backup targets) are severely hampered by our lack of knowledge of the physical properties of dynamically favorable NEOs. In particular, current mission scenarios tend to favor primitive low-albedo objects. For the vast majority of NEOs the albedo is unknown. Here we report new constraints on the size and albedo of 65 NEOs with rendezvous deltaV < 7 km/s. Our results are based on thermal-IR flux data obtained in the framework of our ongoing (2009--2011) ExploreNEOs survey (Trilling et al. 2010) using NASA's "Warm Spitzer" space telescope. As of 2010 July 14, we have results for 293 objects in hand (including the 65 low-deltaV NEOs presented here); before the end of 2011 we expect to have measured the size and albedo of ~ 700 NEOs (including probably ~ 160 low-deltaV NEOs). While there are reasons to believe that primitive volatile-rich materials are universally low in albedo, the converse need not be true: the orbital evolution of some dark objects likely has caused them to lose their volatiles by coming too close to the Sun. For all our targets, we give the closest perihelion distance they are likely to have reached (using orbital integrations from Marchi et al. 2009) and corresponding upper limits on the past surface temperature. Low-deltaV objects for which both albedo and thermal history may suggest a primitive composition include (162998) 2001 SK162, (68372) 2001 PM9, and (100085) 1992 UY4.

We have begun the ExploreNEOs project in which we observe some 700 Near Earth Objects (NEOs) at 3.6 and 4.5 microns with the Spitzer Space Telescope in its Warm Spitzer mode. From these measurements and catalog optical photometry we derive albedos and diameters of the observed targets. The overall goal of our ExploreNEOs program is to study the history of near-Earth space by deriving the physical properties of a large number of NEOs. In this paper we describe both the scientific and technical construction of our ExploreNEOs program. We present our observational, photometric, and thermal modeling techniques. We present results from the first 101 targets observed in this program. We find that the distribution of albedos in this first sample is quite broad, probably indicating a wide range of compositions within the NEO population. Many objects smaller than one kilometer have high albedos (>0.35), but few objects larger than one kilometer have high albedos. This result is consistent with the idea that these larger objects are collisionally older, and therefore possess surfaces that are more space weathered and therefore darker, or are not subject to other surface rejuvenating events as frequently as smaller NEOs.

Thermal emission from Kuiper Belt object (136108) Haumea was measured with Herschel-PACS at 100 and 160 micrometers for almost a full rotation period. Observations clearly indicate a 100-micrometer thermal lightcurve with an amplitude of a factor of ~ 2, which is positively correlated with the optical lightcurve. This confirms that both are primarily due to shape effects. A 160-micrometer lightcurve is marginally detected. Radiometric fits of the mean Herschel- and Spitzer- fluxes indicate an equivalent diameter D ~ 1300 km and a geometric albedo p_v ~ 0.70-0.75. These values agree with inferences from the optical lightcurve, supporting the hydrostatic equilibrium hypothesis. The large amplitude of the 100-micrometer lightcurve suggests that the object has a high projected a/b axis ratio (~ 1.3) and a low thermal inertia as well as possible variable infrared beaming. This may point to fine regolith on the surface, with a lunar-type photometric behavior. The quality of the thermal data is not sufficient to clearly detect the effects of a surface dark spot.

The goal of the Herschel Open Time Key programme "TNOs are Cool!" is to derive the physical and thermal properties for a large sample of Centaurs and trans-Neptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. We present results for seven targets either observed in PACS point-source, or in mini scan-map mode. Spitzer-MIPS observations were included for three objects. The sizes of these targets range from 100 km to almost 1000 km, five have low geometric albedos below 10%, (145480) 2005 TB190 has a higher albedo above 15%. Classical thermal models driven by an intermediate beaming factor of $\eta$=1.2 or $\eta$-values adjusted to the observed colour temperature fit the multi-band observations well in most cases. More sophisticated thermophysical models give very similar diameter and albedo values for thermal inertias in the range 0-25 Jm-2s-0.5K-1, consistent with very low heat conductivities at temperatures far away from the Sun. The early experience with observing and model strategies will allow us to derive physical and thermal properties for our complete Herschel TNO sample of 140 targets as a benchmark for understanding the solar system debris disk, and extra-solar ones as well.

We report on observations of the dust trail of comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2m telescope at 4.7 AU before aphelion, and at 24 micron with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter beta<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of -4.1, and emission speeds for a beta=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>10 micron) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of 60-600 micron.

Dust is an important constituent in cometary comae; its analysis is one of the major objectives of ESA's Rosetta mission to comet 67P/Churyumov-Gerasimenko (C-G). Several instruments aboard Rosetta are dedicated to studying various aspects of dust in the cometary coma, all of which require a certain level of exposure to dust to achieve their goals. At the same time, impacts of dust particles can constitute a hazard to the spacecraft. To conciliate the demands of dust collection instruments and spacecraft safety, it is desirable to assess the dust environment in the coma even before the arrival of Rosetta. We describe the present status of modelling the dust coma of 67P/C-G and predict the speed and flux of dust in the coma, the dust fluence on a spacecraft along sample trajectories, and the radiation environment in the coma. The model will need to be refined when more details of the coma are revealed by observations. An overview of astronomical observations of 67P/C-G is given and model parameters are derived from these data where possible. For quantities not yet measured for 67P/C-G, we use values obtained for other comets. One of the most important and most controversial parameters is the dust mass distribution. We summarise the mass distribution functions derived from the in-situ measurements at comet 1P/Halley in 1986. For 67P/C-G, constraining the mass distribution is currently only possible by the analysis of astronomical images. We find that the results from such analyses are at present rather heterogeneous, and we identify a need to find a model that is reconcilable with all available observations.

We present mid-infrared (8-33 micron) observations of the binary L5-Trojan system (617) Patroclus-Menoetius before, during, and after two shadowing events, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope.F or the first time, we effectively observe changes in asteroid surface temperature in real time, allowing the thermal inertia to be determined very directly. A new detailed binary thermophysical model is presented which accounts for the system's known mutual orbit, arbitrary component shapes, and thermal conduction in the presence of eclipses. We obtain two local thermal-inertia values, representative of the respective shadowed areas: 21+/14 MKS and 6.4+/-1.6 MKS. The average thermal inertia is estimated to be 20+/-15 MKS, potentially with significant surface heterogeneity. This first thermal-inertia measurement for a Trojan asteroid indicates a surface covered in fine regolith. The diameters of Patroclus and Menoetius are 106 +/- 11 and 98+/-10 km, respectively, in agreement with previous findings. Taken together with the system's known total mass, this implies a bulk mass density of 1.08 +/-0.33 g/cm3, significantly below the mass density of L4-Trojan asteroid (624) Hektor and suggesting a bulk composition dominated by water ice.

The Method of Light Curve Simulations is a tool that has been applied to X-ray monitoring observations of Active Galactic Nuclei (AGN) for the characterization of the Power Density Spectrum (PDS) of temporal variability and measurement of associated break frequencies (which appear to be an important diagnostic for the mass of the black hole in these systems as well as their accretion state). It relies on a model for the PDS that is fit to the observed data. The determination of confidence regions on the fitted model parameters is of particular importance, and we show how the Neyman construction based on distributions of estimates may be implemented in the context of light curve simulations. We believe that this procedure offers advantages over the method used in earlier reports on PDS model fits, not least with respect to the correspondence between the size of the confidence region and the precision with which the data constrain the values of the model parameters. We plan to apply the new procedure to existing RXTE and XMM observations of Seyfert I galaxies as well as RXTE observations of the Seyfert II galaxy NGC 4945.

We re-analyze the Chandra ACIS spectrum of the kpc-scale jet in PKS 0637-752 to investigate the possible low energy cutoff in the relativistic electron spectrum producing the non-thermal radiation in the scenario of inverse Compton emission off the cosmic microwave background. This was among the first objects targeted by the Chandra Observatory and gives a unique opportunity to study the low energy X-ray emission free of detector contamination. As previously reported in the literature, the spectrum can be fit by a power law, with the slope predicted by the radio spectrum, modified by low energy absorption through the Galaxy as determined from the spectrum of the quasar core and by HI 21 cm observations. We obtain a marginally better fit with an model of inverse Compton emission produced by an electron population that exhibits a cutoff at (gamma_min delta_10) between about 50 and 80 (assuming Gamma = delta). This range for gamma_min is higher than has previously been assumed in broad-band spectral fits to the jet emission. The observed optical flux can be used to place a lower limit on gamma_min; the constraint is not very strong, but does suggest that gamma_min must be higher than 1 to avoid overproducing the optical emission. We investigate the effect of uncertainties in the column density for galactic absorption as well as the calibration of Chandra for these early observations. Finally, we discuss the implication of these limits on the jet luminosity in this source.

Near Earth Objects (NEOs) are fragments of remnant primitive bodies that date from the era of Solar System formation. At present, the physical properties and origins of NEOs are poorly understood. We have measured thermal emission from three NEOs -- (6037) 1988 EG, 1993 GD, and 2005 GL -- with Spitzer's IRAC instrument at 3.6, 4.5, 5.8, and 8.0 microns (the last object was detected only at 5.8 and 8.0 microns). The diameters of these three objects are 400 m, 180 m, and 160 m, respectively, with uncertainties of around 20% (including both observational and systematic errors). For all three the geometric albedos are around 0.30, in agreement with previous results that most NEOs are S-class asteroids. For the two objects detected at 3.6 and 4.5 microns, diameters and albedos based only on those data agree with the values based on modeling the data in all four bands. This agreement, and the high sensitivity of IRAC, show the promise of the Spitzer Warm Mission for determining the physical parameters for a large number of NEOs.

Thermal inertia determines the temperature distribution over the surface of an asteroid and therefore governs the magnitude the Yarkovsky effect. The latter causes gradual drifting of the orbits of km-sized asteroids and plays an important role in the delivery of near-Earth asteroids (NEAs) from the main belt and in the dynamical spreading of asteroid families. At present, very little is known about the thermal inertia of asteroids in the km size range. Here we show that the average thermal inertia of a sample of NEAs in the km-size range is 200 $\pm$ 40 J m&#8722;2 s&#8722;0.5 K&#8722;1. Furthermore, we identify a trend of increasing thermal inertia with decreasing asteroid diameter, D. This indicates that the dependence of the drift rate of the orbital semimajor axis on the size of asteroids due to the Yarkovsky effect is a more complex function than the generally adopted D^(&#8722;1) dependence, and that the size distribution of objects injected by Yarkovsky-driven orbital mobility into the NEA source regions is less skewed to smaller sizes than generally assumed. We discuss how this fact may help to explain the small difference in the slope of the size distribution of km-sized NEAs and main-belt asteroids.

For light curves sampled on an uneven grid of observation times, the shape of the power density spectrum (PDS) includes severe distortion effects due to the window function, and simulations of light curves are indispensable to recover the true PDS. We present an improved method for comparing light curves generated from a PDS model to the measured data and apply it to a 50-day long RXTE observations of NGC 4945, a Seyfert 2 galaxy with well-determined mass from megamaser observations. The improvements over previously reported investigations include the adjustment of the PDS model normalization for each simulated light curve in order to directly investigate how well the chosen PDS shape describes the source data. We furthermore implement a robust goodness-of-fit measure that does not depend on the form of the variable used to describe the power in the periodogram. We conclude that a knee-type function (smoothly broken power law) describes the data better than a simple power law; the best-fit break frequency is ~10^-6 Hz.