We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

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.

Planet yield calculations may be used to inform the target selection strategy and science operations of space observatories. Forthcoming and proposed NASA missions, such as the Wide-Field Infrared Survey Telescope (WFIRST), the Habitable Exoplanet Imaging Mission (HabEx), and the Large UV/Optical/IR Surveyor (LUVOIR), are expected to be equipped with sensitive coronagraphs and/or starshades. We are developing a suite of numerical simulations to quantify the extent to which ground-based radial velocity (RV) surveys could boost the detection efficiency of direct imaging missions. In this paper, we discuss the first step in the process of estimating planet yields: generating synthetic planetary systems consistent with observed occurrence rates from multiple detection methods. In an attempt to self-consistently populate stars with orbiting planets, it is found that naive extrapolation of occurrence rates (mass, semi-major axis) results in an unrealistically large number-density of Neptune-mass planets beyond the ice-line ($a \gtrsim 5$au), causing dynamic interactions that would destabilize orbits. We impose a stability criterion for multi-planet systems based on mutual Hill radii separation. Considering the influence of compact configurations containing Jovian-mass and Neptune-mass planets results in a marked suppression in the number of terrestrial planets that can exist at large radii. This result has a pronounced impact on planet yield calculations particularly in regions accessible to high-contrast imaging and microlensing. The dynamically compact configurations and occurrence rates that we develop may be incorporated as input into joint RV and direct imaging yield calculations to place meaningful limits on the number of detectable planets with future missions.

Sub-stellar objects exhibit photometric variability, which is believed to be caused by a number of processes, such as magnetically-driven spots or inhomogeneous cloud coverage. Recent models have shown that turbulent flows and waves, including internal gravity waves, may play an important role in cloud evolution. The aim of this paper is to investigate the effect of IGW on dust nucleation and dust growth, and whether observations of the resulting cloud structures could be used to recover atmospheric density information. For a simplified atmosphere in two dimensions, we numerically solved the governing fluid equations to simulate the effect on dust nucleation and mantle growth as a result of the passage of an IGW. Furthermore, we derived an expression that relates the properties of the wave-induced cloud structures to observable parameters in order to deduce the atmospheric density. Numerical simulations show that the $\rho, p, T$ variations caused by gravity waves lead to an increase of the nucleation rate by up to a factor 20, and an increase of the mantle growth rate by up to a factor 1.6, compared to their equilibrium values. An exploration of the wider parameter space shows that in absolute terms, the increase in nucleation due to IGW is stronger in cooler (T dwarfs) and TiO2-rich sub-stellar atmospheres. The relative increase, however, is greater in warmer (L dwarf) and TiO2-poor atmospheres due to conditions less suited for efficient nucleation at equilibrium. These variations lead to banded areas in which dust formation is much more pronounced, similar to the cloud structures observed on Earth. We show that IGW in the atmosphere of sub-stellar objects can produce banded clouds structures similar to that observed on Earth. We propose a method with which potential observations of banded clouds could be used to estimate the atmospheric density of sub-stellar objects.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

We present new high fidelity optical coronagraphic imagery of the inner $\sim$50 au of AU Mic's edge-on debris disk using the BAR5 occulter of the Hubble Space Telescope Imaging Spectrograph (HST/STIS) obtained on 26-27 July 2018. This new imagery reveals that "feature A", residing at a projected stellocentric separation of 14.2 au on SE-side of the disk, exhibits an apparent "loop-like" morphology at the time of our observations. The loop has a projected width of 1.5 au and rises 2.3 au above the disk midplane. We also explored TESS photometric observations of AU Mic that are consistent with evidence of two starspot complexes in the system. The likely co-alignment of the stellar and disk rotational axes breaks degeneracies in detailed spot modeling, indicating that AU Mic's projected magnetic field axis is offset from its rotational axis. We speculate that small grains in AU Mic's disk could be sculpted by a time-dependent wind that is influenced by this offset magnetic field axis, analogous to co-rotating Solar interaction regions that sculpt and influence the inner and outer regions of our own Heliosphere. Alternatively, if the observed spot modulation is indicative of a significant mis-alignment of the stellar and disk rotational axes, we suggest the disk could still be sculpted by the differential equatorial versus polar wind that it sees with every stellar rotation.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

We have obtained Hubble Space Telescope STIS and NICMOS, and Gemini/GPI scattered light images of the HD 191089 debris disk. We identify two spatial components: a ring resembling Kuiper Belt in radial extent (FWHM: ${\sim}$25 au, centered at ${\sim}$46 au), and a halo extending to ${\sim}$640 au. We find that the halo is significantly bluer than the ring, consistent with the scenario that the ring serves as the "birth ring" for the smaller dust in the halo. We measure the scattering phase functions in the 30{\deg}-150{\deg} scattering angle range and find the halo dust is both more forward- and backward-scattering than the ring dust. We measure a surface density power law index of -0.68${\pm}$0.04 for the halo, which indicates the slow-down of the radial outward motion of the dust. Using radiative transfer modeling, we attempt to simultaneously reproduce the (visible) total and (near-infrared) polarized intensity images of the birth ring. Our modeling leads to mutually inconsistent results, indicating that more complex models, such as the inclusion of more realistic aggregate particles, are needed.

Instrumentation techniques in the field of direct imaging of exoplanets have greatly advanced over the last two decades. Two of the four NASA-commissioned large concept studies involve a high-contrast instrument for the imaging and spectral characterization of exo-Earths from space: LUVOIR and HabEx. This whitepaper describes the status of 8 optical testbeds in the US and France currently in operation to experimentally validate the necessary technologies to image exo-Earths from space. They explore two complementary axes of research: (i) coronagraph designs and manufacturing and (ii) active wavefront correction methods and technologies. Several instrument architectures are currently being analyzed in parallel to provide more degrees of freedom for designing the future coronagraphic instruments. The necessary level of performance has already been demonstrated in-laboratory for clear off-axis telescopes (HabEx-like) and important efforts are currently in development to reproduce this accomplishment on segmented and/or on-axis telescopes (LUVOIR-like) over the next two years.

