Search

Estimating the exoplanet scientific productivity of the Habitable Worlds Observatory requires estimating science exposure times. From exoplanet yields to spectral retrievals, exposure times are at the heart of our understanding of the capabilities of this future mission. As such, ensuring accuracy and consistency between different exposure time calculators (ETCs) is critical. We summarize the efforts of the Exoplanet Science Yield sub-Working Group's ETC Calibration Task Group, which conducted a calibration study from March 4 to June 30 of 2024. We compare three commonly-used coronagraphic exposure time calculators. We find that the ETCs use a broad variety of differing methods, assumptions, and inputs that produce variation in the final exposure times at the ~60% level. The causes for the disagreement have largely been identified, flagged for further development efforts, and in some cases retired since the conclusion of this effort. We expect that addressing the flagged efforts will bring the ETCs to within better than ~30% agreement., Comment: 24 pages, 4 figures, 3 tables

A crucial component of the high-contrast instrumental chain in astronomy is the wavefront sensor (WFS). A key property of this component is its sensitivities, which reflect its ability to efficiently use incoming photons to encode the phase aberrations. This paper introduces a new class of highly sensitive wavefront sensors that approach the fundamental sensitivity limits dictated by physics. Assuming a high Strehl regime, we define what linear operator is describing the ideal WFS that would achieve maximum sensitivity. We then show that there is a substantial similarity between this ideal WFS and the second-order ideal coronagraph. Leveraging the exhibited link between ideal wavefront sensing and coronagraphy, we propose a novel WFS concept based on high-performance coronagraphic architecture : the bivortex WFS. This sensor employs charge-2 vortex masks. Simulations for an ideal system demonstrate that this sensor achieves unprecedented sensitivity, even surpassing the highly sensitive Zernike WFS class (especially for low spatial frequencies), while paving the way for new high-contrast architectures integrating simultaneous sensing and coronagraphy.

The Gemini Planet Imager (GPI) is a high contrast imaging instrument that aims to detect and characterize extrasolar planets. GPI is being upgraded to GPI 2.0, with several subsystems receiving a re-design to improve its contrast. To enable observations on fainter targets and increase performance on brighter ones, one of the upgrades is to the adaptive optics system. The current Shack-Hartmann wavefront sensor (WFS) is being replaced by a pyramid WFS with an low-noise electron multiplying CCD (EMCCD). EMCCDs are detectors capable of counting single photon events at high speed and high sensitivity. In this work, we characterize the performance of the HN\"u 240 EMCCD from N\"uv\"u Cameras, which was custom-built for GPI 2.0. Through our performance evaluation we found that the operating mode of the camera had to be changed from inverted-mode (IMO) to non-inverted mode (NIMO) in order to improve charge diffusion features found in the detector's images. Here, we characterize the EMCCD's noise contributors (readout noise, clock-induced charges, dark current) and linearity tests (EM gain, exposure time) before and after the switch to NIMO., Comment: Proceeding of the SPIE Astronomical Telescopes+Instrumentation. 14 pages, 15 figures

The Earth's atmosphere is comprised of turbulent layers that result in speckled and blurry images from ground-based visible and infrared observations. Adaptive Optics (AO) systems are employed to measure the perturbed wavefront with a wavefront sensor (WFS) and correct for these distortions with a deformable mirror. Therefore, understanding and characterising the atmosphere is crucial for the design and functionality of AO systems. One parameter for characterizing the atmosphere is the atmospheric coherence time, which is a function of the effective wind velocity of the atmosphere. This parameter dictates how fast the AO system needs to correct for the atmosphere. If not fast enough, phenomena such as the wind butterfly effect can occur, hindering high-contrast coronographic imaging. This effect is a result of fast, strong, high-altitude turbulent layers. This paper presents two methods for estimating the effective wind velocity, using pseudo-open loop WFS slopes. The first method uses a spatial-temporal covariance map and the second uses the power spectral density of the defocus term. We show both simulated results and preliminary results from the Gemini Planet Imager AO telemetry., Comment: 7 pages, 7 figures, submitted to SPIE Astronomical + Telescopes 2024

The Zernike wavefront sensor (ZWFS) stands out as one of the most sensitive optical systems for measuring the phase of an incoming wavefront, reaching photon efficiencies close to the fundamental limit. This quality, combined with the fact that it can easily measure phase discontinuities, has led to its widespread adoption in various wavefront control applications, both on the ground but also for future space-based instruments. Despite its advantages, the ZWFS faces a significant challenge due to its extremely limited dynamic range, making it particularly challenging for ground-based operations. To address this limitation, one approach is to use the ZWFS after a general adaptive optics (AO) system; however, even in this scenario, the dynamic range remains a concern. This paper investigates two optical configurations of the ZWFS: the conventional setup and its phase-shifted counterpart, which generates two distinct images of the telescope pupil. We assess the performance of various reconstruction techniques for both configurations, spanning from traditional linear reconstructors to gradient-descent-based methods. The evaluation encompasses simulations and experimental tests conducted on the Santa cruz Extreme Adaptive optics Lab (SEAL) bench at UCSC. Our findings demonstrate that certain innovative reconstruction techniques introduced in this study significantly enhance the dynamic range of the ZWFS, particularly when utilizing the phase-shifted version.

The Santa Cruz Extreme AO Lab (SEAL) testbed is an optical bench meant to design and develop new wavefront control techniques for high-contrast imaging for segmented telescopes. These techniques allow for astronomical efficiency in exoplanet imaging and characterization. SEAL consists of several wavefront sensors (WFS) and deformable mirrors (DM) that are currently performing techniques like predictive control or non-linear reconstruction. In this paper, we present the implementation and characterization of a new coronagraphic branch on SEAL and assess the contrast limitations in the testbed. For our coronagraphic branch, we used a vector vortex coronagraph which has high contrast performance. The W. M. Keck Observatory also uses a vortex coronagraph, allowing us to compare the limitations with our own coronagraph. We relied on the testbed and simulations of the vortex coronagraph to compare performance with expected ones. To create a more reliable simulation, we also injected in our numerical model data collected by a Zernike Wavefront sensor (ZWFS) used to perform fine wavefront sensing on the bench. Now that the coronagraphic branch is aligned on SEAL, we will be able to use contrast as a metric for the performance of wavefront control methods on the bench., Comment: 10 pages, 14 figures, SPIE Astronomical Telescopes + Instrumentation: Adaptive Optics Systems IX (Yokohama 2024), Paper 13097-297

