Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay close to their parent bodies. Here we present ALMA 1.3 mm observations of HD 95086, a young early-type star that hosts a directly imaged giant planet b and a massive debris disk with both asteroid- and Kuiper-belt analogs. The location of the Kuiper-belt analog is resolved for the first time. The system can be depicted as a broad ($\Delta R/R \sim$0.84), inclined (30\arcdeg$\pm$3\arcdeg) ring with millimeter emission peaked at 200$\pm$6 au from the star. The 1.3 mm disk emission is consistent with a broad disk with sharp boundaries from 106$\pm$6 to 320$\pm$20 au with a surface density distribution described by a power law with an index of --0.5$\pm$0.2. Our deep ALMA map also reveals a bright source located near the edge of the ring, whose brightness at 1.3 mm and potential spectral energy distribution are consistent with it being a luminous star-forming galaxy at high redshift. We set constraints on the orbital properties of planet b assuming co-planarity with the observed disk.

The architectures of debris disks encode the history of planet formation in these systems. Studies of debris disks via their spectral energy distributions (SEDs) have found infrared excesses arising from cold dust, warm dust, or a combination of the two. The cold outer belts of many systems have been imaged, facilitating their study in great detail. Far less is known about the warm components, including the origin of the dust. The regularity of the disk temperatures indicates an underlying structure that may be linked to the water snow line. If the dust is generated from collisions in an exo-asteroid belt, the dust will likely trace the location of the water snow line in the primordial protoplanetary disk where planetesimal growth was enhanced. If instead the warm dust arises from the inward transport from a reservoir of icy material farther out in the system, the dust location is expected to be set by the current snow line. We analyze the SEDs of a large sample of debris disks with warm components. We find that warm components in single-component systems (those without detectable cold components) follow the primordial snow line rather than the current snow line, so they likely arise from exo-asteroid belts. While the locations of many warm components in two-component systems are also consistent with the primordial snow line, there is more diversity among these systems, suggesting additional effects play a role.

Debris disk morphology is wavelength dependent due to the wide range of particle sizes and size-dependent dynamics influenced by various forces. Resolved images of nearby debris disks reveal complex disk structures that are difficult to distinguish from their spectral energy distributions. Therefore, multi-wavelength resolved images of nearby debris systems provide an essential foundation to understand the intricate interplay between collisional, gravitational, and radiative forces that govern debris disk structures. We present the SOFIA 35 um resolved disk image of epsilon Eri, the closest debris disk around a star similar to the early Sun. Combining with the Spitzer resolved image at 24 um and 15-38 um excess spectrum, we examine two proposed origins of the inner debris in epsilon Eri: (1) in-situ planetesimal belt(s) and (2) dragged-in grains from the cold outer belt. We find that the presence of in-situ dust-producing planetesmial belt(s) is the most likely source of the excess emission in the inner 25 au region. Although a small amount of dragged-in grains from the cold belt could contribute to the excess emission in the inner region, the resolution of the SOFIA data is high enough to rule out the possibility that the entire inner warm excess results from dragged-in grains, but not enough to distinguish one broad inner disk from two narrow belts.

We characterize the first 40 Myr of evolution of circumstellar disks through a unified study of the infrared properties of members of young clusters and associations with ages from 2 Myr up to ~ 40 Myr: NGC 1333, NGC 1960, NGC 2232, NGC 2244, NGC 2362, NGC 2547, IC 348, IC 2395, IC 4665, Chamaeleon I, Orion OB1a and OB1b, Taurus, the \b{eta} Pictoris Moving Group, \r{ho} Ophiuchi, and the associations of Argus, Carina, Columba, Scorpius-Centaurus, and Tucana-Horologium. Our work features: 1.) a filtering technique to flag noisy backgrounds, 2.) a method based on the probability distribution of deflections, P(D), to obtain statistically valid photometry for faint sources, and 3.) use of the evolutionary trend of transitional disks to constrain the overall behavior of bright disks. We find that the fraction of disks three or more times brighter than the stellar photospheres at 24 {\mu}m decays relatively slowly initially and then much more rapidly by ~ 10 Myr. However, there is a continuing component until ~ 35 Myr, probably due primarily to massive clouds of debris generated in giant impacts during the oligarchic/chaotic growth phases of terrestrial planets. If the contribution from primordial disks is excluded, the evolution of the incidence of these oligarchic/chaotic debris disks can be described empirically by a log-normal function with the peak at 12 - 20 Myr, including ~ 13 % of the original population, and with a post-peak mean duration of 10 - 20 Myr.

