# Revised astrometric calibration of the Gemini Planet Imager

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, et al.

► To cite this version:

Robert J. de Rosa, Meiji M. Nguyen, Jeffrey Chilcote, Bruce Macintosh, Marshall D. Perrin, et al. Revised astrometric calibration of the Gemini Planet Imager. Journal of Astronomical Telescopes Instruments and Systems, 2020, 6, 10.1117/1.JATIS.6.1.015006. insu-03705228

HAL Id: insu-03705228

https://insu.hal.science/insu-03705228v1

Submitted on 28 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Revised astrometric calibration of the
Gemini Planet Imager

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

Patrick Ingraham

Paul Kalas

James E. Larkin

Jérôme Maire

Franck Marchis

Mark S. Marley

Christian Marois

Stanimir Metchev

Maxwell A. Millar-Blanchaer

Rebecca Oppenheimer

David Palmer

Jennifer Patience

Lisa Poyneer

Laurent Pueyo

Abhijith Rajan

Fredrik T. Rantakyrö

Jean-Baptiste Ruffio

Dmitry Savransky

Adam C. Schneider

Anand Sivaramakrishnan

Inseok Song

Remi Souummer

Sandrine Thomas

J. Kent Wallace

Kimberly Ward-Duong

Sloane Wiktorowicz

Schuyler Wolff

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, Patrick Ingraham, Paul Kalas, James E. Larkin, Jérôme Maire, Franck Marchis, Mark S. Marley, Christian Marois, Stanimir Metchev, Maxwell A. Millar-Blanchaer, Rebecca Oppenheimer, David Palmer, Jennifer Patience, Lisa Poyneer, Laurent Pueyo, Abhijith Rajan, Fredrik T. Rantakyrö, Jean-Baptiste Ruffio, Dmitry Savransky, Adam C. Schneider, Anand Sivaramakrishnan, Inseok Song, Remi Souummer, Sandrine Thomas, J. Kent Wallace, Kimberly Ward-Duong, Sloane Wiktorowicz, Schuyler Wolff, "Revised astrometric calibration of the Gemini Planet Imager," *J. Astron. Telesc. Instrum. Syst.* **6**(1), 015006 (2020), doi: 10.1117/1.JATIS.6.1.015006
Revised astrometric calibration
of the Gemini Planet Imager

Robert J. De Rosa,<sup>a,b,*</sup> Meiji M. Nguyen,<sup>c</sup> Jeffrey Chilcote,<sup>d</sup>

Bruce Macintosh,<sup>b</sup> Marshall D. Perrin,<sup>e</sup> Quinn Konopacky,<sup>f</sup>

Jason J. Wang,<sup>g</sup> Gaspard Duchêne,<sup>c,h</sup> Eric L. Nielsen,<sup>b</sup> Julien Rameau,<sup>h,i</sup>

S. Mark Ammons,<sup>j</sup> Vanessa P. Bailey,<sup>k</sup> Travis Barman,<sup>l</sup> Joanna Bulger,<sup>m,n</sup>

Tara Cotten,<sup>o</sup> Rene Doyon,<sup>i</sup> Thomas M. Esposito,<sup>c</sup> Michael P. Fitzgerald,<sup>p</sup>

Katherine B. Follette,<sup>q</sup> Benjamin L. Gerard,<sup>r,s</sup> Stephen J. Goodsell,<sup>t</sup>

James R. Graham,<sup>c</sup> Alexandra Z. Greenbaum,<sup>u</sup> Pascale Hibon,<sup>a</sup>

Li-Wei Hung,<sup>v</sup> Patrick Ingraham,<sup>w</sup> Paul Kalas,<sup>c,x</sup> James E. Larkin,<sup>p</sup>

Jérôme Maire,<sup>f</sup> Franck Marchis,<sup>x</sup> Mark S. Marley,<sup>y</sup> Christian Marois,<sup>s,r</sup>

Stanimir Metchev,<sup>z,aa</sup> Maxwell A. Millar-Blanchaer,<sup>k</sup>

Rebecca Oppenheimer,<sup>ab</sup> David Palmer,<sup>j</sup> Jennifer Patience,<sup>ac</sup> Lisa Poyneer,<sup>j</sup>

Laurent Pueyo,<sup>e</sup> Abhijith Rajan,<sup>e</sup> Fredrik T. Rantakyrö,<sup>ad</sup>

Jean-Baptiste Ruffio,<sup>b</sup> Dmitry Savransky,<sup>ae</sup> Adam C. Schneider,<sup>ac</sup>

Anand Sivaramakrishnan,<sup>e</sup> Inseok Song,<sup>o</sup> Remi Souummer,<sup>e</sup>

Sandrine Thomas,<sup>w</sup> J. Kent Wallace,<sup>k</sup> Kimberly Ward-Duong,<sup>q</sup>

Sloane Wiktorowicz,<sup>af</sup> and Schuyler Wolff<sup>ag</sup>

<sup>a</sup>European Southern Observatory, Alonso de Cordova, Vitacura, Santiago, Chile

<sup>b</sup>Stanford University, Kavli Institute for Particle Astrophysics and Cosmology, Stanford, California, United States

<sup>c</sup>University of California, Department of Astronomy, Berkeley, California, United States

<sup>d</sup>University of Notre Dame, Department of Physics, Notre Dame, Indiana, United States

<sup>e</sup>Space Telescope Science Institute, Baltimore, Maryland, United States

<sup>f</sup>University of California San Diego, Center for Astrophysics and Space Science, La Jolla, California, United States

<sup>g</sup>California Institute of Technology, Department of Astronomy, Pasadena, California, United States

<sup>h</sup>Université Grenoble Alpes/Centre National de la Recherche Scientifique, Institut de planétologie et d'astrophysique de Grenoble, Grenoble, France

<sup>i</sup>Université de Montréal, Institut de Recherche sur les Exoplanètes, Département de Physique, Montréal, Québec, Canada

<sup>j</sup>Lawrence Livermore National Laboratory, Livermore, California, United States

<sup>k</sup>Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, United States

<sup>l</sup>University of Arizona, Lunar and Planetary Laboratory, Tucson, Arizona, United States

<sup>m</sup>University of Hawaii, Institute for Astronomy, Honolulu, Hawaii, United States

<sup>n</sup>National Astronomical Observatory of Japan, Subaru Telescope, Hilo, Hawaii, United States

<sup>o</sup>University of Georgia, Department of Physics and Astronomy, Athens, Georgia, United States

<sup>p</sup>University of California, Department of Physics and Astronomy, Los Angeles, California, United States

<sup>q</sup>Amherst College, Physics and Astronomy Department, Amherst, Massachusetts, United States

<sup>r</sup>University of Victoria, Victoria, British Columbia, Canada

<sup>s</sup>National Research Council of Canada Herzberg, Victoria, British Columbia, Canada

<sup>t</sup>Gemini Observatory, Hilo, Hawaii, United States

<sup>u</sup>University of Michigan, Department of Astronomy, Ann Arbor, Michigan, United States

<sup>v</sup>National Park Service, Natural Sounds and Night Skies Division, Fort Collins, Colorado, United States

<sup>w</sup>Large Synoptic Survey Telescope, Tucson, Arizona, United States

<sup>x</sup>SETI Institute, Carl Sagan Center, Mountain View, California, United States

<sup>y</sup>NASA Ames Research Center, Mountain View, California, United States
<sup>2</sup>University of Western Ontario, Centre for Planetary Science and Exploration,  
Department of Physics and Astronomy, London, Ontario, Canada

aaStony Brook University, Department of Physics and Astronomy, Stony Brook, New York, United States

abAmerican Museum of Natural History, Department of Astrophysics, New York, United States

acArizona State University, School of Earth and Space Exploration, Tempe, Arizona, United States

adGemini Observatory, Casilla, La Serena, Chile

aeCornell University, Sibley School of Mechanical and Aerospace Engineering, Ithaca, New York, United States

afThe Aerospace Corporation, El Segundo, California, United States

agLeiden University, Leiden Observatory, Leiden, The Netherlands

**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 the 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. © The Authors. *Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.* [DOI: 10.1117/1.JATIS.6.1.015006]

**Keywords:** high-contrast imaging; astrometric calibration; Gemini Planet Imager; data processing.

Paper 19112 received Oct. 18, 2019; accepted for publication Jan. 21, 2020; published online
Feb. 13, 2020.

**1 Introduction**

The Gemini Planet Imager<sup>1,2</sup> (GPI) is an instrument, currently at the Gemini South telescope, Chile, that was designed to achieve high contrast at small angular separations to resolve planetary-mass companions around nearby, young stars. Many high-contrast imaging observations also require highly precise and accurate astrometry. One of the objectives of the large Gemini Planet Imager Exoplanet Survey<sup>3</sup> (GPIES) was to characterize via relative astrometry the orbits of the brown dwarfs and exoplanets imaged as a part of the campaign.<sup>4</sup> These measurements have been used to investigate the dynamical stability of the multiplanet HR 8799 system,<sup>5</sup> the interactions between substellar companions and circumstellar debris disks,<sup>6,7</sup> and to directly measure the mass of β Pictoris b.<sup>8</sup> Improved astrometric accuracy and precision can reveal systematic discrepancies between instruments that need to be considered when performing orbital fits using astrometric records from multiple instruments. Accurate, precise astrometry can also help with common proper motion confirmation or rejection of detected candidate companions.

Previous work has demonstrated that the location of a faint substellar companion relative
to the host star can be measured within a reduced and postprocessed GPI image to a precision

<footnote>*Address all correspondence to Robert J. De Rosa, E-mail: rderosa@eso.org</footnote>
of ~700'th of a pixel.⁹ Since GPI's science camera is an integral field spectrograph (IFS)/
polarimeter, “pixel” in this context means the spatial pixel sampling set by the IFS lenslet array
rather than of the subsequent Hawaii-2RG detector. Converting these precise measurements of
the relative position of the companion from pixels into an on-sky separation and position angle
(PA) requires a precise and accurate astrometric calibration of the instrument. The plate scale of
the instrument is required to convert from pixels in the reconstructed data cubes into arcseconds
and the angle of north on an image that has been derotated to put north up based on the astro-
metric information within the header. The previous astrometric calibration (a plate scale of
$14.166 \pm 0.007$ mas px⁻¹ and a north offset angle of -0.10 ± 0.13 deg) was based on obser-
vations of calibration binaries and multiple systems obtained during the first two years of oper-
ations of the instrument.⁴,¹⁰

In the course of several investigations using GPI that relied upon astrometric measurements,
over time it became apparent that there were potentially remaining systematic biases after
that calibration, particularly in regard to the north angle correction. This motivated a careful,
thorough calibration effort into GPI astrometry, an effort that eventually grew to include cross
checks of the GPI data processing pipeline, the performance of several Gemini observatory
systems, and a complete reanalysis of all astrometric calibration targets observed with GPI.

This paper presents the findings of those efforts and the resulting improved knowledge of GPI's astrometric calibration. After introducing some background information regarding GPI and the Gemini architecture (Sec. 2), we describe two issues that we identified and fixed in the data reduction pipeline (DRP) (Sec. 3), a retroactive calibration of clock biases affecting some GPI observations (Sec. 4), and a model to calibrate for small apparent PA changes in some observations, at small zenith distances (Sec. 5). With those issues corrected, we revisit the astrometric calibration of GPI based on observations of several calibration binaries and multiple systems (Secs. 6 and 7). Compared to the prior calibration values, we find no significant difference in the plate scale. However, we find a different value for the true north correction by +0.36 deg, along with tentative low-significance evidence for small gradual drifts in that correction over time. Finally, we discuss the effect of the revised astrometric calibration on the astrometric measurements of several substellar companions (Sec. 8).

## 2 GPI and Gemini Systems Architecture Context

### 2.1 GPI Optical Assemblies

The GPI¹,² combines three major optical assemblies (Fig. 1). The adaptive optics (AO) system is mounted on a single thick custom optical bench. The Cassegrain focus of the telescope is located within the AO assembly. On that bench, the beam encounters a linear thin-plate atmospheric dispersion corrector, steerable pupil-alignment fold mirror, an off-axis parabolic (OAP) relay to the first deformable mirror, and an OAP relay to the second deformable mirror. After that, the beam is refocused to f/64. The last optic on the AO bench is a wheel containing microdot-patterned coronagraphic apodizer masks.¹¹,¹² These apodizer masks also include a square grid pattern that induces a regular pattern of diffracted copies of the stellar point spread function (PSF).¹³,¹⁴

The second optical assembly is an infrared wavefront sensor known as the calibration (CAL) system.¹⁵ It contains the focal plane mask component of the coronagraph (a flat mirror with a central hole) and collimating and steering optics.

The third assembly is the IFS.¹⁶,¹⁷ The input collimated beam is refocused onto a grid of
lenslets that serve as the image focal plane of the system. After this, the spectrograph optics relay
and disperse the lenslet images, but since the beam has been segmented, these can no longer
introduce astrometric effects. The lenslet array samples the focal plane and produces a grid of
“spots” or micropupils, each of which is an image of the telescope pupil. The only aberrations
affecting the image quality of the field are from elements in front of the lenslet array.¹⁷

Each of these three assemblies is independently mounted by three bipods. The bipods are supported by a steel truss structure that attaches to a square front mounting plate. The mounting plate attaches to the Gemini Instrument Support Structure (ISS) with large fixed kinematic pins. The ISS is a rotating cube located just above the Cassegrain focus of the telescope.
<figure><img src="image_2.png" /><figcaption>**Fig. 1** (a) Computer aided design (CAD) rendering of the GPI assembly showing the AO, CAL, and IFS optical benches and the supporting truss structure and mounting plate. For scale, the mounting plate is 1.2 m on a side. (Note that this shows an earlier version of the truss while the as-built structure is slightly different.) (b) Schematic showing the light path through the three optical assemblies.</figcaption></figure>

