We present the result of a search for orphan Gamma-Ray Burst (GRB) afterglows in the Pan-STARRS1 (PS1) COSMOS survey. There is extensive theoretical and observational evidence suggesting that GRBs are collimated jets; the direct observation of orphan GRB afterglows would further support this model. An optimal survey strategy is designed by coupling the PS1 survey with the Subaru/Hyper-Suprime-Cam (HSC) survey. The PS1 COSMOS survey, one of the survey fields in the PS1 Medium Deep Survey (PS1/MDS), searches a field of 7 deg$^2$ from December 2011 to January 2014, reaching a limiting magnitude R $\sim$ 23. The dense cadence of PS1/MDS is crucial for identifying transients, and the deep magnitude reached by the HSC survey (R $\sim$ 26) is important for evaluating potential GRB hosts. A transient classification method is employed to select potential orphan GRB afterglow candidates. After a thorough analysis of the transient and host galaxy properties, we conclude that there are no candidates in this survey field. The null result implies that the consideration of jet structures is essential for further orphan GRB afterglow surveys.

The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $\sigma_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $\sigma_e$ of 0.045 and a typical sin(i) x gaussian inclination distribution with a width $\sigma_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 member and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range.

Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

We constrain the rest-frame FUV (1546\AA), NUV (2345\AA) and U-band (3690\AA) luminosity functions (LFs) and luminosity densities (LDs) with unprecedented precision from $z\sim0.2$ to $z\sim3$ (FUV, NUV) and $z\sim2$ (U-band). Our sample of over 4.3 million galaxies, selected from the CFHT Large Area $U$-band Deep Survey (CLAUDS) and HyperSuprime-Cam Subaru Strategic Program (HSC-SSP) data lets us probe the very faint regime (down to $M_\mathrm{FUV},M_\mathrm{NUV},M_\mathrm{U} \simeq -15$ at low redshift) while simultaneously detecting very rare galaxies at the bright end down to comoving densities $<10^{-5}$ Mpc$^{-3}$. Our FUV and NUV LFs are well fitted by single Schechter functions, with faint-end slopes that are very stable up to $z\sim2$. We confirm, but self-consistently and with much better precision than previous studies, that the LDs at all three wavelengths increase rapidly with lookback time to $z\sim1$, and then much more slowly at $1<z<2$--$3$. Evolution of the FUV and NUV LFs and LDs at $z<1$ is driven almost entirely by the fading of the characteristic magnitude, $M^\star_{UV}$, while at $z>1$ it is due to the evolution of both $M^\star_{UV}$ and the characteristic number density $\phi^\star_{UV}$. In contrast, the U-band LF has an excess of faint galaxies and is fitted with a double-Schechter form; $M^\star_\mathrm{U}$, both $\phi^\star_\mathrm{U}$ components, and the bright-end slope evolve throughout $0.2<z<2$, while the faint-end slope is constant over at least the measurable $0.05<z<0.6$. We present tables of our Schechter parameters and LD measurements that can be used for testing theoretical galaxy evolution models and forecasting future observations.

The surface characterization of Trans-Neptunian Binaries (TNBs) is key to understanding the properties of the disk of planetesimals from which these objects formed. In the optical wavelengths, it has been demonstrated that most equal-sized component systems share similar colors, suggesting they have a similar composition. The color homogeneity of binary pairs contrasts with the overall diversity of colors in the Kuiper belt, which was interpreted as evidence that Trans-Neptunian Objects (TNOs) formed from a locally homogeneous and globally heterogeneous protoplanetary disk. In this paradigm, binary pairs must have formed early, before the dynamically hot TNOs were scattered out from their formation location. The latter inferences, however, relied on the assumption that the matching colors of the binary components imply matching composition. Here, we test this assumption by examining the component-resolved photometry of three TNBs found in the Outer Solar System Origins Survey: 505447 (2013 SQ99), 511551 (2014 UD225) and 506121 (2016 BP81), across the visible and J-band near-infrared wavelength range. We report similar colors within 2 sigma for the binary pairs suggestive of similar reflectance spectra and hence surface composition. This advocates for gravitational collapse of pebble clouds as a possible TNO formation route. We however stress that several similarly small TNOs, including at least one binary, have been shown to exhibit substantial spectral variability in the near-infrared, implying color equality of binary pairs is likely to be violated in some cases.

We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0<z<1.4, drawn from the Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) survey. Based on more than 500 hours of Gemini/GMOS spectroscopy, and deep multi-band photometry taken with a range of observatories, we probe the SMFs down to a stellar mass limit of 10^9.7 Msun (10^9.5 Msun for star-forming galaxies). At this early epoch, the fraction of quiescent galaxies is already highly elevated in the clusters compared to the field at the same redshift. The quenched fraction excess (QFE) represents the fraction of galaxies that would be star-forming in the field, but are quenched due to their environment. The QFE is strongly mass dependent, and increases from ~30% at Mstar=10^9.7 Msun, to ~80% at Mstar=10^11.0 Msun. Nonetheless, the shapes of the SMFs of the two individual galaxy types, star-forming and quiescent galaxies, are identical between the clusters and the field - to high statistical precision. Yet, along with the different quiescent fractions is the total galaxy SMF environmentally dependent, with a relative deficit of low-mass galaxies in the clusters. These results are in stark contrast with findings in the local Universe, and thus require a substantially different quenching mode to operate at early times. We discuss these results in the light of several popular quenching models.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$\alpha)\leq50$\AA\ and we find no correlation or anti-correlation between EW (Ly$\alpha$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

We present a two-component Machine Learning (ML) based approach for classifying astronomical images by data-quality via an examination of sources detected in the images and image pixel values from representative sources within those images. The first component, which uses a clustering algorithm, creates a proper and small fraction of the image pixels to determine the quality of the observation. The representative images (and associated tables) are ~800 times smaller than the original images, significantly reducing the time required to train our algorithm. The useful information in the images is preserved, permitting them to be classified in different categories, but the required storage is reduced. The second component, which is a deep neural network model, classifies the representative images. Using ground-based telescope imaging data, we demonstrate that the method can be used to separate 'usable' images from those that present some problems for scientific projects -- such as images that were taken in sub-optimal conditions. This method uses two different data sets as input to a deep model and provides better performance than if we only used the images' pixel information. The method may be used in cases where large and complex data sets should be examined using deep models. Our automated classification approach achieves 97% agreement when compared to classification generated via manual image inspection. We compare our method with traditional results and show that the method improves the results by about 10%, and also presents more comprehensive outcomes.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada-French-Hawaii telescope (CFHT). We use a new set of data extracted from VESTIGE to study the impact of the hostile cluster environment on the star formation process down to the scale of HII regions (~ 50 pc). HII regions are identified and their parameters measured using the HIIphot code on a sample of 114 late-type galaxies spanning a wide range in morphological type (Sa-Sd, Im, BCD), stellar mass (10^6.5 <= M_star <= 10^11 Mo), and star formation activity (10^-3 <= SFR <= 10 Mo yr^-1). Owing to the exquisite average resolution of the VESTIGE data (0.65 arcsec), we detect 11302 HII regions with an Halpha luminosity L(Halpha) >= 10^37 erg s^-1. We show that the typical number of HII regions in gas-stripped objects is significantly lower than in healthy late-types of similar stellar mass. We also show that in these gas-stripped galaxies the number of HII regions significantly drops outside the effective radius, suggesting that the quenching process occurs outside-in, in agreement with other multifrequency observations. These new results consistently confirm that the main mechanism responsible for the decrease of the star formation activity observed in cluster galaxies is ram pressure, allowing us to discard other milder processes such as starvation or strangulation unable to reproduce the observed radially truncated profiles.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies.

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $\sigma_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $\alpha=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history.

We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]$<-1.2$, the algorithm succeeds in identifying more than 70% of the giants in the training/test set, with a dwarf contamination fraction below 30% (with respect to the SDSS/SEGUE dataset). The photometric metallicity estimates have uncertainties better than 0.2 dex when compared with the spectroscopic measurements. The distances estimated by the algorithm are valid out to a distance of at least $\sim 80$ kpc without requiring any prior on the stellar distribution, and have fully independent uncertainities that take into account both random and systematic errors. These advances allow us to estimate these stellar parameters for approximately 12 million stars in the photometric dataset. This will enable studies involving the chemical mapping of the distant outer disc and the stellar halo, including their kinematics using the Gaia proper motions. This type of algorithm can be applied in the Southern hemisphere to the first release of LSST data, thus providing an almost complete view of the external components of our Galaxy out to at least $\sim 80$ kpc. Critical to the success of these efforts will be ensuring well-defined spectroscopic training sets that sample a broad range of stellar parameters with minimal biases. A catalogue containing the training/test set and all relevant parameters within the public footprint of CFIS is available online.

The properties of low-mass galaxies hosting central black holes provide clues about the formation and evolution of the progenitors of supermassive black holes. In this letter, we present HSC-XD 52, a spectroscopically confirmed low-mass active galactic nucleus (AGN) at an intermediate redshift of $z\sim0.56$. We detect this object as a very luminous X-ray source coincident with a galaxy observed by the Hyper Suprime-Cam (HSC) as part of a broader search for low-mass AGN. We constrain its stellar mass through spectral energy distribution modeling to be LMC-like at $M_\star \approx 3 \times 10^9 M_\odot$, placing it in the dwarf regime. We estimate a central black hole mass of $M_\mathrm{BH} \sim 10^{6} M_\odot$. With an average X-ray luminosity of $L_X \approx 3.5 \times 10^{43}~\mathrm{erg}~\mathrm{s}^{-1}$, HSC-XD 52 is among the most luminous X-ray selected AGN in dwarf galaxies. The spectroscopic and photometric properties of HSC-XD 52 indicate that it is an intermediate redshift counterpart to local low-mass AGN.

The occurrence of planetary nebulae (PNe) in globular clusters (GCs) provides an excellent chance to study low-mass stellar evolution in a special (low-metallicity, high stellar density) environment. We report a systematic spectroscopic survey for the [O{\sc iii}] 5007 emission line of PNe in 1469 Virgo GCs and 121 Virgo ultra-compact dwarfs (UCDs), mainly hosted in the giant elliptical galaxies M87, M49, M86, and M84. We detected zero PNe in our UCD sample and discovered one PN ($M_{5007} = -4.1$ mag) associated with an M87 GC. We used the [O{\sc iii}] detection limit for each GC to estimate the luminosity-specific frequency of PNe, $\alpha$, and measured $\alpha$ in the Virgo cluster GCs to be $\alpha \sim 3.9_{-0.7}^{+5.2}\times 10^{-8}\mathrm{PN}/L_\odot$. $\alpha$ in Virgo GCs is among the lowest values reported in any environment, due in part to the large sample size, and is 5--6 times lower than that for the Galactic GCs. We suggest that $\alpha$ decreases towards brighter and more massive clusters, sharing a similar trend as the binary fraction, and the discrepancy between the Virgo and Galactic GCs can be explained by the observational bias in extragalactic surveys toward brighter GCs. This low but non-zero efficiency in forming PNe may highlight the important role played by binary interactions in forming PNe in GCs. We argue that a future survey of less massive Virgo GCs will be able to determine whether PN production in Virgo GCs is governed by internal process (mass, density, binary fraction), or is largely regulated by external environment.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

While collisional families are common in the asteroid belt, only one is known in the Kuiper belt, linked to the dwarf planet Haumea. The characterization of Haumea's family helps to constrain its origin and, more generally, the collisional history of the Kuiper belt. However, the size distribution of the Haumea family is difficult to constrain from the known sample, which is affected by discovery biases. Here, we use the Outer Solar System Origins Survey (OSSOS) Ensemble to look for Haumea family members. In this OSSOS XVI study we report the detection of three candidates with small ejection velocities relative to the family formation centre. The largest discovery, 2013 UQ15, is conclusively a Haumea family member, with a low ejection velocity and neutral surface colours. Although the OSSOS Ensemble is sensitive to Haumea family members to a limiting absolute magnitude (Hr) of 9.5 (inferred diameter of ~90 km), the smallest candidate is significantly larger, Hr = 7.9. The Haumea family members larger than approximately 20 km in diameter must be characterized by a shallow H-distribution slope in order to produce only these three large detections. This shallow size distribution suggests that the family formed in a graze-and-merge scenario, not a catastrophic collision.

We use published models of the early Solar System evolution with a slow, long-range and grainy migration of Neptune to predict the orbital element distributions and the number of modern-day Centaurs. The model distributions are biased by the OSSOS survey simulator and compared with the OSSOS Centaur detections. We find an excellent match to the observed orbital distribution, including the wide range of orbital inclinations which was the most troublesome characteristic to fit in previous models. A dynamical model, in which the original population of outer disk planetesimals was calibrated from Jupiter Trojans, is used to predict that OSSOS should detect 11+/-4 Centaurs with semimajor axis a<30 au, perihelion distance q>7.5 au and diameter D>10 km (absolute magnitude H_r<13.7 for a 6% albedo). This is consistent with 15 actual OSSOS Centaur detections with H_r<13.7. The population of Centaurs is estimated to be 21,000+/-8,000 for D>10 km. The inner scattered disk at 50<a<200 au should contain (2.0+/-0.8)x10^7 D>10 km bodies and the Oort cloud should contain (5.0+/-1.9)x10^8 D>10 km comets. Population estimates for different diameter cutoffs can be obtained from the size distribution of Jupiter Trojans (N(>D) proportional to D^(-2.1) for 5<D<100 km). We discuss model predictions for the Large Synoptic Survey Telescope observations of Centaurs.

The question of whether galaxy mergers are linked to the triggering of active galactic nuclei (AGN) continues to be a topic of considerable debate. The issue can be broken down into two distinct questions: 1) Can galaxy mergers trigger AGN? 2) Are galaxy mergers the dominant AGN triggering mechanism? A complete picture of the AGN-merger connection requires that both of these questions are addressed with the same dataset. In previous work, we have shown that galaxy mergers selected from the Sloan Digital Sky Survey (SDSS) show an excess of both optically-selected, and mid-IR colour-selected AGN, demonstrating that the answer to the first of the above questions is affirmative. Here, we use the same optical and mid-IR AGN selection to address the second question, by quantifying the frequency of morphological disturbances in low surface brightness r-band images from the Canada France Imaging Survey (CFIS). Only ~30 per cent of optical AGN host galaxies are morphologically disturbed, indicating that recent interactions are not the dominant trigger. However, almost 60 per cent of mid-IR AGN hosts show signs of visual disturbance, indicating that interactions play a more significant role in nuclear feeding. Both mid-IR and optically selected AGN have interacting fractions that are a factor of two greater than a mass and redshift matched non-AGN control sample, an excess that increases with both AGN luminosity and host galaxy stellar mass.

The migration of Neptune's resonances through the proto-Kuiper belt has been imprinted in the distribution of small bodies in the outer Solar System. Here we analyze five published Neptune migration models in detail, focusing on the high pericenter distance (high-q) trans-Neptunian Objects (TNOs) near Neptune's 5:2 and 3:1 mean-motion resonances, because they have large resonant populations, are outside the main classical belt, and are relatively isolated from other strong resonances. We compare the observationally biased output from these dynamical models with the detected TNOs from the Outer Solar System Origins Survey, via its Survey Simulator. All of the four new OSSOS detections of high-q non-resonant TNOs are on the Sunward side of the 5:2 and 3:1 resonances. We show that even after accounting for observation biases, this asymmetric distribution cannot be drawn from a uniform distribution of TNOs at 2sigma confidence. As shown by previous work, our analysis here tentatively confirms that the dynamical model that uses grainy slow Neptune migration provides the best match to the real high-q TNO orbital data. However, due to extreme observational biases, we have very few high-q TNO discoveries with which to statistically constrain the models. Thus, this analysis provides a framework for future comparison between the output from detailed, dynamically classified Neptune migration simulations and the TNO discoveries from future well-characterized surveys. We show that a deeper survey (to a limiting r-magnitude of 26.0) with a similar survey area to OSSOS could statistically distinguish between these five Neptune migration models.

