We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

We review the current status of the facility instrumentation for the Large Binocular Telescope (LBT). The LBT has 2x 8.4m primary mirrors on a single mount with an effective collecting area of 11.8m or 23m when interferometrically combined. The facility instruments are: 1) the Large Binocular Cameras (LBCs), each with a 23'x25' field of view (FOV). The blue and red optimized optical LBCs are mounted at the prime focus of the left and right primary mirrors, respectively. The filter suite of the two LBCs covers 0.3-1.1{\mu}m, including the new TiO (0.78{\mu}m) and CN (0.82{\mu}m) filters; 2) the Multi-Object Double Spectrograph (MODS), two identical optical spectrographs each mounted at a straight through f/15 Gregorian mount. MODS-1 & -2 can do imaging with Sloan filters and medium resolution (R~2000) spectroscopy, each with 24 interchangeable masks (multi-object or longslit) over a 6'x6' FOV. Each MODS is capable of blue (0.32-0.6{\mu}m) and red (0.5-1.05{\mu}m) wavelength only coverage or, using a dichroic, 0.32-1.05{\mu}m coverage; and 3) the two LBT Utility Camera in the Infrared instruments (LUCIs), each mounted at a bent-front Gregorian f/15 port. LUCI-1 & 2 are designed for seeing-limited (4'x4'FOV) and AO (0.5'x0.5' FOV) imaging & spectroscopy over 0.95-2.5{\mu}m with spectroscopic resolutions of R~400-11000, including 32 interchangeable cryogenically cooled masks. All facility instruments are on the LBT and, for the first time, have been on-sky for science. We also report on the first science use of "mixed-mode" (differently paired instruments). While both primary mirrors reside on a single fixed mount, they are capable of operating independently within a defined "co-pointing" limit. This provides users with the additional capability to independently dither each mirror or center observations on two different sets of spatial coordinates within this limit. (ABRIDGED)

We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wavelengths of four gravitationally lensed infrared (IR) luminous star-forming galaxies at redshift 1 < z < 3 from the LUCIFER instrument on the Large Binocular Telescope and the Infrared Spectrograph on Spitzer. The sample was selected to represent pure, actively star-forming systems, absent of active galactic nuclei. The large lensing magnifications result in high signal-to-noise spectra that can probe faint IR recombination lines, including Pa-alpha and Br-alpha at high redshifts. The sample was augmented by three lensed galaxies with similar suites of unpublished data and observations from the literature, resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction with H-alpha observations to probe the extinction, Av, of these systems, as well as testing star formation rate (SFR) indicators against the SFR measured by fitting spectral energy distributions to far-IR photometry. Our galaxies occupy a range of Av from ~0 to 5.9 mag, larger than previously known for a similar range of IR luminosities at these redshifts. Thus, estimates of SFR even at z ~ 2 must take careful count of extinction in the most IR luminous galaxies. We also measure extinction by comparing SFR estimates from optical emission lines with those from far-IR measurements. The comparison of results from these two independent methods indicates a large variety of dust distribution scenarios at 1 < z < 3. Without correcting for dust extinction, the H-alpha SFR indicator underestimates the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution scenario. Individual SFR estimates based on the 6.2 micron PAH emission line luminosity do not show a systematic discrepancy with extinction, although a considerable, ~0.2 dex scatter is observed.

We report the discovery of the first new ultra-cool brown dwarf found with the Wide-field Infrared Survey Explorer (WISE). The object's preliminary designation is WISEPC J045853.90+643451.9. Follow-up spectroscopy with the LUCIFER instrument on the Large Binocular Telescope indicates that it is a very late-type T dwarf with a spectral type approximately equal to T9. Fits to an IRTF/SpeX 0.8-2.5 micron spectrum to the model atmospheres of Marley and Saumon indicate an effective temperature of approximately 600 K as well as the presence of vertical mixing in its atmosphere. The new brown dwarf is easily detected by WISE, with a signal-to-noise ratio of ~36 at 4.6 microns. Current estimates place it at a distance of 6 to 10 pc. This object represents the first in what will likely be hundreds of nearby brown dwarfs found by WISE that will be suitable for follow up observations, including those with the James Webb Space Telescope. One of the two primary scientific goals of the WISE mission is to find the coolest, closest stars to our Sun; the discovery of this new brown dwarf proves that WISE is capable of fulfilling this objective.

