Far-infrared & Sub-millimeter Studies of Circumstellar Disks

Circumstellar disks are critical structures in the star and planet formation processes, as they provide a conduit to channel material onto the central object and supply a reservoir of dust and gas to form planets. This thesis focuses on the far-infrared, and sub-millimeter observations of circumstellar disks at two key evolutionary phases; pri- mordial proto-planetary disks, and evolved debris disks – remnants of a system that has undergone a degree of planet formation. Four individual studies of circumstellar disks are presented in this thesis. The results of a 97% complete census of far-infrared emission measured with the Herschel Space Observatory, targeting stars of spectral types M4 and later, in the Taurus molec- ular cloud are presented. This census is the first large-scale survey sensitive to emission across the stellar and sub-stellar boundary. Results from an initial test grid of model spectral energy distributions, generated with the radiative transfer code MCFOST, show that 73% of the observed Class II population are constrained by canonical disks that are viewed from face-on to edge-on inclinations. Sub-millimeter observations with the Caltech Sub-millimeter Observatory are presented for an association of young T-Tauri stars in the Aquila star-forming region. The results of disk frequency and disk mass of this complete census are investigated in this extremely low stellar density environment. A sub-millimeter investigation for two populations of candidate debris disk; warm and cold excess disks is presented. None of the candidate disks were detected in the sub-millimeter despite upper-limits below that expected, based on blackbody model fits to excesses at shorter wavelengths. Several scenarios are investigated in order to identify the null detection rate, such as stellar multiplicity and background-source contamination. Finally, a partially resolved sub-millimeter map of the debris disk around the HR 8799 multiple planetary system is presented. The planet formation history of the system is investigated through the witnessed morphology of the emission.Science & Technology Facilities Counci

The Peculiar Debris Disk of HD 111520 as Resolved by the Gemini Planet Imager

Using the Gemini Planet Imager (GPI), we have resolved the circumstellar debris disk around HD 111520 at a projected range of ~30-100 AU in both total and polarized $H$-band intensity. The disk is seen edge-on at a position angle of ~165$^{\circ}$ along the spine of emission. A slight inclination or asymmetric warping are covariant and alters the interpretation of the observed disk emission. We employ 3 point spread function (PSF) subtraction methods to reduce the stellar glare and instrumental artifacts to confirm that there is a roughly 2:1 brightness asymmetry between the NW and SE extension. This specific feature makes HD 111520 the most extreme examples of asymmetric debris disks observed in scattered light among similar highly inclined systems, such as HD 15115 and HD 106906. We further identify a tentative localized brightness enhancement and scale height enhancement associated with the disk at ~40 AU away from the star on the SE extension. We also find that the fractional polarization rises from 10 to 40% from 0.5" to 0.8" from the star. The combination of large brightness asymmetry and symmetric polarization fraction leads us to believe that an azimuthal dust density variation is causing the observed asymmetry.Comment: 9 pages, 8 Figures, 1 table, Accepted to Ap

HD 165054: An Astrometric Calibration Field for High-Contrast Imagers in Baade\u27s Window

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

Automated data processing architecture for the Gemini Planet Imager Exoplanet Survey

The Gemini Planet Imager Exoplanet Survey (GPIES) is a multi-year direct imaging survey of 600 stars to discover and characterize young Jovian exoplanets and their environments. We have developed an automated data architecture to process and index all data related to the survey uniformly. An automated and flexible data processing framework, which we term the Data Cruncher, combines multiple data reduction pipelines together to process all spectroscopic, polarimetric, and calibration data taken with GPIES. With no human intervention, fully reduced and calibrated data products are available less than an hour after the data are taken to expedite follow-up on potential objects of interest. The Data Cruncher can run on a supercomputer to reprocess all GPIES data in a single day as improvements are made to our data reduction pipelines. A backend MySQL database indexes all files, which are synced to the cloud, and a front-end web server allows for easy browsing of all files associated with GPIES. To help observers, quicklook displays show reduced data as they are processed in real-time, and chatbots on Slack post observing information as well as reduced data products. Together, the GPIES automated data processing architecture reduces our workload, provides real-time data reduction, optimizes our observing strategy, and maintains a homogeneously reduced dataset to study planet occurrence and instrument performance.Comment: 21 pages, 3 figures, accepted in JATI

Dynamical Mass Measurement of the Young Spectroscopic Binary V343 Normae AaAb Resolved With the Gemini Planet Imager

We present new spatially resolved astrometry and photometry from the Gemini Planet Imager of the inner binary of the young multiple star system V343 Normae, which is a member of the beta Pictoris moving group. V343 Normae comprises a K0 and mid-M star in a ~4.5 year orbit (AaAb) and a wide 10" M5 companion (B). By combining these data with archival astrometry and radial velocities we fit the orbit and measure individual masses for both components of M_Aa = 1.10 +/- 0.10 M_sun and M_Ab = 0.290 +/- 0.018 M_sun. Comparing to theoretical isochrones, we find good agreement for the measured masses and JHK band magnitudes of the two components consistent with the age of the beta Pic moving group. We derive a model-dependent age for the beta Pic moving group of 26 +/- 3 Myr by combining our results for V343 Normae with literature measurements for GJ 3305, which is another group member with resolved binary components and dynamical masses.Comment: 12 pages, 7 figures. Accepted to A