The Colours of the Outer Solar System Origins Survey (Col-OSSOS) is acquiring near-simultaneous $g$, $r$, and $J$ photometry of unprecedented precision with the Gemini North Telescope, targeting nearly a hundred trans-Neptunian objects (TNOs) brighter than $m_r=23.6$ mag discovered in the Outer Solar System Origins Survey. Combining the optical and near-infrared photometry with the well-characterized detection efficiency of the Col-OSSOS target sample will provide the first flux-limited compositional dynamical map of the outer Solar System. In this paper, we describe our observing strategy and detail the data reduction processes we employ, including techniques to mitigate the impact of rotational variability. We present optical and near-infrared colors for 35 TNOs. We find two taxonomic groups for the dynamically excited TNOs, the neutral and red classes, which divide at $g-r \simeq 0.75$. Based on simple albedo and orbital distribution assumptions, we find that the neutral class outnumbers the red class, with a ratio of 4:1 and potentially as high as 11:1. Including in our analysis constraints from the cold classical objects, which are known to exhibit unique albedos and $r-z$ colors, we find that within our measurement uncertainty, our observations are consistent with the primordial Solar System protoplanetesimal disk being neutral-class-dominated, with two major compositional divisions in $grJ$ color space.

We present variability measurements and partial light curves of Trans-Neptunian Objects (TNOs) from a two-night pilot study using Hyper Suprime-Cam (HSC) on the Subaru Telescope (Maunakea, Hawai'i, USA). Subaru's large aperture (8-m) and HSC's large field of view (1.77 square degrees) allow us to obtain measurements of multiple objects with a range of magnitudes in each telescope pointing. We observed 65 objects with m_r = 22.6--25.5 mag in just six pointings, allowing 20--24 visits of each pointing over the two nights. Our sample, all discovered in the recent Outer Solar System Origins Survey (OSSOS), span absolute magnitudes H_r = 6.2--10.8 mag and thus investigates smaller objects than previous light curve projects have typically studied. Our data supports the existence of a correlation between light curve amplitude and absolute magnitude seen in other works, but does not support a correlation between amplitude and orbital inclination. Our sample includes a number of objects from different dynamical populations within the trans-Neptunian region, but we do not find any relationship between variability and dynamical class. We were only able to estimate periods for 12 objects in the sample and found that a longer baseline of observations is required for reliable period analysis. We find that 31 objects (just under half of our sample) have variability greater than 0.4 magnitudes during all of the observations; in smaller 1.25 hr, 1.85 hr and 2.45 hr windows, the median variability is 0.13, 0.16 and 0.19 mags, respectively. The fact that variability on this scale is common for small TNOs has important implications for discovery surveys (such as OSSOS or the Large Synoptic Survey Telescope) and color measurements.

We present the first results of our pilot study of 8 photometrically selected Lyman continuum (LyC) emitting galaxy candidates from the COSMOS field and focus on their optical emission line ratios. Observations were performed in the H and K bands using the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) instrument at the Keck Observatory, targeting the [OII], H$\beta$, and [OIII] emission lines. We find that photometrically selected LyC emitting galaxy candidates have high ionization parameters, based on their high [OIII]/[OII] ratios (O32), with an average ratio for our sample of 2.5$\pm$0.2. Preliminary results of our companion Low Resolution Imaging Spectrometer (LRIS) observations, targeting LyC and Ly$\alpha$, show that those galaxies with the largest O32 are typically found to also be Ly$\alpha$ emitters. High O32 galaxies are also found to have tentative non-zero LyC escape fractions ($f_{esc}(LyC)$) based on $u$ band photometric detections. These results are consistent with samples of highly ionized galaxies, including confirmed LyC emitting galaxies from the literature. We also perform a detailed comparison between the observed emission line ratios and simulated line ratios from density bounded H$_{\textrm{II}}$ regions modeled using the photoionization code MAPPINGS V. Estimates of $f_{esc}(LyC)$ for our sample fall in the range from 0.0-0.23 and suggest possible tension with published correlations between O32 and $f_{esc}(LyC)$, adding weight to dichotomy of arguments in the literature. We highlight the possible effects of clumpy geometry and mergers that may account for such tension.

Both physical and dynamical properties must be considered to constrain the origins of the dynamically excited distant Solar System populations. We present high-precision (g-r) colors for 25 small (Hr>5) dynamically excited Trans-Neptunian Objects (TNOs) and centaurs acquired as part of the Colours of the Outer Solar System Origins Survey (Col-OSSOS). We combine our dataset with previously published measurements and consider a set of 229 colors of outer Solar System objects on dynamically excited orbits. The overall color distribution is bimodal and can be decomposed into two distinct classes, termed `gray' and `red', that each has a normal color distribution. The two color classes have different inclination distributions: red objects have lower inclinations than the gray ones. This trend holds for all dynamically excited TNO populations. Even in the worst-case scenario, biases in the discovery surveys cannot account for this trend: it is intrinsic to the TNO population. Considering that TNOs are the precursors of centaurs, and that their inclinations are roughly preserved as they become centaurs, our finding solves the conundrum of centaurs being the only outer Solar System population identified so far to exhibit this property (Tegler et al. 2016). The different orbital distributions of the gray and red dynamically excited TNOs provide strong evidence that their colors are due to different formation locations in a disk of planetesimals with a compositional gradient.

We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and ~ 8 kpc to the south-east, with a couple of plumes of ionised gas, the weakest of which, at ~ 18 kpc from the nucleus, was previously unknown. The analysis of deep optical images taken from the NGVS survey confirms that this gas filament is associated with dust seen in absorption which is now detected up to ~ 2.4 kpc from the nucleus. We have also analysed the physical and kinematical properties of the ionised gas filament using deep IFU MUSE data covering the central 4.8 x 4.8 kpc^2 of the galaxy. The spectroscopic data confirms a perturbed kinematics of the ionised gas, with differences in velocity of ~ 700-800 km s^-1 on scales of < 1 kpc. The analysis of 2D diagnostic diagrams and the observed relationship between the shock-sensitive [OI]/Halpha line ratio and the velocity dispersion of the gas suggest that the gas is shock-ionised.

Using deep, high resolution optical imaging from the Next Generation Virgo Cluster Survey we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses $10^{5}\lesssim M_{*}/M_{\odot} \lesssim10^{12}$. The nucleation fraction reaches a peak value $f_{n}\approx90\%$ for $M_{*} \approx 10^{9} M_{\odot}$ galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as $M_{*} \approx5\times10^{5} M_{\odot}$. Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at $M_{NSC}\approx7\times10^{5} M_{\odot}$, a factor 3-4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to $M_{NSC}/M_{*}\approx3.6\times10^{-3}$ for $M_{*} \approx 5\times10^{9} M_{\odot}$ galaxies. Nuclei in both more and less massive galaxies are much more prominent: $M_{NSC}\propto M_{*}^{0.46}$ at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including a very similar occupation distribution and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

We have observed the late-type peculiar galaxy NGC 4424 during VESTIGE, a blind narrow-band Halpha[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT. The presence of a 110 kpc long HI tail in the S direction indicates that this galaxy is undergoing a ram pressure stripping event. The deep narrow-band image revealed the presence of a low surface brightness ionised gas tail ~10 kpc long extending from the centre of the galaxy to the NW direction, in the direction opposite to the HI tail. Chandra and XMM X-rays data do not show any compact source in the nucleus nor the presence of an extended tail of hot gas, while MUSE spectroscopy indicates that the gas is photo-ionised in the inner regions and shock-ionised in the outer parts. IFU spectroscopy confirms that the ionised gas is kinematically decoupled from the stellar component and indicates the presence of two kinematically distinct structures in the stellar disc. The analysis of the SED of the galaxy indicates that the activity of star formation has been totally quenched in the outer disc ~ 250-280 Myr ago, while only reduced by ~80% in the central regions. All this observational evidence suggests that NGC 4424 is the remnant of an unequal-mass merger occurred <= 500 Myr ago, when the galaxy was already a member of the Virgo cluster, now undergoing a ram pressure stripping event which has removed the gas and quenched the activity of star formation in the outer disc. The tail of ionised gas probably results from the outflow produced by a central starburst fed by the collapse of gas induced by the merging episode. This outflow is sufficiently powerful to overcome the ram pressure induced by the intracluster medium on the disc of the galaxy crossing the cluster. This analysis thus suggests that feedback can participate in the quenching process of galaxies in high-density regions.

Context. Centaurs are icy objects in transition between the transneptunian region and the inner solar system, orbiting the Sun in the giant planet region. Some Centaurs display cometary activity, which cannot be sustained by the sublimation of water ice in this part of the solar system, and has been hypothesized to be due to the crystallization of amorphous water ice. Aims. In this work, we look at Centaurs discovered by the Outer Solar System Origins Survey (OSSOS) and search for cometary activity. Tentative detections would improve understanding of the origins of activity among these objects. Methods. We search for comae and structures by fitting and subtracting both Point Spread Functions (PSF) and Trailed point-Spread Functions (TSF) from the OSSOS images of each Centaur. When available, Col-OSSOS images were used to search for comae too. Results. No cometary activity is detected in the OSSOS sample. We track the recent orbital evolution of each new Centaur to confirm that none would actually be predicted to be active, and we provide size estimates for the objects. Conclusions. The addition of 20 OSSOS objects to the population of 250 known Centaurs is consistent with the currently understood scenario, in which drastic drops in perihelion distance induce changes in the thermal balance prone to trigger cometary activity in the giant planet region.

We present the stellar density profile of the outer halo of the Galaxy traced over a range of Galactocentric radii from $15< R_{GC} < 220$ kpc by blue horizontal branch (BHB) stars. These stars are identified photometrically using deep $u-$band imaging from the new Canada-France-Imaging-Survey (CFIS) that reaches 24.5 mag. This is combined with $griz$ bands from Pan-STARRS 1 and covers a total of $\sim4000$ deg$^2$ of the northern sky. We present a new method to select BHB stars that has low contamination from blue stragglers and high completeness. We use this sample to measure and parameterize the three dimensional density profile of the outer stellar halo. We fit the profile using (i) a simple power-law with a constant flattening (ii) a flattening that varies as a function of Galactocentric radius (iii) a broken power law profile. We find that outer stellar halo traced by the BHB is well modelled by a broken power law with a constant flattening of $q=0.86 \pm 0.02$, with an inner slope of $\gamma=4.24 \pm 0.08$. This is much steeper than the preferred outer profile that has a slope of $\beta=3.21\pm 0.07$ after a break radius of $r_b=41.4^{+2.5}_{-2.4}$ kpc. The outer profile of the stellar halo trace by BHB stars is shallower than that recently measured using RR Lyrae, a surprising result given the broad similarity of the ages of these stellar populations.

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$.

A unique signature for the presence of massive black holes in very dense stellar regions is occasional giant-amplitude outbursts of multiwavelength radiation from tidal disruption and subsequent accretion of stars that make a close approach to the black holes. Previous strong tidal disruption event (TDE) candidates were all associated with the centers of largely isolated galaxies. Here we report the discovery of a luminous X-ray outburst from a massive star cluster at a projected distance of 12.5 kpc from the center of a large lenticular galaxy. The luminosity peaked at ~10^{43} erg/s and decayed systematically over 10 years, approximately following a trend that supports the identification of the event as a TDE. The X-ray spectra were all very soft, with emission confined to be <3.0 keV, and could be described with a standard thermal disk. The disk cooled significantly as the luminosity decreased, a key thermal-state signature often observed in accreting stellar-mass black holes. This thermal-state signature, coupled with very high luminosities, ultrasoft X-ray spectra and the characteristic power-law evolution of the light curve, provides strong evidence that the source contains an intermediate-mass black hole (IMBH) with a mass of a few ten thousand solar mass. This event demonstrates that one of the most effective means to detect IMBHs is through X-ray flares from TDEs in star clusters.

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System.

We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100~au. The detection of these objects implies a 9:1 resonance population of $1.1\times10^4$ objects with $H_r<8.66$ ($D~\gtrsim~100$~km) on similar orbits (95\% confidence range of $\sim0.4-3\times10^4$). Integrations over 4~Gyr of an ensemble of clones spanning these objects' orbit fit uncertainties reveal that they both have median resonance occupation timescales of $\sim1$~Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last $\sim1$~Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000-4500 objects with $H_r<8.66$; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper belt's dynamical history.

NASA's New Horizons spacecraft will conduct a close flyby of the cold classical Kuiper Belt Object (KBO) designated (486958) 2014 MU69 on January 1, 2019. At a heliocentric distance of 44 AU, "MU69" will be the most distant object ever visited by a spacecraft. To enable this flyby, we have developed an extremely high precision orbit fitting and uncertainty processing pipeline, making maximal use of the Hubble Space Telescope's Wide Field Camera 3 (WFC3) and pre-release versions of the ESA Gaia Data Release 2 (DR2) catalog. This pipeline also enabled successful predictions of a stellar occultation by MU69 in July 2017. We describe how we process the WFC3 images to match the Gaia DR2 catalog, extract positional uncertainties for this extremely faint target (typically 140 photons per WFC3 exposure), and translate those uncertainties into probability distribution functions for MU69 at any given time. We also describe how we use these uncertainties to guide New Horizons, plan stellar occultions of MU69, and derive MU69's orbital evolution and long-term stability.

We present Keck/DEIMOS spectroscopy of globular clusters (GCs) around the ultra-diffuse galaxies (UDGs) VLSB-B, VLSB-D, and VCC615 located in the central regions of the Virgo cluster. We spectroscopically identify 4, 12, and 7 GC satellites of these UDGs, respectively. We find that the three UDGs have systemic velocities ($V_{sys}$) consistent with being in the Virgo cluster, and that they span a wide range of velocity dispersions, from $\sim 16$ to $\sim 47$ km/s, and high dynamical mass-to-light ratios within the radius that contains half the number of GCs ($ 407^{+916}_{-407}$, $21^{+15}_{-11}$, $60^{+65}_{-38}$, respectively). VLSB-D shows possible evidence for rotation along the stellar major axis and its $V_{sys}$ is consistent with that of the massive galaxy M84 and the center of the Virgo cluster itself. These findings, in addition to having a dynamically and spatially ($\sim 1$ kpc) off-centered nucleus and being extremely elongated, suggest that VLSB-D could be tidally perturbed. On the contrary, VLSB-B and VCC615 show no signals of tidal deformation. Whereas the dynamics of VLSB-D suggest that it has a less massive dark matter halo than expected for its stellar mass, VLSB-B and VCC615 are consistent with a $\sim 10^{12}$ M$_{\odot}$ dark matter halo. Although our samples of galaxies and GCs are small, these results suggest that UDGs may be a diverse population, with their low surface brightnesses being the result of very early formation, tidal disruption, or a combination of the two.

The scattering trans-Neptunian Objects (TNOs) can be measured to smaller sizes than any other distant small-body population. We use the largest sample yet obtained, 68 discoveries, primarily by the Outer Solar System Origins Survey (OSSOS), to constrain the slope of its luminosity distribution, with sensitivity to much fainter absolute $H$ magnitudes than previous work. Using the analysis technique in Shankman et al. (2016), we confirm that a single slope for the $H$-distribution is not an accurate representation of the scattering TNOs and Centaurs, and that a break in the distribution is required, in support of previous conclusions. A bright-end slope of $\alpha_b=0.9$ transitioning to a faint-end slope $\alpha_f$ of 0.4-0.5 with a differential number contrast $c$ from 1 (a knee) to 10 (a divot) provides an acceptable match to our data. We find that break magnitudes $H_b$ of 7.7 and 8.3, values both previously suggested for dynamically hot Kuiper belt populations, are equally non-rejectable for a range of $\alpha_f$ and $c$ in our statistical analysis. Our preferred divot $H$-distribution transitions to $\alpha_f=0.5$ with a divot of contrast $c=3$ at $H_b=8.3$, while our preferred knee $H$-distribution transitions to $\alpha_f=0.4$ at $H_b=7.7$. The intrinsic population of scattering TNOs required to match the OSSOS detections is $3\times10^6$ for $H_r<12$, and $9\times10^4$ for $H_r<8.66$ ($D\gtrsim100$~km), with Centaurs having an intrinsic population two orders of magnitude smaller.

We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the prospector fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). We find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. In the future, our approach can be applied to a large and diverse sample of UDGs down to faint surface brightness limits, with the goal of constraining their stellar ages, stellar metallicities, and circumstellar and diffuse interstellar dust content.

