We have obtained Gemini Planet Imager (GPI) J-, H-, K1-, and K2-Spec observations of the iconic debris ring around the young, main-sequence star HR 4796A. We applied several point-spread function (PSF) subtraction techniques to the observations (Mask-and-Interpolate, RDI-NMF, RDI-KLIP, and ADI-KLIP) to measure the geometric parameters and the scattering phase function for the disk. To understand the systematic errors associated with PSF subtraction, we also forward-modeled the observations using a Markov Chain Monte Carlo framework and a simple model for the disk. We found that measurements of the disk geometric parameters were robust, with all of our analyses yielding consistent results; however, measurements of the scattering phase function were challenging to reconstruct from PSF-subtracted images, despite extensive testing. As a result, we estimated the scattering phase function using disk modeling. We searched for a dependence of the scattering phase function with respect to the GPI filters but found none. We compared the H-band scattering phase function with that measured by Hubble Space Telescope STIS at visual wavelengths and discovered a blue color at small scattering angles and a red color at large scattering angles, consistent with predictions and laboratory measurements of large grains. Finally, we successfully modeled the SPHERE H2 HR 4796A scattered phase function using a distribution of hollow spheres composed of silicates, carbon, and metallic iron.

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

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

Context. Searching for high-redshift galaxies is a field of intense activity in modern observational cosmology that will continue to grow with future ground-based and sky observatories. Over the last few years, a lot has been learned about the high-z Universe. Aims. Despite extensive Ly-alpha Blobs (LAB) surveys from low to high redshifts, giant LABs over 100 kpc have been found mostly at z~2-4. This redshift range is coincident with the transition epoch of galactic gas-circulation processes from inflows to outflows at z~2.5-3. This suggests that the formation of giant LABs may be related to a combination of gas inflows and outflows. Their extreme youth makes them interesting objects in the study of galaxy formation as they provide insight into some of the youngest known highly star forming galaxies, with only modest time investments using ground-based telescopes. Methods. Systematic narrow-band Ly-alpha nebula surveys are ongoing, but they are limited in their covered redshift range and their comoving volume. This poses a significant problem when searching for such rare sources. To address this problem, we developed a systematic searching tool, ATACAMA (A Tool for seArChing for lArge LyMan Alpha nebulae) designed to find large Ly-alpha nebulae at any redshift within deep multi-wavelength broad-band imaging. Results. We identified a Ly-alpha nebula candidate at zphot~3.3 covering an isophotal area of 29.4sq.arcsec. Its morphology shows a bright core and a faint core which coincides with the morphology of previously known Ly-alpha blobs. A first estimation of the Ly-alpha equivalent width and line flux agree with the values from the study led by several groups.

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

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

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

We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$\alpha$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$\alpha$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$\alpha$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}\alpha}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$\alpha$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $\phi^*=10^{-2.85}$ Mpc$^{-3}$, $\alpha = -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $\phi^*=10^{-2.8}$ Mpc$^{-3}$, $\alpha=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $\phi^*=10^{-2.91}$ Mpc$^{-3}$, $\alpha = -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$\alpha$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $\rho_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $\rho_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $\rho_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$\alpha$ surveys.

While most of the inter-galactic medium (IGM) today is permeated by ionized hydrogen, it was largely filled with neutral hydrogen for the first 700 million years after the Big Bang. The process that ionized the IGM (cosmic reionization) is expected to be spatially inhomogeneous, with fainter galaxies playing a significant role. However, we still have only a few direct constraints on the reionization process. Here we report the first spectroscopic confirmation of two galaxies and very likely a third galaxy in a group (hereafter EGS77) at redshift z = 7.7, merely 680 Myrs after the Big Bang. The physical separation among the three members is < 0.7 Mpc. We estimate the radius of ionized bubble of the brightest galaxy to be about 1.02 Mpc, and show that the individual ionized bubbles formed by all three galaxies likely overlap significantly, forming a large yet localized ionized region, which leads to the spatial inhomogeneity in the reionization process. It is striking that two of three galaxies in EGS77 are quite faint in the continuum, thanks to our selection of reionizing sources using their Lyman-alpha line emission. Indeed, one is the faintest spectroscopically confirmed galaxy yet discovered at such high redshifts. Our observations provide direct constraints in the process of cosmic reionization, and allow us to investigate the properties of sources responsible for reionizing the universe.

We present a catalog of 208 $0.3 < z < 2.1$ Emission Line Galaxies (ELG) selected from 1D slitless spectroscopy obtained using Hubble's WFC3 G102 grism, as part of the Faint Infrared Grism Survey (FIGS). We identify ELG candidates by searching for significant peaks in all continuum-subtracted G102 spectra, and, where possible, confirm candidates by identifying consistent emission lines in other available spectra or with published spectroscopic redshifts. We provide derived emission line fluxes and errors, redshifts, and equivalent widths (EW) for H$\alpha$ $\lambda6563$, [OIII]$\lambda\lambda4959,5007$, and [OII]$\lambda\lambda3727$ emission lines, for emission line galaxies down to AB(F105W) $ > 28$ and $> 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ line flux. We use the resulting line catalog to investigate a possible relationship between line emission and a galaxy's environment. We use 7th-nearest-neighbor distances to investigate the typical surroundings of ELGs compared to non-ELGs, and we find that [OIII] emitters are preferentially found at intermediate galaxy densities near galaxy groups. We characterize these ELGs in terms of the galaxy specific star formation rate (SSFR) versus stellar mass, and find no significant influence of environment on that relation. We calculate star formation rates (SFR), and find no dependence of SFR on local galaxy surface density for $0.3 < z < 0.8$ H$\alpha$ emitters and for $0.8<z<1.3$ [OIII] emitters. We find similar rates of close-pair interaction between ELGs and non-ELGs. For galaxy surface densities $\Sigma \leq 30$ Mpc$^{-2}$, we find no consistent effect of environment on star formation.

We present the largest-ever sample of 79 Ly$\alpha$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$\alpha$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly\alpha}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$\alpha$ luminosity function.

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

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

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

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

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

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

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

Aims. We seek is to identify old and massive galaxies at 0.5<z<2.1 on the basis of the magnesium index MgUV and then study their physical properties. We computed the MgUV index based on the best spectral fitting template of $\sim$3700 galaxies using data from the VLT VIMOS Deep Survey (VVDS) and VIMOS Ultra Deep Survey (VUDS) galaxy redshift surveys. Based on galaxies with the largest signal to noise and the best fit spectra we selected 103 objects with the highest spectral MgUV signature. We performed an independent fit of the photometric data of these galaxies and computed their stellar masses, star formation rates, extinction by dust and age, and we related these quantities to the MgUV index. We find that the MgUV index is a suitable tracer of early-type galaxies at an advanced stage of evolution. Selecting galaxies with the highest MgUV index allows us to choose the most massive, passive, and oldest galaxies at any epoch. The formation epoch t_f computed from the fitted age as a function of the total mass in stars supports the downsizing formation paradigm in which galaxies with the highest mass formed most of their stars at an earlier epoch.

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

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

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.

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

We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency.

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

We derive direct measurement gas-phase metallicities of $7.4 < 12 + \log(O/H) < 8.4$ for 14 low-mass Emission Line Galaxies (ELGs) at $0.3 < z < 0.8$ identified in the Faint Infrared Grism Survey (FIGS). We use deep slitless G102 grism spectroscopy of the Hubble Ultra Deep Field (HUDF), dispersing light from all objects in the field at wavelengths between 0.85 and 1.15 microns. We run an automatic search routine on these spectra to robustly identify 71 emission line sources, using archival data from VLT/MUSE to measure additional lines and confirm redshifts. We identify 14 objects with $0.3 < z < 0.8$ with measurable O[III]$\lambda$4363 \AA\ emission lines in matching VLT/MUSE spectra. For these galaxies, we derive direct electron-temperature gas-phase metallicities with a range of $7.4 < 12 + \log(O/H) < 8.4$. With matching stellar masses in the range of $10^{7.9} M_{\odot} < M_{\star} < 10^{10.4} M_{\odot}$, we construct a mass-metallicity (MZ) relation and find that the relation is offset to lower metallicities compared to metallicities derived from alternative methods (e.g.,$R_{23}$, O3N2, N2O2) and continuum selected samples. Using star formation rates (SFR) derived from the $H\alpha$ emission line, we calculate our galaxies' position on the Fundamental Metallicity Relation (FMR), where we also find an offset toward lower metallicities. This demonstrates that this emission-line-selected sample probes objects of low stellar masses but even lower metallicities than many comparable surveys. We detect a trend suggesting galaxies with higher Specific Star Formation (SSFR) are more likely to have lower metallicity. This could be due to cold accretion of metal-poor gas that drives star formation, or could be because outflows of metal-rich stellar winds and SNe ejecta are more common in galaxies with higher SSFR.

