Showing 1–50 of 87 results for author: Greenbaum, A

Current laboratory performance of starlight suppression systems, and potential pathways to desired Habitable Worlds Observatory exoplanet science capabilities

Authors: Bertrand Mennesson, Ruslan Belikov, Emiel Por, Eugene Serabyn, Garreth Ruane, A. J. Eldorado Riggs, Dan Sirbu, Laurent Pueyo, Remi Soummer, Jeremy Kasdin, Stuart Shaklan, Byoung-Joon Seo, Christopher Stark, Eric Cady, Pin Chen, Brendan Crill, Kevin Fogarty, Alexandra Greenbaum, Olivier Guyon, Roser Juanola-Parramon, Brian Kern, John Krist, Bruce Macintosh, David Marx, Dimitri Mawet , et al. (12 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: 63 pages, 28 pages, submitted to JATIS

The James Webb Interferometer: Space-based interferometric detections of PDS 70 b and c at 4.8 $μ$m

Authors: Dori Blakely, Doug Johnstone, Gabriele Cugno, Anand Sivaramakrishnan, Peter Tuthill, Ruobing Dong, Benjamin J. S. Pope, Loïc Albert, Max Charles, Rachel A. Cooper, Matthew De Furio, Louis Desdoigts, René Doyon, Logan Francis, Alexandra Z. Greenbaum, David Lafrenière, James P. Lloyd, Michael R. Meyer, Laurent Pueyo, Shrishmoy Ray, Joel Sánchez-Bermúdez, Anthony Soulain, Deepashri Thatte, Thomas Vandal

Abstract: We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of… ▽ More We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of both PDS 70 b and c provide evidence for circumplanetary disk emission through fitting SED models to these new measurements and those found in the literature. We also newly detect emission within the disk gap at an SNR of $\sim$4, at a position angle of $207^{+11}_{-10}$ degrees, and an unconstrained separation within $\sim$200 mas. Follow-up observations will be needed to determine the nature of this emission. We place a 5$σ$ upper limit of $Δ$mag = 7.56 on the contrast of the candidate PDS 70 d at 4.8 $μ$m, which indicates that if the previously observed emission at shorter wavelengths is due to a planet, this putative planet has a different atmospheric composition than PDS 70 b or c. Finally, we place upper limits on emission from any additional planets in the disk gap. We find an azimuthally averaged 5$σ$ upper limit of $Δ$mag $\approx$ 7.5 at separations greater than 125 mas. These are the deepest limits to date within $\sim$250 mas at 4.8 $μ$m and the first space-based interferometric observations of this system. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

High-precision atmospheric characterization of a Y dwarf with JWST NIRSpec G395H spectroscopy: isotopologue, C/O ratio, metallicity, and the abundances of six molecular species

Authors: Ben W. P. Lew, Thomas Roellig, Natasha E. Batalha, Michael Line, Thomas Greene, Sagnick Murkherjee, Richard Freedman, Michael Meyer, Charles Beichman, Catarina Alves De Oliveira, Matthew De Furio, Doug Johnstone, Alexandra Z. Greenbaum, Mark Marley, Jonathan J. Fortney, Erick T. Young, Jarron Leisenring, Martha Boyer, Klaus Hodapp, Karl Misselt, John Stansberry, Marcia Rieke

Abstract: The launch of the James Webb Space Telescope (JWST) marks a pivotal moment for precise atmospheric characterization of Y dwarfs, the coldest brown dwarf spectral type. In this study, we leverage moderate spectral resolution observations (R $\sim$ 2700) with the G395H grating of the Near-Infrared Spectrograph (NIRSpec) onboard of JWST to characterize the nearby (9.9 pc) Y dwarf WISEPA J182831.08+26… ▽ More The launch of the James Webb Space Telescope (JWST) marks a pivotal moment for precise atmospheric characterization of Y dwarfs, the coldest brown dwarf spectral type. In this study, we leverage moderate spectral resolution observations (R $\sim$ 2700) with the G395H grating of the Near-Infrared Spectrograph (NIRSpec) onboard of JWST to characterize the nearby (9.9 pc) Y dwarf WISEPA J182831.08+265037.8 (WISE 1828). With the NIRSpec G395H 2.88-5.12 $\mathrmμ$m spectrum, we measure the abundances of CO, CO$_2$, CH$_4$, H$_2$S, NH$_3$, and H$_2$O, which are the major carbon, nitrogen, oxygen, and sulfur bearing species in the atmosphere. Based on the retrieved volume mixing ratios with the atmospheric retrieval framework CHIMERA, we report that the C/O ratio is $0.45 \pm 0.01$, close to the solar C/O value of 0.55, and the metallicity to be +0.30 $\pm$ 0.02 dex. Comparison between the retrieval results with the forward modeling results suggests that the model bias for C/O and metallicity could be as high as 0.03 and 0.97 dex respectively. We also report a lower limit of the $^{12}$CO/$^{13}$CO ratio of $>40 $, being consistent with the nominal solar value of 90. Our results highlight the potential of JWST in measuring the C/O ratios down to percent-level precision and characterizing isotopologues of cold planetary atmospheres similar to WISE 1828. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 18 pages + references, including 11 figures, accepted for publication in The Astronomical Journal

CD-27 11535: Evidence for a Triple System in the $β$ Pictoris Moving Group

Authors: Andrew D. Thomas, Eric L. Nielsen, Robert J. De Rosa, Anne E. Peck, Bruce Macintosh, Jeffrey Chilcote, Paul Kalas, Jason J. Wang, Sarah Blunt, Alexandra Greenbaum, Quinn M. Konopacky, Michael J. Ireland, Peter Tuthill, Kimberly Ward-Duong, Lea A. Hirsch, Ian Czekala, Franck Marchis, Christian Marois, Max A. Millar-Blanchaer, William Roberson, Adam Smith, Hannah Gallamore, Jessica Klusmeyer

Abstract: We present new spatially resolved astrometry and photometry of the CD-27 11535 system, a member of the $β$ 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 discrepa… ▽ More We present new spatially resolved astrometry and photometry of the CD-27 11535 system, a member of the $β$ 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. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 12 pages + references and appendix, 12 figures, 5 tables

Journal ref: AJ 166 246 (2023)

A First Look with JWST Aperture Masking Interferometry (AMI): Resolving Circumstellar Dust around the Wolf-Rayet Binary WR 137 beyond the Rayleigh Limit

Authors: Ryan M. Lau, Matthew J. Hankins, Joel Sanchez-Bermudez, Deepashri Thatte, Anthony Soulain, Rachel A. Cooper, Anand Sivaramakrishnan, Michael F. Corcoran, Alexandra Z. Greenbaum, Theodore R. Gull, Yinuo Han, Olivia C. Jones, Thomas Madura, Anthony F. J. Moffat, Mark R. Morris, Takashi Onaka, Christopher M. P. Russell, Noel D. Richardson, Nathan Smith, Peter Tuthill, Kevin Volk, Gerd Weigelt, Peredur M. Williams

Abstract: We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and Augus… ▽ More We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and August as part of the Director's Discretionary Early Release Science (DD-ERS) program 1349. Interferometric observables (squared visibilities and closure phases) from the WR 137 "interferogram" were extracted and calibrated using three independent software tools: ImPlaneIA, AMICAL, and SAMpip. The analysis of the calibrated observables yielded consistent values except for slightly discrepant closure phases measured by ImPlaneIA. Based on all three sets of calibrated observables, images were reconstructed using three independent software tools: BSMEM, IRBis, and SQUEEZE. All reconstructed image combinations generated consistent images in both F380M and F480M filters. The reconstructed images of WR 137 reveal a bright central core with a $\sim300$ mas linear filament extending to the northwest. A geometric colliding-wind model with dust production constrained to the orbital plane of the binary system and enhanced as the system approaches periapsis provided a general agreement with the interferometric observables and reconstructed images. Based on a colliding-wind dust condensation analysis, we suggest that dust formation within the orbital plane of WR 137 is induced by enhanced equatorial mass-loss from the rapidly rotating O9 companion star, whose axis of rotation is aligned with that of the orbit. △ Less

Submitted 22 December, 2023; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 18 pages, 8 figures, Accepted for publication in ApJ. Updated plotting error in Fig. 2

Searching for Planets Orbiting Fomalhaut with JWST/NIRCam

Authors: Marie Ygouf, Charles Beichman, Jorge Llop-Sayson, Geoffrey Bryden, Jarron Leisenring, Andras Gaspar, John Krist, Marcia Rieke, George Rieke, Schuyler Wolff, Thomas Roellig, Kate Su, Kevin Hainline, Klaus Hodapp, Thomas Greene, Michael Meyer, Doug Kelly, Karl Misselt, John Stansberry, Martha Boyer, Doug Johnstone, Scott Horner, Alexandra Greenbaum

Abstract: We report observations with the JWST/NIRCam coronagraph of the Fomalhaut system. This nearby A star hosts a complex debris disk system discovered by the IRAS satellite. Observations in F444W and F356W filters using the round 430R mask achieve a contrast ratio of ~ 4 x 10-7 at 1'' and ~ 4 x 10-8 outside of 3''. These observations reach a sensitivity limit <1 MJup across most of the disk region. Con… ▽ More We report observations with the JWST/NIRCam coronagraph of the Fomalhaut system. This nearby A star hosts a complex debris disk system discovered by the IRAS satellite. Observations in F444W and F356W filters using the round 430R mask achieve a contrast ratio of ~ 4 x 10-7 at 1'' and ~ 4 x 10-8 outside of 3''. These observations reach a sensitivity limit <1 MJup across most of the disk region. Consistent with the hypothesis that Fomalhaut b is not a massive planet but is a dust cloud from a planetesimal collision, we do not detect it in either F356W or F444W (the latter band where a Jovian-sized planet should be bright). We have reliably detected 10 sources in and around Fomalhaut and its debris disk, all but one of which are coincident with Keck or HST sources seen in earlier coronagraphic imaging; we show them to be background objects, including the "Great Dust Cloud" identified in MIRI data. However, one of the objects, located at the edge of the inner dust disk seen in the MIRI images, has no obvious counterpart in imaging at earlier epochs and has a relatively red [F356W]-[F444W]>0.7 mag (Vega) color. Whether this object is a background galaxy, brown dwarf, or a Jovian mass planet in the Fomalhaut system will be determined by an approved Cycle 2 follow-up program. Finally, we set upper limits to any scattered light from the outer ring, placing a weak limit on the dust albedo at F356W and F444W. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 24 pages, 17 figures

JWST-TST High Contrast: Achieving direct spectroscopy of faint substellar companions next to bright stars with the NIRSpec IFU

Authors: Jean-Baptiste Ruffio, Marshall D. Perrin, Kielan K. W. Hoch, Jens Kammerer, Quinn M. Konopacky, Laurent Pueyo, Alex Madurowicz, Emily Rickman, Christopher A. Theissen, Shubh Agrawal, Alexandra Z. Greenbaum, Brittany E. Miles, Travis S. Barman, William O. Balmer, Jorge Llop-Sayson, Julien H. Girard, Isabel Rebollido, Rémi Soummer, Natalie H. Allen, Jay Anderson, Charles A. Beichman, Andrea Bellini, Geoffrey Bryden, Néstor Espinoza, Ana Glidden , et al. (11 additional authors not shown)

Abstract: The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5 um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starli… ▽ More The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5 um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starlight suppression techniques. We develop specific data reduction strategies to study faint companions around bright stars, and assess the performance of NIRSpec at high contrast. First, we demonstrate an approach to forward model the companion signal and the starlight directly in the detector images, which mitigates the effects of NIRSpec's spatial undersampling. We demonstrate a sensitivity to planets that are 3e-6 fainter than their stars at 1'', or 3e-5 at 0.3''. Then, we implement a reference star point spread function (PSF) subtraction and a spectral extraction that does not require spatially and spectrally regularly sampled spectral cubes. This allows us to extract a moderate resolution (R~2,700) spectrum of the faint T-dwarf companion HD 19467 B from 2.9-5.2 um with signal-to-noise ratio (S/N)~10 per resolution element. Across this wavelength range, HD~19467~B has a flux ratio varying between 1e-5-1e-4 and a separation relative to its star of 1.6''. A companion paper by Hoch et al. more deeply analyzes the atmospheric properties of this companion based on the extracted spectrum. Using the methods developed here, NIRSpec's sensitivity may enable direct detection and spectral characterization of relatively old (~1 Gyr), cool (~250 K), and closely separated (~3-5 au) exoplanets that are less massive than Jupiter. △ Less

Submitted 31 May, 2024; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: Accepted to AJ. The data analysis scripts for this work are published https://github.com/jruffio/HD_19467_B (https://zenodo.org/doi/10.5281/zenodo.11391740). The main revisions of the manuscript are listed in the change history section of the readme

Optimal Polynomial Approximation to Rational Matrix Functions Using the Arnoldi Algorithm

Authors: Tyler Chen, Anne Greenbaum, Natalie Wellen

Abstract: Given an $n$ by $n$ matrix $A$ and an $n$-vector $b$, along with a rational function $R(z) := D(z )^{-1} N(z)$, we show how to find the optimal approximation to $R(A) b$ from the Krylov space, $\mbox{span}( b, Ab, \ldots , A^{k-1} b)$, using the basis vectors produced by the Arnoldi algorithm. To find this optimal approximation requires running $\max \{ \mbox{deg} (D) , \mbox{deg} (N) \} - 1$ extr… ▽ More Given an $n$ by $n$ matrix $A$ and an $n$-vector $b$, along with a rational function $R(z) := D(z )^{-1} N(z)$, we show how to find the optimal approximation to $R(A) b$ from the Krylov space, $\mbox{span}( b, Ab, \ldots , A^{k-1} b)$, using the basis vectors produced by the Arnoldi algorithm. To find this optimal approximation requires running $\max \{ \mbox{deg} (D) , \mbox{deg} (N) \} - 1$ extra Arnoldi steps and solving a $k + \max \{ \mbox{deg} (D) , \mbox{deg} (N) \}$ by $k$ least squares problem. Here {\em optimal} is taken to mean optimal in the $D(A )^{*} D(A)$-norm. Similar to the case for linear systems, we show that eigenvalues alone cannot provide information about the convergence behavior of this algorithm and we discuss other possible error bounds for highly nonnormal matrices. △ Less

Submitted 29 June, 2023; originally announced June 2023.

MSC Class: 65F10 (Primary); 65F60; 47A12; 47A25 (Secondary)

Evaluating the GeoSnap 13-$μ$m Cut-Off HgCdTe Detector for mid-IR ground-based astronomy

Authors: Jarron M. Leisenring, Dani Atkinson, Rory Bowens, Vincent Douence, William F. Hoffmann, Michael R. Meyer, John Auyeung, James Beletic, Mario S. Cabrera, Alexandra Z. Greenbaum, Phil Hinz, Derek Ives, William J. Forrest, Craig W. McMurtry, Judith L. Pipher, Eric Viges

Abstract: New mid-infrared HgCdTe (MCT) detector arrays developed in collaboration with Teledyne Imaging Sensors (TIS) have paved the way for improved 10-$μ$m sensors for space- and ground-based observatories. Building on the successful development of longwave HAWAII-2RGs for space missions such as NEO Surveyor, we characterize the first 13-$μ$m GeoSnap detector manufactured to overcome the challenges of hi… ▽ More New mid-infrared HgCdTe (MCT) detector arrays developed in collaboration with Teledyne Imaging Sensors (TIS) have paved the way for improved 10-$μ$m sensors for space- and ground-based observatories. Building on the successful development of longwave HAWAII-2RGs for space missions such as NEO Surveyor, we characterize the first 13-$μ$m GeoSnap detector manufactured to overcome the challenges of high background rates inherent in ground-based mid-IR astronomy. This test device merges the longwave HgCdTe photosensitive material with Teledyne's 2048x2048 GeoSnap-18 (18-$μ$m pixel) focal plane module, which is equipped with a capacitive transimpedance amplifier (CTIA) readout circuit paired with an onboard 14-bit analog-to-digital converter (ADC). The final assembly yields a mid-IR detector with high QE, fast readout (>85 Hz), large well depth (>1.2 million electrons), and linear readout. Longwave GeoSnap arrays would ideally be deployed on existing ground-based telescopes as well as the next generation of extremely large telescopes. While employing advanced adaptive optics (AO) along with state-of-the-art diffraction suppression techniques, instruments utilizing these detectors could attain background- and diffraction-limited imaging at inner working angles <10 $λ/D$, providing improved contrast-limited performance compared to JWST MIRI while operating at comparable wavelengths. We describe the performance characteristics of the 13-$μ$m GeoSnap array operating between 38 and 45K, including quantum efficiency, well depth, linearity, gain, dark current, and frequency-dependent (1/f) noise profile. △ Less

Submitted 17 July, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 17 pages, 17 figures. Published in special edition of Astronomische Nachrichten / Astronomical Notes as a contribution to SDW2022

Journal ref: Astron.Nachr./AN, e20230103 (2023)

Detecting Exoplanets Closer to Stars with Moderate Spectral Resolution Integral-Field Spectroscopy

Authors: Shubh Agrawal, Jean-Baptiste Ruffio, Quinn M. Konopacky, Bruce Macintosh, Dimitri Mawet, Eric L. Nielsen, Kielan K. W. Hoch, Michael C. Liu, Travis S. Barman, William Thompson, Alexandra Z. Greenbaum, Christian Marois, Jenny Patience

Abstract: 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 spec… ▽ More 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σ$ 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. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Astronomical Journal on May 12, 2023

JWST/NIRCam discovery of the first Y+Y brown dwarf binary: WISE J033605.05$-$014350.4

Authors: Per Calissendorff, Matthew De Furio, Michael Meyer, Loïc Albert, Christian Aganze, Mohamad Ali-Dib, Daniella C. Bardalez Gagliuffi, Frederique Baron, Charles A. Beichman, Adam J. Burgasser, Michael C. Cushing, Jacqueline Kelly Faherty, Clémence Fontanive, Christopher R. Gelino, John E. Gizis, Alexandra Z. Greenbaum, J. Davy Kirkpatrick, Sandy K. Leggett, Frantz Martinache, David Mary, Mamadou N'Diaye, Benjamin J. S. Pope, Thomas L Roellig, Johannes Sahlmann, Anand Sivaramakrishnan , et al. (3 additional authors not shown)

Abstract: We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in… ▽ More We report the discovery of the first brown dwarf binary system with a Y dwarf primary, WISE J033605.05$-$014350.4, observed with NIRCam on JWST with the F150W and F480M filters. We employed an empirical point spread function binary model to identify the companion, located at a projected separation of 84 milliarcseconds, position angle of 295 degrees, and with contrast of 2.8 and 1.8 magnitudes in F150W and F480M, respectively. At a distance of 10$\,$pc based on its Spitzer parallax, and assuming a random inclination distribution, the physical separation is approximately 1$\,$au. Evolutionary models predict for that an age of 1-5 Gyr, the companion mass is about 4-12.5 Jupiter masses around the 7.5-20 Jupiter mass primary, corresponding to a companion-to-host mass fraction of $q=0.61\pm0.05$. Under the assumption of a Keplerian orbit the period for this extreme binary is in the range of 5-9 years. The system joins a small but growing sample of ultracool dwarf binaries with effective temperatures of a few hundreds of Kelvin. Brown dwarf binaries lie at the nexus of importance for understanding the formation mechanisms of these elusive objects, as they allow us to investigate whether the companions formed as stars or as planets in a disk around the primary. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 8 pages, 3 figures, 1 table. Accepted for publication in Astrophysical Journal Letters

Near-Optimal Approximation of Matrix Functions by the Lanczos Method

Authors: Noah Amsel, Tyler Chen, Anne Greenbaum, Cameron Musco, Chris Musco

Abstract: We study the widely used Lanczos method for approximating the action of a matrix function $f(\mathbf{A})$ on a vector $\mathbf{b}$ (Lanczos-FA). For the function $f(\mathbf{A})=\mathbf{A}^{-1}$, it is known that, when $\mathbf{A}$ is positive definite, the $\mathbf{A}$-norm error of Lanczos-FA after $k$ iterations matches the optimal approximation from the Krylov subspace of degree $k$ generated b… ▽ More We study the widely used Lanczos method for approximating the action of a matrix function $f(\mathbf{A})$ on a vector $\mathbf{b}$ (Lanczos-FA). For the function $f(\mathbf{A})=\mathbf{A}^{-1}$, it is known that, when $\mathbf{A}$ is positive definite, the $\mathbf{A}$-norm error of Lanczos-FA after $k$ iterations matches the optimal approximation from the Krylov subspace of degree $k$ generated by $\mathbf A$ and $\mathbf b$. In this work, we ask whether Lanczos-FA also obtains similarly strong optimality guarantees for other functions $f$. We partially answer this question by showing that, up to a multiplicative approximation factor, Lanczos-FA also matches the optimal approximation from the Krylov subspace for any rational function with real poles outside the interval containing $\mathbf{A}$'s eigenvalues. The approximation factor depends on the degree of $f$'s denominator and the condition number of $\mathbf{A}$, but not on the number of iterations $k$. This result provides theoretical justification for the excellent performance of Lanczos-FA on functions that are well approximated by rationals. Additionally, using different techniques, we prove that Lanczos-FA achieves a weaker notion of optimality for the functions $\mathbf A^{1/2}$ and $\mathbf A^{-1/2}$. Experiments confirm that our new bounds more accurately predict the convergence of Lanczos-FA than existing bounds, in particular those based on uniform polynomial approximation. △ Less

Submitted 12 December, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

MSC Class: 65F60; 65F50; 68Q25

GMRES, pseudospectra, and Crouzeix's conjecture for shifted and scaled Ginibre matrices

Authors: Tyler Chen, Anne Greenbaum, Thomas Trogdon

Abstract: We study the GMRES algorithm applied to linear systems of equations involving a scaled and shifted $N\times N$ matrix whose entries are independent complex Gaussians. When the right hand side of this linear system is independent of this random matrix, the $N\to\infty$ behavior of the GMRES residual error can be determined exactly. To handle cases where the right hand side depends on the random mat… ▽ More We study the GMRES algorithm applied to linear systems of equations involving a scaled and shifted $N\times N$ matrix whose entries are independent complex Gaussians. When the right hand side of this linear system is independent of this random matrix, the $N\to\infty$ behavior of the GMRES residual error can be determined exactly. To handle cases where the right hand side depends on the random matrix, we study the pseudospectra and numerical range of Ginibre matrices and prove a restricted version of Crouzeix's conjecture. △ Less

Submitted 3 March, 2023; originally announced March 2023.

MSC Class: 68Q25; 65F35; 15A60

JWST Observations of the Enigmatic Y Dwarf WISE 1828+2650: I. Limits to a Binary Companion

Authors: Matthew De Furio, Ben W. Lew, Charles A. Beichman, Thomas Roellig, Geoffrey Bryden, David R. Ciardi, Michael R. Meyer, Marcia J. Rieke, Alexandra Z. Greenbaum, Jarron Leisenring, Jorge Llop-Sayson, Marie Ygouf, Loïc Albert, Martha L. Boyer, Daniel J. Eisenstein, Klaus W. Hodapp, Scott Horner, Doug Johnstone, Douglas M. Kelly, Karl A. Misselt, George H. Rieke, John A. Stansberry, Erick T. Young

Abstract: The Y-dwarf WISE 1828+2650 is one of the coldest known Brown Dwarfs with an effective temperature of $\sim$300 K. Located at a distance of just 10 pc, previous model-based estimates suggest WISE1828+2650 has a mass of $\sim$5-10 Mj, making it a valuable laboratory for understanding the formation, evolution and physical characteristics of gas giant planets. However, previous photometry and spectros… ▽ More The Y-dwarf WISE 1828+2650 is one of the coldest known Brown Dwarfs with an effective temperature of $\sim$300 K. Located at a distance of just 10 pc, previous model-based estimates suggest WISE1828+2650 has a mass of $\sim$5-10 Mj, making it a valuable laboratory for understanding the formation, evolution and physical characteristics of gas giant planets. However, previous photometry and spectroscopy have presented a puzzle with the near-impossibility of simultaneously fitting both the short (0.9-2.0 microns) and long wavelength (3-5 microns) data. A potential solution to this problem has been the suggestion that WISE 1828+2650 is a binary system whose composite spectrum might provide a better match to the data. Alternatively, new models being developed to fit JWST/NIRSpec and MIRI spectroscopy might provide new insights. This article describes JWST/NIRCam observations of WISE 1828+2650 in 6 filters to address the binarity question and to provide new photometry to be used in model fitting. We also report Adaptive Optics imaging with the Keck 10 m telescope. We find no evidence for multiplicity for a companion beyond 0.5 AU with either JWST or Keck. Companion articles will present low and high resolution spectra of WISE 1828+2650 obtained with both NIRSpec and MIRI. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 15 pages, 9 figures, Accepted by ApJ on Feb. 21 2023

K-Spectral Sets

Authors: Anne Greenbaum, Natalie Wellen

Abstract: We use results in [M. Crouzeix and A. Greenbaum,Spectral sets: numerical range and beyond, SIAM Jour. Matrix Anal. Appl., 40 (2019), pp. 1087-1101] to derive a variety of K-spectral sets and show how they can be used in some applications. We compare the K-values derived here to those that can be derived from a straightforward application of the Cauchy integral formula, by replacing the norm of the… ▽ More We use results in [M. Crouzeix and A. Greenbaum,Spectral sets: numerical range and beyond, SIAM Jour. Matrix Anal. Appl., 40 (2019), pp. 1087-1101] to derive a variety of K-spectral sets and show how they can be used in some applications. We compare the K-values derived here to those that can be derived from a straightforward application of the Cauchy integral formula, by replacing the norm of the integral by the integral of the resolvent norm. While, in some cases, the new upper bounds on the optimal K-value are much tighter than those from the Cauchy integral formula, we show that in many cases of interest, the two values are of the same order of magnitude, with the bounds from the Cauchy integral formula actually being slightly smaller. We give a partial explanation of this in terms of the numerical range of the resolvent at points near an ill-conditioned eigenvalue. △ Less

Submitted 2 November, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 20 pages, 10 figures. For relevant code, see https://github.com/tygris/k-spectral-sets. Abstract refers to arXiv:1803.10904

MSC Class: 47A25; 47A30

Testing the Interaction Between a Substellar Companion and a Debris Disk in the HR 2562 System

Authors: Stella Yimiao Zhang, Gaspard Duchêne, Robert J. De Rosa, Megan Ansdell, Quinn Konopacky, Thomas Esposito, Eugene Chiang, Malena Rice, Brenda Matthews, Paul Kalas, Bruce Macintosh, Franck Marchis, Stan Metchev, Jenny Patience, Julien Rameau, Kimberly Ward-Duong, Schuyler Wolff, Michael P. Fitzgerald, Vanessa P. Bailey, Travis S. Barman, Joanna Bulger, Christine H. Chen, Jeffrey K. Chilcotte, Tara Cotten, René Doyon , et al. (29 additional authors not shown)

Abstract: 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… ▽ More 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$σ$ 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. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 22 pages, 17 figures

First Observations of the Brown Dwarf HD 19467 B with JWST

Authors: Alexandra Z. Greenbaum, Jorge Llop-Sayson, Ben Lew, Geoffrey Bryden, Thomas Roellig, Marie Ygouf, B. J. Fulton, Daniel R. Hey, Daniel Huber, Sagnick Mukherjee, Michael Meyer, Jarron Leisenring, Marcia Rieke, Martha Boyer, Joseph J. Green, Doug Kelly, Karl Misselt, Eugene Serabyn, John Stansberry, Laurie E. U. Chu, Matthew De Furio, Doug Johnstone, Joshua E. Schlieder, Charles Beichman

Abstract: We observed HD 19467 B with JWST's NIRCam in six filters spanning 2.5-4.6 $μm$ with the Long Wavelength Bar coronagraph. The brown dwarf HD 19467 B was initially identified through a long-period trend in the radial velocity of G3V star HD 19467. HD 19467 B was subsequently detected via coronagraphic imaging and spectroscopy, and characterized as a late-T type brown dwarf with approximate temperatu… ▽ More We observed HD 19467 B with JWST's NIRCam in six filters spanning 2.5-4.6 $μm$ with the Long Wavelength Bar coronagraph. The brown dwarf HD 19467 B was initially identified through a long-period trend in the radial velocity of G3V star HD 19467. HD 19467 B was subsequently detected via coronagraphic imaging and spectroscopy, and characterized as a late-T type brown dwarf with approximate temperature $\sim1000$K. We observed HD 19467 B as a part of the NIRCam GTO science program, demonstrating the first use of the NIRCam Long Wavelength Bar coronagraphic mask. The object was detected in all 6 filters (contrast levels of $2\times10^{-4}$ to $2\times10^{-5}$) at a separation of 1.6 arcsec using Angular Differential Imaging (ADI) and Synthetic Reference Differential Imaging (SynRDI). Due to a guidestar failure during acquisition of a pre-selected reference star, no reference star data was available for post-processing. However, RDI was successfully applied using synthetic Point Spread Functions (PSFs) developed from contemporaneous maps of the telescope's optical configuration. Additional radial velocity data (from Keck/HIRES) are used to constrain the orbit of HD 19467 B. Photometric data from TESS are used to constrain the properties of the host star, particularly its age. NIRCam photometry, spectra and photometry from literature, and improved stellar parameters are used in conjunction with recent spectral and evolutionary substellar models to derive physical properties for HD 19467 B. Using an age of 9.4$\pm$0.9 Gyr inferred from spectroscopy, Gaia astrometry, and TESS asteroseismology, we obtain a model-derived mass of 62$\pm 1M_{J}$, which is consistent within 2-$σ$ with the dynamically derived mass of 81$^{+14}_{-12}M_{J}$. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 21 pages, 19 figures. Accepted to AAS Journals

The Near Infrared Imager and Slitless Spectrograph for JWST -- V. Kernel Phase Imaging and Data Analysis

Authors: Jens Kammerer, Rachel A. Cooper, Thomas Vandal, Deepashri Thatte, Frantz Martinache, Anand Sivaramakrishnan, Alexander Chaushev, Tomas Stolker, James P. Lloyd, Loïc Albert, René Doyon, Steph Sallum, Marshall D. Perrin, Laurent Pueyo, Antoine Mérand, Alexandre Gallenne, Alexandra Greenbaum, Joel Sanchez-Bermudez, Dori Blakely, Doug Johnstone, Kevin Volk, Andre Martel, Paul Goudfrooij, Michael R. Meyer, Chris J. Willott , et al. (4 additional authors not shown)

Abstract: Kernel phase imaging (KPI) enables the direct detection of substellar companions and circumstellar dust close to and below the classical (Rayleigh) diffraction limit. We present a kernel phase analysis of JWST NIRISS full pupil images taken during the instrument commissioning and compare the performance to closely related NIRISS aperture masking interferometry (AMI) observations. For this purpose,… ▽ More Kernel phase imaging (KPI) enables the direct detection of substellar companions and circumstellar dust close to and below the classical (Rayleigh) diffraction limit. We present a kernel phase analysis of JWST NIRISS full pupil images taken during the instrument commissioning and compare the performance to closely related NIRISS aperture masking interferometry (AMI) observations. For this purpose, we develop and make publicly available the custom "Kpi3Pipeline" enabling the extraction of kernel phase observables from JWST images. The extracted observables are saved into a new and versatile kernel phase FITS file (KPFITS) data exchange format. Furthermore, we present our new and publicly available "fouriever" toolkit which can be used to search for companions and derive detection limits from KPI, AMI, and long-baseline interferometry observations while accounting for correlated uncertainties in the model fitting process. Among the four KPI targets that were observed during NIRISS instrument commissioning, we discover a low-contrast (~1:5) close-in (~1 $λ/D$) companion candidate around CPD-66~562 and a new high-contrast (~1:170) detection separated by ~1.5 $λ/D$ from 2MASS~J062802.01-663738.0. The 5-$σ$ companion detection limits around the other two targets reach ~6.5 mag at ~200 mas and ~7 mag at ~400 mas. Comparing these limits to those obtained from the NIRISS AMI commissioning observations, we find that KPI and AMI perform similar in the same amount of observing time. Due to its 5.6 times higher throughput if compared to AMI, KPI is beneficial for observing faint targets and superior to AMI at separations >325 mas. At very small separations (<100 mas) and between ~250-325 mas, AMI slightly outperforms KPI which suffers from increased photon noise from the core and the first Airy ring of the point-spread function. △ Less

Submitted 3 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 34 pages, 17 figures, accepted for publication in PASP

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- IV. Aperture Masking Interferometry

Authors: Anand Sivaramakrishnan, Peter Tuthill, James P. Lloyd, Alexandra Z. Greenbaum, Deepashri Thatte, Rachel A. Cooper, Thomas Vandal, Jens Kammerer, Joel Sanchez-Bermudez, Benjamin J. S. Pope, Dori Blakely, Loïc Albert, Neil J. Cook, Doug Johnstone, André R. Martel, Kevin Volk, Anthony Soulain, Étienne Artigau, David Lafrenière, Chris J. Willott, Sébastien Parmentier, K. E. Saavik Ford, Barry McKernan, M. Begoña Vila, Neil Rowlands , et al. (14 additional authors not shown)

Abstract: The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, pres… ▽ More The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, present some sample science cases, explain its operational observing strategies, indicate how AMI proposals can be developed with data simulations, and how AMI data can be analyzed. We also present key results from commissioning AMI. Since the allied Kernel Phase Imaging (KPI) technique benefits from AMI operational strategies, we also cover NIRISS KPI methods and analysis techniques, including a new user-friendly KPI pipeline. The NIRISS KPI bright limit is $\simeq 8$ W2 magnitudes. AMI (and KPI) achieve an inner working angle of $\sim 70$ mas that is well inside the $\sim 400$ mas NIRCam inner working angle for its circular occulter coronagraphs at comparable wavelengths. △ Less

Submitted 7 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 30 pages, 10 figures

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems II: A 1 to 20 Micron Spectrum of the Planetary-Mass Companion VHS 1256-1257 b

Authors: Brittany E. Miles, Beth A. Biller, Polychronis Patapis, Kadin Worthen, Emily Rickman, Kielan K. W. Hoch, Andrew Skemer, Marshall D. Perrin, Niall Whiteford, Christine H. Chen, B. Sargent, Sagnick Mukherjee, Caroline V. Morley, Sarah E. Moran, Mickael Bonnefoy, Simon Petrus, Aarynn L. Carter, Elodie Choquet, Sasha Hinkley, Kimberly Ward-Duong, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Laurent Pueyo, Shrishmoy Ray, Karl R. Stapelfeldt , et al. (79 additional authors not shown)

Abstract: 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… ▽ More 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 $μ$m to 20 $μ$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. △ Less

Submitted 4 July, 2024; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 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

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems I: High Contrast Imaging of the Exoplanet HIP 65426 b from 2-16 $μ$m

Authors: Aarynn L. Carter, Sasha Hinkley, Jens Kammerer, Andrew Skemer, Beth A. Biller, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Simon Petrus, Jordan M. Stone, Kimberly Ward-Duong, Jason J. Wang, Julien H. Girard, Dean C. Hines, Marshall D. Perrin, Laurent Pueyo, William O. Balmer, Mariangela Bonavita, Mickael Bonnefoy, Gael Chauvin, Elodie Choquet, Valentin Christiaens, Camilla Danielski, Grant M. Kennedy, Elisabeth C. Matthews, Brittany E. Miles , et al. (86 additional authors not shown)

Abstract: 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 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$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 exo… ▽ More 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 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$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 $μ$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$σ$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $μ$m and MIRI at 11.3 $μ$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 $μ$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. △ Less

Submitted 3 May, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

Atmospheric Monitoring and Precise Spectroscopy of the HR 8799 Planets with SCExAO/CHARIS

Authors: Jason J. Wang, Peter Gao, Jeffrey Chilcote, Julien Lozi, Olivier Guyon, Christian Marois, Robert J. De Rosa, Ananya Sahoo, Tyler D. Groff, Sebastien Vievard, Nemanja Jovanovic, Alexandra Z. Greenbaum, Bruce Macintosh

Abstract: 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… ▽ More 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. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 21 pages, 12 figures, Accepted to AJ, time series photometry and stacked spectra available as text files

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

Authors: Evan A. Rich, John D. Monnier, Alicia Aarnio, Anna S. E. Laws, Benjamin R. Setterholm, David J. Wilner, Nuria Calvet, Tim Harries, Chris Miller, Claire L. Davies, Fred C. Adams, Sean M. Andrews, Jaehan Bae, Catherine Espaillat, Alexandra Z. Greenbaum, Sasha Hinkley, Stefan Kraus, Lee Hartmann, Andrea Isella, Melissa McClure, Rebecca Oppenheimer, Laura M. Pérez, Zhaohuan Zhu

Abstract: We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of… ▽ More We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 51 pages, 31 figures, 7 tables, accepted to AJ

The JWST Early Release Science Program for the Direct Imaging & Spectroscopy of Exoplanetary Systems

Authors: Sasha Hinkley, Aarynn L. Carter, Shrishmoy Ray, Andrew Skemer, Beth Biller, Elodie Choquet, Maxwell A. Millar-Blanchaer, Stephanie Sallum, Brittany Miles, Niall Whiteford, Polychronis Patapis, Marshall D. Perrin, Laurent Pueyo, Glenn Schneider, Karl Stapelfeldt, Jason Wang, Kimberly Ward-Duong, Brendan P. Bowler, Anthony Boccaletti, Julien H. Girard, Dean Hines, Paul Kalas, Jens Kammerer, Pierre Kervella, Jarron Leisenring , et al. (61 additional authors not shown)

Abstract: The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe e… ▽ More The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$μ$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$μ$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles. △ Less

Submitted 12 September, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 21 pages, 7 figures. Accepted for Publication in PASP

Low-memory Krylov subspace methods for optimal rational matrix function approximation

Authors: Tyler Chen, Anne Greenbaum, Cameron Musco, Christopher Musco

Abstract: We describe a Lanczos-based algorithm for approximating the product of a rational matrix function with a vector. This algorithm, which we call the Lanczos method for optimal rational matrix function approximation (Lanczos-OR), returns the optimal approximation from a given Krylov subspace in a norm depending on the rational function's denominator, and can be computed using the information from a s… ▽ More We describe a Lanczos-based algorithm for approximating the product of a rational matrix function with a vector. This algorithm, which we call the Lanczos method for optimal rational matrix function approximation (Lanczos-OR), returns the optimal approximation from a given Krylov subspace in a norm depending on the rational function's denominator, and can be computed using the information from a slightly larger Krylov subspace. We also provide a low-memory implementation which only requires storing a number of vectors proportional to the denominator degree of the rational function. Finally, we show that Lanczos-OR can be used to derive algorithms for computing other matrix functions, including the matrix sign function and quadrature based rational function approximations. In many cases, it improves on the approximation quality of prior approaches, including the standard Lanczos method, with little additional computational overhead. △ Less

Submitted 31 May, 2023; v1 submitted 22 February, 2022; originally announced February 2022.

Journal ref: SIAM Journal on Matrix Analysis and Applications 2023 44:2, 670-692

Point Source Localization with a Planar Optical Phased Array Compressive Sensor

Authors: Julian A. Brown, Steven J. Spector, Michael Moebius, Lucas Benney, Daniel Vresilovic, Brian Dolle, Alexandra Z. Greenbaum, Alex Huang, Christopher V. Poulton, Michael R. Watts, Robin Dawson, Benjamin F. Lane, J. P. Laine, Kerri Cahoy, Hannah A. Clevenson

Abstract: Compressive sensing has been used to demonstrate scene reconstruction and source localization in a wide variety of devices. To date, optical compressive sensors have not been able to achieve significant volume reduction relative to conventional optics of equivalent angular resolution. Here, we adapt silicon-photonic optical phased array technology to demonstrate, to our knowledge, the first applic… ▽ More Compressive sensing has been used to demonstrate scene reconstruction and source localization in a wide variety of devices. To date, optical compressive sensors have not been able to achieve significant volume reduction relative to conventional optics of equivalent angular resolution. Here, we adapt silicon-photonic optical phased array technology to demonstrate, to our knowledge, the first application of compressive imaging in a photonic-integrated device. Our novel sensor consists of an $8\times 8$ grid of grating couplers with a spacing of $100~μ$m. Path-matched waveguides route to a single multimode interferometer (MMI), which mixes and randomizes the signals into 64 outputs to be used for compressed sensing. Our device is fully passive, having no need for phase shifters, as measurement matrix calibration makes the measurements robust to phase errors. For testing, we use an Amplified Spontaneous Emission (ASE) source with a bandwidth of 40 nm, centered at 1545 nm. We demonstrate simultaneous multi-point (2 sources demonstrated in this work) brightness recovery and localization with better than 10 arcsecond precision in a sub-millimeter thick form-factor. We achieve a single source recovery rate higher than 99.9\% using 10 of the 64 outputs, and a 90\% recovery rate with only 6 outputs, 10 times fewer than the 64 needed for conventional imaging. This planar optical phased array compressive sensor is well-suited for imaging sparse scenes in applications constrained by form factor, volume, or high-cost detectors, with the potential to revolutionize endoscopy, beam locators, and LIDAR. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 10 pages, 6 figures, pre-print to be published in SPIE Optical Engineering and Applications

Error bounds for Lanczos-based matrix function approximation

Authors: Tyler Chen, Anne Greenbaum, Cameron Musco, Christopher Musco

Abstract: We analyze the Lanczos method for matrix function approximation (Lanczos-FA), an iterative algorithm for computing $f(\mathbf{A}) \mathbf{b}$ when $\mathbf{A}$ is a Hermitian matrix and $\mathbf{b}$ is a given vector. Assuming that $f : \mathbb{C} \rightarrow \mathbb{C}$ is piecewise analytic, we give a framework, based on the Cauchy integral formula, which can be used to derive a priori and a pos… ▽ More We analyze the Lanczos method for matrix function approximation (Lanczos-FA), an iterative algorithm for computing $f(\mathbf{A}) \mathbf{b}$ when $\mathbf{A}$ is a Hermitian matrix and $\mathbf{b}$ is a given vector. Assuming that $f : \mathbb{C} \rightarrow \mathbb{C}$ is piecewise analytic, we give a framework, based on the Cauchy integral formula, which can be used to derive a priori and a posteriori error bounds for Lanczos-FA in terms of the error of Lanczos used to solve linear systems. Unlike many error bounds for Lanczos-FA, these bounds account for fine-grained properties of the spectrum of $\mathbf{A}$, such as clustered or isolated eigenvalues. Our results are derived assuming exact arithmetic, but we show that they are easily extended to finite precision computations using existing theory about the Lanczos algorithm in finite precision. We also provide generalized bounds for the Lanczos method used to approximate quadratic forms $\mathbf{b}^\textsf{H} f(\mathbf{A}) \mathbf{b}$, and demonstrate the effectiveness of our bounds with numerical experiments. △ Less

Submitted 18 May, 2022; v1 submitted 17 June, 2021; originally announced June 2021.

Journal ref: SIAM Journal on Matrix Analysis and Applications, 2022, Vol. 43, No. 2 : pp. A3084-A3108

Positive Pressure Testing Booths Development and Deployment In Response To The COVID-19 Outbreak

Authors: Kevin Aroom, Jiawei Ge, Lidia Al-Zogbi, Marcee White, Adrienne Trustman, Adena Greenbaum, Jason Farley, Axel Krieger

Abstract: The COVID-19 pandemic left an unprecedented impact on the general public health, resulting in thousands of deaths in the US alone. Nationwide testing plans were initiated to control the spread, with drive-through being the currently dominant testing approach, which, however, exhausts personal protective equipment supplies, and is unfriendly to individuals not owning a vehicle. Walk-up testing boot… ▽ More The COVID-19 pandemic left an unprecedented impact on the general public health, resulting in thousands of deaths in the US alone. Nationwide testing plans were initiated to control the spread, with drive-through being the currently dominant testing approach, which, however, exhausts personal protective equipment supplies, and is unfriendly to individuals not owning a vehicle. Walk-up testing booths are a safe alternative, but are too prohibitively priced on the market to allow for nationwide deployment. In this paper, we present an accessible, mobile, affordable, and safe version of a positive-pressure COVID-19 testing booth. The booths are manufactured using primarily off-the-shelf components from US vendors with minimized customization. The booths' mobility allows them to be easily transported within local communities to test a larger subset of the population with fewer transportation options. Moreover, the final bill of materials does not surpass USD 3,900, which is about half of the market price. The booths are air conditioned and HEPA filtered to offer healthcare providers a safe and comfortable working environment. The prototype passed required pressure and air exchange tests, and was positively reviewed by two healthcare professionals. Currently, five booths are deployed and used at the Johns Hopkins University School of Nursing, Baltimore City Health Department, and two community health centers in Baltimore. Our design facilitates walk-up testing in the US, as it decreases PPE consumption; reduces the risk of infection; and is accessible to lower-income communities and non-drivers. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Gemini Planet Imager Spectroscopy of the Dusty Substellar Companion HD 206893 B

Authors: K. Ward-Duong, J. Patience, K. Follette, R. J. De Rosa, J. Rameau, M. Marley, D. Saumon, E. L. Nielsen, A. Rajan, A. Z. Greenbaum, J. Lee, J. J. Wang, I. Czekala, G. Duchêne, B. Macintosh, S. Mark Ammons, V. P. Bailey, T. Barman, J. Bulger, C. Chen, J. Chilcote, T. Cotten, R. Doyon, T. M. Esposito, M. P. Fitzgerald , et al. (33 additional authors not shown)

Abstract: We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challengi… ▽ More We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 K to 1800 K for effective temperature and from 3.0 to 5.0 for log($g$), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of sub-micron dust particles, similar to dereddening approaches applied to the peculiar red field L-dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L-dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from four years of astrometric monitoring is consistent with a $\sim$30-year period, orbital inclination of 147$^{\circ}$, and semimajor axis of 10 au, well within the estimated disk inner radius of $\sim$50 au. As one of very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion-disk interaction and provide a benchmark for substellar and planetary atmospheric study. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Accepted for publication in AJ (October 15, 2020). 35 pages, 20 figures, with tables and appendices presented in their entirety

Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e

Authors: P. Mollière, T. Stolker, S. Lacour, G. P. P. L. Otten, J. Shangguan, B. Charnay, T. Molyarova, M. Nowak, Th. Henning, G. -D. Marleau, D. A. Semenov, E. van Dishoeck, F. Eisenhauer, P. Garcia, R. Garcia Lopez, J. H. Girard, A. Z. Greenbaum, S. Hinkley, P. Kervella, L. Kreidberg, A. -L. Maire, E. Nasedkin, L. Pueyo, I. A. G. Snellen, A. Vigan , et al. (3 additional authors not shown)

Abstract: Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our r… ▽ More Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets. △ Less

Submitted 18 June, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: Accepted for publication in A&A, 28 pages, 13 figures, updated author list

Journal ref: A&A 640, A131 (2020)

Multiband Polarimetric Imaging of HR 4796A with the Gemini Planet Imager

Authors: Pauline Arriaga, Michael P. Fitzgerald, Gaspard Duchêne, Paul Kalas, Maxwell A. Millar-Blanchaer, Marshall D. Perrin, Christine H. Chen, Johan Mazoyer, Mark Ammons, Vanessa P. Bailey, Trafis S. Barman, Joanna Bulger, Jeffrey K. Chilcote, Tara Cotten, Robert J. De Rosa, Rene Doyon, Thomas M. Esposito, Katherine B. Follette, Benjamin L. Gerard, Stephen Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Justin Hom, Li-Wei Hung , et al. (27 additional authors not shown)

Abstract: 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 brightne… ▽ More 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. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: Accepted for publication in AJ. 3 tables, 11 figures

Crouzeix's Conjecture and related problems

Authors: Kelly Bickel, Pamela Gorkin, Anne Greenbaum, Thomas Ransford, Felix Schwenninger, Elias Wegert

Abstract: In this paper, we establish several results related to Crouzeix's conjecture. We show that the conjecture holds for contractions with eigenvalues that are sufficiently well-separated. This separation is measured by the so-called separation constant, which is defined in terms of the pseudohyperbolic metric. Moreover, we study general properties of related extremal functions and associated vectors.… ▽ More In this paper, we establish several results related to Crouzeix's conjecture. We show that the conjecture holds for contractions with eigenvalues that are sufficiently well-separated. This separation is measured by the so-called separation constant, which is defined in terms of the pseudohyperbolic metric. Moreover, we study general properties of related extremal functions and associated vectors. Throughout, compressions of the shift serve as illustrating examples which also allow for refined results. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 24 pages

MSC Class: 47A25; 47A12; 15A60; 47A20; 30J10

Journal ref: Comput. Methods Funct. Theory, 2020

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

Authors: Thomas M. Esposito, Paul Kalas, Michael P. Fitzgerald, Maxwell A. Millar-Blanchaer, Gaspard Duchene, Jennifer Patience, Justin Hom, Marshall D. Perrin, Robert J. De Rosa, Eugene Chiang, Ian Czekala, Bruce Macintosh, James R. Graham, Megan Ansdell, Pauline Arriaga, Sebastian Bruzzone, Joanna Bulger, Christine H. Chen, Tara Cotten, Ruobing Dong, Zachary H. Draper, Katherine B. Follette, Li-Wei Hung, Ronald Lopez, Brenda C. Matthews , et al. (40 additional authors not shown)

Abstract: 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. Observ… ▽ More 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 $μ$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 °$). 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. △ Less