The 16 Myr-old A0V star HD 156623 in the Scorpius--Centaurus association hosts a high-fractional-luminosity debris disk, recently resolved in scattered light for the first time by the Gemini Planet Imager (GPI) in polarized intensity. We present new analysis of the GPI H-band polarimetric detection of the HD 156623 debris disk, with particular interest in its unique morphology. This debris disk lacks a visible inner clearing, unlike the majority of low-inclination disks in the GPI sample and in Sco-Cen, and it is known to contain CO gas, positioning it as a candidate ``hybrid'' or ``shielded'' disk. We use radiative transfer models to constrain the geometric parameters of the disk based on scattered light data and thermal models to constrain the unresolved inner radius based on the system's spectral energy distribution (SED). We also compute a measurement of the polarized scattering phase function, adding to the existing sample of empirical phase function measurements. We find that HD 156623's debris disk inner radius is constrained to less than 26.6 AU from scattered light imagery and less than 13.4 AU from SED modeling at a 99.7% confidence interval, and suggest that gas drag may play a role in retaining sub-blowout size dust grains so close to the star., Comment: 18 pages, 10 figures. Accepted for publication in AJ

We summarize the current best polychromatic (10 to 20 % bandwidth) contrast performance demonstrated in the laboratory by different starlight suppression approaches and systems designed to directly characterize exoplanets around nearby stars. We present results obtained by internal coronagraph and external starshade experimental testbeds using entrance apertures equivalent to off-axis or on-axis telescopes, either monolithic or segmented. For a given angular separation and spectral bandwidth, the performance of each starlight suppression system is characterized by the values of raw contrast (before image processing), off-axis (exoplanet) core throughput, and post-calibration contrast (the final 1 sigma detection limit of off-axis point sources, after image processing). To place the current laboratory results in the perspective of the future Habitable Worlds Observatory (HWO) mission, we simulate visible observations of a fiducial Earth/Sun twin system at 12 pc, assuming a 6m (inscribed diameter) collecting aperture and a realistic end-to-end optical throughput. The exposure times required for broadband exoearth detection (20% bandwidth around a wavelength of 0.55 microns) and visible spectroscopic observations (R=70) are then computed assuming various levels of starlight suppression performance, including the values currently demonstrated in the laboratory. Using spectroscopic exposure time as a simple metric, our results point to key starlight suppression system design performance improvements and trades to be conducted in support of HWO exoplanet science capabilities. These trades may be explored via numerical studies, lab experiments, as well as high contrast space-based observations and demonstrations., Comment: 63 pages, 28 pages, submitted to JATIS

The Atmospheric Dispersion Corrector (ADC) of the Gemini Planet Imager (GPI) corrects the chromatic dispersion caused by differential atmospheric refraction (DAR), making it an important optic for exoplanet observation. Despite requiring less than 5 mas of residual DAR to avoid potentially affecting the coronagraph, the GPI ADC averages $\sim7$ and $\sim11$ mas of residual DAR in $H$ and $J$ band respectively. We analyzed GPI data in those bands to find explanations for the underperformance. We found the model GPI uses to predict DAR underestimates humidity's impact on incident DAR, causing on average a 0.54 mas increase in $H$ band residual DAR. Additionally, the GPI ADC consistently undercorrects in $H$ band by about 7 mas, causing almost all the $H$ band residual DAR. $J$ band does not have such an offset. Perpendicular dispersion induced by the GPI ADC, potentially from a misalignment in the prisms' relative orientation, causes 86% of the residual DAR in $J$ band. Correcting these issues could reduce residual DAR, thereby improving exoplanet detection. We also made a new approximation for the index of refraction of air from 0.7 microns to 1.36 microns that more accurately accounts for the effects of humidity., Comment: 21 pages, 10 figures, published in JATIS

The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models., Comment: 32 pages, 16 figures, 6 tables, 2 appendices

We present new spatially resolved astrometry and photometry of the CD-27 11535 system, a member of the $\beta$ Pictoris moving group consisting of two resolved K-type stars on a $\sim$20-year orbit. We fit an orbit to relative astrometry measured from NIRC2, GPI, and archival NaCo images, in addition to literature measurements. However, the total mass inferred from this orbit is significantly discrepant from that inferred from stellar evolutionary models using the luminosity of the two stars. We explore two hypotheses that could explain this discrepant mass sum; a discrepant parallax measurement from Gaia due to variability, and the presence of an additional unresolved companion to one of the two components. We find that the $\sim$20-year orbit could not bias the parallax measurement, but that variability of the components could produce a large amplitude astrometric motion, an effect which cannot be quantified exactly without the individual Gaia measurements. The discrepancy could also be explained by an additional star in the system. We jointly fit the astrometric and photometric measurements of the system to test different binary and triple architectures for the system. Depending on the set of evolutionary models used, we find an improved goodness of fit for a triple system architecture that includes a low-mass ($M=0.177\pm0.055$\,$M_{\odot}$) companion to the primary star. Further studies of this system will be required in order to resolve this discrepancy, either by refining the parallax measurement with a more complex treatment of variability-induced astrometric motion, or by detecting a third companion., Comment: 12 pages + references and appendix, 12 figures, 5 tables

The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet-disk interactions. The multi-wavelength surface brightness, the disk color and geometry permit identification of any asymmetries such as warps or disk offsets from the central star. We find that nineteen of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system are discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power-law, $q$. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications on the evolution of debris disk systems around stars of various spectral types., Comment: 46 pages, 20 figures, 6 tables, accepted for publication in ApJ

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric techniques such as sparse aperture masking. We introduce a nascent Python package, NRM-artist, that we use to design several SCALES masks to be non-redundant and to have uniform coverage in Fourier space. We generate high-fidelity mock SCALES data using the scalessim package for SCALES' low spectral resolution modes across its 2 to 5 micron bandpass. We include realistic noise from astrophysical and instrument sources, including Keck adaptive optics and Poisson noise. We inject planet and disk signals into the mock datasets and subsequently recover them to test the performance of SCALES sparse aperture masking and to determine the sensitivity of various mask designs to different science signals.

SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution. We explore the sensitivity limitations of this mode by simulating planet detections in the presence of realistic noise sources. We use the SCALES simulator scalessim to generate high-fidelity mock observations of planets that include speckle noise from their host stars, as well as other atmospheric and instrumental noise effects. We employ both angular and reference differential imaging as methods of disentangling speckle noise from the injected planet signals. These simulations allow us to assess the feasibility of speckle deconvolution for SCALES medium resolution data, and to test whether one approach outperforms another based on planet angular separations and contrasts.

