Source author record

Christopher C. Stark

Christopher C. Stark appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

ALMA Images the Eccentric HD 53143 Debris Disk

We present ALMA 1.3 mm observations of the HD~53143 debris disk - the first infrared or millimeter image produced of this ~1 Gyr-old solar-analogue. Previous HST STIS coronagraphic imaging did not detect flux along the minor axis of the disk which could suggest a face-on geometry with two 'clumps' of dust. These ALMA observations reveal a disk with a strikingly different structure. In order to fit models to the millimeter visibilities and constrain the uncertainties on the disk parameters, we adopt an MCMC approach. This is the most eccentric debris disk observed to date with a forced eccentricity of $0.21\pm0.02$, nearly twice that of the Fomalhaut debris disk, and also displays apocenter glow. Although this eccentric model fits the outer debris disk well, there are significant interior residuals remaining that may suggest a possible edge-on inner disk, which remains unresolved in these observations. Combined with the observed structure difference between HST and ALMA, these results suggest a potential previous scattering event or dynamical instability in this system. We also note that the stellar flux changes considerably over the course of our observations, suggesting flaring at millimeter wavelengths. Using simultaneous TESS observations, we determine the stellar rotation period to be $9.6\pm0.1$ days.

preprint2022arXiv

ExoVista: A Suite of Planetary System Models for Exoplanet Studies

Studies of future space- and ground-based exoplanet surveys often rely on models of planetary systems to simulate instrument response, estimate scientific yields, perform trade analyses, and study efficient observation strategies. Until now, no planetary system models contained all of the basic physics necessary to enable study with all of the major exoplanet detection methods. Here we introduce a suite of such models generated by a new tool, exoVista. The exoVista tool quickly generates thousands of models of quasi-self-consistent planetary systems around known nearby stars at scattered light wavelengths and efficiently records the position, velocity, spectrum, and physical parameters of all bodies as functions of time. The modeled planetary systems can be used to simulate surveys using the direct imaging, transit, astrometric, and radial velocity techniques, as well as the overlap of these different methods.

preprint2020arXiv

Toward Complete Characterization: Prospects for Directly Imaging Transiting Exoplanets

High contrast direct imaging of exoplanets can provide many important observables, including measurements of the orbit, spectra that probe the lower layers of the atmosphere, and phase variations of the planet, but cannot directly measure planet radius or mass. Our future understanding of directly imaged exoplanets will therefore rely on extrapolated models of planetary atmospheres and bulk composition, which need robust calibration. We estimate the population of extrasolar planets that could serve as calibrators for these models. Critically, this population of "standard planets" must be accessible to both direct imaging and the transit method, allowing for radius measurement. We show that the search volume of a direct imaging mission eventually overcomes the transit probability falloff with semi-major axis, so that as long as cold planets are not exceedingly rare, the population of transiting planets and directly imageable planets overlaps. Using current extrapolations of Kepler occurrence rates, we estimate that ~8 standard planets could be characterized shortward of 800 nm with an ambitious future direct imaging mission like LUVOIR-A and several dozen could be detected at V band. We show the design space that would expand the sample size and discuss the extent to which ground- and space-based surveys could detect this small but crucial population of planets.

preprint2016arXiv

Apodized pupil Lyot coronagraphs for arbitrary apertures. V. Hybrid Shaped Pupil designs for imaging Earth-like planets with future space observatories