As observational datasets become larger and more complex, so too are the questions being asked of these data. Data simulations, i.e., synthetic data with properties (pixelization, noise, PSF, artifacts, etc.) akin to real data, are therefore increasingly required for several purposes, including: (1) testing complicated measurement methods, (2) comparing models and astrophysical simulations to observations in a manner that requires as few assumptions about the data as possible, (3) predicting observational results based on models and astrophysical simulations for, e.g., proposal planning, and (4) mitigating risk for future observatories and missions by effectively priming and testing pipelines. We advocate for an increase in using synthetic data to plan for and interpret real observations as a matter of routine. This will require funding for (1) facilities to provide robust data simulators for their instruments, telescopes, and surveys, and (2) making synthetic data publicly available in archives (much like real data) so as to lower the barrier of entry to all.

Cold debris disks (T$<$200 K) are analogues to the dust in the Solar System's Kuiper belt--dust generated from the evaporation and collision of minor bodies perturbed by planets, our Sun, and the local interstellar medium. Scattered light from debris disks acts as both a signpost for unseen planets as well as a source of contamination for directly imaging terrestrial planets, but many details of these disks are poorly understood. We lay out a critical observational path for the study of nearby debris disks that focuses on defining an empirical relationship between scattered light and thermal emission from a disk, probing the dynamics and properties of debris disks, and directly determining the influence of planets on disks. We endorse the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe. This white paper extends and complements the material presented therein with a focus on debris disks around nearby stars. Separate complementary papers are being submitted regarding the inner warm regions of debris disks (Mennesson et al.), the modeling of debris disk evolution (Gaspar et al.), studies of dust properties (Chen et al.), and thermal emission from disks (Su et al.).

The expected yield of potentially Earth-like planets is a useful metric for designing future exoplanet-imaging missions. Recent yield studies of direct-imaging missions have focused primarily on yield methods and trade studies using "toy" models of missions. Here we increase the fidelity of these calculations substantially, adopting more realistic exoplanet demographics as input, an improved target list, and a realistic distribution of exozodi levels. Most importantly, we define standardized inputs for instrument simulations, use these standards to directly compare the performance of realistic instrument designs, include the sensitivity of coronagraph contrast to stellar diameter, and adopt engineering-based throughputs and detector parameters. We apply these new high-fidelity yield models to study several critical design trades: monolithic vs segmented primary mirrors, on-axis vs off-axis secondary mirrors, and coronagraphs vs starshades. We show that as long as the gap size between segments is sufficiently small, there is no difference in yield for coronagraph-based missions with monolithic off-axis telescopes and segmented off-axis telescopes, assuming that the requisite engineering constraints imposed by the coronagraph can be met in both scenarios. We show that there is currently a factor of ~2 yield penalty for coronagraph-based missions with on-axis telescopes compared to off-axis telescopes, and note that there is room for improvement in coronagraph designs for on-axis telescopes. We also reproduce previous results in higher fidelity showing that the yields of coronagraph-based missions continue to increase with aperture size while the yields of starshade-based missions turnover at large apertures if refueling is not possible. Finally, we provide absolute yield numbers with uncertainties that include all major sources of astrophysical noise to guide future mission design.

A coronagraphic starlight suppression system situated on a future flagship space observatory offers a promising avenue to image Earth-like exoplanets and search for biomarkers in their atmospheric spectra. One NASA mission concept that could serve as the platform to realize this scientific breakthrough is the Large UV/Optical/IR Surveyor (LUVOIR). Such a mission would also address a broad range of topics in astrophysics with a multiwavelength suite of instruments. The apodized pupil Lyot coronagraph (APLC) is one of several coronagraph design families that the community is assessing as part of NASAs Exoplanet Exploration Program Segmented aperture coronagraph design and analysis (SCDA) team. The APLC is a Lyot-style coronagraph that suppresses starlight through a series of amplitude operations on the on-axis field. Given a suite of seven plausible segmented telescope apertures, we have developed an object-oriented software toolkit to automate the exploration of thousands of APLC design parameter combinations. This has enabled us to empirically establish relationships between planet throughput and telescope aperture geometry, inner working angle, bandwidth, and contrast level. In parallel with the parameter space exploration, we have investigated several strategies to improve the robustness of APLC designs to fabrication and alignment errors. We also investigate the combination of APLC with wavefront control or complex focal plane masks to improve inner working angle and throughput. Preliminary scientific yield evaluations based on design reference mission simulations indicate the APLC is a very competitive concept for surveying the local exoEarth population with a mission like LUVOIR.

Future space-based direct imaging missions are poised to search for biosignatures in the atmospheres of potentially habitable planets orbiting nearby AFGKM stars. Although these missions could conduct a survey of high-priority target stars to detect candidate Earth-like planets, conducting a precursor radial velocity (RV) survey will benefit future direct imaging missions in four ways. First, an RV survey capable of detecting signals as small as 8 cm/s over timescales of a few years could discover potentially habitable Earth-mass planets orbiting dozens of nearby GKM stars accessible to space-based direct imaging. Second, RVs will improve scheduling efficiency by reducing the required number of revisits for orbit determination, and revealing when a planet of interest is most observable. Third, RV observations will reveal the masses of gas and ice giants that could be mistaken for Earth-mass planets, thereby reducing the time spent identifying false positives. Fourth, mass measurements from RVs will provide the surface gravities necessary for interpreting atmospheric spectra and potential biosignatures.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

The small class of known stars with unusually warm, dusty debris disks is a key sample to probe in order to understand cascade models and extreme collisions that likely lead to the final configurations of planetary systems. Because of its extreme dustiness and small radius, the disk of BD +20 307 has a short predicted collision time and is therefore an interesting target in which to look for changes in dust quantity and composition over time. To compare with previous ground and Spitzer Space Telescope data, SOFIA photometry and spectroscopy were obtained. The system's 8.8-12.5 $\mu$m infrared emission increased by $10 \pm 2 \%$ over nine years between the SOFIA and earlier Spitzer measurements. In addition to an overall increase in infrared excess, there is a suggestion of a greater increase in flux at shorter wavelengths (less than 10.6 $\mu$m) compared to longer wavelengths (greater than 10.6 $\mu$m). Steady-state collisional cascade models cannot explain the increase in BD +20 307's disk flux over such short timescales. A catastrophic collision between planetary-scale bodies is still the most likely origin for the system's extreme dust; however, the cause for its recent variation requires further investigation.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

