Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s.

We present a response to the 2018 Request for Information (RFI) from the NITRD, NCO, NSF regarding the "Update to the 2016 National Artificial Intelligence Research and Development Strategic Plan." Through this document, we provide a response to the question of whether and how the National Artificial Intelligence Research and Development Strategic Plan (NAIRDSP) should be updated from the perspective of Fermilab, America's premier national laboratory for High Energy Physics (HEP). We believe the NAIRDSP should be extended in light of the rapid pace of development and innovation in the field of Artificial Intelligence (AI) since 2016, and present our recommendations below. AI has profoundly impacted many areas of human life, promising to dramatically reshape society --- e.g., economy, education, science --- in the coming years. We are still early in this process. It is critical to invest now in this technology to ensure it is safe and deployed ethically. Science and society both have a strong need for accuracy, efficiency, transparency, and accountability in algorithms, making investments in scientific AI particularly valuable. Thus far the US has been a leader in AI technologies, and we believe as a national Laboratory it is crucial to help maintain and extend this leadership. Moreover, investments in AI will be important for maintaining US leadership in the physical sciences.

The observation of gravitational waves from compact binary coalescences by LIGO and Virgo has begun a new era in astronomy. A critical challenge in making detections is determining whether loud transient features in the data are caused by gravitational waves or by instrumental or environmental sources. The citizen-science project \emph{Gravity Spy} has been demonstrated as an efficient infrastructure for classifying known types of noise transients (glitches) through a combination of data analysis performed by both citizen volunteers and machine learning. We present the next iteration of this project, using similarity indices to empower citizen scientists to create large data sets of unknown transients, which can then be used to facilitate supervised machine-learning characterization. This new evolution aims to alleviate a persistent challenge that plagues both citizen-science and instrumental detector work: the ability to build large samples of relatively rare events. Using two families of transient noise that appeared unexpectedly during LIGO's second observing run (O2), we demonstrate the impact that the similarity indices could have had on finding these new glitch types in the Gravity Spy program.

Large modern surveys require efficient review of data in order to find transient sources such as supernovae, and to distinguish such sources from artefacts and noise. Much effort has been put into the development of automatic algorithms, but surveys still rely on human review of targets. This paper presents an integrated system for the identification of supernovae in data from Pan-STARRS1, combining classifications from volunteers participating in a citizen science project with those from a convolutional neural network. The unique aspect of this work is the deployment, in combination, of both human and machine classifications for near real-time discovery in an astronomical project. We show that the combination of the two methods outperforms either one used individually. This result has important implications for the future development of transient searches, especially in the era of LSST and other large-throughput surveys.

The Wide-field Infrared Survey Explorer (WISE) is a powerful tool for finding nearby brown dwarfs and searching for new planets in the outer solar system, especially with the incorporation of NEOWISE and NEOWISE-Reactivation data. So far, searches for brown dwarfs in WISE data have yet to take advantage of the full depth of the WISE images. To efficiently search this unexplored space via visual inspection, we have launched a new citizen science project, called "Backyard Worlds: Planet 9," which asks volunteers to examine short animations composed of difference images constructed from time-resolved WISE coadds. We report the discovery of the first new substellar object found by this project, WISEA J110125.95+540052.8, a T5.5 brown dwarf located approximately 34 pc from the Sun with a total proper motion of $\sim$0.7 as yr$^{-1}$. WISEA J110125.95+540052.8 has a WISE $W2$ magnitude of $W2=15.37 \pm 0.09$, this discovery demonstrates the ability of citizen scientists to identify moving objects via visual inspection that are 0.9 magnitudes fainter than the $W2$ single-exposure sensitivity, a threshold that has limited prior motion-based brown dwarf searches with WISE.

(abridged for arXiv) With the first direct detection of gravitational waves, the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) has initiated a new field of astronomy by providing an alternate means of sensing the universe. The extreme sensitivity required to make such detections is achieved through exquisite isolation of all sensitive components of LIGO from non-gravitational-wave disturbances. Nonetheless, LIGO is still susceptible to a variety of instrumental and environmental sources of noise that contaminate the data. Of particular concern are noise features known as glitches, which are transient and non-Gaussian in their nature, and occur at a high enough rate so that accidental coincidence between the two LIGO detectors is non-negligible. In this paper we describe an innovative project that combines crowdsourcing with machine learning to aid in the challenging task of categorizing all of the glitches recorded by the LIGO detectors. Through the Zooniverse platform, we engage and recruit volunteers from the public to categorize images of glitches into pre-identified morphological classes and to discover new classes that appear as the detectors evolve. In addition, machine learning algorithms are used to categorize images after being trained on human-classified examples of the morphological classes. Leveraging the strengths of both classification methods, we create a combined method with the aim of improving the efficiency and accuracy of each individual classifier. The resulting classification and characterization should help LIGO scientists to identify causes of glitches and subsequently eliminate them from the data or the detector entirely, thereby improving the rate and accuracy of gravitational-wave observations. We demonstrate these methods using a small subset of data from LIGO's first observing run.