Deep multicolour surveys are the main tool to explore the formation and evolution of the faint galaxies which are beyond the spectroscopic limit with the present technology. The photometric properties of these faint galaxies are usually compared with current renditions of semianalytical models to provide constraints on the fundamental physical processes involved in galaxy formation and evolution, namely the mass assembly and the star formation. Galaxy counts over large sky areas in the near-UV band are important because they are difficult to obtain given the low efficiency of near-UV instrumentation, even at 8m class telescopes. A large instrumental field of view helps in minimizing the biases due to the cosmic variance. We have obtained deep images in the 360nm U band provided by the blue channel of the Large Binocular Camera at the prime focus of the Large Binocular Telescope. We have derived over an area of ~0.4 sq. deg. the galaxy number counts down to U=27 in the Vega system (corresponding to U=27.86 in the AB system) at a completeness level of 30% reaching the faintest current limit for this wavelength and sky area. The shape of the galaxy counts in the U band can be described by a double power-law, the bright side being consistent with the shape of shallower surveys of comparable or greater areas. The slope bends over significantly at U>23.5 ensuring the convergence of the contribution by star forming galaxies to the EBL in the near-UV band to a value which is more than 70% of the most recent upper limits derived for this band. We have jointly compared our near-UV and K band counts collected from the literature with few selected hierarchical CDM models emphasizing critical issues in the physical description of the galaxy formation and evolution.

We present the characteristics and some early scientific results of the first instrument at the Large Binocular Telescope (LBT), the Large Binocular Camera (LBC). Each LBT telescope unit will be equipped with similar prime focus cameras. The blue channel is optimized for imaging in the UV-B bands and the red channel for imaging in the VRIz bands. The corrected field-of-view of each camera is approximately 30 arcminutes in diameter, and the chip area is equivalent to a 23x23 arcmin2 field. In this paper we also present the commissioning results of the blue channel. The scientific and technical performance of the blue channel was assessed by measurement of the astrometric distortion, flat fielding, ghosts, and photometric calibrations. These measurements were then used as input to a data reduction pipeline applied to science commissioning data. The measurements completed during commissioning show that the technical performance of the blue channel is in agreement with original expectations. Since the red camera is very similar to the blue one we expect similar performance from the commissioning that will be performed in the following months in binocular configuration. Using deep UV image, acquired during the commissioning of the blue camera, we derived faint UV galaxy-counts in a ~500 sq. arcmin. sky area to U(Vega)=26.5. These galaxy counts imply that the blue camera is the most powerful UV imager presently available and in the near future in terms of depth and extent of the field-of-view. We emphasize the potential of the blue camera to increase the robustness of the UGR multicolour selection of Lyman break galaxies at redshift z~3.

We present a 900 sec, wide-field U image of the inner region of the Andromeda galaxy obtained during the commissioning of the blue channel of the Large Binocular Camera mounted on the prime focus of the Large Binocular Telescope. Relative photometry and absolute astrometry of individual sources in the image was obtained along with morphological parameters aimed at discriminating between stars and extended sources, e.g. globular clusters. The image unveils the near-ultraviolet view of the inner ring of star formation recently discovered in the infrared by the Spitzer Space Telescope and shows in great detail the fine structure of the dust lanes associated with the galaxy inner spiral arms. The capabilities of the blue channel of the Large Binocular Camera at the Large Binocular Telescope (LBC-Blue) are probed by direct comparison with ultraviolet GALEX observations of the same region in M31. We discovered 6 new candidate stellar clusters in this high-background region of M31. We also recovered 62 bona-fide globulars and 62 previously known candidates from the Revised Bologna Catalogue of the M31 globular clusters, and firmly established the extended nature of 19 of them.

JHK photometry in the Mauna Kea Observatory (MKO) near-IR system is presented for 115 stars. Of these, 79 are UKIRT standards and 42 are LCO standards. The average brightness is 11.5 mag, with a range of 10 to 15. The average number of nights each star was observed is 4, and the average of the internal error of the final results is 0.011 mag. These JHK data agree with those reported by other groups to 0.02 mag. The measurements are used to derive transformations between the MKO JHK photometric system and the UKIRT, LCO and 2MASS systems. The 2MASS-MKO data scatter by 0.05 mag for redder stars: 2MASS-J includes H2O features in dwarfs and MKO-K includes CO features in giants. Transformations derived for stars whose spectra contain only weak features cannot give accurate transformations for objects with strong absorption features within a filter bandpasses. We find evidence of systematic effects at the 0.02 mag level in the photometry of stars with J<11 and H,K<10.5. This is due to an underestimate of the linearity correction for stars observed with the shortest exposure times; very accurate photometry of stars approaching the saturation limits of infrared detectors which are operated in double-read mode is difficult to obtain. Four stars in the sample, GSPC S705-D, FS 116 (B216-b7), FS 144 (Ser-EC84) and FS 32 (Feige 108), may be variable. 84 stars in the sample have 11< J< 15 and 10.5<H,K<15, are not suspected to be variable, and have magnitudes with an estimated error <0.027 mag; 79 of these have an error of <0.020 mag. These represent the first published high-accuracy JHK stellar photometry in the MKO photometric system; we recommend these objects be employed as primary standards for that system [abridged].