We introduce a new class of solutions for Apodized Pupil Lyot Coronagraphs (APLC) with segmented aperture telescopes to remove broadband diffracted light from a star with a contrast level of $10^{10}$. These new coronagraphs provide a key advance to enabling direct imaging and spectroscopy of Earth twins with future large space missions. Building on shaped pupil (SP) apodization optimizations, our approach enables two-dimensional optimizations of the system to address any aperture features such as central obstruction, support structures or segment gaps. We illustrate the technique with a design that could reach $10^{10}$ contrast level at 34\,mas for a 12\,m segmented telescope over a 10\% bandpass centered at a wavelength $λ_0=$500\,nm. These designs can be optimized specifically for the presence of a resolved star, and in our example, for stellar angular size up to 1.1\,mas. This would allow probing the vicinity of Sun-like stars located beyond 4.4\,pc, therefore fully retiring this concern. If the fraction of stars with Earth-like planets is $η_{\Earth}=0.1$, with 18\% throughput, assuming a perfect, stable wavefront and considering photon noise only, 12.5 exo-Earth candidates could be detected around nearby stars with this design and a 12\,m space telescope during a five-year mission with two years dedicated to exo-Earth detection (one total year of exposure time and another year of overheads). Our new hybrid APLC/SP solutions represent the first numerical solution of a coronagraph based on existing mask technologies and compatible with segmented apertures, and that can provide contrast compatible with detecting and studying Earth-like planets around nearby stars. They represent an important step forward towards enabling these science goals with future large space missions.

preprint2016arXiv

Deep HST/STIS Visible-Light Imaging of Debris Systems around Solar Analog Hosts

We present new Hubble Space Telescope observations of three a priori known starlight-scattering circumstellar debris systems (CDSs) viewed at intermediate inclinations around nearby close-solar analog stars: HD 207129, HD 202628, and HD 202917. Each of these CDSs possesses ring-like components that are more-massive analogs of our solar system's Edgeworth- Kuiper belt. These systems were chosen for follow-up observations to provide higher-fidelity and better sensitivity imaging for the sparse sample of solar-analog CDSs that range over two decades in systemic ages with HD 202628 and HD 202917 (both ~ 2.3 Gyr) currently the oldest CDSs imaged in visible or near-IR light. These deep (10 - 14 ksec) observations, with six-roll point-spread-function template subtracted visible-light coronagraphy using the Space Telescope Imaging Spectrograph, were designed to better reveal their angularly large, diffuse/low surface brightness, debris rings, and for all targets probe their exo-ring environments for starlight-scattering materials that present observational challenges for current ground-based facilities and instruments. Contemporaneously also observing with a narrower occulter position, these observations additionally probe the CDS endo-ring environments seen to be relatively devoid of scatterers. We discuss the morphological, geometrical, and photometric properties of these CDSs also in the context of other FGK-star hosted CDSs we have previously imaged as a homogeneously observed ensemble. From this combined sample we report a general decay in quiescent disk F_disk/F_star optical brightness ~ t^-0.8, similar to what is seen in at thermal IR wavelengths, and CDSs with a significant diversity in scattering phase asymmetries, and spatial distributions of their starlight-scattering grains.

preprint2016arXiv

Discovery of an Inner Disk Component around HD 141569 A

We report the discovery of a scattering component around the HD 141569 A circumstellar debris system, interior to the previously known inner ring. The discovered inner disk component, obtained in broadband optical light with HST/STIS coronagraphy, was imaged with an inner working angle of 0".25, and can be traced from 0".4 (~46 AU) to 1".0 (~116 AU) after deprojection using i=55deg. The inner disk component is seen to forward scatter in a manner similar to the previously known rings, has a pericenter offset of ~6 AU, and break points where the slope of the surface brightness changes. It also has a spiral arm trailing in the same sense as other spiral arms and arcs seen at larger stellocentric distances. The inner disk spatially overlaps with the previously reported warm gas disk seen in thermal emission. We detect no point sources within 2" (~232 AU), in particular in the gap between the inner disk component and the inner ring. Our upper limit of 9+/-3 M_J is augmented by a new dynamical limit on single planetary mass bodies in the gap between the inner disk component and the inner ring of 1 M_J, which is broadly consistent with previous estimates.

preprint2016arXiv

Maximized ExoEarth Candidate Yields for Starshades

