Discovery of extreme asymmetry in the debris disk surrounding HD 15115

We report the first scattered light detection of a dusty debris disk surrounding the F2V star HD 15115 using the Hubble Space Telescope in the optical, and Keck adaptive optics in the near-infrared. The most remarkable property of the HD 15115 disk relative to other debris disks is its extreme length asymmetry. The east side of the disk is detected to ~315 AU radius, whereas the west side of the disk has radius >550 AU. We find a blue optical to near-infrared scattered light color relative to the star that indicates grain scattering properties similar to the AU Mic debris disk. The existence of a large debris disk surrounding HD 15115 adds further evidence for membership in the Beta Pic moving group, which was previously argued based on kinematics alone. Here we hypothesize that the extreme disk asymmetry is due to dynamical perturbations from HIP 12545, an M star 0.5 degrees (0.38 pc) east of HD 15115 that shares a common proper motion vector, heliocentric distance, galactic space velocity, and age.Comment: 4 pages, 3 figures, ApJ Letters, accepte

Feasibility of transit photometry of nearby debris discs

Dust in debris discs is constantly replenished by collisions between larger objects. In this paper, we investigate a method to detect these collisions. We generate models based on recent results on the Fomalhaut debris disc, where we simulate a background star transiting behind the disc, due to the proper motion of Fomalhaut. By simulating the expanding dust clouds caused by the collisions in the debris disc, we investigate whether it is possible to observe changes in the brightness of the background star. We conclude that in the case of the Fomalhaut debris disc, changes in the optical depth can be observed, with values of the optical depth ranging from $10^{-0.5}$ for the densest dust clouds to $10^{-8}$ for the most diffuse clouds with respect to the background optical depth of $\sim1.2\times10^{-3}$.Comment: 19 pages, 15 figures, accepted for publication in MNRA

STIS Coronagraphic Imaging of Fomalhaut: Main Belt Structure and the Orbit of Fomalhaut b

We present new optical coronagraphic data of the bright star Fomalhaut obtained with the HST in 2010/2012 using STIS. Fomalhaut b is recovered at both epochs to high significance. The observations include the discoveries of tenuous nebulosity beyond the main dust belt detected to at least 209 AU projected radius and a ~50 AU wide azimuthal gap in the belt northward of Fom b. The morphology of Fomalhaut b appears elliptical in the STIS detections. We show that residual noise in the processed data can plausibly result in point sources appearing extended. A MCMC analysis demonstrates that the orbit of Fom b is highly eccentric, with e=0.8+/-0.1, a=177+/-68 AU, and q = 32+/-24 AU. Fom b is apsidally aligned with the belt and 90% of allowed orbits have mutual inclination 36 deg or less. Fomalhaut b's orbit is belt-crossing in projection, but only 12% of possible orbits have nodes within a 25 AU wide belt annulus (133-158 AU). The high e invokes a dynamical history where Fom b may have experienced a significant dynamical interaction with a hypothetical planet Fomalhaut c, and the current orbital configuration may be relatively short-lived. The new value for the periastron distance diminishes the Hill radius of Fom b and any weakly bound satellite system surrounding a planet would be sheared and dynamically heated at periapse. We argue that Fom b's minimum mass is that of a dwarf planet in order for a circumplanetary satellite system to remain bound to a sufficient radius from the planet to be consistent with the dust scattered light hypothesis. Fom b may be optically bright because the recent passage through periapse and/or the ascending node has increased the erosion rates of planetary satellites. In the coplanar case, Fomalhaut b will collide with the main belt around 2032, and the subsequent emergent phenomena may help determine its physical nature.Comment: 49 Pages, 33 Figures, 5 Tables; Submitted to ApJ, Dec. 31, 201

Vertical structure of debris discs

