Search

We advocate for the prioritization of high-precision photometric monitoring of exoplanet and brown dwarf companions to detect brightness variability arising from features in their atmospheres. Measurements of photometric variability provide not only an insight into the physical appearances of these companions, but are also a direct probe of their atmospheric structures and dynamics, and yield valuable estimates of their rotation periods. JWST is uniquely capable of monitoring faint exoplanet companions over their full rotation periods, thanks to its inherent stability and powerful high-contrast coronagraphic imaging modes. Rotation period measurements can be further combined with measurements of v sin i obtained using high-resolution spectroscopy to infer the viewing angle of a companion. Photometric monitoring over multiple rotation periods and at multiple epochs will allow both short- and long-term time evolution in variability signals to be traced. Furthermore, the differences between the layers in a companion's atmosphere can be probed by obtaining simultaneous photometric monitoring at different wavelengths through NIRCam dual-band coronagraphy. Overall, JWST will reach the highest sensitivities to variability to date and enable the light curves of substellar companions to be characterised with unprecedented cadence and precision at the sub-percent level., Comment: 4 pages, 2 figures, white paper submitted in response to the call by the Working Group on Strategic Exoplanet Initiatives with HST and JWST (details at https://outerspace.stsci.edu/display/HPR/Strategic+Exoplanet+Initiatives+with+HST+and+JWST & final report at arXiv:2404.02932), adapted to include author list and affiliations

Spatially-resolved images of debris disks are necessary to determine disk morphological properties and the scattering phase function (SPF) which quantifies the brightness of scattered light as a function of phase angle. Current high-contrast imaging instruments have successfully resolved several dozens of debris disks around other stars, but few studies have investigated trends in the scattered-light, resolved population of debris disks in a uniform and consistent manner. We have combined Karhunen-Loeve Image Projection (KLIP) with radiative-transfer disk forward modeling in order to obtain the highest quality image reductions and constrain disk morphological properties of eight debris disks imaged by the Gemini Planet Imager at H-band with a consistent and uniformly-applied approach. In describing the scattering properties of our models, we assume a common SPF informed from solar system dust scattering measurements and apply it to all systems. We identify a diverse range of dust density properties among the sample, including critical radius, radial width, and vertical width. We also identify radially narrow and vertically extended disks that may have resulted from substellar companion perturbations, along with a tentative positive trend in disk eccentricity with relative disk width. We also find that using a common SPF can achieve reasonable model fits for disks that are axisymmetric and asymmetric when fitting models to each side of the disk independently, suggesting that scattering behavior from debris disks may be similar to Solar System dust., Comment: 23+5 pages, 12+6 figures, 15 pages of Online Supplemental Material included; Accepted for publication in MNRAS

The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models., Comment: 32 pages, 16 figures, 6 tables, 2 appendices

The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet-disk interactions. The multi-wavelength surface brightness, the disk color and geometry permit identification of any asymmetries such as warps or disk offsets from the central star. We find that nineteen of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system are discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power-law, $q$. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications on the evolution of debris disk systems around stars of various spectral types., Comment: 46 pages, 20 figures, 6 tables, accepted for publication in ApJ

We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{\mu m}$ as a part of the JWST Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of $0.5\lambda/D$ for an interferometer), which are inaccessible with the classical inner working angles of the JWST coronagraphs. When combined with JWST's unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a $5\sigma$ contrast of $\Delta m{\sim}7.62{\pm}0.13$ mag relative to the host star at separations ${\gtrsim}0.07{"}$, and the contrast deteriorates steeply at separations ${\lesssim}0.07{"}$. However, we detect no additional companions interior to the known companion HIP 65426 b (at separation ${\sim}0.82{"}$ or, $87^{+108}_{-31}\,\rm{au}$). Our observations thus rule out companions more massive than $10{-}12\,\rm{M_{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP 65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on JWST is sensitive to planetary mass companions at close-in separations (${\gtrsim}0.07{"}$), even for thousands of more distant stars at $\sim$100 pc, in addition to the stars in the nearby young moving groups as stated in previous works. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening an essentially unexplored parameter space., Comment: 20 pages, 12 figures, submitted to ApJL

We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim \lambda/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy., Comment: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measurements using precision radial velocity monitoring calibrated with a suite of state-of-the-art absolute and relative wavelength references. MODHIS is the counterpart to HISPEC for the Thirty Meter Telescope and is being developed in parallel with similar scientific goals. In this proceeding, we provide a brief overview of the current design of both instruments, and the requirements for the two spectrographs as guided by the scientific goals for each. We then outline the current science case for HISPEC and MODHIS, with focuses on the science enabled for exoplanet discovery and characterization. We also provide updated sensitivity curves for both instruments, in terms of both signal-to-noise ratio and predicted radial velocity precision., Comment: 25 pages, 9 figures. To appear in the Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