In typical Gemini operations, the ISS rotator operates to keep the sky PA fixed on the science focal plane. High-contrast imaging typically instead tries to fix the telescope pupil on the science instrument to allow angular differential imaging<sup>18</sup> (ADI). In GPI's case, this is always done at a single orientation (corresponding to GPI's vertical axis parallel to the telescope vertical axis). In the simplest case, this would involve stopping all rotator motion. However, as discussed in Sec. 5, in some but not all observations, the observatory software instead tries to maintain the absolute (sky) vertical angle (VA) stationary on the science focal plane, which must be accounted for in astrometric observations.

## 2.2 Software Interface and IFS Operation

The software architecture for GPI and the Gemini South telescope is complex, as is typical for a major observatory. Simple operations often require interactions between several different computers. For example, taking an image with the IFS is a process that involves four separate computer systems; the main Gemini environment that runs the observatory's control software, GPI's top level computer (TLC) that is interfaced with each component of the instrument, the IFS "host" computer that acts as an interface between the UNIX-based TLC and the Windows-based detector software, and the IFS "brick" that interfaces directly with the Hawaii-2RG detector.<sup>17</sup> Three of these four computer systems are responsible for populating the flexible image transport system (FITS)<sup>19</sup> image header keywords appended to each image. The Gemini environment handles telescope-specific quantities such as the telescope mount position, the TLC handles keywords associated with other parts of the instrument such as the AO system, and the IFS brick records detector-specific quantities. Each of these computer systems also maintains its own clock, although only the clock of the Gemini and environment and the IFS brick are relevant for the purposes of this study. These clocks are used when appending various timestamps to FITS headers during the process of obtaining an image. In theory, these clocks should all be synchronized periodically with Gemini's Network Time Protocol (NTP) server.

The IFS camera is controlled by the IFS brick, a computer used to interface with the Teledyne JADE2 electronics and Hawaii-2RG detector. This computer is responsible for commanding the camera, calculating count rates for each pixel based on raw up-the-ramp (UTR) reads,<sup>20</sup> sending completed images back to the observatory computers and providing ancillary metadata including the start and end time of the exposure (UTSTART and UTEND) that are stored in the FITS header. The detector is operated almost exclusively in UTR mode; correlated double sampling (CDS) mode<sup>21</sup> images have been taken in the laboratory, but this mode is not available for a standard observing sequence. The IFS runs at a fixed pixel clocking rate of 1.45479 s for a full read or
reset of the detector. The IFS software allows for multiple exposures to be coadded together prior to writing an FITS file. This mode has lower operational overheads and greater operational efficiency compared to individual exposures, and therefore, is frequently used for short exposures (from 1.5 to 10 s per coadd) but not generally used for long exposures (60 s per coadd) due to field rotation.

## 3 Improvements in the GPI Data Reduction Pipeline

The GPI DRP<sup>22,23</sup> is an open-source pipeline that performs basic reduction steps on data obtained with GPI's IFS to remove a variety of instrumental systematics and produce science-ready spectrophotometrically and astrometrically calibrated data cubes. The DRP corrects for detector dark current, identifies and corrects bad pixels and cosmic ray events, extracts the microspectra in the two-dimensional (2-D) image to construct a three-dimensional (3-D) (x, y, λ) data cube (or x, y, Stokes in polarimetry mode), and corrects for the small geometric distortion measured in the laboratory during the integration of the instrument.<sup>4</sup>

Critically, the DRP calculates the average parallactic angle between the start and end of an exposure, an angle that is used to rotate the reduced data cubes so that the vector toward celestial north is almost aligned with the columns of the image. We have identified and corrected in the latest data pipeline version two issues with the calculation of average parallactic angle that affect a subset of GPI measurements. These issues are most pronounced for observations taken at a very small zenith distance, where the parallactic angle is changing very rapidly. An example dataset showing the combined effect of these issues, and those described in Secs. 4 and 5, on observations of a calibration binary is shown in Fig. 2.

<figure><img src="image_3.png" /></figure>

**Fig. 2** Two GPI images of the binary star HD 6307 demonstrating the error in the calculation of AVPARANG in the previous version of the pipeline. (a), (c) The old reduction and (b), (d) new reduction (right column) for two images (a), (b) ~3 min and (c), (d) <1 min before the target transited the meridian. Each image has been rotated such that north is up based on the value of AVPARANG in the header of the reduced image (white compass). We use the prime symbol to denote the fact that the old reduction does not correctly rotate north up. The original detector coordinate axes are also shown (yellow compass). Note the flip of the x axis due to the odd number of reflections within the instrument. A significant change in the sky PA of the companion is seen between the two images in (a), (c), due to a combination of the errors described in Sec. 3. The PA of the companion is stable after the revisions to the pipeline.
## 3.1 Calculation of Average Parallactic Angle from Precise Exposure Start and End Times

Calculating the time-averaged parallactic angle during the course of an exposure requires accurate and precise knowledge of the exact start and end times of that exposure. We found that the GPI DRP was not originally using a sufficiently precise value for the start time in the case of an exposure with more than one coadd. Doing this correctly requires an understanding of the low-level details of the UTR readout of the Hawaii-2RG detector and the surrounding GPI and Gemini software.

The header of a raw GPI FITS file contains four timestamps saved at various times during the acquisition of an image with the IFS: UT, MJD-OBS, UTSTART, and UTEND. The keywords UT and MJD-OBS contain the time at the moment the header keyword values were queried by the Gemini master process prior to the start of the exposure. UT is reported in the coordinated universal time (UTC) scale, whereas MJD-OBS is reported in the terrestrial time scale, a scale linked to the International Atomic Time that is running ~65 s ahead of UTC. Because these keywords are written during exposure setup by a different computer system, neither is a highly precise metric for the exact exposure time start. The other keywords (UTSTART and UTEND) are generated by the IFS brick upon receipt of the command to execute an exposure and after the final read of the last coadd has completed. These two timestamps are reported in the UTC scale. Because they are written by the same computer that directly controls the readout, these are more accurate values for exposure timing. UTSTART is written when the IFS software receives the command to start an exposure, but since the Hawaii-2RG will be in continuous reset mode between exposures, it must wait some fraction of a read time to complete the current reset before the requested exposure can begin. Thus, the true exposure start time will be some unknown fraction of a read time after UTSTART. The final keyword UTEND is written with negligible delay immediately at the moment the last read of the last pixel is concluded. A schematic diagram of the read and resets of the Hawaii-2RG is shown for two example exposures in Fig. 3.

The pipeline was, therefore, written under the assumption that the UTEND keyword provides the most accurate way to determine the true start and end time of each exposure, which, in turn, is used to calculate the average parallactic angle during the exposure. The effective end time of

<figure><img src="image_4.png" /><figcaption>**Fig. 3** Reads (blue) and resets (red) of the Hawaii-2RG for two example exposures: (a) a single coadd exposure with eleven reads and (b) a five coadd exposure with three reads per coadd. The Hawaii-2RG is in continuous reset mode prior to the start of an exposure. The UTSTART keyword is generated when the exposure is commanded by the IFS software, which can be up to one and a half times the read-out time prior to the start of the exposure. UTEND is generated at the end of the final read. The EXPSTART and EXPEND values are calculated by the pipeline. The erroneous formula for computing EXSPTART for exposures with coadds is shown in red in (b).</figcaption></figure>
the exposure can be calculated as occurring half a read time prior to UTEND, i.e., the time at which half of the detector pixels have been read. The effective start time of the exposure, i.e., when half of the detector pixels have been read for the first time and can be calculated working backward from UTEND toward UTSTART. We do so based on the read time ($t_{read}$), number of reads per coadd ($n_{read}$, where $n_{read} - 1$ multiplied by $t_{read}$ yields the integration time per coadd), and number of coadds ($n_{coadd}$). The pipeline writes two additional keywords to the science extension of the reduced FITS file that stores the calculated effective start (EXPSTART) and end (EXPEND) times of the exposure calculated using UTEND, $t_{read}$, $n_{read}$, and $n_{coadd}$. EXPSTART and EXPEND are then used to calculate the average parallactic angle over the course of the exposure, which is written as keyword AVPARANG.

Inadvertently, versions 1.4 and prior of the GPI pipeline contained an error in this calculation by not correctly accounting for the number of coadds. The total exposure time including overheads was calculated as $t_{exp} = t_{read} \times (n_{read} - 3/2)$, where $n_{read}$ is the number of reads per coadd. Instead, the exposure time is more correctly calculated as

$$t_{exp} = t_{read} \times [n_{coadd} \times (n_{read} + 1) - 2], \quad (1)$$

where the additional terms account for the extra resets that occur between each coadd. The effect of this error was negligible for single-coadd exposures, the most common type of exposures taken with GPI; 89% of on-sky observations were taken with a single coadd. For images with multiple coadds, the effect can be very significant, with the error on the estimated time elapsed during the complete observation of

$$\Delta t = t_{read} \times (n_{coadd} \times n_{read} + n_{coadd} - n_{read} - 1/2). \quad (2)$$

To demonstrate how large this error can get for exposures with multiple coadds, an exposure with an integration time of 1.45 s with 10 coadds has a $\Delta t$ of 40 s, an error equivalent to 98% of the actual time spent exposing (see Fig. 4). A large $\Delta t$ can cause a significant and systematic error in the parallactic angle used to rotate the reduced data cubes north up as EXPSTART and

<figure><img src="image_5.png" /><figcaption>**Fig. 4** Error in the calculated duration of an exposure as a function of the number of reads (approximately equivalent to the integration time per coadd divided by 1.45 s) and the number of coadds. Dashed lines denote contours of $\Delta T = 1, 10, 30, 60$, and 120 s. All unique combinations of $n_{read}$ and $n_{coadd}$ for all on-sky GPI images within the GPIES database are plotted. Combinations with more than 100 images are shown as red circles (size scaled by the number), whereas combinations with less than 100 are shown as small gray circles. The vast majority of GPI exposures are taken with a single coadd, but for some frames with multiple coadds $\Delta T$ exceeded 120 s.</figcaption></figure>
EXPSTOP header keywords are converted into the hour angle at the start and end of the exposure from which the parallactic angle is calculated. This is most pronounced for targets observed at a small zenith distance where the parallactic angle is changing most rapidly. This error not only affects the astrometry of substellar companions, but also the measurement of binaries observed with other instruments that were used to calibrate GPI's true north offset angle.

After this inaccuracy was discovered, the GPI pipeline was updated to perform the correct calculation, as in version 1.5.

## 3.2 Average Parallactic Angle During Transits

A second issue affecting a small number of observations is related to time-averaging during exposures that span transit.

The pipeline computes the average parallactic angle between the start and end of an exposure via Romberg's method. For northern targets that transit during an exposure, the function contains a discontinuity at an hour angle ($H$) of $H = 0$ rad where the parallactic angle jumps from $-\pi$ to $+\pi$. This discontinuity can easily be avoided by performing the integration between $H = H_0$ and $H = 0$ rad, and between $H = 0$ rad and $H = H_1$, where $H_0$ and $H_1$ are the hour angles at the start and end of the exposure. The prior versions 1.4 of the pipeline and earlier contained an error in how this calculation was performed. As an example, the average parallactic angle $p_{avg}$ for an exposure with $|H_0| < H_1$ was calculated as

$$p_{avg} = \frac{1}{H_1 - H_0} \left[ \int_{H_0}^{0} |p(H, \phi, \delta)| dH + \int_{0}^{H_1} p(H, \phi, \delta) dH \right], \quad (3)$$

rather than

$$p_{avg} = \frac{1}{H_1 - H_0} \left[ \int_{H_0}^{0} [p(H, \phi, \delta) + 2\pi] dH + \int_{0}^{H_1} p(H, \phi, \delta) dH \right]. \quad (4)$$

This error only affects sequences where the target star transited the meridian between the start and end of an exposure. The magnitude of this error depends on exactly when transit occurred relative to the start and the end and the declination of the target. The net effect of the error on companion astrometry is small as it will only affect one of ~40 images taken in a typical GPI observing sequence.

This issue has also been corrected as of the latest version of the GPI pipeline.

## 4 Inaccuracies in Some FITS Header Time Information

The pipeline necessarily relies on the accuracy of the FITS header keywords in the data it is processing. However, it been proven that the FITS header keyword time information is not always as reliable as we would like. A review of FITS header timing information allowed us to uncover several periods in which misconfiguration or malfunction of time server software resulted in systematic errors in header keyword information. We were able to reconstruct the past history of such timing drifts sufficiently well as to be able to retroactively calibrate it out when reprocessing older data.

As a reminder, the UTSTART keyword is written by the IFS brick computer. The clock on the IFS brick is, at least in theory, configured to automatically synchronize once per week with Gemini's NTP server. This server provides a master time reference signal to maintain the accurate timings necessary for telescope pointing and control. In order to cause a noticeable error in the average parallactic angle, the IFS brick time stamps would have to be between a few and a few tens of seconds out of sync, depending on the declination of the star (Fig. 5). The regular synchronization of the clock on the IFS brick was intended to be sufficient to prevent it from drifting at such an amplitude relative to the time maintained by Gemini's NTP server.

However, it was eventually discovered that this time synchronization has not always operated as intended, resulting in significant clock offsets for some periods. The history of the offset
<figure><img src="image_6.png" /><figcaption>**Fig. 5** Error in the value of AVPARANG for a 60-s exposure induced by an IFS clock drift of between 1 and 60 s as a function of the hour angle at the start of the exposure for three target declinations. The error is most significant for targets that transit the meridian with a small zenith distance (bottom panel).</figcaption></figure>