The vertical thickness of debris discs is often used as a measure of these systems' dynamical excitation and as clues to the presence of hidden massive perturbers such as planetary embryos. However, this argument could be flawed because the observed dust should be naturally placed on inclined orbits by the combined effect of radiation pressure and mutual collisions. We critically reinvestigate this issue and numerically estimate what the "natural" vertical thickness of a collisionally evolving disc is, in the absence of any additional perturbing body. We use a deterministic collisional code, following the dynamical evolution of a population of indestructible test grains suffering mutual inelastic impacts. Grain differential sizes as well as the effect of radiation pressure are taken into account. We find that, under the coupled effect of radiation pressure and collisions, grains naturally acquire inclinations of a few degrees. The disc is stratified with respect to grain sizes, with the smallest grains having the largest vertical dispersion and the bigger ones clustered closer to the midplane. Debris discs should have a minimum "natural" observed aspect ratio $h_{min}\sim 0.04\pm0.02$ at visible to mid-IR wavelengths where the flux is dominated by the smallest bound grains. These values are comparable to the estimated thicknesses of many vertically resolved debris discs, as is illustrated with the specific example of AU Mic. For all systems with $h \sim h_{min}$, the presence (or absence) of embedded perturbing bodies cannot be inferred from the vertical dispersion of the discComment: accepted for publication in Astronomy and Astrophysics (full abstract in the pdf file

Modeling Self-Subtraction in Angular Differential Imaging: Application to the HD 32297 Debris Disk

We present a new technique for forward-modeling self-subtraction of spatially extended emission in observations processed with angular differential imaging (ADI) algorithms. High-contrast direct imaging of circumstellar disks is limited by quasi-static speckle noise and ADI is commonly used to suppress those speckles. However, the application of ADI can result in self-subtraction of the disk signal due to the disk's finite spatial extent. This signal attenuation varies with radial separation and biases measurements of the disk's surface brightness, thereby compromising inferences regarding the physical processes responsible for the dust distribution. To compensate for this attenuation, we forward-model the disk structure and compute the form of the self-subtraction function at each separation. As a proof of concept, we apply our method to 1.6 and 2.2 micron Keck AO NIRC2 scattered-light observations of the HD 32297 debris disk reduced using a variant of the "locally optimized combination of images" (LOCI) algorithm. We are able to recover disk surface brightness that was otherwise lost to self-subtraction and produce simplified models of the brightness distribution as it appears with and without self-subtraction. From the latter models, we extract radial profiles for the disk's brightness, width, midplane position, and color that are unbiased by self-subtraction. Our analysis of these measurements indicates a break in the brightness profile power law at r~110 AU and a disk width that increases with separation from the star. We also verify disk curvature that displaces the midplane by up to 30 AU towards the northwest relative to a straight fiducial midplane.Comment: Accepted for publication in ApJ, 20 pages, 10 figures, 1 tabl

Close stellar encounters with planetesimal discs: The dynamics of asymmetry in the Beta Pictoris system

We numerically investigate the dynamics of how a close stellar fly-by encounter of a symmetrical circumstellar planetesimal disc can give rise to the many kinds of asymmetries and substructures attributed to the edge-on dusty disc of Beta Pic. In addition we present new optical coronagraphic observations of the outer parts of Beta Pic's disc, and report that the radial extent is significantly greater than was found in previous measurements. The northeasterly extension of the disc's midplane is now measured out to 1835au from the star; the southwesterly component is measured out to 1450au. Hence we use the length asymmetry induced in a distribution of simulation test particles as the principal diagnostic feature when modelling the disc response, in order to constrain fly-by parameters. In particular we favour a low inclination prograde and near-parabolic orbit perturber of mass approximately 0.5 Solar masses. These initial conditions suggest that the perturber could have been physically associated with Beta Pic prior to the encounter. Thus we also consider the possibility that the perturber could be bound to Beta Pic: a consideration also of general interest where dust discs are known to exist in binary star systems. In some of our models, we can relate groupings of perturbed particles to the large-scale structure of the Beta Pic disc. The groupings correspond to: high eccentricity and inclination particles that reach apocentre and maximum height in the southwest, moderately eccentric and low inclination particles that reach apocentre in the northeast, and relatively unperturbed particles inside approximately 200au radius.Comment: Accepted by MNRAS. 15 pages, 19 figures (mainly low resolution). High quality PostScript from http://www.maths.qmw.ac.uk/~jdl