We report direct observational evidence for a latitudinal dependence of dust cloud opacity in ultracool dwarfs, indicating that equatorial latitudes are cloudier than polar latitudes. These results are based on a strong positive correlation between the viewing geometry and the mid-infrared silicate absorption strength in mid-L dwarfs using mid-infrared spectra from the Spitzer Space Telescope and spin axis inclination measurements from available information in the literature. We confirmed that the infrared color anomalies of L dwarfs positively correlate with dust cloud opacity and viewing geometry, where redder objects are inclined equator-on and exhibit more opaque dust clouds while dwarfs viewed at higher latitudes and with more transparent clouds are bluer. These results show the relevance of viewing geometry to explain the appearance of brown dwarfs and provide insight into the spectral diversity observed in substellar and planetary atmospheres. We also find a hint that dust clouds at similar latitudes may have higher opacity in low-surface gravity dwarfs than in higher-gravity objects., Comment: Accepted for publication in ApJL

Analysis of all archival 5--14 micron spectra of field ultracool dwarfs from the Infrared Spectrograph on the Spitzer Space Telescope has shown that absorption by silicates in the 8--11 micron region is seen in most L-type (1300 K to 2200 K) dwarfs. The absorption is caused by silicate-rich clouds in the atmospheres of L dwarfs and is strongest at L4--L6 spectral types. Herein we compare averages of the mid-infrared silicate absorption signatures of L3--L7 dwarfs that have low ($\lesssim$10$^{4.5}$ cm s$^{-2}$) vs.\ high ($\gtrsim$10$^5$ cm s$^{-2}$) surface gravity. We find that the silicate absorption feature is sensitive to surface gravity and indicates a difference in grain size and composition between dust condensates in young and old mid-L dwarfs. The mean silicate absorption profile of low-gravity mid-L dwarfs matches expectations for $\sim$1 micron-sized amorphous iron- and magnesium-bearing pyroxene (Mg$_x$Fe$_{1-x}$SiO$_3$) grains. High-gravity mid-L dwarfs have silicate absorption better represented by smaller ($\lesssim$0.1 $\mu$m) and more volatile amorphous enstatite (MgSiO$_3$) or SiO grains. This is the first direct spectroscopic evidence for gravity-dependent sedimentation of dust condensates in ultracool atmospheres. It confirms theoretical expectations for lower sedimentation efficiencies in low-gravity atmospheres and independently confirms their increased dustiness., Comment: Accepted for publication in MNRAS

We use 1-m class telescopes and the Transiting Exoplanet Survey Satellite (TESS) to explore the photometric variability of all known rapidly rotating ($v\sin{i}\gtrsim30$ km\,s$^{-1}$) ultra-cool ($\geq$M7) dwarfs brighter than $I\approx17.5$ mag. For a sample of 13 M7--L1.5 dwarfs without prior photometric periods, we obtained $I$-band light curves with the SMARTS 1.3m and WIYN 0.9m telescopes and detected rotation-modulated photometric variability in three of them. Seven of our targets were also observed by TESS and six of them show significant periodicities compatible with the estimated rotation periods of the targets. We investigate the potential of TESS to search for rotation-modulated photometric variability in ultra-cool dwarfs and find that its long stare enables $<$80~h periodic variations to be retrieved with $\leq$1\% amplitudes for ultra-cool dwarfs up to a TESS magnitude of 16.5. We combine these results with the periods of all other known photometrically-periodic ultra-cool dwarfs from the literature, and find that the periods of ultra-cool dwarfs range between 1 and 24 h, although the upper limit is likely an observational bias. We also observe that the minimum rotation periods follow a lower envelope that runs from $\approx$2 h at spectral type $\approx$M8 to $\approx$1 h at spectral type T., Comment: Accepted for publication in MNRAS

We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $\mu$m to 20 $\mu$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion., Comment: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $\mu$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $\mu$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $\mu$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$\sigma$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $\mu$m and MIRI at 11.3 $\mu$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $\mu$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail., Comment: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