For the first time in human history, technologies have matured sufficiently to enable a mission capable of discovering and characterizing habitable planets like Earth orbiting sunlike stars other than the Sun. At the same time, such a platform would enable unique science not possible from ground-based facilities. This science is broad and exciting, ranging from new investigations of our own solar system to a full range of astrophysics disciplines. The Habitable Exoplanet Observatory, or HabEx, is one of four studies currently being undertaken by NASA in preparation for the 2020 Astrophysics Decadal Survey. HabEx has been designed to be the Great Observatory of the 2030s, with community involvement through a competed and funded Guest Observer (GO) program. This interim report describes the HabEx baseline concept, which is a space-based 4-meter diameter telescope mission concept with ultraviolet (UV), optical, and near-infrared (near-IR) imaging and spectroscopy capabilities. More information on HabEx can be found at https://www.jpl.nasa.gov/habex

Coronagraph instruments on future space telescopes will enable the direct detection and characterization of Earth-like exoplanets around Sun-like stars for the first time. The quest for the optimal optical coronagraph designs has made rapid progress in recent years thanks to the Segmented Coronagraph Design and Analysis (SCDA) initiative led by the Exoplanet Exploration Program at NASA's Jet Propulsion Laboratory. As a result, several types of high-performance designs have emerged that make use of dual deformable mirrors to (1) correct for optical aberrations and (2) suppress diffracted starlight from obstructions and discontinuities in the telescope pupil. However, the algorithms used to compute the optimal deformable mirror surface tend to be computationally intensive, prohibiting large scale design surveys. Here, we utilize the Fast Linearized Coronagraph Optimizer (FALCO), a tool that allows for rapid optimization of deformable mirror shapes, to explore trade-offs in coronagraph designs for obstructed and segmented space telescopes. We compare designs for representative shaped pupil Lyot and vortex coronagraphs, two of the most promising concepts for the LUVOIR space mission concept. We analyze the optical performance of each design, including their throughput and ability to passively suppress light from partially resolved stars in the presence of low-order aberrations. Our main result is that deformable mirror based apodization can sufficiently suppress diffraction from support struts and inter-segment gaps whose widths are on the order of $\sim$0.1% of the primary mirror diameter to detect Earth-sized planets within a few tens of milliarcseconds from the star.

We present the most distant detection of cosmic voids ($z \sim 2.3$) and the first detection of three-dimensional voids in the Lyman-$\alpha$ forest. We used a 3D tomographic map of the absorption with effective comoving spatial resolution of $2.5\,h^{-1}\mathrm{Mpc}$ and volume of $3.15\times 10^5\,h^{-3}\mathrm{Mpc}^3$, which was reconstructed from moderate-resolution Keck-I/LRIS spectra of 240 background Lyman-break galaxies and quasars in a $0.16\,\mathrm{deg}^2$ footprint in the COSMOS field. Voids were detected using a spherical overdensity finder calibrated from hydrodynamical simulations of the intergalactic medium. This allows us to identify voids in the IGM corresponding to voids in the underlying matter density field, yielding a consistent volume fraction of voids in both data (19.5%) and simulations (18.2%). We fit excursion set models to the void radius function and compare the radially-averaged stacked profiles of large voids ($r > 5$ $h^{-1}$ Mpc) to stacked voids in mock observations and the simulated density field. Comparing with 432 coeval galaxies with spectroscopic redshifts in the same volume as the tomographic map, we find that the tomography-identified voids are underdense in galaxies by 5.95$\sigma$ compared to random cells.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

Future NASA concept missions that are currently under study, like Habitable Exoplanet Imaging Mission (HabEx) & Large Ultra-Violet Optical Infra Red (LUVOIR) Surveyor, would discover a large diversity of exoplanets. We propose here a classification scheme that distinguishes exoplanets into different categories based on their size and incident stellar flux, for the purpose of providing the expected number of exoplanets observed (yield) with direct imaging missions. The boundaries of this classification can be computed using the known chemical behavior of gases and condensates at different pressures and temperatures in a planetary atmosphere. In this study, we initially focus on condensation curves for sphalerite ZnS, H2O, CO2 and CH4. The order in which these species condense in a planetary atmosphere define the boundaries between different classes of planets. Broadly, the planets are divided into rocky (0.5 - 1.0RE), super-Earths (1.0- 1.75RE), sub-Neptunes (1.75-3.5RE), sub-Jovians (3.5 - 6.0RE) and Jovians (6-14.3RE) based on their planet sizes, and 'hot', 'warm' and 'cold' based on the incident stellar flux. We then calculate planet occurrence rates within these boundaries for different kinds of exoplanets, \eta_{planet}, using the community co-ordinated results of NASA's Exoplanet Program Analysis Group's Science Analysis Group-13 (SAG-13). These occurrence rate estimates are in turn used to estimate the expected exoplanet yields for direct imaging missions of different telescope diameter.

We present the first scattered-light images of two debris disks around the F8 star HD 104860 and the F0V star HD 192758, respectively $\sim45$ and $\sim67$ pc away. We detected these systems in the F110W and F160W filters through our re-analysis of archival Hubble Space Telescope NICMOS data with modern starlight subtraction techniques. Our image of HD 104860 confirms the morphology previously observed by Herschel in thermal emission with a well-defined ring at radius $\sim114$ au inclined $\sim58$ degrees. Although the outer edge profile is consistent with dynamical evolution models, the sharp inner edge suggests sculpting by unseen perturbers. Our images of HD 192758 reveal a disk at radius $\sim95$ au inclined by $\sim59$ degrees, never resolved so far. These disks have low scattering albedos of 10% and 13% respectively, inconsistent with water ice grain compositions. They are reminiscent of several other disks with similar inclination and scattering albedos: Fomalhaut, HD 92945, HD 202628, and HD 207129. They are also very distinct from brighter disks in the same inclination bin, which point to different compositions between these two populations. Varying scattering albedo values can be explained by different grain porosities, chemical compositions, or grain size distributions, which may indicate distinct formation mechanisms or dynamical processes at work in these systems. Finally, these faint disks with large infrared excesses may be representative of an underlying population of systems with low albedo values. Searches with more sensitive instruments on HST or on the James Webb Space Telescope and using state-of-the art starlight-subtraction methods may help discover more of such faint systems.

