We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

We present $L$'-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wavefront sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semi-major axis of $20^{+3}_{-4}$~au and PDS 70 c to have a semi-major axis of $34^{+12}_{-6}$~au (95\% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2-3~$R_{Jup}$. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 $M_{\textrm{Jup}}$ and a mean mass accretion rate between $3 \times 10^{-7}$ and $8 \times 10^{-7}~M_{\textrm{Jup}}/\textrm{yr}$. For PDS 70 c, we computed a mass between 1 and 3 $M_{\textrm{Jup}}$ and mean mass accretion rate between $1 \times 10^{-7}$ and $5 \times~10^{-7} M_{\textrm{Jup}}/\textrm{yr}$. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets' SEDs.

We present the results of an independent search of all ~200,000 stars observed over the four year Kepler mission (Q1-Q17) for multiplanet systems, using a three-transit minimum detection criteria to search orbital periods up to hundreds of days. We incorporate both automated and manual triage, and provide estimates of the completeness and reliability of our vetting pipeline. Our search returned 17 planet candidates (PCs) in addition to thousands of known Kepler Objects of Interest (KOIs), with a 98.8% recovery rate of already confirmed planets. We highlight the discovery of one candidate, KIC-7340288 b, that is both rocky (radius < 1.6 $R_{\bigoplus}$) and in the Habitable Zone (insolation between 0.25 and 2.2 times the Earth's insolation). Another candidate is an addition to the already known KOI-4509 system. We also present adaptive optics imaging follow-up for six of our new PCs, two of which reveal a line-of-sight stellar companion within 4''.

Currently, postdoctoral fellow (PDF) researchers in Canada face challenges due to the precarious nature of their employment and their overall low compensation and benefits coverage. This report presents three themes, written as statements of need, to support an inclusive and thriving PDF community. These themes are the need for better terms of employment and conditions, the need for access to grants by non-permanent research staff, and the need for a sustainable PDF hiring model that considers the outcomes for the PDFs. We make six recommendations: R1. PDFs should be hired and compensated as skilled experts in their areas, not as trainees. R2. Standard PDF hiring practices should be revised to be more inclusive of different life circumstances. - R2.1 Allow PDFs the option of part-time employment. - R2.2 Remove years-since-PhD time limits from PDF jobs. - R2.3 Financially support PDF hires for relocation and visa expenses. R3. CASCA should form a committee to advocate for and provide support to astronomy PDFs in Canada. R4. CASCA should encourage universities to create offices dedicated to their PDFs. R5. PDFs and other PhD-holding term researchers with a host institution should be able to compete for and win grants to self-fund their own research. R6. Astronomy in Canada should hire general-purpose continuing support scientist positions instead of term PDFs to fill project or mission-specific requirements. In short, we ask for prioritization of people over production of papers.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

orbitize! is an open-source, object-oriented software package for fitting the orbits of directly-imaged objects. It packages the Orbits for the Impatient (OFTI) algorithm and a parallel-tempered Markov Chain Monte Carlo (MCMC) algorithm into a consistent and intuitive Python API. orbitize! makes it easy to run standard astrometric orbit fits; in less than 10 lines of code, users can read in data, perform one fit using OFTI and another using MCMC, and make two publication-ready figures. Extensive pedagogical tutorials, intended to be navigable by both orbit-fitting novices and seasoned experts, are available on our documentation page. We have designed the orbitize! API to be flexible and easy to use/modify for unique cases. orbitize! was designed by members of the exoplanet imaging community to be a central repository for algorithms, techniques, and know-how developed by this community. We intend for it to continue to expand and change as the field progresses and new techniques are developed, and call for community involvement in this process. Complete and up-to-date documentation is available at orbitize.info.

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

The Kepler, K2 and TESS transit surveys are revolutionizing our understanding of planets orbiting close to their host stars and our understanding of exoplanet systems in general, but there remains a gap in our understanding of wide-orbit planets. This gap in our understanding must be filled if we are to understand planet formation and how it affects exoplanet habitability. We summarize current and planned exoplanet detection programs using a variety of methods: microlensing (including WFIRST), radial velocities, Gaia astrometry, and direct imaging. Finally, we discuss the prospects for joint analyses using results from multiple methods and obstacles that could hinder such analyses. We endorse the findings and recommendations published in the 2018 National Academy report on Exoplanet Science Strategy. This white paper extends and complements the material presented therein.

Microlensing is the only method that can detect and measure mass of wide orbit, low mass, solar system analog exoplanets. Mass measurements of such planets would yield massive science on planet formation, exoplanet demographics, free floating planets, planet frequencies towards the galaxy. High res follow-up observations of past microlens targets provide a mass measurement of microlens planets and hosts at an uncertainty of <20%. This will be primary method for mass measurement with WFIRST. We advocate for the fact that high resolution observations with AO, HST and JWST(in future) remain necessary in coming decade to develop the methods, to determine the field and filter selection, understand the systematics and to develop a robust pipeline to release high quality data products from WFIRST microlensing survey such that the astronomy community can promptly engage in the science. We also support future high res obs with US ELTs with advanced Laser AO systems in context of enhancing the science return of WFIRST microlensing survey. We endorse the 2018 Exoplanet Science Strategy report published by the National Academy. This white paper extends and complements the material presented therein. In particular, this white paper supports the recommendation of the National Academy Exoplanet Science Strategy report that: NASA should launch WFIRST to conduct its microlensing survey of distant planets and to demonstrate the technique of coronagraphic spectroscopy on exoplanet targets. This white paper also supports to the finding from that report which states "A number of activities, including precursor and concurrent observations using ground- and space-based facilities, would optimize the scientific yield of the WFIRST microlensing survey."

We use radial velocity observations to search for long-period gas giant companions in systems hosting inner super-Earth (1-4 R_Earth, 1-10 M_Earth) planets to constrain formation and migration scenarios for this population. We consistently re-fit published RV datasets for 65 stars and find 9 systems with statistically significant trends indicating the presence of an outer companion. We combine these RV data with AO images to constrain the masses and semi-major axes of these companions. We quantify our sensitivity to the presence of long-period companions by fitting the sample with a power law distribution and find an occurrence rate of 39+/-7% for companions 0.5-20 M_Jup and 1-20 AU. Half of our systems were discovered by the transit method and half were discovered by the RV method. While differences in RV baselines and number of data points between the two samples lead to different sensitivities to distant companions, we find that occurrence rates of gas giant companions in each sample are consistent at the 0.5$\sigma$ level. We compare the frequency of Jupiter analogs in these systems to the equivalent rate from field star surveys and find that Jupiter analogs are more common around stars hosting super-Earths. We conclude that the presence of outer gas giants does not suppress the formation of inner super-Earths, and that these two populations of planets instead appear to be correlated. We also find that the stellar metallicities of systems with gas giant companions are higher than those without companions, in agreement with the well-established metallicity correlation from RV surveys of field stars.

Reference star differential imaging (RDI) is a powerful strategy for high contrast imaging. Using example observations taken with the vortex coronagraph mode of Keck/NIRC2 in $L^\prime$ band, we demonstrate that RDI provides improved sensitivity to point sources at small angular separations compared to angular differential imaging (ADI). Applying RDI to images of the low-mass stellar companions HIP 79124 C (192 mas separation, $\Delta L^\prime$=4.01) and HIP 78233 B (141 mas separation, $\Delta L^\prime$=4.78), the latter a first imaging detection, increases the significance of their detections by up to a factor of 5 with respect to ADI. We compare methods for reference frames selection and find that pre-selection of frames improves detection significance of point sources by up to a factor of 3. In addition, we use observations of the circumstellar disks around MWC 758 and 2MASS J16042165-2130284 to show that RDI allows for accurate mapping of scattered light distributions without self-subtraction artifacts.

The Coronagraph Instrument (CGI) for NASA's Wide Field Infrared Survey Telescope (WFIRST) will constitute a dramatic step forward for high-contrast imaging, integral field spectroscopy, and polarimetry of exoplanets and circumstellar disks, aiming to improve upon the sensitivity of current ground-based direct imaging facilities by 2-3 orders of magnitude. Furthermore, CGI will serve as a pathfinder for future exo-Earth imaging and characterization missions by demonstrating wavefront control, coronagraphy, and spectral retrieval in a new contrast regime, and by validating instrument and telescope models at unprecedented levels of precision. To achieve this jump in performance, it is critical to draw on the experience of ground-based high-contrast facilities. We discuss several areas of relevant commonalities, including: wavefront control, post-processing of integral field unit data, and calibration and observing strategies.

The NIRC2 vortex coronagraph is an instrument on Keck II designed to directly image exoplanets and circumstellar disks at mid-infrared bands $L^\prime$ (3.4-4.1 $\mu$m) and $M_s$ (4.55-4.8 $\mu$m). We analyze imaging data and corresponding adaptive optics telemetry, observing conditions, and other metadata over a three year time period to characterize the performance of the instrument and predict the detection limits of future observations. We systematically process images from 359 observations of 304 unique stars to subtract residual starlight (i.e., the coronagraphic point spread function) of the target star using two methods: angular differential imaging (ADI) and reference star differential imaging (RDI). We find that for the typical parallactic angle (PA) rotation of our dataset ($\sim$10$^{\circ}$), RDI provides gains over ADI for angular separations smaller than 0.25$^{\prime\prime}$. Furthermore, we find a power-law relation between the angular separation from the host star and the minimum PA rotation required for ADI to outperform RDI, with a power-law index of -1.18$\pm$0.08. Finally, we use random forest models to estimate ADI and RDI post-processed detection limits a priori. These models, which we provide publicly on a website, explain 70%-80% of the variance in ADI detection limits and 30%-50% of the variance in RDI detection limits. Averaged over a range of angular separations, our models predict both ADI and RDI contrast to within a factor of 2. These results illuminate important factors in high-contrast imaging observations with the NIRC2 vortex coronagraph, help improve observing strategies, and inform future upgrades to the hardware.

We have discovered a doubly eclipsing, bound, quadruple star system in the field of K2 Campaign 7. EPIC 219217635 is a stellar image with $Kp = 12.7$ that contains an eclipsing binary (`EB') with $P_A = 3.59470$ d and a second EB with $P_B = 0.61825$ d. We have obtained followup radial-velocity (`RV') spectroscopy observations, adaptive optics imaging, as well as ground-based photometric observations. From our analysis of all the observations, we derive good estimates for a number of the system parameters. We conclude that (1) both binaries are bound in a quadruple star system; (2) a linear trend to the RV curve of binary A is found over a 2-year interval, corresponding to an acceleration, $\dot \gamma = 0.0024 \pm 0.0007$ cm s$^{-2}$; (3) small irregular variations are seen in the eclipse-timing variations (`ETVs') detected over the same interval; (4) the orbital separation of the quadruple system is probably in the range of 8-25 AU; and (5) the orbital planes of the two binaries must be inclined with respect to each other by at least 25$^\circ$. In addition, we find that binary B is evolved, and the cooler and currently less massive star has transferred much of its envelope to the currently more massive star. We have also demonstrated that the system is sufficiently bright that the eclipses can be followed using small ground-based telescopes, and that this system may be profitably studied over the next decade when the outer orbit of the quadruple is expected to manifest itself in the ETV and/or RV curves.

Our Keck/NIRC2 imaging survey searches for stellar companions around 144 systems with radial velocity (RV) detected giant planets to determine whether stellar binaries influence the planets' orbital parameters. This survey, the largest of its kind to date, finds eight confirmed binary systems and three confirmed triple systems. These include three new multi-stellar systems (HD 30856, HD 86081, and HD 207832) and three multi-stellar systems with newly confirmed common proper motion (HD 43691, HD 116029, and HD 164509). We combine these systems with seven RV planet-hosting multi-stellar systems from the literature in order to test for differences in the properties of planets with semimajor axes ranging between 0.1-5 au in single vs multi-stellar systems. We find no evidence that the presence or absence of stellar companions alters the distribution of planet properties in these systems. Although the observed stellar companions might influence the orbits of more distant planetary companions in these systems, our RV observations currently provide only weak constraints on the masses and orbital properties of planets beyond 5 au. In order to aid future efforts to characterize long period RV companions in these systems, we publish our contrast curves for all 144 targets. Using four years of astrometry for six hierarchical triple star systems hosting giant planets, we fit the orbits of the stellar companions in order to characterize the orbital architecture in these systems. We find that the orbital plane of the secondary and tertiary companions are inconsistent with an edge-on orbit in four out of six cases.

We announce the discovery of KELT-16b, a highly irradiated, ultra-short period hot Jupiter transiting the relatively bright ($V = 11.7$) star TYC 2688-1839-1. A global analysis of the system shows KELT-16 to be an F7V star with $T_\textrm{eff} = 6236\pm54$ K, $\log{g_\star} = 4.253_{-0.036}^{+0.031}$, [Fe/H] = -0.002$_{-0.085}^{+0.086}$, $M_\star = 1.211_{-0.046}^{+0.043} M_\odot$, and $R_\star = 1.360_{-0.053}^{+0.064} R_\odot$. The planet is a relatively high mass inflated gas giant with $M_\textrm{P} = 2.75_{-0.15}^{+0.16} M_\textrm{J}$, $R_\textrm{P} = 1.415_{-0.067}^{+0.084} R_\textrm{J}$, density $\rho_\textrm{P} = 1.20\pm0.18$ g cm$^{-3}$, surface gravity $\log{g_\textrm{P}} = 3.530_{-0.049}^{+0.042}$, and $T_\textrm{eq} = 2453_{-47}^{+55}$ K. The best-fitting linear ephemeris is $T_\textrm{C} = 2457247.24791\pm0.00019$ BJD$_{tdb}$ and $P = 0.9689951 \pm 0.0000024$ d. KELT-16b joins WASP-18b, -19b, -43b, -103b, and HATS-18b as the only giant transiting planets with $P < 1$ day. Its ultra-short period and high irradiation make it a benchmark target for atmospheric studies by HST, Spitzer, and eventually JWST. For example, as a hotter, higher mass analog of WASP-43b, KELT-16b may feature an atmospheric temperature-pressure inversion and day-to-night temperature swing extreme enough for TiO to rain out at the terminator. KELT-16b could also join WASP-43b in extending tests of the observed mass-metallicity relation of the Solar System gas giants to higher masses. KELT-16b currently orbits at a mere $\sim$ 1.7 Roche radii from its host star, and could be tidally disrupted in as little as a few $\times 10^{5}$ years (for a stellar tidal quality factor of $Q_*' = 10^5$). Finally, the likely existence of a widely separated bound stellar companion in the KELT-16 system makes it possible that Kozai-Lidov oscillations played a role in driving KELT-16b inward to its current precarious orbit.

A vortex coronagraph is now available for high contrast observations with the Keck/NIRC2 instrument at L band. Reaching the optimal performance of the coronagraph requires fine control of the wavefront incident on the phase mask. In particular, centering errors can lead to significant stellar light leakage that degrades the contrast performance and prevents the observation of faint planetary companions around the observed stars. It is thus critical to correct for the possible slow drift of the star image from the phase mask center, generally due to mechanical flexures induced by temperature and/or gravity field variation, or to misalignment between the optics that rotate in pupil tracking mode. A control loop based on the QACITS algorithm for the vortex coronagraph has thus been developed and deployed for the Keck/NIRC2 instrument. This algorithm executes the entire observing sequence, including the calibration steps, initial centering of the star on the vortex center and stabilisation during the acquisition of science frames. On-sky data show that the QACITS control loop stabilizes the position of the star image down to 2.4 mas rms at a frequency of about 0.02 Hz. However, the accuracy of the estimator is probably limited by a systematic error due to a misalignment of the Lyot stop with respect to the entrance pupil, estimated to be on the order of 4.5 mas. A method to reduce the amplitude of this bias down to 1 mas is proposed. The QACITS control loop has been successfully implemented and provides a robust method to center and stabilize the star image on the vortex mask. In addition, QACITS ensures a repeatable pointing quality and significantly improves the observing efficiency compared to manual operations. It is now routinely used for vortex coronagraph observations at Keck/NIRC2, providing contrast and angular resolution capabilities suited for exoplanet and disk imaging.

We present a strongly interacting quadruple system associated with the K2 target EPIC 220204960. The K2 target itself is a Kp = 12.7 magnitude star at Teff ~ 6100 K which we designate as "B-N" (blue northerly image). The host of the quadruple system, however, is a Kp = 17 magnitude star with a composite M-star spectrum, which we designate as "R-S" (red southerly image). With a 3.2" separation and similar radial velocities and photometric distances, 'B-N' is likely physically associated with 'R-S', making this a quintuple system, but that is incidental to our main claim of a strongly interacting quadruple system in 'R-S'. The two binaries in 'R-S' have orbital periods of 13.27 d and 14.41 d, respectively, and each has an inclination angle of >89 degrees. From our analysis of radial velocity measurements, and of the photometric lightcurve, we conclude that all four stars are very similar with masses close to 0.4 Msun. Both of the binaries exhibit significant ETVs where those of the primary and secondary eclipses 'diverge' by 0.05 days over the course of the 80-day observations. Via a systematic set of numerical simulations of quadruple systems consisting of two interacting binaries, we conclude that the outer orbital period is very likely to be between 300 and 500 days. If sufficient time is devoted to RV studies of this faint target, the outer orbit should be measurable within a year.

We present here our observations and analysis of the dayside emission spectrum of the hot Jupiter WASP-103b. We observed WASP-103b during secondary eclipse using two visits of the Hubble Space Telescope with the G141 grism on Wide Field Camera 3 in spatial scan mode. We generated secondary eclipse light curves of the planet in both blended white-light and spectrally binned wavechannels from 1.1-1.7 micron and corrected the light curves for flux contamination from a nearby companion star. We modeled the detector systematics and secondary eclipse spectrum using Gaussian process regression and found that the near-IR emission spectrum of WASP-103b is featureless across the observed near-IR region to down to a sensitivity of 175 ppm, and shows a shallow slope towards the red. The atmosphere has a single brightness temperature of T_B = 2890 K across this wavelength range. This region of the spectrum is indistinguishable from isothermal, but may not manifest from a physically isothermal system, i.e. pseudo-isothermal. A Solar-metallicity profile with a thermal inversion layer at 10^-2 bar fits WASP-103b's spectrum with high confidence, as do an isothermal profile with Solar metallicity and a monotonically decreasing atmosphere with C/O>1. The data rule out a monotonically decreasing atmospheric profile with Solar composition, and we rule out a low-metallicity decreasing profile as non-physical for this system. The pseudo-isothermal profile could be explained by a thermal inversion layer just above the layer probed by our observations, or by clouds or haze in the upper atmosphere. Transmission spectra at optical wavelengths would allow us to better differentiate between potential atmospheric models.

This paper presents the development of an in situ measurement system known as the Debris Resistive Acoustic Grid Orbital Navy/NASA Sensor (DRAGONS). The DRAGONS system is designed to detect impacts caused by particles ranging from 50 micrometers to 1 mm at both low-earth and geostationary orbits. DRAGONS utilizes a combination of low-cost sensor technologies to facilitate accurate measurements and approximations of the size, velocity, and angle of impacting micrometeoroids and orbital debris (MMOD). Two thin layers of kapton sheets with resistive traces are used to detect the changes in resistance that are directly proportional to the impacting force caused by the fast traveling particles. Four polyvinylidene fluoride-based sensors are positioned in the back of each kapton sheet to measure acoustic strain caused by an impact. The electronic hardware module that controls all operations employs a low-power, modular, and compact design that enables it to be installed as a low-resource load on a host satellite. Laboratory results demonstrate that in addition to having the ability to detect an impact event, the DRAGONS system can determine impact location, speed, and angle of impact with a mean error of 1.4 cm, 0.2 km/s, and 5{\deg}. The DRAGONS system could be deployed as an add-on subsystem of a payload to enable a real-time, in-depth study of the properties of MMOD.

We report the discovery of an extremely red planetary-mass companion to 2MASS J22362452+4751425, a $\approx$0.6 $M_{\odot}$ late-K dwarf likely belonging to the $\sim$120 Myr AB Doradus moving group. 2M2236+4751 b was identified in multi-epoch NIRC2 adaptive optics imaging at Keck Observatory at a separation of 3.7$"$, or 230 $\pm$ 20 AU in projection at the kinematic distance of 63 $\pm$ 5 pc to its host star. Assuming membership in the AB Dor group, as suggested from its kinematics, the inferred mass of 2M2236+4751 b is 11-14 $M_\mathrm{Jup}$. Follow-up Keck/OSIRIS $K$-band spectroscopy of the companion reveals strong CO absorption similar to other faint red L dwarfs and lacks signs of methane absorption despite having an effective temperature of $\approx$900-1200 K. With a ($J$-$K$)$_\mathrm{MKO}$ color of 2.69 $\pm$ 0.12 mag, the near-infrared slope of 2M2236+4751 b is redder than all of the HR 8799 planets and instead resembles the $\approx$23 Myr isolated planetary-mass object PSO J318.5-22, implying that similarly thick photospheric clouds can persist in the atmospheres of giant planets at ages beyond 100 Myr. In near-infrared color-magnitude diagrams, 2M2236+4751 b is located at the tip of the red L dwarf sequence and appears to define the "elbow" of the AB Dor substellar isochrone separating low-gravity L dwarfs from the cooler young T dwarf track. 2M2236+4751 b is the reddest substellar companion to a star and will be a valuable benchmark to study the shared atmospheric properties of young low-mass brown dwarfs and extrasolar giant planets.

Small inner working angle coronagraphs are essential to benefit from the full potential of large and future extremely large ground-based telescopes, especially in the context of the detection and characterization of exoplanets. Among existing solutions, the vortex coronagraph stands as one of the most effective and promising solutions. However, for focal-plane coronagraph, a small inner working angle comes necessarily at the cost of a high sensitivity to pointing errors. This is the reason why a pointing control system is imperative to stabilize the star on the vortex center against pointing drifts due to mechanical flexures, that generally occur during observation due for instance to temperature and/or gravity variations. We have therefore developed a technique called QACITS (Quadrant Analysis of Coronagraphic Images for Tip-tilt Sensing), which is based on the analysis of the coronagraphic image shape to infer the amount of pointing error. It has been shown that the flux gradient in the image is directly related to the amount of tip-tilt affecting the beam. The main advantage of this technique is that it does not require any additional setup and can thus be easily implemented on all current facilities equipped with a vortex phase mask. In this paper, we focus on the implementation of the QACITS sensor at Keck/NIRC2, where an L-band AGPM has been recently commissioned (June and October 2015), successfully validating the QACITS estimator in the case of a centrally obstructed pupil. The algorithm has been designed to be easily handled by any user observing in vortex mode, which is available for science in shared risk mode since 2016B.

Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that $47\%\pm7\%$ of hot Jupiter systems have stellar companions with semi-major axes between 50-2000 AU. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of $4.4\sigma$. In the 1-50AU range, only $3.9^{+4.5}_{-2.0}\%$ of hot Jupiters host stellar companions compared to the field star value of $16.4\%\pm0.7\%$, which is a $2.7\sigma$ difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai-Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai-Lidov oscillations assuming initial semi-major axes between 1-5 AU, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

We present a quintuple star system that contains two eclipsing binaries. The unusual architecture includes two stellar images separated by 11" on the sky: EPIC 212651213 and EPIC 212651234. The more easterly image (212651213) actually hosts both eclipsing binaries which are resolved within that image at 0.09", while the westerly image (212651234) appears to be single in adaptive optics (AO), speckle imaging, and radial velocity (RV) studies. The 'A' binary is circular with a 5.1-day period, while the 'B' binary is eccentric with a 13.1-day period. The gamma velocities of the A and B binaries are different by ~10 km/s. That, coupled with their resolved projected separation of 0.09", indicates that the orbital period and separation of the 'C' binary (consisting of A orbiting B) are ~65 years and ~25 AU, respectively, under the simplifying assumption of a circular orbit. Motion within the C orbit should be discernible via future RV, AO, and speckle imaging studies within a couple of years. The C system (i.e., 212651213) has a radial velocity and proper motion that differ from that of 212651234 by only ~1.4 km/s and ~3 mas/yr. This set of similar space velocities in 3 dimensions strongly implies that these two objects are also physically bound, making this at least a quintuple star system.

We report the discovery of EPIC201702477b, a transiting brown dwarf in a long period (40.73691 +/- 0.00037 day) and eccentric (e=0.2281 +/- 0.0026) orbit. This system was initially reported as a planetary candidate based on two transit events seen in K2 Campaign 1 photometry and later validated as an exoplanet. We confirm the transit and refine the ephemeris with two subsequent ground-based detections of the transit using the LCOGT 1m telescope network. We rule out any transit timing variations above the level of 30s. Using high precision radial velocity measurements from HARPS and SOPHIE we identify the transiting companion as a brown dwarf with a mass, radius, and bulk density of 66.9 +/- 1.7 M$_J$, 0.757 +/- 0.065 R$_J$, and 191+/-51 g.cm$^{-3}$ respectively. EPIC201702477b is the smallest radius brown dwarf yet discovered, with a mass just below the H-burning limit. It has the highest density of any planet, substellar mass object or main-sequence star discovered so far. We find evidence in the set of known transiting brown dwarfs for two populations of objects - high mass brown dwarfs and low mass brown dwarfs. The higher-mass population have radii in very close agreement to theoretical models, and show a lower-mass limit around 60 M$_J$. This may be the signature of mass-dependent ejection of systems during the formation process.

We conducted a Doppler survey at Keck combined with NIRC2 K-band AO imaging to search for massive, long-period companions to 123 known exoplanet systems with one or two planets detected using the radial velocity (RV) method. Our survey is sensitive to Jupiter mass planets out to 20 AU for a majority of stars in our sample, and we report the discovery of eight new long-period planets, in addition to 20 systems with statistically significant RV trends indicating the presence of an outer companion beyond 5 AU. We combine our RV observations with AO imaging to determine the range of allowed masses and orbital separations for these companions, and account for variations in our sensitivity to companions among stars in our sample. We estimate the total occurrence rate of companions in our sample to be 52 +/- 5% over the range 1 - 20 M_Jup and 5 - 20 AU. Our data also suggest a declining frequency for gas giant planets in these systems beyond 3-10 AU, in contrast to earlier studies that found a rising frequency for giant planets in the range 0.01-3 AU. This suggests either that the frequency of gas giant planets peaks between 3-10 AU, or that outer companions in these systems have a different semi-major axis distribution than the overall gas giant planet population. Our results also suggest that hot gas giants may be more likely to have an outer companion than cold gas giants. We find that planets with an outer companion have higher average eccentricities than their single counterparts, suggesting that dynamical interactions between planets may play an important role in these systems.

KIC 7177553 was observed by the Kepler satellite to be an eclipsing eccentric binary star system with an 18-day orbital period. Recently, an eclipse timing study of the Kepler binaries has revealed eclipse timing variations in this object with an amplitude of about 100 sec, and an outer period of 529 days. The implied mass of the third body is that of a superJupiter, but below the mass of a brown dwarf. We therefore embarked on a radial velocity study of this binary to determine its system configuration and to check the hypothesis that it hosts a giant planet. From the radial velocity measurements, it became immediately obvious that the same Kepler target contains another eccentric binary, this one with a 16.5-day orbital period. Direct imaging using adaptive optics reveals that the two binaries are separated by 0.4 arcsec (about 167 AU), and have nearly the same magnitude (to within 2%). The close angular proximity of the two binaries, and very similar Gamma velocities, strongly suggest that KIC 7177553 is one of the rare SB4 systems consisting of two eccentric binaries where at least one system is eclipsing. Both systems consist of slowly rotating, non-evolved, solar-like stars of comparable masses. From the orbital separation and the small difference in Gamma velocity, we infer that the period of the outer orbit most likely lies in the range 1000 to 3000 years. New images taken over the next few years, as well as the high-precision astrometry of the Gaia satellite mission, will allow us to set much narrower constraints on the system geometry. Finally, we note that the observed eclipse timing variations in the Kepler data cannot be produced by the second binary. Further spectroscopic observations on a longer time scale will be required to prove the existence of the massive planet.

Surveys of nearby field stars indicate that stellar binaries are common, yet little is known about the effects that these companions may have on planet formation and evolution. The Friends of Hot Jupiters project uses three complementary techniques to search for stellar companions to known planet-hosting stars: radial velocity monitoring, adaptive optics imaging, and near-infrared spectroscopy. In this paper, we examine high-resolution K band infrared spectra of fifty stars hosting gas giant planets on short-period orbits. We use spectral fitting to search for blended lines due to the presence of cool stellar companions in the spectra of our target stars, where we are sensitive to companions with temperatures between 3500-5000 K and projected separations less than 100 AU in most systems. We identify eight systems with candidate low-mass companions, including one companion that was independently detected in our AO imaging survey. For systems with radial velocity accelerations, a spectroscopic non-detection rules out scenarios involving a stellar companion in a high inclination orbit. We use these data to place an upper limit on the stellar binary fraction at small projected separations, and show that the observed population of candidate companions is consistent with that of field stars and also with the population of wide-separation companions detected in our previous AO survey. We find no evidence that spectroscopic stellar companions are preferentially located in systems with short-period gas giant planets on eccentric and/or misaligned orbits.

Multi-star systems are common, yet little is known about a stellar companion's influence on the formation and evolution of planetary systems. For instance, stellar companions may have facilitated the inward migration of hot Jupiters towards to their present day positions. Many observed short period gas giant planets also have orbits that are misaligned with respect to their star's spin axis, which has also been attributed to the presence of a massive outer companion on a non-coplanar orbit. We present the results of a multi-band direct imaging survey using Keck NIRC2 to measure the fraction of short period gas giant planets found in multi-star systems. Over three years, we completed a survey of 50 targets ("Friends of Hot Jupiters") with 27 targets showing some signature of multi-body interaction (misaligned or eccentric orbits) and 23 targets in a control sample (well-aligned and circular orbits). We report the masses, projected separations, and confirmed common proper motion for the 19 stellar companions found around 17 stars. Correcting for survey incompleteness, we report companion fractions of $48\%\pm9\%$, $47\%\pm12\%$, and $51\%\pm13\%$ in our total, misaligned/eccentric, and control samples, respectively. This total stellar companion fraction is $2.8\,\sigma$ larger than the fraction of field stars with companions approximately $50-2000\,$AU. We observe no correlation between misaligned/eccentric hot Jupiter systems and the incidence of stellar companions. Combining this result with our previous radial velocity survey, we determine that $72\% \pm 16\%$ of hot Jupiters are part of multi-planet and/or multi-star systems.

We report secondary eclipse photometry of the hot Jupiter HAT-P-32Ab, taken with Hale/WIRC in H and Ks bands and with Spitzer/IRAC at 3.6 and 4.5 micron. We carried out adaptive optics imaging of the planet host star HAT-P-32A and its companion HAT-P-32B in the near-IR and the visible. We clearly resolve the two stars from each other and find a separation of 2.923" +/- 0. 004" and a position angle 110.64 deg +/- 0.12 deg. We measure the flux ratios of the binary in g' r' i' z' and H & Ks bands, and determine Teff = 3565 +/- 82 K for the companion star, corresponding to an M1.5 dwarf. We use PHOENIX stellar atmosphere models to correct the dilution of the secondary eclipse depths of the hot Jupiter due to the presence of the M1.5 companion. We also improve the secondary eclipse photometry by accounting for the non-classical, flux-dependent nonlinearity of the WIRC IR detector in the H band. We measure planet-to-star flux ratios of 0.090 +/- 0.033%, 0.178 +/- 0.057%, 0.364 +/- 0.016%, and 0.438 +/- 0.020% in the H, Ks, 3.6 and 4.5 micron bands, respectively. We compare these with planetary atmospheric models, and find they prefer an atmosphere with a temperature inversion and inefficient heat redistribution. However, we also find that the data are equally well-described by a blackbody model for the planet with Tp = 2042 +/- 50 K. Finally, we measure a secondary eclipse timing offset of 0.3 +/- 1.3 min from the predicted mid-eclipse time, which constrains e = 0.0072 +0.0700/-0.0064 when combined with RV data and is more consistent with a circular orbit.

We present high spatial resolution images that demonstrate the hot Jupiters WASP-12b and HAT-P-8b orbit the primary star of hierarchical triple star systems. In each case, two distant companions with colors and brightness consistent with M dwarfs co-orbit the planet host as well as one another. Our adaptive optics images spatially resolve the secondary around WASP-12, previously identified by Bergfors et al. 2011 and Crossfield et al. 2012, into two distinct sources separated by 84.3+/-0.6 mas (21 +/- 3 AU). We find that the secondary to HAT-P-8, also identified by Bergfors et al. 2011, is in fact composed of two stars separated by 65.3+/-0.5 mas (15+/-1 AU). Our follow-up observations demonstrate physical association through common proper-motion. HAT-P-8 C has a particularly low mass, which we estimate to be 0.18+/-0.02Msun using photometry. Due to their hierarchy, WASP-12 BC and HAT-P-8 BC will enable the first dynamical mass determination for hot Jupiter stellar companions. These previously well-studied planet hosts now represent higher-order multi-star systems with potentially complex dynamics, underscoring the importance of diffraction-limited imaging and providing additional context for understanding the migrant population of transiting hot Jupiters.

In this paper we search for distant massive companions to known transiting gas giant planets that may have influenced the dynamical evolution of these systems. We present new radial velocity observations for a sample of 51 planets obtained using the Keck HIRES instrument, and find statistically significant accelerations in fifteen systems. Six of these systems have no previously reported accelerations in the published literature: HAT-P-10, HAT-P-22, HAT-P-29, HAT-P-32, WASP-10, and XO-2. We combine our radial velocity fits with Keck NIRC2 adaptive optics (AO) imaging data to place constraints on the allowed masses and orbital periods of the companions responsible for the detected accelerations. The estimated masses of the companions range between 1-500 M_Jup, with orbital semi-major axes typically between 1-75 AU. A significant majority of the companions detected by our survey are constrained to have minimum masses comparable to or larger than those of the transiting planets in these systems, making them candidates for influencing the orbital evolution of the inner gas giant. We estimate a total occurrence rate of 51 +/- 10% for companions with masses between 1-13 M_Jup and orbital semi-major axes between 1-20 AU in our sample. We find no statistically significant difference between the frequency of companions to transiting planets with misaligned or eccentric orbits and those with well-aligned, circular orbits. We combine our expanded sample of radial velocity measurements with constraints from transit and secondary eclipse observations to provide improved measurements of the physical and orbital characteristics of all of the planets included in our survey.

We present observations at 250, 350, and 500 um of the nearby galaxy cluster Abell 3112 (z=0.075) carried out with BLAST, the Balloon-borne Large Aperture Submillimeter Telescope. Five cluster members are individually detected as bright submillimetre sources. Their far-infrared SEDs and optical colours identify them as normal star-forming galaxies of high mass, with globally evolved stellar populations. They all have B-R colours of 1.38+/-0.08, transitional between the blue, active population and the red, evolved galaxies that dominate the cluster core. We stack to determine the mean submillimetre emission from all cluster members, which is determined to be 16.6+/-2.5, 6.1+/-1.9, and 1.5+/-1.3 mJy at 250, 350, and 500 um, respectively. Stacking analyses of the submillimetre emission of cluster members reveal trends in the mean far-infrared luminosity with respect to cluster-centric radius and Ks-band magnitude. We find that a large fraction of submillimetre emission comes from the boundary of the inner, virialized region of the cluster, at cluster-centric distances around R_500. Stacking also shows that the bulk of the submillimetre emission arises in intermediate-mass galaxies (L<L*), with Ks magnitude ~1 mag fainter than the giant ellipticals. The results and constraints obtained in this work will provide a useful reference for the forthcoming surveys to be conducted on galaxy clusters by Herschel.

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has made one square degree, deep, confusion limited maps at three different bands, centered on the Great Observatories Origins Deep Survey South field. By calculating the covariance of these maps with catalogs of 24 micron sources from the Far-Infrared Deep Extragalactic Legacy Survey (FIDEL), we have determined that the total submillimeter intensities are 8.60 +/- 0.59, 4.93 +/- 0.34, and 2.27 +/- 0.20 nW m^2 sr^(-1) at 250, 350, and 500 micron, respectively. These numbers are more precise than previous estimates of the cosmic infrared background (CIB) and are consistent with 24 micron-selected galaxies generating the full intensity of the CIB. We find that the fraction of the CIB that originates from sources at z >= 1.2 increases with wavelength, with 60% from high redshift sources at 500 micron. At all BLAST wavelengths, the relative intensity of high-z sources is higher for 24 micron-faint sources than it is for 24 micron-bright sources. Galaxies identified as active galactic nuclei (AGN) by their Infrared Array Camera (IRAC) colors are 1.6-2.6 times brighter than the average population at 250-500 micron, consistent with what is found for X-ray-selected AGN. BzK-selected galaxies are found to be moderately brighter than typical 24 micron-selected galaxies in the BLAST bands. These data provide high precision constraints for models of the evolution of the number density and intensity of star forming galaxies at high redshift.

We directly measure redshift evolution in the mean physical properties (far-infrared luminosity, temperature, and mass) of the galaxies that produce the cosmic infrared background (CIB), using measurements from the Balloon-borne Large Aperture Sub-millimeter Telescope (BLAST), and Spitzer which constrain the CIB emission peak. This sample is known to produce a surface brightness in the BLAST bands consistent with the full CIB, and photometric redshifts are identified for all of the objects. We find that most of the 70 micron background is generated at z <~ 1 and the 500 micron background generated at z >~ 1. A significant growth is observed in the mean luminosity from ~ 10^9 - 10^12 L_sun, and in the mean temperature by 10 K, from redshifts 0< z < 3. However, there is only weak positive evolution in the comoving dust mass in these galaxies across the same redshift range. We also measure the evolution of the far-infrared luminosity density, and the star-formation rate history for these objects, finding good agreement with other infrared studies up to z ~1, exceeding the contribution attributed to optically-selected galaxies.

Submillimetre surveys during the past decade have discovered a population of luminous, high-redshift, dusty starburst galaxies. In the redshift range 1 <= z <= 4, these massive submillimetre galaxies go through a phase characterized by optically obscured star formation at rates several hundred times that in the local Universe. Half of the starlight from this highly energetic process is absorbed and thermally re-radiated by clouds of dust at temperatures near 30 K with spectral energy distributions peaking at 100 microns in the rest frame. At 1 <= z <= 4, the peak is redshifted to wavelengths between 200 and 500 microns. The cumulative effect of these galaxies is to yield extragalactic optical and far-infrared backgrounds with approximately equal energy densities. Since the initial detection of the far-infrared background (FIRB), higher-resolution experiments have sought to decompose this integrated radiation into the contributions from individual galaxies. Here we report the results of an extragalactic survey at 250, 350 and 500 microns. Combining our results at 500 microns with those at 24 microns, we determine that all of the FIRB comes from individual galaxies, with galaxies at z >= 1.2 accounting for 70 per cent of it. As expected, at the longest wavelengths the signal is dominated by ultraluminous galaxies at z > 1.