Submitted 23 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: arXiv resubmission with typographical corrections. Accepted for publication in AJ. 19 figures, 7 tables

The Gemini Planet Imager view of the HD 32297 debris disk

Authors: Gaspard Duchene, Malena Rice, Justin Hom, Joseph Zalesky, Thomas M. Esposito, Maxwell A. Millar-Blanchaer, Bin Ren, Paul Kalas, Michael Fitzgerald, Pauline Arriaga, Sebastian Bruzzone, Joanna Bulger, Christine H. Chen, Eugene Chiang, Tara Cotten, Ian Czekala, Robert J. De Rosa, Ruobing Dong, Zachary H. Draper, Katherine B. Follette, James R. Graham, Li-Wei Hung, Ronald Lopez, Bruce Macintosh, Brenda C. Matthews , et al. (38 additional authors not shown)

Abstract: We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We ana… ▽ More We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: Accepted for publication in the Astronomical Journal

HD 165054: an astrometric calibration field for high-contrast imagers in Baade's Window

Authors: Meiji M. Nguyen, Robert J. De Rosa, Jason J. Wang, Thomas M. Esposito, Paul Kalas, James R. Graham, Bruce Macintosh, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Gaspard Duchêne, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, Alexandra Z. Greenbaum, Pascale Hibon, Justin Hom, Li-Wei Hung, Patrick Ingraham, Quinn Konopacky, James E. Larkin , et al. (29 additional authors not shown)