In contemporary sub-stellar model atmospheres, dust growth occurs through neutral gas-phase surface chemistry. Recently, there has been a growing body of theoretical and observational evidence suggesting that ionisation processes can also occur. As a result, atmospheres are populated by regions composed of plasma, gas and dust, and the consequent influence of plasma processes on dust evolution is enhanced. This paper aims to introduce a new model of dust growth and destruction in sub-stellar atmospheres via plasma deposition and plasma sputtering. Using example sub-stellar atmospheres from Drift-Phoenix, we have compared plasma deposition and sputtering timescales to those from neutral gas-phase surface chemistry to ascertain their regimes of influence. We calculated the plasma sputtering yield and discuss the circumstances where plasma sputtering dominates over deposition. Within the highest dust density cloud regions, plasma deposition and sputtering dominates over neutral gas-phase surface chemistry if the degree of ionisation is $\gtrsim10^{-4}$. Loosely bound grains with surface binding energies of the order of $0.1-1$ eV are susceptible to destruction through plasma sputtering for feasible degrees of ionisation and electron temperatures; whereas, strong crystalline grains with binding energies of the order $10$ eV are resistant to sputtering. The mathematical framework outlined sets the foundation for the inclusion of plasma deposition and plasma sputtering in global dust cloud formation models of sub-stellar atmospheres.

The segmented coronagraph design and analysis (SCDA) study is a coordinated effort, led by Stuart Shaklan (JPL) and supported by NASA's Exoplanet Exploration Program (ExEP), to provide efficient coronagraph design concepts for exoplanet imaging with future segmented aperture space telescopes. This document serves as an update on the apodized vortex coronagraph designs devised by the Caltech/JPL SCDA team. Apodized vortex coronagraphs come in two flavors, where the apodization is achieved either by use of 1) a gray-scale semi-transparent pupil mask or 2) a pair of deformable mirrors in series. Each approach has attractive benefits. This document presents a comprehensive review of the former type. Future theoretical investigations will further explore the use of deformable mirrors for apodization.

We present two state-of-the-art models of the solar system, one corresponding to the present day and one to the Archean Eon 3.5 billion years ago. Each model contains spatial and spectral information for the star, the planets, and the interplanetary dust, extending to 50 AU from the sun and covering the wavelength range 0.3 to 2.5 micron. In addition, we created a spectral image cube representative of the astronomical backgrounds that will be seen behind deep observations of extrasolar planetary systems, including galaxies and Milky Way stars. These models are intended as inputs to high-fidelity simulations of direct observations of exoplanetary systems using telescopes equipped with high-contrast capability. They will help improve the realism of observation and instrument parameters that are required inputs to statistical observatory yield calculations, as well as guide development of post-processing algorithms for telescopes capable of directly imaging Earth-like planets.

We have conducted a mini-survey for low-frequency radio emission from some of the closest brown dwarfs to the Sun with rapid rotation rates: SIMP J013656.5+093347, WISEPC J150649.97+702736.0, and WISEPA J174124.26+255319.5. We have placed robust 3-sigma upper limits on the flux density in the 111 - 169 MHz frequency range for these targets: WISE 1506: < 0.72 mJy; WISE 1741: < 0.87 mJy; SIMP 0136: < 0.66 mJy. At 8 hours of integration per target to achieve these limits, we find that systematic and detailed study of this class of object at LOFAR frequencies will require a substantial dedication of resources.

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.

An increasing number of observations have shown that gaseous debris discs are not an exception. However, until now we only knew of cases around A stars. Here we present the first detection of 12CO (2-1) disc emission around an F star, HD 181327, obtained with ALMA observations at 1.3 mm. The continuum and CO emission are resolved into an axisymmetric disc with ring-like morphology. Using a Markov chain Monte Carlo method coupled with radiative transfer calculations we study the dust and CO mass distribution. We find the dust is distributed in a ring with a radius of 86.0 +- 0.4 AU and a radial width of 23.2 +- 1.0 AU. At this frequency the ring radius is smaller than in the optical, revealing grain size segregation expected due to radiation pressure. We also report on the detection of low level continuum emission beyond the main ring out to ~200 AU. We model the CO emission in the non-LTE regime and we find that the CO is co-located with the dust, with a total CO gas mass ranging between 1.2x10^-6 Mearth and 2.9x10^-6 Mearth, depending on the gas kinetic temperature and collisional partners densities. The CO densities and location suggest a secondary origin, i.e. released from icy planetesimals in the ring. We derive a CO+CO2 cometary composition that is consistent with Solar system comets. Due to the low gas densities it is unlikely that the gas is shaping the dust distribution.

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.

ATLAST is a particular realization of LUVOIR, the Large Ultraviolet Optical Infrared telescope, a ~10 m diameter space telescope being defined for consideration in the 2020 Decadal Review of astronomy and astrophysics. ATLAST/LUVOIR is generally thought of as an ambient temperature (~300 K) system, and little consideration has been given to using it at infrared wavelengths longward of ~2 {\mu}m. We assess the scientific and technical benefits of operating such a telescope further into the infrared, with particular emphasis on the study of exoplanets, which is a major science theme for ATLAST/LUVOIR. For the study of exoplanet atmospheres, the capability to work at least out to 5.0 {\mu}m is highly desirable. Such an extension of the long wavelength limit of ATLAST would greatly increase its capabilities for studies of exoplanet atmospheres and provide powerful capabilities for the study of a wide range of astrophysical questions. We present a concept for a fiber-fed grating spectrometer which would enable R = 200 spectroscopy on ATLAST with minimal impact on the other focal planet instruments. We conclude that it is technically feasible and highly desirable scientifically to extend the wavelength range of ATLAST to at least 5 {\mu}m.