The design and scale of a future mission to directly image and characterize potentially Earth-like planets will be impacted, to some degree, by the expected yield of such planets. Recent efforts to increase the estimated yields, by creating observation plans optimized for the detection and characterization of Earth-twins, have focused solely on coronagraphic instruments; starshade-based missions could benefit from a similar analysis. Here we explore how to prioritize observations for a starshade given the limiting resources of both fuel and time, present analytic expressions to estimate fuel use, and provide efficient numerical techniques for maximizing the yield of starshades. We implemented these techniques to create an approximate design reference mission code for starshades and used this code to investigate how exoEarth candidate yield responds to changes in mission, instrument, and astrophysical parameters for missions with a single starshade. We find that a starshade mission operates most efficiently somewhere between the fuel- and exposure-time limited regimes, and as a result, is less sensitive to photometric noise sources as well as parameters controlling the photon collection rate in comparison to a coronagraph. We produced optimistic yield curves for starshades, assuming our optimized observation plans are schedulable and future starshades are not thrust-limited. Given these yield curves, detecting and characterizing several dozen exoEarth candidates requires either multiple starshades or an eta_Earth > ~0.3.

preprint2015arXiv

ATLAST detector needs for direct spectroscopic biosignature characterization in the visible and near-IR

Are we alone? Answering this ageless question will be a major focus for astrophysics in coming decades. Our tools will include unprecedentedly large UV-Optical-IR space telescopes working with advanced coronagraphs and starshades. Yet, these facilities will not live up to their full potential without better detectors than we have today. To inform detector development, this paper provides an overview of visible and near-IR (VISIR; $λ=0.4-1.8~μ\textrm{m}$) detector needs for the Advanced Technology Large Aperture Space Telescope (ATLAST), specifically for spectroscopic characterization of atmospheric biosignature gasses. We also provide a brief status update on some promising detector technologies for meeting these needs in the context of a passively cooled ATLAST.

preprint2015arXiv

First images of debris disks around TWA 7, TWA 25, HD 35650, and HD 377

We present the first images of four debris disks observed in scattered light around the young (4--250 Myr old) M dwarfs TWA 7 and TWA 25, the K6 star HD 35650, and the G2 star HD 377. We obtained these images by reprocessing archival Hubble Space Telescope NICMOS coronagraph data with modern post-processing techniques as part of the Archival Legacy Investigation of Circumstellar Environments (ALICE) program. All four disks appear faint and compact compared with other debris disks resolved in scattered light. The disks around TWA 25, HD 35650, and HD 377 appear very inclined, while TWA 7's disk is viewed nearly face-on. The surface brightness of HD 35650's disk is strongly asymmetric. These new detections raise the number of disks resolved in scattered light around M and late-K stars from one (the AU Mic system) to four. This new sample of resolved disks enables comparative studies of heretofore scarce debris disks around low-mass stars relative to solar-type stars.

preprint2015arXiv

Lower Limits on Aperture Size for an ExoEarth-Detecting Coronagraphic Mission

The yield of Earth-like planets will likely be a primary science metric for future space-based missions that will drive telescope aperture size. Maximizing the exoEarth candidate yield is therefore critical to minimizing the required aperture. Here we describe a method for exoEarth candidate yield maximization that simultaneously optimizes, for the first time, the targets chosen for observation, the number of visits to each target, the delay time between visits, and the exposure time of every observation. This code calculates both the detection time and multi-wavelength spectral characterization time required for planets. We also refine the astrophysical assumptions used as inputs to these calculations, relying on published estimates of planetary occurrence rates as well as theoretical and observational constraints on terrestrial planet sizes and classical habitable zones. Given these astrophysical assumptions, optimistic telescope and instrument assumptions, and our new completeness code that produces the highest yields to date, we suggest lower limits on the aperture size required to detect and characterize a statistically-motivated sample of exoEarths.

preprint2015arXiv

The Pseudo-zodi Problem for Edge-on Planetary Systems

Future direct observations of extrasolar Earth-sized planets in the habitable zone could be hampered by a worrisome source of noise, starlight-reflecting exozodiacal dust. Mid-infrared surveys are currently underway to constrain the amount of exozodiacal dust in the habitable zones around nearby stars. However, at visible wavelengths another source of dust, invisible to these surveys, may dominate over exozodiacal dust. For systems observed near edge-on, a cloud of dust with face-on optical depth 10^-7 beyond ~5 AU can mimic the surface brightness of a cloud of exozodiacal dust with equal optical depth if the dust grains are sufficiently forward-scattering. We posit that dust migrating inward from cold debris belts via Poynting-Robertson drag could produce this "pseudo-zodiacal" effect, potentially making it ~50% as common as exozodiacal clouds. We place constraints on the disk radii and scattering phase function required to produce the effect.