We present the results of an optical spectroscopic survey of a sample of 40 candidate obscured quasars identified on the basis of their mid-infrared emission detected by the Wide-Field Infrared Survey Explorer (WISE). Optical spectra for this survey were obtained using the Robert Stobie Spectrograph (RSS) on the Southern African Large Telescope (SALT). Our sample was selected with WISE colors characteristic of AGNs, as well as red optical to mid-IR colors indicating that the optical/UV AGN continuum is obscured by dust. We obtain secure redshifts for the majority of the objects that comprise our sample (35/40), and find that sources that are bright in the WISE W4 (22$\mu$m) band are typically at moderate redshift (<z> = 0.35$) while sources fainter in W4 are at higher redshifts (<z> = 0.73$). The majority of the sources have narrow emission lines, with optical colors and emission line ratios of our WISE-selected sources that are consistent with the locus of AGN on the rest-frame $g-z$ color vs. [NeIII]$\lambda$3869 / [OII]$\lambda\lambda$3726+3729 line ratio diagnostic diagram. We also use empirical AGN and galaxy templates to model the spectral energy distributions (SEDs) for the objects in our sample, and find that while there is significant variation in the observed SEDs for these objects, the majority require a strong AGN component. Finally, we use the results from our analysis of the optical spectra and the SEDs to compare our selection criteria to alternate criteria presented in the literature. These results verify the efficacy of selecting luminous obscured AGNs based on their WISE colors.

Recent cosmological modeling efforts have shown that a local underdensity on scales of a few hundred Mpc (out to z ~ 0.1), could produce the apparent acceleration of the expansion of the universe observed via type Ia supernovae. Several studies of galaxy counts in the near-infrared (NIR) have found that the local universe appears under-dense by ~25-50% compared with regions a few hundred Mpc distant. Galaxy counts at low redshifts sample primarily L ~ L* galaxies. Thus, if the local universe is under-dense, then the normalization of the NIR galaxy luminosity function (LF) at z>0.1 should be higher than that measured for z<0.1. Here we present a highly complete (> 90%) spectroscopic sample of 1436 galaxies selected in the H-band to study the normalization of the NIR LF at 0.1<z<0.3 and address the question of whether or not we reside in a large local underdensity. We find that for the combination of our six fields, the product phi* L* at 0.1 < z < 0.3 is ~ 30% higher than that measured at lower redshifts. While our statistical errors in this measurement are on the ~10% level, we find the systematics due to cosmic variance may be larger still. We investigate the effects of cosmic variance on our measurement using the COSMOS cone mock catalogs from the Millennium simulation and recent empirical estimates. We find that our survey is subject to systematic uncertainties due to cosmic variance at the 15% level ($1 sigma), representing an improvement by a factor of ~ 2 over previous studies in this redshift range. We conclude that observations cannot yet rule out the possibility that the local universe is under-dense at z<0.1.

The Baldwin, Phillips, and Terlevich emission-line ratio diagnostic ([OIII]/H{\beta} versus [NII]/H{\alpha}, hereafter BPT diagram) efficiently separates galaxies whose signal is dominated by star formation <BPT-SF> from those dominated by AGN activity (BPT-AGN). Yet this BPT diagram is limited to z < 0.5, the redshift at which [NII]{\lambda}6584 leaves the optical spectral window. Using the Sloan Digital Sky Survey (SDSS), we construct a new diagnostic, or TBT diagram, that is based on rest-frame g-z color, [NeIII]{\lambda}3869, and [OII]{\lambda}{\lambda}3726 + 3729 and can be used for galaxies out to z < 1.4. The TBT diagram identifies 98.7% of the SDSS BPT-AGN as TBT-AGN and 97% of the SDSS BPT-SF as TBT-SF. Furthermore, it identifies 97% of the OPTX Chandra X-ray selected AGNs as TBT-AGN. This is in contrast to the BPT diagram, which misidentifies 20% of X-ray selected AGNs as BPT-SF. We use the GOODS-N and Lockman Hole galaxy samples, with their accompanying deep Chandra imaging, to perform X-ray and infrared stacking analyses to further validate our TBT-AGN and TBT-SF selections; that is, we verify the dominance of AGN activity in the former and star formation activity in the latter. Finally, we address the inclusion of the majority of the BPT-comp (sources lying between the BPT-SF and BPT-AGN regimes) in our TBT-AGN regime. We find that the stacked BPT-comp source is X-ray hard (<{\Gamma}eff> = 1.0 +/-0.4) and has a high X-ray luminosity to total infrared luminosity ratio. This suggests that, on average, the X-ray signal in BPT-comp is dominated by obscured or low accretion rate AGN activity rather than by star formation, supporting their inclusion in the TBT-AGN regime.