(Abr.) Laser guide stars employed at astronomical observatories provide artificial wavefront reference sources to help correct (in part) the impact of atmospheric turbulence on astrophysical observations. Following the recent commissioning of the 4 Laser Guide Star Facility (4LGSF) on UT4 at the VLT, we characterize the spectral signature of the uplink beams from the 22W lasers to assess the impact of laser scattering from the 4LGSF on science observations. We use the MUSE optical integral field spectrograph to acquire spectra at a resolution of R~3000 of the uplink laser beams over the wavelength range of 4750\AA\ to 9350\AA. We report the first detection of laser-induced Raman scattering by N2, O2, CO2, H2O and (tentatively) CH4 molecules in the atmosphere above the astronomical observatory of Cerro Paranal. In particular, our observations reveal the characteristic spectral signature of laser photons -- but 480\AA\ to 2210\AA\ redder than the original laser wavelength of 5889.959\AA\ -- landing on the 8.2m primary mirror of UT4 after being Raman-scattered on their way up to the sodium layer. Laser-induced Raman scattering is not unique to the observatory of Cerro Paranal, but common to any astronomical telescope employing a laser-guide-star (LGS) system. It is thus essential for any optical spectrograph coupled to a LGS system to handle thoroughly the possibility of a Raman spectral contamination via a proper baffling of the instrument and suitable calibrations procedures. These considerations are particularly applicable for the HARMONI optical spectrograph on the upcoming Extremely Large Telescope. At sites hosting multiple telescopes, laser collision prediction tools also ought to account for the presence of Raman emission from the uplink laser beam(s) to avoid the unintentional contamination of observations acquired with telescopes in the vicinity of a LGS system.

The first observations of laser guide star photons Raman-scattered by air molecules above the Very Large Telescope (VLT) were reported in June 2017. The initial detection came from the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph, following the installation of the 4 Laser Guide Star Facility (4LGSF) on the Unit Telescope 4 (UT4) of the VLT. In this Letter, we delve further into the symbiotic relationship between the 4LGSF laser guide star system, the UT4 telescope, and MUSE by monitoring the spectral contamination of MUSE observations by Raman photons over a 27 month period. This dataset reveals that dust particles deposited on the primary and tertiary mirrors of UT4 -- responsible for a reflectivity loss of ~8% at 6000{\AA} -- contribute (60$\pm5)% to the laser line fluxes detected by MUSE. The flux of Raman lines, contaminating scientific observations acquired with optical spectrographs, thus provides a new, non-invasive means to monitor the evolving scatter properties of the mirrors of astronomical telescopes equipped with laser guide star systems.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

We present new high resolution imaging of a light-scattering dust ring and halo around the young star HD 35841. Using spectroscopic and polarimetric data from the Gemini Planet Imager in H-band (1.6 microns), we detect the highly inclined (i=85 deg) ring of debris down to a projected separation of ~12 au (~0.12") for the first time. Optical imaging from HST/STIS shows a smooth dust halo extending outward from the ring to >140 au (>1.4"). We measure the ring's scattering phase function and polarization fraction over scattering angles of 22-125 deg, showing a preference for forward scattering and a polarization fraction that peaks at ~30% near the ansae. Modeling of the scattered-light disk indicates that the ring spans radii of ~60-220 au, has a vertical thickness similar to that of other resolved dust rings, and contains grains as small as 1.5 microns in diameter. These models also suggest the grains have a low porosity, are more likely to consist of carbon than astrosilicates, and contain significant water ice. The halo has a surface brightness profile consistent with that expected from grains pushed by radiation pressure from the main ring onto highly eccentric but still bound orbits. We also briefly investigate arrangements of a possible inner disk component implied by our spectral energy distribution models, and speculate about the limitations of Mie theory for doing detailed analyses of debris disk dust populations.

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.

We improve the accuracy of photometric redshifts by including low-resolution spectral data from the G102 grism on the Hubble Space Telescope, which assists in redshift determination by further constraining the shape of the broadband Spectral Energy Disribution (SED) and identifying spectral features. The photometry used in the redshift fits includes near-IR photometry from FIGS+CANDELS, as well as optical data from ground-based surveys and HST ACS, and mid-IR data from Spitzer. We calculated the redshifts through the comparison of measured photometry with template galaxy models, using the EAZY photometric redshift code. For objects with F105W $< 26.5$ AB mag with a redshift range of $0 < z < 6$, we find a typical error of $\Delta z = 0.03 * (1+z)$ for the purely photometric redshifts; with the addition of FIGS spectra, these become $\Delta z = 0.02 * (1+z)$, an improvement of 50\%. Addition of grism data also reduces the outlier rate from 8\% to 7\% across all fields. With the more-accurate spectrophotometric redshifts (SPZs), we searched the FIGS fields for galaxy overdensities. We identified 24 overdensities across the 4 fields. The strongest overdensity, matching a spectroscopically identified cluster at $z=0.85$, has 28 potential member galaxies, of which 8 have previous spectroscopic confirmation, and features a corresponding X-ray signal. Another corresponding to a cluster at $z=1.84$ has 22 members, 18 of which are spectroscopically confirmed. Additionally, we find 4 overdensities that are detected at an equal or higher significance in at least one metric to the two confirmed clusters.

The Gemini/HST Galaxy Cluster Project (GCP) covers 14 z=0.2-1.0 clusters with X-ray luminosity of L_500 >= 10^44 ergs/s in the 0.1-2.4 keV band. In this paper we provide homogeneously calibrated X-ray luminosities, masses and radii, and we present the complete catalog of the ground-based photometry for the GCP clusters. The clusters were observed with Gemini North or South in three or four of the optical passbands g', r', i' and z'. The photometric catalog includes consistently calibrated total magnitudes, colors, and geometrical parameters. The photometry reaches ~25 mag in the passband closest to rest frame B-band. We summarize comparisons of our photometry with data from the Sloan Digital Sky Survey. We describe the sample selection for our spectroscopic observations, and establish the calibrations to obtain rest frame magnitudes and colors. Finally, we derive the color-magnitude relations for the clusters and briefly discuss these in the context of evolution with redshift. Consistent with our results based on spectroscopic data, the color-magnitude relations support passive evolution of the red-sequence galaxies. The absence of change in the slope with redshift, constrains the allowable age variation along the red sequence to <0.05 dex between the brightest cluster galaxies and those four magnitudes fainter. The paper serves as the main reference for the GCP cluster and galaxy selection, X-ray data and ground-based photometry.

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.

We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf.

We constrain the quasar contribution to cosmic reionization based on our deep optical survey of z~6 quasars down to z_R=24.15 using Subaru/Suprime-Cam in three UKIDSS-DXS fields covering 6.5 deg^2. In Kashikawa et al. (2015), we select 17 quasar candidates and report our initial discovery of two low-luminosity quasars (M_1450~ -23) from seven targets, one of which might be a Lyman alpha emitting galaxy. From an additional optical spectroscopy, none of the four candidates out of the remaining ten turn out to be genuine quasars. Moreover, the deeper optical photometry provided by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) shows that, unlike the two already-known quasars, the i-z and z-y colors of the last six candidates are consistent with M- or L-type brown dwarfs. Therefore, the quasar luminosity function (QLF) in the previous paper is confirmed. Compiling QLF measurements from the literature over a wide magnitude range, including an extremely faint AGN candidate from Parsa et al. (2017}, to fit them with a double power-law, we find that the best-fit faint-end slope is alpha=-2.04^+0.33_-0.18 (-1.98^+0.48_-0.21) and characteristic magnitude is M_1450^*=-25.8^+1.1_-1.9 (-25.7^+1.0_-1.8) in the case of two (one) quasar detection. Our result suggests that, if the QLF is integrated down to M_1450=-18, quasars produce ~1-12% of the ionizing photons required to ionize the whole universe at z~6 with 2sigma confidence level, assuming that the escape fraction is f_esc=1 and the IGM clumpy factor is C=3. Even when the systematic uncertainties are taken into account, our result supports the scenario that quasars are the minor contributors of reionization.

Narrowband imaging is a highly successful approach for finding large numbers of high redshift Lya emitting galaxies (LAEs) up to z~6.6. However, at z>~7 there are as yet only 3 narrowband selected LAEs with spectroscopic confirmations (two at z~6.9-7.0, one at z~7.3), which hinders extensive studies on cosmic reionization and galaxy evolution at this key epoch. We have selected 23 candidate z~6.9 LAEs in COSMOS field with the large area narrowband survey LAGER (Lyman-Alpha Galaxies at the End of Reionization). In this work we present spectroscopic followup observations of 12 candidates using IMACS on Magellan. For 9 of these, the observations are sufficiently deep to detect the expected lines. Lya emission lines are identified in six sources (yielding a success rate of 2/3), including 3 luminous LAEs with Lya luminosities of L(Lya) ~ 10^{43.5} erg/s, the highest among known spectroscopically confirmed galaxies at >~7.0. This triples the sample size of spectroscopically confirmed narrowband selected LAEs at z>~7, and confirms the bright end bump in the Lya luminosity function we previously derived based on the photometric sample, supporting a patchy reionization scenario. Two luminous LAEs appear physically linked with projected distance of 1.1 pMpc and velocity difference of ~ 170 km/s. They likely sit in a common ionized bubble produced by themselves or with close neighbors, which reduces the IGM attenuation of Lya. A tentative narrow NV${\lambda}$1240 line is seen in one source, hinting at activity of a central massive black hole with metal rich line emitting gas.

The Faint Infrared Grism Survey (FIGS) is a deep Hubble Space Telescope (HST) WFC3/IR (Wide Field Camera 3 Infrared) slitless spectroscopic survey of four deep fields. Two fields are located in the Great Observatories Origins Deep Survey-North (GOODS-N) area and two fields are located in the Great Observatories Origins Deep Survey-South (GOODS-S) area. One of the southern fields selected is the Hubble Ultra Deep Field. Each of these four fields were observed using the WFC3/G102 grism (0.8$\mu m$-1.15$\mu m$ continuous coverage) with a total exposure time of 40 orbits (~ 100 kilo-seconds) per field. This reaches a 3 sigma continuum depth of ~26 AB magnitudes and probes emission lines to $\approx 10^{-17}\ erg\ s^{-1} \ cm^{-2}$. This paper details the four FIGS fields and the overall observational strategy of the project. A detailed description of the Simulation Based Extraction (SBE) method used to extract and combine over 10000 spectra of over 2000 distinct sources brighter than m_F105W=26.5 mag is provided. High fidelity simulations of the observations is shown to significantly improve the background subtraction process, the spectral contamination estimates, and the final flux calibration. This allows for the combination of multiple spectra to produce a final high quality, deep, 1D-spectra for each object in the survey.