We present Gemini South/IGRINS observations of the 1060 K T6 dwarf 2MASS J08173001$-$6155158 with unprecedented resolution ($R\equiv\lambda/\Delta\lambda=45\,000$) and signal-to-noise ratio (SNR > 200) for a late-type T dwarf. We use this benchmark observation to test the reliability of molecular line lists used up-to-date atmospheric models. We determine which spectroscopic regions should be used to estimate the parameters of cold brown dwarfs and, by extension, exoplanets. We present a detailed spectroscopic atlas with molecular identifications across the $H$ and $K$ bands of the near-infrared. We find that water (H$_2$O) line lists are overall reliable. We find the most discrepancies amongst older methane (CH$_4$) line lists, and that the most up-to-date CH$_4$ line lists correct many of these issues. We identify individual ammonia (NH$_3$) lines, a hydrogen sulfide (H$_2$S) feature at 1.5900 $\mu$m, and a molecular hydrogen (H$_2$) feature at 2.1218 $\mu$m. These are the first unambiguous detections of H$_2$S and H$_2$ absorption features in an extra-solar atmosphere. With the H$_2$ detection, we place an upper limit on the atmospheric dust concentration of this T6 dwarf: at least 500 times less than the interstellar value, implying that the atmosphere is effectively dust-free. We additionally identify several features that do not appear in the model spectra. Our assessment of the line lists is valuable for atmospheric model applications to high-dispersion, low-SNR, high-background spectra, such as an exoplanet around a star. We demonstrate a significant enhancement in the detection of the CH$_4$ absorption signal in this T6 dwarf with the most up-to-date line lists., Comment: 45 pages, 63 figures, 5 tables, Accepted to MNRAS

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

We present a uniform analysis of all mid-infrared $R\approx90$ spectra of field M5--T9 dwarfs obtained with the Spitzer Infrared Spectrograph (IRS). The sample contains 113 spectra out of which 12 belong to late-M dwarfs, 69 to L dwarfs, and 32 to T dwarfs. Sixty-eight of these spectra are presented for the first time. We measure strengths of the main absorption bands in the IRS spectra, namely H$_2$O at 6.25 $\mu$m, CH$_4$ at 7.65 $\mu$m, NH$_3$ at 10.5 $\mu$m, and silicates over 8--11 $\mu$m. Water absorption is present in all spectra and strengthens with spectral type. The onset of methane and ammonia occurs at the L8 and T2.5 types, respectively, although ammonia can be detectable as early as T1.5. Silicate absorption sets in at spectral type L2, is on average the strongest in L4--L6 dwarfs, and disappears past L8. However, silicate absorption can also be absent from the spectra at any L subtype. We find a positive correlation between the silicate absorption strength and the excess (deviation from median) near-infrared colour at a given L subtype, which supports the idea that variations of silicate cloud thickness produce the observed colour scatter in L dwarfs. We also find that variable L3--L7 dwarfs are twice more likely to have above-average silicate absorption than non-variables. The ensemble of results solidifies the evidence for silicate condensate clouds in the atmospheres of L dwarfs, and for the first time observationally establishes their emergence and sedimentation between effective temperatures of $\approx$2000 K and $\approx$1300 K, respectively., Comment: Accepted for publication in MNRAS

We present a survey for photometric variability in young, low-mass brown dwarfs with the Spitzer Space Telescope. The 23 objects in our sample show robust signatures of youth and share properties with directly-imaged exoplanets. We present three new young objects: 2MASS J03492367$+$0635078, 2MASS J09512690 $-$8023553 and 2MASS J07180871$-$6415310. We detect variability in 13 young objects, and find that young brown dwarfs are highly likely to display variability across the L2--T4 spectral type range. In contrast, the field dwarf variability occurrence rate drops for spectral types $>$L9. We examine the variability amplitudes of young objects and find an enhancement in maximum amplitudes compared to field dwarfs. We speculate that the observed range of amplitudes within a spectral type may be influenced by secondary effects such as viewing inclination and/or rotation period. We combine our new rotation periods with the literature to investigate the effects of mass on angular momentum evolution. While high mass brown dwarfs ($>30 M_{\mathrm{Jup}}$) spin up over time, the same trend is not apparent for lower mass objects ($<30 M_{\mathrm{Jup}}$), likely due to the small number of measured periods for old, low-mass objects. The rotation periods of companion brown dwarfs and planetary-mass objects are consistent with those of isolated objects with similar ages and masses, suggesting similar angular momentum histories. Within the AB Doradus group, we find a high variability occurrence rate and evidence for common angular momentum evolution. The results are encouraging for future variability searches in directly-imaged exoplanets with facilities such as the James Webb Space Telescope and 30-meter telescopes.