between the IFS brick clock and UTC cannot be recovered directly from the various logs and
headers generated by the IFS. Instead, we can use the difference between the UT and UTSTART
header keywords as a proxy. The first timestamp is generated when the command to execute an
observation is issued by Gemini's Sequence Executor (SeqExec) and is assumed to be accurate;
a significant offset in the observatory's clock would quickly become apparent when attempting
to guide the telescope. The second timestamp is generated when the IFS brick receives the
command to start an exposure from the GPI TLC. The difference between these two timestamps,
UTSTART-UT, should be small and relatively stable, as there have not been any significant
changes to these software components since the instrument was commissioned in 2014, and
we show below that this time difference does prove to be stable for the majority of GPI data.

We, therefore, data mined all available GPI data to determine the time evolution of the offset between UT and UTSTART during the entire time GPI has been operational. We queried the GPIES Structured Query Language database,²⁴,²⁵ which contains the header information for all images obtained in the GPIES Campaign programs, selected guest observer (GO) programs whose principal investigators have contributed their data into this database and all public calibration programs. We augmented this with all GO programs that were publicly accessible in the Gemini Observatory Science Archive when this analysis was performed. We excluded engineering frames—images that are obtained via GPI’s interactive data language interface—as the UT keyword is populated via a different process for these types of frames. A total of 99,695 measurements of the UT to UTSTART offset spanning the previous six years were obtained, including 93,575 from the GPIES database and 6120 from other GO programs not included within the database.

The evolution of this offset between the installation of the instrument at Gemini South and
now is shown in Fig. 6. We identified several periods of time, two quite extended, where the IFS
clock was not correctly synchronized with the Gemini NTP. From the initial commissioning of
the instrument until the end of 2014, the offset varied significantly, from about 8 s slow to up to
<figure><img src="image_7.png" /><figcaption>**Fig. 6** Difference between UT and UTSTART in the header of coronagraphic (black), polarimetric (blue) images taken during the GPIES campaign, and for GO images (red). The green dashed line denotes the nominal offset between these keywords of 3.38 s (see Fig. 7). The median offset calculated as the median of all frames within a 12-h window (yellow solid line) was used to identify dates where the clock drift was significant (yellow shaded region).</figcaption></figure>

30 s fast. The causes of these variations are not fully known, but we point out that during this first year, GPI was still in commission and shared-risk science verification and software was still significantly in flux. In several instances, negative shifts in the offset are correlated with dates on which the IFS brick was used after having been restarted but prior to the periodic time synchronization occurred. The gradual negative drifts in offset observed at several points imply that the IFS clock was running too fast, gaining time at a rate of ~1 s per day over this period. Later, other small excursions in April 2016, August 2018, and August 2019 were also apparently caused by the IFS brick being used after an extended time powered off but prior to the scheduled weekly time synchronization. It would, of course, have been better had the time synchronization occurred automatically immediately after each reboot, but that was not the case.

A second long period with a significant offset, between June 2015 and March 2016, was caused by the IFS brick being synchronized to the wrong time server; it was tracking the Global Positioning System (GPS) time scale rather than UTC, and therefore, ran 18 s ahead
<figure><img src="image_8.png" /><figcaption>**Fig. 7** Histogram of the offset between UT and UTSTART for spectroscopic (black) and polarimetric (blue) observations taken between 2016.5 and 2018.5, when there were no clock biases. The width of the distribution is narrow relative to the 20 to 30 s offsets shown in the prior figure, supporting the notion that we can use drifts in UT-UTSTART to track the clock biases affecting UTSTART and UTEND.</figcaption></figure>

of UTC. An extended drift in the offset from April into May 2019 was caused by a failure in the NTP daemon running on a computer intermediate to the IFS brick and Gemini's NTP server. The drift was noticeably slower than in the 2013 to 2014 period, with the IFS brick gaining time at a rate of only one quarter of a second per day.

Improved systems administration can prevent such drifts in the future, but in order to properly calibrate the available data, we must model out the drifts that occurred in the past. The offset between UT and UTSTART remained relatively stable from mid-2016 through mid-2018 and was independent of the observing mode. We measured the median offset value between 2016.5 and 2018.5 as -3.38 s and defined this as the nominal UT to UTSTART offset (Fig. 7). We used a rolling median with a width of 12 h to calculate the value of the offset at a resolution of 1 h between late 2013 and 2019. A lookup table was created that the pipeline queries when reducing an IFS image so that it can apply a correction to UTSTART and UTEND if the observation was taken during a period identified as having a significant offset (Fig. 6).

# 5 Modeling Apparent Image Rotation at Gemini's Cassegrain Port

Recall from Sec. 2.1 that GPI always operates in ADI mode, with its pupil fixed or nearly fixed relative to the telescope pupil. GPI is attached to Gemini's ISS, which itself is mounted on the Cassegrain port of the telescope. A Cassegrain instrument rotator is used to maintain a fixed PA between the columns on an instrument's detector and either celestial north or the zenith. For an ideal altitude-azimuth telescope with an elevation axis perfectly aligned with local vertical and with an azimuth platform perpendicular to vertical, an instrument mounted on the Cassegrain port would observe the north angle changing with the parallactic angle as the telescope tracked a star through the meridian. The angle between the columns on the instrument detector and the direction of vertical would remain fixed (Fig. 8). Differences between true vertical and the vertical axis of the telescope cause this angle to vary slightly, an effect most pronounced for stars observed near the meridian with a small zenith distance (≤5 deg). When enabled, Gemini South's instrument rotator compensates for this motion, keeping the VA fixed on the detector (Fig. 9).

Due to difficulties maintaining the AO guide loops for targets with a very small zenith distance, it became common for some operators to keep the instrument rotator drive disabled while GPI was in operation, regardless of the target elevation. However, this practice was inconsistently applied. The drive was disabled and the rotator was kept at a nominal home position for 99 of the 317 nights on which GPI was used over the last 6 years. For data taken on these nights, a small correction needs to be applied to the parallactic angle in the header to compensate for this small motion of the VA as a star is tracked through the meridian.
<figure><img src="image_9.png" /><figcaption>**Fig. 8** Image (black) and sky (red) coordinate systems for observations taken (a), (d) before, (b), (e) during, and (c), (f) after meridian transit for (a)–(c) northerly and (d)–(f) southerly targets. The angle of the vertical vector (green) remains fixed relative to the image coordinate system for an ideal altitude-azimuth telescope, here at an angle of ~23.5 deg from the x axis within a reduced GPI data cube. Any offset between true vertical and the vertical axis of the telescope will cause the vertical vector within a reduced image to move slightly as the target crosses the meridian, the magnitude of which would be imperceptible in this diagram for a small offset as is the case for Gemini South, but significant relative to the precision of astrometric measurements made with GPI.</figcaption></figure>

Such a correction relies on precise knowledge of the telescope mount alignment. Sufficiently precise information on the Gemini South telescope mount is not publicly available. We, therefore, derived *post facto* knowledge of the Gemini South telescope mount based on the behavior of the Cassegrain rotator on nights when it was activated.

<figure><img src="image_10.png" /><figcaption>**Fig. 9** Angle of the instrument rotator as a function of hour angle for GPI observations where the rotator drive was enabled. The color of the symbol denotes the declination of the target. The instrument rotator angle has a different behavior for northern and southern targets due to the nonperpendicularity of the Gemini South telescope.</figcaption></figure>
We constructed a simple model to predict the correction to the parallactic angle caused by the nonperpendicular nature of the telescope.²⁶ For a perfect telescope, the parallactic angle of a source *p* is calculated as (Fig. 10)