Abstract: We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square… ▽ More We present a study of the HD 165054 astrometric calibration field that has been periodically observed with the Gemini Planet Imager. HD 165054 is a bright star within Baade's Window, a region of the galactic plane with relatively low extinction from interstellar dust. HD 165054 was selected as a calibrator target due to the high number density of stars within this region ($\sim 3$ stars per square arcsecond with $H<22$), necessary because of the small field-of-view of the Gemini Planet Imager. Using nine epochs spanning over five years, we have fit a standard five-parameter astrometric model to the astrometry of seven background stars within close proximity to HD 165054 (angular separation $< 2$ arcsec). We achieved a proper motion precision of $\sim 0.3$ mas/yr, and constrained the parallax of each star to be $\lesssim 1$ mas. Our measured proper motions and parallax limits are consistent with the background stars being a part of the galactic bulge. Using these measurements we find no evidence of any systematic trend of either the plate scale or the north angle offset of GPI between 2014 and 2019. We compared our model describing the motions of the seven background stars to observations of the same field in 2014 and 2018 obtained with Keck/NIRC2, an instrument with an excellent astrometric calibration. We find that predicted position of the background sources is consistent with that measured by NIRC2, within the uncertainties of the calibration of the two instruments. In the future, we will use this field as a standard astrometric calibrator for the upgrade of GPI and potentially for other high-contrast imagers. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 25 pages, 12 figures