We present the first results from the ongoing LAGER project (Lyman Alpha Galaxies in the Epoch of Reionization), which is the largest narrowband survey for $z \sim$ 7 galaxies to date. Using a specially built narrowband filter NB964 for the superb large-area Dark-Energy Camera (DECam) on the NOAO/CTIO 4m Blanco telescope, LAGER has collected 34 hours NB964 narrowband imaging data in the 3 deg$^2$ COSMOS field. We have identified 23 Lyman Alpha Emitter (LAE) candidates at $z$ = 6.9 in the central 2-deg$^2$ region, where DECam and public COSMOS multi-band images exist. The resulting luminosity function can be described as a Schechter function modified by a significant excess at the bright end (4 galaxies with $L_{Ly\alpha} \sim $ 10$^{43.4\pm0.2}$ erg s$^{-1}$). The number density at $L_{Ly\alpha}\sim$ 10$^{43.4\pm0.2}$ erg s$^{-1}$ is little changed from z= 6.6, while at fainter $L_{Ly\alpha}$ it is substantially reduced. Overall, we see a fourfold reduction in Ly$\alpha$ luminosity density from $z$ = 5.7 to 6.9. Combined with a more modest evolution of the continuum UV luminosity density, this suggests a factor of $\sim 3$ suppression of Ly$\alpha$ by radiative transfer through the $z \sim$ 7 intergalactic medium (IGM). It indicates an IGM neutral fraction $x_{HI}$ $\sim$ 0.4--0.6 (assuming Ly$\alpha$ velocity offsets of 100-200 km s$^{-1}$). The changing shape of the Ly$\alpha$ luminosity function between $z\lesssim 6.6$ and $z=6.9$ supports the hypothesis of ionized bubbles in a patchy reionization at $z\sim$ 7.

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.