We investigate whether varying the dust composition (described by the optical constants) can solve a persistent problem in debris disk modeling--the inability to fit the thermal emission without over-predicting the scattered light. We model five images of the beta Pictoris disk: two in scattered light from HST/STIS at 0.58 microns and HST/WFC3 at 1.16 microns, and three in thermal emission from Spitzer/MIPS at 24 microns, Herschel/PACS at 70 microns, and ALMA at 870 microns. The WFC3 and MIPS data are published here for the first time. We focus our modeling on the outer part of this disk, consisting of a parent body ring and a halo of small grains. First, we confirm that a model using astronomical silicates cannot simultaneously fit the thermal and scattered light data. Next, we use a simple, generic function for the optical constants to show that varying the dust composition can improve the fit substantially. Finally, we model the dust as a mixture of the most plausible debris constituents: astronomical silicates, water ice, organic refractory material, and vacuum. We achieve a good fit to all datasets with grains composed predominantly of silicates and organics, while ice and vacuum are, at most, present in small amounts. This composition is similar to one derived from previous work on the HR 4796A disk. Our model also fits the thermal SED, scattered light colors, and high-resolution mid-IR data from T-ReCS for this disk. Additionally, we show that sub-blowout grains are a necessary component of the halo.

Fomalhaut plays an important role in the study of debris disks and small bodies in other planetary systems. The proximity and luminosity of the star make key features of its debris, like the water ice-line, accessible. Here we present ALMA cycle 1, 870 \mu m (345 GHz) observations targeted at the inner part of the Fomalhaut system with a synthesized beam of 0.45"x0.37" (~3 AU linear resolution at the distance of Fomalhaut) and a rms of 26 \mu Jy/beam. The high angular resolution and sensitivity of the ALMA data enable us to place strong constraints on the nature of the warm excess revealed by Spitzer and Herschel observations. We detect a point source at the star position with a total flux consistent with thermal emission from the stellar photosphere. No structures that are brighter than 3\sigma\ are detected in the central 15 AU x 15 AU region. Modeling the spectral energy distribution using parameters expected for a dust-producing planetesimal belt indicates a radial location in the range ~8-15 AU. This is consistent with the location where ice sublimates in Fomalhaut, i.e., an asteroid-belt analog. The 3\sigma\ upper limit for such a belt is <1.3 mJy at 870 \mu m. We also interpret the 2 and 8-13 \mu m interferometric measurements to reveal the structure in the inner 10 AU region as dust naturally connected to this proposed asteroid belt by Poynting-Robertson drag, dust sublimation, and magnetically trapped nano grains.

The final assembly of terrestrial planets occurs via massive collisions, which can launch copious clouds of dust that are warmed by the star and glow in the infrared. We report the real-time detection of a debris-producing impact in the terrestrial planet zone around a 35-million year-old solar analog star. We observed a substantial brightening of the debris disk at 3-5 {\mu}m, followed by a decay over a year, with quasi-periodic modulations of the disk flux. The behavior is consistent with the occurrence of a violent impact that produced vapor out of which a thick cloud of silicate spherules condensed that were ground into dust by collisions. These results demonstrate how the time domain can become a new dimension for the study of terrestrial planet formation.

Luminous debris disks of warm dust in the terrestrial planet zones around solar-like stars are recently found to vary, indicative of ongoing large-scale collisions of rocky objects. We use Spitzer 3.6 and 4.5 {\mu}m time-series observations in 2012 and 2013 (extended to 2014 in one case) to monitor 5 more debris disks with unusually high fractional luminosities ("extreme debris disk"), including P1121 in the open cluster M47 (80 Myr), HD 15407A in the AB Dor moving group (80 Myr), HD 23514 in the Pleiades (120 Myr), HD 145263 in the Upper Sco Association (10 Myr), and the field star BD+20 307 (>1 Gyr). Together with the published results for ID8 in NGC 2547 (35 Myr), this makes the first systematic time-domain investigation of planetary impacts outside the solar system. Significant variations with timescales shorter than a year are detected in five out of the six extreme debris disks we have monitored. However, different systems show diverse sets of characteristics in the time domain, including long-term decay or growth, disk temperature variations, and possible periodicity.