We present an atlas of 88 z~5.7 and 30 z~6.5 Ly alpha emitters obtained from a wide-field narrowband survey. We combined deep narrowband imaging in 120A bandpass filters centered at 8150A and 9140A with deep BVRIz broadband imaging to select high-redshift galaxy candidates over an area of 4180 square arcmin. The goal was to obtain a uniform selection of comparable depth over the 7 targeted fields in the two filters. For the GOODS-N region of the HDF-N field, we also selected candidates using a 120A filter centered at 9210A. We made spectroscopic observations with Keck DEIMOS of nearly all the candidates to obtain the final sample of Ly alpha emitters. At the 3.3A resolution of the DEIMOS observations the asymmetric profile for Ly alpha emission with its steep blue fall-off can be clearly seen in the spectra of nearly all the galaxies. We show that the spectral profiles are surprisingly similar for many of the galaxies and that the composite spectral profiles are nearly identical at z=5.7 and z=6.5. We analyze the distributions of line widths and Ly alpha equivalent widths and find that the lines are marginally narrower at the higher redshift, with median values of 0.77A at z=6.5 and 0.92A at z=5.7. The line widths have a dependence on the Ly alpha luminosity of the form L(L alpha)^(0.3). We compare the surface densities and the luminosity functions at the two redshifts and find that there is a multiplicative factor of 2 decrease in the number density of bright Ly alpha emitters from z=5.7 to z=6.5, while the characteristic luminosity is unchanged.

We compare optical and hard X-ray identifications of AGNs using a uniformly selected (above a flux limit of f_2-8 keV = 3.5e-15 erg/cm2/s) and highly optically spectroscopically complete ( > 80% for f_2-8 keV > 1e-14 erg/cm2/s and > 60% below) 2-8 keV sample observed in three Chandra fields (CLANS, CLASXS, and the CDF-N). We find that empirical emission-line ratio diagnostic diagrams misidentify 20-50% of the X-ray selected AGNs that can be put on these diagrams as star formers, depending on which division is used. We confirm that there is a large (2 orders in magnitude) dispersion in the log ratio of the [OIII]5007A to hard X-ray luminosities for the non-broad line AGNs, even after applying reddening corrections to the [OIII] luminosities. We find that the dispersion is similar for the broad-line AGNs, where there is not expected to be much X-ray absorption from an obscuring torus around the AGN nor much obscuration from the galaxy along the line-of-sight if the AGN is aligned with the galaxy. We postulate that the X-ray selected AGNs that are misidentified by the diagnostic diagrams have low [OIII] luminosities due to the complexity of the structure of the narrow-line region, which causes many ionizing photons from the AGN not to be absorbed. This would mean that the [OIII] luminosity can only be used to predict the X-ray luminosity to within a factor of ~3 (one sigma). Despite selection effects, we show that the shapes and normalizations of the [OIII] and transformed hard X-ray luminosity functions show reasonable agreement, suggesting that the [OIII] samples are not finding substantially more AGNs at low redshifts than hard X-ray samples.

We present a deep, wide-field near-infrared (NIR) survey over five widely separated fields at high Galactic latitude covering a total of ~ 3 deg^2 in J, H, and Ks. The deepest areas of the data (~ 0.25 deg^2) extend to a 5 sigma limiting magnitude of JHKs > 24 in the AB magnitude system. Although depth and area vary from field to field, the overall depth and large area of this dataset make it one of the deepest wide-field NIR imaging surveys to date. This paper discusses the observations, data reduction, and bright galaxy counts in these fields. We compare the slope of the bright galaxy counts with the Two Micron All Sky Survey (2MASS) and other counts from the literature and explore the relationship between slope and supergalactic latitude. The slope near the supergalactic equator is sub- Euclidean on average pointing to the possibility of a decreasing average space density of galaxies by ~ 10-15% over scales of ~ 250-350 Mpc. On the contrary, the slope at high supergalactic latitudes is strongly super-Euclidean on average suggesting an increase in the space density of galaxies as one moves from the voids just above and below the supergalactic plane out to distances of ~ 250-350 Mpc. These results suggest that local large scale structure could be responsible for large discrepancies in the measured slope between different studies in the past. In addition, the local universe away from the supergalactic plane appears to be underdense by ~ 25-100% relative to the space densities of a few hundred megaparsecs distant. Subject headings: cosmology: observations and large scale structure of universe-galaxies: fundamental parameters (counts)-infrared: galaxies

