Engineering Papers⌕ Search

Engineering topics

Doyon, Rene

Publications and source records attributed to Doyon, Rene.

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

The Gemini Planet Imager View of the HD 32297 Debris Disk

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

79 ASTRONOMY AND ASTROPHYSICS↗

Revised astrometric calibration of the Gemini Planet Imager

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.

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

An Updated Visual Orbit of the Directly Imaged Exoplanet 51 Eridani b and Prospects for a Dynamical Mass Measurement with Gaia

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

79 ASTRONOMY AND ASTROPHYSICS↗

Detection of a Low-mass Stellar Companion to the Accelerating A2IV Star HR 1645

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

79 ASTRONOMY AND ASTROPHYSICS↗

Asymmetries in adaptive optics point spread functions

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

79 ASTRONOMY AND ASTROPHYSICS↗

An Exo–Kuiper Belt with an Extended Halo around HD 191089 in Scattered Light

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

79 ASTRONOMY AND ASTROPHYSICS↗

GPI Spectra of HR8799 C, D, and E in H-K Bands with KLIP Forward Modeling

We demonstrate KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR8799, using PyKLIP. We report new and re-reduced spectrophotometry of HR8799 c, d, and e from H-K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting fake sources and recovering them over a range of parameters. The K1/K2 spectra for planets c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR8799e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We compare planets c, d, and e with M, L, and T-type field objects. All objects are consistent with low gravity mid-to-late L dwarfs, however, a lack of standard spectra for low gravity late L-type objects lead to poor fit for gravity. We place our results in context of atmospheric models presented in previous publications and discuss differences in the spectra of the three planets.

planets and satellites: gaseous planets { stars: i↗

Spectroscopy of luminous infrared galaxies at 2 microns: 1. The ultraluminous galaxies (L(sub IR) approximately greater than 10 (exp 12) solar luminosity)

We present high-quality spectra covering the K window at a resolving power of 340 for a sample of 13 ultraluminous (L(sub IR) approximately greater than 10(exp 12) solar luminosity) infrared-selected galaxies, and line fluxes for a comparison sample of 24 lower luminosity galaxies. The 2 micrometers spectra of 10 of the ultraluminous galaxies are characterized by emission and absorption features commonly associated with stars and star formation; two others have the red power-law spectra and Br gamma line widths of Seyfert 1 galaxies; the final galaxy has strong emission from hot dust. We have found no broad-line active nuclei not already known from optical observations, despite the fact that the extinction at 2 micrometers is 1/10 that at optical wavelengths; any putative Seyfert 1 nuclei must be deeply buried. Powerful continua and emission lines from H2 and Br gamma are detected in all the ultraluminous galaxies. Comparing the H2 1-0 S(1), Br gamma, and 2 micrometers and far-infrared luminosities to those of the lower luminosity galaxies yields several major results. First, the dereddened Br gamma emission, relative to the far-infrared luminosity is significantly depressed in the ultraluminous sample, when compared to the lower luminosity galaxies. Five of the ultraluminous galaxies have L(sub Br gamma)L(sub IR) ratios lower than for any of the comparison objects. Second, the H2 1-0 S(1) luminosity is also responsible, directly or indirectly, for producing the excited H2, and that the H2 apparently comes from optically thin regions in both classes of objects. Third, eight of the 13 ultraluminous systems have lower 2 micrometers/far-infrared luminosity ratios than any of the lower luminosity galaxies, and five of these are the galaxies also deficient in Br gamma. These three findings may be understood if the the H2, Br gamma, and 2 mircometers continua in the ultraluminous galaxies arise from spatially distinct regions, with the continuum and Br gamma largely coming from volumes optically thick even at 2 micrometers, and obscured in such a fashion that the extinctions measured using optical spectroscopy do not properly measure the true optical depths. If this is the case, then even near-infrared spectroscopy may be unable to exclude the presence of undetected powerful active galactive nuclei in the ultraluminous galaxies.

Goldader, Jeffrey D.↗

The He I 2.06 microns/Br-gamma ratio in starburst galaxies - An objective constraint on the upper mass limit to the initial mass function

The use of the He I 2.06 microns/Br-gamma ratio as a constraint on the massive stellar population in star-forming galaxies is developed. A theoretical relationship between the He I 2.06 microns/Br-gamma ratio and the effective temperature of the exciting star in H II regions is derived. The effects of collisional excitation and dust within the nebula on the ratio are also considered. It is shown that the He I 2.06 microns/Br-gamma ratio is a steep function of the effective temperature, a property which can be used to determine the upper mass limit of the initial mass function (IMF) in galaxies. This technique is reliable for upper mass limits less than about 40 solar masses. New near-infrared spectra of starburst galaxies are presented. The He I 2.06 microns/Br-gamma ratios observed imply a range of upper mass limits from 27 to over 40 solar masses. There is also evidence that the upper mass limit is spatially dependent within a given galaxy. These results suggest that the upper mass limit is not a uniquely defined parameter of the IMF and probably varies with local physical conditions.

Doyon, Rene↗