We present our first results from a high-contrast imaging search for planetary mass companions around stars with gapped debris disks, as inferred from the stars' bright infrared excesses. For the six considered stars, we model the disks' unresolved infrared spectral energy distributions (SEDs) in order to derive the temperature and location of the disk components. With VLT/NaCo Apodizing Phase Plate coronagraphic L'-band imaging, we search for planetary mass companions that may be sculpting the disks. We detect neither disks nor companions in this sample, confirmed by comparing plausible point sources with archival data. In order to calculate our mass sensitivity limit, we revisit the stellar age estimates. One target, HD 17848, at 540$\pm$100 Myr old is significantly older than previously estimated. We then discuss our high-contrast imaging results with respect to the disk properties.

HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing 5% of the total dust mass. We also present detailed analysis of the disk SED and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (175 K) belt, a cold (55 K) disk, and an extended disk halo (up to 800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass 5 MJ, could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system.

We report the discovery of a planetary-mass companion, HD 106906 b, with the new Magellan Adaptive Optics (MagAO) + Clio2 system. The companion is detected with Clio2 in three bands: $J$, $K_S$, and $L^\prime$, and lies at a projected separation of 7.1'' (650 AU). It is confirmed to be comoving with its $13\pm2$ Myr-old F5 host using Hubble Space Telescope/Advanced Camera for Surveys astrometry over a time baseline of 8.3 yr. DUSTY and COND evolutionary models predict the companion's luminosity corresponds to a mass of $11\pm2 M_{Jup}$, making it one of the most widely separated planetary-mass companions known. We classify its Magellan/Folded-Port InfraRed Echellette $J/H/K$ spectrum as L$2.5\pm1$; the triangular $H$-band morphology suggests an intermediate surface gravity. HD 106906 A, a pre-main-sequence Lower Centaurus Crux member, was initially targeted because it hosts a massive debris disk detected via infrared excess emission in unresolved Spitzer imaging and spectroscopy. The disk emission is best fit by a single component at 95 K, corresponding to an inner edge of 15-20 AU and an outer edge of up to 120 AU. If the companion is on an eccentric ($e>0.65$) orbit, it could be interacting with the outer edge of the disk. Close-in, planet-like formation followed by scattering to the current location would likely disrupt the disk and is disfavored. Furthermore, we find no additional companions, though we could detect similar-mass objects at projected separations $>35$ AU. In situ formation in a binary-star-like process is more probable, although the companion-to-primary mass ratio, at $<1\%$, is unusually small.

Cold debris disks trace the limits of planet formation or migration in the outer regions of planetary systems, and thus have the potential to answer many of the outstanding questions in wide-orbit planet formation and evolution. We characterized the infrared excess spectral energy distributions of 174 cold debris disks around 546 main-sequence stars observed by both Spitzer IRS and MIPS. We found a trend between the temperature of the inner edges of cold debris disks and the stellar type of the stars they orbit. This argues against the importance of strictly temperature-dependent processes (e.g. non-water ice lines) in setting the dimensions of cold debris disks. Also, we found no evidence that delayed stirring causes the trend. The trend may result from outward planet migration that traces the extent of the primordial protoplanetary disk, or it may result from planet formation that halts at an orbital radius limited by the efficiency of core accretion.

We review the nearby debris disk structures revealed by multi-wavelength images from Spitzer and Herschel, and complemented with detailed spectral energy distribution modeling. Similar to the definition of habitable zones around stars, debris disk structures should be identified and characterized in terms of dust temperatures rather than physical distances so that the heating power of different spectral type of stars is taken into account and common features in disks can be discussed and compared directly. Common features, such as warm (~150 K) dust belts near the water-ice line and cold (~50 K) Kuiper-belt analogs, give rise to our emerging understanding of the levels of order in debris disk structures and illuminate various processes about the formation and evolution of exoplanetary systems. In light of the disk structures in the debris disk twins (Vega and Fomalhaut), and the current limits on the masses of planetary objects, we suggest that the large gap between the warm and cold dust belts is the best signpost for multiple (low-mass) planets beyond the water-ice line.

We present Herschel spatially resolved images of the debris disc orbiting the subgiant Kappa CrB. Not only are these the first resolved images of a debris disc orbiting a subgiant, but Kappa CrB is a rare example of an intermediate mass star where a detailed study of the structure of the planetary system can be made, including both planets and planetesimal belt(s). The only way to discover planets around such stars using the radial velocity technique is to observe 'retired' A stars, which are cooler and slower rotators compared to their main-sequence counterparts. A planetary companion has already been detected orbiting the subgiant Kappa CrB, with revised parameters of m sin i = 2.1MJ and apl = 2.8AU (Johnson et al. 2008a). We present additional Keck I HIRES radial velocity measurements that provide evidence for a second planetary companion, alongside Keck II AO imaging that places an upper limit on the mass of this companion. Modelling of our Herschel images shows that the dust is broadly distributed, but cannot distinguish between a single wide belt (from 20 to 220AU) or two narrow dust belts (at around 40 and 165AU). Given the existence of a second planetary companion beyond approximately 3AU it is possible that the absence of dust within approximately 20AU is caused by dynamical depletion, although the observations are not inconsistent with depletion of these regions by collisional erosion, which occurs at higher rates closer to the star.

We present a 3-5um LBT/MMT adaptive optics imaging study of three Upper Scorpius stars with brown dwarf (BD) companions with very low-masses/mass ratios (M_BD < 25M_Jup; M_BD / M_star ~ 1-2%), and wide separations (300-700 AU): GSC 06214, 1RXS 1609, and HIP 78530. We combine these new thermal IR data with existing 1-4um and 24um photometry to constrain the properties of the BDs and identify evidence for circumprimary/secondary disks in these unusual systems. We confirm that GSC 06214B is surrounded by a disk, further showing this disk produces a broadband IR excess due to small dust near the dust sublimation radius. An unresolved 24um excess in the system may be explained by the contribution from this disk. 1RXS 1609B exhibits no 3-4um excess, nor does its primary; however, the system as a whole has a modest 24um excess, which may come from warm dust around the primary and/or BD. Neither object in the HIP 78530 system exhibits near- to mid-IR excesses. We additionally find that the 1-4um colors of HIP 78530B match a spectral type of M3+-2, inconsistent with the M8 spectral type assigned based on its near-IR spectrum, indicating it may be a low-mass star rather than a BD. We present new upper limits on additional low-mass companions in the system (<5M_Jup beyond 175AU). Finally, we examine the utility of circumsecondary disks as probes of the formation histories of wide BD companions, finding that the presence of a disk may disfavor BD formation near the primary with subsequent outward scattering.

Vega and Fomalhaut, are similar in terms of mass, ages, and global debris disk properties; therefore, they are often referred as "debris disk twins". We present Spitzer 10-35 um spectroscopic data centered at both stars, and identify warm, unresolved excess emission in the close vicinity of Vega for the first time. The properties of the warm excess in Vega are further characterized with ancillary photometry in the mid infrared and resolved images in the far-infrared and submillimeter wavelengths. The Vega warm excess shares many similar properties with the one found around Fomalhaut. The emission shortward of ~30 um from both warm components is well described as a blackbody emission of ~170 K. Interestingly, two other systems, eps Eri and HR 8799, also show such an unresolved warm dust using the same approach. These warm components may be analogous to the solar system's zodiacal dust cloud, but of far greater. The dust temperature and tentative detections in the submillimeter suggest the warm excess arises from dust associated with a planetesimal ring located near the water-frost line and presumably created by processes occurring at similar locations in other debris systems as well. We also review the properties of the 2 um hot excess around Vega and Fomalhaut, showing that the dust responsible for the hot excess is not spatially associated with the dust we detected in the warm belt. We suggest it may arise from hot nano grains trapped in the magnetic field of the star. Finally, the separation between the warm and cold belt is rather large with an orbital ratio >~10 in all four systems. In light of the current upper limits on the masses of planetary objects and the large gap, we discuss the possible implications for their underlying planetary architecture, and suggest that multiple, low-mass planets likely reside between the two belts in Vega and Fomalhaut.

We present observations of the debris disk around gamma Doradus, an F1V star, from the Herschel Key Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/Submillimetre). The disk is well-resolved at 70, 100 and 160 micron, resolved along its major axis at 250 micron, detected but not resolved at 350 micron, and confused with a background source at 500 micron. It is one of our best resolved targets and we find it to have a radially broad dust distribution. The modelling of the resolved images cannot distinguish between two configurations: an arrangement of a warm inner ring at several AU (best-fit 4 AU) and a cool outer belt extending from ~55 to 400 AU or an arrangement of two cool, narrow rings at ~70 AU and ~190 AU. This suggests that any configuration between these two is also possible. Both models have a total fractional luminosity of ~10^{-5} and are consistent with the disk being aligned with the stellar equator. The inner edge of either possible configuration suggests that the most likely region to find planets in this system would be within ~55 AU of the star. A transient event is not needed to explain the warm dust's fractional luminosity.