We present spectro-photometry spanning 1-5 $\mu$m of 51 Eridani b, a 2-10 M$_\text{Jup}$ planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new $K1$ (1.90-2.19 $\mu$m) and $K2$ (2.10-2.40 $\mu$m) spectra taken with the Gemini Planet Imager as well as an updated $L_P$ (3.76 $\mu$m) and new $M_S$ (4.67 $\mu$m) photometry from the NIRC2 Narrow camera. The new data were combined with $J$ (1.13-1.35 $\mu$m) and $H$ (1.50-1.80 $\mu$m) spectra from the discovery epoch with the goal of better characterizing the planet properties. 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4-T8) and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud free and patchy/intermediate opacity clouds. Model fits suggest that 51 Eri b has an effective temperature ranging between 605-737 K, a solar metallicity, a surface gravity of $\log$(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the SED. From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 ($\log L/L_{\odot}$), leaving 51 Eri b in the unique position as being one of the only directly imaged planet consistent with having formed via cold-start scenario. Comparisons of the planet SED against warm-start models indicates that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M$_{\oplus}$.

We present H band spectroscopic and Halpha photometric observations of HD 100546 obtained with GPI and MagAO. We detect H band emission at the location of the protoplanet HD 100546b, but show that choice of data processing parameters strongly affects the morphology of this source. It appears point-like in some aggressive reductions, but rejoins an extended disk structure in the majority of the others. Furthermore, we demonstrate that this emission appears stationary on a timescale of 4.6 yrs, inconsistent at the 2sigma level with a Keplerian clockwise orbit at 59 au in the disk plane. The H band spectrum of the emission is inconsistent with any type of low effective temperature object or accreting protoplanetary disk. It strongly suggests a scattered light origin, as it is consistent with the spectrum of the star and the spectra extracted at other locations in the disk. A non detection at the 5sigma level of HD 100546b in differential Halpha imaging places an upper limit, assuming the protoplanet lies in a gap free of extinction, on the accretion luminosity and accretion rate of 1.7E-4 Lsun and MMdot<6.4E-7Mjup^2/yr for 1Rjup. These limits are comparable to the accretion luminosity and rate of TTauri-stars or LkCa 15b. Taken together, these lines of evidence suggest that the H band source at the location of HD 100546b is not emitted by a planetary photosphere or an accreting circumplanetary disk but is a disk feature enhanced by the PSF subtraction process. This non-detection is consistent with the non-detection in the K band reported in an earlier study but does not exclude the possibility that HD 100546b is deeply embedded.

We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at H{\alpha}. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.

Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.

We present the first scattered-light images of the debris disk around 49 ceti, a ~40 Myr A1 main sequence star at 59 pc, famous for hosting two massive dust belts as well as large quantities of atomic and molecular gas. The outer disk is revealed in reprocessed archival Hubble Space Telescope NICMOS F110W images, as well as new coronagraphic H band images from the Very Large Telescope SPHERE instrument. The disk extends from 1.1" (65 AU) to 4.6" (250 AU), and is seen at an inclination of 73degr, which refines previous measurements at lower angular resolution. We also report no companion detection larger than 3 M_Jup at projected separations beyond 20 AU from the star (0.34"). Comparison between the F110W and H-band images is consistent with a grey color of 49 ceti's dust, indicating grains larger than >2microns. Our photometric measurements indicate a scattering efficiency / infrared excess ratio of 0.2-0.4, relatively low compared to other characterized debris disks. We find that 49 ceti presents morphological and scattering properties very similar to the gas-rich HD 131835 system. From our constraint on the disk inclination we find that the atomic gas previously detected in absorption must extend to the inner disk, and that the latter must be depleted of CO gas. Building on previous studies, we propose a schematic view of the system describing the dust and gas structure around 49 ceti and hypothetic scenarios for the gas nature and origin.

Uncovering the ingredients and the architecture of planetary systems is a very active field of research that has fuelled many new theories on giant planet formation, migration, composition, and interaction with the circumstellar environment. We aim at discovering and studying new such systems, to further expand our knowledge of how low-mass companions form and evolve. We obtained high-contrast H-band images of the circumstellar environment of the F5V star HD206893, known to host a debris disc never detected in scattered light. These observations are part of the SPHERE High Angular Resolution Debris Disc Survey (SHARDDS) using the InfraRed Dual-band Imager and Spectrograph (IRDIS) installed on VLT/SPHERE. We report the detection of a source with a contrast of 3.6x10^{-5} in the H-band, orbiting at a projected separation of 270 milliarcsecond or 10 au, corresponding to a mass in the range 24 to 73 Mjup for an age of the system in the range 0.2 to 2 Gyr. The detection was confirmed ten months later with VLT/NaCo, ruling out a background object with no proper motion. A faint extended emission compatible with the disc scattered light signal is also observed. The detection of a low-mass companion inside a massive debris disc makes this system an analog of other young planetary systems such as beta Pictoris, HR8799 or HD95086 and requires now further characterisation of both components to understand their interactions.

Galaxies at high redshifts provide a valuable tool to study cosmic dawn, and therefore it is crucial to reliably identify these galaxies. Here, we present an unambiguous and first simultaneous detection of both the Lyman-alpha emission and the Lyman break from a z = 7.512+/- 0.004 galaxy, observed in the Faint Infrared Grism Survey (FIGS). These spectra, taken with G102 grism on Hubble Space Telescope (HST), show a significant emission line detection (6 sigma) in multiple observational position angles (PA), with total integrated Ly{\alpha} line flux of 1.06+/- 0.12 e10-17erg s-1cm-2. The line flux is nearly a factor of four higher than the previous MOSFIRE spectroscopic observations of faint Ly{\alpha} emission at {\lambda} = 1.0347{\mu}m, yielding z = 7.5078+/- 0.0004. This is consistent with other recent observations implying that ground-based near-infrared spectroscopy underestimates total emission line fluxes, and if confirmed, can have strong implications for reionization studies that are based on ground-based Lyman-{\alpha} measurements. A 4-{\sigma} detection of the NV line in one PA also suggests a weak Active Galactic Nucleus (AGN), potentially making this source the highest-redshift AGN yet found. Thus, this observation from the Hubble Space Telescope clearly demonstrates the sensitivity of the FIGS survey, and the capability of grism spectroscopy to study the epoch of reionization.

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.

We present $H$-band near-infrared polarimetric imaging observations of the F5V star HD~157587 obtained with the Gemini Planet Imager (GPI) that reveal the debris disk as a bright ring structure at a separation of $\sim$80$-$100~AU. The new GPI data complement recent HST/STIS observations that show the disk extending out to over 500~AU. The GPI image displays a strong asymmetry along the projected minor axis as well as a fainter asymmetry along the projected major axis. We associate the minor and major axis asymmetries with polarized forward scattering and a possible stellocentric offset, respectively. To constrain the disk geometry we fit two separate disk models to the polarized image, each using a different scattering phase function. Both models favor a disk inclination of $\sim 70\degr$ and a $1.5\pm0.6$ AU stellar offset in the plane of the sky along the projected major axis of the disk. We find that the stellar offset in the disk plane, perpendicular to the projected major axis is degenerate with the form of the scattering phase function and remains poorly constrained. The disk is not recovered in total intensity due in part to strong adaptive optics residuals, but we recover three point sources. Considering the system's proximity to the galactic plane and the point sources' positions relative to the disk, we consider it likely that they are background objects and unrelated to the disk's offset from the star.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

A principal scientific goal of the Gemini Planet Imager (GPI) is obtaining milliarcsecond astrometry to constrain exoplanet orbits. However, astrometry of directly imaged exoplanets is subject to biases, systematic errors, and speckle noise. Here we describe an analytical procedure to forward model the signal of an exoplanet that accounts for both the observing strategy (angular and spectral differential imaging) and the data reduction method (Karhunen-Lo\`eve Image Projection algorithm). We use this forward model to measure the position of an exoplanet in a Bayesian framework employing Gaussian processes and Markov chain Monte Carlo (MCMC) to account for correlated noise. In the case of GPI data on $\beta$ Pic b, this technique, which we call Bayesian KLIP-FM Astrometry (BKA), outperforms previous techniques and yields 1$\sigma$-errors at or below the one milliarcsecond level. We validate BKA by fitting a Keplerian orbit to twelve GPI observations along with previous astrometry from other instruments. The statistical properties of the residuals confirm that BKA is accurate and correctly estimates astrometric errors. Our constraints on the orbit of $\beta$ Pic b firmly rule out the possibility of a transit of the planet at 10-$\sigma$ significance. However, we confirm that the Hill sphere of $\beta$ Pic b will transit, giving us a rare chance to probe the circumplanetary environment of a young, evolving exoplanet. We provide an ephemeris for photometric monitoring of the Hill sphere transit event, which will begin at the start of April in 2017 and finish at the end of January in 2018.

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

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.

We present astrometric monitoring of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager between 2013 and 2016. A small but significant position angle change is detected at constant separation; the orbital motion is confirmed with literature measurements. Efficient Monte Carlo techniques place preliminary constraints on the orbital parameters of HD 95086 b. With 68% confidence, a semimajor axis of 61.7^{+20.7}_{-8.4} au and an inclination of 153.0^{+9.7}_{-13.5} deg are favored, with eccentricity less than 0.21. Under the assumption of a co-planar planet-disk system, the periastron of HD 95086 b is beyond 51 au with 68% confidence. Therefore HD 95086 b cannot carve the entire gap inferred from the measured infrared excess in the SED of HD 95086. We use our sensitivity to additional planets to discuss specific scenarios presented in the literature to explain the geometry of the debris belts. We suggest that either two planets on moderately eccentric orbits or three to four planets with inhomogeneous masses and orbital properties are possible. The sensitivity to additional planetary companions within the observations presented in this study can be used to help further constrain future dynamical simulations of the planet-disk system.

During the commissioning of the Gemini MCAO System (GeMS), we had the opportunity to obtain data with the Gemini Multi-Object Spectrograph (GMOS), the most utilised instrument at Gemini South Observatory, in March and May 2012. Several globular clusters were observed in imaging mode that allowed us to study the performance of this new and untested combination. GMOS is a visible instrument, hence pushing MCAO toward the visible.We report here on the results with the GMOS instruments, derive photometric performance in term of Full Width Half Maximum (FWHM) and throughput. In most of the cases, we obtained an improvement factor of at least 2 against the natural seeing. This result also depends on the Natural Guide Star constellation selected for the observations and we then study the impact of the guide star selection on the FWHM performance.We also derive a first astrometric analysis showing that the GeMS+GMOS system provide an absolute astrometric precision better than 8mas and a relative astrometric precision lower than 50 mas.

We present new $H$ (1.5-1.8 $\mu$m) photometric and $K_1$ (1.9-2.2 $\mu$m) spectroscopic observations of the young exoplanet HD 95086 b obtained with the Gemini Planet Imager. The $H$-band magnitude has been significantly improved relative to previous measurements, whereas the low resolution $K_1$ ($\lambda/\delta\lambda \approx 66$) spectrum is featureless within the measurement uncertainties, and presents a monotonically increasing pseudo-continuum consistent with a cloudy atmosphere. By combining these new measurements with literature $L^{\prime}$ photometry, we compare the spectral energy distribution of the planet to other young planetary-mass companions, field brown dwarfs, and to the predictions of grids of model atmospheres. HD 95086 b is over a magnitude redder in $K_1-L^{\prime}$ color than 2MASS J12073346-3932539 b and HR 8799 c and d, despite having a similar $L^{\prime}$ magnitude. Considering only the near-infrared measurements, HD 95086 b is most analogous to the brown dwarfs 2MASS J2244316+204343 and 2MASS J21481633+4003594, both of which are thought to have dusty atmospheres. Morphologically, the spectral energy distribution of HD 95086 b is best fit by low temperature ($T_{\rm eff} =$ 800-1300 K), low surface gravity spectra from models which simulate high photospheric dust content. This range of effective temperatures is consistent with field L/T transition objects, but the spectral type of HD 95086 b is poorly constrained between early L and late T due to its unusual position the color-magnitude diagram, demonstrating the difficulty in spectral typing young, low surface gravity substellar objects. As one of the reddest such objects, HD 95086 b represents an important empirical benchmark against which our current understanding of the atmospheric properties of young extrasolar planets can be tested.

We present the Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole (SEP) field. The large area covered (7.7 deg$^2$), together with one of the lowest Galactic cirrus emissions in the entire sky and a very extensive coverage by Spitzer, Herschel, Akari, and GALEX, make the SIMES field ideal for extragalactic studies. The elongated geometry of the SIMES area ($\approx$4:1), allowing for a significant cosmic variance reduction, further improves the quality of statistical studies in this field. Here we present the reduction and photometric measurements of the Spitzer/IRAC data. The survey reaches a depth of 1.93 and 1.75 $\mu$Jy (1$\sigma$) at 3.6 and 4.5 $\mu$m, respectively. We discuss the multiwavelength IRAC--based catalog, completed with optical, mid-- and far--IR observations. We detect 341,000 sources with F$_{3.6\mu m} \geq 3\sigma$. Of these, 10% have an associated 24 $\mu$m counterpart, while 2.7% have an associated SPIRE source. We release the catalog through the NASA/IPAC Infrared Science Archive (IRSA). Two scientific applications of these IRAC data are presented in this paper: first we compute integral number counts at 3.6 $\mu$m. Second, we use the [3.6]--[4.5] color index to identify galaxy clusters at z$>$1.3. We select 27 clusters in the full area, a result consistent with previous studies at similar depth.

We present H and K band imaging polarimetry for the PDS 66 circumstellar disk obtained during the commissioning of the Gemini Planet Imager (GPI). Polarization images reveal a clear detection of the disk in to the 0.12'' inner working angle (IWA) in H band, almost 3 times as close to the star as the previous HST observations with NICMOS and STIS (0.35'' effective IWA). The centro-symmetric polarization vectors confirm that the bright inner disk detection is due to circumstellar scattered light. A more diffuse disk extends to a bright outer ring centered at 80 AU. We discuss several physical mechanisms capable of producing the observed ring + gap structure. GPI data confirm enhanced scattering on the East side of the disk which is inferred to be nearer to us. We also detect a lateral asymmetry in the South possibly due to shadowing from material within the inner working angle. This likely corresponds to a temporally variable azimuthal asymmetry observed in HST/STIS coronagraphic imaging.

We present the first scattered-light image of the debris disk around HD 131835 in $H$ band using the Gemini Planet Imager. HD 131835 is a $\sim$15 Myr old A2IV star at a distance of $\sim$120 pc in the Sco-Cen OB association. We detect the disk only in polarized light and place an upper limit on the peak total intensity. No point sources resembling exoplanets were identified. Compared to its mid-infrared thermal emission, the disk in scattered light shows similar orientation but different morphology. The scattered-light disk extends from $\sim$75 to $\sim$210 AU in the disk plane with roughly flat surface density. Our Monte Carlo radiative transfer model can well describe the observations with a model disk composed of a mixture of silicates and amorphous carbon. In addition to the obvious brightness asymmetry due to stronger forward scattering, we discover a weak brightness asymmetry along the major axis with the northeast side being 1.3 times brighter than the southwest side at a 3-{\sigma} level.

We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$\deg$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet.

We present $H$-band observations of $\beta$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $\beta$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($\phi=85.27{\deg}^{+0.26}_{-0.19}$) with a position angle $\theta_{PA}=30.35{\deg}^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $\theta_{PA}\approx29{\deg}$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $\beta$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $\Omega=31.75{\deg}\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $\beta$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $\beta$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $\beta$ Pic b; 2) the mutual inclination of the inner disk and $\beta$ Pic $b$ is $4{\deg}$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $\beta$ Pic $b$ or self-stirring by the disk's parent bodies.

We present Gemini Planet Imager (GPI) observations of AU Microscopii, a young M dwarf with an edge-on, dusty debris disk. Integral field spectroscopy and broadband imaging polarimetry were obtained during the commissioning of GPI. In our broadband imaging polarimetry observations, we detect the disk only in total intensity and find asymmetries in the morphology of the disk between the southeast and northwest sides. The southeast side of the disk exhibits a bump at 1$''$ (10 AU projected separation) that is three times more vertically extended and three times fainter in peak surface brightness than the northwest side at similar separations. This part of the disk is also vertically offset by 69$\pm$30 mas to the northeast at 1$''$ when compared to the established disk mid-plane and consistent with prior ALMA and Hubble Space Telescope/STIS observations. We see hints that the southeast bump might be a result of detecting a horizontal sliver feature above the main disk that could be the disk backside. Alternatively when including the morphology of the northwest side, where the disk mid-plane is offset in the opposite direction $\sim$50 mas between 0$.''$4 and 1$.''$2, the asymmetries suggest a warp-like feature. Using our integral field spectroscopy data to search for planets, we are 50% complete for $\sim$4 $M_\mathrm{Jup}$ planets at 4 AU. We detect a source, resolved only along the disk plane, that could either be a candidate planetary mass companion or a compact clump in the disk.

Context. Lyman Break Galaxies (LBG) and Narrow Band (NB) surveys have been successful at detecting large samples of high-redshift galaxies. Both methods are subject to contamination from low-redshift interlopers. Aims. In this paper, our aim is to investigate the nature of low-redshift interlopers in NB Lyman-$\alpha$ emitters (LAE) searches. Methods. From previous HAWK-I NB imaging at z $\sim$7.7 we identify three objects that would have been selected as high-redshift LAEs had our optical data been one magnitude shallower (but still one to two magnitudes fainter than the near infrared data). We follow-up these objects in spectroscopy with XSHOOTER at the VLT. Results. Despite low quality data due to bad weather conditions, for each of the three objects we identify one, and only one emission line, in the spectra of the objects, that we identify as the O[III]5007A line. This result combined to spectral energy density fitting and tests based on line ratios of several populations of galaxies we infer that the 3 objects are ultrastrong line emitters at redshifts $\sim$1.1. Conclusions. From this work and the literature we remark that the O[III] line appears to be a common source of contamination in high-redshift LBG and LAE samples and we suggest that efforts be put to characterize with high accuracy the O[III] luminosity function out to redshift $\sim$3 or higher.

We present the discovery of one or two extremely faint z~6 quasars in 6.5 deg^2 utilizing a unique capability of the wide-field imaging of the Subaru/Suprime-Cam. The quasar selection was made in (i'-z_B) and (z_B-z_R) colors, where z_B and z_R are bandpasses with central wavelengths of 8842A and 9841A, respectively. The color selection can effectively isolate quasars at z~6 from M/L/T dwarfs without the J-band photometry down to z_R<24.0, which is 3.5 mag. deeper than SDSS. We have selected 17 promising quasar candidates. The follow-up spectroscopy for seven targets identified one apparent quasar at z=6.156 with M_1450=-23.10. We also identified one possible quasar at z=6.041 with a faint continuum of M_1450=-22.58 and a narrow Lyman-alpha emission with HWHM=427 km/s, which cannot be distinguished from Lyman-alpha emitters. We derive the quasar luminosity function at z~6 by combining our faint quasar sample with the bright quasar samples by SDSS and CFHQS. Including our data points invokes a higher number density in the faintest bin of the quasar luminosity function than the previous estimate employed. This suggests a steeper faint-end slope than lower-z, though it is yet uncertain based on a small number of spectroscopically identified faint quasars and several quasar candidates are still remain to be diagnosed. The steepening of the quasar luminosity function at the faint-end does increase the expected emission rate of the ionizing photon, however, it only changes by a factor of ~2-6. This was found to be still insufficient for the required photon budget of reionization at z~6.

During the first-light run of the Gemini Planet Imager (GPI) we obtained K-band spectra of exoplanets HR 8799 c and d. Analysis of the spectra indicates that planet d may be warmer than planet c. Comparisons to recent patchy cloud models and previously obtained observations over multiple wavelengths confirm that thick clouds combined with horizontal variation in the cloud cover generally reproduce the planets' spectral energy distributions. When combined with the 3 to 4 um photometric data points, the observations provide strong constraints on the atmospheric methane content for both planets. The data also provide further evidence that future modeling efforts must include cloud opacity, possibly including cloud holes, disequilibrium chemistry, and super-solar metallicity.

The Gemini Planet Imager (GPI) entered on-sky commissioning and had its first-light at the Gemini South (GS) telescope in November 2013. GPI is an extreme adaptive optics (XAO), high-contrast imager and integral-field spectrograph dedicated to the direct detection of hot exo-planets down to a Jupiter mass. The performance of the apodized pupil Lyot coronagraph depends critically upon the residual wavefront error (design goal of 60 nm RMS with 5 mas RMS tip/tilt), and therefore is most sensitive to vibration (internal or external) of Gemini's instrument suite. Excess vibration can be mitigated by a variety of methods such as passive or active dampening at the instrument or telescope structure or Kalman filtering of specific frequencies with the AO control loop. Understanding the sources, magnitudes and impact of vibration is key to mitigation. This paper gives an overview of related investigations based on instrument data (GPI AO module) as well as external data from accelerometer sensors placed at different locations on the GS telescope structure. We report the status of related mitigation efforts, and present corresponding results.

The Gemini Planet Imager (GPI) entered on-sky commissioning phase, and had its First Light at the Gemini South telescope in November 2013. Meanwhile, the fast loops for atmospheric correction of the Extreme Adaptive Optics (XAO) system have been closed on many dozen stars at different magnitudes (I=4-8), elevation angles and a variety of seeing conditions, and a stable loop performance was achieved from the beginning. Ultimate contrast performance requires a very low residual wavefront error (design goal 60 nm RMS), and optimization of the planet finding instrument on different ends has just begun to deepen and widen its dark hole region. Laboratory raw contrast benchmarks are in the order of 10^-6 or smaller. In the telescope environment and in standard operations new challenges are faced (changing gravity, temperature, vibrations) that are tackled by a variety of techniques such as Kalman filtering, open-loop models to keep alignment to within 5 mas, speckle nulling, and a calibration unit (CAL). The CAL unit was especially designed by the Jet Propulsion Laboratory to control slowly varying wavefront errors at the focal plane of the apodized Lyot coronagraph by the means of two wavefront sensors: 1) a 7x7 low order Shack-Hartmann SH wavefront sensor (LOWFS), and 2) a special Mach-Zehnder interferometer for mid-order spatial frequencies (HOWFS) - atypical in that the beam is split in the focal plane via a pinhole but recombined in the pupil plane with a beamsplitter. The original design goal aimed for sensing and correcting on a level of a few nm which is extremely challenging in a telescope environment. This paper focuses on non-common path low order wavefront correction as achieved through the CAL unit on sky. We will present the obtained results as well as explain challenges that we are facing.

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

We present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses a new integral field polarimetry architecture to provide high contrast linear polarimetry with minimal systematic biases between the orthogonal polarizations. We describe the design, data reduction methods, and performance of polarimetry with GPI. Point spread function subtraction via differential polarimetry suppresses unpolarized starlight by a factor of over 100, and provides sensitivity to circumstellar dust reaching the photon noise limit for these observations. In the case of the circumstellar disk around HR 4796A, GPI's advanced adaptive optics system reveals the disk clearly even prior to PSF subtraction. In polarized light, the disk is seen all the way in to its semi-minor axis for the first time. The disk exhibits surprisingly strong asymmetry in polarized intensity, with the west side >9 times brighter than the east side despite the fact that the east side is slightly brighter in total intensity. Based on a synthesis of the total and polarized intensities, we now believe that the west side is closer to us, contrary to most prior interpretations. Forward scattering by relatively large silicate dust particles leads to the strong polarized intensity on the west side, and the ring must be slightly optically thick in order to explain the lower brightness in total intensity there. These findings suggest that the ring is geometrically narrow and dynamically cold, perhaps shepherded by larger bodies in the same manner as Saturn's F ring.

The Gemini Planet Imager (GPI) has as its science instrument an infrared integral field spectrograph/polarimeter (IFS). Integral field spectrographs are scientificially powerful but require sophisticated data reduction systems. For GPI to achieve its scientific goals of exoplanet and disk characterization, IFS data must be reconstructed into high quality astrometrically and photometrically accurate datacubes in both spectral and polarization modes, via flexible software that is usable by the broad Gemini community. The data reduction pipeline developed by the GPI instrument team to meet these needs is now publicly available following GPI's commissioning. This paper, the first of a series, provides a broad overview of GPI data reduction, summarizes key steps, and presents the overall software framework and implementation. Subsequent papers describe in more detail the algorithms necessary for calibrating GPI data. The GPI data reduction pipeline is open source, available from planetimager.org, and will continue to be enhanced throughout the life of the instrument. It implements an extensive suite of task primitives that can be assembled into reduction recipes to produce calibrated datasets ready for scientific analysis. Angular, spectral, and polarimetric differential imaging are supported. Graphical tools automate the production and editing of recipes, an integrated calibration database manages reference files, and an interactive data viewer customized for high contrast imaging allows for exploration and manipulation of data.

The Gemini Planet Imager instrument's adaptive optics (AO) subsystem was designed specifically to facilitate high-contrast imaging. It features several new technologies, including computationally efficient wavefront reconstruction with the Fourier transform, modal gain optimization every 8 seconds, and the spatially filtered wavefront sensor. It also uses a Linear-Quadratic-Gaussian (LQG) controller (aka Kalman filter) for both pointing and focus. We present on-sky performance results from verification and commissioning runs from December 2013 through May 2014. The efficient reconstruction and modal gain optimization are working as designed. The LQG controllers effectively notch out vibrations. The spatial filter can remove aliases, but we typically use it oversized by about 60% due to stability problems.

The Gemini Planet Imager (GPI) is a new facility, extreme adaptive optics (AO), coronagraphic instrument, currently being integrated onto the 8-meter Gemini South telescope, with the ultimate goal of directly imaging extrasolar planets. To achieve the contrast required for the desired science, it is necessary to quantify and mitigate wavefront error (WFE). A large source of potential static WFE arises from the primary AO wavefront sensor (WFS) detector's use of multiple readout segments with independent signal chains including on-chip preamplifiers and external amplifiers. Temperature changes within GPI's electronics cause drifts in readout segments' bias levels, inducing an RMS WFE of 1.1 nm and 41.9 nm over 4.44 degrees Celsius, for magnitude 4 and 11 stars, respectively. With a goal of $<$2 nm of static WFE, these are significant enough to require remedial action. Simulations imply a requirement to take fresh WFS darks every 2 degrees Celsius of temperature change, for a magnitude 6 star; similarly, for a magnitude 7 star, every 1 degree Celsius of temperature change. For sufficiently dim stars, bias drifts exceed the signal, causing a large initial WFE, and the former periodic requirement practically becomes an instantaneous/continuous one, making the goal of $<$2 nm of static WFE very difficult for stars of magnitude 9 or fainter. In extreme cases, this can cause the AO loops to destabilize due to perceived nonphysical wavefronts, as some of the WFS's Shack-Hartmann quadcells are split between multiple readout segments. Presented here is GPI's AO WFS geometry, along with detailed steps in the simulation used to quantify bias drift related WFE, followed by laboratory and on sky results, and concluded with possible methods of remediation.

An Atmospheric Dispersion Corrector (ADC) uses a double-prism arrangement to nullify the vertical chromatic dispersion introduced by the atmosphere at non-zero zenith distances. The ADC installed in the Gemini Planet Imager (GPI) was first tested in August 2012 while the instrument was in the laboratory. GPI was installed at the Gemini South telescope in August 2013 and first light occurred later that year on November 11th. In this paper, we give an overview of the characterizations and performance of this ADC unit obtained in the laboratory and on sky, as well as the structure of its control software.

The Gemini Planet Imager (GPI) is a complex optical system designed to directly detect the self-emission of young planets within two arcseconds of their host stars. After suppressing the starlight with an advanced AO system and apodized coronagraph, the dominant residual contamination in the focal plane are speckles from the atmosphere and optical surfaces. Since speckles are diffractive in nature their positions in the field are strongly wavelength dependent, while an actual companion planet will remain at fixed separation. By comparing multiple images at different wavelengths taken simultaneously, we can freeze the speckle pattern and extract the planet light adding an order of magnitude of contrast. To achieve a bandpass of 20%, sufficient to perform speckle suppression, and to observe the entire two arcsecond field of view at diffraction limited sampling, we designed and built an integral field spectrograph with extremely low wavefront error and almost no chromatic aberration. The spectrograph is fully cryogenic and operates in the wavelength range 1 to 2.4 microns with five selectable filters. A prism is used to produce a spectral resolution of 45 in the primary detection band and maintain high throughput. Based on the OSIRIS spectrograph at Keck, we selected to use a lenslet-based spectrograph to achieve an rms wavefront error of approximately 25 nm. Over 36,000 spectra are taken simultaneously and reassembled into image cubes that have roughly 192x192 spatial elements and contain between 11 and 20 spectral channels. The primary dispersion prism can be replaced with a Wollaston prism for dual polarization measurements. The spectrograph also has a pupil-viewing mode for alignment and calibration.

We present on-sky polarimetric observations with the Gemini Planet Imager (GPI) obtained at straight Cassegrain focus on the Gemini South 8-m telescope. Observations of polarimetric calibrator stars, ranging from nearly unpolarized to strongly polarized, enable determination of the combined telescope and instrumental polarization. We find the conversion of Stokes $I$ to linear and circular instrumental polarization in the instrument frame to be $I \rightarrow (Q_{\rm IP}, U_{\rm IP}, P_{\rm IP}, V_{\rm IP}) = (-0.037 \pm 0.010\%, +0.4338 \pm 0.0075\%, 0.4354 \pm 0.0075\%, -6.64 \pm 0.56\%)$. Such precise measurement of instrumental polarization enables $\sim 0.1\%$ absolute accuracy in measurements of linear polarization, which together with GPI's high contrast will allow GPI to explore scattered light from circumstellar disk in unprecedented detail, conduct observations of a range of other astronomical bodies, and potentially even study polarized thermal emission from young exoplanets. Observations of unpolarized standard stars also let us quantify how well GPI's differential polarimetry mode can suppress the stellar PSF halo. We show that GPI polarimetry achieves cancellation of unpolarized starlight by factors of 100-200, reaching the photon noise limit for sensitivity to circumstellar scattered light for all but the smallest separations at which the calibration for instrumental polarization currently sets the limit.

The Gemini Planet Imager (GPI) is a next-generation, facility instrument currently being commissioned at the Gemini South observatory. GPI combines an extreme adaptive optics system and integral field spectrograph (IFS) with an apodized-pupil Lyot coronagraph (APLC) producing an unprecedented capability for directly imaging and spectroscopically characterizing extrasolar planets. GPI's operating goal of $10^{-7}$ contrast requires very precise alignments between the various elements of the coronagraph (two pupil masks and one focal plane mask) and active control of the beam path throughout the instrument. Here, we describe the techniques used to automatically align GPI and maintain the alignment throughout the course of science observations. We discuss the particular challenges of maintaining precision alignments on a Cassegrain mounted instrument and strategies that we have developed that allow GPI to achieve high contrast even in poor seeing conditions.

The Gemini Planet Imager is an extreme AO instrument with an integral field spectrograph (IFS) operating in Y, J, H, and K bands. Both the Gemini telescope and the GPI instrument are very complex systems. Our goal is that the combined telescope and instrument system may be run by one observer operating the instrument, and one operator controlling the telescope and the acquisition of light to the instrument. This requires a smooth integration between the two systems and easily operated control interfaces. We discuss the definition of the software and hardware interfaces, their implementation and testing, and the integration of the instrument with the telescope environment.

The Gemini Planet Imager (GPI) Extreme Adaptive Optics Coronograph contains an interferometric mode: a 10-hole non-redundant mask (NRM) in its pupil wheel. GPI operates at $Y, J, H$, and $K$ bands, using an integral field unit spectrograph (IFS) to obtain spectral data at every image pixel. NRM on GPI is capable of imaging with a half resolution element inner working angle at moderate contrast, probing the region behind the coronagraphic spot. The fine features of the NRM PSF can provide a reliable check on the plate scale, while also acting as an attenuator for spectral standard calibrators that would otherwise saturate the full pupil. NRM commissioning data provides details about wavefront error in the optics as well as operations of adaptive optics control without pointing control from the calibration system. We compare lab and on-sky results to evaluate systematic instrument properties and examine the stability data in consecutive exposures. We discuss early on-sky performance, comparing images from integration and tests with the first on-sky images, and demonstrate resolving a known binary. We discuss the status of NRM and implications for future science with this mode.

The Gemini Planet Imager (GPI) combines extreme adaptive optics, an integral field spectrograph, and a high performance coronagraph to directly image extrasolar planets in the near-infrared. Because the coronagraph blocks most of the light from the star, it prevents the properties of the host star from being measured directly. Instead, satellite spots, which are created by diffraction from a square grid in the pupil plane, can be used to locate the star and extract its spectrum. We describe the techniques implemented into the GPI Data Reduction Pipeline to measure the properties of the satellite spots and discuss the precision of the reconstructed astrometry and spectrophotometry of the occulted star. We find the astrometric precision of the satellite spots in an $H$-band datacube to be $0.05$ pixels and is best when individual satellite spots have a signal to noise ratio (SNR) of $> 20$. In regards to satellite spot spectrophotometry, we find that the total flux from the satellite spots is stable to $\sim 7\%$ and scales linearly with central star brightness and that the shape of the satellite spot spectrum varies on the $2\%$ level.

We present the wavelength calibration for the lenslet-based Integral Field Spectrograph (IFS) that serves as the science instrument for the Gemini Planet Imager (GPI). The GPI IFS features a 2.7" x 2.7" field of view and a 190 x 190 lenslet array (14.3 mas/lenslet) operating in $Y$, $J$, $H$, and $K$ bands with spectral resolving power ranging from $R$ $\sim$ 35 to 78. Due to variations across the field of view, a unique wavelength solution is determined for each lenslet characterized by a two-dimensional position, the spectral dispersion, and the rotation of the spectrum with respect to the detector axes. The four free parameters are fit using a constrained Levenberg-Marquardt least-squares minimization algorithm, which compares an individual lenslet's arc lamp spectrum to a simulated arc lamp spectrum. This method enables measurement of spectral positions to better than 1/10th of a pixel on the GPI IFS detector using Gemini's facility calibration lamp unit GCAL, improving spectral extraction accuracy compared to earlier approaches. Using such wavelength calibrations we have measured how internal flexure of the spectrograph with changing zenith angle shifts spectra on the detector. We describe the methods used to compensate for these shifts when assembling datacubes from on-sky observations using GPI.

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

The Gemini Planet Imager (GPI) is currently in its final Acceptance & Testing stages. GPI is an XAO system based on a tweeter & woofer architecture (43 & 9 actuators respectively across the pupil), with the tweeter being a Boston Michromachines $64^2$ MEMS device. The XAO AO system is tightly integrated with a Lyot apodizing coronagraph. Acceptance testing started in February 2013 at the University of California, Santa Cruz. A conclusive acceptance review was held in July 2013 and the instrument was found ready for shipment to the Gemini South telescope on Cerro Pachon, Chile. Commissioning at the telescope will take place by the end of 2013, matching the summer window of the southern hemisphere. According to current estimates the 3 year planet finding campaign (890 allocated hours) might discover, image, and spectroscopically analyze 20 to 40 new exo-planets. Final acceptance testing of the integrated instrument can always bring up surprises when using cold chamber and flexure rig installations. The latest developments are reported. Also, we will give an overview of GPI's lab performance, the interplay between subsystems such as the calibration unit (CAL) with the AO bench. We report on-going optimizations on the AO controller loop to filter vibrations and last but not least achieved contrast performance applying speckle nulling. Furthermore, we will give an outlook of possible but challenging future upgrades as the implementation of a predictive controller or exchanging the conventional 48x48 SH WFS with a pyramid. With the ELT era arising, GPI will proof as a versatile and path-finding testbed for AO technologies on the next generation of ground-based telescopes.

We present a spectroscopically confirmed sample of Lyman alpha emitting galaxies (LAEs) at z ~ 4.5 in the Extended Chandra Deep Field South (ECDFS), which we combine with a sample of z ~ 4.5 LAEs from the Large Area Lyman Alpha (LALA) survey to build a unified Lya luminosity function (LF). We spectroscopically observed 64 candidate LAEs in the ECDFS, confirming 46 objects as z~4.5 LAEs. We did not detect significant flux from neither C_iv 1549\AA\ nor the He_ii 1640\AA\ emission in individual LAE spectra, even with a coadded spectrum. With the coadded line ratio of He_ii to Lya constraining the Population III star formation rate (SFR) to be <0.3% of the total SFR, and <1.25% of the observed SFR (both at the 2-$\sigma$ level). Only one LAE was detected in both the X-ray and radio, while the other objects remained undetected, even when stacked. The Lya LF in our two deepest narrowband filters in the ECDFS differ at >2$\sigma$ significance, and the product $L^*\Phi^*$ differs by a factor of >3. Similar LF differences have been used to infer evolution in the neutral gas fraction in the intergalactic medium at z>6, yet here the difference is likely due to cosmic variance, given that the two samples are from adjoining line-of-sight volumes. Combining our new sample of LAEs with those from previous LALA narrowband surveys at z = 4.5, we obtain one of the best measured Lya LFs to date of L* = 42.83 $\pm$ 0.06 and $\Phi^*$ = -3.48 $\pm$ 0.09. We compare our new LF to others from the literature to study the evolution of the Lya luminosity density at 0 < z < 7. We find tentative evidence for evolution in the product $L^* \Phi^*$, which approximately tracks the cosmic SFR density, but since field-to-field and survey-to-survey variations are in some cases as large as the possible evolution, some caution is needed in interpreting this trend.

We present spectroscopic measurements of the [OIII] emission line from two subregions of strong Lyman-alpha emission in a radio-quiet Lyman-alpha blob (LAB). The blob under study is LAB1 (Steidel et al. 2000) at z ~ 3.1, and the [OIII] detections are from the two Lyman break galaxies embedded in the blob halo. The [OIII] measurements were made with LUCIFER on the 8.4m Large Binocular Telescope and NIRSPEC on 10m Keck Telescope. Comparing the redshift of the [OIII] measurements to Lyman-alpha redshifts from SAURON (Weijmans et al. 2010) allows us to take a step towards understanding the kinematics of the gas in the blob. Using both LUCIFER and NIRSPEC we find velocity offsets between the [OIII] and Lyman-alpha redshifts that are modestly negative or consistent with 0 km/s in both subregions studied (ranging from -72 +/- 42 -- +6 +/- 33 km/s). A negative offset means Lyman-alpha is blueshifted with respect to [OIII], a positive offset then implies Lyman-alpha is redshifted with respect to [OIII]. These results may imply that outflows are not primarily responsible for Lyman alpha escape in this LAB, since outflows are generally expected to produce a positive velocity offset (McLinden et al. 2011). In addition, we present an [OIII] line flux upper limit on a third region of LAB1, a region that is unassociated with any underlying galaxy. We find that the [OIII] upper limit from the galaxy-unassociated region of the blob is at least 1.4 -- 2.5 times fainter than the [OIII] flux from one of the LBG-associated regions and has an [OIII] to Lyman-alpha ratio measured at least 1.9 -- 3.4 times smaller than the same ratio measured from one of the LBGs.

We present the results of a search for the most luminous star-forming galaxies at redshifts z~6 based on CFHT Legacy Survey data. We identify a sample of 40 Lyman break galaxies brighter than magnitude z'=25.3 across an area of almost 4 square degrees. Sensitive spectroscopic observations of seven galaxies provide redshifts for four, of which only two have moderate to strong Lyman alpha emission lines. All four have clear continuum breaks in their spectra. Approximately half of the Lyman break galaxies are spatially resolved in 0.7 arcsec seeing images, indicating larger sizes than lower luminosity galaxies discovered with the Hubble Space Telescope, possibly due to on-going mergers. The stacked optical and infrared photometry is consistent with a galaxy model with stellar mass ~ 10^{10} solar masses. There is strong evidence for substantial dust reddening with a best-fit A_V=0.7 and A_V>0.48 at 2 sigma confidence, in contrast to the typical dust-free galaxies of lower luminosity at this epoch. The spatial extent and spectral energy distribution suggest that the most luminous z~6 galaxies are undergoing merger-induced starbursts. The luminosity function of z=5.9 star-forming galaxies is derived. This agrees well with previous work and shows strong evidence for an exponential decline at the bright end, indicating that the feedback processes which govern the shape of the bright end are occurring effectively at this epoch.

Lyman-alpha emitting galaxies can be used to study cosmological reionization, because a neutral intergalactic medium scatters Lyman-alpha photons into diffuse halos whose surface brightness falls below typical survey detection limits. Here we present the Lyman-alpha emitting galaxy LAE J095950.99+021219.1, identified at redshift z=6.944 in the COSMOS field using narrowband imaging and followup spectroscopy with the IMACS instrument on the Magellan I Baade telescope. With a single object spectroscopically confirmed so far, our survey remains consistent with a wide range of IGM neutral fraction at redshift seven, but further observations are planned and will help clarify the situation. Meantime, the object we present here is only the third Lyman-alpha selected galaxy to be spectroscopically confirmed at redshift seven, and is 2--3 times fainter than the previously confirmed redshift seven Lyman alpha galaxies.

This paper is a response to a call for white papers solicited by Gemini Observatory and its Science and Technology Advisory Committee, to help define the science case and requirements for a new Gemini instrument, envisaged to consist of a single-object spectrograph at medium resolution simultaneously covering optical and near-infrared wavelengths. In this white paper we discuss the science case for an alternative new instrument, consisting instead of a multi-object, medium-resolution, high-throughput spectrograph, covering simultaneously the optical and near-infrared slices of the electromagnetic spectrum. We argue that combination of wide wavelength coverage at medium resolution with moderate multiplexing power is an innovative path that will enable the pursuit of fundamental science questions in a variety of astrophysical topics, without compromise of the science goals achievable by single-object spectroscopy on a wide baseline. We present a brief qualitative discussion of the main features of a notional hardware design that could conceivably make such an instrument viable.

The study of Ly-alpha emission in the high-redshift universe is a useful probe of the epoch of reionization, as the Ly-alpha line should be attenuated by the intergalactic medium (IGM) at low to moderate neutral hydrogen fractions. Here we present the results of a deep and wide imaging search for Ly-alpha emitters in the COSMOS field. We have used two ultra-narrowband filters (filter width of ~8-9 {\deg}A) on the NEWFIRM camera, installed on the Mayall 4m telescope at Kitt Peak National Observatory, in order to isolate Ly-alpha emitters at z = 7.7; such ultra-narrowband imaging searches have proved to be excellent at detecting Ly-alpha emitters. We found 5-sigma detections of four candidate Ly-alpha emitters in a survey volume of 2.8 x 10^4 Mpc^3 (total survey area ~760 arcmin^2). Each candidate has a line flux greater than 8 x 10^-18 erg s^-1 cm^-2. Using these results to construct a luminosity function and comparing to previously established Ly-alpha luminosity functions at z = 5.7 and z = 6.5, we find no conclusive evidence for evolution of the luminosity function between z = 5.7 and z = 7.7. Statistical Monte Carlo simulations suggest that half of these candidates are real z = 7.7 targets, and spectroscopic follow-up will be required to verify the redshift of these candidates. However, our results are consistent with no strong evolution in the neutral hydrogen fraction of the IGM between z = 5.7 and z = 7.7, even if only one or two of the z = 7.7 candidates are spectroscopically confirmed.

Aims. Ly-alpha emitters (LAEs) can be detected out to very high redshifts during the epoch of reionization. The evolution of the LAE luminosity function with redshift is a direct probe of the Ly-alpha transmission of the intergalactic medium (IGM), and therefore of the IGM neutral-hydrogen fraction. Measuring the Ly-alpha luminosity function (LF) of LAEs at redshift z = 7.7 therefore allows us to constrain the ionizing state of the Universe at this redshift. Methods. We observed three 7.5'x7.5' fields with the HAWK-I instrument at the VLT with a narrow band filter centred at 1.06 $\mu$m and targeting LAEs at redshift z ~ 7.7. The fields were chosen for the availability of multiwavelength data. One field is a galaxy cluster, the Bullet Cluster, which allowed us to use gravitational amplification to probe luminosities that are fainter than in the field. The two other fields are subareas of the GOODS Chandra Deep Field South and CFHTLS-D4 deep field. We selected z=7.7 LAE candidates from a variety of colour criteria, in particular from the absence of detection in the optical bands. Results. We do not find any LAE candidates at z = 7.7 in ~2.4 x 10^4 Mpc^3 down to a narrow band AB magnitude of ~ 26, which allows us to infer robust constraints on the Ly-alpha LAE luminosity function at this redshift. Conclusions. The predicted mean number of objects at z = 6.5, derived from somewhat different LFs of Hu et al. (2010), Ouchi et al. (2010), and Kashikawa et al. (2011) are 2.5, 13.7, and 11.6, respectively. Depending on which of these LFs we refer to, we exclude a scenario with no evolution from z = 6.5 to z = 7.7 at 85% confidence without requiring a strong change in the IGM Ly-alpha transmission, or at 99% confidence with a significant quenching of the IGM Ly-alpha transmission, possibly from a strong increase in the high neutral-hydrogen fraction between these two redshifts.

We report a search for z=7 Ly-alpha emitters (LAEs) using a custom-made Narrow-Band filter, centered at 9755 Angstroms, with the instrument Suprime-Cam installed at the Subaru telescope. We observed two different fields and obtained two sample of 7 Ly-alpha emitters of which 4 are robust in each field. We are covering the luminosity range of 9.10^{42} - 2.10^{43} erg/s in comoving volumes of ~ 4.10^{5} and 4.3.10^{5} Mpc^{3}. From this result, we derived possible z~7 Ly-alpha luminosity functions for the full samples and for a subsample of 4 objects in each field. We do not observe, in each case, any strong evolution between the z=6.5 and z~7 Ly-alpha luminosity functions. Spectroscopic confirmation for these candidate samples is required to establish a definitive measure of the luminosity function at z~7.

We report a search for z~6.96 Ly-alpha emitters (LAEs) using a Narrow-Band filter, centered at 9680 Angstroms, with the IMACS instrument on the Magellan telescope at Las Campanas Observatory. We obtain a sample of 6 Ly-alpha emitter candidates of luminosity ~10^42 erg/s in a total area of 465 square arcmin corresponding to a comoving volume of ~ 72000 Mpc^3. From this result, we derive a Ly-alpha luminosity function (LF) at z~6.96 and compare our sample with the only z~6.96 Ly-alpha emitter spectroscopically confirmed to date (Iye et al. 2006). We find no evolution between the z=5.7 and z~7 Ly-alpha luminosity functions, if a majority of our candidates are confirmed. Spectroscopic confirmation for this sample will enable more robust conclusions.

We present the first spectroscopic measurements of the [OIII] 5007 A line in two z ~ 3.1 Lyman-alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFER1 on the 8.4m Large Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used to identify these Lya emitting galaxies. Using the [OIII] line we have measured accurate systemic redshifts for these two galaxies, and discovered a velocity offset between the [OIII] and Ly-alpha lines in both, with the Lya line peaking 342 and 125 km/s redward of the systemic velocity. These velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the transmission of Lya photons through the interstellar medium and intergalactic medium around the galaxies. By measuring these offsets directly, we can refine both Lya-based tests for reionization, and Lya luminosity function measurements where the Lya forest affects the blue wing of the line. Our work also provides the first direct constraints on the strength of the [OIII] line in high-redshift LAEs. We find [OIII] fluxes of 7 and 36 x 10^-17 erg s^-1 cm^-2 in two z ~ 3.1 LAEs. These lines are strong enough to dominate broad-band flux measurements that include the line (in thiscase, K_s band photometry). Spectral energy distribution fits that do not account for the lines would therefore overestimate the 4000 A (and/or Balmer) break strength in such galaxies, and hence also the ages and stellar masses of such high-z galaxies.

Lyman alpha (Lya) emission lines should be attenuated in a neutral intergalactic medium (IGM). Therefore the visibility of Lya emitters at high redshifts can serve as a valuable probe of reionization at about the 50% level. We present an imaging search for z=7.7 Lya emitting galaxies using an ultra-narrowband filter (filter width= 9A) on the NEWFIRM imager at the Kitt Peak National Observatory. We found four candidate Lya emitters in a survey volume of 1.4 x 10^4 Mpc^3, with a line flux brighter than 6x10^-18 erg/cm^2/s (5 sigma in 2" aperture). We also performed a detailed Monte-Carlo simulation incorporating the instrumental effects to estimate the expected number of Lya emitters in our survey, and found that we should expect to detect one Lya emitter, assuming a non-evolving Lya luminosity function (LF) between z=6.5 and z=7.7. Even if one of the present candidates is spectroscopically confirmed as a z~8 Lya emitter, it would indicate that there is no significant evolution of the Lya LF from z=3.1 to z~8. While firm conclusions would need both spectroscopic confirmations and larger surveys to boost the number counts of galaxies, we successfully demonstrate the feasibility of sensitive near-infrared (1.06 um) narrow-band searches using custom filters designed to avoid the OH emission lines that make up most of the sky background.

We present the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog, which is based upon the SDSS Seventh Data Release. The catalog, which contains 105,783 spectroscopically confirmed quasars, represents the conclusion of the SDSS-I and SDSS-II quasar survey. The catalog consists of the SDSS objects that have luminosities larger than M_i = -22.0 (in a cosmology with H_0 = 70 km/s/Mpc Omega_M = 0.3, and Omega_Lambda = 0.7) have at least one emission line with FWHM larger than 1000 km/s or have interesting/complex absorption features, are fainter than i > 15.0 and have highly reliable redshifts. The catalog covers an area of 9380 deg^2. The quasar redshifts range from 0.065 to 5.46, with a median value of 1.49; the catalog includes 1248 quasars at redshifts greater than four, of which 56 are at redshifts greater than five. The catalog contains 9210 quasars with i < 18; slightly over half of the entries have i< 19. For each object the catalog presents positions accurate to better than 0.1" rms per coordinate, five-band (ugriz) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains radio, near-infrared, and X-ray emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3800-9200 Ang. at a spectral resolution R = 2000 the spectra can be retrieved from the SDSS public database using the information provided in the catalog. Over 96% of the objects in the catalog were discovered by the SDSS. We also include a supplemental list of an additional 207 quasars with SDSS spectra whose archive photometric information is incomplete.

The Ly-alpha luminosity function (LF) of high-redshift Ly-alpha emitters (LAEs) is one of the few observables of the re-ionization epoch accessible to date with 8-10 m class telescopes. The evolution with redshift allows one to constrain the evolution of LAEs and their role in re-ionizing the Universe at the end of the Dark Ages. We have performed a narrow-band imaging program at 1.06 microns at the CFHT, targeting Ly-alpha emitters at redshift z ~ 7.7 in the CFHT-LS D1 field. From these observations we have derived a photometric sample of 7 LAE candidates at z ~ 7.7. We derive luminosity functions for the full sample of seven objects and for sub-samples of four objects. If the brightest objects in our sample are real, we infer a luminosity function which would be difficult to reconcile with previous work at lower redshift. More definitive conclusions will require spectroscopic confirmation.

We present the discovery of a candidate of giant radio-quiet Lyman-alpha (Lya) blob (RQLAB) in a large-scale structure around a high-redshift radio galaxy (HzRG) lying in a giant Lya halo, B3 J2330+3927 at redshift z=3.087. We obtained Lya imaging around B3 J2330+3927 with Subaru/Suprime-Cam to search for Lya emitters (LAEs) and absorbers (LAAs) at redshift z=3.09+-0.03. We detected candidate 127 LAEs and 26 LAAs in the field of view of 31' x 24'. We found that B3 J2330+3927 is surrounded by a 130 kpc Lya halo and a large-scale (60 x 20 comoving Mpc) filamentary structure. The large-scale structure contains one prominent local density peak with an overdensity of greater than 5, which is 8' (15 comoving Mpc) away from B3 J2330+3927. In this peak, we discovered a candidate 100 kpc RQLAB. The existence of both types of Lya nebulae in the same large-scale structure suggests that giant Lya nebulae need special large-scale environments to form. On smaller scales, however, the location of B3 J2330+3927 is not a significant local density peak in this structure, in contrast to the RQLAB. There are two possible interpretations of the difference of the local environments of these two Lya nebulae. Firstly, RQLAB may need a prominent (delta ~ 5) density peak of galaxies to form through intense star-bursts due to frequent galaxy interactions/mergers and/or continuous gas accretion in an overdense environment. On the other hand, Lya halo around HzRG may not always need a prominent density peak to form if the surrounding Lya halo is mainly powered by its radio and AGN activities. Alternatively, both RQLAB and Lya halo around HzRG may need prominent density peaks to form but we could not completely trace the density of galaxies because we missed evolved and dusty galaxies in this survey.

Aims: The first star forming galaxies in the early universe should be copious Ly alpha emitters, and may play a significant role in ionizing the intergalactic medium (IGM). It has been proposed that the luminosity function of Lya emitting galaxies beyond z~6 may be used to constrain the neutral fraction of the IGM during this epoch. In this work we report on a search for Ly alpha emitters at redshift 8.8. Methods: We performed a narrow band imaging programme using ISAAC at the ESO VLT. Seven fields, covering a total area of 31sq. arcmin and for which optical and broad band infra-red images have been obtained in the GOODS survey, were imaged to a limiting flux (respectively luminosity) of ~ 1.3 x 10^{-17} ergs.s^{-1}.cm^{-2} (respectively ~ 1.3 x 10^{43} ergs.s^{-1} in a narrow band filter centered in a region of low OH sky emission at 1.19 micron. Candidate Lyman alpha emitters are objects that are detected in the ISAAC NB images and undetected in the visible broad band images. Results: No z=8.8 Ly alpha emitting galaxies were detected to a limit approaching recent estimates of the luminosity function at z ~ 6. Our results do suggest, however, that detections or substantial constraints could be achieved by this method in the near future with larger field instruments planned for various telescopes.