Near-IR photometry for 20 radio-loud z>3 quasars, 16 of which are radio- selected, are presented. These data sample the rest-frame optical/UV continuum, which is commonly interpreted as emission from an accretion disk. In a previous study, we compared the rest-frame optical/UV continuum shapes of 15 optically bright (V<17.5) z>3 quasars with those of 27 low redshift (z~0.1) ones that were matched to the high redshift sample in evolved luminosity (i.e. having luminosities ranging from 1-7 times the characteristic luminosity, L*, where L*~(1+z)^{~3}) to look for signs of evolution in the central engines. We found the continuum shapes at z~0.1 and z>3 similar, consistent with no significant change in the ratio mdot/M, where mdot is the accretion rate with respect to the Eddington rate and M is the black hole mass. This study expands our earlier high redshift sample to lower luminosity, away from extreme objects and towards a luminosity overlap with lower redshift samples. The distribution of rest-frame optical/UV continuum shapes for this fainter sample is broader, extending further to the red than that of the brighter z>3 one. Three quasars from this fainter sample, two radio-selected and one optically-selected, have optical continuum slopes alpha<-1 (F_{nu}~nu^{alpha}). The optically-selected one, LBQS0056+0125, appears to be reddened by dust along the line of sight or in the host galaxy, whereas the radio-selected ones, PKS2215+02 and TXS2358+189, could derive their red continua from the contribution of a relatively strong synchrotron component to the rest-frame optical. These objects may represent a bridge to a population of very red high redshift quasars to which ongoing or future near-IR, optical and deep X-ray surveys will be sensitive.

We present L'and M' photometry, obtained at UKIRT using the Mauna Kea Observatories Near-IR filter set, for 46 and 31 standard stars, respectively. The L' standards include 25 from the UKIRT in-house "Bright Standards" with magnitudes deriving from Elias et al. (1982) and observations at the IRTF in the early 1980s, and 21 fainter stars. The M' magnitudes derive from the results of Sinton & Tittemore (1984). We estimate the average external error to be 0.015 mag for the bright L' standards and 0.025 mag for the fainter L' standards, and 0.026 mag for the M' standards. The new results provide a network of homogeneously observed standards, and establish reference stars for the MKO system, in these bands. They also extend the available standards to magnitudes which should be faint enough to be accessible for observations with modern detectors on large and very large telescopes.

For 15 bright (V<17.5), high redshift (z>3) quasars, we have obtained IR spectra and photometry, and optical spectrophotometry and photometry, which we use to construct their rest-frame 1285-5100 AA spectral energy distributions (SEDs). High resolution spectroscopy for 7 and L' detections of 4 extend their SEDs shortwards of Ly_alpha or redwards of 7500 AA. The average luminosities within a set of narrow, line-free, windows are computed. Power law fits to those from 1285 to 5100 AA, but excluding the 2000-4000 AA region of the FeII+BaC `small bump', adequately characterize the continuum shapes of most of the objects and yield optical/UV spectral indices (aouv: fnu ~ nu^{aouv}). To look for signs of evolution, we compare the distribution of the aouv for the z>3 quasars with that for a set of 27 z~0.1 quasars which are matched to the high redshift(z) ones in evolved luminosity. The mean (median) aouv for the high and low z samples are -0.32(-0.29) and -0.38(-0.40), respectively. Neither the distributions nor means differ significantly. The distributions of aouv for both samples span a wide range however, with Delta~1, which could be understood as the result of both a diversity in the emitted continua themselves and in the amounts of intrinsic extinction undergone. A clear difference between the high and low z samples occurs in the region of `small bump'. The power law fit residuals for the low z sample show a systematic excess from 2200-3000 AA ; this feature is weak or absent in the high z sample. Further study is needed, but the `small bump' evolution could reflect differences in iron abundance or FeII energy source, or alternatively, an intrinsic turnover in the continuum itself which is present at low but not at high z. (Abridged)

Optical red spectra of a set of 18 bright barred spiral galaxies are presented. The study is aimed at determining the local kinematics, and the physical conditions of ionized gas in the compact nucleus (inside a diameter of 5") and in the circumnuclear regions (inside a diameter of 20"). Only 8 galaxies showed bright emission from their east and west side of the nucleus. The spectrum of each region was analized separately. In other 10 galaxies the line emission was so weak that we were only able to obtain an average spectrum of the central emission. No emission was detected in the remaining 8 galaxies. An estimate of the dynamical mass is presented based on the observed velocities in the circumnuclear regions. In NGC 4314 and NGC 6951, that show H_alpha emission distributed in circumnuclear ring structures, we determine the [NII]/H_alpha and [SII]/H_alpha ratios for the eastern and western regions of the rings. The velocity difference for the two sides is used to derive the rotation velocity of the gas around the compact nucleus. The ratio [NII]6583/H_alpha is a factor of 2 larger in the compact nucleus of NGC 6951 than in its western side. The electron gas densities have been estimated from the [SII] lines ratio.