We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{\mu m}$ as a part of the JWST Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of $0.5\lambda/D$ for an interferometer), which are inaccessible with the classical inner working angles of the JWST coronagraphs. When combined with JWST's unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a $5\sigma$ contrast of $\Delta m{\sim}7.62{\pm}0.13$ mag relative to the host star at separations ${\gtrsim}0.07{"}$, and the contrast deteriorates steeply at separations ${\lesssim}0.07{"}$. However, we detect no additional companions interior to the known companion HIP 65426 b (at separation ${\sim}0.82{"}$ or, $87^{+108}_{-31}\,\rm{au}$). Our observations thus rule out companions more massive than $10{-}12\,\rm{M_{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP 65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on JWST is sensitive to planetary mass companions at close-in separations (${\gtrsim}0.07{"}$), even for thousands of more distant stars at $\sim$100 pc, in addition to the stars in the nearby young moving groups as stated in previous works. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening an essentially unexplored parameter space., Comment: 23 pages, 14 figures, accepted for publication in ApJL

We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim \lambda/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy., Comment: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~\mu$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~\mu$m) high-contrast spectrographs. This will lead to new progress in exoplanet atmospheric studies, including detailed characterization of benchmark systems that will advance the state of the art of atmospheric modeling. SCALES' unique modes, while designed specifically for direct exoplanet characterization, will enable a broader range of novel (exo)planetary observations as well as galactic and extragalactic studies. Here we present the science cases that drive the design of SCALES. We describe an end-to-end instrument simulator that we use to track requirements, and show simulations of expected science yields for each driving science case. We conclude with a discussion of preparations for early science when the instrument sees first light in $\sim2025$., Comment: 10 pages, 16 figures, submitted to Proceedings of the SPIE

The Gemini Planet Imager (GPI) is a high contrast imaging instrument that aims to detect and characterize extrasolar planets. GPI is being upgraded to GPI 2.0, with several subsystems receiving a re-design to improve the instrument's contrast. To enable observations on fainter targets and increase stability on brighter ones, one of the upgrades is to the adaptive optics system. The current Shack-Hartmann wavefront sensor (WFS) is being replaced by a pyramid WFS with an low-noise electron multiplying CCD (EMCCD). EMCCDs are detectors capable of counting single photon events at high speed and high sensitivity. In this work, we characterize the performance of the HN\"u 240 EMCCD from N\"uv\"u Cameras, which was custom-built for GPI 2.0. The HN\"u 240 EMCCD's characteristics make it well suited for extreme AO: it has low dark current ($<$ 0.01 e-/pix/fr), low readout noise (0.1 e-/pix/fr at a gain of 5000), high quantum efficiency ( 90% at wavelengths from 600-800 nm; 70% from 800-900 nm), and fast readout (up to 3000 fps full frame). Here we present test results on the EMCCD's noise contributors, such as the readout noise, pixel-to-pixel variability and CCD bias. We also tested the linearity and EM gain calibration of the detector. All camera tests were conducted before its integration into the GPI 2.0 PWFS system., Comment: 16 pages, 14 figures. Conference Proceedings for AO4ELT7, held in June 2023 in Avignon, France

Adaptive optics (AO) is a technique to improve the resolution of ground-based telescopes by correcting, in real-time, optical aberrations due to atmospheric turbulence and the telescope itself. With the rise of Giant Segmented Mirror Telescopes (GSMT), AO is needed more than ever to reach the full potential of these future observatories. One of the main performance drivers of an AO system is the wavefront sensing operation, consisting of measuring the shape of the above mentioned optical aberrations. Aims. The non-modulated pyramid wavefront sensor (nPWFS) is a wavefront sensor with high sensitivity, allowing the limits of AO systems to be pushed. The high sensitivity comes at the expense of its dynamic range, which makes it a highly non-linear sensor. We propose here a novel way to invert nPWFS signals by using the principle of reciprocity of light propagation and the Gerchberg-Saxton (GS) algorithm. We test the performance of this reconstructor in two steps: the technique is first implemented in simulations, where some of its basic properties are studied. Then, the GS reconstructor is tested on the Santa Cruz Extreme Adaptive optics Laboratory (SEAL) testbed located at the University of California Santa Cruz. This new way to invert the nPWFS measurements allows us to drastically increase the dynamic range of the reconstruction for the nPWFS, pushing the dynamics close to a modulated PWFS. The reconstructor is an iterative algorithm requiring heavy computational burden, which could be an issue for real-time purposes in its current implementation. However, this new reconstructor could still be helpful in the case of many wavefront control operations. This reconstruction technique has also been successfully tested on the Santa Cruz Extreme AO Laboratory (SEAL) bench where it is now used as the standard way to invert nPWFS signal.

The Nancy Grace Roman Space Telescope Coronagraph Instrument is a critical technology demonstrator for NASA's Habitable Worlds Observatory. With a predicted visible-light flux ratio detection limit of 1E-8 or better, it will be capable of reaching new areas of parameter space for both gas giant exoplanets and circumstellar disks. It is in the final stages of integration and test at the Jet Propulsion Laboratory, with an anticipated delivery to payload integration in the coming year. This paper will review the instrument systems, observing modes, potential observing applications, and overall progress toward instrument integration and test., Comment: Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

While radial velocity surveys have demonstrated that the population of gas giants peaks around $3~\text{au}$, the most recent high-contrast imaging surveys have only been sensitive to planets beyond $\sim~10~\text{au}$. Sensitivity at small angular separations from stars is currently limited by the variability of the point spread function. We demonstrate how moderate-resolution integral field spectrographs can detect planets at smaller separations ($\lesssim~0.3$ arcseconds) by detecting the distinct spectral signature of planets compared to the host star. Using OSIRIS ($R$ $\approx$ 4000) at the W. M. Keck Observatory, we present the results of a planet search via this methodology around 20 young targets in the Ophiuchus and Taurus star-forming regions. We show that OSIRIS can outperform high-contrast coronagraphic instruments equipped with extreme adaptive optics and non-redundant masking in the $0.05-0.3$ arcsecond regime. As a proof of concept, we present the $34\sigma$ detection of a high-contrast M dwarf companion at $\approx0.1$" with a flux ratio of $\approx0.92\%$ around the field F2 star HD 148352. We developed an open-source Python package, breads, for the analysis of moderate-resolution integral field spectroscopy data in which the planet and the host star signal are jointly modeled. The diffracted starlight continuum is forward-modeled using a spline model, which removes the need for prior high-pass filtering or continuum normalization. The code allows for analytic marginalization of linear hyperparameters, simplifying posterior sampling of other parameters (e.g., radial velocity, effective temperature). This technique could prove very powerful when applied to integral field spectrographs like NIRSpec on the JWST and other upcoming first-light instruments on the future Extremely Large Telescopes., Comment: Accepted for publication in the Astronomical Journal on May 12, 2023