$$ \tan p = \frac{-\cos \phi \sin A}{\sin \phi \cos E - \cos \phi \sin E \cos A'} \quad (5) $$

where *A* and *E* are the topocentric horizontal coordinates of the target, i.e., azimuth and elevation. If the telescope's azimuth platform is tilted at an angle of θ with an azimuth of Ω, the difference between the true *p* and apparent *p'* parallactic angles is²⁷

$$ p' - p = \Delta p = -\arctan\left(\frac{\cos \omega \sin \theta}{\cos E \cos \theta + \sin E \sin \theta \sin \omega}\right), \quad (6) $$

where ω = Ω - π/2 - A. A tilt in the elevation axis of θE within the plane connecting A ± π/2 causes an additional modification of

$$ \Delta p = -\arcsin(\sin \theta_E / \cos E). \quad (7) $$

These tilts will lead to a slight difference in the elevation and azimuth (E', A') of the telescope mount versus the topocentric elevation and azimuth (E, A) of the target. The telescope elevation and azimuth modified by the azimuth tilt are calculated as

$$ \sin E' = (\sin E \cos \theta - \cos E \sin \theta \sin \omega) \\ A' = \Omega - \arctan\left(\frac{\cos \omega \cos E}{-\cos \theta \sin \omega \cos E - \sin \theta \sin E}\right), \quad (8) $$

and due to an elevation tilt as

$$ \sin E' = \sin E / \cos \theta_E \quad A' = A - \arcsin(\tan E \tan \theta_E). \quad (9) $$

To construct a model of the tilt of the azimuth and elevation axes of the Gemini South telescope, we assumed that the instrument rotator was only compensating for the change in parallactic angle induced by these tilts. We collected measurements of the telescope elevation and azimuth and instrument rotator position on the 207 nights where GPI observations were taken with the rotator drive enabled. As the header stores the mechanical position of the telescope, we inverted the previous equations to compute the topocentric elevation and azimuth. Using these,

<figure><img src="image_11.png" /><figcaption>**Fig. 10** A tilt of the (a) azimuth and (b) elevation axes of the telescope can cause a significant change in the apparent parallactic angle (*p*) of a target (*S*). The magnitude of the tilt in azimuth (*θ*) and elevation (*θ*<sub>E</sub>) axes has been grossly exaggerated for the purpose of this diagram.</figcaption></figure>
<table><caption>**Table 1** Tilt model fit parameters.</caption><thead><tr><th>Start date (UT)</th><th>End date (UT)</th><th>θ (arc sec)</th><th>Ω (deg)</th><th>θ<sub>E</sub> (arc sec)</th><th>N<sub>frames</sub></th></tr></thead><tbody><tr><td>—</td><td>2014-06-30</td><td>27.3</td><td>38.9</td><td>16.8</td><td>3406</td></tr><tr><td>2014-07-01</td><td>2014-12-31</td><td>27.3</td><td>42.8</td><td>15.4</td><td>2787</td></tr><tr><td>2015-01-01</td><td>2015-09-16</td><td>29.5</td><td>43.8</td><td>20.1</td><td>5013</td></tr><tr><td>2015-09-17</td><td>2016-06-30</td><td>29.2</td><td>45.5</td><td>16.8</td><td>6323</td></tr><tr><td>2016-07-01</td><td>2016-12-31</td><td>28.1</td><td>40.5</td><td>16.5</td><td>1806</td></tr><tr><td>2017-01-01</td><td>2017-07-02</td><td>27.5</td><td>50.2</td><td>18.1</td><td>1641</td></tr><tr><td>2017-07-03</td><td>2017-12-31</td><td>26.4</td><td>37.6</td><td>19.7</td><td>2751</td></tr><tr><td>2018-01-01</td><td>2018-06-30</td><td>29.0</td><td>44.5</td><td>15.2</td><td>4682</td></tr><tr><td>2018-07-01</td><td>2019-01-18</td><td>29.5</td><td>45.4</td><td>16.1</td><td>3001</td></tr><tr><td>2019-01-19</td><td>—</td><td>31.1</td><td>51.8</td><td>19.1</td><td>1234</td></tr></tbody></table>

we predicted the change in parallactic angle, and thus the position that the instrument rotator would need to be at to compensate for nonperpendicularity for a given set of tilt parameters ($\theta$, $\Omega$, and $\theta_E$). We performed a least squares minimization to determine the set of tilt parameters that best reproduce the instrument rotator position for 10 roughly 6-month periods over the last 5 years. The break points were chosen arbitrarily to be at the start and midpoint of each year except for years in which a major earthquake occurred near Cerro Pachon (September 17, 2015, and January 19, 2019), and when a break point coincided with a period in which GPI was being used.

The tilt model parameters that best fit the measured instrument rotator positions are given in Table 1. A comparison between the model and data on the night of May 6, 2015, UT is shown in Fig. 11. The model is able to reproduce the commanded rotator positions with residuals smaller than the north calibration uncertainty (discussed below) in all but a handful of the images, specifically those taken at elevations ≥88 deg (Fig. 12).

<figure><img src="image_12.png" /><figcaption>**Fig. 11** Comparison between the sensed rotator angle (black) and that predicted by our simple telescope model (red) for (a) observations taken on May 06, 2015, and (b) the corresponding residuals. The model is able to reproduce the sensed angle well for targets at low elevation (small values of $\Delta p$), but performs worse at very high elevations.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. 12** Residuals between the sensed rotator angle and that predicted by the model for all observations in the GPIES database where the rotator drive was enabled, plotted (a) as a function of elevation, and (b) as a marginalized histogram on a logarithmic scale. The residuals are significant for observations taken at an elevation of $E > 88$ deg; 33 of the 32,644 images in the database have a residual >0.05 deg.</figcaption></figure>

We identified all GPI images to which we had access that were taken with the instrument rotator drive disabled. We used the tilt model parameters in Table 1 and the telescope elevation and azimuth within the header to calculate the correction to apply to the parallactic angle to compensate for the slight change in the angle of vertical on the detector. We created a lookup table with these corrections using the DATALAB header keyword to uniquely assign a correction to a specific GPI observation taken with the rotator drive disabled. Files with DATALAB values not in the lookup table do not have a correction applied. This lookup table contains all GPI observations taken with the drive disabled that were accessible at the time of this study, including GPIES campaign data, GO program data that are ingested into the GPIES database, and GO program data that were public at the time of the analysis.

# 6 North Angle Calibration

The corrections to the GPI DRP described in Secs. 3, 4, and 5 necessitated a revision of GPI's astrometric calibration, specifically the true north angle. The north angle offset is defined as the angle between IFS pixel columns and north in an image that has been rotated to put north up based on the average parallactic angle during the exposure. Here, we define the direction of the north angle offset as $\theta_{\text{true}} - \theta_{\text{observed}}$, a correction that would need to be added to a PA measured in images reduced with the GPI DRP (after correcting for the x axis flip) to recover the true PA of a companion.

We calibrate true north in GPI data based on observations of astrometric reference targets on sky. The small field of view (2.8 arc sec $\times$ 2.8 arc sec) and relatively bright limiting magnitude ($I < 10$) of GPI exclude many of the typical astrometric calibration fields used by other instruments (e.g., M15 and M92). Instead, we rely on periodic observations of a set of calibration binaries that have near-contemporaneous measurements with the well-calibrated NIRC2 camera on the Keck II telescope.²⁸,²⁹

## 6.1 Gemini South/GPI Observations

We have observed nine binary or multiple star systems since the start of routine operations in 2014. A summary of all these observations is given in Table 2. These observations were obtained with GPI's H band filter ($\lambda_{\text{eff}}$ = 1.64 µm) for all except two sequences taken with the K1 filter
<table><caption>Table 2 GPI observing log.</caption><thead><tr><th>Target</th><th>UT date</th><th>Mode</th><th>Filter</th><th>t<sub>int</sub> (s)</th><th>n<sub>coadd</sub></th><th>n<sub>exp</sub></th><th>ρ (px)</th><th>θ (deg)</th></tr></thead><tbody><tr><td>HD 1620</td><td>2015-08-30</td><td>C</td><td>H</td><td>1.45</td><td>10</td><td>23</td><td>41.342 ± 0.027</td><td>181.529 ± 0.040</td></tr><tr><td>HD 1620</td><td>2015-11-05</td><td>D</td><td>H</td><td>1.45</td><td>10</td><td>14</td><td>41.289 ± 0.055</td><td>181.202 ± 0.069</td></tr><tr><td>HD 1620</td><td>2018-07-21</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>10</td><td>41.024 ± 0.043</td><td>178.457 ± 0.058</td></tr><tr><td>HD 1620</td><td>2018-08-09</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>24</td><td>41.001 ± 0.020</td><td>178.434 ± 0.037</td></tr><tr><td>HD 1620</td><td>2018-09-21</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>24</td><td>40.981 ± 0.023</td><td>178.334 ± 0.039</td></tr><tr><td>HD 1620</td><td>2018-11-18</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>25</td><td>40.983 ± 0.027</td><td>178.272 ± 0.069</td></tr><tr><td>HD 1620</td><td>2018-12-20</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>17</td><td>40.982 ± 0.019</td><td>178.058 ± 0.068</td></tr><tr><td>HD 1620</td><td>2019-08-10</td><td>U</td><td>H</td><td>1.45</td><td>1</td><td>9</td><td>40.900 ± 0.048</td><td>(177.400 ± 0.086)</td></tr><tr><td>HD 1620</td><td>2019-08-10</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>7</td><td>40.878 ± 0.019</td><td>(177.377 ± 0.038)</td></tr><tr><td>HD 1620</td><td>2019-08-10</td><td>U</td><td>H</td><td>1.45</td><td>1, 10</td><td>16</td><td>40.882 ± 0.039</td><td>177.389 ± 0.070<sup>a</sup></td></tr><tr><td>HD 6307</td><td>2015-09-01</td><td>D</td><td>H</td><td>1.45</td><td>10</td><td>19</td><td>59.915 ± 0.029</td><td>237.081 ± 0.066</td></tr><tr><td>HD 6307</td><td>2019-08-10</td><td>U</td><td>H</td><td>4.36</td><td>1</td><td>19</td><td>60.365 ± 0.035</td><td>236.659 ± 0.032</td></tr><tr><td>HD 157516</td><td>2015-07-01</td><td>D</td><td>H</td><td>1.45</td><td>10</td><td>13</td><td>48.773 ± 0.025</td><td>142.511 ± 0.019</td></tr><tr><td>HD 157516</td><td>2015-07-29</td><td>D</td><td>H</td><td>1.45</td><td>10</td><td>7</td><td>48.759 ± 0.058</td><td>142.514 ± 0.052</td></tr><tr><td>HD 157516</td><td>2015-07-30</td><td>D</td><td>H</td><td>1.45</td><td>10</td><td>20</td><td>48.788 ± 0.041</td><td>142.457 ± 0.027</td></tr><tr><td>HD 158614</td><td>2019-08-11</td><td>ND</td><td>H</td><td>5.82</td><td>5</td><td>14</td><td>26.238 ± 0.017</td><td>127.714 ± 0.046</td></tr><tr><td>HIP 43947</td><td>2014-05-14</td><td>D</td><td>H</td><td>1.45</td><td>5</td><td>9</td><td>29.994 ± 0.010</td><td>(260.908 ± 0.020)</td></tr><tr><td>HIP 43947</td><td>2014-05-14</td><td>D</td><td>H</td><td>1.45</td><td>1</td><td>20</td><td>29.994 ± 0.056</td><td>(260.909 ± 0.059)</td></tr><tr><td>HIP 43947</td><td>2014-05-14</td><td>D</td><td>H</td><td>1.45</td><td>1, 5</td><td>29</td><td>29.994 ± 0.047</td><td>260.908 ± 0.051<sup>b</sup></td></tr><tr><td>HIP 43947</td><td>2015-01-24</td><td>D</td><td>H</td><td>1.45</td><td>5</td><td>12</td><td>29.991 ± 0.015</td><td>260.872 ± 0.013</td></tr><tr><td>HIP 43947</td><td>2015-04-02</td><td>D</td><td>H</td><td>1.45</td><td>5</td><td>12</td><td>29.982 ± 0.012</td><td>260.878 ± 0.016</td></tr><tr><td>HIP 43947</td><td>2015-04-23</td><td>D</td><td>H</td><td>1.45</td><td>5</td><td>12</td><td>29.978 ± 0.016</td><td>261.036 ± 0.027</td></tr><tr><td>HIP 44804</td><td>2014-03-23</td><td>D</td><td>K1</td><td>1.45</td><td>10</td><td>4</td><td>32.159 ± 0.009</td><td>306.035 ± 0.027</td></tr><tr><td>HIP 44804</td><td>2014-05-14</td><td>D</td><td>H</td><td>1.45</td><td>5</td><td>14</td><td>32.096 ± 0.069</td><td>305.866 ± 0.076</td></tr><tr><td>HIP 80628</td><td>2019-04-27</td><td>ND</td><td>H</td><td>8.73</td><td>3</td><td>12</td><td>69.367 ± 0.014</td><td>55.733 ± 0.016</td></tr><tr><td>HIP 80628</td><td>2019-08-10</td><td>ND</td><td>H</td><td>8.73</td><td>3</td><td>9</td><td>69.770 ± 0.012</td><td>56.362 ± 0.015</td></tr><tr><td>HR 7668</td><td>2016-09-21</td><td>U</td><td>H</td><td>1.45</td><td>10</td><td>5</td><td>37.337 ± 0.012</td><td>114.342 ± 0.017</td></tr><tr><td>HR 7668</td><td>2016-09-21</td><td>U</td><td>K1</td><td>1.45</td><td>10</td><td>15</td><td>37.332 ± 0.021</td><td>(114.308 ± 0.052)<sup>c</sup></td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2014-09-12</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>12</td><td>8.143 ± 0.059</td><td>222.881 ± 0.440</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2014-11-11</td><td>C</td><td>H</td><td>14.55</td><td>2</td><td>13</td><td>8.124 ± 0.021</td><td>223.718 ± 0.165</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2014-12-17</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.067 ± 0.028</td><td>224.067 ± 0.166</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2015-01-31</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.107 ± 0.016</td><td>223.816 ± 0.103</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2015-04-06</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>8</td><td>8.085 ± 0.031</td><td>223.923 ± 0.208</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2015-12-01</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.106 ± 0.035</td><td>224.920 ± 0.196</td></tr></tbody></table>
<table><caption>**Table 2 (Continued).**</caption><thead><tr><th>Target</th><th>UT date</th><th>Mode</th><th>Filter</th><th>t<sub>int</sub> (s)</th><th>n<sub>coadd</sub></th><th>n<sub>exp</sub></th><th>ρ (px)</th><th>θ (deg)</th></tr></thead><tbody><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2015-12-19</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.114 ± 0.024</td><td>224.894 ± 0.156</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-01-21</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.084 ± 0.020</td><td>225.076 ± 0.135</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-02-26</td><td>C</td><td>H</td><td>8.73</td><td>1</td><td>15</td><td>8.113 ± 0.037</td><td>224.837 ± 0.330</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-03-18</td><td>C</td><td>H</td><td>8.73</td><td>3</td><td>7</td><td>8.088 ± 0.017</td><td>225.059 ± 0.122</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-09-19</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.102 ± 0.020</td><td>225.845 ± 0.146</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-11-17</td><td>C</td><td>H</td><td>29.10</td><td>1</td><td>10</td><td>8.073 ± 0.026</td><td>226.106 ± 0.162</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-12-21</td><td>C</td><td>H</td><td>8.73</td><td>3</td><td>9</td><td>8.096 ± 0.023</td><td>225.936 ± 0.135</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2017-02-13</td><td>C</td><td>H</td><td>8.73</td><td>3</td><td>10</td><td>8.070 ± 0.018</td><td>226.289 ± 0.141</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2017-04-20</td><td>C</td><td>H</td><td>8.73</td><td>3</td><td>10</td><td>8.072 ± 0.027</td><td>226.430 ± 0.158</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2017-11-06</td><td>C</td><td>H</td><td>8.73</td><td>3</td><td>10</td><td>8.087 ± 0.026</td><td>226.947 ± 0.176</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2017-11-10</td><td>C</td><td>H</td><td>8.73</td><td>6</td><td>3</td><td>8.076 ± 0.019</td><td>226.860 ± 0.089</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-01-06</td><td>C</td><td>H</td><td>8.73</td><td>6</td><td>7</td><td>8.085 ± 0.010</td><td>227.008 ± 0.074</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-01-29</td><td>C</td><td>H</td><td>8.73</td><td>6</td><td>7</td><td>8.062 ± 0.028</td><td>227.410 ± 0.262</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-03-08</td><td>C</td><td>H</td><td>8.73</td><td>6</td><td>7</td><td>8.078 ± 0.029</td><td>227.127 ± 0.195</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-03-24</td><td>C</td><td>H</td><td>24.73</td><td>2</td><td>11</td><td>8.083 ± 0.032</td><td>227.303 ± 0.177</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-03-26</td><td>C</td><td>H</td><td>14.55</td><td>4</td><td>7</td><td>8.071 ± 0.011</td><td>227.481 ± 0.090</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-04-07</td><td>C</td><td>H</td><td>8.73</td><td>6</td><td>2</td><td>8.127 ± 0.060</td><td>227.125 ± 0.366</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-11-19</td><td>C</td><td>H</td><td>14.55</td><td>4</td><td>7</td><td>8.091 ± 0.016</td><td>228.055 ± 0.072</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2019-08-10</td><td>C</td><td>H</td><td>14.55</td><td>4</td><td>6</td><td>8.070 ± 0.032</td><td>228.762 ± 0.182</td></tr></tbody></table>

Note: C, coronagraphic; D, direct; ND, neutral density; and U, unblocked.

(λ<sub>eff</sub> = 2.06 µm); note that since the spectral filter in the GPI IFS is after the spatial pixellation at the lenslet array, change of filter cannot affect the astrometric calibration. The majority of the observations were obtained in GPI's "direct" mode, a configuration where the various coronagraphic components are removed from the optical path. Some were obtained in "unblocked" mode, which includes the Lyot mask and pupil plane apodizer in the optical path to reduce instrument throughput, preventing saturation from brighter stars. The addition of a neutral density filter in 2017 allowed us to observe calibrator binaries that were significantly brighter than the nominal H-band saturation limit of the IFS in either direct or unblocked mode. Observations of the θ<sup>1</sup> Ori B multiple system were taken in the coronagraphic mode, the typical mode for planet search observations, allowing for a high signal-to-noise ratio (SNR) detection of the fainter stellar components B2, B3, and B4 that all lie within an arcsecond of the primary star.

We do not expect the coronagraph optics to have a significant effect on astrometric measurements, except for those made for objects extremely close to the edge of the focal plane mask, which is not relevant here. The three coronagraph optics are in pupil and focal planes only, so cannot individually introduce distortions. By effectively weighting the beam profile across the pupil, they could, in principle, cause the beam to sample a different portion of any intermediate optics if those optics have polishing errors that could cause a slight field-dependent photocenter

<footnote><sup>a</sup>Calculated using all images obtained on 2019-08-10.</footnote>

<footnote><sup>b</sup>Calculated using all images obtained on 2014-05-14.</footnote>

<footnote><sup>c</sup>These data are not used for deriving the plate scale and north angle in Sec. 7.</footnote>
shift. However, this effect should be negligible. The intermediate optics (see Fig. 1) are small, located in a slow beam, and superpolished to ~1 nm rms wavefront error. Measured distortions are ~3 mas across the field of view⁴ and completely dominated by the geometric effects of the telephoto relay inside the spectrograph, with no evidence for a polishing-error component.

These observations were processed using version 1.5 (revision e0ea9f5) of the GPI DRP, incorporating the changes described in Secs. 3, 4, and 5. The data were all processed using the same DRP recipe with standard processing steps. The raw images were dark subtracted and corrected for bad pixels using both a static bad pixel map and outlier identification. The individual microspectra in each 2-D image were reassembled into a 3-D data cube (x, y, λ) using a wavelength solution derived from observations of a calibration argon arc lamp. An additional outlier identification and rejection step was performed on the individual slices of the data cubes. A distortion correction was then applied to each slice based on measurements of a pinhole mask taken during the commissioning of the instrument.⁴

## 6.2 Keck II/NIRC2 Observations

The same nine multiple systems have been observed with the NIRC2 instrument in conjunction with the facility AO system on the Keck II telescope. The isolated calibration binaries have between one and six NIRC2 epochs between 2014 and 2019. The Trapezium cluster that contains $\theta^1$ Ori B has been observed periodically with NIRC2 as an astrometric calibrator field by multiple different teams, with archival measurements extending as far back as December 2001. The observations were taken in a variety of instrument configurations and filters. A summary of these observations is given in Table 3. Datasets were taken in either PA mode, where north remains fixed at a given angle on the detector, or VA mode, where the VA remains fixed and north varies with the parallactic angle of the target.

We reduced these data using a typical near-infrared imaging DRP; correction for nonlinearity,³⁰ dark subtraction, flat fielding, and bad pixel identification and correction. Reduced images were corrected for geometric distortion using the appropriate distortion map.²⁸,²⁹ For observations taken using a subarray of the NIRC2 detector, we zero-padded the images prior to applying the distortion correction as the distortion correction script is hard-coded for 1024 × 1024 px images.³¹ The astrometric calibration of NIRC2 was derived from analyses of globular cluster observations and has been validated with measurements of the locations of SiO masers in the galactic center that were determined precisely using very long baseline radio interferometry measurements.²⁸,²⁹ We used a plate scale of 9.952 ± 0.002 mas px⁻¹ and a north angle offset of -0.252 ± 0.009 deg for data taken prior to April 13, 2015,²⁸ and 9.971 ± 0.005 mas px⁻¹ and a north angle offset of -0.262 ± 0.020 deg for data taken after.²⁹

## 6.3 Relative Astrometry

We used PSF fitting to measure the position of the companion relative to the primary. For the calibration binaries other than $\theta^1$ Ori B, we estimated the location of the primary star within each image (or wavelength slice) by fitting a 2-D Gaussian to a small 7 × 7 pixel stamp centered on an initial estimate of the primary star. The five parameters (x, y, $\sigma_x$, $\sigma_y$, and amplitude A) were allowed to vary except for the NIRC2 data obtained on 2019-04-25 (HIP 80628) and 2019-05-23 (HIP 44804), where $\sigma_x$ and $\sigma_y$ were fixed due to a strongly asymmetric PSF and the proximity of the companion. This process was repeated using the output of the first iteration as the initial guess for the second. We extracted a 15 × 15 px stamp centered on the fitted position of the primary to use as a template to fit the location of the secondary. We used the Nelder-Mead downhill simplex algorithm to determine the pixel offset and flux ratio between the primary and secondary stars by minimizing the squared residuals within a $2\lambda/D$ radius aperture surrounding the secondary. We estimated the uncertainty in the centroid of each fit as the full-width-at-half-maximum divided by the SNR measured as the peak pixel value divided by the standard deviation of pixel values within an annulus $15\lambda/D$ from the star. We corrected differential atmospheric refraction caused by the different zenith angle of the two stars using the model described in Ref. 32. We used the simplifying assumption that the observations were monochromatic at the central wavelength of the filter, negating any stellar color dependence on the effective
<table><caption>Table 3 NIRC2 observing log.</caption><thead><tr><th>Target</th><th>UT date</th><th>Filter</th><th>Rot. + mode</th><th>t<sub>int</sub> (s)</th><th>n<sub>coadd</sub></th><th>n<sub>exp</sub></th><th>ρ (mas)</th><th>θ (deg)</th></tr></thead><tbody><tr><td>HD 1620</td><td>2015-08-02</td><td>H2<sub>2-1</sub></td><td>PA</td><td>0.18</td><td>50</td><td>9</td><td>585.93 ± 0.41</td><td>181.740 ± 0.030</td></tr><tr><td>HD 6307</td><td>2015-08-02</td><td>K'</td><td>PA</td><td>0.181</td><td>100</td><td>9</td><td>848.36 ± 0.48</td><td>237.136 ± 0.031</td></tr><tr><td>HD 157516</td><td>2015-05-11</td><td>K<sub>cont</sub></td><td>VA</td><td>1.0</td><td>15</td><td>3</td><td>690.57 ± 0.50</td><td>142.678 ± 0.035</td></tr><tr><td>HD 158614</td><td>2014-05-13</td><td>Bry</td><td>PA</td><td>0.053</td><td>100</td><td>12</td><td>787.50 ± 0.35</td><td>147.406 ± 0.017</td></tr><tr><td>HD 158614</td><td>2019-08-17</td><td>Bry</td><td>PA</td><td>0.053</td><td>100</td><td>45</td><td>369.64 ± 0.22</td><td>128.014 ± 0.026</td></tr><tr><td>HD 158614</td><td>2019-08-26</td><td>Bry</td><td>PA</td><td>0.053</td><td>100</td><td>42</td><td>367.36 ± 0.22</td><td>127.799 ± 0.023</td></tr><tr><td>HIP 43947</td><td>2014-03-13</td><td>K'</td><td>VA</td><td>1.0</td><td>1</td><td>4</td><td>424.70 ± 0.46</td><td>260.948 ± 0.039</td></tr><tr><td>HIP 44804</td><td>2014-03-13</td><td>Bry</td><td>VA</td><td>0.5</td><td>5</td><td>4</td><td>455.19 ± 0.68</td><td>306.325 ± 0.021</td></tr><tr><td>HIP 44804</td><td>2019-05-23</td><td>H<sub>cont</sub></td><td>VA</td><td>2.0</td><td>10</td><td>16</td><td>444.73 ± 0.31</td><td>297.145 ± 0.049</td></tr><tr><td>HIP 80628</td><td>2014-03-13</td><td>Bry</td><td>VA</td><td>0.181</td><td>1</td><td>4</td><td>859.09 ± 0.39</td><td>43.640 ± 0.024</td></tr><tr><td>HIP 80628</td><td>2019-04-25</td><td>H<sub>cont</sub></td><td>PA</td><td>0.1</td><td>50</td><td>9</td><td>982.01 ± 0.54</td><td>56.187 ± 0.025</td></tr><tr><td>HIP 80628</td><td>2019-05-15</td><td>H<sub>cont</sub></td><td>PA</td><td>0.017</td><td>100</td><td>14</td><td>983.13 ± 0.59</td><td>56.327 ± 0.026</td></tr><tr><td>HIP 80628</td><td>2019-05-23</td><td>H<sub>cont</sub></td><td>VA</td><td>0.01</td><td>100</td><td>10</td><td>983.64 ± 0.63</td><td>56.242 ± 0.023</td></tr><tr><td>HIP 80628</td><td>2019-08-17</td><td>Bry</td><td>PA</td><td>0.0176</td><td>100</td><td>33</td><td>988.58 ± 0.65</td><td>56.855 ± 0.027</td></tr><tr><td>HIP 80628</td><td>2019-08-26</td><td>Bry</td><td>PA</td><td>0.0176</td><td>100</td><td>42</td><td>989.20 ± 0.56</td><td>56.881 ± 0.024</td></tr><tr><td>HIP 80628</td><td>2019-08-26</td><td>Bry</td><td>VA</td><td>0.0176</td><td>100</td><td>42</td><td>988.90 ± 0.63</td><td>56.863 ± 0.024</td></tr><tr><td>HR 7668</td><td>2016-07-22</td><td>Bry</td><td>PA</td><td>1.0</td><td>10</td><td>9</td><td>528.55 ± 0.41</td><td>114.725 ± 0.034</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2001-12-20</td><td>NB<sub>2.108</sub></td><td>PA</td><td>0.2</td><td>25</td><td>6</td><td>115.69 ± 0.40</td><td>209.32 ± 0.20</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2004-10-03</td><td>Bry</td><td>PA</td><td>0.2</td><td>100</td><td>2</td><td>116.97 ± 0.77</td><td>212.17 ± 0.38</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2005-02-16</td><td>NB<sub>2.108</sub></td><td>PA</td><td>0.2</td><td>50</td><td>3</td><td>116.34 ± 0.45</td><td>212.70 ± 0.22</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2005-02-25</td><td>Bry</td><td>PA</td><td>0.2</td><td>50</td><td>3</td><td>116.93 ± 0.30</td><td>212.94 ± 0.15</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2011-02-06</td><td>Bry</td><td>VA</td><td>0.726</td><td>1</td><td>6</td><td>114.97 ± 0.89</td><td>219.47 ± 0.44</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2011-02-06</td><td>Bry</td><td>PA</td><td>0.726</td><td>1</td><td>9</td><td>116.03 ± 0.71</td><td>219.35 ± 0.35</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2014-09-03</td><td>K'</td><td>VA</td><td>0.032</td><td>100</td><td>6</td><td>115.12 ± 0.14</td><td>223.90 ± 0.07</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2014-12-06</td><td>H</td><td>VA</td><td>0.053</td><td>100</td><td>15</td><td>115.41 ± 0.28</td><td>223.99 ± 0.14</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2015-10-27</td><td>Bry</td><td>PA</td><td>0.75</td><td>30</td><td>9</td><td>115.07 ± 0.23</td><td>224.93 ± 0.11</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-01-18</td><td>Bry</td><td>PA</td><td>0.75</td><td>30</td><td>10</td><td>115.52 ± 0.20</td><td>225.08 ± 0.10</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-02-04</td><td>Bry</td><td>VA</td><td>0.75</td><td>30</td><td>6</td><td>114.88 ± 0.16</td><td>225.14 ± 0.08</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-02-21</td><td>Bry</td><td>PA</td><td>0.181</td><td>300</td><td>9</td><td>115.17 ± 0.19</td><td>225.13 ± 0.09</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2016-08-20</td><td>K<sub>s</sub></td><td>PA</td><td>0.181</td><td>1</td><td>4</td><td>115.52 ± 0.44</td><td>226.23 ± 0.22</td></tr><tr><td>θ<sup>1</sup> Ori B2-B3</td><td>2018-02-13</td><td>Bry</td><td>PA</td><td>0.75</td><td>1</td><td>11</td><td>115.31 ± 0.18</td><td>227.59 ± 0.08</td></tr></tbody></table>
wavelength. This effect causes a reduction in the separation of a binary star along the elevation axis and was typically very small; at most 0.3 mas for the NIRC2 observations of HIP 80628 taken at an elevation of ~35 deg. PAs measured in datasets taken in VA mode were corrected by the parallactic angle at the middle of the exposure such that they were effectively measured relative from north.

The small angular separation between the two components of the $\theta^1$ Ori B2-B3 binary required us to use either $\theta^1$ Ori B1 for the NIRC2 observations or $\theta^1$ Ori B4 for the GPI observations as a reference PSF. We used this template PSF to simultaneously fit the location and fluxes of the two components of the B2-B3 binary following a similar procedure. We used a Fourier high-pass filter to subtract the seeing halo from B1 that was introducing a background signal for both B4 and the B2-B3 binary. The relative astrometries are listed in Table 2 for GPI and in Table 3 for NIRC2. We did not apply any correction for the differential atmospheric refraction for these observations given the extremely small difference in zenith angle between the two stars. We did not use the relative astrometry of B1-B2, B1-B3, or B1-B4 as B1 was obscured by GPI's focal plane mask, nor did we use B2-B4 or B3-B4 as the relative motion of these three stars cannot be described using a simple Keplerian model.

As a verification of the relative astrometry presented here, we performed an independent analysis of a subset of both the GPI and NIRC2 observations using the procedure described in Ref. 4. The GPI data were reduced with the same version of the DRP, whereas the NIRC2 data were reduced with a separate pipeline that performed the same functions as described in Sec. 6.2. Once the data were reduced, relative astrometry was performed using *StarFinder*.$^{33}$ For this subset of observations, we measured consistent separations and PAs to the values reported in Tables 2 and 3.

## 6.4 Accounting for Orbital Motion

Orbital motion of the calibration binaries between the NIRC2 and GPI epochs can introduce a significant bias in the north angle offset measurement. We fit Keplerian orbits to each of the calibration binaries using the NIRC2 astrometry presented in Table 3. These fits allowed us to simulate NIRC2 measurements on the same epoch as the GPI observations listed in Table 2, mitigating the bias induced by orbital motion. We use the parallel-tempered affine invariant Markov chain Monte Carlo (MCMC) package `emcee`$^{34}$ to sample the posterior distributions of the Campbell elements describing the visual orbit and of the system parallax. A complete description of the fitting procedure as applied to the 51 Eridani system can be found in Ref. 35. We used prior distributions for the system mass based on the blended spectral type and flux ratios of the components, and for the system parallax using measurements from either Hipparcos$^{36}$ or Gaia.$^{37}$ We used a parallax of $2.41 \pm 0.03$ mas for $\theta^1$ Ori B2-B3.$^{38}$ We also fitted the radial velocity measurements of both components of the HD 158614 binary$^{39}$ to help further constrain its orbital parameters. We purposely excluded astrometric measurements from other instruments and assumed that the NIRC2 astrometric calibration was stable before and after the realignment procedure in mid-2015.

The PA of the visual orbit and corresponding residuals are shown in Fig. 13 for the nine calibration binaries. We simulated NIRC2 measurements at the epoch of the GPI observations by drawing 10,000 orbits at random from MCMC chains and converting the orbital elements into separations and PAs at the desired epoch. We used the median of the resulting distribution of separations and PAs as the simulated measurement and the standard deviation as the uncertainty. These simulated measurements are reported in Table 4. The small semimajor axis of the HIP 43947 binary led to a significant uncertainty on the simulated NIRC2 observation despite the short 50-day baseline between the NIRC2 and GPI observations, precluding a measurement of the north offset angle with this binary. This was also the case for all but one epoch of both the HD 1620 and HD 6307 systems. Additional observations of these systems with NIRC2 to reduce the orbital uncertainties will be required for more precise predictions at these epochs. The remaining binaries (HD 157516, HD 158614, HIP 44804, HIP 80628, HR 7668, and $\theta^1$ Ori B2-B3) either had enough NIRC2 measurements to sufficiently constrain the orbit at the GPI epochs or were close enough in time that the orbital motion between the NIRC2 and GPI epochs was smaller than the measurement uncertainties.
<figure><img src="image_14.png" /><figcaption>**Fig. 13** (a) PA and (b) residuals of the orbits (blue lines) consistent with the NIRC2 astrometry in Table 3 (squares). The dates of GPI observations are highlighted; green dashed lines denote epochs that were used for the astrometric calibration, and red dotted lines denote epochs where the orbital motion is significant relative to the GPI measurement uncertainties. In a subset of the plots in (b), the date range has been restricted to focus on the dates of the GPI observations.</figcaption></figure>

# 7 Revised Astrometric Calibration

## 7.1 GPI Plate Scale

The plate scale for GPI was measured using the predicted separations in angular units from the orbit fit to the NIRC2 measurements and the pixel separations measured in the reduced GPI images (Table 4). We saw no evidence of a variation in the plate scale with time (Fig. 14) and adopted a single value of $14.161 \pm 0.021 \text{ mas px}^{-1}$. This measurement is consistent with the previous plate scale of $14.166 \pm 0.007 \text{ mas px}^{-1}$,⁴,¹⁰ but with a larger uncertainty. The pipeline changes described in Secs. 3, 4, and 5 have no impact on the separation of two stars within a reduced GPI image. The slight difference in the inferred plate scale can instead be ascribed to
<table><caption>Table 4 GPI plate scale and north offset angle.</caption><thead><tr><th>UT date</th><th>Target</th><th>ρ<sub>orbit</sub> (mas)</th><th>θ<sub>orbit</sub> (deg)</th><th>ρ<sub>orbit</sub>/ρ<sub>GPI</sub> (mas px<sup>-1</sup>)</th><th>θ<sub>orbit</sub> - θ<sub>GPI</sub> (deg)</th></tr></thead><tbody><tr><td>2014-03-23</td><td>HIP 44804</td><td>455.26 ± 0.63</td><td>306.274 ± 0.020</td><td>14.157 ± 0.020</td><td>0.239 ± 0.034</td></tr><tr><td>2014-05-14</td><td>HIP 43947</td><td>424.81 ± 12.42</td><td>260.872 ± 2.154</td><td>(14.163 ± 0.415)</td><td>(-0.036 ± 2.155)</td></tr><tr><td>2014-05-14</td><td>HIP 44804</td><td>455.06 ± 0.64</td><td>306.028 ± 0.020</td><td>14.178 ± 0.036</td><td>0.162 ± 0.079</td></tr><tr><td colspan="6">Weighted mean (2013-11-11 to 2014-09-08): 0.23 ± 0.11 deg</td></tr><tr><td>2014-09-12</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.88 ± 0.14</td><td>223.640 ± 0.052</td><td>14.123 ± 0.132</td><td>0.239 ± 0.414</td></tr><tr><td>2014-11-11</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.87 ± 0.14</td><td>223.826 ± 0.052</td><td>14.132 ± 0.044</td><td>0.189 ± 0.198</td></tr><tr><td>2014-12-17</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.86 ± 0.13</td><td>223.938 ± 0.052</td><td>14.178 ± 0.056</td><td>-0.143 ± 0.185</td></tr><tr><td>2015-01-24</td><td>HIP 43947</td><td>418.33 ± 65.15</td><td>260.540 ± 11.398</td><td>(13.949 ± 2.172)</td><td>(-0.332 ± 11.398)</td></tr><tr><td>2015-01-31</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.86 ± 0.13</td><td>224.078 ± 0.052</td><td>14.138 ± 0.032</td><td>0.241 ± 0.124</td></tr><tr><td>2015-04-02</td><td>HIP 43947</td><td>414.66 ± 80.02</td><td>260.435 ± 14.097</td><td>(13.830 ± 2.669)</td><td>(-0.443 ± 14.097)</td></tr><tr><td>2015-04-06</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.85 ± 0.12</td><td>224.283 ± 0.053</td><td>14.168 ± 0.055</td><td>0.308 ± 0.225</td></tr><tr><td>2015-04-23</td><td>HIP 43947</td><td>413.34 ± 84.75</td><td>260.404 ± 14.978</td><td>(13.788 ± 2.827)</td><td>(-0.632 ± 14.978)</td></tr><tr><td>2015-07-01</td><td>HD 157516</td><td>690.60 ± 0.59</td><td>142.678 ± 0.053</td><td>14.159 ± 0.014</td><td>0.167 ± 0.056</td></tr><tr><td>2015-07-29</td><td>HD 157516</td><td>690.60 ± 0.74</td><td>142.678 ± 0.071</td><td>14.164 ± 0.023</td><td>0.164 ± 0.088</td></tr><tr><td>2015-07-30</td><td>HD 157516</td><td>690.60 ± 0.74</td><td>142.678 ± 0.072</td><td>14.155 ± 0.019</td><td>0.221 ± 0.077</td></tr><tr><td>2015-08-30</td><td>HD 1620</td><td>585.89 ± 0.90</td><td>181.740 ± 0.105</td><td>14.172 ± 0.024</td><td>0.211 ± 0.112</td></tr><tr><td>2015-09-01</td><td>HD 6307</td><td>848.39 ± 0.78</td><td>237.138 ± 0.062</td><td>14.160 ± 0.015</td><td>0.057 ± 0.091</td></tr><tr><td colspan="6">Weighted mean (2014-09-08 to 2015-10-31): 0.17 ± 0.14 deg</td></tr><tr><td>2015-11-05</td><td>HD 1620</td><td>585.88 ± 2.82</td><td>181.727 ± 0.346</td><td>(14.190 ± 0.071)</td><td>(0.525 ± 0.353)</td></tr><tr><td>2015-12-01</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.84 ± 0.11</td><td>225.024 ± 0.055</td><td>14.135 ± 0.073</td><td>0.052 ± 0.200</td></tr><tr><td>2015-12-19</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.83 ± 0.11</td><td>225.080 ± 0.055</td><td>14.142 ± 0.049</td><td>0.204 ± 0.214</td></tr><tr><td>2016-01-21</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.83 ± 0.11</td><td>225.183 ± 0.055</td><td>14.179 ± 0.040</td><td>0.197 ± 0.177</td></tr><tr><td>2016-02-26</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.83 ± 0.11</td><td>225.295 ± 0.056</td><td>14.139 ± 0.082</td><td>0.319 ± 0.491</td></tr><tr><td>2016-03-18</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.83 ± 0.11</td><td>225.360 ± 0.056</td><td>14.171 ± 0.042</td><td>0.321 ± 0.155</td></tr><tr><td colspan="6">Weighted mean (2015-10-31 to 2016-09-05): 0.21 ± 0.23 deg</td></tr><tr><td>2016-09-19</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.82 ± 0.13</td><td>225.938 ± 0.059</td><td>14.135 ± 0.037</td><td>0.166 ± 0.144</td></tr><tr><td>2016-09-21</td><td>HR 7668</td><td>528.57 ± 0.52</td><td>114.727 ± 0.055</td><td>14.157 ± 0.015</td><td>0.385 ± 0.058</td></tr><tr><td>2016-11-17</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.82 ± 0.14</td><td>226.121 ± 0.060</td><td>14.126 ± 0.051</td><td>0.099 ± 0.190</td></tr><tr><td>2016-12-21</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.15</td><td>226.227 ± 0.061</td><td>14.164 ± 0.047</td><td>0.287 ± 0.163</td></tr><tr><td>2017-02-13</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.16</td><td>226.394 ± 0.062</td><td>14.217 ± 0.044</td><td>0.211 ± 0.184</td></tr><tr><td>2017-04-20</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.17</td><td>226.603 ± 0.064</td><td>14.194 ± 0.051</td><td>0.223 ± 0.176</td></tr><tr><td colspan="6">Weighted mean (2016-09-05 to 2017-10-13): 0.32 ± 0.15 deg</td></tr><tr><td>2017-11-06</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.22</td><td>227.224 ± 0.069</td><td>14.195 ± 0.052</td><td>0.254 ± 0.220</td></tr><tr><td>2017-11-10</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.22</td><td>227.237 ± 0.069</td><td>14.175 ± 0.045</td><td>0.303 ± 0.131</td></tr></tbody></table>
<table><caption>**Table 4 (Continued).**</caption><thead><tr><th>UT date</th><th>Target</th><th>ρ<sub>orbit</sub> (mas)</th><th>θ<sub>orbit</sub> (deg)</th><th>ρ<sub>orbit</sub>/ρ<sub>GPI</sub> (mas px<sup>-1</sup>)</th><th>θ<sub>orbit</sub> - θ<sub>GPI</sub> (deg)</th></tr></thead><tbody><tr><td>2018-01-06</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.23</td><td>227.414 ± 0.071</td><td>14.181 ± 0.035</td><td>0.377 ± 0.109</td></tr><tr><td>2018-01-29</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.24</td><td>227.485 ± 0.072</td><td>14.200 ± 0.057</td><td>0.179 ± 0.275</td></tr><tr><td>2018-03-08</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.25</td><td>227.603 ± 0.073</td><td>14.161 ± 0.063</td><td>0.369 ± 0.202</td></tr><tr><td>2018-03-24</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.26</td><td>227.653 ± 0.073</td><td>14.201 ± 0.071</td><td>0.234 ± 0.218</td></tr><tr><td>2018-03-26</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.26</td><td>227.660 ± 0.073</td><td>14.181 ± 0.039</td><td>0.162 ± 0.114</td></tr><tr><td>2018-04-07</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.26</td><td>227.700 ± 0.074</td><td>14.078 ± 0.116</td><td>0.179 ± 0.534</td></tr><tr><td>2018-07-21</td><td>HD 1620</td><td>584.69 ± 32.23</td><td>181.531 ± 3.961</td><td>(14.252 ± 0.786)</td><td>(3.074 ± 3.961)</td></tr><tr><td>2018-08-09</td><td>HD 1620</td><td>584.64 ± 32.79</td><td>181.528 ± 4.030</td><td>(14.259 ± 0.800)</td><td>(3.094 ± 4.030)</td></tr><tr><td colspan="6">Weighted mean (2017-10-13 to 2018-09-01): 0.28 ± 0.19 deg</td></tr><tr><td>2018-09-21</td><td>HD 1620</td><td>584.52 ± 34.08</td><td>181.520 ± 4.189</td><td>(14.263 ± 0.832)</td><td>(3.186 ± 4.189)</td></tr><tr><td>2018-11-18</td><td>HD 1620</td><td>584.41 ± 35.81</td><td>181.509 ± 4.402</td><td>(14.260 ± 0.874)</td><td>(3.237 ± 4.403)</td></tr><tr><td>2018-11-19</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.81 ± 0.33</td><td>228.402 ± 0.081</td><td>14.171 ± 0.050</td><td>0.267 ± 0.113</td></tr><tr><td>2018-12-20</td><td>HD 1620</td><td>584.28 ± 36.77</td><td>181.504 ± 4.520</td><td>(14.257 ± 0.897)</td><td>(3.446 ± 4.521)</td></tr><tr><td>2019-04-27</td><td>HIP 80628</td><td>982.22 ± 0.43</td><td>56.215 ± 0.020</td><td>14.160 ± 0.007</td><td>0.482 ± 0.026</td></tr><tr><td>2019-08-10</td><td>θ<sup>1</sup> Ori B2-B3</td><td>114.83 ± 0.42</td><td>229.224 ± 0.091</td><td>14.186 ± 0.079</td><td>0.299 ± 0.220</td></tr><tr><td>2019-08-10</td><td>HD 1620</td><td>583.63 ± 43.77</td><td>181.462 ± 5.385</td><td>(14.276 ± 1.071)</td><td>(4.073 ± 5.385)</td></tr><tr><td>2019-08-10</td><td>HD 6307</td><td>847.62 ± 30.81</td><td>237.156 ± 2.673</td><td>(14.042 ± 0.510)</td><td>(0.497 ± 2.673)</td></tr><tr><td>2019-08-10</td><td>HIP 80628</td><td>988.21 ± 0.38</td><td>56.801 ± 0.016</td><td>14.164 ± 0.006</td><td>0.439 ± 0.022</td></tr><tr><td>2019-08-11</td><td>HD 158614</td><td>371.35 ± 0.19</td><td>128.153 ± 0.017</td><td>14.153 ± 0.012</td><td>0.439 ± 0.049</td></tr><tr><td colspan="6">Weighted mean (2018-09-01 to 2019-08-27): 0.45 ± 0.11 deg</td></tr><tr><td colspan="6">Weighted mean (all): 0.36 ± 0.12 deg</td></tr><tr><td colspan="6">14.161 ± 0.021 mas px<sup>-1</sup></td></tr></tbody></table>

Note: measurements in parentheses are not included in weighted mean.

changes in the way the relative positions of the two components of each calibration binary were measured, the greater number of measurements, or simply to measurement uncertainties.

## 7.2 GPI North Offset Angle

The north offset angle for GPI was measured by taking the difference of the PA of the companion predicted from the NIRC2-only orbit fit ($\theta_{orbit}$) and the measured PA within the reduced GPI data cubes ($\theta_{GPI}$). This difference is reported in Table 4 for each calibration binary measurement. We calculated a weighted mean of $0.36 \pm 0.12$ deg for the full set of measurements, with the error calculated assuming that they were not independent. The aforementioned 0.12-deg uncertainty includes a 0.1-deg uncertainty that was added in quadrature to account for systematics and uncertainties in the relative astrometry. Measurements with large uncertainties in the predicted PA ($\theta_{orbit}$) were excluded. The measured offsets and the best fit model are plotted in Fig. 15(a). Although the model is consistent with the measurements given the sizes of the uncertainties on both the measurements and the model ($\chi^2 = 1.2, \nu = 36$), there does appear to be a slight
<figure><img src="image_15.png" /></figure>

**Fig. 14** Measurements of the plate scale of GPI derived from calibration binaries (red circles) and the $\theta^1$ Ori B2-B3 binary (black squares). The mean and standard deviation (blue solid line and shaded region) were calculated using a weighted mean and assuming that the measurements were not independent. The previous astrometric calibration is overplotted for reference (gray dashed line and shaded region).

trend of increasing north offset angle over the course of 6 years when comparing the calibration binary measurements in early-2014 and mid-2019.

One plausible cause of a rotation of the instrument with respect to the telescope is the annual shutdown of the telescope when both the instrument and ISS are removed to perform maintenance. We fit a variable north offset angle that remains static between the dates of telescope shutdowns. A series of weighted means were calculated using measurements between each shutdown, as listed in Table 4 and plotted in Fig. 15. This model reproduces the trend of increasing north offset angle during the previous 6 years and is an improved fit ($\chi^2 = 0.4, \nu = 31$) relative to the single-valued model. We opted to use this variable north offset angle model for the final astrometric calibration of the instrument.

<figure><img src="image_16.png" /></figure>

**Fig. 15** Measurements of the north offset angle of GPI derived from calibration binaries (red circles) and the $\theta^1$ Ori B2-B3 binary (black squares). We fit the north angle assuming it is either (a) a constant calibration for the entire date ranges or (b) that it varies between telescope shutdowns. The mean and standard deviation (blue solid line and shaded region) are calculated as in Fig. 14. The previous astrometric calibration is overplotted for reference (gray dashed line and shaded region).
<figure><img src="image_17.png" /><figcaption>**Fig. 16** (a) One wavelength slice of a reduced GPI data cube for a postalignment image taken using GPI's internal source on November 12, 2014. The four satellite spots generated by the grid on the pupil plane apodizer are clearly visible. The PA between the bottom left (S1) and top right (S2) satellite spot, measured from S1 to S2 counter-clockwise from vertical, (b) plotted as a function of date for each postalignment image taken since the instrument was commissioned.</figcaption></figure>

## 7.3 Instrument Stability

The cause of the change of the north offset angle over time is not known. In principle, a movement of the IFS or the CAL system on their bipod mounts could produce a clocking of the focal plane with respect to the telescope, although a movement of 5 mm would be required. We excluded rotations internal to the instrument by measuring the angle between two of the satellite spots within a postalignment image taken routinely before instrument operation. These satellite spots are generated by a periodic wire grid on the pupil plane apodizer,¹³,¹⁴ located on the AO bench (Fig. 1). A physical rotation of the IFS relative to the apodizer would manifest itself as a rotation of the satellite spots within the focal plane as recorded by the IFS. We measured the angle between the bottom left and top right satellite spots in 406 postalignment images taken between late-2014 and mid-2019 using the satellite spot finding algorithm that is a part of the GPI DRP. We find no significant trend in this angle over the past 5 years (Fig. 16), although a significant offset of ~0.1 deg is seen for a few months at the start of 2016 that coincides with mechanical difficulties with the wheel containing the pupil plane apodizers. Excluding this period, we find an angle between these two satellite spots of 335.96 ± 0.02 deg. The stability of this angle implies that the change in the north offset angle seen in Fig. 15 is caused by a mechanical rotation upstream of the pupil plane mechanism containing the apodizer. The GPI optics upstream of this are all rigidly mounted in a single plane onto a thick optical bench and are extremely unlikely to produce such a rotation. In principle, a rotation of the outer truss structure holding all three assemblies with respect to the mounting plate could rotate the focal plane, but again that would have to be on the order of 5 mm, essentially impossible. GPI has an extremely rigid truss structure supporting various subcomponents. Integrated finite element analysis/optical modeling shows that flexure motions of any component relative to the optical axis are <25 µm over the operating range of gravity vectors.⁴⁰ Although we did not explicitly model rotation, if any hypothetical rotation component involves displacements on the same scale, the angular rotation would be on the order of 0.01 deg. The pins that locate GPI onto the ISS face have much more precise tolerances than that as well (<0.23 mm).

## 8 Revised Astrometry for Substellar Companions

The changes to the pipeline described in Secs. 3, 4, and 5 and the revised astrometric calibration of the instrument described in Sec. 7 both necessitate a revision of previously published relative astrometry of substellar companions measured using GPI observations. Revisions for $\beta$ Pictoris b,⁸ 51 Eridani b,³⁵ and HD 206893 B⁴¹ are presented in other works. Here, we present corrections
<figure><img src="image_18.png" /><figcaption>**Fig. 17** The change in the average parallactic angle header keyword (AVPARANG) due to the revisions to the pipeline described in Secs. 3, 4, and 5 for four GPIES targets that have published astrometry of substellar companions. The change in parallactic angle varies for each epoch, and for HD 984 varies significantly within a single epoch. The error in the integration described in Sec. 3.2 is apparent in several epochs, most noticeably for the February 29, 2016, dataset on HD 95086 (middle row in left column).</figcaption></figure>

to the astrometry for the exoplanets in the HR 8799<sup>5</sup> and HD 95086<sup>7</sup> systems, and the brown dwarfs HR 2562 B<sup>42</sup> and HD 984 B<sup>43</sup> that correct for the changes to the pipeline and the revised astrometric calibration of the instrument. We reduced the same images used in the previous studies with the latest version of the GPI DRP. The revisions described in Secs. 3, 4, and 5 all affect the AVPARANG header keyword. The change in this value is plotted as a function of frame number for each observing sequence in Fig. 17. Δ AVPARANG is typically small and static, only changing by at most ~0.05 deg between the start and end of the J-band sequence on HD 984 taken on August 30, 2015. The effect of the parallactic angle integration error described in Sec. 3.2 is apparent in several epochs.

The median Δ AVPARANG was used in conjunction with the revised north offset angle described in Sec. 7 to revise the previously published astrometry. We assumed that a single offset
<table><caption>**Table 5** Revised companion astrometry.</caption><thead><tr><th>Object</th><th>Date (UT)</th><th>Band</th><th>ρ<sub>original</sub> (mas)</th><th>θ<sub>original</sub> (deg)</th><th>ρ<sub>revised</sub> (mas)</th><th>θ<sub>revised</sub> (deg)</th></tr></thead><tbody><tr><td>HR 8799 c</td><td>2013-11-17</td><td>K1</td><td>949.5 ± 0.5</td><td>325.18 ± 0.14</td><td>949.1 ± 1.4</td><td>325.51 ± 0.12</td></tr><tr><td>HR 8799 d</td><td>2013-11-17</td><td>K1</td><td>654.6 ± 0.9</td><td>214.15 ± 0.15</td><td>654.3 ± 1.3</td><td>214.48 ± 0.13</td></tr><tr><td>HR 8799 e</td><td>2013-11-17</td><td>K1</td><td>382.6 ± 2.1</td><td>265.13 ± 0.24</td><td>382.4 ± 2.2</td><td>265.46 ± 0.23</td></tr><tr><td>HR 8799 b</td><td>2014-09-12</td><td>H</td><td>1721.2 ± 1.4</td><td>65.46 ± 0.14</td><td>1720.5 ± 2.8</td><td>65.74 ± 0.15</td></tr><tr><td>HR 8799 c</td><td>2014-09-12</td><td>H</td><td>949.0 ± 1.1</td><td>326.53 ± 0.14</td><td>948.6 ± 1.7</td><td>326.81 ± 0.15</td></tr><tr><td>HR 8799 d</td><td>2014-09-12</td><td>H</td><td>662.5 ± 1.3</td><td>216.57 ± 0.17</td><td>662.2 ± 1.6</td><td>216.85 ± 0.18</td></tr><tr><td>HR 8799 c</td><td>2016-09-19</td><td>H</td><td>944.2 ± 1.0</td><td>330.01 ± 0.14</td><td>943.8 ± 1.7</td><td>330.43 ± 0.16</td></tr><tr><td>HR 8799 d</td><td>2016-09-19</td><td>H</td><td>674.5 ± 1.0</td><td>221.81 ± 0.15</td><td>674.2 ± 1.4</td><td>222.23 ± 0.17</td></tr><tr><td>HR 8799 e</td><td>2016-09-19</td><td>H</td><td>384.8 ± 1.7</td><td>281.68 ± 0.25</td><td>384.6 ± 1.8</td><td>282.10 ± 0.26</td></tr><tr><td>HD 95086 b</td><td>2013-12-10</td><td>K1</td><td>619.0 ± 5.0</td><td>150.90 ± 0.50</td><td>618.9 ± 4.9</td><td>151.10 ± 0.44</td></tr><tr><td>HD 95086 b</td><td>2013-12-11</td><td>H</td><td>618.0 ± 11.0</td><td>150.30 ± 1.10</td><td>617.8 ± 11.1</td><td>150.45 ± 1.11</td></tr><tr><td>HD 95086 b</td><td>2014-05-13</td><td>K1</td><td>618.0 ± 8.0</td><td>150.20 ± 0.70</td><td>617.7 ± 8.0</td><td>150.55 ± 0.71</td></tr><tr><td>HD 95086 b</td><td>2015-04-06</td><td>K1</td><td>622.0 ± 7.0</td><td>148.80 ± 0.60</td><td>621.9 ± 7.3</td><td>149.06 ± 0.64</td></tr><tr><td>HD 95086 b</td><td>2015-04-08</td><td>K1</td><td>622.0 ± 4.0</td><td>149.00 ± 0.40</td><td>621.7 ± 4.1</td><td>149.25 ± 0.39</td></tr><tr><td>HD 95086 b</td><td>2016-02-29</td><td>H</td><td>621.0 ± 5.0</td><td>147.80 ± 0.50</td><td>620.3 ± 4.8</td><td>148.09 ± 0.57</td></tr><tr><td>HD 95086 b</td><td>2016-03-06</td><td>H</td><td>620.0 ± 5.0</td><td>147.20 ± 0.50</td><td>619.8 ± 4.8</td><td>147.50 ± 0.57</td></tr><tr><td>HR 2562 B</td><td>2016-01-25</td><td>H</td><td>619.0 ± 3.0</td><td>297.56 ± 0.35</td><td>618.8 ± 3.0</td><td>297.76 ± 0.40</td></tr><tr><td>HR 2562 B</td><td>2016-01-28</td><td>K1</td><td>618.0 ± 5.0</td><td>297.40 ± 0.25</td><td>617.8 ± 5.1</td><td>297.50 ± 0.30</td></tr><tr><td>HR 2562 B</td><td>2016-01-28</td><td>K2</td><td>618.0 ± 4.0</td><td>297.76 ± 0.37</td><td>618.0 ± 4.1</td><td>297.88 ± 0.42</td></tr><tr><td>HR 2562 B</td><td>2016-02-25</td><td>K2</td><td>619.0 ± 2.0</td><td>297.50 ± 0.25</td><td>618.9 ± 2.1</td><td>297.58 ± 0.31</td></tr><tr><td>HR 2562 B</td><td>2016-02-28</td><td>J</td><td>620.0 ± 3.0</td><td>297.90 ± 0.25</td><td>620.2 ± 3.0</td><td>298.11 ± 0.32</td></tr><tr><td>HD 984 B</td><td>2015-08-30</td><td>H</td><td>216.3 ± 1.0</td><td>83.30 ± 0.30</td><td>216.2 ± 1.0</td><td>83.76 ± 0.30</td></tr><tr><td>HD 984 B</td><td>2015-08-30</td><td>J</td><td>217.9 ± 0.7</td><td>83.60 ± 0.20</td><td>217.8 ± 0.8</td><td>84.00 ± 0.21</td></tr></tbody></table>

to the measured PA of a companion accurately describes the effect of the change to the parallactic angle for each frame within a sequence. As the maximum change in Δ AVPARANG over a sequence was 0.05 deg, the effect on the companion astrometry is likely on this order or smaller. For the majority of cases, Δ AVPARANG changes by <100' of a degree over the course of a full observing sequence. The previous and revised astrometries for each published epoch are given in Table 5. We find small but not significant changes in the measured separations, and significant changes in the measured PAs due to the significant change in the north offset angle described in Sec. 7.

## 9 Discussion/Conclusion

We have identified and corrected several issues with the GPI DRP that affected astrometric measurements of both calibration binaries and substellar objects whose orbital motion was being monitored. We reprocessed the calibration data after implementing these fixes into the pipeline and revised the astrometric calibration of the instrument. The most significant change was to the
north offset angle; changing from −0.10 ± 0.13 deg to between 0.17 ± 0.14 deg and 0.45 ± 0.11 deg, depending on the date. The plate scale of the instrument was also remeasured as 14.161 ± 0.021 mas px⁻¹, consistent with the previous calibration albeit with a larger uncertainty.

Although the change to the astrometric calibration of the instrument is significant relative to the stated uncertainties, the impact should be limited to studies that combine GPI astrometry with that from instruments of similar precision. The revised calibration should not have a significant impact on the results and interpretation of studies that used GPI astrometry either solely or in conjunction with astrometry from instruments with significantly worse astrometric precision;⁶,⁹,¹⁰ an offset in the north angle will simply change the PA of the orbit on the sky (Ω). A more significant effect might be seen for orbit fits that combined astrometry from GPI with astrometry of a similar precision from other instruments.⁷,⁴⁴ The magnitude of the effect on the derived orbital parameters is likely small. All but one of the substellar companions studied with GPI have a small fraction of their complete orbits measured, and so the change of the shape of the posterior distributions describing the orbital elements is likely not statistically significant. The precision of astrometric measurements made with GPI is currently limited by measurement uncertainties except for widely separated companions such as the HR 8799 bcd, and the highest SNR measurements of β Pic b made in 2013 when the projected separation was ~430 mas, where the north angle uncertainty dominates the PA error budget. Lower SNR measurements of faint companions such as 51 Eri b are less affected, with the north angle uncertainty being between a factor of two and five smaller than the measurement uncertainty.

Future studies using archival GPI data will need to account for both the changes to the pipeline and the revision to the astrometric calibration. The updated pipeline is publicly available on the GPI instrument website⁴⁵ and on GitHub.⁴⁶ All users wishing to perform precision astrometry will have to reduce their data using the latest version of the pipeline, especially those obtained on the highlighted dates in Fig. 6, and apply the revised astrometric calibration presented in Sec. 7. The measurements presented here demonstrate the importance of continued astrometric calibration, especially for instruments on the Cassegrain mount of a telescope. Improvements to the limiting magnitude of GPI’s AO system as it is moved to Gemini North will allow us to use globular clusters as astrometric calibrations instead of isolated binaries, allowing for a more precise determination of the north angle via a comparison to both archival Hubble Space Telescope and contemporaneous Keck/NIRC2 observations.

This study also demonstrates the importance of precise and accurate astrometric calibration of instruments designed for high-contrast imaging of extrasolar planets. Instruments equipped with IFS necessarily have a small field of view, challenging for astrometric calibration that typically relies on images of globular clusters extending over several to tens of arcseconds. These results also demonstrate the importance of accounting for orbital motion, either between the two components of a calibration binary and/or the photocenter motion of one of the components if one of the components is itself a tight binary. A similar problem arises with the use of SiO masers near the Galactic Center;²⁸ the location of the infrared source is not necessarily coincident with that of the radio emission that the infrared astrometric reference frame is tied.⁴⁷ Precise and accurate astrometric calibration of future instruments with very narrow fields of view such as the Coronagraphic Instrument on the Wide Field Infrared Survey Telescope⁴⁸ will require a careful calibration strategy to mitigate the effects of these and other biases.

## Acknowledgments

The authors wish to thank Brian Chinn, Carlos Quiroz, Ignacio Arriagada, Thomas Hayward, and Carlos Alvarez for their useful discussions relating to this work. Supported by NSF under Grant Nos. AST-1411868 (R. D. R., E. L. N., K. B. F., B. M., and J. P.); AST-141378 (G. D.); AST-1518332 (R. D. R., J. J. W., T. M. E., J. R. G., and P. G. K.); and AST1411868 (J. H. and J. P.). Supported by NASA under Grant Nos. NNX14AJ80G (R. D. R., E. L. N., B. M., F. M., and M. P.); NSSC17K0535 (R. D. R, E. L. N., B. M., and J. B. R.); NNX15AC89G and NNX15AD95G (R. D. R., B. M., J. E. W., T. M. E., G. D., J. R. G., and P. G. K.). This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344. This work benefited from
NASA's Nexus for Exoplanet System Science (NExSS) research coordination network spon-
sored by NASA's Science Mission Directorate. J. R. was supported by the French National
Research Agency in the framework of the Investissements d'Avenir Program (Grant
No. ANR-15-IDEX-02). Based on observations obtained at the Gemini Observatory, which
is operated by the Association of Universities for Research in Astronomy, Inc., under a co-
operative agreement with the NSF on behalf of the Gemini partnership: the National Science
Foundation (United States), National Research Council (Canada), CONICYT (Chile), Ministerio
de Ciencia, Tecnología e Innovación Productiva (Argentina), Ministério da Ciência, Tecnologia e
Inovação (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). Some
of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a
scientific partnership among the California Institute of Technology, the University of California,
and the National Aeronautics and Space Administration. The observatory was made possible by
the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and
acknowledge the very significant cultural role and reverence that the summit of Maunakea has
always had within the indigenous Hawaiian community. We are most fortunate to have the
opportunity to conduct observations from this mountain. This work has made use of data from
the European Space Agency (ESA) mission Gaia,⁴⁹ processed by the Gaia Data Processing and
Analysis Consortium (DPAC).⁵⁰ Funding for the DPAC has been provided by the national
institutions, in particular the institutions participating in the Gaia Multilateral Agreement.
This research has made use of the SIMBAD database and the VizieR catalog access tool, both
operated at the CDS, Strasbourg, France. This research has made use of the Washington Double
Star Catalog maintained at the U.S. Naval Observatory. The authors have no relevant financial
interests and no other potential conflicts of interest to disclose.