During pilot observations of the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT, we have observed the spiral galaxy NGC 4254 (M99). Deep Halpha+[NII] narrow-band and GALEX UV images revealed the presence of 60 compact (70-500 pc radius) star forming regions up to ~ 20 kpc outside the optical disc of the galaxy. These regions are located along a tail of HI gas stripped from the disc of the galaxy after a rapid gravitational encounter with another Virgo cluster member that simulations indicate occurred 280-750 Myr ago. We have combined the VESTIGE data with multifrequency data from the UV to the far-infrared to characterise the stellar populations of these regions and study the star formation process in an extreme environment such as the tails of stripped gas embedded in the hot intracluster medium. The colour, spectral energy distribution (SED), and linear size consistently indicate that these regions are coeval and have been formed after a single burst of star formation that occurred ~< 100 Myr ago. These regions might become free floating objects within the cluster potential well, and be the local analogues of compact sources produced after the interaction of gas-rich systems that occurred during the early formation of clusters.

Substructure in globular cluster (GC) populations around large galaxies is expected in galaxy formation scenarios that involve accretion or merger events, and it has been searched for using direct associations between GCs and structure in the diffuse galaxy light, or with GC kinematics. Here, we present a search for candidate substructures in the GC population around the Virgo cD galaxy M87 through the analysis of the spatial distribution of the GC colors.~The study is based on a sample of $\sim\!1800$ bright GCs with high-quality $u,g,r,i,z,K_s$ photometry, selected to ensure a low contamination by foreground stars or background galaxies.~The spectral energy distributions of the GCs are associated with formal estimates of age and metallicity, which are representative of its position in a 4-D color-space relative to standard single stellar population models.~Dividing the sample into broad bins based on the relative formal ages, we observe inhomogeneities which reveal signatures of GC substructures.~The most significant of these is a spatial overdensity of GCs with relatively young age labels, of diameter $\sim\!0.1$\,deg ($\sim\!30\,$kpc), located to the south of M87.~The significance of this detection is larger than about 5$\sigma$ after accounting for estimates of random and systematic errors.~Surprisingly, no large Virgo galaxy is present in this area, that could potentially host these GCs.~But candidate substructures in the M87 halo with equally elusive hosts have been described based on kinematic studies in the past.~The number of GC spectra available around M87 is currently insufficient to clarify the nature of the new candidate substructure.

We describe a program to measure surface brightness fluctuation (SBF) distances to galaxies observed in the Next Generation Virgo Cluster Survey (NGVS), a photometric imaging survey covering $104~deg^2$ of the Virgo cluster in the ${u}^*,g,i,z$ bandpasses with the Canada-France Hawaii Telescope. We describe the selection of the sample galaxies, the procedures for measuring the apparent $i$-band SBF magnitude $\bar{i}$, and the calibration of the absolute $\bar{M}_i$ as a function of observed stellar population properties. The multi-band NGVS data set provides multiple options for calibrating the SBF distances, and we explore various calibrations involving individual color indices as well as combinations of two different colors. Within the color range of the present sample, the two-color calibrations do not significantly improve the scatter with respect to wide-baseline, single-color calibrations involving $u^{*}$. We adopt the ${u}^*{-}z$ calibration as reference for the present galaxy sample, with an observed scatter of 0.11 mag. For a few cases that lack good ${u}^*$ photometry, we use an alternative relation based on a combination of $g{-}i$ and $g{-}z$ colors, with only a slightly larger observed scatter of 0.12 mag. The agreement of our measurements with the best existing distance estimates provides confidence that our measurements are accurate. We present a preliminary catalog of distances for 89 galaxies brighter than $B_T\approx13.0$ mag within the survey footprint, including members of the background M and W Clouds at roughly twice the distance of the main body of the Virgo cluster. The extension of the present work to fainter and bluer galaxies is in progress.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band H$\alpha$+[NII] imaging survey carried out with MegaCam at the CFHT. During pilot observations we have observed NGC 4330, an intermediate mass, edge-on star forming spiral currently falling into the core of the Virgo cluster. New deep observations revealed a low surface brightness 10 kpc tail exhibiting a peculiar filamentary structure. The filaments are remarkably parallel one another and clearly indicate the direction of motion of the galaxy in the Virgo potential. Motivated by the detection of these features, indicating ongoing gas stripping, we collected literature photometry in 15 bands from the far-UV to the far-IR and deep optical long slit spectroscopy using the FORS2 instrument at the ESO Very Large Telescope. Using a newly developed Monte Carlo code that jointly fits spectroscopy and photometry, we reconstructed the star formation histories in apertures along the major axis of the galaxy. Our results have been validated against the output of CIGALE, a fitting code which has been previously used for similar studies. We found a clear outside-in gradient with radius of the time when the quenching event started: the outermost radii have been stripped 500 Myr ago, while the stripping has reached the inner 5 kpc from the center in the last 100 Myr. Regions at even smaller radii are currently still forming stars fueled by the presence of HI and H2 gas. When compared to statistical studies of the quenching timescales in the local Universe we find that ram pressure stripping of the cold gas is an effective mechanism to reduce the transformation times for galaxies falling into massive clusters. Future systematic studies of all the active galaxies observed by VESTIGE in the Virgo cluster will extend these results to a robust statistical framework.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.

We revisit the proposed extended Schmidt law (Shi et al. 2011) which points that the star formation efficiency in galaxies depends on the stellar mass surface density, by investigating spatially-resolved star formation rates (SFRs), gas masses and stellar masses of star formation regions in a vast range of galactic environments, from the outer disks of dwarf galaxies to spiral disks and to merging galaxies as well as individual molecular clouds in M33. We find that these regions are distributed in a tight power-law as Sigma_SFR ~(Sigma_star^0.5 Sigma_gas )^1.09, which is also valid for the integrated measurements of disk and merging galaxies at high-z. Interestingly, we show that star formation regions in the outer disks of dwarf galaxies with Sigma_SFR down to 10^(-5) Msun/yr/kpc^2, which are outliers of both Kennicutt-Schmidt and Silk-Elmegreen law, also follow the extended Schmidt law. Other outliers in the Kennicutt-Schmidt law, such as extremely-metal poor star-formation regions, also show significantly reduced deviations from the extended Schmidt law. These results suggest an important role for existing stars in helping to regulate star formation through the effect of their gravity on the mid-plane pressure in a wide range of galactic environments.

We measured stacked weak lensing cluster masses for a sample of 1325 galaxy clusters detected by the RedGOLD algorithm in the Canada-France-Hawaii Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at $0.2<z<0.5$, in the optical richness range $10<\lambda<70$. After a selection of our best richness subsample ($20<\lambda<50$), this is the most comprehensive lensing study of a $\sim 100\%$ complete and $\sim 90\%$ pure optical cluster catalogue in this redshift range, with a total of 346 clusters in $\sim164~deg^2$. We test three different mass models, and our best model includes a basic halo model, with a Navarro Frenk and White profile, and correction terms that take into account cluster miscentering, non-weak shear, the two-halo term, the contribution of the Brightest Cluster Galaxy, and an a posteriori correction for the intrinsic scatter in the mass-richness relation. With this model, we obtain a mass-richness relation of $\log{M_{\rm 200}/M_{\odot}}=(14.48\pm0.04)+(1.14\pm0.23)\log{(\lambda/40)}$ (statistical uncertainties). This result is consistent with other published lensing mass-richness relations. When compared to X-ray masses and mass proxies, we find that on average weak lensing masses are $\sim 10\%$ higher than those derived in the X-ray in the range $2\times10^{13}M_{\rm \odot}<E(z) M^{X}_{\rm 200}<2\times10^{14}M_{\rm \odot}$, in agreement with most previous results and simulations. We also give the coefficients of the scaling relations between the lensing mass and X-ray mass proxies, $L_X$ and $T_X$, and compare them with previous results.

We present an analysis of 39 nuclei and their early-type hosts in the Virgo Cluster using ten broadband filters: F300W, F475W, F850LP, F160W, $u^*griz$, and $K_s$. We describe the Virgo Redux program, which provides high-resolution UV and NIR imaging. Combining this data with optical and NIR imaging from the ACS Virgo Cluster Survey and the Next Generation Virgo Cluster Survey, we estimate masses, metallicities and ages using simple stellar population (SSP) models. For 19 nuclei, we compare to SSP parameters derived from Keck and Gemini spectra and find reasonable agreement between the photometric and spectroscopic metallicity: the RMS scatter is 0.3 dex. We reproduce the nucleus-galaxy mass fraction of $0.33^{+0.09}_{-0.07}$ percent for galaxy stellar masses $10^{8.4}-10^{10.3} M_\odot$ with a typical precision of $\sim$35% for the nuclei masses. Based on available model predictions, there is no single preferred formation scenario for nuclei, suggesting that nuclei are formed stochastically through a mix of processes. Nuclei metallicities are statistically identical to those of their hosts, appearing $0.07 \pm 0.3$ dex more metal-rich on average -- although, omitting galaxies with unusual origins, nuclei are $0.20\pm0.28$ dex more metal-rich. Nuclei appear to be $0.56 \pm 0.12$ dex more metal rich than ultra-compact dwarf galaxies (UCDs) at fixed mass. We find no clear age difference between nuclei and their galaxies, with nuclei displaying a broad range of ages. Interestingly, we find that the most massive nuclei may be flatter and more closely aligned with the semi-major axes of their hosts, suggesting that they formed through predominantly dissipative processes.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

Several different classes of trans-Neptunian objects (TNOs) have been identified based on their optical and near-infrared colors. As part of the Colours of the Outer Solar System Origins Survey, we have obtained $g$, $r$, and $z$ band photometry of 26 TNOs using Subaru and Gemini Observatories. Previous color surveys have not utilized $z$ band reflectance, and the inclusion of this band reveals significant surface reflectance variations between sub-populations. The colors of TNOs in $g-r$ and $r-z$ show obvious structure, and appear consistent with the previously measured bi-modality in $g-r$. The distribution of colors of the two dynamically excited surface types can be modeled using the two-component mixing models from Fraser \& Brown (2012). With the combination of $g-r$ and $r-z$, the dynamically excited classes can be separated cleanly into red and neutral surface classes. In $g - r$ and $r - z$, the two dynamically excited surface groups are also clearly distinct from the cold classical TNO surfaces, which are red, with $g-r\gtrsim$0.85 and $r-z\lesssim$0.6, while all dynamically excited objects with similar $g-r$ colors exhibit redder $r-z$ colors. The $z$ band photometry makes it possible for the first time to differentiate the red excited TNO surfaces from the red cold classical TNO surfaces. The discovery of different $r-z$ colors for these cold classical TNOs makes it possible to search for cold classical surfaces in other regions of the Kuiper belt and to completely separate cold classical TNOs from the dynamically excited population, which overlaps in orbital parameter space.

We use three different techniques to identify hundreds of white dwarf (WD) candidates in the Next Generation Virgo Cluster Survey (NGVS) based on photometry from the NGVS and GUViCS, and proper motions derived from the NGVS and the Sloan Digital Sky Survey (SDSS). Photometric distances for these candidates are calculated using theoretical color-absolute magnitude relations while effective temperatures are measured by fitting their spectral energy distributions. Disk and halo WD candidates are separated using a tangential velocity cut of 200 km~s$^{-1}$ in a reduced proper motion diagram, which leads to a sample of six halo WD candidates. Cooling ages, calculated for an assumed WD mass of 0.6$M_{\odot}$, range between 60 Myr and 6 Gyr, although these estimates depend sensitively on the adopted mass. Luminosity functions for the disk and halo subsamples are constructed and compared to previous results from the SDSS and SuperCOSMOS survey. We compute a number density of (2.81 $\pm$ 0.52) $\times 10^{-3}$~pc$^{-3}$ for the disk WD population--- consistent with previous measurements. We find (7.85 $\pm$ 4.55) $\times 10^{-6}$~pc$^{-3}$ for the halo, or 0.3\% of the disk. Observed stellar counts are also compared to predictions made by the TRILEGAL and Besan\c{c}on stellar population synthesis models. The comparison suggests that the TRILEGAL model overpredicts the total number of WDs. The WD counts predicted by the Besan\c{c}on model agree with the observations, although a discrepancy arises when comparing the predicted and observed halo WD populations; the difference is likely due to the WD masses in the adopted model halo.

We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.

The accumulating, but small, set of large semi-major axis trans-Neptunian objects (TNOs) shows an apparent clustering in the orientations of their orbits. This clustering must either be representative of the intrinsic distribution of these TNOs, or else arise as a result of observation biases and/or statistically expected variations for such a small set of detected objects. The clustered TNOs were detected across different and independent surveys, which has led to claims that the detections are therefore free of observational bias. This apparent clustering has led to the so-called "Planet 9" hypothesis that a super-Earth currently resides in the distant solar system and causes this clustering. The Outer Solar System Origins Survey (OSSOS) is a large program that ran on the Canada-France-Hawaii Telescope from 2013--2017, discovering more than 800 new TNOs. One of the primary design goals of OSSOS was the careful determination of observational biases that would manifest within the detected sample. We demonstrate the striking and non-intuitive biases that exist for the detection of TNOs with large semi-major axes. The eight large semi-major axis OSSOS detections are an independent dataset, of comparable size to the conglomerate samples used in previous studies. We conclude that the orbital distribution of the OSSOS sample is consistent with being detected from a uniform underlying angular distribution.

Large samples of globular clusters (GC) with precise multi-wavelength photometry are becoming increasingly available and can be used to constrain the formation history of galaxies. We present the results of an analysis of Milky Way (MW) and Virgo core GCs based on five optical-near-infrared colors and ten synthetic stellar population models. For the MW GCs, the models tend to agree on photometric ages and metallicities, with values similar to those obtained with previous studies. When used with Virgo core GCs, for which photometry is provided by the Next Generation Virgo cluster Survey (NGVS), the same models generically return younger ages. This is a consequence of the systematic differences observed between the locus occupied by Virgo core GCs and models in panchromatic color space. Only extreme fine-tuning of the adjustable parameters available to us can make the majority of the best-fit ages old. Although we cannot exclude that the formation history of the Virgo core may lead to more conspicuous populations of relatively young GCs than in other environments, we emphasize that the intrinsic properties of the Virgo GCs are likely to differ systematically from those assumed in the models. Thus, the large wavelength coverage and photometric quality of modern GC samples, such as used here, is not by itself sufficient to better constrain the GC formation histories. Models matching the environment-dependent characteristics of GCs in multi-dimensional color space are needed to improve the situation.

We present the Pristine survey, a new narrow-band photometric survey focused on the metallicity-sensitive Ca H & K lines and conducted in the northern hemisphere with the wide-field imager MegaCam on the Canada-France-Hawaii Telescope (CFHT). This paper reviews our overall survey strategy and discusses the data processing and metallicity calibration. Additionally we review the application of these data to the main aims of the survey, which are to gather a large sample of the most metal-poor stars in the Galaxy, to further characterise the faintest Milky Way satellites, and to map the (metal-poor) substructure in the Galactic halo. The current Pristine footprint comprises over 1,000 deg2 in the Galactic halo ranging from b~30 to 78 and covers many known stellar substructures. We demonstrate that, for SDSS stellar objects, we can calibrate the photometry at the 0.02-magnitude level. The comparison with existing spectroscopic metallicities from SDSS/SEGUE and LAMOST shows that, when combined with SDSS broad-band g and i photometry, we can use the CaHK photometry to infer photometric metallicities with an accuracy of ~0.2 dex from [Fe/H]=-0.5 down to the extremely metal-poor regime ([Fe/H]<-3.0). After the removal of various contaminants, we can efficiently select metal-poor stars and build a very complete sample with high purity. The success rate of uncovering [Fe/H]SEGUE<-3.0 stars among [Fe/H]Pristine<-3.0 selected stars is 24% and 85% of the remaining candidates are still very metal poor ([Fe/H]<-2.0). We further demonstrate that Pristine is well suited to identify the very rare and pristine Galactic stars with [Fe/H]<-4.0, which can teach us valuable lessons about the early Universe.