preprint2014arXiv

Maximizing the ExoEarth Candidate Yield from a Future Direct Imaging Mission

ExoEarth yield is a critical science metric for future exoplanet imaging missions. Here we estimate exoEarth candidate yield using single visit completeness for a variety of mission design and astrophysical parameters. We review the methods used in previous yield calculations and show that the method choice can significantly impact yield estimates as well as how the yield responds to mission parameters. We introduce a method, called Altruistic Yield Optimization, that optimizes the target list and exposure times to maximize mission yield, adapts maximally to changes in mission parameters, and increases exoEarth candidate yield by up to 100% compared to previous methods. We use Altruistic Yield Optimization to estimate exoEarth candidate yield for a large suite of mission and astrophysical parameters using single visit completeness. We find that exoEarth candidate yield is most sensitive to telescope diameter, followed by coronagraph inner working angle, followed by coronagraph contrast, and finally coronagraph contrast noise floor. We find a surprisingly weak dependence of exoEarth candidate yield on exozodi level. Additionally, we provide a quantitative approach to defining a yield goal for future exoEarth-imaging missions.

preprint2014arXiv

On the Morphology and Chemical Composition of the HR 4796A Debris Disk

[abridged] We present resolved images of the HR 4796A debris disk using the Magellan adaptive optics system paired with Clio-2 and VisAO. We detect the disk at 0.77 \microns, 0.91 \microns, 0.99 \microns, 2.15 \microns, 3.1 \microns, 3.3 \microns, and 3.8 \microns. We find that the deprojected center of the ring is offset from the star by 4.76$\pm$1.6 AU and that the deprojected eccentricity is 0.06$\pm$0.02, in general agreement with previous studies. We find that the average width of the ring is 14$^{+3}_{-2}%$, also comparable to previous measurements. Such a narrow ring precludes the existence of shepherding planets more massive than \about 4 \mj, comparable to hot-start planets we could have detected beyond \about 60 AU in projected separation. Combining our new scattered light data with archival HST/STIS and HST/NICMOS data at \about 0.5-2 \microns, along with previously unpublished Spitzer/MIPS thermal emission data and all other literature thermal data, we set out to constrain the chemical composition of the dust grains. After testing 19 individual root compositions and more than 8,400 unique mixtures of these compositions, we find that good fits to the scattered light alone and thermal emission alone are discrepant, suggesting that caution should be exercised if fitting to only one or the other. When we fit to both the scattered light and thermal emission simultaneously, we find mediocre fits (reduced chi-square \about 2). In general, however, we find that silicates and organics are the most favored, and that water ice is usually not favored. These results suggest that the common constituents of both interstellar dust and solar system comets also may reside around HR 4796A, though improved modeling is necessary to place better constraints on the exact chemical composition of the dust.

preprint2014arXiv

Probing for Exoplanets Hiding in Dusty Debris Disks: Disk Imaging, Characterization, and Exploration with HST/STIS Multi-Roll Coronagraphy

Spatially resolved scattered-light images of circumstellar (CS) debris in exoplanetary systems constrain the physical properties and orbits of the dust particles in these systems. They also inform on co-orbiting (but unseen) planets, systemic architectures, and forces perturbing starlight-scattering CS material. Using HST/STIS optical coronagraphy, we have completed the observational phase of a program to study the spatial distribution of dust in ten CS debris systems, and one "mature" protoplanetrary disk all with HST pedigree, using PSF-subtracted multi-roll coronagraphy. These observations probe stellocentric distances > 5 AU for the nearest stars, and simultaneously resolve disk substructures well beyond, corresponding to the giant planet and Kuiper belt regions in our Solar System. They also disclose diffuse very low-surface brightness dust at larger stellocentric distances. We present new results inclusive of fainter disks such as HD92945 confirming, and better revealing, the existence of a narrow inner debris ring within a larger diffuse dust disk. Other disks with ring-like sub-structures, significant asymmetries and complex morphologies include: HD181327 with a posited spray of ejecta from a recent massive collision in an exo-Kuiper belt; HD61005 suggested interacting with the local ISM; HD15115 & HD32297, discussed also in the context of environmental interactions. These disks, and HD15745, suggest debris system evolution cannot be treated in isolation. For AU Mic's edge-on disk, out-of-plane surface brightness asymmetries at > 5 AU may implicate one or more planetary perturbers. Time resolved images of the MP Mus proto-planetary disk provide spatially resolved temporal variability in the disk illumination. These and other new images from our program enable direct inter-comparison of the architectures of these exoplanetary debris systems in the context of our own Solar System.