Brown dwarfs and giant gas extrasolar planets have cold atmospheres with a rich chemical compositions from which mineral cloud particles form. Their properties, like particle sizes and material composition, vary with height, and the mineral cloud particles are charged due to triboelectric processes in such dynamic atmospheres. The dynamics of the atmospheric gas is driven by the irradiating host star and/or by the rotation of the objects that changes during its lifetime. Thermal gas ionisation in these ultra-cool but dense atmospheres allows electrostatic interactions and magnetic coupling of a substantial atmosphere volume. Combined with a strong magnetic field $\gg B_{\rm Earth}$, a chromosphere and aurorae might form as suggested by radio and X-ray observations of brown dwarfs. Non-equilibrium processes like cosmic ray ionisation and discharge processes in clouds will increase the local pool of free electrons in the gas. Cosmic rays and lighting discharges also alter the composition of the local atmospheric gas such that tracer molecules might be identified. Cosmic rays affect the atmosphere through air showers which was modelled with a 3D Monte Carlo radiative transfer code to be able to visualise their spacial extent. Given a certain degree of thermal ionisation of the atmospheric gas, we suggest that electron attachment to charge mineral cloud particles is too inefficient to cause an electrostatic disruption of the cloud particles. Cloud particles will therefore not be destroyed by Coulomb explosion for the local temperature in the collisional dominated brown dwarf and giant gas planet atmospheres. However, the cloud particles are destroyed electrostatically in regions with strong gas ionisation. The potential size of such cloud holes would, however, be too small and might occur too far inside the cloud to mimic the effect of, e.g., magnetic field induced star spots.

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.

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 $\lambda_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 $\eta_{\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.

This is a joint summary of the reports from the three Astrophysics Program Analysis Groups (PAGs) in response to the "Planning for the 2020 Decadal Survey" charge given by the Astrophysics Division Director Paul Hertz. This joint executive summary contains points of consensus across all three PAGs. Additional findings specific to the individual PAGs are reported separately in the individual reports. The PAGs concur that all four large mission concepts identified in the white paper as candidates for maturation prior to the 2020 Decadal Survey should be studied in detail. These include the Far-IR Surveyor, the Habitable-Exoplanet Imaging Mission, the UV/Optical/IR Surveyor, and the X-ray Surveyor. This finding is predicated upon assumptions outlined in the white paper and subsequent charge, namely that 1) major development of future large flagship missions under consideration are to follow the implementation phases of JWST and WFIRST; 2) NASA will partner with the European Space Agency on its L3 Gravitational Wave Surveyor; 3) the Inflation Probe be classified as a probe-class mission to be developed according to the 2010 Decadal Survey report. If these key assumptions were to change, this PAG finding would need to be re-evaluated. The PAGs find that there is strong community support for the second phase of this activity - maturation of the four proposed mission concepts via Science and Technology Definition Teams (STDTs). The PAGs find that there is strong consensus that all of the STDTs contain broad and interdisciplinary representation of the science community. Finally, the PAGs find that there is community support for a line of Probe-class missions within the Astrophysics mission portfolio (condensed).

Results are presented of a first study of collisionless magnetic reconnection starting from a recently found exact nonlinear force-free Vlasov-Maxwell equilibrium. The initial state has a Harris sheet magnetic field profile in one direction and a non-uniform guide field in a second direction, resulting in a spatially constant magnetic field strength as well as a constant initial plasma density and plasma pressure. It is found that the reconnection process initially resembles guide field reconnection, but that a gradual transition to anti-parallel reconnection happens as the system evolves. The time evolution of a number of plasma parameters is investigated, and the results are compared with simulations starting from a Harris sheet equilibrium and a Harris sheet plus constant guide field equilibrium.