The cold classical Kuiper belt objects have low inclinations and eccentricities and are the only Kuiper belt population suspected to have formed in situ. Compared with the dynamically excited populations, which exhibit a broad range of colours and a low binary fraction of ~10% cold classical Kuiper belt objects typically have red optical colours with ~30% of the population found in binary pairs; the origin of these differences remains unclear. We report the detection of a population of blue-coloured, tenuously bound binaries residing among the cold classical Kuiper belt objects. Here we show that widely separated binaries could have survived push-out into the cold classical region during the early phases of Neptune's migration. The blue binaries may be contaminants, originating at ~38 au, and could provide a unique probe of the formative conditions in a region now nearly devoid of objects. The idea that the blue objects, which are predominantly binary, are the products of push-out requires that the planetesimals formed entirely as multiples. Plausible formation routes include planetesimal formation via pebble accretion and subsequent binary production through dynamic friction and binary formation during the collapse of a cloud of solids.

We report the discovery of the minor planet 2013 SY$_{99}$, on an exceptionally distant, highly eccentric orbit. With a perihelion of 50.0 au, 2013 SY$_{99}$'s orbit has a semi-major axis of $730 \pm 40$ au, the largest known for a high-perihelion trans-Neptunian object (TNO), well beyond those of (90377) Sedna and 2012 VP$_{113}$. Yet, with an aphelion of $1420 \pm 90$ au, 2013 SY$_{99}$'s orbit is interior to the region influenced by Galactic tides. Such TNOs are not thought to be produced in the current known planetary architecture of the Solar System, and they have informed the recent debate on the existence of a distant giant planet. Photometry from the Canada-France-Hawaii Telescope, Gemini North and Subaru indicate 2013 SY$_{99}$ is $\sim 250$ km in diameter and moderately red in colour, similar to other dynamically excited TNOs. Our dynamical simulations show that Neptune's weak influence during 2013 SY$_{99}$'s perihelia encounters drives diffusion in its semi-major axis of hundreds of astronomical units over 4 Gyr. The overall symmetry of random walks in semi-major axis allow diffusion to populate 2013 SY$_{99}$'s orbital parameter space from the 1000-2000 au inner fringe of the Oort cloud. Diffusion affects other known TNOs on orbits with perihelia of 45 to 49 au and semi-major axes beyond 250 au, providing a formation mechanism that implies an extended population, gently cycling into and returning from the inner fringe of the Oort cloud.

We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process.

We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.

Multiwavelength flares from tidal disruption and accretion of stars can be used to find and study otherwise dormant massive black holes in galactic nuclei. Previous well-monitored candidate flares are short-lived, with most emission confined to within ~1 year. Here we report the discovery of a well observed super-long (>11 years) luminous soft X-ray flare from the nuclear region of a dwarf starburst galaxy. After an apparently fast rise within ~4 months a decade ago, the X-ray luminosity, though showing a weak trend of decay, has been persistently high at around the Eddington limit (when the radiation pressure balances the gravitational force). The X-ray spectra are generally soft (steeply declining towards higher energies) and can be described with Comptonized emission from an optically thick low-temperature corona, a super-Eddington accretion signature often observed in accreting stellar-mass black holes. Dramatic spectral softening was also caught in one recent observation, implying either a temporary transition from the super-Eddington accretion state to the standard thermal state or the presence of a transient highly blueshifted (~0.36c) warm absorber. All these properties in concert suggest a tidal disruption event of an unusually long super-Eddington accretion phase that has never been observed before.

We aim to explore the relationship between globular cluster total number, $N_{\rm GC}$, and central black hole mass, $M_\bullet$, in spiral galaxies, and compare it with that recently reported for ellipticals. We present results for the Sbc galaxy NGC 4258, from Canada France Hawaii Telescope data. Thanks to water masers with Keplerian rotation in a circumnuclear disk, NGC 4258 has the most precisely measured extragalactic distance and supermassive black hole mass to date. The globular cluster (GC) candidate selection is based on the ($u^*\ -\ i^\prime$) vs. ($i^\prime\ -\ K_s$) diagram, which is a superb tool to distinguish GCs from foreground stars, background galaxies, and young stellar clusters, and hence can provide the best number counts of GCs from photometry alone, virtually free of contamination, even if the galaxy is not completely edge-on. The mean optical and optical-near infrared colors of the clusters are consistent with those of the Milky Way and M 31, after extinction is taken into account. We directly identify 39 GC candidates; after completeness correction, GC luminosity function extrapolation and correction for spatial coverage, we calculate a total $N_{\rm GC} = 144\pm31^{+38}_{-36}$ (random and systematic uncertainties, respectively). We have thus increased to 6 the sample of spiral galaxies with measurements of both $M_\bullet$ and $N_{\rm GC}$. NGC 4258 has a specific frequency $S_{\rm N} = 0.4\pm0.1$ (random uncertainty), and is consistent within 2$\sigma$ with the $N_{\rm GC}$ vs. $M_\bullet$ correlation followed by elliptical galaxies. The Milky Way continues to be the only spiral that deviates significantly from the relation.

Deep imaging of the diffuse light emitted by the stellar fine structures and outer halos around galaxies is now often used to probe their past mass assembly. Because the extended halos survive longer than the relatively fragile tidal features, they trace more ancient mergers. We use images reaching surface brightness limits as low as 28.5-29 mag.arcsec-2 (g-band) to obtain light and color profiles up to 5-10 effective radii of a sample of nearby early-type galaxies. They were acquired with MegaCam as part of the CFHT MATLAS large programme. These profiles may be compared to those produced by simulations of galaxy formation and evolution, once corrected for instrumental effects. Indeed they can be heavily contaminated by the scattered light caused by internal reflections within the instrument. In particular, the nucleus of galaxies generates artificial flux in the outer halo, which has to be precisely subtracted. We present a deconvolution technique to remove the artificial halos that makes use of very large kernels. The technique based on PyOperators is more time efficient than the model-convolution methods also used for that purpose. This is especially the case for galaxies with complex structures that are hard to model. Having a good knowledge of the Point Spread Function (PSF), including its outer wings, is critical for the method. A database of MegaCam PSF models corresponding to different seeing conditions and bands was generated directly from the deep images. It is shown that the difference in the PSFs in different bands causes artificial changes in the color profiles, in particular a reddening of the outskirts of galaxies having a bright nucleus. The method is validated with a set of simulated images and applied to three representative test cases: NGC 3599, NGC 3489, and NGC 4274, and exhibiting for two of them a prominent ghost halo. The method successfully removes it.

Globular clusters (GCs) are some of the most visible tracers of the merging and accretion history of galaxy halos. Metal-poor GCs, in particular, are thought to arrive in massive galaxies largely through dry, minor merging events, but it is rare to see a direct connection between GCs and visible stellar streams. NGC 474 is a post-merger early-type galaxy with dramatic fine structures made of concentric shells and radial streams that have been more clearly revealed by deep imaging. We present a study of GCs in NGC 474 to better establish the relationship between merger-induced fine structure and the GC system. We find that many GCs are superimposed on visible streams and shells, and about 35% of GCs outside $3R_{\rm e,galaxy}$ are located in regions of fine structure. The spatial correlation between the GCs and fine structure is significant at the 99.9% level, showing that this correlation is not coincidental. The colors of the GCs on the fine structures are mostly blue, and we also find an intermediate-color population that is dominant in the central region, and which will likely passively evolve to have colors consistent with a traditional metal-rich GC population. The association of the blue GCs with fine structures is direct confirmation that many metal-poor GCs are accreted onto massive galaxy halos through merging events, and that progenitors of these mergers are sub-L* galaxies.

The Next Generation Virgo Cluster Survey is a deep (with a $2\sigma$ detection limit $\mu_g$ = 29~mag~arcsec$^{-2}$ in the $g-$band) optical panchromatic survey targeting the Virgo cluster from its core to virial radius, for a total areal coverage of 104 square degrees. As such, the survey is well suited for the study of galaxies' outskirts, haloes and low surface brightness features that arise from dynamical interactions within the cluster environment. We report the discovery of extremely faint ($\mu_g$ $>$ 25 mag arcsec$^{-2}$) shells in three Virgo cluster early-type dwarf galaxies, VCC~1361, VCC~1447 and VCC~1668. Among them, VCC~1447 has an absolute magnitude M$_{g}$ = -11.71 mag and is {\it the least massive galaxy with a shell system discovered to date}. We present a detailed study of these low surface brightness features. We detect between three and four shells in each of our galaxies. Within the uncertainties, we find no evidence of a color difference between the galaxy main body and shell features. The observed arcs of the shells are located upto several effective radii of the galaxies. We further explore the origin of these low surface brightness features with the help of idealized numerical simulations. We find that a near equal mass merger is best able to reproduce the main properties of the shells, including their quite symmetric appearance and their alignment along the major axis of the galaxy. The simulations provide support for a formation scenario in which a recent merger, between two near-equal mass, gas-free dwarf galaxies forms the observed shell systems.

We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening.

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

Low surface brightness galaxies (LSBGs) represent a significant percentage of local galaxies but their formation and evolution remain elusive. They may hold crucial information for our understanding of many key issues (i.e., census of baryonic and dark matter, star formation in the low density regime, mass function). The most massive examples - the so called giant LSBGs - can be as massive as the Milky Way, but with this mass being distributed in a much larger disk. Malin 1 is an iconic giant LSBG, perhaps the largest disk galaxy known. We attempt to bring new insights on its structure and evolution on the basis of new images covering a wide range in wavelength. We have computed surface brightness profiles (and average surface brightnesses in 16 regions of interest), in six photometric bands (FUV, NUV, u, g, i, z). We compared these data to various models, testing a variety of assumptions concerning the formation and evolution of Malin 1. We find that the surface brightness and color profiles can be reproduced by a long and quiet star-formation history due to the low surface density; no significant event, such as a collision, is necessary. Such quiet star formation across the giant disk is obtained in a disk model calibrated for the Milky Way, but with an angular momentum approximately 20 times larger. Signs of small variations of the star-formation history are indicated by the diversity of ages found when different regions within the galaxy are intercompared.For the first time, panchromatic images of Malin 1 are used to constrain the stellar populations and the history of this iconic example among giant LSBGs. Based on our model, the extreme disk of Malin 1 is found to have a long history of relatively low star formation (about 2 Msun/yr). Our model allows us to make predictions on its stellar mass and metallicity.

We present an analysis of high-quality photometry for globular clusters (GCs) in the Virgo cluster core region, based on data from the Next Generation Virgo Cluster Survey (NGVS) pilot field, and in the Milky Way (MW) based on VLT/X-Shooter spectrophotometry. We find significant discrepancies in color-color diagrams between sub-samples from different environments, confirming that the environment has a strong influence on the integrated colors of GCs. GC color distributions along a single color are not sufficient to capture the differences we observe in color-color space. While the average photometric colors become bluer with increasing radial distance to the cD galaxy M87, we also find a relation between the environment and the slope and intercept of the color-color relations. A denser environment seems to produce a larger dynamic range in certain color indices. We argue that these results are not due solely to differential extinction, IMF variations, calibration uncertainties, or overall age/metallicity variations. We therefore suggest that the relation between the environment and GC colors is, at least in part, due to chemical abundance variations, which affect stellar spectra and stellar evolution tracks. Our results demonstrate that stellar population diagnostics derived from model predictions which are calibrated on one particular sample of GCs may not be appropriate for all extragalactic GCs. These results advocate a more complex model of the assembly history of GC systems in massive galaxies that goes beyond the simple bimodality found in previous decades.

We present the Red-sequence Cluster Lensing Survey (RCSLenS), an application of the methods developed for the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) to the ~785deg$^2$, multi-band imaging data of the Red-sequence Cluster Survey 2 (RCS2). This project represents the largest public, sub-arcsecond seeing, multi-band survey to date that is suited for weak gravitational lensing measurements. With a careful assessment of systematic errors in shape measurements and photometric redshifts we extend the use of this data set to allow cross-correlation analyses between weak lensing observables and other data sets. We describe the imaging data, the data reduction, masking, multi-colour photometry, photometric redshifts, shape measurements, tests for systematic errors, and a blinding scheme to allow for more objective measurements. In total we analyse 761 pointings with r-band coverage, which constitutes our lensing sample. Residual large-scale B-mode systematics prevent the use of this shear catalogue for cosmic shear science. The effective number density of lensing sources over an unmasked area of 571.7deg$^2$ and down to a magnitude limit of r~24.5 is 8.1 galaxies per arcmin$^2$ (weighted: 5.5 arcmin$^{-2}$) distributed over 14 patches on the sky. Photometric redshifts based on 4-band griz data are available for 513 pointings covering an unmasked area of 383.5 deg$^2$ We present weak lensing mass reconstructions of some example clusters as well as the full survey representing the largest areas that have been mapped in this way. All our data products are publicly available through CADC at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/community/rcslens/query.html in a format very similar to the CFHTLenS data release.

The central region of the Virgo cluster of galaxies contains thousands of globular clusters (GCs), an order of magnitude more than the numbers found in the Local Group. Relics of early star formation epochs in the universe, these GCs also provide ideal targets to test our understanding of the Spectral Energy Distributions (SEDs) of old stellar populations. Based on photometric data from the Next Generation Virgo cluster Survey (NGVS) and its near-infrared counterpart NGVS-IR, we select a robust sample of 1846 GCs with excellent photometry and spanning the full range of colors present in the Virgo core. The selection exploits the well defined locus of GCs in the uiK diagram and the fact that the globular clusters are marginally resolved in the images. We show that the GCs define a narrow sequence in 5-dimensional color space, with limited but real dispersion around the mean sequence. The comparison of these SEDs with the predictions of eleven widely used population synthesis models highlights differences between models, and also shows that no single model adequately matches the data in all colors. We discuss possible causes for some of these discrepancies. Forthcoming papers of this series will examine how best to estimate photometric metallicities in this context, and compare the Virgo globular cluster colors with those in other environments.

Diffuse Galactic light has been observed in the optical since the 1930s. We propose that, when observed in the optical with deep imaging surveys, it can be used as a tracer of the turbulent cascade in the diffuse interstellar medium (ISM), down to scales of about 1 arcsec. Here we present a power spectrum analysis of the dust column density of a diffuse cirrus at high Galactic latitude (l ~ 198 deg, b ~ 32 deg) as derived from the combination of a MegaCam g-band image, obtained as part of the MATLAS large programme at the CFHT, with Planck radiance and Wise 12 micron data. The combination of these three datasets have allowed us to compute the density power spectrum of the HI over scales of more than three orders of magnitude. We found that the density field is well described by a single power law over scales ranging from 0.01 to 50 pc. The exponent of the power spectrum, gamma=-2.9 +- 0.1, is compatible with what is expected for thermally bi-stable and turbulent HI. We did not find any steepening of the power spectrum at small scales indicating that the typical scale at which turbulent energy is dissipated in this medium is smaller than 0.01pc. The ambipolar diffusion scenario that is usually proposed as the main dissipative agent, is consistent with our data only if the density of the cloud observed is higher than the typical values assumed for the cold neutral medium gas. We discuss the new avenue offered by deep optical imaging surveys for the study of the low density ISM structure and turbulence.

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic Plane Survey (CFEPS; Gladman et al. 2012), with which we find broad agreement. We confirm that the 5:2 resonance is more populated than models of the outer Solar System's dynamical history predict; our minimum population estimate shows that the high eccentricity (e>0.35) portion of the resonance is at least as populous as the 2:1 and possibly as populated as the 3:2 resonance. One OSSOS block was well-suited to detecting objects trapped at low libration amplitudes in Neptune's 3:2 resonance, a population of interest in testing the origins of resonant TNOs. We detected three 3:2 objects with libration amplitudes below the cutoff modeled by CFEPS; OSSOS thus offers new constraints on this distribution. The OSSOS detections confirm that the 2:1 resonance has a dynamically colder inclination distribution than either the 3:2 or 5:2 resonances. Using the combined OSSOS and CFEPS 2:1 detections, we constrain the fraction of 2:1 objects in the symmetric mode of libration to be 0.2-0.85; we also constrain the fraction of leading vs. trailing asymmetric librators, which has been theoretically predicted to vary depending on Neptune's migration history, to be 0.05-0.8. Future OSSOS blocks will improve these constraints.

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $\alpha=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $\alpha = -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $\alpha=-1.21\pm0.05$.