preprint2014arXiv

Revealing Asymmetries in the HD 181327 Debris Disk: A Recent Massive Collision or ISM Warping

New multi-roll coronagraphic images of the HD 181327 debris disk obtained using the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST) reveal the debris ring in its entirety at high S/N and unprecedented spatial resolution. We present and apply a new multi-roll image processing routine to identify and further remove quasi-static PSF-subtraction residuals and quantify systematic uncertainties. We also use a new iterative image deprojection technique to constrain the true disk geometry and aggressively remove any surface brightness asymmetries that can be explained without invoking dust density enhancements/deficits. The measured empirical scattering phase function for the disk is more forward scattering than previously thought and is not well-fit by a Henyey-Greenstein function. The empirical scattering phase function varies with stellocentric distance, consistent with the expected radiation pressured-induced size segregation exterior to the belt. Within the belt, the empirical scattering phase function contradicts unperturbed debris ring models, suggesting the presence of an unseen planet. The radial profile of the flux density is degenerate with a radially-varying scattering phase function; therefore estimates of the ring's true width and edge slope may be highly uncertain. We detect large scale asymmetries in the disk, consistent with either the recent catastrophic disruption of a body with mass >1% the mass of Pluto, or disk warping due to strong interactions with the interstellar medium (ISM).

preprint2013arXiv

A Search for Exozodiacal Clouds with Kepler

Planets embedded within dust disks may drive the formation of large scale clumpy dust structures by trapping dust into resonant orbits. Detection and subsequent modeling of the dust structures would help constrain the mass and orbit of the planet and the disk architecture, give clues to the history of the planetary system, and provide a statistical estimate of disk asymmetry for future exoEarth-imaging missions. Here we present the first search for these resonant structures in the inner regions of planetary systems by analyzing the light curves of hot Jupiter planetary candidates identified by the Kepler mission. We detect only one candidate disk structure associated with KOI 838.01 at the 3-sigma confidence level, but subsequent radial velocity measurements reveal that KOI 838.01 is a grazing eclipsing binary and the candidate disk structure is a false positive. Using our null result, we place an upper limit on the frequency of dense exozodi structures created by hot Jupiters. We find that at the 90% confidence level, less than 21% of Kepler hot Jupiters create resonant dust clumps that lead and trail the planet by ~90 degrees with optical depths >~5*10^-6, which corresponds to the resonant structure expected for a lone hot Jupiter perturbing a dynamically cold dust disk 50 times as dense as the zodiacal cloud.

preprint2013arXiv

A Survey for Very Short-Period Planets in the Kepler Data

We conducted a search for very short-period transiting objects in the publicly available Kepler dataset. Our preliminary survey has revealed four planetary candidates, all with orbital periods less than twelve hours. We have analyzed the data for these candidates using photometric models that include transit light curves, ellipsoidal variations, and secondary eclipses to constrain the candidates' radii, masses, and effective temperatures. Even with masses of only a few Earth masses, the candidates' short periods mean they may induce stellar radial velocity signals (a few m/s) detectable by currently operating facilities. The origins of such short-period planets are unclear, but we discuss the possibility that they may be the remnants of disrupted hot Jupiters. Whatever their origins, if confirmed as planets, these candidates would be among the shortest-period planets ever discovered. Such planets would be particularly amenable to discovery by the planned TESS mission.

preprint2011arXiv

The Transit Light Curve of an Exozodiacal Dust Cloud