The majority of debris discs discovered so far have only been detected through infrared excess emission above stellar photospheres. While disc properties can be inferred from unresolved photometry alone under various assumptions for the physical properties of dust grains, there is a degeneracy between disc radius and dust temperature that depends on the grain size distribution and optical properties. By resolving the disc we can measure the actual location of the dust. The launch of Herschel, with an angular resolution superior to previous far-infrared telescopes, allows us to spatially resolve more discs and locate the dust directly. Here we present the nine resolved discs around A stars between 20 and 40 pc observed by the DEBRIS survey. We use these data to investigate the disc radii by fitting narrow ring models to images at 70, 100 and 160 {\mu}m and by fitting blackbodies to full spectral energy distributions. We do this with the aim of finding an improved way of estimating disc radii for unresolved systems. The ratio between the resolved and blackbody radii varies between 1 and 2.5. This ratio is inversely correlated with luminosity and any remaining discrepancies are most likely explained by differences to the minimum size of grain in the size distribution or differences in composition. We find that three of the systems are well fit by a narrow ring, two systems are borderline cases and the other four likely require wider or multiple rings to fully explain the observations, reflecting the diversity of planetary systems.

We have obtained Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) 24 {\mu}m and 70 {\mu}m observations of 215 nearby, Hipparcos B- and A-type common proper motion single and binary systems in the nearest OB association, Scorpius-Centaurus. Combining our MIPS observations with those of other ScoCen stars in the literature, we estimate 24 {\mu}m B+A-type disk fractions of 17/67 (25+6%), 36/131 (27+4%), and 23/95 (24+5%) for Upper Scorpius (\sim11 Myr), Upper Centaurus Lupus (\sim15 Myr), and Lower Centaurus Crux (\sim17 Myr), respectively, somewhat smaller disk fractions than previously obtained for F- and G-type members. We confirm previous IRAS excess detections and present new discoveries of 51 protoplanetary and debris disk systems, with fractional infrared luminosities ranging from LIR/L\ast = 1e-6 to 1e-2 and grain temperatures ranging from Tgr = 40 - 300 K. In addition, we confirm that the 24 {\mu}m and 70 {\mu}m excesses (or fractional infrared luminosities) around B+A type stars are smaller than those measured toward F+G type stars and hypothesize that the observed disk property dependence on stellar mass may be the result of a higher stellar companion fraction around B- and A-type stars at 10 - 200 AU and/or the presence of Jupiter-mass companions in the disks around F- and G- type stars. Finally, we note that the majority of the ScoCen 24 {\mu}m excess sources also possess 12 {\mu}m excess, indicating that Earth-like planets may be forming via collisions in the terrestrial planet zone at \sim10 - 100 Myr.