Imaging the 44 AU Kuiper Belt-analogue debris ring around HD 141569A with GPI polarimetry

Authors: J. S. Bruzzone, S. Metchev, G. Duchene, M. A. Millar-Blanchaer, R. Dong, J. J. Wang, J. R. Graham, J. Mazoyer, S. Wolff, S. M. Ammons, A. C. Schneider, A. Z. Greenbaum, B. C. Matthews, P. Arriaga, V. P. Bailey, T. Barman, J. Bulger, J. Chilcote, T. Cotten, R. J. De Rosa, R. Doyon, M. P. Fitzgerald, K. B. Follette, B. L. Gerard, S. J. Goodsell , et al. (31 additional authors not shown)

Abstract: We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east… ▽ More We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: 16 pages, 8 figures, accepted at AJ

The Gemini Planet Imager Exoplanet Survey: Dynamical Mass of the Exoplanet beta Pictoris b from Combined Direct Imaging and Astrometry

Authors: Eric L. Nielsen, Robert J. De Rosa, Jason J. Wang, Johannes Sahlmann, Paul Kalas, Gaspard Duchene, Julien Rameau, Mark S. Marley, Didier Saumon, Bruce Macintosh, Maxwell A. Millar-Blanchaer, Meiji M. Nguyen, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham , et al. (29 additional authors not shown)