References

1. B. A. Macintosh et al., “The Gemini Planet Imager: from science to design to construction,” *Proc. SPIE* **7015**, 701518 (2008).

2. B. Macintosh et al., “First light of the Gemini Planet Imager,” *Proc. Natl. Acad. Sci. U. S. A.* **111**, 12661–12666 (2014).

3. E. L. Nielsen et al., “The Gemini Planet Imager Exoplanet Survey: giant planet and brown dwarf demographics from 10 to 100 au,” *Astron. J.* **158**, 13 (2019).

4. Q. M. Konopacky et al., “Gemini Planet Imager observational calibrations V: astrometry and distortion,” *Proc. SPIE* **9147**, 914784 (2014).

5. J. J. Wang et al., “Dynamical constraints on the HR 8799 Planets with GPI,” *Astron. J.* **156**, 192 (2018).

6. M. A. Millar-Blanchaer et al., “β Pictoris’ inner disk in polarized light and new orbital parameters for β pictoris b,” *Astrophys. J.* **811**, 18 (2015).

7. J. Rameau et al., “Constraints on the architecture of the HD 95086 planetary system with the Gemini Planet Imager,” *Astrophys. J.* **822**, L29 (2016).

8. E. L. Nielsen et al., “The gemini planet imager exoplanet survey: dynamical mass of the exoplanet β Pictoris b from combined direct imaging and astrometry,” *Astron. J.* **159**(2) (2020).