We present $Spitzer$ IRS 5--14 $\mu$m spectra and 16 $\mu$m and 22 $\mu$m photometry of the T2.5 companion to the $\sim$300 Myr-old G0V star HN Peg. We incorporate previous 0.8--5 $\mu$m observations to obtain the most comprehensive spectral energy distribution of an intermediate-gravity L/T-transition dwarf which, together with an accurate Gaia EDR3 parallax of the primary, enable us to derive precise fundamental parameters. We find that young ($\approx$0.1--0.3 Gyr) early-T dwarfs on average have $\approx$140 K lower effective temperatures, $\approx$20% larger radii, and similar bolometric luminosities compared to $\gtrsim$1 Gyr-old field dwarfs with similar spectral types. Our accurate infrared spectrophotometry offers new detail at wavelengths where the dominant carbon-bearing molecules have their strongest transitions: at 3.4 $\mu$m for methane and at 4.6 $\mu$m for carbon monoxide. We assess the performance of various widely available photospheric models and find that models with condensates and/or clouds better reproduce the full SED of this moderately young early-T dwarf. However, cloud-free models incorporating a more general convective instability treatment reproduce at least the low-resolution near-IR spectrum similarly well. Our analysis of $R\approx2300$ $J$-band spectra shows that the near-infrared potassium absorption lines in HN Peg B have similar strengths to those seen in both younger and older T2-T3 dwarfs. We conclude that while alkali lines are well-established as surface gravity indicators for L-type or warmer stars, they are insensitive to surface gravity in early-T dwarfs, Comment: Accepted for publication in ApJ