We study the statistics of the Lyman-$\alpha$ forest in a flat LCDM cosmology with the N-body + Eulerian hydrodynamics code Nyx. We produce a suite of simulations, covering the observationally relevant redshift range $2 \leq z \leq 4$. We find that a grid resolution of 20 kpc/h is required to produce one percent convergence of Lyman-$\alpha$ flux statistics, up to k = 10 h/Mpc. In addition to establishing resolution requirements, we study the effects of missing modes in these simulations, and find that box sizes of L > 40 Mpc/h are needed to suppress numerical errors to a sub-percent level. Our optically-thin simulations with the ionizing background prescription of Haardt & Madau (2012) reproduce an IGM equation of state with $T_0 \approx 10^4 K$ and $\gamma \approx 1.55$ at z=2, with a mean transmitted flux close to the observed values. When using the ionizing background prescription of Faucher-Giguere et al. (2009), the mean flux is 10-15 per cent below observed values at z=2, and a factor of 2 too small at z = 4. We show the effects of the common practice of rescaling optical depths to the observed mean flux and how it affects convergence rates. We also investigate the common practice of `splicing' results from a number of different simulations to estimate the 1D flux power spectrum and show it is accurate at the 10 percent level. Finally, we find that collisional heating of the gas from dark matter particles is negligible in modern cosmological simulations.

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.

The well-studied M9 dwarf TVLM 513-46546 is a rapid rotator (P_rot ~ 2 hr) hosting a stable, dipolar magnetic field of ~3 kG surface strength. Here we report its detection with ALMA at 95 GHz at a mean flux density of $56 \pm 12$ uJy, making it the first ultracool dwarf detected in the millimeter band, excluding young, disk-bearing objects. We also report flux density measurements from unpublished archival VLA data and new optical monitoring data from the Liverpool Telescope. The ALMA data are consistent with a power-law radio spectrum that extends continuously between centimeter and millimeter wavelengths. We argue that the emission is due to the synchrotron process, excluding thermal, free-free, and electron cyclotron maser emission as possible sources. During the interval of the ALMA observation that phases with the maximum of the object's optical variability, the flux density is higher at a ~1.8 sigma significance level. These early results show how ALMA opens a new window for studying the magnetic activity of ultracool dwarfs, particularly shedding light on the particle acceleration mechanism operating in their immediate surroundings.

The appearance of debris disks around distant stars depends upon the scattering/phase function (SPF) of the material in the disk. However, characterizing the SPFs of these extrasolar debris disks is challenging because only a limited range of scattering angles are visible to Earth-based observers. By contrast, Saturn's tenuous rings can be observed over a much broader range of geometries, so their SPFs can be much better constrained. Since these rings are composed of small particles released from the surfaces of larger bodies, they are reasonable analogs to debris disks and so their SPFs can provide insights into the plausible scattering properties of debris disks. This work examines two of Saturn's dusty rings: the G ring (at 167,500 km from Saturn's center) and the D68 ringlet (at 67,600 km). Using data from the cameras onboard the Cassini spacecraft, we are able to estimate the rings' brightnesses at scattering angles ranging from 170 to 0.5 degrees. We find that both of the rings exhibit extremely strong forward-scattering peaks, but for scattering angles above 60 degrees their brightnesses are nearly constant. These SPFs can be well approximated by a linear combination of three Henyey-Greenstein functions, and are roughly consistent with the SPFs of irregular particles from laboratory measurements. Comparing these data to Fraunhofer and Mie models highlights several challenges involved in extracting information about particle compositions and size distributions from SPFs alone. The SPFs of these rings also indicate that the degree of forward scattering in debris disks may be greatly underestimated.

Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars at $z\sim 2.3-3$ within a $11.5'\times13.5'$ area of the COSMOS field ($\sim 1200\,\mathrm{deg}^2$ projected area density or $\sim 2.4\,h^{-1}\,\mathrm{Mpc}$ mean transverse separation), we reconstruct a 3D tomographic map of the foreground Ly$\alpha$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $\sigma_{3d}\approx3.5\,h^{-1}\,\mathrm{Mpc}$ comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the tomographic map. We find an extended IGM overdensity with deep absorption troughs at $z=2.45$ associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be $M_{\rm dm} (z=2.45) = (9\pm4)\times 10^{13}\,h^{-1}\,\mathrm{M_\odot}$, and argue based on this mass and absorption strength that it will form at least one $z\sim0$ galaxy cluster with $M(z=0) = (3\pm 2) \times 10^{14}\,h^{-1}\mathrm{M_\odot}$, although its elongated nature suggests that it will likely collapse into two separate clusters. We also point out a compact overdensity of six MOSDEF galaxies at $z=2.30$ within a $r\sim 1\,h^{-1}\,\mathrm{Mpc}$ radius and $\Delta z\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$\alpha$ forest tomography on larger volumes to study galaxy properties as a function of environment, as well as revealing the large-scale IGM overdensities associated with protoclusters and other features of large-scale structure.

Radio and X-ray emission from brown dwarfs suggest that an ionised gas and a magnetic field with a sufficient flux density must be present. We perform a reference study for late M-dwarfs, brown dwarfs and giant gas planet to identify which ultra-cool objects are most susceptible to plasma and magnetic processes. Only thermal ionisation is considered. We utilise the {\sc Drift-Phoenix} model grid where the local atmospheric structure is determined by the global parameters T$_{\rm eff}$, $\log(g)$ and [M/H]. Our results show that it is not unreasonable to expect H$_{\alpha}$ or radio emission to origin from Brown Dwarf atmospheres as in particular the rarefied upper parts of the atmospheres can be magnetically coupled despite having low degrees of thermal gas ionisation. Such ultra-cool atmospheres could therefore drive auroral emission without the need for a companion's wind or an outgassing moon. The minimum threshold for the magnetic flux density required for electrons and ions to be magnetised is well above typical values of the global magnetic field of a brown dwarf and a giant gas planet. Na$^{+}$, K$^{+}$ and Ca$^{+}$ are the dominating electron donors in low-density atmospheres (low log(g), solar metallicity) independent of T$_{\rm eff}$. Mg$^{+}$ and Fe$^{+}$ dominate the thermal ionisation in the inner parts of M-dwarf atmospheres. Molecules remain unimportant for thermal ionisation. Chemical processes (e.g. cloud formation) affecting the most abundant electron donors, Mg and Fe, will have a direct impact on the state of ionisation in ultra-cool atmospheres.

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; $\lambda=0.4-1.8~\mu\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.

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.

Recent observations of brown dwarf spectroscopic variability in the infrared infer the presence of patchy cloud cover. This paper proposes a mechanism for producing inhomogeneous cloud coverage due to the depletion of cloud particles through the Coulomb explosion of dust in atmospheric plasma regions. Charged dust grains Coulomb-explode when the electrostatic stress of the grain exceeds its mechanical tensile stress, which results in grains below a critical radius $a<a^{\rm Coul}_{\rm crit}$ being broken up. This work outlines the criteria required for the Coulomb explosion of dust clouds in substellar atmospheres, the effect on the dust particle size distribution function, and the resulting radiative properties of the atmospheric regions. Our results show that for an atmospheric plasma region with an electron temperature of $T_{e}=10$~eV ($\approx10^{5}$~K), the critical grain radius varies from $10^{-7}$ to $10^{-4}$~cm, depending on the grains' tensile strength. Higher critical radii up to $10^{-3}$~cm are attainable for higher electron temperatures. We find that the process produces a bimodal particle size distribution composed of stable nanoscale seed particles and dust particles with $a\geq a^{\rm Coul}_{\rm crit}$, with the intervening particle sizes defining a region devoid of dust. As a result, the dust population is depleted, and the clouds become optically thin in the wavelength range $0.1-10~\mu$m, with a characteristic peak that shifts to higher wavelengths as more sub-micrometer particles are destroyed. In an atmosphere populated with a distribution of plasma volumes, this will yield regions of contrasting radiative properties, thereby giving a source of inhomogeneous cloud coverage. The results presented here may also be relevant for dust in supernova remnants and protoplanetary disks.

We present a new method of finding cosmic voids using tomographic maps of Ly{\alpha} forest flux. We identify cosmological voids with radii of 2 - 12 $h^{-1}$Mpc in a large N-body simulation at $z = 2.5$, and characterize the signal of the high-redshift voids in density and Ly{\alpha} forest flux. The void properties are similar to what has been found at lower redshifts, but they are smaller and have steeper radial density profiles. Similarly to what has been found for low-redshift voids, the radial velocity profiles have little scatter and agree very well with the linear theory prediction. We run the same void finder on an ideal Ly{\alpha} flux field and tomographic reconstructions at various spatial samplings. We compare the tomographic map void catalogs to the density void catalog and find good agreement even with modest-sized voids ($r > 6 \, h^{-1}$Mpc). Using our simple void-finding method, the configuration of the ongoing CLAMATO survey covering 1 deg$^2$ would provide a sample of about 100 high-redshift voids. We also provide void-finding forecasts for larger area surveys, and discuss how these void samples can be used to test modified gravity models, study high-redshift void galaxies, and to make an Alcock-Paczynski measurement. To aid future work in this area, we provide public access to our simulation products, catalogs, and sample tomographic flux maps.

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.

We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.

We present a new method of finding protoclusters using tomographic maps of Ly$\alpha$ Forest flux. We review our method of creating tomographic flux maps and discuss our new high performance implementation, which makes large reconstructions computationally feasible. Using a large N-body simulation, we illustrate how protoclusters create large-scale flux decrements, roughly 10 $h^{-1}$Mpc across, and how we can use this signal to find them in flux maps. We test the performance of our protocluster finding method by running it on the ideal, noiseless map and tomographic reconstructions from mock surveys, and comparing to the halo catalog. Using the noiseless map, we find protocluster candidates with about 90% purity, and recover about 75% of the protoclusters that form massive clusters ($> 3 \times 10^{14} \, h^{-1} M_{\odot}$). We construct mock surveys similar to the ongoing COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) survey. While the existing data has an average sightline separation of 2.3 $h^{-1}$Mpc, we test separations of 2 - 6 $h^{-1}$Mpc to see what can be tolerated for our application. Using reconstructed maps from small separation mock surveys, the protocluster candidate purity and completeness are very close what was found in the noiseless case. As the sightline separation increases, the purity and completeness decrease, although they remain much higher than we initially expected. We extended our test cases to mock surveys with an average separation of 15 $h^{-1}$Mpc, meant to reproduce high source density areas of the BOSS survey. We find that even with such a large sightline separation, the method can still be used to find some of the largest protoclusters.

[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.

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.

We present the first observations of foreground Lyman-$\alpha$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$\alpha$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$\alpha$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

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.

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).

We explore the feasibility of using the atmosphere of Jupiter to detect Ultra-High-Energy Cosmic Rays (UHECR's). The large surface area of Jupiter allows us to probe cosmic rays of higher energies than previously accessible. Cosmic ray extensive air showers in Jupiter's atmosphere could in principle be detected by the Large Area Telescope (LAT) on the Fermi observatory. In order to be observed, these air showers would need to be oriented toward the Earth, and would need to occur sufficiently high in the atmosphere that the gamma rays can penetrate. We demonstrate that, under these assumptions, Jupiter provides an effective cosmic ray "detector" area of $3.3 \times 10^7$ km$^2$. We predict that Fermi-LAT should be able to detect events of energy $>10^{21}$ eV with fluence $10^{-7}$ erg cm$^{-2}$ at a rate of about one per month. The observed number of air showers may provide an indirect measure of the flux of cosmic rays $\gtrsim 10^{20}$ eV. Extensive air showers also produce a synchrotron signature that may be measurable by ALMA. Simultaneous observations of Jupiter with ALMA and Fermi-LAT could be used to provide broad constraints on the energies of the initiating cosmic rays.

Mineral clouds in substellar atmospheres play a special role as a catalyst for a variety of charge processes. If clouds are charged, the surrounding environment becomes electrically activated, and ensembles of charged grains are electrically discharging (e.g. by lightning), which significantly infuences the local chemistry creating conditions similar to those thought responsible for life in early planetary atmospheres. We note that such lightning discharges contribute also to the ionisation state of the atmosphere. We apply scaling laws for electrical discharge processes from laboratory measurements and numerical experiments to Drift-Phoenix model atmosphere results to model the discharge's propagation downwards (as lightning) and upwards (as sprites) through the atmospheric clouds. We evaluate the spatial extent and energetics of lightning discharges. The atmospheric volume affected (e.g. by increase of temperature or electron number) is larger in a brown dwarf atmosphere ($10^8 -~10^{10}$m$^3$) than in a giant gas planet's ($10^4 -~10^{6}$m$^3$). Our results suggest that the total dissipated energy in one event is $<10^{12}$ J for all models of initial solar metallicity. First attempts to show the infuence of lightning on the local gas phase indicate an increase of small carbohydrate molecules like CH and CH$_2$ at the expense of CO and CH$_4$. Dust forming molecules are destroyed and the cloud particle properties are frozen-in unless enough time is available for complete evaporation. We summarise instruments potentially suitable to observe lightning on extrasolar objects.

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.

Charged dust grains in the atmospheres of exoplanets may play a key role in the formation of prebiotic molecules, necessary to the origin of life. Dust grains submerged in an atmospheric plasma become negatively charged and attract a flux of ions that are accelerated from the plasma. The energy of the ions upon reaching the grain surface may be sufficient to overcome the activation energy of particular chemical reactions that would be unattainable via ion and neutral bombardment from classical, thermal excitation. As a result, prebiotic molecules or their precursors could be synthesised on the surface of dust grains that form clouds in exoplanetary atmospheres. This paper investigates the energization of the plasma ions, and the dependence on the plasma electron temperature, in the atmospheres of substellar objects such as gas giant planets. Calculations show that modest electron temperatures of $\approx 1$ eV ($\approx 10^{4}$ K) are enough to accelerate ions to sufficient energies that exceed the activation energies required for the formation of formaldehyde, ammonia, hydrogen cyanide and the amino acid glycine.

Observations of continuous radio and sporadic X-ray emission from low-mass objects suggest they harbour localized plasmas in their atmospheric environments. For low-mass objects, the degree of thermal ionization is insufficient to qualify the ionized component as a plasma, posing the question: what ionization processes can efficiently produce the required plasma that is the source of the radiation? We propose Alfv\'{e}n ionization as a mechanism for producing localized pockets of ionized gas in the atmosphere, having sufficient degrees of ionization ($\geq10^{-7}$) that they constitute plasmas. We outline the criteria required for Alfv\'{e}n ionization and demonstrate it's applicability in the atmospheres of low-mass objects such as giant gas planets, brown dwarfs and M-dwarfs for both solar and sub-solar metallicities. We find that Alfv\'{e}n ionization is most efficient at mid to low atmospheric pressures where a seed plasma is easier to magnetize and the pressure gradients needed to drive the required neutral flows are the smallest. For the model atmospheres considered, our results show that degrees of ionization of $10^{-6}-1$ can be obtained as a result of Alfv\'{e}n ionization. Observable consequences include continuum Bremsstrahlung emission, superimposed with spectral lines from the plasma ion species (e.g. He, Mg, H$_{2}$ or CO lines). Forbidden lines are also expected from the metastable population. The presence of an atmospheric plasma opens the door to a multitude of plasma and chemical processes not yet considered in current atmospheric models. The occurrence of Alfv\'{e}n ionization may also be applicable to other astrophysical environments such as protoplanetary disks.

Being able to accurately estimate stellar parameters based on spectral observations is important not only for understanding the stars themselves but it is also vital for the determination of exoplanet parameters. M dwarfs are discussed as targets for planet detection as these stars are less massive, less luminous and have smaller radii making it possible to detect smaller and lighter planets. Therefore M-dwarfs could prove to be a valuable source for examining the lower mass end of planet distribution, but in order to do that, one must first take care to understand the characteristics of the host stars well enough. Up to date, there are several families of stellar model atmospheres. We focus on the ATLAS9, MARCS and Drift-Phoenix families in the M-dwarf parameter space. We examine the differences in the (Tgas, pgas) structures, synthetic photometric fluxes and related colour indices.We find discrepancies in the hotter regions of the stellar atmosphere between the ATLAS and MARCS models. The MARCS and Drift-Phoenix models appear to agree to a better extend with variances of less than 300K. We have compiled the broad-band synthetic photometric fluxes of all models for the Johnson UBVRI and 2MASS JHKs. The fluxes of MARCS differ from both ATLAS and Drift-Phoenix models in the optical range.

Electric discharges were detected directly in the cloudy atmospheres of Earth, Jupiter and Saturn, are debatable for Venus, and indirectly inferred for Neptune and Uranus in our solar system. Sprites (and other types of transient luminous events) have been detected only on Earth, and are theoretically predicted for Jupiter, Saturn and Venus. Cloud formation is a common phenomenon in ultra-cool atmospheres such as in Brown Dwarf and extrasolar planetary atmospheres. Cloud particles can be expected to carry considerable charges which may trigger discharge events via small-scale processes between individual cloud particles (intra-cloud discharges) or large-scale processes between clouds (inter-cloud discharges). We investigate electrostatic breakdown characteristics, like critical field strengths and critical charge densities per surface, to demonstrate under which conditions mineral clouds undergo electric discharge events which may trigger or be responsible for sporadic X-ray emission. We apply results from our kinetic dust cloud formation model that is part of the Drift-Phoenix model atmosphere simulations. We present a first investigation of the dependence of the breakdown conditions in Brown Dwarf and giant gas exoplanets on the local gas-phase chemistry, the effective temperature and primordial gas-phase metallicity. Our results suggest that different intra-cloud discharge processes dominate at different heights inside mineral clouds: local coronal (point discharges) and small-scale sparks at the bottom region of the cloud where the gas density is high, and flow discharges and large-scale sparks near, and maybe above, the cloud top. The comparison of the thermal degree of ionisation and the number density of cloud particles allows us to suggest the efficiency with which discharges will occur in planetary atmospheres.

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.

It has long been suspected that metal polluted white dwarfs (types DAZ, DBZ, and DZ) and white dwarfs with dusty disks possess planetary systems, but a specific physical mechanism by which planetesimals are perturbed close to a white dwarf has not yet been fully posited. In this paper we demonstrate that mass loss from a central star during post main sequence evolution can sweep planetesimals into interior mean motion resonances with a single giant planet. These planetesimals are slowly removed through chaotic excursions of eccentricity that in time create radial orbits capable of tidally disrupting the planetesimal. Numerical N-body simulations of the Solar System show that a sufficient number of planetesimals are perturbed to explain white dwarfs with both dust and metal pollution, provided other white dwarfs have more massive relic asteroid belts. Our scenario requires only one Jupiter-sized planet and a sufficient number of asteroids near its 2:1 interior mean motion resonance. Finally, we show that once a planetesimal is perturbed into a tidal crossing orbit, it will become disrupted after the first pass of the white dwarf, where a highly eccentric stream of debris forms the main reservoir for dust producing collisions. These simulations, in concert with observations of white dwarfs, place interesting limits on the frequency of planetary systems around main sequence stars, the frequency of planetesimal belts, and the probability that dust may obscure future terrestrial planet finding missions.

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.

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.

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.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

Earth-sized planets around nearby stars are being detected for the first time by ground-based radial velocity and space-based transit surveys. This milestone is opening the path towards the definition of missions able to directly detect the light from these planets, with the identification of bio-signatures as one of the main objectives. In that respect, both ESA and NASA have identified nulling interferometry as one of the most promising techniques. The ability to study distant planets will however depend on exozodiacal dust clouds surrounding the target stars. In this paper, we assess the impact of exozodiacal dust clouds on the performance of an infrared nulling interferometer in the Emma X-array configuration. For the nominal mission architecture with 2-m aperture telescopes, we found that point-symmetric exozodiacal dust discs about 100 times denser than the solar zodiacal cloud can be tolerated in order to survey at least 150 targets during the mission lifetime. Considering modeled resonant structures created by an Earth-like planet orbiting at 1 AU around a Sun-like star, we show that the tolerable dust density for planet detection goes down to about 15 times the solar zodiacal density for face-on systems and decreases with the disc inclination. The upper limits on the tolerable exozodiacal dust density derived in this study must be considered as rather pessimistic, but still give a realistic estimation of the typical sensitivity that we will need to reach on exozodiacal discs in order to prepare the scientific programme of future Earth-like planet characterisation missions.

We present a novel description of how energetic electrons may be ejected from the pulsar interior into the atmosphere, based on the collective electrostatic oscillations of interior electrons confined to move parallel to the magnetic field. The size of the interior magnetic field influences the interior plasma frequency, via the associated matter density compression. The plasma oscillations occur close to the regions of maximum magnetic field curvature, that is, close to the magnetic poles where the majority of magnetic flux emerges. Given that these oscillations have a density-dependent maximum amplitude before wave-breaking occurs, such waves can eject energetic electrons using only the self-field of the electron population in the interior. Moreover, photons emitted by electrons in the bulk of the oscillation can escape along the field lines by virtue of the lower opacity there (and the fact that they are emitted predominantly in this direction), leading to features in the spectra of pulsars.

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.

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.

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.

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}}/\beta$, where $a_{\rm p}$ is the planet's semi-major axis and $\beta$ 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.

The inner Solar System contains a cloud of small (1-100 micron) dust grains created when small bodies-asteroids, comets, and Kuiper belt objects-collide and outgas. This dust cloud, the zodiacal cloud probably has extrasolar analogs, exozodiacal clouds. Exozodiacal clouds are related to debris disks, clouds of rocks and dust orbiting main sequence stars thought to represent the debris left over from planet formation. Some debris disks appear to contain distinct inner clouds that could be considered massive exozodiacal clouds (e.g. Koerner et al. 1998, Absil et al. 2006). This white paper addresses the need for future theoretical work on the dynamics of exozodiacal clouds. This theoretical work should help us discover new planets and understand exozodiacal clouds as astrophysical noise. So far, observations of nearby stars have not provided good constraints on the structures of exozodiacal clouds. But future observations probably will demand a better theoretical understanding of these systems.

Electrostatic oscillations in cold electron-positron plasmas can be coupled to a propagating electromagnetic mode if the background magnetic field is inhomogeneous. Previous work considered this coupling in the quasi-linear regime, successfully simulating the electromagnetic mode. Here we present a stability analysis of the non-linear problem, perturbed from dynamical equilibrium, in order to gain some insight into the modes present in the system.