9. J. J. Wang et al., “The orbit and transit prospects for β pictoris b constrained with one milliarcsecond astrometry,” *Astron. J.* **152**, 97 (2016).

10. R. J. De Rosa et al., “Astrometric confirmation and preliminary orbital parameters of the young exoplanet 51 Eridani b with the Gemini Planet Imager,” *Astrophys. J. Lett.* **814**, L3 (2015).

11. R. Soummer et al., “The Gemini Planet Imager coronagraph testbed,” *Proc. SPIE* **7440**, 74400R (2009).

12. R. Soummer et al., “Apodized pupil Lyot coronagraphs for arbitrary apertures. III. Quasi-achromatic solutions,” *Astrophys. J.* **729**, 144 (2011).

13. C. Marois et al., “Accurate astrometry and photometry of saturated and coronagraphic point spread functions,” *Astrophys. J.* **647**, 612–619 (2006).

14. A. Sivaramakrishnan and B. R. Oppenheimer, “Astrometry and photometry with coronagraphs,” *Astrophys. J.* **647**, 620–629 (2006).
15. J. K. Wallace et al., “The Gemini Planet Imager calibration wavefront sensor instrument,” *Proc. SPIE* **7736**, 77365D (2010).

16. J. K. Chilcote et al., “Performance of the integral field spectrograph for the Gemini Planet Imager,” *Proc. SPIE* **8446**, 84468W (2012).

