from the sizes of protoplanetary disks , we expect planets to form out to tens of au .
the many giant planets found on small orbits by radial velocity and transit techniques are believed to have formed beyond the water ice lines of their stars and migrated inward .
inward migration of massive planets can have a profound ( and usually destructive ) effect on the smaller objects in a planetary system , such as earth - sized planets within the habitable zones .
it is therefore important to determine how many systems retain their giant planets on wide orbits .
direct detection of these planets is limited to the most massive examples .
indirect detection through the use of debris disks addresses this limitation .
a debris disk consists of the circumstellar solid material that remains after the protoplanetary disk gas has dispersed and giant planets might have formed .
although most of the mass in a debris disk is harbored by the parent bodies ( planetesimals ) , dust generated in their collisions accounts for the majority of the disk s surface area , and observations of debris disks via their thermal emission or reflected stellar radiation trace this dust .
while tens of debris disks have been spatially resolved , the majority are detected only as an infrared excess in the spectral energy distribution ( sed ) of their star .
typical debris disk temperatures are tens to a few hundred kelvins , emitting as modified blackbodies that peak in the mid to far infrared .
these wavelengths are well - suited for study with the _ spitzer space telescope _ @xcite infrared spectrograph ( irs ;
@xcite ) and multiband imaging photometer for _ spitzer _ ( mips ; @xcite ) . for a recent review of debris disks ,
see @xcite .
debris disks often appear constrained to one or two discrete rings .
this is evident from images of resolved disks such as fomalhaut @xcite , and from seds of unresolved debris disks that are fit well by one or two blackbody functions , corresponding to dust at one or two distinct radial locations .
of the 28 disks without strong emission features presented by @xcite , all but one were fit better with a blackbody function than with a continuous disk model .
the exception was hr8799 , which was later determined to be best modelled by two blackbodies @xcite .
@xcite originally fit some excess seds with a power law , representing a continuous radial distribution of dust ; however , @xcite argued that these power law fits require the optical depth of the disk to increase with orbital radius ( which is theoretically implausible ) , and they found that these targets can be fit well by two blackbodies instead . these  warm " and  cold " debris disk components may be analogous to the asteroid belt and kuiper belt in the solar system .
why do rings form , and what sets their location ?
ice lines are one possibility . during the protoplanetary disk phase , a radial pressure gradient in the gas
partially counteracts the gravitational force on the gas , allowing it to rotate at a sub - keplerian velocity .
solid particles orbit at keplerian rates , and thus experience a head wind that slows their rotation and makes them spiral inwards . when these solids reach the ice line , the volatile component sublimates , producing a local pressure increase that counteracts the overall pressure gradient .
this creates a zone where the particles can settle without a headwind .
thus , there is a tendency to have a planetesimal belt near the ice line , and this belt can then produce grains that assume a specific temperature .
@xcite found similar warm disk temperatures ( 190 k ) around stars of different stellar types .
this can potentially be explained by the presence of the water ice line at 150 - 170 k. @xcite argue that ice lines of other species ( e.g. co , @xmath0 , @xmath1 ) at lower temperatures may play an important role in planet formation in the outer solar system .
if ice lines are also responsible for setting the location of the cold components , we would expect these components to have similar temperatures over a range of stellar types .
planets may also cause the discrete ring structures .
planetesimals will be scattered away once a planet is massive enough to dominate the gravity in its vicinity .
observations support the expected relation between planets and disk structure : the four giant planets imaged around hr 8799 @xcite are located in the gap between the warm and cold debris disk components @xcite ; the imaged planet orbiting @xmath2 pic appears to sculpt the inner edge of its debris disk @xcite ; the well - resolved cold debris ring around fomalhaut is likely confined by planets @xcite ; and in the solar system , neptune sculpts the inner edge of the kuiper belt @xcite .
if the locations of debris disks are set by planets , we can use cold debris disks to investigate wide - orbit planets , which are not easily studied by other means .
of the over 850 confirmed exoplanets , 35 have orbits larger than 5 au , and only 17 have orbits larger than 10 au ( nasa exoplanet archive ) .
this is likely due to the observational biases of the radial velocity and transit detection techniques .
direct imaging can detect wide - orbit planets , but the current technology is only sensitive to very massive planets around young , nearby stars .
uranus and neptune , for instance , could not be directly detected from outside the solar system .
the planet formation processes beyond 5 - 10 au are poorly understood .
for example , did the planets around hr 8799 form at their current locations or did they ( or at least their cores ) form on smaller orbits , then migrate outwards via scattering interactions with other planets or planetesimals ?
cold debris disks provide an observational test of these planet formation and migration theories .
planet formation by core accretion becomes less efficient farther from the star due to a radial increase in the dynamical timescale and decline in the surface density of solids in the protoplanetary disk @xcite .
therefore , the inner edge of a cold debris disk may represent the outer limit of efficient core accretion @xcite .
a planet migrating outwards into a debris disk will push the inner edge of the disk outward as well .
uranus and neptune may have formed on smaller orbits and migrated into the kuiper belt @xcite . in this alternate scenario ,
the inner edge of a cold debris disk may represent the limits of outward migration .
in this paper , we focus on _ spitzer _ measurements of cold debris disks , showing how their temperatures vary with the temperature of their central star and what this implies about planet formation and migration on wide - orbits .
first , we describe the selection of our target stars (  [ sec : targetselection ] ) .
then , we outline our photometric (  [ sec : photometry ] ) and irs (  [ sec : irsdatareduction ] ) data acquisition / reduction .
next , we detail our modeling of the stellar photosphere sed (  [ sec : photospheremodel ] ) , our derivation of the infrared excess , and our fitting of blackbodies to the excess (  [ sec : blackbodyfitting ] ) .
finally , we analyze the results (  [ sec : results ] ) and discuss the implications for wide - orbit planet formation and migration (  [ sec : discussion ] ) , before offering a summary and concluding remarks (  [ sec : conclusions ] ) .
we searched the _ spitzer _ observers log for main - sequence stars that were observed with the irs long low ( ll ) module ( both orders ) in staring mode and with mips at 24 @xmath3 and 70 @xmath3 , and we accumulated a sample of 546 targets .
the stellar properties of our target list are summarized in table [ table : stellarproperties ] . after reducing and analyzing the data , we refined the sample as described in
[ sec : blackbodyfitting ] .
this yielded 225 stars with significant excess ( of which 174 had cold components ) , which are detailed in table [ table : excessdata ] .
it is important to note that our sample comprised stars from a variety of _ spitzer _ observational programs , each having targets selected in a different manner ( e.g. nearby stars , stars with previously detected infrared excess , stars with rv - detected planets , etc . ) .
hence , our sample was not selected to be statistically representative in any one sense ; rather , it was designed to include as many relevant targets as possible .
we calculated the stellar temperature , @xmath4 , from the star s v-@xmath5 color using tabulated values for main - sequence stars ( or interpolations between these values ) from @xcite , and we estimated a 200 k one - sigma uncertainty on these values .
@xmath4 is listed for targets with significant excess in table [ table : excessdata ] . over the range of spectral types and metallicities of our sample of stars ,
v-@xmath5 is an accurate indicator of stellar effective temperature , with a peak - to - peak scatter of less than @xmath6 1.5% and weak metallicity dependence @xcite .
we use @xmath4 to parametrize stellar type , and this degree of accuracy is sufficient for our purpose . some of the uncertainty in @xmath4 reflects the effect of interstellar reddening on v-@xmath5 , although the great majority of stars in this sample are within the local bubble and therefore should not be strongly reddened ( of the 174 targets for which we detect cold components , only 29 are beyond 100 pc , 19 are beyond 120 pc , and 4 are beyond 150 pc )
.    ages for these stars were estimated from a combination of chromospheric activity measurements , x - ray emission , placement on the hr diagram , surface gravity , membership in clusters and associations , and gyrochronology collected from the literature .
@xcite describe the intercomparison of these methods and how they are applied , and the quoted ages are on the scale calibrated by @xcite .
the ages are listed in table [ table : stellarproperties ] , along with references for the measurements used to derive the ages and a quality flag for the age accuracy ( ranging from 0 if no age could be determined to 3 if there were three or more age measurements with good agreement ) .
cccccccccc hip110 & hd224873 & & g5 & 49.6 & 0.42 & 2 & 2,20 & 5346816 & not included + hip345 & hd225200 & hr9102 & a0v & 124.8 & 0.1 & 2 & 31 & 12720128 & warm + cold + hip394 & hd225239 & hr9107 & g2v & 39.2 & 0.25 & 1 & 2 & 4028672 & not included + hip490 & hd105 & & g0v & 39.4 & 0.17 & 3 & 2,3,4,5,6,9,12 & 5295616 & cold + hip522 & hd142 & hr6 & f7v & 25.7 & 4.8 & 2 & 3,17,27 & 22296064 & not included + hip544 & hd166 & hr8 & k0ve & 13.7 & 0.24 & 3 & 1,2,6,15,23 & 12717824 & warm + cold + hip560 & hd203 & hr9 & f3vn & 39.4 & 0.01 & 2 & 38 & 14983424 & cold + hip682 &
hd377 & & g2v & 39.1 & 0.22 & 3 & 2,4,9,12,20,32 & 5268736 & warm + cold + hip910 & hd693 & hr33 & f8vfe-08h-05 & 18.7 & 3 & 3 & 3,4,23,24 & 14494720 & not included + & & & & & & & & 4008448 & + hip919 & hd691 & & k0v & 34.2 & 0.24 & 3 & 2,9,12,14,18,20,39 & 5345280 & not included + hip1031 & hd870 & & k0v & 20.2 & 2.3 & 2 & 3,5 & 16026112 & cold + hip1134 & hd984 & & f5 & 47.1 & 0.25 & 3 & 2,4,9,20 & 5271808 & not included + hip1292 & hd1237 & & g8.5vk : & 17.5 & 0.3 & 3 & 2,3,4,5,6 & 22290176 & not included + hip1368 & & & k7 & 15 & 0.9 & 3 & 1,2,39 & 25674240 & cold + hip1473 & hd1404 & hr68 & a2v & 41.3 & 0.45 & 1 & 7 & 14160384 & warm +      although we built our target list around the _ spitzer _ irs spectra , we supported these spectra in our analysis with a suite of photometric data .
the properties of the photometric systems we used are summarized in table [ table : photometrybands ] and the photometric data for the targets with significant excess ( see  [ sec : blackbodyfitting ] ) are given in table [ table : excessdata ] .
mips photometry at 24 @xmath3 provided an important calibration reference for the irs data (  [ sec : irsdatareduction ] ) , while 70 @xmath3 photometry provided a crucial constraint on the temperature of cold debris disks (  [ sec : blackbodyfitting ] ) .
we used our in - house debris disk pipeline to reduce and extract photometry for the mips data as part of the effort to preserve the legacy of _ spitzer _ measurements on debris disk studies ( @xcite ; k. su et al . 2013 , in prep ) .
basic reduction ( up to the post - bcd mosaics ) and calibrations of the mips data follow the descriptions by @xcite and @xcite .
in addition , the extraction of the 24 @xmath3 photometry was briefly described in @xcite , where the source position at 24 @xmath3 is determined by a combination of psf fitting and 2d gaussian fitting .
both psf fitting and aperture photometry aperture photometry used an aperture radius of 6.255@xmath7 , a sky annulus of 19.92 - 29.88@xmath7 , and an aperture correction factor of 1.6994 .
] were performed , and we preferentially used the aperture photometry at 24 @xmath3 because this was used to calibrate the photosphere model at 24 @xmath3 , as described in  [ sec : photospheremodel ] and the appendix ( the psf and aperture photometry agree to within a few percent ) .
we then used the 24 @xmath3 source position to perform psf fitting for data at 70 @xmath3 by minimizing the residual signal at the source position ( for details , see k. su et al .
2013 , in prep ) . for faint sources located in areas with structured background ,
the resultant 70 @xmath3 photometry can be negative , which reflects non - detection .
the 1@xmath8 photometry uncertainty ( listed in table [ table : excessdata ] ) includes the pixel - to - pixel variation near the source of interest , and detector repeatability ( 1% and 5% of the source flux at 24 and 70 @xmath3 , respectively ) .
mips photometry data for all of our targets can be found in @xcite , @xcite , and k. su et al .
( 2013 , in prep ) .
v , j , h , and @xmath5 photometry were used to model the stellar photosphere sed (  [ sec : photospheremodel ] ) and to estimate @xmath4 (  [ sec : targetselection ] ) .
we obtained _ hipparcos _ v band @xcite and 2mass j , h , and k band photometry @xcite from the vizier online database .
many of the targets in the sample are nearby stars and are severely saturated in the 2mass data . to overcome this , we used heritage aperture photometry , which we transformed to match the 2mass system .
when both 2mass and transformed heritage photometry of high quality were available , we averaged them .
the references for the heritage photometry are in table [ table : excessdata ] .
irac photometry ( 3.6 @xmath3 ) was also used to model the stellar photosphere sed (  [ sec : photospheremodel ] ) .
these data were taken in _ spitzer _
cycle 7 ( pid 70076 , pi : su ) .
all the irac data were taken in subarray mode with four dithered positions to avoid saturation .
we used the basic calibrated data ( bcd ) products provided by the _
science center ( pipeline version s18.18 ) , and performed the necessary steps ( pixel solid angle correction and pixel phase correction ; k. su et al .
2013 , in prep ) to extract the photometry .
aperture photometry was used for each individual data frame ( 64 frames per dithered position ) ; and the final quoted flux is the median value of all measurements per star .
wise data @xcite were obtained from vizier , but are not included in table [ table : excessdata ] . because wise photometry is less accurate than _ spitzer _ photometry , we did not use it quantitatively .
instead , we used it as a qualitative confirmation of our _ spitzer _ data .
our irs reduction started with the level 1 bcd products , downloaded from the _ spitzer _ heritage archive .
in addition to irs ll data , we also reduced short low ( sl ) module data , when available , and included them in our analysis .
the astronomical observation requests ( aors ) used for each target are listed in table [ table : stellarproperties ] .
the basic reduction was performed using the spectroscopic modeling analysis and reduction tool ( smart ) software package @xcite , scripted into a series of automated routines . for each spectral order of each irs aor , three files ( 2d spectra , uncertainty , and mask )
were combined into a single 3-plane file .
bad and rogue pixels were removed by the routine irsclean with the clean parameter set to 4096 ( pixels with this value or higher were included in the clean ) .
next , when available , multiple data collection events ( dces ) for the same nod position were combined , and then the background was removed from each 2d spectrum by subtraction of the opposite nod .
the 2d spectra were then converted into 1d spectra using smart s optimal 2 nod extraction @xcite .
the 1d spectra from both nods were combined .
the result was a wavelength , flux , and uncertainty vector of each spectral order for each aor .
the  bonus " third order data were not used .
the remainder of the data reduction and analysis was performed with the matlab software package .
irs data at wavelengths near the order edges often were of poor quality , so data longward of 38 @xmath3 and shortward of 20.75 @xmath3 were discarded from the ll first order , data shortward of 14.3 @xmath3 were discarded from the ll second order , data longward of 14.7 @xmath3 and shortward of 7.55 @xmath3 were discarded from the sl first order , and data shortward of 5.25 @xmath3 were discarded from the sl second order .
systematic offsets in flux existed between irs orders , which we fixed by applying a multiplicative correction factor to the ll first order spectrum and to both sl order spectra ( if they existed ) , in order to align them with the ll second order flux .
initially , we determined the order correction factors using an automated routine , but due to the variety of shapes of the spectra and the presence of outlying data points , we found that fine - tuning these factors by eye was more reliable .
the data from all available orders and modules were then combined into a single spectrum .
next , we cut outlying data points from the spectrum in an iterative process .
the spectrum was fit with a polynomial and the standard deviation of residual flux values around the fit was calculated .
points lying more than three standard deviations from the fit were discarded .
this process was iterated six times .
because the scatter of the data generally increased towards longer wavelengths , we applied this process separately to the short and long ends of our spectrum ; if sl data were available , the two sections were divided at 14 @xmath3 and each section fit with a fourth degree polynomial , whereas if no sl data were available , the sections were divided at 25 @xmath3 and each section was fit with a second degree polynomial .
this procedure , on average , cut 10 data points from each spectrum , with the first iteration typically cutting five points , and the sixth iteration only cutting zero , one , or two points ( nearly 90% of spectra had no points cut in the sixth iteration ) . before cutting ,
the spectra had 181 , 296 , or 373 points , depending on the available orders .
after the outliers were cut , the spectrum was smoothed by binning to a wavelength resolution of 0.7 @xmath3 .
the absolute calibration of irs data is only accurate to @xmath910% , whereas the mips calibration is accurate to @xmath92% @xcite . to reduce systematic errors in our irs data
, we multiplicatively scaled our spectrum to be consistent with the measured mips 24 @xmath3 flux density for each target .
the irs spectra were calibrated as point sources , so if a source was slightly extended , it would have an incorrect slit - loss correction .
the stellar photospheres ( point - like ) are dominant in the irs spectra for the majority of sources in the sample ; therefore , this has no significant impact on the result .
if the mips photometry for a source was contaminated by background emission , this contamination would be passed to the irs spectrum of the source through this scaling ( the irs data may or may not have picked up the contamination , depending on the slit orientation ) .
nevertheless , an additional multiplicative factor was later applied to the irs data ( see  [ sec : photospheremodel ] ) , which corrected any lingering errors from this process .
some targets were observed with more than one irs aor .
we combined the spectra from these aors , interleaving their data points , then re - smoothing to a spectral resolution of 0.7 @xmath3 .
although the spectrum of a main - sequence star peaks in the visible or near - ir , the stellar photosphere can still be the dominant source of flux density in the mid - ir .
thus , to extract the thermal emission of the debris disk ( the infrared excess ) , an accurate model spectrum of the stellar photosphere must be generated and subtracted from the data .
it is common to use photosphere spectra from detailed numerical models of stellar structures , such as kurucz / atlas9 @xcite .
however , these models are not well - tested in the mid- and far - infrared .
@xcite compared various families of models ( kurucz / atlas9 , marcs , and nextgen / phoenix ) and found that the choice of model family can affect whether a star is determined to have an infrared excess or not .
we opted for a simpler model , a rayleigh - jeans relation given by @xmath10 , where rj is an amplitude scale factor .
this model is appropriate for wavelengths beyond 5 @xmath11 m because stars in this study , mostly a through k dwarfs , have virtually no spectral features in the mid - ir , behaving as blackbodies in the rayleigh - jeans regime . for example , in the range of most interest for this study , the departure of a reference a0 star spectrum @xcite from a rayleigh - jeans fit between 15 and 30 @xmath3 is no more than @xmath6 0.3% , and the atlas9 @xcite solar model follows rayleigh - jeans behavior to within @xmath6 0.24% from 20 to 40 @xmath3 .
however , the theoretical spectra show minor systematic differences from the observations @xcite , so use of them to improve the knowledge of the seds would be risky .    modeling
the photosphere thus came down to finding the appropriate scale factor ( rj ) for each star in our sample .
we used a set of color relations ( derived from a sample of stars known to have no ir excess ) to predict the mips 24 @xmath3 magnitude of the photosphere , [ 24 ] , from measured v , j , h , @xmath5 , and irac magnitudes
bands where ir excesses due to very close - in dust are very rare .
the details of these relations are described in the appendix .
this procedure estimates the 24 @xmath11 m photospheric outputs to about 2% rms , based solely on data at wavelengths short of 4 @xmath11 m .
we then converted [ 24 ] to flux density units and used it to find the scaling factor for the photosphere model , which is given in table [ table : excessdata ] for targets with significant excess ( see  [ sec : blackbodyfitting ] ) .
we finally applied a small multiplicative factor to the irs data if it clearly looked slightly offset from the photosphere model .
then , we subtracted the photosphere model from the irs spectra and mips 70 @xmath3 point to obtain the infrared excess , @xmath12 .
the photosphere model was assumed to have an uncertainty of 2% , thus the uncertainty in the excess flux density was @xmath13 ^ 2}.\ ] ]       data are omitted from the right panels , as these points were not used to constrain the fits . ]    to interpret the excesses in a general way , we needed to assign them fiducial temperatures .
the emission of disk grains is a function of their optical constants ; however , as discussed in
[ sec : introduction ] , debris disks radiate like blackbodies at one or two temperatures .
the simplest possible description of the emission is in terms of these temperatures . in this section ,
we describe how we determined if each target had significant excess and whether the excess was best described by one or two blackbodies .
we fit @xmath14 with a one component blackbody model @xmath15 where @xmath16 is the planck function .
the best - fit parameters , @xmath17 and @xmath18 were found by minimizing the reduced chi - squared , @xmath19 ^ 2 } { \sigma_\text{excess}\left(\lambda_i\right)^2 } \nonumber\\ & + 28\frac { \left[f_{\nu\,\text{excess}}\left(70 \micron\right)-f_{\nu\,\text{model}}\left(70 \micron\right)\right]^2 } { \sigma_\text{excess}\left(70 \micron\right)^2 } \bigg\},\end{aligned}\ ] ] where
@xmath20 , n is the number of data points in the irs excess data , and n is the number of free parameters in the fit ( n=2 ) .
the mips 70 @xmath3 data point was weighted in the fit as 28 irs data points , equivalent to the number of irs data points ( 0.7 @xmath3 resolution ) that would fit inside the equivalent width of the mips 70 @xmath3 spectral response function ( 19.65 @xmath3 ) .
the fit was performed using matlab s lsqcurvefit algorithm .
@xmath18 was constrained to between 0 and 500 k.    while @xmath17 represents the amplitude of the debris disk emission , a more useful measure of a debris disk s brightness is the fractional excess , @xmath21 .
we calculated @xmath22 according to @xmath23 from equation 2 of @xcite . @xmath24 and @xmath25 are the peak flux density values of the disk emission and stellar photosphere emission , respectively , which occur at wavelengths @xmath26 and @xmath27 .
while the peak flux densities for the disk components were easily calculated from our best fit model , the rayleigh - jeans stellar photosphere model had no maximum .
to overcome this , we created a blackbody function using the temperature of the star from table [ table : stellarproperties ] , and scaled it to match the flux density of our rayleigh - jeans model at 24 @xmath3 .
from this representation of the photosphere , we found the maximum flux density , and calculated @xmath22 .
data are omitted from the right panels , as these points were not used to constrain the fits . ]    with our best fits in hand , we refined our sample to only those targets with statistically significant excess . first , we inspected all fits by eye and discarded targets with poor quality data that resulted in clearly unrealistic fits .
second , we discarded targets with excess that was too faint , @xmath28 .
third , we inspected the fits and identified all targets whose excess relied solely on the mips 70 @xmath3 data point ( i.e. there was no excess in the irs data ) .
for these cases , we required this mips point to represent a significant excess , so we discarded targets where @xmath29 or where @xmath30 .
this process resulted in 321 of the original 546 stars having no significant excess .    of the remaining 225 targets with significant excess
, we next determined whether the excess consisted of one or two components . to do this
, we fit the excess sed of each target with a model consisting of the sum of two blackbodies , @xmath31 finding the optimal set of the four parameters @xmath32 , @xmath33 , @xmath34 , and @xmath35 was again done by minimizing the reduced chi - squared ( equation [ eq : reducedchisq ] ) , now with n=4 .
data are omitted from the right panels , as these points were not used to constrain the fits . ]
our definition of  cold " was the coldest well - detected component of the excess that was below 130 k. @xcite found a continuous distribution of warm components above @xmath9130 k , centered around 190 k , whose temperatures were set by the water ice line .
thus , any component above 130 k was considered
warm " for the purposes of this study .
a component with a temperature of 110 k , for example , would be considered warm if another , colder temperature component was also detected , but would be considered cold if no colder component was detected . to implement this
, we fit each target once using 100 k as the division between warm and cold , and again with the division at 130 k. we then selected the better of these two cases ( based on reduced chi - squared ) to represent the best two - component model . for many targets , these two cases yielded identical fits . in both cases ,
the minimum cold temperature allowed by the fit was 0 k , and the maximum allowed warm temperature was 500 k.    deciding if each target warranted a two - component model was a two - step process .
first , we required that the reduced chi - squared of the two - component fit be at least three times better than that of the one - component fit ( i.e. all targets where @xmath36/@xmath37 were deemed to be better fit by a single component ) .
second , for the remaining targets , we identified those for which the presence of the cold component relied entirely on the mips 70 @xmath3 point , meaning the irs excess was entirely fit by the warm component .
for these targets to have significant cold components , we required their mips 70 @xmath3 excess to fall more than 3@xmath8 above the excess predicted by the warm component .
that is , if one of these targets had @xmath38 we concluded that it was best fit by one component .
this step was analogous to the third step we performed when deciding if each target had excess or not .
the targets that passed both of these criteria were deemed to have two components .
this process identified 100 single - component cold disks , 51 single - component warm disks , and 74 two - component disks .
we found only one target ( hip45585 ) that required two warm components ( both @xmath39130 k ) to fit properly , and we discarded this target from our results as it had no cold components ( it is counted in the 321 targets that we discarded ) .
we found no targets requiring two components colder than 100 k. our methods were insensitive to possible additional cold components below @xmath950 k. we concluded that our fitting procedure allowed us to find all cold components above this limit .
equivalently , our method reliably fit the inner edges of the cold disks .
for the two - component disks , we calculated the fractional excesses for the warm and cold components , @xmath40 and @xmath41 , using equation [ eq : fwyatt2 ] .
one - component excesses were deemed warm or cold depending if they had temperatures above or below 130 k. the verdicts for all targets are given in table [ table : stellarproperties ] and the parameters of our best fits for targets with significantly - detected components are listed in table [ table : excessdata ] .
we estimated that the uncertainty in @xmath33 was 10 k. this estimate was conservative ; warmer cold components had more irs data points in excess of the photosphere , so their fits were more constrained , with a temperature uncertainty of 3 to 7 k.    we were unable to accurately constrain the temperatures of very cold components , which were detected only at 70 @xmath11 m . with only one data point , blackbodies at a wide range of cold temperatures could be given the amplitude necessary to match the point , so a degeneracy existed between t and @xmath22 in blackbody fits to the excess sed .
we made this distinction at 45 k ; above this temperature there was enough information in the irs spectra to constrain the temperature .
cold temperatures above 45 k were considered true detections , while cold temperatures below 45 k were considered unconstrained and were replaced with upper limits at 45 k. by using the 70 @xmath3 data to determine if such targets had significant excesses , we ensured that these disks were truly very cold , rather than warm but very faint .
of the 174 cold components in our sample , 25 had temperature upper limits .
some irs spectra exhibit mineralogical emission features ( e.g. hip41081 , hip57971 ) .
our model did not explicitly account for these features , but they peak sharply at @xmath910 @xmath3 and do not resemble blackbody functions .
these features did not influence our fits , except to raise the value of the minimum reduced chi - squared .
the data and best fits for a small sample of our targets are shown in figures [ fig : twocompgallery ] , [ fig : coldcompgallery ] , and [ fig : warmcompgallery ] ( for targets with two components , one cold component , and one warm component , respectively ) .
the left panels in all of these figures show the measured data with the photosphere and blackbody models , while the right panels illustrate the excess data and models with the photosphere subtracted .
histograms of the significant warm and cold debris disk temperatures are shown in figure [ fig : temphistogram ] .
45 k around cooler stars .
the green line is the best fit trend to the data , determined by a bayesian linear regression , @xmath42 .
a representative error bar is in the lower right . ]
@xcite studied the behavior of the warm disk components in detail , but their sample ( restricted to stars with ages less than 1 gyr ) had too few 70 @xmath3 detections of solar - type stars ( 9 ) to determine any trends between the cold component temperature and stellar type .
to look for such a trend in our larger sample , we plot the cold component temperature versus stellar temperature in figure [ fig : tdisktstarscatter ] .
well - constrained disk temperatures are black circles ( open circles are young systems with age less than 25 myr ) , and upper limits are downward facing blue triangles .
although there is substantial scatter in the cold component temperatures , a positive correlation between cold debris disk temperature and stellar temperature is evident .    to quantify this trend ,
we fit a linear relation to the data by bayesian analysis using the routine linmix_err .
the routine properly handles upper - limits , and it uses uncertainties in both x and y directions ( we assumed an uncertainty of 10 k in @xmath33 and 200 k in @xmath4 ) . in its model , the routine also includes the normally - distributed intrinsic scatter of the data around the trend .
the posterior distributions of the slope , intercept , and intrinsic scatter from the bayesian regression are show in figure [ fig : trendposteriors ] .
the slope of the linear regression was @xmath43 , thus the existence of a trend between @xmath33 and @xmath4 is significant at level greater than 4.5@xmath8 .
the 1@xmath8 intrinsic scatter around the trend was @xmath44 k. the best - fit trend line from the bayesian analysis , @xmath42 , is plotted in green in figure [ fig : tdisktstarscatter ] . the choice of a linear fit is not physically motivated ; its purpose is to show that  to first order  there is a correlation between cold debris disk temperature and stellar type .    a significant source of scatter around the trend results from the 10 targets with @xmath45 k and @xmath46 k. six of these stars ( hip560 , hip64184 , hip64995 , hip65875 , hip67497 , and hip78663 ) are less than 25 myr old ( most are part of the scorpius - centaurus association ) .
although we show in  [ sec : discussion ] that there is no significant trend in @xmath33 with age across the broad range of stellar ages in our sample , the excesses of these targets may be a product of their relatively young ages .
they represent a period ( 5 - 50 myr ) when considerable collisional activity may still be occurring as a result of planet building @xcite .     vs. @xmath4 trend slope ,
intercept , and intrinsic scatter , generated by the linmix_err routine . ]
we now consider what determines the temperature of cold debris disks and how a trend with stellar type might arise .    as discussed in  1 , @xcite found similar warm disk temperatures ( 190 k ) around stars of different stellar types , and they attributed the effect to a particle trap at the water ice line .
if ice lines of other species set the location of cold debris disks , we would expect cold disks around different type stars to have a common temperature .
however , the trend of cold component temperature with stellar type ( figure [ fig : tdisktstarscatter ] ) is inconsistent with any strictly temperature - dependent mechanism for setting the location of cold debris disks .
@xcite show that during the protoplanetary disk phase , the water ice line location is time dependent , and predicting it requires accounting for viscous heating in the disk , as well as the evolution of the pre - main - sequence stellar luminosity .
such considerations have not been applied to potential ice lines of species other than water , and we do not consider this level of detail here
.    delayed stirring , in which the radial location of dust in a debris disk moves outwards with time , is a postulated mechanism for debris disk evolution @xcite .
this could occur if the parent bodies are distributed in a broad ring .
particles in a debris disk collide and grind into dust faster on smaller orbits where the dynamical timescale is shorter .
thus , the location of the emitting dust moves outward with time , and , therefore , becomes cooler with time . because late type stars have longer lifetimes than early type stars , the late type stars in our sample are generally older than the early type stars .
if delayed stirring occurs , then this age bias in our sample could explain the observed trend in @xmath33 with @xmath4 . to test this hypothesis , we plot @xmath33 against the age of the system ( if available ) for targets in three @xmath4 bins ( 5000 to 6000 k , 6000 to 7000 k , and 7500 to 9500 k ) , shown in figure [ fig : agevstdisk ] .
we see no trend in @xmath33 with age in any bin , suggesting that delayed stirring does not produce the trend of disk temperature with stellar type .
is plotted against the age of the system for targets in three @xmath4 bins .
no trend is seen with age , suggesting age does not play a confounding role in our discovered @xmath33 vs. @xmath4 trend .
this also argues against the occurrence of delayed stirring . ]
perhaps cold debris disks are all at roughly the same orbital distance , regardless of stellar type . assuming the disk is heated to its equilibrium temperature , the relation between the disk s temperature , location , and stellar type is given by @xmath47 the second and third lines of equation [ eq : equilibriumtemp ] are derived assuming @xmath48 and @xmath49 ( valid in the roughly solar mass range ) .
so with @xmath50 constant , early type stars would host warmer disks .
@xcite perform detailed simulations of the evolution of debris disks with inner and outer edges at 30 and 150 au , respectively , around stars with mass ranging from 1 to 3 @xmath51 .
their models output the disk emission at 24 and 70 @xmath3 , which show that the ( color ) temperature of these cold debris disks does increase with stellar mass , in agreement with this expectation .
but is the orbital location of cold debris disks truly constant with stellar type ? if the size of cold debris disks traces the size of the original protoplanetary disks , we can use observations of the latter to address this question .
@xcite find that protoplanetary disk mass and radius are related by @xmath52 furthermore , observations @xcite show that @xmath53 combining these relations yields @xmath54 so protoplanetary disk size does increase with earlier stellar type , albeit slowly .
how quickly must the disk size increase for it to maintain a constant temperature ? from equation [ eq : equilibriumtemp ] , @xmath55 although the size of disks does increase with earlier spectral type , it does so more slowly than required to maintain a constant temperature , thus disks are expected to be warmer around early type stars , consistent with our findings . substituting equation [ eq : ppradius ] into equation [ eq : equilibriumtemp ]
reveals how the disk temperature would vary with spectral type in this case : @xmath56 new results from @xcite suggest that equation [ eq : ppandstarmass ] may hold only for @xmath57 , with disk mass constant or possibly even decreasing as @xmath58 for higher mass stars .
if this were true and translated into a flat or decreasing @xmath50 with @xmath59 , then the disk temperature would increase even faster with spectral type than derived in equation [ eq : migresult ] .
is it plausible that the inner edge of cold debris disks scales with the size of the original protoplanetary disk ?
some mechanism must clear the material inside cold debris disks and set their inner edge , and planets are a common explanation .
an outwardly migrating planet would set an inner edge , but migration from planet - planet scattering can be a chaotic and unpredictable phenomenon ; simulations show that the planet s final location depends sensitively on the initial conditions of the system @xcite .
this suggests that a trend with stellar temperature would not arise .
outward migration via scattering through a smooth disk of planetesimals would proceed in a more orderly manner , and the planet would halt its migration when the surface density of planetesimals decreased below a certain threshold , at a location that scales with the size of the original protoplanetary disk .
planet formation by core accretion ( without migration ) would also clear debris inside the cold component and create its inner edge .
planetesimals would be either incorporated into the planets as they formed , or scattered away once the planets became massive enough to dominate the gravitational field in their vicinity . as mentioned in
[ sec : introduction ] , core accretion efficiency declines with increasing orbital radius , leaving an outer zone of planetesimals beyond the planets . the timescale for the formation of a planet with a given mass scales as @xmath60 where p is the orbital period and @xmath61 is the surface density of solids . substituting kepler s third law , @xmath62 ( @xmath63 is the orbital distance ) and the typical protoplanetary disk structure , @xmath64 @xcite into equation [ eq : formationtime ] gives @xmath65 after a time ( @xmath66 ) planets will have formed out to a given orbital location ( @xmath63 ) , which sets the inner edge of the cold debris disk .
so setting @xmath66 constant and @xmath67 predicts that the size of the disk would scale with spectral type as @xmath68 by comparing this with equation [ eq : equilibriumradius ] , we see that in this case as well , the size of the disk grows more slowly with stellar type than required to maintain a constant equilibrium temperature , consistent with our observational results .
substituting equation [ eq : formradius ] into equation [ eq : equilibriumtemp ] shows how the disk temperature would vary with stellar type if its inner edge were set by the limits of planet formation : @xmath69
we studied the circumstellar environment of 546 main - sequence stars via the mid infrared emission of their debris disks , as measured by the _ spitzer _ space telescope . after subtracting a model of the flux expected from the stellar photosphere
, we obtained an sed of the infrared excess for each target .
we found 225 targets with significant excess : 100 with a single cold component , 51 with a single warm component , and 74 with two components . examining the results revealed a trend between the temperature of the inner edge of the cold debris disk component and that of its host star .
photosphere magnitude for a sample of 1037 stars .
the blue curve shows the gaussian fit to the distribution with @xmath70 and @xmath71 . ]
this trend is inconsistent with theories that predict the location of cold debris disks to be strictly temperature - dependent , i.e. we rule out the dominance of ice lines in sculpting the outer regions of planetary systems .
we also rule out delayed stirring as the source of this trend .
the trend can potentially be explained if the outward migration of planets traces the extent of the primordial protoplanetary disk , which tends to be limited to warmer equilibrium temperatures for hotter stars .
the trend can also be explained if planets form in situ out to a distance where the core accretion efficiency drops below a certain threshold , leaving a cold debris disk that is warmer around earlier type stars .
we thank scott kenyon and benjamin bromley for helpful science discussions , and the referee for many useful comments .
we also thank vianney lebouteiller for help in scripting the irs reduction pipeline .
this work was supported by contract 1255094 from caltech / jpl to the university of arizona .
we used the simbad database and the vizier catalogue access tool , operated at cds , strasbourg , france .
we also used the nasa / ipac infrared science archive , which is operated by the jet propulsion laboratory , california institute of technology , under contract with nasa .
this work is based on observations made with the _ spitzer _ space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa .
and we used data products from the two micron all sky survey , which is a joint project of the university of massachusetts and the infrared processing and analysis center / california institute of technology , funded by nasa and the nsf .
finally , we used data products from the wide - field infrared survey explorer , which is a joint project of the university of california , los angeles , and the jet propulsion laboratory / california institute of technology , funded by nasa .
here we describe how we used photometry to derive a model photosphere magnitude at 24 @xmath3 .
the relation among v , @xmath5 , and [ irac ] used to predict [ 24 ] is equation [ eq : mips24fromivk ] , but a number of other relations were used beforehand to improve the accuracy of the result .
first , we preprocessed @xmath5 , h , and j photometry to remove small nonlinearities in the relations .
these small corrections to the 2mass photometry were derived by comparison with our all - sky warm mission measurements in irac band 1 .
the preprocessing replaced @xmath5 with @xmath72 according to @xmath73 h with h according to @xmath74 and j with j according to @xmath75    equation [ eq : mips24fromivk ] utilizes @xmath76 , a combination of measured and derived @xmath5 magnitudes .
one of these , @xmath77 , was derived by solving the relation among v , @xmath77 , and j , given by @xmath78 where @xmath79 .
this relation is a fit to the behavior of more than 1000 main - sequence stars , as are the other relations that are given below .
@xmath80 was derived ( if h was measured ) by solving the relation among v , @xmath80 , and h , given by @xmath81 where @xmath82 .
@xmath76 was then calculated by averaging @xmath72 , @xmath77 , and @xmath80 with relative weights of 1 , 0.75 , and 0.47 , respectively .
the weights were determined from typical 2mass errors and by minimizing the residuals in comparing the projected photospheric levels with the 24 @xmath11 m measurements .
equation [ eq : mips24fromivk ] also uses @xmath83_\text{super}$ ] , a combination of measured and derived irac magnitudes .
the derived irac magnitude , @xmath83_1 $ ] , was calculated by solving @xmath84_1 = k_{s\,\text{super } } + 0.0012x^4 - 0.01448x^3 + 0.0419x^2 - 0.056x + 0.0295,\ ] ] where @xmath85 . the measured [ irac ] and derived @xmath83_1 $ ] were averaged together with relative weights of 1 and 0.5 , respectively , yielding @xmath83_\text{super}$ ] .
finally , [ 24 ] was calculated with @xmath86 = [ \text{irac}]_\text{super } + 0.001x^5 - 0.0097x^4 + 0.0179x^3 + 0.0305x^2 - 0.0821x + 0.0312,\ ] ] where , again , @xmath85 .
this process required v and at least one of j , h , or @xmath5 to proceed . targets that did not have this minimum photometry available were not included in our sample , although nearly all targets in our sample had a complete suite of v , j , h , and @xmath5 photometry available ( table [ table : excessdata ] ) .    the accuracy and precision of this method is illustrated in figure [ fig:24accuracy ] , which shows the distribution of the difference between the observed and predicted [ 24 ] for over 1000 stars .
a gaussian fit to this distribution ( the blue curve ) has a mean of 0.00158 and a standard deviation of 0.0224 .
cccccccccccccccc hip345 & 6.39 & 6.28 & 6.25 & 6.26 & & 8843 & & 35.7 @xmath6 0.38 & 97.24 @xmath6 5.35 & 11.9 & 94.9 @xmath6 5.35 & 53 & 185 & 9.42 & 4.56 + hip490 & 7.53 & 6.46 & 6.19 & 6.12 & & 5923 & & 28.39 @xmath6 0.3 & 152.7 @xmath6 9.71 & 14.2 & 149.9 @xmath6 9.71 & 48 & & 42 & + hip544 & 6.13 & 4.73 & 4.63 & 4.31 & & 5405 & & 159.1 @xmath6 1.58 & 105.8 @xmath6 6.32 & 72.6 & 91.57 @xmath6 6.33 & 50 & 126 & 4.99 & 4.78 + hip560 & 6.19 & 5.45 & 5.33 & 5.24 & & 6635 & & 117.2 @xmath6 1.12 & 67.08 @xmath6 7.06 & 31.6 & 60.89 @xmath6 7.06 & 129 & & 18.4 & + hip682 & 7.59 & 6.42 & 6.15 & 6.12 & & 5809 & & 36.49 @xmath6 0.39 & 170.6 @xmath6 10.62 & 14.7 & 167.7 @xmath6 10.6 & @xmath8745 & 119 & 47.3 & 11.6 + hip1031 & 7.23 & 5.85 & 5.47 & 5.38 & & 5372 & 1.9 @xmath6 0.03 & 50.67 @xmath6 0.51 & 22.51 @xmath6 4.98 & 27.8 & 17.07 @xmath6 4.99 & 50 & & 3.29 & + hip1368 & 9 & 6.38 & 5.75 & 5.58 & & 3903 & 1.7 @xmath6 0.03 & 46.6 @xmath6 0.48 & 20.74 @xmath6 5.15 & 26.5 & 15.55 @xmath6 5.15 & @xmath8745 & & 9.99 & + hip1473 & 4.52 & 4.34 & 4.42 & 4.46 & & 9164 & & 154.4 @xmath6 1.56 & 43.82 @xmath6 6.5 & 65.7 & 30.94 @xmath6 6.51 & & 133 & & 1.38 + hip1481 & 7.46 & 6.46 & 6.25 & 6.15 & & 6258 & & 34.76 @xmath6 0.36 & -0.83 @xmath6 2.98 & 13.7 & -3.518 @xmath6 2.98 & & 217 & & 11.9 + hip1499 & 6.46 & 5.33 & 5.04 & 4.89 & & 5692 & 2.95 @xmath6 0.03 & 80.94 @xmath6 0.81 & 61.43 @xmath6 6.81 & 43.1 & 52.98 @xmath6 6.81 & 54 & & 5.54 & + hip2072 & 3.9 & 3.74 & 3.58 & 3.51 & & 8131 & & 311.2 @xmath6 3.1 & 72.97 @xmath6 4.43 & 143 & 44.94 @xmath6 4.47 & @xmath8745 & 169 & 0.577 & 1.47 + hip2472 & 4.77 & 4.67 & 4.77 & 4.7 & & 9096 & & 112.5 @xmath6 1.12 & 76.98 @xmath6 6.54 & 49.5 & 67.27 @xmath6 6.54 & 69 & 192 & 1.47 & 1.5 + hip2578 & 5.07 & 5.06 & 5.16 & 4.99 & & 9028 & & 232 @xmath6 2.31 & 56.48 @xmath6 4.13 & 35.5 & 49.52
@xmath6 4.13 & & 194 & & 15.6 + hip2710 & 6.91 & 6.04 & 5.85 & 5.75 & & 6473 & 1.33 @xmath6 0.02 & 40.62 @xmath6 0.44 & 104.7 @xmath6 6.5 & 19.8 & 100.8 @xmath6 6.5 & @xmath8745 & 110 & 16.5 & 2.94 + hip2843 & 6.71 & 5.84 & 5.66 & 5.59 & & 6510 & 1.57 @xmath6 0.03 & 42.13 @xmath6 0.46 & 21.54 @xmath6 3.76 & 23.3 & 16.98 @xmath6 3.76 & 64 & & 2.43 & +    ccc v & 0.5423 & 3729 + 2mass j & 1.241 & 1623 + 2mass h & 1.6513 & 1075 + 2mass @xmath5 & 2.1657 & 676 + wise 1 & 3.35 & 309.54 + wise 2 & 4.60 & 171.79 + wise 3 & 11.56 & 31.67 + wise 4 & 22.09 & 8.36 + irac 1 & 3.6 & 269.53 + mips 24 & 23.675 & 7.17 + mips 70 & 71.42 & 0.778 +