We reanalyze near-infrared spectra of the young extrasolar giant planet 51 Eridani b which was originally presented in (Macintosh et al. 2015) and (Rajan et al. 2017) using modern atmospheric models which include a self-consistent treatment of disequilibrium chemistry due to turbulent vertical mixing. In addition, we investigate the possibility that significant opacity from micrometeors or other impactors in the planet's atmosphere may be responsible for shaping the observed spectral energy distribution (SED). We find that disequilibrium chemistry is useful for describing the mid-infrared colors of the planet's spectra, especially in regards to photometric data at M band around 4.5 $\mu$m which is the result of super-equilibrium abundances of carbon monoxide, while the micrometeors are unlikely to play a pivotal role in shaping the SED. The best-fitting, micrometeroid-dust-free, disequilibrium chemistry, patchy cloud model has the following parameters: effective temperature $T_\textrm{eff} = 681$ K with clouds (or without clouds, i.e. the grid temperature $T_\textrm{grid}$ = 900 K), surface gravity $g$ = 1000 m/s$^2$, sedimentation efficiency $f_\textrm{sed}$ = 10, vertical eddy diffusion coefficient $K_\textrm{zz}$ = 10$^3$ cm$^2$/s, cloud hole fraction $f_\textrm{hole}$ = 0.2, and planet radius $R_\textrm{planet}$ = 1.0 R$_\textrm{Jup}$., Comment: 22 pages, 14 figures, Accepted to AJ

The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the orbital fit and, in combination with absolute astrometry from GAIA, place a 3$\sigma$ upper limit of 18.5 $M_J$ on the companion's mass. To interpret the ALMA observations, we used radiative transfer modeling to determine the disk properties. We find that the disk is well resolved and nearly edge on. While the misalignment angle between the disk and the orbit is weakly constrained due to the short orbital arc available, the data strongly support a (near) coplanar geometry for the system. Furthermore, we find that the models that describe the ALMA data best have an inner radius that is close to the companion's semi-major axis. Including a posteriori knowledge of the system's SED further narrows the constraints on the disk's inner radius and place it at a location that is in reasonable agreement with, possibly interior to, predictions from existing dynamical models of disk truncation by an interior substellar companion. HR\,2562 has the potential over the next few years to become a new testbed for dynamical interaction between a debris disk and a substellar companion., Comment: 22 pages, 17 figures

The carbon-to-oxygen (C/O) ratio in an exoplanet atmosphere has been suggested as a potential diagnostic of planet formation. Now that a number of exoplanets have measured C/O ratios, it is possible to examine this diagnostic at a population level. Here, we present an analysis of currently measured C/O ratios of directly imaged and transit/eclipse planets. First, we derive atmospheric parameters for the substellar companion HD 284149 AB b using data taken with the OSIRIS integral field spectrograph at the W.M. Keck Observatory and report two non-detections from our ongoing imaging spectroscopy survey with Keck/OSIRIS. We find an effective temperature of $T_\mathrm{eff} = 2502$ K, with a range of 2291-2624 K, $\log g=4.52$, with a range of 4.38-4.91, and [M/H] = 0.37, with a range of 0.10-0.55. We derive a C/O of 0.59$^{+0.15}_{-0.30}$ for HD 284149 AB b. We add this measurement to the list of C/O ratios for directly imaged planets and compare them with those from a sample of transit/eclipse planets. We also derive the first dynamical mass estimate for HD 284149 AB b, finding a mass of $\sim$28 $M_\mathrm{Jup}$. There is a trend in C/O ratio with companion mass ($M_{\mathrm{Jup}}$), with a break seen around 4$M_{\mathrm{Jup}}$. We run a Kolmogorov-Smirnov and an Anderson-Darling test on planets above and below this mass boundary, and find that they are two distinct populations. This could be additional evidence of two distinct populations possibly having two different formation pathways, with companion mass as an indicator of most likely formation scenario., Comment: 26 pages, 16 figures

The Earth's turbulent atmosphere results in speckled and blurred images of astronomical objects when observed by ground based visible and near-infrared telescopes. Adaptive optics (AO) systems are employed to reduce these atmospheric effects by using wavefront sensors (WFS) and deformable mirrors. Some AO systems are not fast enough to correct for strong, fast, high turbulence wind layers leading to the wind butterfly effect, or wind-driven halo, reducing contrast capabilities in coronagraphic images. Estimating the effective wind speed of the atmosphere allows us to calculate the atmospheric coherence time. This is not only an important parameter to understand for site characterization but could be used to help remove the wind butterfly in post processing. Here we present a method for estimating the atmospheric effective wind speed from spatio-temporal covariance maps generated from pseudo open-loop (POL) WFS data. POL WFS data is used as it aims to reconstruct the full wavefront information when operating in closed-loop. The covariance maps show how different atmospheric turbulent layers traverse the telescope. Our method successfully recovered the effective wind speed from simulated WFS data generated with the soapy python library. The simulated atmospheric turbulence profiles consist of two turbulent layers of ranging strengths and velocities. The method has also been applied to Gemini Planet Imager (GPI) AO WFS data. This gives insight into how the effective wind speed can affect the wind-driven halo seen in the AO image point spread function. In this paper, we will present results from simulated and GPI WFS data., Comment: 7 pages, 6 figures

Accreting protoplanets represent a window into planet formation processes. We report H{\alpha} differential imaging results from the deepest and most comprehensive accreting protoplanet survey to date, acquired with the Magellan Adaptive Optics (MagAO) system's VisAO camera. The fourteen transitional disks targeted are ideal candidates for protoplanet discovery due to their wide, heavily depleted central cavities, wealth of non-axisymmetric circumstellar disk features evocative of ongoing planet formation, and ongoing stellar accretion. To address the twin challenges of morphological complexity in the target systems and PSF instability, we develop novel approaches for frame selection and optimization of the Karhounen-Loeve Image Processing algorithm pyKLIP. We detect one new candidate protoplanet, CS Cha "c", at a separation of 75mas and a Delta mag of 5.1 and robustly recover the HD142527 B and HD100453 B low mass stellar companions across multiple epochs. Though we cannot rule out a substantial scattered light contribution to its emission, we also recover LkCa 15 b. Its presence inside of the cleared disk cavity and consistency with a forward-modeled point source suggest that it remains a viable protoplanet candidate. The protoplanet PDS 70 c was marginally recovered under our conservative general methodology. However, through targeted optimization in H-alpha} imagery, we tentatively recover PDS 70 c in three epochs and PDS 70 b in one epoch. Of the many other previously-reported companions and companion candidates around objects in the sample, we do not recover any additional robust candidates. However, lack of recovery at moderate H-alpha contrast does not rule out the presence of protoplanets at these locations, and we report limiting H-alpha contrasts in such cases., Comment: Final version, in press at AJ