We analyze the kinematics of six Virgo cluster dwarf early-type galaxies (dEs) from their globular cluster (GC) systems. We present new Keck/DEIMOS spectroscopy for three of them and reanalyze the data found in the literature for the remaining three. We use two independent methods to estimate the rotation amplitude (Vmax) and velocity dispersion (sigma_GC) of the GC systems and evaluate their statistical significance by simulating non-rotating GC systems with the same number of GC satellites and velocity uncertainties. Our measured kinematics agree with the published values for the three galaxies from the literature and, in all cases, some rotation is measured. However, our simulations show that the null hypothesis of being non-rotating GC systems cannot be ruled out. In the case of VCC1861, the measured Vmax and the simulations indicate that it is not rotating. In the case of VCC1528, the null hypothesis can be marginally ruled out, thus, it might be rotating although further confirmation is needed. In our analysis, we find that, in general, the measured Vmax tends to be overestimated and the measured sigma_GC tends to be underestimated by amounts that depend on the intrinsic Vmax/sigma_GC, the number of observed GCs (N_GC), and the velocity uncertainties. The bias is negligible when N_GC>~20. In those cases where a large N_GC is not available, it is imperative to obtain data with small velocity uncertainties. For instance, errors of <2km/s lead to Vmax<10km/s for a system that is intrinsically not rotating.

(Abridged) We investigate the intrinsic shapes of low-luminosity galaxies in the central 300 kpc of the Virgo cluster using deep imaging obtained as part of the NGVS. We build a sample of nearly 300 red-sequence cluster members in the yet unexplored $-14 < M_{g} < -8$ magnitude range. The observed distribution of apparent axis ratios is then fit by families of triaxial models with normally-distributed intrinsic ellipticities and triaxialities. We develop a Bayesian framework to explore the posterior distribution of the model parameters, which allows us to work directly on discrete data, and to account for individual, surface brightness-dependent axis ratio uncertainties. For this population we infer a mean intrinsic ellipticity E=0.43, and a mean triaxiality T=0.16. This implies that faint Virgo galaxies are best described as a family of thick, nearly oblate spheroids with mean intrinsic axis ratios 1:0.94:0.57. We additionally attempt a study of the intrinsic shapes of Local Group satellites of similar luminosities. For the LG population we infer a slightly larger mean intrinsic ellipticity E=0.51, and the paucity of objects with round apparent shapes translates into more triaxial mean shapes, 1:0.76:0.49. We finally compare the intrinsic shapes of NGVS low-mass galaxies with samples of more massive quiescent systems, and with field, star-forming galaxies of similar luminosities. We find that the intrinsic flattening in this low-luminosity regime is almost independent of the environment in which the galaxy resides--but there is a hint that objects may be slightly rounder in denser environments. The comparable flattening distributions of low-luminosity galaxies that have experienced very different degrees of environmental effects suggests that internal processes are the main drivers of galaxy structure at low masses--with external mechanisms playing a secondary role.

We obtained using MegaCam at the CFHT a deep narrow band Halpha+[NII] wide field image of NGC 4569, the brightest late-type galaxy in the Virgo cluster. The image reveals the presence of long tails of diffuse ionised gas without any associated stellar component extending from the disc of the galaxy up to ~ 80 kpc (projected distance) with a typical surface brightness of a few 10^-18 erg s-1 cm-2 arcsec-2. These features provide direct evidence that NGC 4569 is undergoing a ram presure stripping event. The image also shows a prominent 8 kpc spur of ionised gas associated to the nucleus that spectroscopic data identify as an outflow. With some assumptions on the 3D distribution of the gas, we use the Halpha surface brightness of these extended low surface brightness features to derive the density and the mass of the gas stripped during the interaction of the galaxy with the ICM. The comparison with ad-hoc chemo-spectrophotometric models of galaxy evolution indicates that the mass of the Halpha emitting gas in the tail is comparable to that of the cold phase stripped from the disc, suggesting that the gas is ionised within the tail during the stripping process. The lack of star forming regions suggests that mechanisms other than photoionisation are responsible for the excitation of the gas (shocks, heat conduction, magneto hydrodynamic waves). This analysis indicates that ram pressure stripping is efficient in massive (M_star ~ 10^10.5 Mo) galaxies located in intermediate mass (~ 10^14 Mo) clusters under formation. It also shows that the mass of gas expelled by the nuclear outflow is ~ 1 % than that removed during the ram pressure stripping event. All together these results indicate that ram pressure stripping, rather than starvation through nuclear feedback, can be the dominant mechanism responsible for the quenching of the star formation activity of galaxies in high density environments.

The Next Generation Virgo Cluster Survey (NGVS) is a deep $u^*giz$ survey targeting the Virgo cluster of galaxies at 16.5~Mpc. This survey provides high-quality photometry over an $\sim$ 100 deg$^2$ region straddling the constellations of Virgo and Coma Berenices. This sightline through the Milky Way is noteworthy in that it intersects two of the most prominent substructures in the Galactic halo: the Virgo Over-Density (VOD) and Sagittarius stellar stream (close to its bifurcation point). In this paper, we use deep $u^*gi$ imaging from the NGVS to perform tomography of the VOD and Sagittarius stream using main-sequence turnoff (MSTO) stars as a halo tracer population. The VOD, whose centroid is known to lie at somewhat lower declinations ($\alpha \sim 190^\circ$, $\delta \sim -5^\circ$) than is covered by the NGVS, is nevertheless clearly detected in the NGVS footprint at distances between $\sim$ 8 and 25~kpc. By contrast, the Sagittarius stream is found to slice directly across the NGVS field at distances between 25 and 40~kpc, with a density maximum at $\simeq$~35~kpc. No evidence is found for new substructures beyond the Sagittarius stream, at least out to a distance of $\sim$~90 kpc --- the largest distance to which we can reliably trace the halo using MSTO stars. We find clear evidence for a distance gradient in the Sagittarius stream across the $\sim 30$~deg of sky covered by the NGVS and its flanking fields. We compare our distance measurements along the stream to those predicted by leading stream models.

Context. The elliptical galaxy NGC 3923 is known to be surrounded by a number of stellar shells, probable remnants of an accreted galaxy. Despite its uniqueness, the deepest images of its outskirts come from the 1980s. On the basis of the modified Newtonian dynamics (MOND), it has recently been predicted that a new shell lies in this region. Aims. We obtain the deepest image ever of the galaxy, map the tidal features in it, and search for the predicted shell. Methods. The image of the galaxy was taken by the MegaCam camera at the Canada-France-Hawaii Telescope in the g' band. It reached the surface-brightness limit of 29 mag/arcsec2. In addition, we reanalyzed an archival HST image of the galaxy. Results. We detected up to 42 shells in NGC 3923. This is by far the highest number among all shell galaxies. We present the description of the shells and other tidal features in the galaxy. A probable progenitor of some of these features was discovered. The shell system likely originates from two or more progenitors. The predicted shell was not detected, but the new image revealed that the prediction was based on incorrect assumptions and poor data.

We measure the absolute magnitude, $H$, distribution, $dN(H) \propto 10^{\alpha H}$ of the scattering Trans-Neptunian Objects (TNOs) as a proxy for their size-frequency distribution. We show that the H-distribution of the scattering TNOs is not consistent with a single-slope distribution, but must transition around $H_g \sim 9$ to either a knee with a shallow slope or to a divot, which is a differential drop followed by second exponential distribution. Our analysis is based on a sample of 22 scattering TNOs drawn from three different TNO surveys, the Canada-France Ecliptic Plane Survey (CFEPS, Petit et al. 2011), Alexandersen et al. (2014), and the Outer Solar System Origins Survey (OSSOS, Bannister et al. 2016), all of which provide well characterized detection thresholds, combined with a cosmogonic model for the formation of the scattering TNO population. Our measured absolute magnitude distribution result is independent of the choice of cosmogonic model. Based on our analysis, we estimate that number of scattering TNOs is (2.4-8.3)$\times 10^5$ for $H_r < 12$. A divot $H$-distribution is seen in a variety of formation scenarios and may explain several puzzles in Kuiper Belt science. We find that a divot $H$-distribution simultaneously explains the observed scattering TNO, Neptune Trojan, Plutino, and Centaur $H$-distributions while simultaneously predicting a large enough scattering TNO population to act as the sole supply of the Jupiter-Family Comets.

We use the GALEX data of the GUViCS survey to construct the NUV luminosity function of the Virgo cluster over ~ 300 deg.2, an area covering the cluster and its surrounding regions up to ~ 1.8 virial radii. The NUV luminosity function is also determined for galaxies of different morphological type and NUV-i colour, and for the different substructures within the cluster. These luminosity functions are robust vs. statistical corrections since based on a sample of 833 galaxies mainly identified as cluster members with spectroscopic redshift (808) or high-quality optical scaling relations (10). We fit these luminosity functions with a Schechter function, and compare the fitted parameters with those determined for other nearby clusters and for the field. The faint end slope of the Virgo NUV luminosity function (alpha = -1.19), here sampled down to ~ NUV = -11.5 mag, is significantly flatter than the one measured in other nearby clusters and similar to the field one. Similarly M* = -17.56 is one-to-two magnitudes fainter than measured in Coma, A1367, the Shapley supercluster, and the field. These differences seem due to the quite uncertain statistical corrections and the small range in absolute magnitude sampled in these clusters. We do not observe strong systematic differences in the overall NUV luminosity function of the core of the cluster with respect to that of its periphery. We notice, however, that the relative contribution of red-to-blue galaxies at the faint end is inverted, with red quiescent objects dominating the core of the cluster and star forming galaxies dominating beyond one virial radius. These observational evidences are discussed in the framework of galaxy evolution in dense environments.

We report on the discovery of an ultrasoft X-ray transient source, 3XMM J152130.7+074916. It was serendipitously detected in an XMM-Newton observation on 2000 August 23, and its location is consistent with the center of the galaxy SDSS J152130.72+074916.5 (z=0.17901 and d_L=866 Mpc). The high-quality X-ray spectrum can be fitted with a thermal disk with an apparent inner disk temperature of 0.17 keV and a rest-frame 0.24-11.8 keV unabsorbed luminosity of ~5e43 erg/s, subject to a fast-moving warm absorber. Short-term variability was also clearly observed, with the spectrum being softer at lower flux. The source was covered but not detected in a Chandra observation on 2000 April 3, a Swift observation on 2005 September 10, and a second XMM-Newton observation on 2014 January 19, implying a large variability (>260) of the X-ray flux. The optical spectrum of the candidate host galaxy, taken ~11 yrs after the XMM-Newton detection, shows no sign of nuclear activity. This, combined with its transient and ultrasoft properties, leads us to explain the source as tidal disruption of a star by the supermassive black hole in the galactic center. We attribute the fast-moving warm absorber detected in the first XMM-Newton observation to the super-Eddington outflow associated with the event and the short-term variability to a disk instability that caused fast change of the inner disk radius at a constant mass accretion rate.

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$\alpha$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

We use imaging from the Next Generation Virgo cluster Survey (NGVS) to present a comparative study of ultra-compact dwarf (UCD) galaxies associated with three prominent Virgo sub-clusters: those centered on the massive, red-sequence galaxies M87, M49 and M60. We show how UCDs can be selected with high completeness using a combination of half-light radius and location in color-color diagrams ($u^*iK_s$ or $u^*gz$). Although the central galaxies in each of these sub-clusters have nearly identical luminosities and stellar masses, we find large differences in the sizes of their UCD populations, with M87 containing ~3.5 and 7.8 times more UCDs than M49 and M60, respectively. The relative abundance of UCDs in the three regions scales in proportion to sub-cluster mass, as traced by X-ray gas mass, total gravitating mass, number of globular clusters, and number of nearby galaxies. We find that the UCDs are predominantly blue in color, with ~85% of the UCDs having colors similar to blue GCs and stellar nuclei of dwarf galaxies. We present evidence that UCDs surrounding M87 and M49 may follow a morphological sequence ordered by the prominence of their outer, low surface brightness envelope, ultimately merging with the sequence of nucleated low-mass galaxies, and that envelope prominence correlates with distance from either galaxy. Our analysis provides evidence that tidal stripping of nucleated galaxies is an important process in the formation of UCDs.

We report the discovery of three large (R29 >~ 1 arcminute) extremely low surface brightness (mu_(V,0) ~ 27.0) galaxies identified using our deep, wide-field imaging of the Virgo Cluster from the Burrell Schmidt telescope. Complementary data from the Next Generation Virgo Cluster Survey do not resolve red giant branch stars in these objects down to i=24, yielding a lower distance limit of 2.5 Mpc. At the Virgo distance, these objects have half-light radii 3-10 kpc and luminosities L_V=2-9x10^7 Lsun. These galaxies are comparable in size but lower in surface brightness than the large ultradiffuse LSB galaxies recently identified in the Coma cluster, and are located well within Virgo's virial radius; two are projected directly on the cluster core. One object appears to be a nucleated LSB in the process of being tidally stripped to form a new Virgo ultracompact dwarf galaxy. The others show no sign of tidal disruption, despite the fact that such objects should be most vulnerable to tidal destruction in the cluster environment. The relative proximity of Virgo makes these objects amenable to detailed studies of their structural properties and stellar populations. They thus provide an important new window onto the connection between cluster environment and galaxy evolution at the extremes.

The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$.

A fraction of the early-type galaxy population hosts a prominent dust lane. Methods to quantify the dust content of these systems based on optical imaging data usually yield dust masses which are an order of magnitude lower than dust masses derived from the observed FIR emission. High-quality optical data from the Next Generation Virgo cluster Survey (NGVS) and FIR/submm observations from the Herschel Virgo Cluster Survey (HeViCS) allow us to revisit previous methods to determine the dust content in galaxies and explore new ones. We aim to derive the dust mass in NGC 4370 from both optical and FIR data, and investigate the need to invoke a putative diffuse dust component. We create color and attenuation maps, which are converted to approximate dust mass maps based on simple dust geometries. Dust masses are also derived from SED fits to FIR/submm observations. Finally, inverse radiative transfer fitting is performed to investigate more complex dust geometries. The empirical methods applied to the optical data yield lower limits of 3.4e5 solar masses, an order of magnitude below the total dust masses derived from SED fitting. In contrast, radiative transfer models yield dust masses which are slightly lower, but fully consistent with the FIR-derived mass. Dust is more likely to be distributed in a ring around the centre of NGC 4370 as opposed to an exponential disc or a simple foreground screen. Moreover, using inverse radiative transfer fitting, we are able to constrain most of the parameters describing these geometries. The resulting dust masses are high enough to account for the dust observed at FIR/submm wavelengths, so that no diffuse dust component needs to be invoked. We furthermore caution for the interpretation of dust masses and optical depths based on optical data alone, using overly simplistic star-dust geometries and the neglect of scattering effects. [ABRIDGED]

