Abstract
High-contrast imaging has afforded astronomers the opportunity to study light directly emitted by adolescent (tens of megayears) and “proto” (<10 Myr) planets still undergoing formation. ...Direct detection of these planets is enabled by empirical point-spread function (PSF) modeling and removal algorithms. The computational intensity of such algorithms, as well as their multiplicity of tunable input parameters, has led to the prevalence of ad hoc optimization approaches to high-contrast imaging results. In this work, we present a new, systematic approach to optimization vetted using data of the high-contrast stellar companion HD 142527 B from the Magellan Adaptive Optics Giant Accreting Protoplanet Survey (GAPlanetS). More specifically, we present a grid search technique designed to explore three influential parameters of the PSF subtraction algorithm
pyKLIP
: annuli, movement, and KL modes. We consider multiple metrics for postprocessed image quality in order to optimally recover at H
α
(656 nm) synthetic planets injected into contemporaneous continuum (643 nm) images. These metrics include peak (single-pixel) signal-to-noise ratio (S/N), average (multipixel average) S/N, 5
σ
contrast, and false-positive fraction. We apply continuum-optimized KLIP reduction parameters to six H
α
direct detections of the low-mass stellar companion HD 142527 B and recover the companion at a range of separations. Relative to a single-informed, nonoptimized set of KLIP parameters applied to all data sets uniformly, our multimetric grid search optimization led to improvements in companion S/N of up to 1.2
σ
, with an average improvement of 0.6
σ
. Since many direct imaging detections lie close to the canonical 5
σ
threshold, even such modest improvements may result in higher yields in future imaging surveys.
High-contrast imaging observations are fundamentally limited by the spatially and temporally correlated noise source called speckles. Suppression of speckle noise is the key goal of wavefront control ...and adaptive optics (AO), coronagraphy, and a host of post-processing techniques. Speckles average at a rate set by the statistical speckle lifetime, and speckle-limited integration time in long exposures is directly proportional to this lifetime. As progress continues in post-coronagraph wavefront control, residual atmospheric speckles will become the limiting noise source in high-contrast imaging, so a complete understanding of their statistical behavior is crucial to optimizing high-contrast imaging instruments. Here we present a novel power spectral density method for calculating the lifetime, and develop a semi-analytic method for predicting intensity PSDs behind a coronagraph. Considering a frozen-flow turbulence model, we analyze the residual atmosphere speckle lifetimes in a MagAO-X-like AO system as well as 25–39 m giant segmented mirror telescope scale systems. We find that standard AO control shortens atmospheric speckle lifetime from ∼130 ms to ∼50 ms, and predictive control will further shorten the lifetime to ∼20 ms on 6.5 m MagAO-X. We find that speckle lifetimes vary with diameter, wind speed, seeing, and location within the AO control region. On bright stars lifetimes remain within a rough range of ∼20 ms to ∼100 ms. Due to control system dynamics there are no simple scaling laws which apply across a wide range of system characteristics. Finally, we use these results to argue that telemetry-based post-processing should enable ground-based telescopes to achieve the photon-noise limit in high-contrast imaging.
Abstract
We present resolved images of the inner disk component around HD 141569A using the Magellan adaptive optics system with the Clio2 1–5
μ
m camera, offering a glimpse of a complex system ...thought to be in a short evolutionary phase between protoplanetary and debris disk stages. We use a reference star along with the Karhunen–Loéve image projection (KLIP) algorithm for point-spread function subtraction to detect the disk inward to about 0.″24 (∼25 au assuming a distance of 111 pc) at high signal-to-noise ratios at
L
′
(3.8
μ
m),
Ls
(3.3
μ
m), and narrowband
Ice
(3.1
μ
m). We identify an arc or spiral arm structure at the southeast extremity, consistent with previous studies. We implement forward modeling with a simple disk model within the framework of a Markov Chain Monte Carlo sampler to better constrain the geometrical attributes and photometry using our KLIP-reduced disk images. We then leverage these modeling results to facilitate a comparison of the measured brightness in each passband to find a reduction in scattered light from the disk in the
Ice
filter, implying significant absorption due to water ice in the dust. Additionally, our best-fit disk models exhibit peak brightness in the southwestern, back-scattering region of the disk, which we suggest to be possible evidence of 3.3
μ
m polycyclic aromatic hydrocarbon emission. However, we point out the need for additional observations with bluer filters and more complex modeling to confirm these hypotheses.
Abstract
Starlight subtraction algorithms based on the method of Karhunen–Loève eigenimages have proved invaluable to exoplanet direct imaging. However, they scale poorly in runtime when paired with ...differential imaging techniques. In such observations, reference frames and frames from which starlight is to be subtracted are drawn from the same set of data, requiring a new subset of references (and eigenimages) for each frame processed to avoid self-subtraction of the signal of interest. The data rates of extreme adaptive optics instruments are such that the only way to make this computationally feasible has been to downsample the data. We develop a technique that updates a precomputed singular value decomposition of the full data set to remove frames (i.e., a “downdate”) without a full recomputation, yielding the modified eigenimages. This not only enables analysis of much larger data volumes in the same amount of time, but also exhibits near-linear scaling in runtime as the number of observations increases. We apply this technique to archival data and investigate its scaling behavior for very large numbers of frames
N
. The resulting algorithm provides speed improvements of 2.6× (for 200 eigenimages at
N
= 300) to 140× (at
N
= 10
4
) with the advantage only increasing as
N
grows. This algorithm has allowed us to substantially accelerate Karhunen–Loève image projection (KLIP) even for modest
N
, and will let us quickly explore how KLIP parameters affect exoplanet characterization in large-
N
data sets.
Abstract
A lower bound on unbiased estimates of wave front errors (WFEs) is presented for the linear regime of small perturbation and active control of a high-contrast region (dark hole). Analytical ...approximations and algorithms for computing the closed-loop covariance of the WFE modes are provided for discrete- and continuous-time linear WFE dynamics. Our analysis applies to both image-plane and non-common-path wave front sensing (WFS) with Poisson-distributed measurements and noise sources (i.e., photon-counting mode). Under this assumption, we show that recursive estimation benefits from infinitesimally short exposure times, is more accurate than batch estimation and, for high-order WFE drift dynamical processes, scales better than batch estimation with amplitude and star brightness. These newly derived contrast scaling laws are a generalization of previously known theoretical and numerical results for turbulence-driven adaptive optics. For space-based coronagraphs, we propose a scheme for combining models of WFE drift, low-order non-common-path WFS (LOWFS) and high-order image-plane WFS (HOWFS) into closed-loop contrast estimates. We also analyze the impact of residual low-order WFE, sensor noise, and other sources incoherent with the star, on closed-loop dark hole maintenance and the resulting contrast. As an application example, our model suggests that the Roman Space Telescope might operate in a regime that is dominated by incoherent sources rather than WFE drift, where the WFE drift can be actively rejected throughout the observations with residuals significantly dimmer than the incoherent sources. The models proposed in this paper make possible the assessment of the closed-loop contrast of coronagraphs with combined LOWFS and HOWFS capabilities, and thus help estimate WFE stability requirements of future instruments.
Abstract
Multi-wavelength observations provide a complementary view of the formation of young, directly imaged planet-mass companions. We report the ALMA 1.3 mm and
Magellan
adaptive optics H
α
,
,
, ...and
Y
S
observations of the GQ Lup system, a classical T Tauri star with a
substellar companion at ∼110 au projected separation. We estimate the accretion rates for both components from the observed H
α
fluxes. In our ∼0.″05 resolution ALMA map, we resolve GQ Lup A’s disk in the dust continuum, but no signal is found from the companion. The disk is compact, with a radius of ∼22 au, a dust mass of ∼6
M
⊕
, an inclination angle of ∼56°, and a very flat surface density profile indicative of a radial variation in dust grain sizes. No gaps or inner cavity are found in the disk, so there is unlikely a massive inner companion to scatter GQ Lup B outward. Thus, GQ Lup B might have formed in situ via disk fragmentation or prestellar core collapse. We also show that GQ Lup A’s disk is misaligned with its spin axis, and possibly with GQ Lup B’s orbit. Our analysis on the tidal truncation radius of GQ Lup A’s disk suggests that GQ Lup B’s orbit might have a low eccentricity.
Abstract
Companions embedded in the cavities of transitional circumstellar disks have been observed to exhibit excess luminosity at H
α
, an indication that they are actively accreting. We report 5 ...yr (2013–2018) of monitoring of the position and H
α
excess luminosity of the embedded, accreting low-mass stellar companion HD 142527 B from the MagAO/VisAO instrument. We use
pyklip
, a Python implementation of the Karhunen–Loeve Image Processing algorithm, to detect the companion. Using
pyklip
forward modeling, we constrain the relative astrometry to 1–2 mas precision and achieve sufficient photometric precision (±0.2 mag, 3% error) to detect changes in the H
α
contrast of the companion over time. In order to accurately determine the relative astrometry of the companion, we conduct an astrometric calibration of the MagAO/VisAO camera against 20 yr of Keck/NIRC2 images of the Trapezium cluster. We demonstrate agreement of our VisAO astrometry with other published positions for HD 142527 B, and use
orbitize!
to generate a posterior distribution of orbits fit to the relative astrometry of HD 142527 B. Our data suggest that the companion is close to periastron passage, on an orbit significantly misaligned with respect to both the wide circumbinary disk and the recently observed inner disk encircling HD 142527 A. We translate observed H
α
contrasts for HD 142527 B into mass accretion rate estimates on the order of 4–9 × 10
−10
M
⊙
yr
−1
. Photometric variation in the H
α
excess of the companion suggests that the accretion rate onto the companion is variable. This work represents a significant step toward observing accretion-driven variability onto protoplanets, such as PDS 70 b&c.
Abstract
Accreting protoplanets are windows into planet formation processes, and high-contrast differential imaging is an effective way to identify them. We report results from the Giant Accreting ...Protoplanet Survey (GAPlanetS), which collected H
α
differential imagery of 14 transitional disk host stars with the Magellan Adaptive Optics System. To address the twin challenges of morphological complexity and point-spread function instability, GAPlanetS required novel approaches for frame selection and optimization of the Karhounen–Loéve Image Processing algorithm
pyKLIP
. We detect one new candidate, CS Cha “c,” at a separation of 68 mas and a modest Δmag of 2.3. We recover the HD 142527 B and HD 100453 B accreting stellar companions in several epochs, and the protoplanet PDS 70 c in 2017 imagery, extending its astrometric record by nine months. Though we cannot rule out scattered light structure, we also recover LkCa 15 “b,” at H
α
; its presence inside the disk cavity, absence in Continuum imagery, and consistency with a forward-modeled point source suggest that it remains a viable protoplanet candidate. Through targeted optimization, we tentatively recover PDS 70 c at two additional epochs and PDS 70 b in one epoch. Despite numerous previously reported companion candidates around GAplanetS targets, we recover no additional point sources. Our moderate H
α
contrasts do not preclude most protoplanets, and we report limiting H
α
contrasts at unrecovered candidate locations. We find an overall detection rate of ∼36
−
22
+
26
%
, considerably higher than most direct imaging surveys, speaking to both GAPlanetS’s highly targeted nature and the promise of H
α
differential imaging for protoplanet identification.