The HR 8799 system hosts four massive planets orbiting 15 and 80 AU. Studies of the system's orbital stability and its outer debris disk open the possibility of additional planets, both interior to and exterior to the known system. Reaching a sufficient sensitivity to search for interior planets is very challenging due to the combination of bright quasi static speckle noise close to the stellar diffraction core and relatively fast orbital motion. In this work, we present a deep L-band imaging campaign using NIRC2 at Keck comprising 14 observing sequences. We further re-reduce archival data for a total of 16.75 hours, one of the largest uniform datasets of a single direct imaging target. Using a Bayesian modeling technique for detecting planets in images while compensating for plausible orbital motion, we then present deep limits on the existence of additional planets in the HR 8799 system. The final combination shows a tentative candidate, consistent with 4-7 $M_{jup}$ at 4-5 AU, detected with an equivalent false alarm probability better than $3\sigma$. This analysis technique is widely applicable to archival data and to new observations from upcoming missions that revisit targets at multiple epochs., Comment: 26 pages, 18 figures. Accepted to AJ. Data available: DOI 10.5281/zenodo.6823071. Updated with corrected citations and higher resolution stability plots

The Gemini Planet Imager (GPI) is an integral field spectrograph (IFS) and coronagraph that is one of the few current generation instruments optimized for high-contrast direct imaging of substellar companions. The instrument is in the process of being upgraded and moved from its current mount on the Gemini South Observatory in Cerro Pachon, Chile, to its twin observatory, Gemini North, on Mauna Kea (a process colloquially dubbed 'GPI 2.0'). We present the designs that have been developed for the part of GPI 2.0 that pertains to upgrading various optical components of the GPI coronagraphic system. More specifically, we present new designs for the apodizer and Lyot stop (LS) that achieve better raw contrast at the inner working angle of the dark zone as well as improved core throughput while retaining a similar level of robustness to LS misalignment. To generate these upgraded designs, we use our own publicly available software package called APLC-Optimization that combines a commercial linear solver (Gurobi) with a high contrast imaging simulation package (HCIPy) in order to iteratively propagate light through a simulated model of an apodized phase lyot coronagraph (APLC), optimizing for the best coronagraph performance metrics. The designs have recently finished being lithographically printed by a commercial manufacturer and will be ready for use when GPI 2.0 goes on-sky in 2023., Comment: Conference proceedings published in SPIE Astronomical Telescopes + Instrumentation (July 2022), 12 pages, 7 figures

We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $\mu$m to 20 $\mu$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion., Comment: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $\mu$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $\mu$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $\mu$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$\sigma$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $\mu$m and MIRI at 11.3 $\mu$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $\mu$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail., Comment: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

Direct imaging characterization of extrasolar planets is often done at low spectral resolution. We model the spectrograph for the Gemini Planet Imager upgrade (GPI 2.0) and assess the instrument's potential for allowing observers to constrain exoplanet properties through analysis of near-infrared spectra. We simulated noisy observations followed by calculations of posterior distributions from maximum likelihood comparison with the Sonora 2018 model grid. Preliminary results suggest that GPI 2.0 should allow observers to constrain temperature with sufficient accuracy, but gravity remains largely uncertain. We also explore the effects of incorporating convolution with the instrument line spread function into our simulation and compare the results with our preliminary findings., Comment: 11 pages, 10 figures, to be published in proceedings for SPIE Astronomical Telescopes + Instrumentation 2022

The atmospheres of gas giant planets are thought to be inhomogeneous due to weather and patchy clouds. We present two full nights of coronagraphic observations of the HR 8799 planets using the CHARIS integral field spectrograph behind the SCExAO adaptive optics system on the Subaru Telescope to search for spectrophomometric variability. We did not detect significant variability signals, but placed the lowest variability upper limits for HR 8799 c and d. Based on injection-recovery tests, we expected to have a 50% chance to detect signals down to 10% H-band photometric variability for HR 8799 c and down to 30% H-band variability for HR 8799 d. We also investigated spectral variability and expected a 50% chance to recovery 20% variability in the H/K flux ratio for HR 8799 c. We combined all the data from the two nights to obtain some of the most precise spectra obtained for HR 8799 c, d, and e. Using a grid of cloudy radiative-convective-thermochemical equilibrium models, we found all three planets prefer supersolar metallicity with effective temperatures of ~1100 K. However, our high signal-to-noise spectra show that HR 8799 d has a distinct spectrum from HR 8799 c, possibly preferring more vertically extended and uniform clouds and indicating that the planets are not identical., Comment: 21 pages, 12 figures, Accepted to AJ, time series photometry and stacked spectra available as text files

We present moderate-resolution ($R\sim4000$) $K$ band spectra of the planetary-mass companion VHS 1256 b. The data were taken with the OSIRIS integral field spectrograph at the W.M. Keck Observatory. The spectra reveal resolved molecular lines from H$_{2}$O and CO. The spectra are compared to custom $PHOENIX$ atmosphere model grids appropriate for young, substellar objects. We fit the data using a Markov Chain Monte Carlo forward modeling method. Using a combination of our moderate-resolution spectrum and low-resolution, broadband data from the literature, we derive an effective temperature of 1240 K, with a range of 1200--1300 K, a surface gravity of $\log{g}=$ 3.25, with a range of 3.25--3.75 and a cloud parameter of $\log P_{cloud}=$ 6, with a range of 6.0--6.6. These values are consistent with previous studies, regardless of the new, larger system distance from GAIA EDR3 (22.2$^{+1.1}_{-1.2}$ pc). We derive a C/O ratio of 0.590$_{-0.354}^{+0.280}$ for VHS 1256b. Both our OSIRIS data and spectra from the literature are best modeled when using a larger 3 $\mu$m grain size for the clouds than used for hotter objects, consistent with other sources in the L/T transition region. VHS 1256 b offers an opportunity to look for systematics in the modeling process that may lead to the incorrect derivation of properties like C/O ratio in the high contrast regime., Comment: 16 pages, 10 figures. arXiv admin note: text overlap with arXiv:2009.08959