We present diffraction-limited \ks band and \lprime adaptive optics images of the edge-on debris disk around the nearby F2 star HD 15115, obtained with a single 8.4 m primary mirror at the Large Binocular Telescope. At \ks band the disk is detected at signal-to-noise per resolution element (SNRE) \about 3-8 from \about 1-2\fasec 5 (45-113 AU) on the western side, and from \about 1.2-2\fasec 1 (63-90 AU) on the east. At \lprime the disk is detected at SNRE \about 2.5 from \about 1-1\fasec 45 (45-90 AU) on both sides, implying more symmetric disk structure at 3.8 \microns . At both wavelengths the disk has a bow-like shape and is offset from the star to the north by a few AU. A surface brightness asymmetry exists between the two sides of the disk at \ks band, but not at \lprime . The surface brightness at \ks band declines inside 1\asec (\about 45 AU), which may be indicative of a gap in the disk near 1\asec. The \ks - \lprime disk color, after removal of the stellar color, is mostly grey for both sides of the disk. This suggests that scattered light is coming from large dust grains, with 3-10 \microns -sized grains on the east side and 1-10 \microns dust grains on the west. This may suggest that the west side is composed of smaller dust grains than the east side, which would support the interpretation that the disk is being dynamically affected by interactions with the local interstellar medium.