Improving and Assessing Planet Sensitivity of the GPI Exoplanet Survey with a Forward Model Matched Filter

We present a new matched filter algorithm for direct detection of point sources in the immediate vicinity of bright stars. The stellar Point Spread Function (PSF) is first subtracted using a Karhunen-Lo\'eve Image Processing (KLIP) algorithm with Angular and Spectral Differential Imaging (ADI and SDI). The KLIP-induced distortion of the astrophysical signal is included in the matched filter template by computing a forward model of the PSF at every position in the image. To optimize the performance of the algorithm, we conduct extensive planet injection and recovery tests and tune the exoplanet spectra template and KLIP reduction aggressiveness to maximize the Signal-to-Noise Ratio (SNR) of the recovered planets. We show that only two spectral templates are necessary to recover any young Jovian exoplanets with minimal SNR loss. We also developed a complete pipeline for the automated detection of point source candidates, the calculation of Receiver Operating Characteristics (ROC), false positives based contrast curves, and completeness contours. We process in a uniform manner more than 330 datasets from the Gemini Planet Imager Exoplanet Survey (GPIES) and assess GPI typical sensitivity as a function of the star and the hypothetical companion spectral type. This work allows for the first time a comparison of different detection algorithms at a survey scale accounting for both planet completeness and false positive rate. We show that the new forward model matched filter allows the detection of $50\%$ fainter objects than a conventional cross-correlation technique with a Gaussian PSF template for the same false positive rate.Comment: ApJ accepte

The MALATANG Survey : The L GAS-L IR Correlation on Sub-kiloparsec Scale in Six Nearby Star-forming Galaxies as Traced by HCN J = 4 → 3 and HCO + J = 4 → 3

This is an author-created, un-copyedited version of an article published in The Astrophysical Journal. The Version of Record is available online at https://doi.org/10.3847/1538-4357/aac512.We present HCN J = 4→3 and HCO+ J = 4→3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2×2 region at 14 (FWHM) resolution (corresponding to linear scales of ∼0.2-1.0 kpc). The LIR-Ldense relation, where the dense gas is traced by the HCN J = 4→3 and the HCO+ J = 4→3 emission, measured in our sample of spatially resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, LIR/Ldense, shows systematic variations with LIR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between LIR/Ldense ratio and the warm-dust temperature gauged by the 70 μm/100 μm flux ratio. We find that the luminosity ratios of IR/HCN (4-3) and IR/HCO+ (4-3), which can be taken as a proxy for the star formation efficiency (SFE) in the dense molecular gas (SFE dense), appear to be nearly independent of the dense gas fraction ( f dense) for our sample of galaxies. The SFE of the total molecular gas (SFEmol) is found to increase substantially with f dense when combining our data with those on local (ultra)luminous infrared galaxies and high-z quasars. The mean LHCN(4-3) LHCO+(4-3) line ratio measured for the six targeted galaxies is 0.9±0.6. No significant correlation is found for the L'HCN(4-3) L'HCO+(4-3) ratio with the star formation rate as traced by L IR, nor with the warm-dust temperature, for the different populations of galaxies.Peer reviewe

GPI spectra of HR 8799 c, d, and e from 1.5 to 2.4μ\mum with KLIP Forward Modeling

We explore KLIP forward modeling spectral extraction on Gemini Planet Imager coronagraphic data of HR 8799, using PyKLIP and show algorithm stability with varying KLIP parameters. We report new and re-reduced spectrophotometry of HR 8799 c, d, and e in H & K bands. We discuss a strategy for choosing optimal KLIP PSF subtraction parameters by injecting simulated sources and recovering them over a range of parameters. The K1/K2 spectra for HR 8799 c and d are similar to previously published results from the same dataset. We also present a K band spectrum of HR 8799 e for the first time and show that our H-band spectra agree well with previously published spectra from the VLT/SPHERE instrument. We show that HR 8799 c and d show significant differences in their H & K spectra, but do not find any conclusive differences between d and e or c and e, likely due to large error bars in the recovered spectrum of e. Compared to M, L, and T-type field brown dwarfs, all three planets are most consistent with mid and late L spectral types. All objects are consistent with low gravity but a lack of standard spectra for low gravity limit the ability to fit the best spectral type. We discuss how dedicated modeling efforts can better fit HR 8799 planets' near-IR flux and discuss how differences between the properties of these planets can be further explored.Comment: Accepted to AJ, 25 pages, 16 Figure

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

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