The Nancy Grace Roman Space Telescope's (Roman) Coronagraph Instrument is a technology demonstration equipped to achieve flux contrast levels of up to 10$^{-9}$. This precision depends upon the quality of observations and their resultant on-sky corrections via an absolute flux calibration (AFC). Our plan utilizes 10 dim and 4 bright standard photometric calibrator stars from Hubble Space Telescope's (HST) CALSPEC catalog to yield a final AFC error of 1.94\% and total observation time of $\sim$22 minutes. Percent error accounts for systematic uncertainties (filters, upstream optics, quantum efficiency) in Roman component instrumentation along with shot noise for a signal to noise ratio (SNR) of 500., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2022

The prospect of combining integral field spectroscopy with the solar gravitational lens (SGL) to spectrally and spatially resolve the surfaces and atmospheres of extrasolar planets is investigated. The properties of hyperbolic orbits visiting the focal region of the SGL are calculated analytically, demonstrating trade offs between departure velocity and time of arrival, as well as gravity assist maneuvers and heliocentric angular velocity. Numerical integration of the solar barycentric motion demonstrates that navigational acceleration of $\textrm{d}v \lesssim 80 \frac{\textrm{m}}{\textrm{s}} + 6.7 \frac{\textrm{m}}{\textrm{s}} \frac{t}{\textrm{year}}$ is needed to obtain and maintain alignment. Obtaining target ephemerides of sufficient precision is an open problem. The optical properties of an oblate gravitational lens are reviewed, including calculations of the magnification and the point-spread function that forms inside a telescope. Image formation for extended, incoherent sources is discussed when the projected image is smaller than, approximately equal to, and larger than the critical caustic. Sources of contamination which limit observational SNR are considered in detail, including the sun, the solar corona, the host star, and potential background objects. A noise mitigation strategy of spectrally and spatially separating the light using integral field spectroscopy is emphasized. A pseudoinverse-based image reconstruction scheme demonstrates that direct reconstruction of an Earth-like source from \textit{single} measurements of the Einstein ring is possible when the critical caustic and observed SNR are sufficiently large. In this arrangement, a mission would not require multiple telescopes or navigational symmetry breaking, enabling continuous monitoring of the atmospheric composition and dynamics on other planets., Comment: 34 pages, 21 figures, accepted to ApJ

We evaluate the potential for gravitational-wave (GW) detection in the frequency band from 10 nHz to 1 $\mu$Hz using extremely high-precision astrometry of a small number of stars. In particular, we argue that non-magnetic, photometrically stable hot white dwarfs (WD) located at $\sim$ kpc distances may be optimal targets for this approach. Previous studies of astrometric GW detection have focused on the potential for less precise surveys of large numbers of stars; our work provides an alternative optimization approach to this problem. Interesting GW sources in this band are expected at characteristic strains around $h_c \sim 10^{-17} \times \left(\mu\text{Hz}/f_{\text{GW}}\right)$. The astrometric angular precision required to see these sources is $\Delta \theta \sim h_c$ after integrating for a time $T \sim 1/f_{\text{GW}}$. We show that jitter in the photometric center of WD of this type due to starspots is bounded to be small enough to permit this high-precision, small-$N$ approach. We discuss possible noise arising from stellar reflex motion induced by orbiting objects and show how it can be mitigated. The only plausible technology able to achieve the requisite astrometric precision is a space-based stellar interferometer. Such a future mission with few-meter-scale collecting dishes and baselines of $\mathcal{O}(100\text{ km})$ is sufficient to achieve the target precision. This collector size is broadly in line with the collectors proposed for some formation-flown, space-based astrometer or optical synthetic-aperture imaging-array concepts proposed for other science reasons. The proposed baseline is however somewhat larger than the km-scale baselines discussed for those concepts, but we see no fundamental technical obstacle to utilizing such baselines. A mission of this type thus also holds the promise of being one of the few ways to access interesting GW sources in this band., Comment: 19 pages. Published version

NASA's next flagship mission, the Nancy Grace Roman Space Telescope, is a 2.4-meter observatory set to launch no later than May 2027. Roman features two instruments: the Wide Field Imager and the Coronagraph Instrument. Roman's Coronagraph is a Technology Demonstration that will push the current capabilities of direct imaging to smaller contrast ratios ($\sim$10$^{-9}$) and inner-working angles (3~$\lambda$/D). In order to achieve this high precision, Roman Coronagraph data must be calibrated to remove as many potential sources of error as possible. Here we present a detailed overview of the Nancy Grace Roman Space Telescope Coronagraph Instrument Observation Calibration Plan including identifying potential sources of error and how they will be mitigated via on-sky calibrations., Comment: Posting for public information on the current status of the Roman Coronagraph Observation Calibration Plan; latest updates as of July 29, 2022

Exozodiacal dust, warm debris from comets and asteroids in and near the habitable zone of stellar systems, reveals the physical processes that shape planetary systems. Scattered light from this dust is also a source of background flux which must be overcome by future missions to image Earthlike planets. This study quantifies the sensitivity of the Nancy Grace Roman Space Telescope Coronagraph to light scattered by exozodi, the zodiacal dust around other stars. Using a sample of 149 nearby stars, previously selected for optimum detection of habitable exoplanets by space observatories, we find the maximum number of exozodiacal disks with observable \textit{inner} habitable zone boundaries is six and the number of observable outer habitable boundaries is 74. One zodi was defined as the visible-light surface brightness of 22 $m_{\rm V}\ $arcsec$^{-2}$ around a solar-mass star, approximating the scattered light brightness in visible light at the Earth-equivalent insolation. In the speckle limited case, where the signal-to-noise ratio is limited by speckle temporal stability rather than shot noise, the median $5\sigma$ sensitivity to habitable zone exozodi is 12 zodi per resolution element. This estimate is calculated at the inner-working angle of the coronagraph, for the current best estimate performance, neglecting margins on the uncertainty in instrument performance and including a post-processing speckle suppression factor. For an log-norm distribution of exozodi levels with a median exozodi of 3$\times$ the solar zodi, we find that the Roman Coronagraph would be able to make 5$\sigma$ detections of exozodiacal disks in scattered light from 13 systems with a 95\% confidence interval spanning 7-20 systems. This sensitivity allows Roman Coronagraph to complement ground-based measurements of exozodiacal thermal emission and constrain dust albedos., Comment: 30 pages, 8 figures, 2 tables, accepted to PASP