Planets embedded within debris disks gravitationally perturb nearby dust and can create clumpy, azimuthally asymmetric circumstellar ring structures that rotate in lock with the planet. The Earth creates one such structure in the solar zodiacal dust cloud. In an edge-on system, the dust "clumps" periodically pass in front of the star as the planet orbits, occulting and forward-scattering starlight. In this paper, we predict the shape and magnitude of the corresponding transit signal. To do so, we model the dust distributions of collisional, steady-state exozodiacal clouds perturbed by planetary companions. We examine disks with dusty ring structures formed by the planet's resonant trapping of in-spiraling dust for a range of planet masses and semi-major axes, dust properties, and disk masses. We synthesize edge-on images of these models and calculate the transit signatures of the resonant ring structures. The transit light curves created by dusty resonant ring structures typically exhibit two broad transit minima that lead and trail the planetary transit. We find that Jupiter-mass planets embedded within disks hundreds of times denser than our zodiacal cloud can create resonant ring structures with transit depths up to $\sim10^{-4}$, possibly detectable with \emph{Kepler}. Resonant rings produced by planets more or less massive than Jupiter produce smaller transit depths. Observations of these transit signals may provide upper limits on the degree of asymmetry in exozodiacal clouds.

preprint2010arXiv

Collisional Grooming Models of the Kuiper Belt Dust Cloud

We modeled the 3-D structure of the Kuiper Belt dust cloud at four different dust production rates, incorporating both planet-dust interactions and grain-grain collisions using the collisional grooming algorithm. Simulated images of a model with a face-on optical depth of ~10^-4 primarily show an azimuthally-symmetric ring at 40-47 AU in submillimeter and infrared wavelengths; this ring is associated with the cold classical Kuiper Belt. For models with lower optical depths (10^-6 and 10^-7), synthetic infrared images show that the ring widens and a gap opens in the ring at the location of of Neptune; this feature is caused by trapping of dust grains in Neptune's mean motion resonances. At low optical depths, a secondary ring also appears associated with the hole cleared in the center of the disk by Saturn. Our simulations, which incorporate 25 different grain sizes, illustrate that grain-grain collisions are important in sculpting today's Kuiper Belt dust, and probably other aspects of the Solar System dust complex; collisions erase all signs of azimuthal asymmetry from the submillimeter image of the disk at every dust level we considered. The model images switch from being dominated by resonantly-trapped small grains ("transport dominated") to being dominated by the birth ring ("collision dominated") when the optical depth reaches a critical value of tau ~ v/c, where v is the local Keplerian speed.

preprint2010arXiv

The cold origin of the warm dust around epsilon Eridani

Context: The K2V star eps Eri hosts one known inner planet, an outer Kuiper belt analog, and an inner disk of warm dust. Spitzer/IRS measurements indicate that the warm dust is present at distances as close as a few AU from the star. Its origin is puzzling, since an "asteroid belt" that could produce this dust would be unstable because of the known inner planet. Aims: Here we test the hypothesis that the observed warm dust is generated by collisions in the outer belt and is transported inward by Poynting-Robertson (P-R) drag and strong stellar winds. Methods: We simulated a steady-state distribution of dust particles outside 10AU with a collisional code and in the inner region (r<10AU) with single-particle numerical integrations. By assuming homogeneous spherical dust grains composed of water ice and silicate, we calculated the thermal emission of the dust and compared it with observations. We investigated two different orbital configurations for the inner planet inferred from RV measurements, one with a highly eccentric orbit of e=0.7 and another one with a moderate one of e=0.25. We also produced a simulation without a planet. Results: Our models can reproduce the shape and magnitude of the observed SED from mid-IR to sub-mm wavelengths, as well as the Spitzer/MIPS radial brightness profiles. The best-fit dust composition includes both ice and silicates. The results are similar for the two possible planetary orbits and without a planet. Conclusions: The observed warm dust in the system can indeed stem from the outer belt and be transported inward by P-R and stellar wind drag. The inner planet has little effect on the distribution of dust, so that the planetary orbit could not be constrained. Reasonable agreement between the model and observations can only be achieved by relaxing the assumption of purely silicate dust and assuming a mixture of silicate and ice in comparable amounts.

preprint2009arXiv