Abstract: We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit… ▽ More We present new observations of the planet beta Pictoris b from 2018 with GPI, the first GPI observations following conjunction. Based on these new measurements, we perform a joint orbit fit to the available relative astrometry from ground-based imaging, the Hipparcos Intermediate Astrometric Data (IAD), and the Gaia DR2 position, and demonstrate how to incorporate the IAD into direct imaging orbit fits. We find a mass consistent with predictions of hot-start evolutionary models and previous works following similar methods, though with larger uncertainties: 12.8 [+5.3, -3.2] M_Jup. Our eccentricity determination of 0.12 [+0.04, -0.03] disfavors circular orbits. We consider orbit fits to several different imaging datasets, and find generally similar posteriors on the mass for each combination of imaging data. Our analysis underscores the importance of performing joint fits to the absolute and relative astrometry simultaneously, given the strong covariance between orbital elements. Time of conjunction is well constrained within 2.8 days of 2017 September 13, with the star behind the planet's Hill sphere between 2017 April 11 and 2018 February 16 (+/- 18 days). Following the recent radial velocity detection of a second planet in the system, beta Pic c, we perform additional two-planet fits combining relative astrometry, absolute astrometry, and stellar radial velocities. These joint fits find a significantly smaller mass for the imaged planet beta Pic b, of 8.0 +/- 2.6 M_Jup, in a somewhat more circular orbit. We expect future ground-based observations to further constrain the visual orbit and mass of the planet in advance of the release of Gaia DR4. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: 29 Pages, 23 Figures. Accepted to AJ