We present Gemini GMOS-IFU data of eight compact low-mass early-type galaxies (ETGs) in the Virgo cluster. We analyse their stellar kinematics, stellar population, and present two-dimensional maps of these properties covering the central 5"x 7" region. We find a large variety of kinematics: from non- to highly-rotating objects, often associated with underlying disky isophotes revealed by deep images from the Next Generation Virgo Cluster Survey. In half of our objects, we find a centrally-concentrated younger and more metal-rich stellar population. We analyze the specific stellar angular momentum through the lambdaR parameter and find six fast-rotators and two slow-rotators, one having a thin counter-rotating disk. We compare the local galaxy density and stellar populations of our objects with those of 39 more extended low-mass Virgo ETGs from the SMAKCED survey and 260 massive ($M>10^{10}$\Msun) ETGs from the A3D sample. The compact low-mass ETGs in our sample are located in high density regions, often close to a massive galaxy and have, on average, older and more metal-rich stellar populations than less compact low-mass galaxies. We find that the stellar population parameters follow lines of constant velocity dispersion in the mass-size plane, smoothly extending the comparable trends found for massive ETGs. Our study supports a scenario where low-mass compact ETGs have experienced long-lived interactions with their environment, including ram-pressure stripping and gravitational tidal forces, that may be responsible for their compact nature.

We present the results of deep \chandra\ imaging of the central region of the Extended Groth Strip, the AEGIS-X Deep (AEGIS-XD) survey. When combined with previous \chandra\ observations of a wider area of the strip, AEGIS-X Wide (AEGIS-XW; Laird et~al. 2009), these provide data to a nominal exposure depth of 800ks in the three central ACIS-I fields, a region of approximately $0.29$~deg$^{2}$. This is currently the third deepest X-ray survey in existence, a factor $\sim 2-3$ shallower than the Chandra Deep Fields (CDFs) but over an area $\sim 3$ times greater than each CDF. We present a catalogue of 937 point sources detected in the deep \chandra\ observations. We present identifications of our X-ray sources from deep ground-based, Spitzer, GALEX and HST imaging. Using a likelihood ratio analysis, we associate multi band counterparts for 929/937 of our X-ray sources, with an estimated 95~\% reliability, making the identification completeness approximately 94~\% in a statistical sense. Reliable spectroscopic redshifts for 353 of our X-ray sources are provided predominantly from Keck (DEEP2/3) and MMT Hectospec, so the current spectroscopic completeness is $\sim 38$~per cent. For the remainder of the X-ray sources, we compute photometric redshifts based on multi-band photometry in up to 35 bands from the UV to mid-IR. Particular attention is given to the fact that the vast majority the X-ray sources are AGN and require hybrid templates. Our photometric redshifts have mean accuracy of $\sigma=0.04$ and an outlier fraction of approximately 5\%, reaching $\sigma=0.03$ with less than 4\% outliers in the area covered by CANDELS . The X-ray, multi-wavelength photometry and redshift catalogues are made publicly available.

The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.

The Trans-Neptunian Objects (TNOs) may preserve evidence of planet building in their orbital and size-distributions. While all populations show steep size-distributions for large objects, recently relative deficit of Neptunian Trojans and scattering objects with diameters $D<100\,\mathrm{km}$ were detected. We have investigated this deficit with a 32 square degree survey, detecting 77 TNOs to a limiting $r$-band magnitude of 24.6. Our Plutinos sample (18 objects in 3:2 mean motion resonance with Neptune) also shows a deficit of $D<100\,\mathrm{km}$ objects. We reject a single power-law size-distribution and find that the Plutinos favour a divot. The Plutinos are thus added the list of populations with a deficit of $D<100\,\mathrm{km}$ objects. The fact that three independent samples of three different populations show this trend suggests that it is a real feature, possibly shared by all hot TNO populations as a remnant of "born big" planetesimal formation processes. We surmise the existence of $9000\pm3000$ Plutinos with $H_r\leq8.66$ and $37000^{+12000}_{-10000}$ Plutinos with $H_r\leq10.0$. Our survey also discovered one temporary Uranian Trojan, one temporary Neptunian Trojans and one stable Neptunian Trojan, from which we derive populations of $110^{+500}_{-100}$, $210^{+900}_{-200}$ and $150^{+600}_{-140}$ , respectively, with $H_r\leq10.0$. The Neptunian Trojans are thus less numerous than the main belt asteroids, which has over 700 asteroids with $H_r\leq10.0$. With such numbers, the temporary Neptunian Trojans cannot be previously stable Trojans that happen to be escaping the resonance now; they must be captured from another reservoir. With three 3:1 and one 4:1 resonators, we add to the growing evidence that the high-order resonances are more populated than expected.

The New Horizons spacecraft, launched by NASA in 2006, will arrive in the Pluto-Charon system on July 14, 2015. There, it will spend a few hours imaging Pluto and its moons. It will then have a small amount of reserve propellant which will be used to direct the probe on to a second, yet to be discovered object in the Kuiper Belt. Data from the MegaPrime camera on CFHT was used to build a precise, high density astrometric reference frame for both the final approach into the Pluto system and the search for the secondary target. Pluto currently lies in the galactic plane. This is a hindrance in that there are potential problems with confusion. However, it is also a benefit, since it allows the use of the UCAC4 astrometric reference catalog, which is normally too sparse for use with MegaCam images. The astrometric accuracy of the final catalogs, as measured by the residuals, is 0.02 arcseconds.

The Next Generation Virgo Cluster Survey is an optical imaging survey covering 104 deg^2 centered on the Virgo cluster. Currently, the complete survey area has been observed in the u*giz-bands and one third in the r-band. We present the photometric redshift estimation for the NGVS background sources. After a dedicated data reduction, we perform accurate photometry, with special attention to precise color measurements through point spread function-homogenization. We then estimate the photometric redshifts with the Le Phare and BPZ codes. We add a new prior which extends to iAB = 12.5 mag. When using the u*griz-bands, our photometric redshifts for 15.5 \le i \lesssim 23 mag or zphot \lesssim 1 galaxies have a bias |\Delta z| < 0.02, less than 5% outliers, and a scatter \sigma_{outl.rej.} and an individual error on zphot that increase with magnitude (from 0.02 to 0.05 and from 0.03 to 0.10, respectively). When using the u*giz-bands over the same magnitude and redshift range, the lack of the r-band increases the uncertainties in the 0.3 \lesssim zphot \lesssim 0.8 range (-0.05 < \Delta z < -0.02, \sigma_{outl.rej} ~ 0.06, 10-15% outliers, and zphot.err. ~ 0.15). We also present a joint analysis of the photometric redshift accuracy as a function of redshift and magnitude. We assess the quality of our photometric redshifts by comparison to spectroscopic samples and by verifying that the angular auto- and cross-correlation function w(\theta) of the entire NGVS photometric redshift sample across redshift bins is in agreement with the expectations.

We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution. We present the clearest evidence to date showing the difference in concentration between red and blue GCs over the extent of the cluster, where the red (metal-rich) GCs are largely located around the massive early-type galaxies, whilst the blue (metal-poor) GCs have a more extended spatial distribution, with significant populations present beyond 83' (215 kpc) along the major axes of M49 and M87. The GC distribution around M87 and M49 shows remarkable agreement with the shape, ellipticity and boxiness of the diffuse light surrounding both galaxies. We find evidence for spatial enhancements of GCs surrounding M87 that may be indicative of recent interactions or an ongoing merger history. We compare the GC map to the locations of Virgo galaxies and the intracluster X-ray gas, and find good agreement between these baryonic structures. The Virgo cluster contains a total population of 67300$\pm$14400 GCs, of which 35% are located in M87 and M49 alone. We compute a cluster-wide specific frequency S_N,CL=$2.8\pm0.7$, including Virgo's diffuse light. The GC-to-baryonic mass fraction is e_b=$5.7\pm1.1\times10^{-4} $and the GC-to-total cluster mass formation efficiency is e_t=$2.9\pm0.5\times10^{-5}$, values slightly lower than, but consistent with, those derived for individual galactic halos. Our results show that the production of the complex structures in the unrelaxed Virgo cluster core (including the diffuse intracluster light) is an ongoing process.(abridged)

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

In this paper we introduce the deepest and most extensive ultraviolet extragalactic source catalogs of the Virgo Cluster area to date. Archival and targeted GALEX imaging is compiled and combined to provide the deepest possible coverage over ~120 deg^2 in the NUV (lambda_eff=2316 angstroms) and ~40 deg^2 in the FUV (lambda_eff=1539 angstroms) between 180 deg <= R.A. <= 195 deg and 0 deg <= Decl. <= 20 deg. We measure the integrated photometry of 1770 extended UV sources of all galaxy types and use GALEX pipeline photometry for 1,230,855 point-like sources in the foreground, within, and behind the cluster. Extended source magnitudes are reliable to m_UV ~22, showing ~0.01 sigma difference from their asymptotic magnitudes. Point-like source magnitudes have a 1 sigma standard deviation within ~0.2 mag down to m_uv ~23. The point-like source catalog is cross-matched with large optical databases and surveys including the SDSS DR9 (> 1 million Virgo Cluster sources), the Next Generation Virgo Cluster Survey (NGVS; >13 million Virgo Cluster sources), and the NED (~30,000 sources in the Virgo Cluster). We find 69% of the entire UV point-like source catalog has a unique optical counterpart, 11% of which are stars and 129 are Virgo cluster members neither in the VCC nor part of the bright CGCG galaxy catalog (i.e., m_pg < 14.5). These data are collected in four catalogs containing the UV extended sources, the UV point-like sources, and two catalogs each containing the most relevant optical parameters of UV-optically matched point-like sources for further studies from SDSS and NGVS. The GUViCS catalogs provide a unique set of data for future works on UV and multiwavelength studies in the cluster and background environments.

The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5\sigma limiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with \mu_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies.

Trojan objects share a planet's orbit, never straying far from the triangular Lagrangian points, 60 degrees ahead of (L4) or behind (L5) the planet. We report the detection of a Uranian Trojan; in our numerical integrations, 2011 QF99 oscillates around the Uranian L4 Lagrange point for >70 kyr and remains co-orbital for ~1 Myr before becoming a Centaur. We constructed a Centaur model, supplied from the transneptunian region, to estimate temporary co-orbital capture frequency and duration (to factor of two accuracy), finding that at any time 0.4% and 2.8% of the population will be Uranian and Neptunian co-orbitals, respectively. The co-orbital fraction (~2.4%) among Centaurs in the IAU Minor Planet Centre database is thus as expected under transneptunian supply.