Using a novel wide-slit, multi-object approach with the GMOS spectrograph on the 8-meter Gemini South telescope, we have obtained precise time-series spectrophotometry of the binary brown dwarf Luhman 16 at optical wavelengths over two full nights. The B component of this binary system is known to be variable in the red optical and near-infrared with a period of 5 hr and an amplitude of 5--20%. Our observations probe its spectrally-resolved variability in the 6000--10000 Angstrom range. At wavelengths affected by the extremely strong, broadened spectral lines of the neutral alkali metals (the potassium doublet centered near 7682 Angstroms and the sodium doublet at 5893 Angstroms), we see photometric variations that differ strikingly from the those of the 8000--10000 Angstrom `red continuum' that dominates our detected flux. On UT 2014 February 24, these variations are anticorrelated with the red continuum, while on Feb 25 they have a large relative phase shift. The extent to which the wavelength-dependent photometric behavior diverges from that of the red continuum appears to correlate with the strength of the alkali absorption. We consider but ultimately reject models in which our observations are explained by lightning or auroral activity. A more likely cause is cloud-correlated, altitude-dependent variations in the gas-phase abundances of sodium and potassium, which are in chemical equilibrium with their chlorides in brown dwarf atmospheres. Clouds could influence these chemical equilibria by changing the atmospheric temperature profile and/or through cloud particles acting as chemical catalysts., Comment: 26 pages, 14 figures, accepted by ApJ

We present the discovery of rapid photometric variability in three ultra-cool dwarfs from long-duration monitoring with the Spitzer Space Telescope. The T7, L3.5, and L8 dwarfs have the shortest photometric periods known to date: ${1.080}^{+0.004}_{-0.005}$ h, ${1.14}^{+0.03}_{-0.01}$ h, and ${1.23}^{+0.01}_{-0.01}$ h, respectively. We confirm the rapid rotation through moderate-resolution infrared spectroscopy, which reveals projected rotational velocities between 79 and 104 km s$^{-1}$. We compare the near-infrared spectra to photospheric models to determine the objects' fundamental parameters and radial velocities. We find that the equatorial rotational velocities for all three objects are $\gtrsim$100 km s$^{-1}$. The three L and T dwarfs reported here are the most rapidly spinning and likely the most oblate field ultra-cool dwarfs known to date. Correspondingly, all three are excellent candidates for seeking auroral radio emission and net optical/infrared polarization. As of this writing, 78 L-, T-, and Y-dwarf rotation periods have now been measured. The clustering of the shortest rotation periods near 1 h suggests that brown dwarfs are unlikely to spin much faster., Comment: 35 pages, 13 figures, 6 tables, accepted for publication in the Astronomical Journal, typos corrected in updated version

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., Comment: Accepted for publication in AJ. 3 tables, 11 figures

We present \textit{Spitzer Space Telescope} variability monitoring observations of three low-gravity L dwarfs with previous detections of variability in the near-IR, 2MASS J0045+16, 2MASS J0501-00 and 2MASS J1425-36. We detect significant, periodic variability in two of our targets, 2MASS J0045+16 and 2MASS J0501-00. We do not detect variability in 2MASS J1425-36. Combining our new rotation periods with rotational velocities, we calculate inclination angles of $22\pm1^{\circ}$, ${60^{+13 }_{-8}} ^{\circ}$ and $52^{+19}_{-13}~^{\circ}$ for 2MASS J0045+16, 2MASS J0501-00 and 2MASS J1425-36 respectively. Our three new objects are consistent with the tentative relations between inclination, amplitude and color anomaly previously reported. Objects with the highest variability amplitudes are inclined equator-on, while the maximum observed amplitude decreases as the inclination angle decreases. We also find a correlation between the inclination angle and $(J-K)_{\mathrm{2MASS}}$ color anomaly for the sample of objects with measured inclinations. Compiling the entire sample of brown dwarfs with \textit{Spitzer} variability detections, we find no enhancement in amplitude for young, early-L dwarfs compared to the field dwarf population. We find a possible enhancement in amplitude of low-gravity late-L dwarfs at $4.5~\mu$m. We do not find a correlation between amplitude ratio and spectral type for field dwarfs or for the young population. Finally, we compile the rotation periods of a large sample of brown dwarfs with ages 1 Myr to 1 Gyr and compare the rotation rates predicted by evolutionary models assuming angular momentum conservation. We find that the rotation rates of the current sample of brown dwarfs fall within the expected range set by evolutionary models and breakup limits., Comment: Accepted for publication in AJ

We present new parallax measurements from the CFHT Infrared Parallax Program and the Pan-STARRS 3$\pi$ Steradian Survey for the young ($\approx150-300$ Myr) triple system VHS J125601.92$-$125723.9. This system is composed of a nearly equal-flux binary ("AB") and a wide, possibly planetary-mass companion ("b"). The system's published parallactic distance ($12.7\pm1.0$ pc) implies absolute magnitudes unusually faint compared to known young objects and is in tension with the spectrophotometric distance for the central binary ($17.2\pm2.6$ pc). Our CFHT and Pan-STARRS parallaxes are consistent, and the more precise CFHT result places VHS J1256-1257 at $22.2^{+1.1}_{-1.2}$ pc. Our new distance results in higher values for the companion's mass ($19\pm5$ M$_{\rm Jup}$) and temperature ($1240\pm50$ K), and also brings the absolute magnitudes of all three components into better agreement with known young objects., Comment: Research Notes of the American Astronomical Society, in press. Arxiv version has additional information about the astrometry (text & figure)

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., Comment: 25 pages, 12 figures

HD106906b is an ~11$M_{\mathrm{Jup}}$, ~15Myr old directly-imaged exoplanet orbiting at an extremely large distance from its host star. The wide separation (7.11 arcsec) between HD106906b and its host star greatly reduces the difficulty in direct-imaging observations, making it one of the most favorable directly-imaged exoplanets for detailed characterization. In this paper, we present HST/WFC3/IR time-resolved observations of HD106906b in the F127M, F139M, and F153M bands. We have achieved ~1% precision in the lightcurves in all three bands. The F127M lightcurve demonstrates marginally-detectable ($2.7\sigma$ significance) variability with a best-fitting period of 4 hr, while the lightcurves in the other two bands are consistent with flat lines. We construct primary-subtracted deep images and use these images to exclude additional companions to HD106906 that are more massive than 4$M_{\mathrm{Jup}}$ and locate at projected distances of more than ~500 au. We measure the astrometry of HD106906b in two HST/WFC3 epochs and achieve precisions better than 2.5 mas. The position angle and separation measurements do not deviate from those in the 2004 HST/ACS/HRC images for more than $1\sigma$ uncertainty. We provide the HST/WFC3 astrometric results for 25 background stars that can be used as reference sources in future precision astrometry studies. Our observations also provide the first 1.4-micron water band photometric measurement for HD106906b. HD106906b's spectral energy distribution and the best-fitting BT-Settl model have an inconsistency in the 1.4-micron water absorption band, which highlights the challenges in modeling atmospheres of young planetary-mass objects., Comment: Accepted for publication to AJ

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., Comment: 29 Pages, 23 Figures. Accepted to AJ

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., Comment: 20 pages, 9 figures

This white paper proposes leveraging high-quality Gaia data available to the worldwide scientific community and complement it with support from Canadian-related facilities to place Canada as a leader in the fields of stellar associations and exoplanet science, and to train Canadian highly qualified personnel through graduate and post-graduate research grants. Gaia has sparked a new era in the study of stellar kinematics by measuring precise distances and proper motions for 1.3 billion stars. These data have already generated more than 1700 scientific papers and are guaranteed to remain the source of many more papers for the upcoming decades. More than 900 new age-calibrated young low-mass stars have already been discovered as a direct consequence of the second Gaia data release. Some of these may already be host stars to known exoplanet systems or may become so with the progress of the TESS mission that is expected to discover 10,000 nearby transiting exoplanets in the upcoming decade. This places Canada in a strategic position to leverage Gaia data because it has access to several high-resolution spectrometers on 1-4 m class telescopes (e.g. The ESPaDOnS, SPIRou and NIRPS), that would allow to quickly characterize this large number of low-mass stars and their exoplanet systems. This white paper describes the opportunity in such scientific projects that could place Canada as a leader in the fields of stellar associations and exoplanets., Comment: Submitted to CASCA Long-Range Plan 2020

Digital tracking enables telescopes to detect asteroids several times fainter than conventional techniques. We describe our optimized methodology to acquire, process, and interpret digital tracking observations, and we apply it to probe the apparent magnitude distribution of main belt asteroids fainter than any previously detected from the ground. All-night integrations with the Dark Energy Camera (DECam) yield 95\% completeness at $R$ magnitude 25.0, and useful sensitivity to $R=25.6$ mag when we use an analytical detection model to correct flux overestimation bias. In a single DECam field observed over two nights, we detect a total of 3234 distinct asteroids, of which 3123 are confirmed on both nights. At opposition from the Sun, we find a sky density of $697 \pm 15$ asteroids per square degree brighter than $R = 25.0$ mag, and $1031 \pm 23$ brighter than $R = 25.6$ mag. We agree with published results for the sky density and apparent magnitude distribution of asteroids brighter than $R=23$ mag. For a power law defined by $dN/dR \propto 10^{\alpha R}$, we find marginally acceptable fits with a constant slope $\alpha = 0.28 \pm 0.02$ from $R=20$ to 25.6 mag. Better fits are obtained for a broken power law with $\alpha=0.218 \pm 0.026$ for $R=20$ to 23.5 mag, steepening to $\alpha=0.340 \pm 0.025$ for $R = 23.5$ to 25.6 mag. The constant or steepening power law indicates asteroids fainter than $R = 23.5$ mag are abundant, contrary to some previous claims but consistent with theory., Comment: 40 pages, 15 figures, accepted by AJ

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., Comment: 14 pages, 9 figures. Accepted for publication in the Astronomical Journal

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}$)., Comment: 17 pages, 11 figures. Accepted for publication in the Astronomical Journal

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., Comment: 36 pages, 17 figures, accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems. Updated version includes revisions made during the referee process

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., Comment: 28 Pages, 13 Figures, Accepted to JATIS

We use the Wide Field Camera 3 on the {\sl Hubble Space Telescope} to spectrophotometrically monitor the young L7.5 companion HD~203030B. Our time series reveal photometric variability at 1.27\,$\mu$m and 1.39\,$\mu$m on time scales compatible with rotation. We find a rotation period of $7.5^{+0.6}_{-0.5}$ h: comparable to those observed in other brown dwarfs and planetary-mass companions younger than 300 Myr. We measure variability amplitudes of $1.1\pm0.3\%$ (1.27\,$\mu$m) and $1.7\pm0.4\%$ (1.39\,$\mu$m), and a phase lag of 56$^\circ\pm$28$^\circ$ between the two light curves. We attribute the difference in photometric amplitudes and phases to a patchy cloud layer that is sinking below the level where water vapor becomes opaque. HD 203030B and the few other known variable young late-L dwarfs are unlike warmer (earlier-type and/or older) L dwarfs, for which variability is much less wavelength-dependent across the 1.1--1.7$\mu$m region. We further suggest that a sinking of the top-most cloud deck below the level where water or carbon monoxide gas become opaque may also explain the often enhanced variability amplitudes of even earlier-type low-gravity L dwarfs. Because these condensate and gas opacity levels are already well-differentiated in T dwarfs, we do not expect the same variability amplitude enhancement in young vs.\ old T dwarfs., Comment: Accepted for publication in The Astrophysical Journal

We present individual dynamical masses for the nearby M9.5+T5.5 binary WISE J072003.20$-$084651.2AB, a.k.a. Scholz's star. Combining high-precision CFHT/WIRCam photocenter astrometry and Keck adaptive optics resolved imaging, we measure the first high-quality parallactic distance ($6.80_{-0.06}^{+0.05}$ pc) and orbit ($8.06_{-0.25}^{+0.24}$ yr period) for this system composed of a low-mass star and brown dwarf. We find a moderately eccentric orbit ($e = 0.240_{-0.010}^{+0.009}$), incompatible with previous work based on less data, and dynamical masses of $99\pm6$ $M_{\rm Jup}$ and $66\pm4$ $M_{\rm Jup}$ for the two components. The primary mass is marginally inconsistent (2.1$\sigma$) with the empirical mass$-$magnitude$-$metallicity relation and models of main-sequence stars. The relatively high mass of the cold ($T_{\rm eff} = 1250\pm40$ K) brown dwarf companion indicates an age older than a few Gyr, in accord with age estimates for the primary star, and is consistent with our recent estimate of $\approx$70 $M_{\rm Jup}$ for the stellar/substellar boundary among the field population. Our improved parallax and proper motion, as well as an orbit-corrected system velocity, improve the accuracy of the system's close encounter with the solar system by an order of magnitude. WISE J0720$-$0846AB passed within $68.7\pm2.0$ kAU of the Sun $80.5\pm0.7$ kyr ago, passing through the outer Oort cloud where comets can have stable orbits., Comment: accepted to AJ

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., Comment: 29 pages, 18 figures, ApJ accepted

Context. The interaction between low-mass companions and the debris discs they reside in is still not fully understood. A debris disc can evolve due to self-stirring, a process in which planetesimals can excite their neighbours to the point of destructive collisions. In addition, the presence of a companion could further stir the disc (companion-stirring). Additional information is necessary to understand this fundamental step in the formation and evolution of a planetary system, and at the moment of writing only a handful of systems are known where a companion and a debris disc have both been detected and studied at the same time. Aims. Our primary goal is to augment the sample of these systems and to understand the relative importance between self-stirring and companion-stirring. Methods. In the course of the VLT/NaCo Imaging Survey for Planets around Young stars (ISPY), we observed HD 193571, an A0 debris disc hosting star at a distance of 68 pc with an age between 60 and 170 Myr. We obtained two sets of observations in L' band and a third epoch in H band using the GPI instrument at Gemini-South. Results. A companion was detected in all three epochs at a projected separation of 11 au (0.17 arcsec), and co-motion was confirmed through proper motion analysis. Given the inferred disc size of 120 au, the companion appears to reside within the gap between the host star and the disc. Comparison between the L' and H band magnitude and evolutionary tracks suggests a mass of 0.31 - 0.39 solar masses . Conclusions. We discovered a previously unknown M-dwarf companion around HD 193571, making it the third low-mass stellar object discovered within a debris disc. A comparison to self- and companion-stirring models suggests that the companion is likely responsible for the stirring of the disc., Comment: 10 pages, 8 figures, accepted for publication in Astronomy&Astrophysics

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

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., Comment: 52 pages, 18 figures. AJ in press

High resolution spectroscopy of the lowest-mass stars and brown dwarfs reveals their origins, multiplicity, compositions and physical properties, with implications for the star formation and chemical evolution history of the Milky Way. We motivate the need for high-resolution, infrared spectroscopic surveys to reach these faint sources., Comment: 9 pages, submitted to Astro2020 Science White Paper Call

Directly-imaged planetary-mass companions offer unique opportunities in atmospheric studies of exoplanets. They share characteristics of both brown dwarfs and transiting exoplanets, therefore, are critical for connecting atmospheric characterizations for these objects. Rotational phase mapping is a powerful technique to constrain the condensate cloud properties in ultra-cool atmospheres. Applying this technique to directly-imaged planetary-mass companions will be extremely valuable for constraining cloud models in low mass and surface gravity atmospheres and for determining the rotation rate and angular momentum of substellar companions. Here, we present Hubble Space Telescope Wide Field Camera 3 near-infrared time-resolved photometry for three planetary-mass companions, AB Pic B, 2M0122B, and 2M1207b. Using two-roll differential imaging and hybrid point spread function modeling, we achieve sub-percent photometric precision for all three observations. We find tentative modulations ($<\!\!2\sigma$) for AB Pic B and 2M0122B but cannot reach conclusive results on 2M1207b due to strong systematics. The relatively low significance of the modulation measurements cannot rule out the hypothesis that these planetary-mass companions have the same vertical cloud structures as brown dwarfs. Our rotation rate measurements, combined with archival period measurements of planetary-mass companions and brown dwarfs do not support a universal mass-rotation relation. The high precision of our observations and the high occurrence rates of variable low-surface gravity objects encourage high-contrast time-resolved observations with the James Webb Space Telescope., Comment: Accepted for publication in AAS Journal

We report the results of a $J$-band survey for photometric variability in a sample of young, low-gravity objects using the New Technology Telescope (NTT) and the United Kingdom InfraRed Telescope (UKIRT). Surface gravity is a key parameter in the atmospheric properties of brown dwarfs and this is the first large survey that aims to test the gravity dependence of variability properties. We do a full analysis of the spectral signatures of youth and assess the group membership probability of each target using membership tools from the literature. This results in a 30 object sample of young low-gravity brown dwarfs. Since we are lacking in objects with spectral types later than L9, we focus our statistical analysis on the L0-L8.5 objects. We find that the variability occurrence rate of L0-L8.5 low-gravity brown dwarfs in this survey is $30^{+16}_{-8}\%$. We reanalyse the results of Radigan 2014 and find that the field dwarfs with spectral types L0-L8.5 have a variability occurrence rate of $11^{+13}_{-4}\%$. We determine a probability of $98\%$ that the samples are drawn from different distributions. This is the first quantitative indication that the low-gravity objects are more likely to be variable than the field dwarf population. Furthermore, we present follow-up $J_S$ and $K_S$ observations of the young, planetary-mass variable object PSO 318.5-22 over three consecutive nights. We find no evidence of phase shifts between the $J_S$ and $K_S$ bands and find higher $J_S$ amplitudes. We use the $J_S$ lightcurves to measure a rotational period of $8.45\pm0.05~$hr for PSO 318.5-22., Comment: accepted for publication in MNRAS

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., Comment: 21 pages, 10 figures, accepted to AJ

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., Comment: Accepted to AJ. 22 pages, 11 figures, 3 tables

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

Time-resolved observations of brown dwarfs' rotational modulations provide powerful insights into the properties of condensate clouds in ultra-cool atmospheres. Multi-wavelength light curves reveal cloud vertical structures, condensate particle sizes, and cloud morphology, which directly constrain condensate cloud and atmospheric circulation models. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared G141 taken in six consecutive orbits observations of HN Peg B, an L/T transition brown dwarf companion to a G0V type star. The best-fit sine wave to the $1.1-1.7\mu$m broadband light curve has the amplitude of $1.206\pm0.025\%$ and period of $15.4\pm0.5$ hr. The modulation amplitude has no detectable wavelength dependence except in the 1.4 $\mu$m water absorption band, indicating that the characteristic condensate particle sizes are large ($>1\mu$m). We detect significantly ($4.4\sigma$) lower modulation amplitude in the 1.4$\mu$m water absorption band, and find that HN Peg B's spectral modulation resembles those of early T type brown dwarfs. We also describe a new empirical interpolation method to remove spectral contamination from the bright host star. This method may be applied in other high-contrast time-resolved observations with WFC3., Comment: Accepted for publication in AAS Journal

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

We have discovered a potential T0 $\pm$ 1 subdwarf from a search for sources in the AllWISE2 Motion Survey that do not have counterparts in surveys at shorter wavelengths. With a tangential velocity of $\sim$170 kms, this object - WISE J071121.36-573634.2 - has kinematics that are consistent with the thick-disk population of the Milky Way. Spectral fits suggest a low-metallicity for this object but also allow for the possibility of unresolved multiplicity. If WISE J0711-5736 is indeed an sdT0 dwarf, it would be only the second early-T subdwarf discovered to date., Comment: 11 pages, 6 figures, accepted to AJ

We report results from the preliminary trials of Colibri, a dedicated fast-photometry array for the detection of small Kuiper belt objects through serendipitous stellar occultations. Colibri's novel data processing pipeline analyzed 4000 star hours with two overlapping-field EMCCD cameras, detecting no Kuiper belt objects and one false positive occultation event in a high ecliptic latitude field. No occultations would be expected at these latitudes, allowing these results to provide a control sample for the upcoming main Colibri campaign. The empirical false positive rate found by the processing pipeline is consistent with the 0.002% simulation-determined false positive rate. We also describe Colibri's software design, kernel sets for modeling stellar occultations, and method for retrieving occultation parameters from noisy diffraction curves. Colibri's main campaign will begin in mid-2018, operating at a 40 Hz sampling rate., Comment: Accepted for publication by PASP, Oct 31st 2017

Upon its discovery in 2006, the young L7.5 companion to the solar analog HD 203030 was found to be unusual in being $\approx$200 K cooler than older late-L dwarfs. HD 203030B offered the first clear indication that the effective temperature at the L-to-T spectral type transition depends on surface gravity: now a well-known characteristic of low-gravity ultra-cool dwarfs. An initial age analysis of the G8V primary star indicated that the system was 130--400 Myr old, and so the companion between 12--31 $M_{\rm Jup}$. Using moderate resolution near-infrared spectra of HD 203030B, we now find features of very low gravity comparable to those of 10--150 Myr-old L7--L8 dwarfs. We also obtained more accurate near infrared and {\sl Spitzer}/IRAC photometry, and find a $(J-K)_{\rm MKO}$ color of $2.56\pm0.13$ mag---comparable to those observed in other young planetary-mass objects---and a luminosity of log($L_{\rm bol}/L_{\odot}$)$\,=\,-4.75\pm0.04$ dex. We further reassess the evidence for the young age of the host star, HD 203030, with a more comprehensive analysis of the photometry and updated stellar activity measurements and age calibrations. Summarizing the age diagnostics for both components of the binary, we adopt an age of 100 Myr for HD 203030B and an age range of 30--150 Myr. Using cloudy evolutionary models, the new companion age range and luminosity result in a mass of 11 $M_{\rm Jup}$ with a range of 8--15 $M_{\rm Jup}$, and an effective temperature of $1040\pm50$ K., Comment: 12 pages, 7 figures, accepted for publication in AJ