First Resolved Scattered-Light Images of Four Debris Disks in Scorpius-Centaurus with the Gemini Planet Imager

Authors: Justin Hom, Jennifer Patience, Thomas M. Esposito, Gaspard Duchêne, Kadin Worthen, Paul Kalas, Hannah Jang-Condell, Kezman Saboi, Pauline Arriaga, Johan Mazoyer, Schuyler Wolff, Maxwell A. Millar-Blanchaer, Michael P. Fitzgerald, Marshall D. Perrin, Christine H. Chen, Bruce Macintosh, Brenda C. Matthews, Jason J. Wang, James R. Graham, Franck Marchis, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey K. Chilcote , et al. (35 additional authors not shown)

Abstract: We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, a… ▽ More We present the first spatially resolved scattered-light images of four debris disks around members of the Scorpius-Centaurus (Sco-Cen) OB Association with high-contrast imaging and polarimetry using the Gemini Planet Imager (GPI). All four disks are resolved for the first time in polarized light and one disk is also detected in total intensity. The three disks imaged around HD 111161, HD 143675, and HD 145560 are symmetric in both morphology and brightness distribution. The three systems span a range of inclinations and radial extents. The disk imaged around HD 98363 shows indications of asymmetries in morphology and brightness distribution, with some structural similarities to the HD 106906 planet-disk system. Uniquely, HD 98363 has a wide co-moving stellar companion Wray 15-788 with a recently resolved disk with very different morphological properties. HD 98363 A/B is the first binary debris disk system with two spatially resolved disks. All four targets have been observed with ALMA, and their continuum fluxes range from one non-detection to one of the brightest disks in the region. With the new results, a total of 15 A/F-stars in Sco-Cen have resolved scattered light debris disks, and approximately half of these systems exhibit some form of asymmetry. Combining the GPI disk structure results with information from the literature on millimeter fluxes and imaged planets reveals a diversity of disk properties in this young population. Overall, the four newly resolved disks contribute to the census of disk structures measured around A/F-stars at this important stage in the development of planetary systems. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: 20 pages, 9 figures