Astrometry has long been a promising technique for exoplanet detection. At the theoretical limits, astrometry would allow for the detection of smaller planets than previously seen by current exoplanet search methods, but stellar activity may make these theoretical limits unreachable. Astrometric jitter of a Sun-like star due to magnetic activity in its photosphere induces apparent variability in the photocenter of order $0.5\ \textrm{m}R_\odot$. This jitter creates a fundamental astrophysical noise floor preventing detection of lower-mass planets in a single spectral band. By injecting planet orbits into simulated solar data at five different passbands, we investigate mitigation of this fundamental astrometric noise using correlations across passbands. For a true solar analog and a planet at 1 au semimajor axis, the $6\sigma$ detection limit set by stellar activity for an ideal telescope at the best single passband is $0.01$ Earth masses. We found that pairs of passbands with highly correlated astrometric jitter due to stellar activity, but with less motion in the redder band, enable higher-precision measurements of the common signal from the planet. Using this method improves detectable planet masses at 1 au by up to a factor of $10$, corresponding to at best $0.005$ Earth masses for a Sun-like star with a perfect telescope. Given these results, we recommend that future astrometry missions consider proceeding with two or more passbands to reduce noise due to stellar activity., Comment: 11 pages, 9 figures. Published in AJ

The four directly imaged planets orbiting the star HR 8799 are an ideal laboratory to probe atmospheric physics and formation models. We present more than a decade's worth of Keck/OSIRIS observations of these planets, which represent the most detailed look at their atmospheres to-date by its resolution and signal to noise ratio. We present the first direct detection of HR 8799 d, the second-closest known planet to the star, at moderate spectral resolution with Keck/OSIRIS (K-band; R~4,000). Additionally, we uniformly analyze new and archival OSIRIS data (H and K band) of HR 8799 b, c, and d. First, we show detections of water (H2O) and carbon monoxide (CO) in the three planets and discuss the ambiguous case of methane (CH4) in the atmosphere of HR 8799b. Then, we report radial velocity (RV) measurements for each of the three planets. The RV measurement of HR 8799 d is consistent with predictions made assuming coplanarity and orbital stability of the HR 8799 planetary system. Finally, we perform a uniform atmospheric analysis on the OSIRIS data, published photometric points, and low resolution spectra. We do not infer any significant deviation from to the stellar value of the carbon to oxygen ratio (C/O) of the three planets, which therefore does not yet yield definitive information about the location or method of formation. However, constraining the C/O ratio for all the HR 8799 planets is a milestone for any multiplanet system, and particularly important for large, widely separated gas giants with uncertain formation processes., Comment: Accepted to AJ

The Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope will demonstrate the high-contrast technology necessary for visible-light exoplanet imaging and spectroscopy from space via direct imaging of Jupiter-size planets and debris disks. This in-space experience is a critical step toward future, larger missions targeted at direct imaging of Earth-like planets in the habitable zones of nearby stars. This paper presents an overview of the current instrument design and requirements, highlighting the critical hardware, algorithms, and operations being demonstrated. We also describe several exoplanet and circumstellar disk science cases enabled by these capabilities. A competitively selected Community Participation Program team will be an integral part of the technology demonstration and could perform additional CGI observations beyond the initial tech demo if the instrument performance warrants it., Comment: SPIE Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave, 114431U (15 December 2020)

Launching a starshade to rendezvous with the Nancy Grace Roman Space Telescope would provide the first opportunity to directly image the habitable zones of nearby sunlike stars in the coming decade. A report on the science and feasibility of such a mission was recently submitted to NASA as a probe study concept. The driving objective of the concept is to determine whether Earth-like exoplanets exist in the habitable zones of the nearest sunlike stars and have biosignature gases in their atmospheres. With the sensitivity provided by this telescope, it is possible to measure the brightness of zodiacal dust disks around the nearest sunlike stars and establish how their population compares to our own. In addition, known gas-giant exoplanets can be targeted to measure their atmospheric metallicity and thereby determine if the correlation with planet mass follows the trend observed in the Solar System and hinted at by exoplanet transit spectroscopy data. In this paper we provide the details of the calculations used to estimate the sensitivity of Roman with a starshade and describe the publicly available Python-based source code used to make these calculations. Given the fixed capability of Roman and the constrained observing windows inherent for the starshade, we calculate the sensitivity of the combined observatory to detect these three types of targets and we present an overall observing strategy that enables us to achieve these objectives., Comment: Accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems (JATIS)

We explore the impact of outer stellar companions on the occurrence rate of giant planets detected with radial velocities. We searched for stellar and planetary companions to a volume-limited sample of solar-type stars within 25 pc. Using adaptive optics imaging from the Lick 3m and Palomar 200" Telescopes, we characterized the multiplicity of our sample stars, down to the bottom of the main sequence. With these data, we confirm field star multiplicity statistics from previous surveys. We combined three decades of radial velocity data from the California Planet Search with new RV data from Keck/HIRES and APF/Levy to search for planets in the same systems. Using an updated catalog of both stellar and planetary companions and injection/recovery tests to determine our sensitivity, we measured the occurrence rate of planets among the single and multiple star systems. We found that planets with masses of 0.1-10 $M_{Jup}$ and semi-major axes of 0.1-10 AU have an occurrence rate of $0.18^{+0.04}_{-0.03}$ planets per single star, and $0.12\pm0.04$ planets per binary primary. Only one planet-hosting binary system in our sample had a binary separation $<100$ AU, and none had a separation $<50$ AU. We found planet occurrence rates of $0.20^{+0.07}_{-0.06}$ planets per star for binaries with separation $a_B > 100$ AU, and $0.04^{+0.04}_{-0.02}$ planets per star for binaries with separation $a_B<100$ AU. The similarity in the planet occurrence rate around single stars and wide primaries implies that wide binary systems should host more planets than single star systems, since they have more potential host stars. We estimated a system-wide planet occurrence rate of 0.3 planets per wide binary system for binaries with separations $a_B > 100$ AU. Finally, we found evidence that giant planets in binary systems have a different semi-major axis distribution than their counterparts in single star systems., Comment: 58 pages, 13 figures, accepted to The Astronomical Journal

During its move from the mountaintop of Cerro Pachon in Chile to the peak of Mauna Kea in Hawaii, the Gemini Planet Imager will make a pit stop to receive various upgrades, including a pyramid wavefront sensor. As a highly non-linear sensor, a standard approach to linearize the response of the pyramid is induce a rapid circular modulation of the beam around the pyramid tip, trading off sensitivity for robustness during high turbulence. Using high temporal resolution Fourier Optics based simulations, we investigate phase reconstruction approaches that attempt to optimize the performance of the sensor with a dynamically adjustable modulation parameter. We have studied the linearity and gain stability of the sensor under different modulation and seeing conditions, and the ability of the sensor to correct non-common-path errors. We will also show performance estimates which includes a comparative analysis of the atmospheric columns above the two mountains, as well as the Error Transfer Functions of the two systems., Comment: 20 pages, 19 figures. To appear in Proceedings of SPIE Adaptive Optics Systems VII