51 Ophiuchus: A Possible Beta Pictoris Analog Measured with the Keck Interferometer Nuller

We present observations of the 51 Ophiuchi circumstellar disk made with the Keck interferometer operating in nulling mode at N-band. We model these data simultaneously with VLTI-MIDI visibility data and a Spitzer IRS spectrum using a variety of optically-thin dust cloud models and an edge-on optically-thick disk model. We find that single-component optically-thin disk models and optically-thick disk models are inadequate to reproduce the observations, but an optically-thin two-component disk model can reproduce all of the major spectral and interferometric features. Our preferred disk model consists of an inner disk of blackbody grains extending to ~4 AU and an outer disk of small silicate grains extending out to ~1200 AU. Our model is consistent with an inner "birth" disk of continually colliding parent bodies producing an extended envelope of ejected small grains. This picture resembles the disks around Vega, AU Microscopii, and Beta Pictoris, supporting the idea that 51 Ophiuchius may be a Beta Pictoris analog.

preprint2009arXiv

A New Algorithm for Self-Consistent 3-D Modeling of Collisions in Dusty Debris Disks

We present a new "collisional grooming" algorithm that enables us to model images of debris disks where the collision time is less than the Poynting Robertson time for the dominant grain size. Our algorithm uses the output of a collisionless disk simulation to iteratively solve the mass flux equation for the density distribution of a collisional disk containing planets in 3 dimensions. The algorithm can be run on a single processor in ~1 hour. Our preliminary models of disks with resonant ring structures caused by terrestrial mass planets show that the collision rate for background particles in a ring structure is enhanced by a factor of a few compared to the rest of the disk, and that dust grains in or near resonance have even higher collision rates. We show how collisions can alter the morphology of a resonant ring structure by reducing the sharpness of a resonant ring's inner edge and by smearing out azimuthal structure. We implement a simple prescription for particle fragmentation and show how Poynting-Robertson drag and fragmentation sort particles by size, producing smaller dust grains at smaller circumstellar distances. This mechanism could cause a disk to look different at different wavelengths, and may explain the warm component of dust interior to Fomalhaut's outer dust ring seen in the resolved 24 micron Spitzer image of this system.

preprint2009arXiv

Collisional Grooming of Debris Disks

Debris disk images show clumps, rings, warps, and other structures, many of which have been interpreted as perturbations from hidden planets. But so far, no models of these structures have properly accounted for collisions between dust grains. We have developed new steady-state 3D models of debris disks that self-consistently incorporate grain-grain collisions. We summarize our algorithm and use it to illustrate how collisions interact with resonant trapping in the presence of a planet.

preprint2008arXiv

The Detectability of Exo-Earths and Super-Earths Via Resonant Signatures in Exozodiacal Clouds

Directly imaging extrasolar terrestrial planets necessarily means contending with the astrophysical noise of exozodiacal dust and the resonant structures created by these planets in exozodiacal clouds. Using a custom tailored hybrid symplectic integrator we have constructed 120 models of resonant structures created by exo-Earths and super-Earths on circular orbits interacting with collisionless steady-state dust clouds around a Sun-like star. Our models include enough particles to overcome the limitations of previous simulations that were often dominated by a handful of long-lived particles, allowing us to quantitatively study the contrast of the resulting ring structures. We found that in the case of a planet on a circular orbit, for a given star and dust source distribution, the morphology and contrast of the resonant structures depend on only two parameters: planet mass and $\sqrt{a_{\rm p}}/β$, where $a_{\rm p}$ is the planet's semi-major axis and $β$ is the ratio of radiation pressure force to gravitational force on a grain. We constructed multiple-grain-size models of 25,000 particles each and showed that in a collisionless cloud, a Dohnanyi crushing law yields a resonant ring whose optical depth is dominated by the largest grains in the distribution, not the smallest. We used these models to estimate the mass of the lowest-mass planet that can be detected through observations of a resonant ring for a variety of assumptions about the dust cloud and the planet's orbit. Our simulations suggest that planets with mass as small as a few times Mar's mass may produce detectable signatures in debris disks for semi-major axes greater than 10 AU.