Detection of a low-mass stellar companion to the accelerating A2IV star HR 1645

Authors: Robert J. De Rosa, Eric L. Nielsen, Julien Rameau, Gaspard Duchêne, Alexandra Z. Greenbaum, Jason J. Wang, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Pascale Hibon, Justin Hom, Li-Wei Hung, Patrick Ingraham, Paul Kalas Quinn Konopacky, James E. Larkin , et al. (28 additional authors not shown)

Abstract: The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to… ▽ More The $\sim500$\, Myr A2IV star HR 1645 has one of the most significant low-amplitude accelerations of nearby early-type stars measured from a comparison of the {\it Hipparcos} and {\it Gaia} astrometric catalogues. This signal is consistent with either a stellar companion with a moderate mass ratio ($q\sim0.5$) on a short period ($P<1$\,yr), or a substellar companion at a separation wide enough to be resolved with ground-based high contrast imaging instruments; long-period equal mass ratio stellar companions that are also consistent with the measured acceleration are excluded with previous imaging observations. The small but significant amplitude of the acceleration made HR 1645 a promising candidate for targeted searches for brown dwarf and planetary-mass companions around nearby, young stars. In this paper we explore the origin of the astrometric acceleration by modelling the signal induced by a wide-orbit M8 companion discovered with the Gemini Planet Imager, as well as the effects of an inner short-period spectroscopic companion discovered a century ago but not since followed-up. We present the first constraints on the orbit of the inner companion, and demonstrate that it is a plausible cause of the astrometric acceleration. This result demonstrates the importance of vetting of targets with measured astrometric acceleration for short-period stellar companions prior to conducting targeted direct imaging surveys for wide-orbit substellar companions. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 14 pages, 9 figures. Accepted for publication in the Astronomical Journal

An updated visual orbit of the directly-imaged exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia

Authors: Robert J. De Rosa, Eric L. Nielsen, Jason J. Wang, S. Mark Ammons, Gaspard Duchêne, Bruce Macintosh, Meiji M. Nguyen, Julien Rameau, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Justin Hom, Li-Wei Hung, Patrick Ingraham , et al. (30 additional authors not shown)

Abstract: We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional… ▽ More We present a revision to the visual orbit of the young, directly-imaged exoplanet 51 Eridani b using four years of observations with the Gemini Planet Imager. The relative astrometry is consistent with an eccentric ($e=0.53_{-0.13}^{+0.09}$) orbit at an intermediate inclination ($i=136_{-11}^{+10}$\,deg), although circular orbits cannot be excluded due to the complex shape of the multidimensional posterior distribution. We find a semi-major axis of $11.1_{-1.3}^{+4.2}$\,au and a period of $28.1_{-4.9}^{+17.2}$\,yr, assuming a mass of 1.75\,M$_{\odot}$ for the host star. We find consistent values with a recent analysis of VLT/SPHERE data covering a similar baseline. We investigated the potential of using absolute astrometry of the host star to obtain a dynamical mass constraint for the planet. The astrometric acceleration of 51~Eri derived from a comparison of the {\it Hipparcos} and {\it Gaia} catalogues was found to be inconsistent at the 2--3$σ$ level with the predicted reflex motion induced by the orbiting planet. Potential sources of this inconsistency include a combination of random and systematic errors between the two astrometric catalogs or the signature of an additional companion within the system interior to current detection limits. We also explored the potential of using {\it Gaia} astrometry alone for a dynamical mass measurement of the planet by simulating {\it Gaia} measurements of the motion of the photocenter of the system over the course of the extended eight-year mission. We find that such a measurement is only possible ($>98$\% probability) given the most optimistic predictions for the {\it Gaia} scan astrometric uncertainties for bright stars, and a high mass for the planet ($\gtrsim3.6$\,M$_{\rm Jup}$). △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 17 pages, 11 figures. Accepted for publication in the Astronomical Journal

Revised Astrometric Calibration of the Gemini Planet Imager

Authors: Robert J. De Rosa, Meiji M. Nguyen, Jeffrey Chilcote, Bruce Macintosh, Marshall D. Perrin, Quinn Konopacky, Jason J. Wang, Gaspard Duchêne, Eric L. Nielsen, Julien Rameau, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Tara Cotten, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Li-Wei Hung , et al. (27 additional authors not shown)

Abstract: We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in redu… ▽ More We present a revision to the astrometric calibration of the Gemini Planet Imager (GPI), an instrument designed to achieve the high contrast at small angular separations necessary to image substellar and planetary-mass companions around nearby, young stars. We identified several issues with the GPI Data Reduction Pipeline (DRP) that significantly affected the determination of angle of north in reduced GPI images. As well as introducing a small error in position angle measurements for targets observed at small zenith distances, this error led to a significant error in the previous astrometric calibration that has affected all subsequent astrometric measurements. We present a detailed description of these issues, and how they were corrected. We reduced GPI observations of calibration binaries taken periodically since the instrument was commissioned in 2014 using an updated version of the DRP. These measurements were compared to observations obtained with the NIRC2 instrument on Keck II, an instrument with an excellent astrometric calibration, allowing us to derive an updated plate scale and north offset angle for GPI. This revised astrometric calibration should be used to calibrate all measurements obtained with GPI for the purposes of precision astrometry. △ Less

Submitted 8 April, 2020; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 36 pages, 17 figures, accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems. Updated version includes revisions made during the referee process

A Search for Intermediate Separation Low Mass Binaries in the Orion Nebula Cluster

Authors: Matthew De Furio, Megan Reiter, Michael Meyer, Alexandra Greenbaum, Trent Dupuy, Adam Kraus

Abstract: We present the results of a binary population study in the Orion Nebula Cluster (ONC) using archival Hubble Space Telescope (HST) data obtained with the Advanced Camera for Surveys (ACS) in Johnson V filter (HST Proposal 10246, PI M. Robberto). Young clusters and associations hold clues to the origin and properties of multiple star systems. Binaries with separations $< 100 $ AU are useful as trace… ▽ More We present the results of a binary population study in the Orion Nebula Cluster (ONC) using archival Hubble Space Telescope (HST) data obtained with the Advanced Camera for Surveys (ACS) in Johnson V filter (HST Proposal 10246, PI M. Robberto). Young clusters and associations hold clues to the origin and properties of multiple star systems. Binaries with separations $< 100 $ AU are useful as tracers of the initial binary population since they are not as likely to be destroyed through dynamical interactions. Low mass, low stellar density star-forming regions such as Taurus-Auriga, reveal an excess of multiples compared to the Galactic Field. Studying the binary population of higher mass, higher stellar density star-forming regions like the ONC provides useful information concerning the origin of the Galactic Field star population. In this survey, we characterize the previously unexplored (and incomplete) separation parameter space of binaries in the ONC (15 - 160 AU) by fitting a double-PSF model built from empirical PSFs. We identified 14 candidate binaries (11 new detections) and find that 8$_{-2\%}^{+4\%}$ of our observed sample are in binary systems, complete over mass ratios and separations of 0.6 $< $ q $< $ 1.0 and 30 $< $ a $< $ 160 AU. This is consistent with the Galactic Field M-dwarf population over the same parameter ranges, 6.5\% $\pm$ 3\%. Therefore, high mass star forming regions like the ONC would not require further dynamical evolution for their binary population to resemble the Galactic Field, as some models have hypothesized for young clusters. △ Less

Submitted 4 October, 2019; originally announced October 2019.

Comments: 14 pages, 3 figures

Asymmetries in adaptive optics point spread functions

Authors: Alexander Madurowicz, Bruce Macintosh, Vanessa P. Bailey, Jeffrey Chilcote, Marshall Perrin, Lisa Poyneer, Laurent Pueyo, Jean-Baptiste Ruffio, Travis Barman, Joanna Bulger, Tara Cotten, Robert J. De Rosa, Rene Doyon, Gaspard Duchêne, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon, Li-Wei Hung, Patrick Ingraham, Paul Kalas , et al. (23 additional authors not shown)