17. J. E. Larkin et al., “The integral field spectrograph for the Gemini Planet Imager,” *Proc. SPIE* **9147**, 91471K (2014).

18. C. Marois et al., “Angular differential imaging: a powerful high-contrast imaging technique,” *Astrophys. J.* **641**, 556–564 (2006).

19. W. D. Pence et al., “Definition of the flexible image transport system (FITS), version 3.0,” *Astron. Astrophys.* **524**, A42 (2010).

20. A. M. Fowler and I. Gatley, “Noise reduction strategy for hybrid IR focal-plane arrays,” *Proc. SPIE* **1541**, 127–133 (1991).

21. I. S. McLean, *Electronic Imaging in Astronomy: Detectors and Instrumentation*, 2nd ed., Springer-Verlag, Berlin, Heidelberg (2008).

22. M. D. Perrin et al., “Gemini Planet Imager observational calibrations I: overview of the GPI data reduction pipeline,” *Proc. SPIE* **9147**, 91473J (2014).

23. M. D. Perrin et al., “Gemini Planet Imager observational calibrations XI: pipeline improvements and enhanced calibrations after two years on sky,” *Proc. SPIE* **9908**, 990837 (2016).

24. J. J. Wang et al., “The automated data processing architecture for the GPI exoplanet survey,” *Proc. SPIE* **10400**, 1040026 (2017).