We report the discovery of 2010 GB$_{174}$, a likely new member of the Inner Oort Cloud (IOC). 2010 GB$_{174}$ is one of 91 Trans Neptunian Objects (TNOs) and Centaurs discovered in a 76 deg$^2$ contiguous region imaged as part of the Next Generation Virgo Cluster Survey (NGVS) --- a moderate ecliptic latitude survey reaching a mean limiting magnitude of $g^\prime \simeq 25.5$ --- using MegaPrime on the 3.6m Canada France Hawaii Telescope. 2010 GB$_{174}$ is found to have an orbit with semi-major axis $a\simeq350.8$ AU, inclination $i \simeq 21.6^\circ$ and pericentre $q\sim48.5$ AU. This is the second largest perihelion distance among known solar system objects. Based on the sky coverage and depth of the NGVS, we estimate the number of IOC members with sizes larger than 300 km ($H_V \le 6.2$ mag) to be $\simeq 11\,000$. A comparison of the detection rate from the NGVS and the PDSSS (a characterized survey that `re-discovered' the IOC object Sedna) gives, for an assumed a power-law LF for IOC objects, a slope of $\alpha \simeq 0.7 \pm 0.2$, with only two detections in this region this slope estimate is highly uncertain.

Here we present high cadence photometry taken by the Acquisition Camera on Gemini South, of a close passage by the $\sim540$ km radius Kuiper Belt Object, (50000) Quaoar, of a r'=20.2 background star. Observations before and after the event show that the apparent impact parameter of the event was $0.019\pm0.004$", corresponding to a close approach of $580\pm120$ km to the centre of Quaoar. No signatures of occultation by either Quaoar's limb or its potential atmosphere are detectable in the relative photometry of Quaoar and the target star, which were unresolved during closest approach. From this photometry we are able to put constraints on any potential atmosphere Quaoar might have. Using a Markov chain Monte Carlo and likelihood approach, we place pressure upper limits on sublimation supported, isothermal atmospheres of pure N$_2$, CO, and CH$_4$. For N$_2$ and CO, the upper limit surface pressures are 1 and 0.7 $\mu{bar}$ respectively. The surface temperature required for such low sublimation pressures is $\sim33$ K, much lower than Quaoar's mean temperature of $\sim44$ K measured by others. We conclude that Quaoar cannot have an isothermal N$_2$ or CO atmosphere. We cannot eliminate the possibility of a CH$_4$ atmosphere, but place upper surface pressure and mean temperature limits of $\sim138$ nbar and $\sim44$ K respectively.

Here we present observations of 7 large Kuiper Belt Objects. From these observations, we extract a point source catalog with $\sim0.01"$ precision, and astrometry of our target Kuiper Belt Objects with $0.04-0.08"$ precision within that catalog. We have developed a new technique to predict the future occurrence of stellar occultations by Kuiper Belt Objects. The technique makes use of a maximum likelihood approach which determines the best-fit adjustment to cataloged orbital elements of an object. Using simulations of a theoretical object, we discuss the merits and weaknesses of this technique compared to the commonly adopted ephemeris offset approach. We demonstrate that both methods suffer from separate weaknesses, and thus, together provide a fair assessment of the true uncertainty in a particular prediction. We present occultation predictions made by both methods for the 7 tracked objects, with dates as late as 2015. Finally, we discuss observations of three separate close passages of Quaoar to field stars, which reveal the accuracy of the element adjustment approach, and which also demonstrate the necessity of considering the uncertainty in stellar position when assessing potential occultations.

We present an investigation into the origins of a series of interlaced narrow filamentary stellar structures, loops and plumes in the vicinity of the Virgo Cluster, edge-on spiral galaxy, NGC 4216 that were previously identified by the Blackbird Telescope. Using the deeper, higher-resolution and precisely calibrated optical CFHT/MegaCam images obtained as part of the Next Generation Virgo Cluster Survey (NGVS), we confirm the previously identified features and identify a few additional structures. The NGVS data allowed us to make a physical study of these low-surface brightness features and investigate their origin. The likely progenitors of the structures were identified as either already catalogued VCC dwarfs or newly discovered satellites caught in the act of being destroyed. They have the same g-i color index and likely contain similar stellar populations. The alignment of three dwarfs along an apparently single stream is intriguing, and we cannot totally exclude that these are second-generation dwarf galaxies being born inside the filament from the debris of an original dwarf. The observed complex structures, including in particular a stream apparently emanating from a satellite of a satellite, point to a high rate of ongoing dwarf destruction/accretion in the region of the Virgo Cluster where NGC 4216 is located. We discuss the age of the interactions and whether they occurred in a group that is just falling into the cluster and shows signs of so-called "pre-processing" before it gets affected by the cluster environment, or in a group which already ventured towards the central regions of Virgo Cluster.

We present the discovery of a long-term stable L5 (trailing) Neptune Trojan in data acquired to search for candidate Trans-Neptunian objects for the New Horizons spacecraft to fly by during an extended post-Pluto mission. This Neptune Trojan, 2011 HM102, has the highest inclination (29.4 degrees) of any known member of this population. It is intrinsically brighter than any single L5 Jupiter Trojan at H~8.18. We have determined its gri colors (a first for any L5 Neptune Trojan), which we find to be similar to the moderately red colors of the L4 Neptune Trojans, suggesting similar surface properties for members of both Trojan clouds. We also present colors derived from archival data for two L4 Neptune Trojans (2006 RJ103 and 2007 VL305), better refining the overall color distribution of the population. In this document we describe the discovery circumstances, our physical characterization of 2011 HM102, and this object's implications for the Neptune Trojan population overall. Finally, we discuss the prospects for detecting 2011 HM102 from the New Horizons spacecraft during their close approach in mid- to late-2013.

We present the results of a search for extended X-ray sources and their corresponding galaxy groups from 800-ks Chandra coverage of the All-wavelength Extended Groth Strip International Survey (AEGIS). This yields one of the largest X-ray selected galaxy group catalogs from a blind survey to date. The red-sequence technique and spectroscopic redshifts allow us to identify 100$%$ of reliable sources, leading to a catalog of 52 galaxy groups. The groups span the redshift range $z\sim0.066-1.544$ and virial mass range $M_{200}\sim1.34\times 10^{13}-1.33\times 10^{14}M_\odot$. For the 49 extended sources which lie within DEEP2 and DEEP3 Galaxy Redshift Survey coverage, we identify spectroscopic counterparts and determine velocity dispersions. We select member galaxies by applying different cuts along the line of sight or in projected spatial coordinates. A constant cut along the line of sight can cause a large scatter in scaling relations in low-mass or high-mass systems depending on the size of cut. A velocity dispersion based virial radius can more overestimate velocity dispersion in comparison to X-ray based virial radius for low mass systems. There is no significant difference between these two radial cuts for more massive systems. Independent of radial cut, overestimation of velocity dispersion can be created in case of existence of significant substructure and also compactness in X-ray emission which mostly occur in low mass systems. We also present a comparison between X-ray galaxy groups and optical galaxy groups detected using the Voronoi-Delaunay method (VDM) for DEEP2 data in this field.

This paper describes a new catalog that supplements the existing DEEP2 Galaxy Redshift Survey photometric and spectroscopic catalogs with ugriz photometry from two other surveys; the Canada-France-Hawaii Legacy Survey (CFHTLS) and the Sloan Digital Sky Survey (SDSS). Each catalog is cross-matched by position on the sky in order to assign ugriz photometry to objects in the DEEP2 catalogs. We have recalibrated the CFHTLS photometry where it overlaps DEEP2 in order to provide a more uniform dataset. We have also used this improved photometry to predict DEEP2 BRI photometry in regions where only poorer measurements were available previously. In addition, we have included improved astrometry tied to SDSS rather than USNO-A2.0 for all DEEP2 objects. In total this catalog contains ~27,000 objects with full ugriz photometry as well as robust spectroscopic redshift measurements, 64% of which have r > 23. By combining the secure and accurate redshifts of the DEEP2 Galaxy Redshift Survey with ugriz photometry, we have created a catalog that can be used as an excellent testbed for future photo-z studies, including tests of algorithms for surveys such as LSST and DES.

(Abridged) Distant galaxy clusters provide important tests of the growth of large scale structure in addition to highlighting the process of galaxy evolution in a consistently defined environment at large look back time. We present a sample of 22 distant (z>0.8) galaxy clusters and cluster candidates selected from the 9 deg2 footprint of the overlapping X-ray Multi Mirror (XMM) Large Scale Structure (LSS), CFHTLS Wide and Spitzer SWIRE surveys. Clusters are selected as extended X-ray sources with an accompanying overdensity of galaxies displaying optical to mid-infrared photometry consistent with z>0.8. Nine clusters have confirmed spectroscopic redshifts in the interval 0.8<z<1.2, four of which are presented here for the first time. A further 11 candidate clusters have between 8 and 10 band photometric redshifts in the interval 0.8<z<2.2, while the remaining two candidates do not have information in sufficient wavebands to generate a reliable photometric redshift. All of the candidate clusters reported in this paper are presented for the first time. Those confirmed and candidate clusters with available near infrared photometry display evidence for a red sequence galaxy population, determined either individually or via a stacking analysis, whose colour is consistent with the expectation of an old, coeval stellar population observed at the cluster redshift. We further note that the sample displays a large range of red fraction values indicating that the clusters may be at different stages of red sequence assembly. We compare the observed X-ray emission to the flux expected from a suite of model clusters and find that the sample displays an effective mass limit M200 ~ 1e14 Msolar with all clusters displaying masses consistent with M200 < 5e14 Msolar. This XMM distant cluster study represents a complete sample of X-ray selected z>0.8 clusters.

We present data products from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). CFHTLenS is based on the Wide component of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). It encompasses 154 deg^2 of deep, optical, high-quality, sub-arcsecond imaging data in the five optical filters u^*g'r'i'z'. The article presents our data processing of the complete CFHTLenS data set. We were able to obtain a data set with very good image quality and high-quality astrometric and photometric calibration. Our external astrometric accuracy is between 60-70 mas with respect to SDSS data and the internal alignment in all filters is around 30 mas. Our average photometric calibration shows a dispersion on the order of 0.01 to 0.03 mag for g'r'i'z' and about 0.04 mag for u^* with respect to SDSS sources down to i <= 21. In the spirit of the CFHTLS all our data products are released to the astronomical community via the Canadian Astronomy Data Centre. We give a description and how-to manuals of the public products which include image pixel data, source catalogues with photometric redshift estimates and all relevant quantities to perform weak lensing studies.

Context: Several galaxies in the Virgo cluster are known to have large HI gas tails related to a recent ram-pressure stripping event. The Virgo cluster has been extensively observed at 1539 A in the far-ultraviolet for the GALEX Ultraviolet Virgo Cluster Survey (GUViCS), and in the optical for the Next Generation Virgo Survey (NGVS), allowing a study of the stellar emission potentially associated with the gas tails of 8 cluster members. On the theoretical side, models of ram-pressure stripping events have started to include the physics of star formation. Aim: We aim to provide quantitative constraints on the amount of star formation taking place in the ram-pressure stripped gas, mainly on the basis of the far-UV emission found in the GUViCS images in relation with the gas content of the tails. Methods: We have performed three comparisons of the young stars emission with the gas column density: visual, pixel-by-pixel and global. We have compared our results to other observational and theoretical studies. Results: We find that the level of star formation taking place in the gas stripped from galaxies by ram-pressure is low with respect to the available amount of gas. Star formation is lower by at least a factor 10 compared to the predictions of the Schmidt Law as determined in regular spiral galaxy disks. It is also lower than measured in dwarfs galaxies and the outer regions of spirals, and than predicted by some numerical simulations. We provide constraints on the star formation efficiency in the ram-pressure stripped gas tails, and compare these with current models.

Context: We study the peculiar interacting galaxy system of VCC1249/M49 located in the core of the Virgo B subcluster. Owing to a recent interaction between the dwarf galaxy VCC1249 and the halo gas of the gE M49, neutral hydrogen has been displaced from the interstellar medium of this dwarf into the Virgo ICM. Observations also reveal multiple compact star-forming regions that are embedded in this HI cloud, with a projected separation up to 13 kpc from VCC1249 in the northwest direction. Aims: Motivated by recent NUV imaging from GUViCS of the VCC1249/M49 system that shows significant ongoing/recent star formation in the compact regions, we aim to constrain the origin of these outlying HII regions with a multi-wavelength approach. Methods: Using deep optical (u, g, i, z) imaging from NGVS and new Halpha imaging obtained at the San Pedro Martir observatory together with Keck long-slit spectroscopy, we characterize the SFR, ages, and metallicity of VCC1249 and its outlying compact regions. Moreover, we analyze the color and luminosity profile of the galaxy to investigate its recent interaction with M49. Results: Our new observations indicate that VCC1249 underwent a recent interaction with M49 in which both ram-pressure stripping and tidal interaction occured. The joint action of the two mechanisms led to the removal of the HI gas from the ISM of VCC1249, while the gravitational tides triggered the stellar tail and counter-tail of VCC1249. Our SED analysis reveals that the star formation in this galaxy was truncated around 200 Myr ago and that the outlying HII regions were born in situ about 10 Myr ago out of pre-enriched gas removed from the dwarf galaxy. These observations also reveal that interactions between central and satellite galaxies similar to the one between VCC1249/M49 may be an effective way of dispersing metals into the halos of massive galaxies.

It is very difficult to find archival images of solar system objects. While regular archive searches can find images at a fixed location, they cannot find images of moving targets. Archival images have become increasingly useful to galactic and stellar astronomers the last few years but, until now, solar system researchers have been at a disadvantage in this respect. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for images containing that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, several ESO instruments and the SDSS for a total of 6.5 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos

While regular archive searches can find images at a fixed location, they cannot find images of moving targets such as asteroids. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for observations of that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, and several ESO instruments for a total of 1.6 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos

This paper describes the image stacks and catalogs of the CFHT Legacy Survey produced using using the MegaPipe data pipeline at the Canadian Astronomy Data Centre. The Legacy Survey is divided into to two parts: The Deep Survey consists of 4 fields each of 1 square degree, with magnitude limits of u=27.5, g=27.9, r=27.7, i=27.4, and z=26.2, and contains 1.6 million sources. The Wide Survey consists of 171 square degree pointings split over 4 fields, with magnitude limits of u=26.0, g=26.5, r=25.9, i=25.7, and z=24.6, and contains 30 million sources. This paper describes the calibration, image stacking and catalog generation process. The images and catalogs are available on the web through several interfaces: normal image and text file catalog downloads, a Google Sky interface, an image cutout service, and a catalog database query service.

The mass assembly of galaxies leaves various imprints on their surroundings, such as shells, streams and tidal tails. The frequency and properties of these fine structures depend on the mechanism driving the mass assembly: e.g. a monolithic collapse, rapid cold-gas accretion followed by violent disk instabilities, minor mergers or major dry / wet mergers. Therefore, by studying the outskirts of galaxies, one can learn about their main formation mechanism. I present here our on-going work to characterize the outskirts of Early-Type Galaxies (ETGs), which are powerful probes at low redshift of the hierarchical mass assembly of galaxies. This work relies on ultra-deep optical images obtained at CFHT with the wide-field of view MegaCam camera of field and cluster ETGs obtained as part of the Atlas-3D and NGVS projects. State of the art numerical simulations are used to interpret the data. The images reveal a wealth of unknown faint structures at levels as faint as 29 mag arcsec-2 in the g-band. Initial results for two galaxies are presented here.

The Next Generation Virgo Cluster Survey (NGVS) is a CFHT Large Program that is using the wide field of view capabilities of the MegaCam camera to map the entire Virgo Cluster from its core to virial radius. The observing strategy has been optimized to detect very low surface brightness structures in the cluster, including intracluster stellar streams and faint dwarf spheroidal galaxies. We present here the current status of this ongoing survey, with an emphasis on the detection and analysis of the very low-mass galaxies in the cluster that have been revealed by the NGVS.

We have combined the large SN Ia database of the Canada-France-Hawaii Telescope Supernova Legacy Survey and catalogs of galaxies with photometric redshifts, VLA 1.4 GHz radio sources, and Spitzer infrared sources. We present eight SNe Ia in early-type host galaxies which have counterparts in the radio and infrared source catalogs. We find the SN Ia rate in subsets of radio and infrared early-type galaxies is ~1-5 times the rate in all early-type galaxies, and that any enhancement is always <~ 2 sigma. Rates in these subsets are consistent with predictions of the two component "A+B" SN Ia rate model. Since infrared properties of radio SN Ia hosts indicate dust obscured star formation, we incorporate infrared star formation rates into the "A+B" model. We also show the properties of SNe Ia in radio and infrared galaxies suggest the hosts contain dust and support a continuum of delay time distributions for SNe Ia, although other delay time distributions cannot be ruled out based on our data.

We present the results of our Thousand Asteroid Light Curve Survey (TALCS) conducted with the Canada-France-Hawaii Telescope in September 2006. Our untargeted survey detected 828 Main Belt asteroids to a limiting magnitude of g'~22.5 corresponding to a diameter range of 0.4 km <= D <= 10 km. Of these, 278 objects had photometry of sufficient quality to perform rotation period fits. We debiased the observations and light curve fitting process to determine the true distribution of rotation periods and light curve amplitudes of Main Belt asteroids. We confirm a previously reported excess in the fraction of fast rotators but find a much larger excess of slow rotating asteroids (~15% of our sample). A few percent of objects in the TALCS size range have large light curve amplitudes of ~1 mag. Fits to the debiased distribution of light curve amplitudes indicate that the distribution of triaxial ellipsoid asteroid shapes is proportional to the square of the axis-ratio, (b/a)^2, and may be bi-modal. Finally, we find six objects with rotation periods that may be less than 2 hours with diameters between 400 m and 1.5 km, well above the break-up limit for a gravitationally-bound aggregate. Our debiased data indicate that this population represents <4% of the Main Belt in the 1-10 km size range.

(Abridged) We present new optical and near-infrared imaging for a sample of 98 spectroscopically-selected galaxy groups at 0.25<z<0.55. We measure accurate colours for group members and the surrounding field population, statistically complete above a stellar mass limit of M=1E10 Msun. The overall colour distribution is bimodal in both the field and group samples; but at fixed luminosity the fraction of group galaxies populating the red peak is larger, by 20+/-7 per cent, than that of the field. In particular, group members with early-type morphologies, as identified in HST imaging, exhibit a tight red sequence, similar to that seen for more massive clusters. We show that approximately 20-30 per cent of galaxies on the red sequence may be dust-reddened galaxies with non-negligible star formation and early-spiral morphologies. This is true of both the field and group sample, and shows little dependence on near infrared luminosity. Thus, the fraction of bright group members with no sign of star formation or AGN activity is 54+/-6 per cent. Our field sample, which includes galaxies in all environments, contains 35+/-3 per cent of such inactive galaxies, consistent with the amount expected if all such galaxies are located in groups and clusters. This reinforces our earlier conclusions, that dense environments at z<0.5 are associated with a premature cessation of star formation in some galaxies; in particular we find no evidence for significantly enhanced star formation in these environments. Simple galaxy formation models predict a quenching of star formation in groups that is too efficient, overpopulating the red sequence. Attempts to fix this by increasing the timescale of this quenching equally for all group members distorts the colour distribution in a way that is inconsistent with observations.

This paper describes the MegaPipe image processing pipeline at the Canadian Astronomical Data Centre (CADC). The pipeline takes multiple images from the MegaCam mosaic camera on CFHT and combines them into a single output image. MegaPipe takes as input detrended MegaCam images and does a careful astrometric and photometric calibration on them. The calibrated images are then resampled and combined into image stacks. MegaPipe is run on PI data by request, data from large surveys (the CFHT Legacy Survey and the Next Generation Virgo Survey) and all non-proprietary MegaCam data in the CFHT archive. The stacked images and catalogs derived from these images are available through the CADC website. Currently, 1500 square degrees have been processed.

The formation and growth of galaxy disks over cosmic time is crucial to our understanding of galaxy formation. Despite steady improvements in the size and quality of disk samples over the last decade, many aspects of galaxy disk evolution remain unclear. Using two square degrees of deep, wide-field i'-band imaging from the Canada-France-Hawaii Telescope Legacy Survey, we compute size functions for 6000 disks from z=0.2 to z=1 and explore luminosity and number density evolution scenarios with an emphasis on the importance of selection effects on the interpretation of the data. We also compute the size function of a very large sample of disks from the Sloan Digital Sky Survey to use as a local (z ~ 0.1) comparison. CFHTLS size functions computed with the same fixed luminosity-size selection window at all redshifts exhibit evolution that appears to be best modelled by a pure number density evolution. The z = 0.3 size function is an excellent match to the z = 0.9 one if disks at the highest redshift are a factor of 2.5 more abundant than in the local universe. The SDSS size function would also match the z = 0.9 CFHTLS size function very well with a similar change in number density. On the other hand, the CFHTLS size functions computed with a varying luminosity-size selection window with redshift remain constant if the selection window is shifted by 1.0$-$1.5 mag towards fainter magnitudes with decreasing redshift. There is a weak dependence on disk scale length with smaller ($h \lesssim $ 4 kpc) disks requiring more luminosity evolution than larger ones. Given that changes in number density are primarily due to mergers and that current estimates of merger rates below z = 1 are low, luminosity evolution appears to be a more plausible scenario to explain the observations.

The Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS) has created a large homogeneous database of intermediate redshift (0.2 < z < 1.0) type Ia supernovae (SNe Ia). The SNLS team has shown that correlations exist between SN Ia rates, properties, and host galaxy star formation rates. The SNLS SN Ia database has now been combined with a photometric redshift galaxy catalog and an optical galaxy cluster catalog to investigate the possible influence of galaxy clustering on the SN Ia rate, over and above the expected effect due to the dependence of SFR on clustering through the morphology-density relation. We identify three cluster SNe Ia, plus three additional possible cluster SNe Ia, and find the SN Ia rate per unit mass in clusters at intermediate redshifts is consistent with the rate per unit mass in field early-type galaxies and the SN Ia cluster rate from low redshift cluster targeted surveys. We also find the number of SNe Ia in cluster environments to be within a factor of two of expectations from the two component SNIa rate model.

This paper describes the MegaPipe image processing pipeline at the Canadian Astronomical Data Centre. The pipeline combines multiple images from the MegaCam mosaic camera on CFHT and combines them into a single output image. MegaPipe takes as input detrended MegaCam images and does a careful astrometric and photometric calibration on them. The calibrated images are then resampled and combined into image stacks. The astrometric calibration of the output images is accurate to within 0.15 arcseconds relative to external reference frames and 0.04 arcseconds internally. The photometric calibration is good to within 0.03 magnitudes. The stacked images and catalogues derived from these images are available through the CADC website:

We present the quantitative rest-frame B morphological evolution and galaxy merger fractions at 0.2 < z < 1.2 as observed by the All-wavelength Extended Groth Strip International Survey (AEGIS). We use the Gini coefficent and M_20 to identify major mergers and classify galaxy morphology for a volume-limited sample of 3009 galaxies brighter than 0.4 L_B^*, assuming pure luminosity evolution of 1.3 M_B per unit redshift. We find that the merger fraction remains roughly constant at 10 +/- 2% for 0.2 < z < 1.2. The fraction of E/S0/Sa increases from 21+/- 3% at z ~ 1.1 to 44 +/- 9% at z ~ 0.3, while the fraction of Sb-Ir decreases from 64 +/- 6% at z ~ 1.1 to 47 +/- 9% at z ~ 0.3. The majority of z < 1.2 Spitzer MIPS 24 micron sources with L(IR) > 10^11 L_sun are disk galaxies, and only ~ 15% are classified as major merger candidates. Edge-on and dusty disk galaxies (Sb-Ir) are almost a third of the red sequence at z ~ 1.1, while E/S0/Sa makeup over 90% of the red sequence at z ~ 0.3. Approximately 2% of our full sample are red mergers. We conclude (1) the galaxy merger rate does not evolve strongly between 0.2 < z < 1.2; (2) the decrease in the volume-averaged star-formation rate density since z ~ 1 is a result of declining star-formation in disk galaxies rather than a disappearing population of major mergers; (3) the build-up of the red sequence at z < 1 can be explained by a doubling in the number of spheroidal galaxies since z ~ 1.2.

We present optical measurements of the faint end of the luminosity function in the core of the Coma cluster. Dwarf galaxies are detected down to a limiting magnitude of R approx. 25.75 in images taken with the Hubble Space Telescope. This represents the faintest determination of the Coma luminosity function to date. With the assumption that errors due to cosmic variance are small, evidence is found for a steep faint end slope with alpha < -2. Such a value is expected in theories in which reionization and other feedback processes are dependent on density.

We visually and quantitatively determine the host galaxy morphologies of 94 intermediate redshift (0.2 < z < 1.2) active galactic nuclei (AGN), selected using Chandra X-ray and Spitzer mid-infrared data in the Extended Groth Strip. Using recently developed morphology measures, the second-order moment of the brightest 20% of a galaxy's flux (M_{20}) and the Gini coefficient, we find that X-ray-selected AGN mostly reside in E/S0/Sa galaxies (53^{+11}_{-10}%), while IR-selected AGN show no clear preference for host morphology. X-ray-selected AGN hosts are members of close pairs more often than the field population by a factor of 3.3+/-1.4, but most of these pair members appear to be undisturbed early-type galaxies and do not tend to show direct evidence of gravitational perturbations or interactions. Thus, the activation mechanism for AGN activity remains unknown, even for pair members.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

We present the first cosmic shear measurements obtained from the T0001 release of the Canada-France-Hawaii Telescope Legacy Survey. The data set covers three uncorrelated patches (D1, D3 and D4) of one square degree each observed in u*, g', r', i' and z' bands, out to i'=25.5. The depth and the multicolored observations done in deep fields enable several data quality controls. The lensing signal is detected in both r' and i' bands and shows similar amplitude and slope in both filters. B-modes are found to be statistically zero at all scales. Using multi-color information, we derived a photometric redshift for each galaxy and separate the sample into medium and high-z galaxies. A stronger shear signal is detected from the high-z subsample than from the low-z subsample, as expected from weak lensing tomography. While further work is needed to model the effects of errors in the photometric redshifts, this results suggests that it will be possible to obtain constraints on the growth of dark matter fluctuations with lensing wide field surveys. The various quality tests and analysis discussed in this work demonstrate that MegaPrime/Megacam instrument produces excellent quality data. The combined Deep and Wide surveys give sigma_8= 0.89 pm 0.06 assuming the Peacock & Dodds non-linear scheme and sigma_8=0.86 pm 0.05 for the halo fitting model and Omega_m=0.3. We assumed a Cold Dark Matter model with flat geometry. Systematics, Hubble constant and redshift uncertainties have been marginalized over. Using only data from the Deep survey, the 1 sigma upper bound for w_0, the constant equation of state parameter is w_0 < -0.8.

We discuss the optical and X-ray spectral properties of the sources detected in a single 200ks Chandra pointing in the Groth-Westphal Strip region. Optical identifications and spectroscopic redshifts are primarily from the DEEP2 survey. This is complemented with deeper (r~26mag) multi-waveband data (ugriz) from the Canada France Hawaii Legacy Survey to estimate photometric redshifts and to optically identify sources fainter than the DEEP2 magnitude limit (R(AB)~24.5mag). We focus our study on the 2-10keV selected sample comprising 97 sources to the limit ~8e-16erg/s/cm2, this being the most complete in terms of optical identification rate (86%) and redshift determination fraction (63%; both spectroscopic and photometric). We first construct the redshift distribution of the sample which shows a peak at z~1. This is in broad agreement with models where less luminous AGNs evolve out to z~1 with powerful QSOs peaking at higher redshift, z~2. Evolution similar to that of broad-line QSOs applied to the entire AGN population (both type-I and II) does not fit the data. We also explore the observed N_H distribution of the sample and estimate a fraction of obscured AGN (N_H>1e22) of ~48%. This is found to be consistent with both a luminosity dependent intrinsic N_H distribution, where less luminous systems comprise a higher fraction of type-II AGNs, and models with a fixed ratio 2:1 between type-I and II AGNs. We further compare our results with those obtained in deeper and shallower surveys. We argue that a luminosity dependent parametrisation of the intrinsic N_H distribution is required to account for the fraction of obscured AGN observed in different samples over a wide range of fluxes.

We measure the build-up of the stellar mass of galaxies from z=6 to z=1. Using 15 band multicolour imaging data in the NICMOS Ultra Deep Field we derive photometric redshifts and masses for 796 galaxies down to H_AB=26.5. The derived evolution of the global stellar mass density of galaxies is consistent with previous star formation rate density measurements over the observed range of redshifts. Ongoing research in the CFHTLS Deep Fields confirms this result at lower redshifts. Further, if the sample is split by morphological type, a substantial increase is seen in the number of bulge dominated galaxies relative to disk-dominated galaxies since z=1.

We present preliminary results of a new quasar survey being undertaken with multi-colour optical imaging from the Canada-France-Hawaii Telescope. The current data consists of 3.83 sq. deg. of imaging in the i' and z' filters to a 10 sigma limit of z'<23.35. Near-infrared photometry of 24 candidate 5.7<z<6.4 quasars confirms them all to be low mass stars including two T dwarfs and four or five L dwarfs. Photometric estimates of the spectral type of the two T dwarfs are T3 and T6. We use the lack of high-redshift quasars in this survey volume to constrain the z=6 quasar luminosity function. For reasonable values of the break absolute magnitude M*_1450 and faint-end slope alpha, we determine that the bright-end slope beta>-3.2 at 95% confidence. We find that the comoving space-density of quasars brighter than M_1450=-23.5 declines by a factor >25 from z=2 to z=6, mirroring the decline observed for high-luminosity quasars. We consider the contribution of the quasar population to the ionizing photon density at z=6 and the implications for reionization. We show that the current constraints on the quasar population give an ionizing photon density <<30% that of the star-forming galaxy population. We conclude that active galactic nuclei make a negligible contribution to the reionization of hydrogen at z~6.

We study the relationships between galaxy total luminosity (M_g), morphology, color and environment as a function of redshift. We use a magnitude-limited sample of 65,624 galaxies in the redshift range 0<z<1.3 taken from one of the Canada-France-Hawaii Telescope Legacy Survey Deep Fields. We parametrize galaxy morphology according to the Sersic index n. Our n>2 number fractions at z=0.1 agree well with those based on SDSS galaxies. We find that the n>2 galaxy number fraction is constant with redshift in the field. However, for overdense environments this fraction is larger and increases towards lower redshifts, higher densities and higher luminosities. Rest-frame color-magnitude diagrams show that the color distribution is bimodal out to our redshift limit of z~1 with a prominent red-sequence of galaxies at 0.2<z<0.4 and a large blue-peak dominance at 0.8<z<1. We use this bimodality to define a red galaxy fraction (rest-frame u-g>1). For all environments, this fraction increases towards lower redshifts and higher luminosities. The red fraction within cluster-like regions changes 60% faster with redshift as compared to the field for M_g<-19.5. Using, for the first time, observations across many cluster-field interfaces distributed over a single, large volume, we trace the large-scale morphology-density relation and the Butcher-Oemler effect over a period of almost 8 Gyr.

We measure the build-up of the stellar mass of galaxies from z=6 to z=1. Using 15 band multicolor imaging data in the NICMOS Ultra Deep Field we derive photometric redshifts and masses for 796 galaxies down to H(AB)=26.5. The derived evolution of the global stellar mass density of galaxies is consistent with previous star formation rate density measurements over the observed range of redshifts. Beyond the observed range, maintaining consistency between the global stellar mass and the observed star formation rate suggests the epoch of galaxy formation was z=16.

(abridged) In this paper we describe in detail the reduction, preparation and reliability of the photometric catalogues which comprise the 1.2 deg^2 CFH12K-VIRMOS deep field. The survey reaches a limiting magnitude of BAB~26.5, VAB~26.2, RAB~25.9 IAB~25.0 and contains 90,729 extended sources in the magnitude range 18.0<IAB<24.0. We demonstrate our catalogues are free from systematic biases and are complete and reliable down these limits. We estimate that the upper limit on bin-to-bin systematic photometric errors for the I- limited sample is ~10% in this magnitude range. We estimate that 68% of the catalogues sources have absolute per co-ordinate astrometric uncertainties less than ~0.38" and ~0.32" (alpha,delta). Our internal (filter-to-filter) per co-ordinate astrometric uncertainties are 0.08" and 0.08" (alpha,delta). We quantify the completeness of our survey in the joint space defined by object total magnitude and peak surface brightness. Finally, we present numerous comparisons between our catalogues and published literature data: galaxy and star counts, galaxy and stellar colours, and the clustering of both point-like and extended populations. In all cases our measurements are in excellent agreement with literature data to IAB<24.0. This combination of depth and areal coverage makes this multi-colour catalogue a solid foundation to select galaxies for follow-up spectroscopy with VIMOS on the ESO-VLT and a unique database to study the formation and evolution of the faint galaxy population to z~1 and beyond.

This paper presents the CFH12K-VIRMOS survey: a deep B, V, R and I imaging survey in four fields totalling more than 17 deg^2, conducted with the 30x40 arcmin^2 field CFH-12K camera. The survey is intended to be a multi-purpose survey used for a variety of science goals, including surveys of very high redshift galaxies and weak lensing studies. Four high galactic latitude fields, each 2x2 deg^2, have been selected along the celestial equator: 0226-04, 1003+01, 1400+05, and 2217+00. The 16 deg^2 of the "wide" survey are covered with exposure times of 2h, 1.5h, 1h, 1h, while the 1.3x1 deg^2 area of the "deep" survey at the center of the 0226-04 field is covered with exposure times of 7h, 4.5h, 3h, 3h, in B,V,R and I respectively. The data is pipeline processed at the Terapix facility at the Institut d'Astrophysique de Paris to produce large mosaic images. The catalogs produced contain the positions, shape, total and aperture magnitudes for the 2.175 million objects. The depth measured (3sigma in a 3 arc-second aperture) is I_{AB}=24.8 in the ``Wide'' areas, and I_{AB}=25.3 in the deep area. Careful quality control has been applied on the data as described in joint papers. These catalogs are used to select targets for the VIRMOS-VLT Deep Survey, a large spectroscopic survey of the distant universe (Le F\`evre et al., 2003). First results from the CFH12K-VIRMOS survey have been published on weak lensing (e.g. van Waerbeke & Mellier 2003). Catalogs and images are available through the VIRMOS database environment under Oracle ({\tt http://www.oamp.fr/virmos}). They will be open for general use on July 1st, 2003.

Photometric redshifts have been calculated for the Hubble Deep Fields. The redshift distributions of the fields differ; there is a large excess of galaxies in the HDF-North in the redshift range 0.4<z<1.2. The difference is consistent with the presence of a weak cluster in the HDFN and is only slightly larger than the cosmic variance in other fields of similar depths and with models of large scale structure.

A tangential distortion of background source galaxies around foreground lens galaxies in the Hubble Deep Field is detected at the 99.3% confidence level. An important element of our analysis is the use of photometric redshifts to determine distances of lens and source galaxies and rest-frame B-band luminosities of the lens galaxies. The lens galaxy halos obey a Tully-Fisher relation between halo circular velocity and luminosity; the typical lens galaxy, at a redshift z = 0.6, has a circular velocity of 210 +/-40 km/s at M_B = -18.5, if q_0 = 0.5. Control tests, in which lens and source positions and source ellipticities are randomized, confirm the significance level of the detection quoted above. Furthermore, a marginal signal is also detected from an independent, fainter sample of source galaxies without photometric redshifts. Potential systematic effects, such as contamination by aligned satellite galaxies, the distortion of source shapes by the light of the foreground galaxies, PSF anisotropies, and contributions from mass distributed on the scale of galaxy groups are shown to be negligible. A comparison of our result with the local Tully-Fisher relation indicates that intermediate-redshift galaxies are fainter than local spirals by 1.0 +/- 0.6 B mag at a fixed circular velocity. This is consistent with some spectroscopic studies of the rotation curves of intermediate-redshift galaxies. This result suggests that the strong increase in the global luminosity density with redshift is dominated by evolution in the galaxy number density.

Results of a blind test of photometric redshift predictions against spectroscopic galaxy redshifts obtained in the Hubble Deep Field with the Keck Telescope are presented. The best photometric redshift schemes predict spectroscopic redshifts with a redshift accuracy of |Delta-z|<0.1 for more than 68 percent of sources and with |Delta-z|<0.3 for 100 percent, when single-feature spectroscopic redshifts are removed from consideration. This test shows that photometric redshift schemes work well at least when the photometric data are of high quality and when the sources are at moderate redshifts.

Photometric redshifts have been determined for the galaxies in the Hubble Deep Field. The resulting redshift distribution shows two peaks: one at z~0.6 and one at z~2.2. Luminosity functions derived from the redshifts show strong luminosity evolution as a function of redshift. This evolution is consistent with the Babul & Rees (1992) scenario wherein massive galaxies form stars at high redshift while star formation in dwarf galaxies is delayed until after z=1.