We compare the optical spectral types with the X-ray spectral properties for a uniformly selected (sources with fluxes greater than the 3 sigma level and above a flux limit of f_2-8 keV > 3.5x10^-15 erg/cm2/s), highly spectroscopically complete (>80% for f_2-8 keV > 10^-14 erg/cm2/s and >60% below) 2-8 keV X-ray sample observed in three Chandra fields (CLANS, CLASXS, and the CDF-N) that cover ~1.2 deg^2. For our sample of 645 spectroscopically observed sources, we confirm that there is significant overlap of the X-ray spectral properties, as determined by the effective photon indices, Geff, obtained from the ratios of the 0.5-2 keV to 2-8 keV counts, for the different optical spectral types. For example, of the broad-line AGNs (non-broad-line AGNs), 20% +/- 3% (33% +/- 4%) have Geff<1.2 (Geff > 1.2). Thus, one cannot use the X-ray spectral classifications and the optical spectral classifications equivalently. Since it is not understood how X-ray and optical classifications relate to the obscuration of the central engine, we strongly advise against a mixed classification scheme, as it can only complicate the interpretation of X-ray AGN samples. We confirm the dependence of optical spectral type on X-ray luminosity, and for z<1, we find a similar luminosity dependence of Geff. However, this dependence breaks down at higher redshifts due to the highly redshift-dependent nature of Geff. We therefore also caution that any classification scheme which depends on Geff is likely to suffer from serious redshift bias.

We use the largest, most uniform, and most spectroscopically complete to faint X-ray flux limits Chandra sample to date to construct hard 2-8 keV rest-frame X-ray luminosity functions (HXLFs) of spectroscopically identified active galactic nuclei (AGNs) to z~5. In addition, we use a new 2-8 keV local sample selected by the very hard (14-195 keV) SWIFT 9-month Burst Alert Telescope (BAT) survey to construct the local 2-8 keV HXLF. We do maximum likelihood fits of the combined distant plus local sample (as well as of the distant sample alone) over the redshift intervals 0<z<1.2, 0<z<3, and 0<z<5 using a variety of analytic forms, which we compare with the HXLFs. We recommend using our luminosity dependent density evolution (LDDE) model fits of the combined distant plus local sample over 0<z<3 for all the spectroscopically identified sources and for the broad-line AGNs.

We use a wide-field (0.9 square degree) X-ray sample with optical and GALEX ultraviolet observations to measure the contribution of Active Galactic Nuclei (AGNs) to the ionizing flux as a function of redshift. Our analysis shows that the AGN contribution to the metagalactic ionizing background peaks around z=2. The measured values of the ionizing background from the AGNs are lower than previous estimates and confirm that ionization from AGNs is insufficient to maintain the observed ionization of the intergalactic medium (IGM) at z>3. We show that only sources with broad lines in their optical spectra have detectable ionizing flux and that the ionizing flux seen in an AGN is not correlated with its X-ray color. We also use the GALEX observations of the GOODS-N region to place a 2-sigma upper limit of 0.008 on the average ionization fraction fnu(700 A)/fnu(1500 A) for 626 UV selected galaxies in the redshift range z=0.9-1.4. We then use this limit to estimate an upper bound to the galaxy contribution in the redshift range z=0-5. If the z~1.15 ionization fraction is appropriate for higher redshift galaxies, then contributions from the galaxy population are also too low to account for the IGM ionization at the highest redshifts (z>4).

We present the redshift catalogs for the X-ray sources detected in the Chandra Deep Field North (CDF-N), the Chandra Large Area Synoptic X-ray Survey (CLASXS), and the Chandra Lockman Area North Survey (CLANS). The catalogs for the CDF-N and CLASXS fields include redshifts from previous work, while the redshifts for the CLANS field are all new. For fluxes above 10^-14 ergs cm^-2 s^-1 (2-8 keV) we have redshifts for 76% of the sources. We extend the redshift information for the full sample using photometric redshifts. The goal of the OPTX Project is to use these three surveys, which are among the most spectroscopically complete surveys to date, to analyze the effect of spectral type on the shape and evolution of the X-ray luminosity functions and to compare the optical spectral types with the X-ray spectral properties. We also present the CLANS X-ray catalog. The nine ACIS-I fields cover a solid angle of ~0.6 square degrees and reach fluxes of 7x10^-16 ergs cm^-2 s^-1 (0.5-2 keV) and 3.5x10^-15 ergs cm^-2 s^-1 (2-8 keV). We find a total of 761 X-ray point sources. Additionally, we present the optical and infrared photometric catalog for the CLANS X-ray sources, as well as updated optical and infrared photometric catalogs for the X-ray sources in the CLASXS and CDF-N fields. The CLANS and CLASXS surveys bridge the gap between the ultradeep pencil-beam surveys, such as the CDFs, and the shallower, very large-area surveys. As a result, they probe the X-ray sources that contribute the bulk of the 2-8 keV X-ray background and cover the flux range of the observed break in the logN-logS distribution. We construct differential number counts for each individual field and for the full sample.