25. D. Savransky et al., “Mining the GPIES database,” *Proc. SPIE* **10703**, 107030H (2018).

26. R. M. Green, *Spherical Astronomy*, p. 533, Cambridge University Press, Cambridge, New York (1985).

27. M. Shepherd, “Pointing the CBI,” 2001, http://www.astro.caltech.edu/~mcs/CBI/pointing/ (accessed 31 January 2020).

28. S. Yelda et al., “Improving Galactic Center Astrometry by reducing the effects of geometric distortion,” *Astrophys. J.* **725**, 331–352 (2010).

29. M. Service et al., “A new distortion solution for NIRC2 on the Keck II telescope,” *Publ. Astron. Soc. Pac.* **128**, 095004 (2016).

30. S. A. Metchev and L. A. Hillenbrand, “The Palomar/Keck adaptive optics survey of young solar analogs: evidence for a universal companion mass function,” *Astron. Astrophys. Suppl. Ser.* **181**, 62–109 (2009).

31. J. Lu, “nirc2_distortion,” 2016, https://github.com/jluastro/nirc2_distortion (accessed 31 January 2020).

32. J. Gubler and D. Tytler, “Differential atmospheric refraction and limitations on the relative astrometric accuracy of large telescopes,” *Publ. Astron. Soc. Pac.* **110**, 738–746 (1998).

33. E. Diolaiti et al., “Analysis of isoplanatic high resolution stellar fields by the StarFinder code,” *Astron. Astrophys. Suppl. Ser.* **147**, 335–346 (2000).

34. D. Foreman-Mackey et al., “emcee: the MCMC Hammer,” *Publ. Astron. Soc. Pac.* **125**, 306–312 (2013).

35. R. J. De Rosa et al., “An updated visual orbit of the directly imaged Exoplanet 51 Eridani b and prospects for a dynamical mass measurement with Gaia,” *Astron. J.* **159**, 1 (2020).

36. F. van Leeuwen, “Validation of the new Hipparcos reduction,” *Astron. Astrophys.* **474**, 653–664 (2007).

37. Gaia Collaboration et al., “GaiaData release 2,” *Astron. Astrophys.* **616**, A1 (2018).

38. M. J. Reid et al., “Trigonometric parallaxes of high mass star forming regions: the structure and kinematics of the Milky Way,” *Astrophys. J.* **783**, 130 (2014).

39. D. Pourbaix et al., “SB9: the ninth catalogue of spectroscopic binary orbits,” *Astron. Astrophys.* **424**, 727–732 (2004).

40. D. Erickson et al., “Integrated modeling for the GPi flexure sensitive structure,” *Proc. SPIE* **7017**, 701711 (2008).

41. Ward-Duong et al., “Gemini planet imager spectroscopy of the dusty substellar companion HD 206893 B,” *Astron. J.* submitted (2020).

42. Q. M. Konopacky et al., “Discovery of a substellar companion to the nearby debris disk host HR 2562,” *Astrophys. J.* **829**, L4 (2016).
43. M. Johnson-Groh et al., “Integral field spectroscopy of the low-mass companion HD 984 B with the Gemini Planet Imager,” *Astrophys. J.* **153**, 190 (2017).

44. A. L. Maire et al., “Hint of curvature in the orbital motion of the exoplanet 51 Eridani b using 3 yr of VLT/SPHERE monitoring,” *Astron. Astrophys.* **624**, A118 (2019).

45. GPI Collaboration, “GPI Data Pipeline Documentation,” 2017, http://docs.planetimager.org/pipeline/ (accessed 31 January 2020).

46. GPI Collaboration, “gpi_pipeLine,” 2019, https://github.com/geminiplanetimager/gpi_pipeLine/ (accessed 31 January 2020).

47. S. Sakai et al., “The Galactic Center: an improved astrometric reference frame for stellar orbits around the supermassive black hole,” *Astrophys. J.* **873**, 65 (2019).

48. D. Spergel et al., “Wide-field infrared survey telescope-astrophysics focused telescope assets WFIRST-AFTA 2015 report,” arXiv:1503.03757 (2015).

49. Gaia Consortium, “Gaia,” 2020, https://www.cosmos.esa.int/web/gaia (accessed 31 January 2020).

50. Gaia Consortium, “The DPAC Consortium,” 2020, https://www.cosmos.esa.int/web/gaia/dpac/consortium (accessed 31 January 2020).

Biographies of the authors are not available.