Direct imaging of exoplanets requires establishing and maintaining a high-contrast dark field (DF) within the science image to a high degree of precision (10(-10)). Current approaches aimed at establishing the DF, such as electric field conjugation (EFC), have been demonstrated in the lab and have proven capable of high-contrast DF generation. The same approaches have been considered for the maintenance of the DF as well. However, these methods rely on phase diversity measurements, which require field modulation; this interrupts the DF and consequently competes with the science acquisition. We introduce and demonstrate spatial linear dark field control (LDFC) as an alternative technique by which the high-contrast DF can be maintained without modulation. Once the DF has been established by conventional EFC, spatial LDFC locks the high-contrast state of the DF by operating a closed loop around the linear response of the bright field (BF) to wavefront variations that modify both the BF and the DF. We describe the fundamental operating principles of spatial LDFC and provide numerical simulations of its operation as a DF stabilization technique that is capable of wavefront correction within the DF without interrupting science acquisition. (c) The Authors.

We review astronomical results in the visible (lambda<1 mu m) with adaptive optics. Other than a brief period in the early 1990s, there has been little (<1 paper/yr) night-time astronomical science published with AO in the visible from 2000-2013 (outside of the solar or Space Surveillance Astronomy communities where visible AO is the norm, but not the topic of this invited review). However, since mid-2013 there has been a rapid increase visible AO with over 50 refereed science papers published in just similar to 2.5 years (visible AO is experiencing a rapid growth rate very similar to that of NIR AO science from 1997-2000; Close 2000). Currently the most productive small (D < 2 m) visible light AO telescope is the UV-LGS Robo-AO system (Baranec, et al. 2016) on the robotic Palomar D=1.5 m telescope (currently relocated to the Kitt Peak 1.8m; Salama et al. 2016). Robo-AO uniquely offers the ability to target >15 objects/hr, which has enabled large (>3000 discrete targets) companion star surveys and has resulted in 23 refereed science publications. The most productive large telescope visible AO system is the D=6.5m Magellan telescope AO system (MagAO). MagAO is an advanced Adaptive Secondary Mirror (ASM) AO system at the Magellan 6.5m in Chile (Morzinski et al. 2016). This ASM secondary has 585 actuators with < 1 msec response times (0.7 ms typically). MagAO utilizes a 1 kHz pyramid wavefront sensor. The relatively small actuator pitch (similar to 22 cm/subap) allows moderate Strehls to be obtained in the visible (0.63-1.05 microns). Long exposures (60s) achieve <30mas resolutions, 30% Strehls at 0.62 microns (r') with the VisAO camera in 0.5" seeing with bright R <= 9 mag stars. These capabilities have led to over 22 MagAO refereed science publications in the visible. The largest (D=8m) telescope to achieve regular visible AO science is SPHERE/ZIMPOL. ZIMPOL is a polarimeter fed by the similar to 1.2 kHz SPHERE ExAO system (Fusco et al. 2016). ZIMPOL's ability to differentiate scattered polarized light from starlight allows the sensitive detection of circumstellar disks, stellar surfaces, and envelopes of evolved AGB stars. Here we review the key steps to having good performance in the visible and review the exciting new AO visible science opportunities and science results in the fields of: exoplanet detection; circumstellar and protoplanetary disks; young stars; AGB stars; emission line jets; and stellar surfaces. The recent rapid increase in the scientific publications and power of visible AO is due to the maturity of the next-generation of AO systems and our new ability probe circumstellar regions with very high (10-30 mas) spatial resolutions that would otherwise require much larger (> 10m) diameter telescopes in the infrared.

We present measurements of the evolutionary timescales of speckles around adaptive optics-corrected PSFs. We placed a SELEX SAPHIRA HgCdTe detector behind the SCExA0 instrument at Subaru Telescope. We analyzed the behavior of speckles at radial distances of 2-8 A/D away from the diffraction-limited PSF in H-band (-1.6 m) images collected at 1 kHz framerates. Speckles evolve with a variety of timescales, and these have not previously been studied at near-infrared wavelengths. Ultimately we would like to image reflected-light exoplanets, which necessitates a fast speckle control loop. Based on our measurements, we calculate the parameters of an optimized control loop that would enable such observations.

The Keck Planet Imager and Characterizer (KPIC) is a cost-effective upgrade path to the W.M. Keck observatory (WMKO) adaptive optics (AO) system, building on the lessons learned from first and second-generation extreme AO (ExA0) coronagraphs. KPIC will explore new scientific niches in exoplanet science, while maturing critical technologies and systems for future ground-based (TMT, FELT, GMT) and space-based planet imagers (HabEx, LUVOIR). The advent of fast low-noise IR cameras (IR-APD, MKIDS, electron injectors), the rapid maturing of efficient wavefront sensing (WFS) techniques (Pyramid, Zernike), small inner working angle (IWA) coronagraphs (e.g., vortex) and associated low-order wavefront sensors (LOWFS), as well as recent breakthroughs in high contrast high resolution spectroscopy, open new direct exoplanet exploration avenues that are complementary to planet imagers such as VLT-SPHERE and the Gemini Planet Imager (GPI). For instance, the search and detailed characterization of planetary systems on solar-system scales around late-type stars, mostly beyond SPHERE and GPI's reaches, can be initiated now at WMKO.

The characterization of exozodiacal light emission is both important for the understanding of planetary systems evolution and for the preparation of future space missions aiming to characterize low mass planets in the habitable zone of nearby main sequence stars. The Large Binocular Telescope Interferometer (LBTI) exozodi survey aims at providing a ten-fold improvement over current state of the art, measuring dust emission levels down to a typical accuracy of similar to 12 zodis per star, for a representative ensemble of similar to 30+ high priority targets. Such measurements promise to yield a final accuracy of about 2 zodis on the median exozodi level of the targets sample. Reaching a 1. measurement uncertainty of 12 zodis per star corresponds to measuring interferometric cancellation ("null") levels, i.e visibilities at the few 100 ppm uncertainty level. We discuss here the challenges posed by making such high accuracy mid-infrared visibility measurements from the ground and present the methodology we developed for achieving current best levels of 500 ppm or so. We also discuss current limitations and plans for enhanced exozodi observations over the next few years at LBTI.

Transition disks, protoplanetary disks with inner clearings, are promising objects in which to directly image forming planets. The high contrast imaging technique of non-redundant masking is well posed to detect planetary mass companions at several to tens of AU in nearby transition disks. We present non-redundant masking observations of the T Cha and LkCa 15 transition disks, both of which host posited sub-stellar mass companions. However, due to a loss of information intrinsic to the technique, observations of extended sources (e.g. scattered light from disks) can be misinterpreted as moving companions. We discuss tests to distinguish between these two scenarios, with applications to the T Cha and LkCa 15 observations. We argue that a static, forward-scattering disk can explain the T Cha data, while LkCa 15 is best explained by multiple orbiting companions.

In this dissertation, we develop multiscale finite element methods and uncertainty quantification technique for Richards' equation, a mathematical model to describe fluid flow in unsaturated porous media. Both coarse-level and fine-level numerical computation techniques are presented. To develop an accurate coarse-scale numerical method, we need to construct an effective multiscale map that is able to capture the multiscale features of the large-scale solution without resolving the small scale details. With a careful choice of the coarse spaces for multiscale finite element methods, we can significantly reduce errors. We introduce several methods to construct coarse spaces for multiscale finite element methods. A coarse space based on local spectral problems is also presented. The construction of coarse spaces begins with an initial choice of multiscale basis functions supported in coarse regions. These basis functions are complemented using weighted local spectral eigenfunctions. These newly constructed basis functions can capture the small scale features of the solution within a coarse-grid block and give us an accurate coarse-scale solution. However, it is expensive to compute the local basis functions for each parameter value for a nonlinear equation. To overcome this difficulty, local reduced basis method is discussed, which provides smaller dimension spaces with which to compute the basis functions. Robust solution techniques for Richards' equation at a fine scale are discussed. We construct iterative solvers for Richards' equation, whose number of iterations is independent of the contrast. We employ two-level domain decomposition pre-conditioners to solve linear systems arising in approximation of problems with high contrast. We show that, by using the local spectral coarse space for the preconditioners, the number of iterations for these solvers is independent of the physical properties of the media. Several numerical experiments are given to support the theoretical results. Last, we present numerical methods for uncertainty quantification applications for Richards' equation. Numerical methods combined with stochastic solution techniques are proposed to sample conductivities of porous media given in integrated data. Our proposed algorithm is based on upscaling techniques and the Markov chain Monte Carlo method. Sampling results are presented to prove the efficiency and accuracy of our algorithm.

A dual channel polarimeter was incorporated into the Hokupa'a adaptive optics system mounted on the Gemini North telescope to enhance sensitivity to detecting the light scattered by circumstellar material. The technique suppressed noise introduced by non-repeatable variations of the point spread function which limit the sensitivity of non-simultaneous adaptive optics imaging. Polarimetric images of the classical T-Tauri star environments around GG Tauri Aab, TW Hydrae, LkCa 15, LkHα 242, GM Aurigae, and SR24 N/S were observed to establish the instrument's sensitivity. A survey of nearby ( d < 25 pc), young ( age < 1 Gyr), solar-analog stars was undertaken with the polarimeter to search for collisionally active debris disks analogous to our young solar system. Of the 24 stars sampled, none were found to have obvious scattered light signatures. Isotropic and Mie scattering model images of debris disks were used to constrain the amount of material around the survey stars to no more than M dust ∼ 10 -2 M Moon of 1-10μ m sized dust contained between 5-50 AU from the sample stars. Particle lifetimes under the influence of the Poynting Robertson Drag, radiation pressure, and solar wind drag are calculated as a function of central star spectral type. The corpuscular drag from stellar winds shorten dust lifetimes by an amount inversely proportional to the stellar wind mass-loss rate. This translates into dust lifetimes 100-1000 times shorter around young solar analog stars compared to the present day. This effect, cam significantly reduce the near-IR detectability of debris disks around these chromospherically active stars.

For decades, astronomers have used observations of mature planetary systems to constrain planet formation theories, beginning with our own solar system and now the thousands of known exoplanets. Recent advances in instrumentation have given us a direct view of some steps in the planet formation process, such as large-scale protostar and protoplanetary disk features and evolution. However, understanding the details of how planets accrete and interact with their environment requires direct observations of protoplanets themselves. Transition disks, protoplanetary disks with inner clearings that may be caused by forming planets, are the best targets for these studies. Their large distances, compared to the stars normally targeted for direct imaging of exoplanets, make protoplanet detection difficult and necessitate novel imaging techniques. In this dissertation, I describe the results of using non-redundant masking (NRM) to search for forming planets in transition disk clearings. I first present a data reduction pipeline that I wrote to this end, using example datasets and simulations to demonstrate reduction and imaging optimizations. I discuss two transition disk NRM case studies: T Cha and LkCa 15. In the case of T Cha, while we detect significant asymmetries, the data cannot be explained by orbiting companions. The fluxes and orbital motion of the LkCa 15 companion signals, however, can be naturally explained by protoplanets in the disk clearing. I use these datasets and simulated observations to illustrate the effects of scattered light from transition disk material on NRM protoplanet searches. I then demonstrate the utility of the dual-aperture Large Binocular Telescope Interferometer's NRM mode on the bright B[e] star MWC 349A. I discuss the implications of this work for planet formation studies as well as future prospects for NRM and related techniques on next generation instruments.

On March 2015 an L'-band vortex coronagraph based on an Annular Groove Phase Mask made up of a diamond sub-wavelength grating was installed on NIRC2 as a demonstration project. This vortex coronagraph operates in the L' band not only in order to take advantage from the favorable star/planet contrast ratio when observing beyond the K band, but also to exploit the fact that the Keck II Adaptive Optics (AO) system delivers nearly extreme adaptive optics image quality (Strehl ratios values near 90%) at 3.7 mu m. We describe the hardware installation of the vortex phase mask during a routine NIRC2 service mission. The success of the project depends on extensive software development which has allowed the achievement of exquisite real-time pointing control as well as further contrast improvements by using speckle nulling to mitigate the effect of static speckles. First light of the new coronagraphic mode was on June 2015 with already very good initial results. Subsequent commissioning nights were interlaced with science nights by members of the VORTEX team with their respective scientific programs. The new capability and excellent results so far have motivated the VORTEX team and the Keck Science Steering Committee (KSSC) to offer the new mode in shared risk mode for 2016B.