Abstract: An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with… ▽ More An explanation for the origin of asymmetry along the preferential axis of the PSF of an AO system is developed. When phase errors from high altitude turbulence scintillate due to Fresnel propagation, wavefront amplitude errors may be spatially offset from residual phase errors. These correlated errors appear as asymmetry in the image plane under the Fraunhofer condition. In an analytic model with an open-loop AO system, the strength of the asymmetry is calculated for a single mode of phase aberration, which generalizes to two dimensions under a Fourier decomposition of the complex illumination. Other parameters included are the spatial offset of the AO correction, which is the wind velocity in the frozen flow regime multiplied by the effective AO time delay, and propagation distance or altitude of the turbulent layer. In this model, the asymmetry is strongest when the wind is slow and nearest to the coronagraphic mask when the turbulent layer is far away, such as when the telescope is pointing low towards the horizon. A great emphasis is made about the fact that the brighter asymmetric lobe of the PSF points in the opposite direction as the wind, which is consistent analytically with the clarification that the image plane electric field distribution is actually the inverse Fourier transform of the aperture plane. Validation of this understanding is made with observations taken from the Gemini Planet Imager, as well as being reproducible in end-to-end AO simulations. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 28 Pages, 13 Figures, Accepted to JATIS

Kernel-phase Detection Limits : Hypothesis Testing and the Example of JWST NIRISS Full Pupil Images

Authors: Alban Ceau, David Mary, Alexandra Greenbaum, Frantz Martinache, Anand Sivaramakrishnan, Romain Laugier, Mamadou N'Diaye

Abstract: The James Webb Space Telescope will offer high-angular resolution observing capability in the near-infrared with masking interferometry on NIRISS, and coronagraphic imaging on NIRCam & MIRI. Full aperture kernel-phase based interferometry complements these observing modes, probing for companions at small separations while preserving the telescope throughput. Our goal is to derive both theoretica… ▽ More The James Webb Space Telescope will offer high-angular resolution observing capability in the near-infrared with masking interferometry on NIRISS, and coronagraphic imaging on NIRCam & MIRI. Full aperture kernel-phase based interferometry complements these observing modes, probing for companions at small separations while preserving the telescope throughput. Our goal is to derive both theoretical and operational contrast detection limits for the kernel-phase analysis of JWST NIRISS full-pupil observations by using tools from hypothesis testing theory, applied to observations of faint brown dwarfs with this instrument, but the tools and methods introduced here are applicable in a wide variety of contexts. We construct a statistically independent set of observables from aberration-robust kernel phases. Three detection tests based on these observable quantities are designed and analysed, all guaranteeing a constant false alarm rate for small phase aberrations. One of these tests, the Likelihood Ratio or Neyman-Pearson test, provides a theoretical performance bound for any detection test. The operational detection method considered here is shown to exhibit only marginal power loss with respect to the theoretical bound. In principle, for the test set to a false alarm probability of 1%, companion at contrasts reaching 10^3 at separations of 200 mas around objects of magnitude 14.1 are detectable. With JWST NIRISS, contrasts of up to 10^4 at separations of 200 mas could be ultimately achieved, barring significant wavefront drift. The proposed detection method is close to the ultimate bound and offers guarantees over the probability of making a false detection for binaries, as well as over the error bars for the estimated parameters of the binaries detectable by JWST NIRISS. This method is not only applicable to JWST NIRISS but to any imaging system with adequate sampling. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 11 pages, 11 figures

Journal ref: A&A 630, A120 (2019)

An Exo-Kuiper Belt and An Extended Halo around HD 191089 in Scattered Light

Authors: Bin Ren, Élodie Choquet, Marshall D. Perrin, Gaspard Duchêne, John H. Debes, Laurent Pueyo, Malena Rice, Christine Chen, Glenn Schneider, Thomas M. Esposito, Charles A. Poteet, Jason J. Wang, S. Mark Ammons, Megan Ansdell, Pauline Arriaga, Vanessa P. Bailey, Travis Barman, Juan Sebastián Bruzzone, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Robert J. De Rosa, Rene Doyon, Michael P. Fitzgerald, Katherine B. Follette , et al. (48 additional authors not shown)

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

Submitted 31 July, 2019; originally announced August 2019.

Comments: 29 pages, 18 figures, ApJ accepted

On the Convergence Rate of Variants of the Conjugate Gradient Algorithm in Finite Precision Arithmetic

Authors: Anne Greenbaum, Hexuan Liu, Tyler Chen

Abstract: We consider three mathematically equivalent variants of the conjugate gradient (CG) algorithm and how they perform in finite precision arithmetic. It was shown in [{\em Behavior of slightly perturbed Lanczos and conjugate-gradient recurrences}, Lin.~Alg.~Appl., 113 (1989), pp.~7-63] that under certain conditions the convergence of a slightly perturbed CG computation is like that of exact CG for a… ▽ More We consider three mathematically equivalent variants of the conjugate gradient (CG) algorithm and how they perform in finite precision arithmetic. It was shown in [{\em Behavior of slightly perturbed Lanczos and conjugate-gradient recurrences}, Lin.~Alg.~Appl., 113 (1989), pp.~7-63] that under certain conditions the convergence of a slightly perturbed CG computation is like that of exact CG for a matrix with many eigenvalues distributed throughout tiny intervals about the eigenvalues of the given matrix, the size of the intervals being determined by how closely these conditions are satisfied. We determine to what extent each of these variants satisfies the desired conditions, using a set of test problems and show that there is significant correlation between how well these conditions are satisfied and how well the finite precision computation converges before reaching its ultimately attainable accuracy. We show that for problems where the width of the intervals containing the eigenvalues of the associated exact CG matrix makes a significant difference in the behavior of exact CG, the different CG variants behave differently in finite precision arithmetic. For problems where the interval width makes little difference or where the convergence of exact CG is essentially governed by the upper bound based on the square root of the condition number of the matrix, the different CG variants converge similarly in finite precision arithmetic until the ultimate level of accuracy is achieved, although this ultimate level of accuracy may be different for the different variants. This points to the need for testing new CG variants on problems that are especially sensitive to rounding errors. △ Less

Submitted 16 July, 2021; v1 submitted 14 May, 2019; originally announced May 2019.

MSC Class: 65F10; 65Y05

Performance of the Gemini Planet Imager Non-Redundant Mask and spectroscopy of two close-separation binaries HR 2690 and HD 142527

Authors: Alexandra Z. Greenbaum, Anthony Cheetham, Anand Sivaramakrishnan, Fredrik T. Rantakyrö, Gaspard Duchêne, Peter Tuthill, Robert J. De Rosa, Rebecca Oppenheimer, Bruce Macintosh, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Joanna Bulger, Andrew Cardwell, Jeffrey Chilcote, Tara Cotten, Rene Doyon, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Pascale Hibon, Li-Wei Hung, Patrick Ingraham , et al. (29 additional authors not shown)

Abstract: The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectr… ▽ More The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$σ$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\simλ/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $σ\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: Accepted to AJ, 22 pages, 14 figures

The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics From 10-100 AU

Authors: Eric L. Nielsen, Robert J. De Rosa, Bruce Macintosh, Jason J. Wang, Jean-Baptiste Ruffio, Eugene Chiang, Mark S. Marley, Didier Saumon, Dmitry Savransky, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Celia Blain, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Ian Czekala, Rene Doyon, Gaspard Duchene, Thomas M. Esposito, Daniel Fabrycky, Michael P. Fitzgerald, Katherine B. Follette, Jonathan J. Fortney, Benjamin L. Gerard , et al. (40 additional authors not shown)

Abstract: We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation be… ▽ More We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^αa^β$, finding $α$ = -2.4 $\pm$ 0.8 and $β$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: 52 pages, 18 figures. AJ in press

First direct detection of an exoplanet by optical interferometry; Astrometry and K-band spectroscopy of HR8799 e

Authors: S. Lacour, M. Nowak, J. Wang, O. Pfuhl, F. Eisenhauer, R. Abuter, A. Amorim, N. Anugu, M. Benisty, J. P. Berger, H. Beust, N. Blind, M. Bonnefoy, H. Bonnet, P. Bourget, W. Brandner, A. Buron, C. Collin, B. Charnay, F. Chapron, Y. Clénet, V. Coudé du Foresto, P. T. de Zeeuw, C. Deen, R. Dembet , et al. (63 additional authors not shown)

Abstract: To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes o… ▽ More To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,μ$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: published in A&A

Substellar Multiplicity Throughout the Ages

Authors: Daniella Bardalez Gagliuffi, Kimberly Ward-Duong, Jacqueline Faherty, Alex Greenbaum, Federico Marocco, Adam Burgasser, Matthew Bate, Trent Dupuy, Christopher Gelino, Johannes Sahlmann, Frantz Martinache, Michael Meyer, Quinn Konopacky, Denise Stephens

Abstract: Substellar multiplicity is a key outcome of the formation process. The biggest challenge for the next decade will be to distinguish between the formation history, environmental conditions, and dynamical evolution leading to the least massive brown dwarfs and the most massive planets at the tail ends of their mass functions. In this white paper, we advocate for a comprehensive characterization of b… ▽ More Substellar multiplicity is a key outcome of the formation process. The biggest challenge for the next decade will be to distinguish between the formation history, environmental conditions, and dynamical evolution leading to the least massive brown dwarfs and the most massive planets at the tail ends of their mass functions. In this white paper, we advocate for a comprehensive characterization of both the statistical distributions of the population of ultracool dwarf multiple systems and the fundamental properties of their individual components as a function of age. A space-based precision astrometry mission in near-infrared wavelengths would provide the necessary measurements to identify and characterize age-calibrated populations of multiple systems. △ Less

Submitted 15 March, 2019; originally announced March 2019.