Debris disks with extremely large infrared excesses (fractional luminosities $> 10^{-2}$) are rare. Those with ages between 30 and 130 Myr are of interest because their evolution has progressed well beyond that of protoplanetary disks (which dissipate with a timescale of order 3 Myr), yet they represent a period when dynamical models suggest that terrestrial planet building may still be progressing through large, violent collisions that could yield large amounts of debris and large infrared excesses. For example, our Moon was formed through a violent collision of two large proto-planets during this age range. We report two disks around the solar-like stars ID8 and HD 23514 in this age range where the 24 {\mu}m infrared excesses vary on timescales of a few years, even though the stars are not variable in the optical. Variations this rapid are difficult to understand if the debris is produced by collisional cascades, as it is for most debris disks. It is possible that the debris in these two systems arises in part from condensates from silicate-rich vapor produced in a series of violent collisions among relatively large bodies. If their evolution is rapid, the rate of detection of extreme excesses would indicate that major collisions may be relatively common in this age range.

We describe a search for dust created in collisions between the Saturnian irregular satellites using archival \emph{Spitzer} MIPS observations. Although we detected a degree scale Saturn-centric excess that might be attributed to an irregular satellite dust cloud, we attribute it to the far-field wings of the PSF due to nearby Saturn. The Spitzer PSF is poorly characterised at such radial distances, and we expect PSF characterisation to be the main issue for future observations that aim to detect such dust. The observations place an upper limit on the level of dust in the outer reaches of the Saturnian system, and constrain how the size distribution extrapolates from the smallest known (few km) size irregulars down to micron-size dust. Because the size distribution is indicative of the strength properties of irregulars, we show how our derived upper limit implies irregular satellite strengths more akin to comets than asteroids. This conclusion is consistent with their presumed capture from the outer regions of the Solar System.

We consider the limiting case for orbital stability of the companions to HR 8799. This case is only consistent with ages for the system of ~100 Myr, not with the 1 Gyr age proposed from astroseismology. The discrepancy probably arises because the inclination of the star is smaller than assumed in analyzing the astroseismology data. Given this young age, the best estimates of the companion masses place them by a small margin on the planet side of the division between planets and brown dwarfs.

In addition to the Sun, six other stars are known to harbor multiple planets and debris disks: HD 69830, HD 38529, HD 128311, HD 202206, HD 82943 and HR 8799. In this paper we set constraints on the location of the dust-producing planetesimals around the latter four systems. We use a radiative transfer model to analyze the spectral energy distributions of the dust disks (including two new Spitzer IRS spectra presented in this paper), and a dynamical model to assess the long-term stability of the planetesimals' orbits. As members of a small group of stars that show evidence of harboring a multiple planets and planetesimals, their study can help us learn about the diversity of planetary systems.

We have carried out a spatial-kinematical study of three proto-planetary nebulae, IRAS 16594-4656, Hen 3-401, and Rob 22. High-resolution H2 images were obtained with NICMOS on the HST and high-resolution spectra were obtained with the Phoenix spectrograph on Gemini-South. IRAS 16594-4656 shows a "peanut-shaped" bipolar structure with H2 emission from the walls and from two pairs of more distant, point-symmetric faint blobs. The velocity structure shows the polar axis to be in the plane of the sky, contrary to the impression given by the more complex visual image and the visibility of the central star, with an ellipsoidal velocity structure. Hen 3-401 shows the H2 emission coming from the walls of the very elongated, open-ended lobes seen in visible light, along with a possible small disk around the star. The bipolar lobes appear to be tilted 10-15 deg with respect to the plane of the sky and their kinematics display a Hubble-like flow. In Rob 22, the H2 appears in the form of an "S" shape, approximately tracing out the similar pattern seen in the visible. H2 is especially seen at the ends of the lobes and at two opposite regions close to the unseen central star. The axis of the lobes is nearly in the plane of the sky. Expansion ages of the lobes are calculated to be approximately 1600 yr (IRAS 16594-4656), 1100 yr (Hen 3-401), and 640 yr (Rob 22), based upon approximate distances.