The Gemini Planet Imager (GPI) is a high-contrast imaging instrument designed to directly image and characterize exoplanets. GPI is currently undergoing several upgrades to improve performance. In this paper, we discuss the upgrades to the GPI IFS. This primarily focuses on the design and performance improvements of new prisms and filters. This includes an improved high-resolution prism which will provide more evenly dispersed spectra across y, J, H and K-bands. Additionally, we discuss the design and implementation of a new low-resolution mode and prism which allow for imaging of all four bands (y, J, H and K-bands) simultaneously at R=10. We explore the possibility of using a multiband filter which would block the light between the four spectral bands. We discuss possible performance improvements from the multiband filter, if implemented. Finally we explore the possibility of making small changes to the optical design to improve the IFS's performance near the edge of the field of view., Comment: 8 pages, 5 figures, Proc. of SPIE Paper No. 11447-410

We present moderate-resolution ($R\sim4000$) $K$ band spectra of the "super-Jupiter," $\kappa$ Andromedae b. The data were taken with the OSIRIS integral field spectrograph at Keck Observatory. The spectra reveal resolved molecular lines from H$_{2}$O and CO. The spectra are compared to a custom $PHOENIX$ atmosphere model grid appropriate for young planetary-mass objects. We fit the data using a Markov Chain Monte Carlo forward modeling method. Using a combination of our moderate-resolution spectrum and low-resolution, broadband data from the literature, we derive an effective temperature of $T_\mathrm{eff}$ = 1950 - 2150 K, a surface gravity of $\log g=3.5 - 4.5$, and a metallicity of [M/H] = $-0.2 - 0.0$. These values are consistent with previous estimates from atmospheric modeling and the currently favored young age of the system ($<$50 Myr). We derive a C/O ratio of 0.70$_{-0.24}^{+0.09}$ for the source, broadly consistent with the solar C/O ratio. This, coupled with the slightly subsolar metallicity, implies a composition consistent with that of the host star, and is suggestive of formation by a rapid process. The subsolar metallicity of $\kappa$ Andromedae b is also consistent with predictions of formation via gravitational instability. Further constraints on formation of the companion will require measurement of the C/O ratio of $\kappa$ Andromedae A. We also measure the radial velocity of $\kappa$ Andromedae b for the first time, with a value of $-1.4\pm0.9\,\mathrm{km}\,\mathrm{s}^{-1}$ relative to the host star. We find that the derived radial velocity is consistent with the estimated high eccentricity of $\kappa$ Andromedae b., Comment: 23 pages, 18 figures. Accepted for publication in AJ

Comparing chemical abundances of a planet and the host star reveals the origin and formation path. Stellar abundance is measured with high-resolution spectroscopy. Planet abundance, on the other hand, is usually inferred from low-resolution data. For directly imaged exoplanets, the data are available from a slew of high-contrast imaging/spectroscopy instruments. Here, we study the chemical abundance of HR 8799 and its planet c. We measure stellar abundance using LBT/PEPSI (R=120,000) and archival HARPS data: stellar [C/H], [O/H], and C/O are 0.11$\pm$0.12, 0.12$\pm$0.14, and 0.54$^{+0.12}_{-0.09}$, all consistent with solar values. We conduct atmospheric retrieval using newly obtained Subaru/CHARIS data together with archival Gemini/GPI and Keck/OSIRIS data. We model the planet spectrum with petitRADTRANS and conduct retrieval using PyMultiNest. Retrieved planetary abundance can vary by $\sim$0.5 dex, from sub-stellar to stellar C and O abundances. The variation depends on whether strong priors are chosen to ensure a reasonable planet mass. Moreover, comparison with previous works also reveals inconsistency in abundance measurements. We discuss potential issues that can cause the inconsistency, e.g., systematics in individual data sets and different assumptions in the physics and chemistry in retrieval. We conclude that no robust retrieval can be obtained unless the issues are fully resolved., Comment: Published on AJ

HR4796A hosts a well-studied debris disk with a long history due to its high fractional luminosity and favorable inclination lending itself well to both unresolved and resolved observations. We present new J- and K1-band images of the resolved debris disk HR4796A taken in the polarimetric mode of the Gemini Planet Imager (GPI). The polarized intensity features a strongly forward scattered brightness distribution and is undetected at the far side of the disk. The total intensity is detected at all scattering angles and also exhibits a strong forward scattering peak. We use a forward modelled geometric disk in order to extract geometric parameters, polarized fraction and total intensity scattering phase functions for these data as well as H-band data previously taken by GPI. We find the polarized phase function becomes increasingly more forward scattering as wavelength increases. We fit Mie and distribution of hollow spheres grain (DHS) models to the extracted functions. We find that while it is possible to describe generate a satisfactory model for the total intensity using a DHS model, but not with a Mie model. We find that no single grain population of DHS or Mie grains of arbitrary composition can simultaneously reproduce the polarized fraction and total intensity scattering phase functions, indicating the need for more sophisticated grain models., Comment: Accepted for publication in AJ. 3 tables, 11 figures

Age is a fundamental parameter of stars, yet in many cases ages of individual stars are presented without robust estimates of the uncertainty. We have developed a Bayesian framework, BAFFLES, to produce the age posterior for a star from its calcium emission strength (log($R'_{HK}$)) or lithium abundance (Li EW) and $B-V$ color. We empirically determine the likelihood functions for calcium and lithium as functions of age from literature measurements of stars in benchmark clusters with well-determined ages. We use a uniform prior on age which reflects a uniform star formation rate. The age posteriors we derive for several test cases are consistent with literature ages found from other methods. BAFFLES represents a robust method to determine the age posterior probability distribution for any field star with $0.45 \leq B-V \leq 0.9$ and a measurement of $R'_{HK}$ and/or $0.35 \leq B-V \leq 1.9$ and measured Li EW. We compile colors, $R'_{HK}$, and Li EW from over 2630 nearby field stars from the literature and present the derived BAFFLES age posterior for each star., Comment: Accepted to ApJ, 80 pages, 14 figures, 3 tables

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $\mu$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 \deg$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present., Comment: arXiv resubmission with typographical corrections. Accepted for publication in AJ. 19 figures, 7 tables