In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$\alpha$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

Quasars at early redshifts ($z > 6$) with companion galaxies offer unique insights into the growth and evolution of the first supermassive black holes. Here, we report on a 150 ks Chandra observation of PSO J308.0416$-$21.2339, a $z=6.23$ quasar with a merging companion galaxy identified in [C II] and rest-frame UV emission. With $72.3^{+9.6}_{-8.6}$ net counts, we find that PSO J308.0416$-$21.2339 is powerful ($L_X = 2.31^{+1.14}_{-0.76} \times 10^{45}\ \textrm{erg}\,\textrm{s}^{-1}\,\textrm{cm}^{-2}$ in rest-frame $2.0-10.0$ keV) yet soft (spectral power-law index $\Gamma=2.39^{+0.37}_{-0.36}$ and optical-to-X-ray slope $\alpha_{\rm OX} = -1.41 \pm 0.11$). In addition, we detect three hard-energy photons $2.''0$ to the west of the main quasar, cospatial with the brightest UV emission of the merging companion. As no soft-energy photons are detected in the same area, this is potentially indicative of a highly-obscured source. With conservative assumptions, and accounting for both background fluctuations and the extended wings of the quasar's emission, these photons only have a probability $P=0.021$ of happening by chance. If confirmed by deeper observations, this system is the first high redshift quasar and companion individually detected in X-rays and is likely a dual AGN.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

We use MUSE/VLT to conduct a survey of $z\sim3$ physical quasar pairs at close separation with a fast observation strategy. Our aim is twofold: (i) explore the Ly$\alpha$ glow around the faint-end of the quasar population; (ii) take advantage of the combined illumination of a quasar pair to unveil large-scale intergalactic structures extending between the two quasars. Here, we report the results for a quasar pair ($z=3.020,3.008$; $i=21.84,22.15$), separated by 11.6 arcsec (or 89 projected kpc). MUSE reveals filamentary Ly$\alpha$ structures extending between the two quasars with an average surface brightness of SB$_{\rm Ly\alpha}=1.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Photoionization models of the constraints in the Ly$\alpha$, HeII, and CIV line emissions show that the emitting structures are intergalactic bridges with an extent between $\sim89$ and up to $\sim600$ kpc. Our models rule out the possibility that the structure extends for $\sim 2.9$ Mpc, i.e., the separation inferred from the uncertain systemic redshift difference of the quasars if the difference was only due to the Hubble flow. At the current spatial resolution and surface brightness limit, the average projected width of an individual bridge is about 35 kpc. We also detect a strong absorption in HI, NV, and CIV along the background sight-line at higher $z$, which we interpret as due to at least two components of cool, metal enriched, and relatively ionized CGM or IGM surrounding the quasar pair. Two additional HI absorbers are detected along both quasar sight-lines at $\sim -900$ and $-2800$ km s$^{-1}$ from the system, with the latter having associated CIV absorption only along the foreground quasar sight-line. The absence of galaxies in the MUSE field of view at the redshifts of these two absorbers suggests that they trace large-scale structures or expanding shells in front of the quasar pair.

We present ALMA 0.28"x0.20" (1.4kpc x 1.0kpc) resolution observations of the [CII] 158um line and underlying dust continuum in the host galaxy of the most distant quasar currently known, ULAS J1342+0928 at z=7.5413. Both the [CII] and continuum are detected and spatially resolved. The [CII] line is ~1.5 times more extended than the continuum emission, showing an elongated and complex structure of approximately 3.2kpc x 6.4kpc. Two separate peaks are clearly seen (>6sigma each) in three 100 km/s width [CII] channel maps. The velocity field of the [CII] gas does not show evidence of a coherent rotation field but rather chaotic motion, reminiscent of an on-going merger.

We present a survey for metal absorption systems traced by neutral oxygen over $3.2 < z < 6.5$. Our survey uses Keck/ESI and VLT/X-Shooter spectra of 199 QSOs with redshifts up to 6.6. In total we detect 74 O I absorbers, of which 57 are separated from the background QSO by more than 5000 km/s. We use a maximum likelihood approach to fit the distribution of O I $\lambda$1302 equivalent widths in bins of redshift, and from this determine the evolution in number density of absorbers with $W_{1302} > 0.05$ angstroms. We find that the number density does not monotonically increase with decreasing redshift, as would naively be expected from the buildup of metal-enriched circumgalactic gas with time. The number density over $4.9 < z < 5.7$ is a factor of 1.7--4.1 lower (68% confidence) than over $5.7 < z < 6.5$, with a lower value at $z < 5.7$ favored with 99% confidence. This decrease suggests that the fraction of metals in a low-ionization phase is larger at $z \sim 6$ than at lower redshifts. Absorption from highly ionized metals traced by C IV is also weaker in higher-redshift O I systems, supporting this picture. The evolution of O I absorbers implies that metal-enriched circumgalactic gas at $z \sim 6$ is undergoing an ionization transition driven by a strengthening ultraviolet background. This in turn suggests that the reionization of the diffuse intergalactic medium may still be ongoing at or only recently ended by this epoch.

We report on ~0.35"(~2 kpc) resolution observations of the [CII] and dust continuum emission from five z>6 quasar host-companion galaxy pairs obtained with the Atacama Large Millimeter/submillimeter Array. The [CII] emission is resolved in all galaxies, with physical extents of 3.2-5.4 kpc. The dust continuum is on-average 40% more compact, which results in larger [CII] deficits in the center of the galaxies. However, the measured [CII] deficits are fully consistent with those found at lower redshifts. Four of the galaxies show [CII] velocity fields that are consistent with ordered rotation, while the remaining six galaxies show no clear velocity gradient. All galaxies have high (~80-200 km/s) velocity dispersions, consistent with the interpretation that the interstellar medium (ISM) of these high redshift galaxies is turbulent. By fitting the galaxies with kinematic models, we estimate the dynamical mass of these systems, which range between (0.3 -> 5.4) x 1E10 Msun. For the three closest separation galaxy pairs, we observe dust and [CII] emission from gas in between and surrounding the galaxies, which is an indication that tidal interactions are disturbing the gas in these systems. Although gas exchange in these tidal interactions could power luminous quasars, the existence of quasars in host galaxies without nearby companions suggests that tidal interactions are not the only viable method for fueling their active centers. These observations corroborate the assertion that accreting supermassive black holes do not substantially contribute to the [CII] and dust continuum emission of the quasar host galaxies, and showcase the diverse ISM properties of galaxies when the universe was less than one billion years old.

Massive, quiescent galaxies are already observed at redshift $z\sim4$, i.e. $\sim$1.5 Gyr after the Big Bang. Current models predict them to be formed via massive, gas-rich mergers at $z>6$. Recent ALMA observations of the cool gas and dust in $z\gtrsim$6 quasars have discovered [CII]- and far infrared-bright galaxies adjacent to several quasars. In this work, we present sensitive imaging and spectroscopic follow-up observations, with HST/WFC3, Spitzer/IRAC, VLT/MUSE, Magellan/FIRE and LBT/LUCI-MODS, of ALMA-detected, dust-rich companion galaxies of four quasars at $z\gtrsim 6$, specifically acquired to probe their stellar content and unobscured star formation rate. Three companion galaxies do not show significant emission in the observed optical/IR wavelength range. The photometric limits suggest that these galaxies are highly dust-enshrouded, with unobscured star formation rates SFR$_{\rm UV}<$few M$_{\odot}$ yr$^{-1}$, and a stellar content of $M_{*}<$10$^{10}$ M$_{\odot}$ yr$^{-1}$. However, the companion to PJ167-13 shows bright rest-frame UV emission (F140W AB = 25.48). Its SED resembles that of a star-forming galaxy with a total SFR$\sim$50 M$_{\odot}$ yr$^{-1}$ and $M_{*}\sim 9 \times 10^{9}$ M$_{\odot}$. All the companion sources are consistent with residing on the galaxy main sequence at $z\sim$6. Additional, deeper data from future facilities, such as JWST, are needed in order to characterize these gas-rich sources in the first Gyr of cosmic history.

We use ALMA observations of the host galaxy of the quasar ULAS-J1342+0928 at z=7.54 to study the dust continuum and far infrared lines emitted from its interstellar medium. The Rayleigh-Jeans tail of the dust continuum is well sampled with eight different spectral setups, and from a modified black body fit we obtain an emissivity coefficient of beta=1.85+-0.3. Assuming a standard dust temperature of 47 K we derive a dust mass of Mdust=0.35x10^8 Msol and a star formation rate of 150+-30 Msol/yr. We have >4sigma detections of the [CII]_158, [OIII]_88 and [NII]_122 atomic fine structure lines and limits on the [CI]_369, [OI]_146 and [NII]_205 emission. We also report multiple limits of CO rotational lines with Jup >= 7, as well as a tentative 3.3sigma detection of the stack of four CO lines (Jup=11, 10, 8, 7). We find line deficits that are in agreement with local ultra luminous infrared galaxies. Comparison of the [NII]_122 and [CII]_158 lines indicates that the [CII]_158 emission arises predominantly from the neutral medium, and we estimate that the photo-disassociation regions in J1342+0928 have densities <~5x10^4cm^-3. The data suggest that ~16% of hydrogen is in ionized form and that the HII regions have high electron densities of n_e>180 cm^-3. Our observations favor a low gas-to-dust ratio of <100 and a metallicity of the interstellar medium comparable to the Solar value. All the measurements presented here suggest that the host galaxy of J1342+0928 is highly enriched in metal and dust, despite being observed just 680 Myr after the Big Bang.

We present deep MUSE observations of five quasars within the first Gyr of the Universe ($z\gtrsim6$), four of which display extended Ly$\alpha$ halos. After PSF-subtraction, we reveal halos surrounding two quasars for the first time, as well as confirming the presence of two more halos for which tentative detections exist in long-slit spectroscopic observations and narrow-band imaging. The four Ly$\alpha$ halos presented here are diverse in morphology and size, they each display spatial asymmetry, and none are centred on the position of the quasar. Spectra of the diffuse halos demonstrate that none are dramatically offset in velocity from the systemic redshift of the quasars ($\Delta$ v $< 200$ kms$^{-1}$), however each halo shows a broad Ly$\alpha$ line, with a velocity width of order $\sim1000$ kms$^{-1}$. Total Ly$\alpha$ luminosities range between $\sim$ $2 \times 10^{43}$ erg s$^{-1}$ and $\sim$ $2 \times 10^{44}$ erg s$^{-1}$, reaching maximum radial extents of $13 - 30$ pkpc from the quasar positions. We find larger sizes and higher Ly$\alpha$ luminosities than previous literature results at this redshift, but find no correlation between the quasar properties and the Ly$\alpha$ halo, suggesting that the detected emission is most closely related to the physical properties of the circum-galactic medium

We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $\mu$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

Luminous high-redshift quasars can be used to probe of the intergalactic medium (IGM) in the early universe because their UV light is absorbed by the neutral hydrogen along the line of sight. They help us to measure the neutral hydrogen fraction of the high-z universe, shedding light on the end of reionization epoch. In this paper, we present a discovery of a new quasar (PSO J006.1240+39.2219) at redshift $z=6.61\pm0.02$ from Panoramic Survey Telescope & Rapid Response System 1. Including this quasar, there are nine quasars above $z>6.5$ up to date. The estimated continuum brightness is $M_\text{1450}$=$-25.96\pm0.08$. PSO J006.1240+39.2219 has a strong Ly~$\alpha$ emission compared with typical low-redshift quasars, but the measured near-zone region size is $R_\text{NZ}=3.2\pm1.1$ proper megaparsecs, which is consistent with other quasars at z$\sim$6.

We present ALMA band 8 observations of the [OIII] 88um line and the underlying thermal infrared continuum emission in the z=6.08 quasar CFHQS J2100-1715 and its dust-obscured starburst companion galaxy (projected distance: ~60 kpc). Each galaxy hosts dust-obscured star formation at rates > 100 M_sun/yr, but only the quasar shows evidence for an accreting 10^9 M_sun black hole. Therefore we can compare the properties of the interstellar medium in distinct galactic environments in two physically associated objects, ~1 Gyr after the Big Bang. Bright [OIII] 88um emission from ionized gas is detected in both systems; the positions and line-widths are consistent with earlier [CII] measurements, indicating that both lines trace the same gravitational potential on galactic scales. The [OIII] 88um/FIR luminosity ratios in both sources fall in the upper range observed in local luminous infrared galaxies of similar dust temperature, although the ratio of the quasar is smaller than in the companion. This suggests that gas ionization by the quasar (expected to lead to strong optical [OIII] 5008A emission) does not dominantly determine the quasar's FIR [OIII] 88um luminosity. Both the inferred number of photons needed for the creation of O++ and the typical line ratios can be accounted for without invoking extreme (top-heavy) stellar initial mass functions in the starbursts of both sources.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

We present Atacama Large Millimeter Array 1mm observations of the rest-frame far-infrared (FIR) dust continuum in 27 quasars at redshifts 6.0 < z < 6.7. We detect FIR emission at >3sigma in all quasar host galaxies with flux densities at ~1900GHz in the rest-frame of 0.12 < S_rest,1900GHz < 5.9mJy, with a median (mean) flux density of 0.88mJy (1.59mJy). The implied FIR luminosities range from L_FIR = (0.27-13)x10^12 L_sun, with 74% of our quasar hosts having L_FIR > 10^12 L_sun. The estimated dust masses are M_dust = 10^7-10^9 M_sun. If the dust is heated only by star formation, then the star formation rates in the quasar host galaxies are between 50 and 2700 M_sun/yr. In the framework of the host galaxy-black hole coevolution model a correlation between ongoing black hole growth and star formation in the quasar host galaxy would be expected. However, combined with results from the literature to create a luminosity-limited quasar sample, we do not find a strong correlation between quasar UV luminosity (a proxy for ongoing black hole growth) and FIR luminosity (star formation in the host galaxy). The absence of such a correlation in our data does not necessarily rule out the coevolution model, and could be due to a variety of effects (including different timescales for black hole accretion and FIR emission).

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

Bright high-redshift quasars ($z>6$), hosting supermassive black holes ($\rm M_{BH} > 10^8 M_{\odot}$), are expected to reside in massive host galaxies embedded within some of the earliest and most massive galaxy overdensities. We analyze 1.2\,mm ALMA dust continuum maps of 35 bright quasars at $6<z<7$ and search the primary beam for excess dust continuum emission from sources with L$_{\rm IR} \gtrsim{ 10^{12}}$\,L$_{\odot}$ as evidence for early protoclusters. We compare the detection rates of continuum sources at $\geq5\sigma$ significance in the fields surrounding the quasars (A$_{\rm eff}$ = 4.3\,arcmin$^2$) with millimeter number counts in blank field surveys. We discover 15 millimeter sources in the fields excluding the quasars themselves, corresponding to an overdensity $\delta_{\rm gal} \equiv (N_{\rm gal} - N_{\rm exp})/N_{\rm exp} = -0.07\pm0.56$, consistent with no detected overdensity of dusty galaxies within 140 physical kpc of the quasars. However, the apparent lack of continuum overdensity does not negate the hypothesis that quasars live in overdense environments, as evidenced by strong [CII] overdensities found on the same scales to similarly-selected quasars. The small field of view of ALMA could miss a true overdensity if it exists on scales larger than 1\,cMpc, if the quasar is not centered in the overdensity, or if quasar feedback plays a role close to the quasar, but it is most likely that the large line of sight volume probed by a continuum survey will wash out a true overdensity signal. We discuss the necessary factors in determining the bias with which dusty star-forming galaxies trace true dark matter overdensities in order to improve upon overdensity searches in dust continuum.

Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission). These quasars have a median redshfit of z=3.17, absolute $i$ magnitude in the range $-29.67\leq M_i(z=2)\leq-27.03$, and different levels of radio-loudness. This sample unveils diverse specimens of Ly$\alpha$ nebulosities extending for tens of kiloparsecs around these quasars (on average out to a maximum projected distance of 80 kpc) above a surface brightness SB$>8.8\times10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ ($2\sigma$). The bulk of the extended Ly$\alpha$ emission is within R<50 kpc, and is characterized by relatively quiescent kinematics, with average velocity dispersions of $\langle \sigma_{\rm Ly\alpha}\rangle < 400$~km~s$^{-1}$. Therefore, the motions within all these Ly$\alpha$ nebulosities have amplitudes consistent with gravitational motions expected in dark matter halos hosting quasars at these redshifts, possibly reflecting the complexity in propagating a fast wind on large scales. Our current data suggest a combination of photoionization and resonant scattering as powering mechanisms of the Ly$\alpha$ emission. We discover the first $z\sim3$ ELAN, which confirms a very low probability ($\sim1\%$) of occurrence of such extreme systems at these cosmic epochs. Finally, we discuss the redshift evolution currently seen in extended Ly$\alpha$ emission around radio-quiet quasars from $z\sim3$ to $z\sim2$, concluding that it is possibly linked to a decrease of cool gas mass within the quasars' CGM from $z\sim3$ to $z\sim2$, and thus to the balance of cool vs hot media.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

We present the discovery of the radio-loud quasar PSO J352.4034-15.3373 at z=5.84 pm 0.02. This quasar is the radio brightest source known, by an order of magnitude, at z~6 with a flux density in the range of 8-100 mJy from 3GHz to 230MHz and a radio loudness parameter R>~1000. This source provides an unprecedented opportunity to study powerful jets and radio-mode feedback at the highest redshifts, and presents the first real chance to probe deep into the neutral intergalactic medium by detecting 21 cm absorption at the end of cosmic reionization.

We present a 45 ks Chandra observation of the quasar ULAS J1342+0928 at z=7.54. We detect 14.0^{+4.8}_{-3.7} counts from the quasar in the observed-frame energy range 0.5-7.0 keV (6-sigma detection), representing the most distant non-transient astronomical source identified in X-rays to date. The present data are sufficient only to infer rough constraints on the spectral parameters. We find an X-ray hardness ratio of HR = -0.51^{+0.26}_{-0.28} between the 0.5-2.0 keV and 2.0-7.0 keV ranges and derive a power-law photon index of Gamma = 1.95^{+0.55}_{-0.53}. Assuming a typical value for high-redshift quasars of Gamma = 1.9, ULAS J1342+0928 has a 2-10 keV rest-frame X-ray luminosity of L_{2-10} = 11.6^{+4.3}_{-3.5} x 10^{44} erg/s. Its X-ray-to-optical power-law slope is alpha_{OX}=-1.67^{+0.16}_{-0.10}, consistent with the general trend indicating that the X-ray emission in the most bolometrically powerful quasars is weaker relative to their optical emission.

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

BL Lac objects emitting in the Very High Energy (VHE) regime are unique tools to peer into the properties of the Extragalactic Background Light (EBL). However, due to the typical absence of features in their spectra, the determination of their redshifts has proven challenging. In this work we exploit the superb spatial resolution delivered by the new Advanced Rayleigh guided Ground layer adaptive Optics System (ARGOS) at the Large Binocular Telescope to detect the host galaxy of HESS J1943+213, a VHE emitting BL Lac shining through the Galaxy. Deep H-band imaging collected during the ARGOS commissioning allowed us to separate the contribution of the nuclear emission and to unveil the properties of the host galaxy with unprecedented detail. The host galaxy is well fitted by a S\'ersic profile with index of n~2 and total magnitude of H~16.15 mag. Under the assumption that BL Lac host galaxies are standard candles, we infer a redshift of z~0.21. In the framework of the current model for the EBL, this value is in agreement with the observed dimming of the VHE spectrum due to the scatter of energetic photons on the EBL.

Measuring the proximity effect and the damping wing of intergalactic neutral hydrogen in quasar spectra during the epoch of reionization requires an estimate of the intrinsic continuum at rest-frame wavelengths $\lambda_{\rm rest}\sim1200$-$1260$ {\AA}. In contrast to previous works which used composite spectra with matched spectral properties or explored correlations between parameters of broad emission lines, we opted for a non-parametric predictive approach based on principal component analysis (PCA) to predict the intrinsic spectrum from the spectral properties at redder (i.e. unabsorbed) wavelengths. We decomposed a sample of $12764$ spectra of $z\sim2$-$2.5$ quasars from SDSS/BOSS into 10 red-side ($1280$ {\AA} $<\lambda_{\rm rest}<2900$ {\AA}) and 6 blue-side ($1180$ {\AA} $<\lambda_{\rm rest}<1280$ {\AA}) PCA basis spectra, and constructed a projection matrix to predict the blue-side coefficients from a fit to the red-side spectrum. We found that our method predicts the blue-side continuum with $\sim6$-$12\%$ precision and $\lesssim1\%$ bias by testing on the full training set sample. We then computed predictions for the blue-side continua of the two quasars currently known at $z>7$: ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$). Both of these quasars are known to exhibit extreme emission line properties, so we individually calibrated the precision of the continuum predictions from similar quasars in the training set. We find that both $z>7$ quasars, and in particular ULAS J1342+0928, show signs of damping wing-like absorption at wavelengths redward of Ly$\alpha$.

We present a survey of the [CII] 158 $\mu$m line and underlying far-infrared (FIR) dust continuum emission in a sample of 27 z>6 quasars using the Atacama Large Millimeter Array (ALMA) at ~1" resolution. The [CII] line was significantly detected (at >5-sigma) in 23 sources (85%). We find typical line luminosities of $L_{\rm [CII]}=10^{9-10}$ L$_\odot$, and an average line width of ~385 km/s. The [CII]-to-far-infrared luminosity ratio ([CII]/FIR) in our sources span one order of magnitude, highlighting a variety of conditions in the star-forming medium. Four quasar host galaxies are clearly resolved in their [CII] emission on a few kpc scales. Basic estimates of the dynamical masses of the host galaxies give masses between $2\times10^{10}$ and $2\times10^{11}$ M$_\odot$, i.e., more than an order of magnitude below what is expected from local scaling relations, given the available limits on the masses of the central black holes ($>3\times10^8$ M$_\odot$, assuming Eddington-limited accretion). In stacked ALMA [CII] spectra of individual sources in our sample, we find no evidence of a deviation from a single Gaussian profile. The quasar luminosity does not strongly correlate with either the [CII] luminosity or equivalent width. This survey (with typical on-source integration times of 8 min) showcases the unparalleled sensitivity of ALMA at millimeter wavelengths, and offers a unique reference sample for the study of the first massive galaxies in the universe.

The properties of circum-galactic gas in the halo of quasar host galaxies are investigated analyzing Mg II 2800 and C IV 1540 absorption-line systems along the line of sight close to quasars. We used optical spectroscopy of closely aligned pairs of quasars (projected distance $\leq$ 200 kpc, but at very different redshift) obtained at the VLT and Gran Telescopio Canarias to investigate the distribution of the absorbing gas for a sample of quasars at z$\sim$1. Absorption systems of EW $\geq$ 0.3 $\rm{\AA}$ associated with the foreground quasars are revealed up to 200 kpc from the centre of the host galaxy, showing that the structure of the absorbing gas is patchy with a covering fraction quickly decreasing beyond 100 kpc. In this contribution we use optical and near-IR images obtained at VLT to investigate the relations between the properties of the circum-galactic medium of the host galaxies and of the large scale galaxy environments of the foreground quasars.

We present IRAM/NOEMA and JVLA observations of the quasar J1342+0928 at z=7.54 and report detections of copious amounts of dust and [CII] emission in the interstellar medium (ISM) of its host galaxy. At this redshift, the age of the universe is 690 Myr, about 10% younger than the redshift of the previous quasar record holder. Yet, the ISM of this new quasar host galaxy is significantly enriched by metals, as evidenced by the detection of the [CII] 158micron cooling line and the underlying far-infrared (FIR) dust continuum emission. To the first order, the FIR properties of this quasar host are similar to those found at a slightly lower redshift (z~6), making this source by far the FIR-brightest galaxy known at z>7.5. The [CII] emission is spatially unresolved, with an upper limit on the diameter of 7 kpc. Together with the measured FWHM of the [CII] line, this yields a dynamical mass of the host of <1.5x10^11 M_sun. Using standard assumptions about the dust temperature and emissivity, the NOEMA measurements give a dust mass of (0.6-4.3)x10^8 M_sun. The brightness of the [CII] luminosity, together with the high dust mass, imply active ongoing star formation in the quasar host. Using [CII]-SFR scaling relations, we derive star formation rates of 85-545 M_sun/yr in the host, consistent with the values derived from the dust continuum. Indeed, an episode of such past high star formation is needed to explain the presence of ~10^8 M_sun of dust implied by the observations.

We present final statistics from a survey for intervening MgII absorption towards 100 quasars with emission redshifts between $z=3.55$ and $z=7.08$. Using infrared spectra from Magellan/FIRE, we detect 279 cosmological MgII absorbers, and confirm that the incidence rate of $W_r>0.3 \AA$ MgII absorption per comoving path length does not evolve measurably between $z=0.25$ and $z=7$. This is consistent with our detection of seven new MgII systems at $z>6$, a redshift range that was not covered in prior searches. Restricting to relatively strong MgII systems ($W_r>1$\AA), there is significant evidence for redshift evolution. These systems roughly double in number density between $z=0$ and $z=2$-$3$, but decline by an order of magnitude from this peak by $z\sim 6$. This evolution mirrors that of the global star formation rate density, which could reflect a connection between star formation feedback and strong MgII absorbers. We compared our results to the Illustris cosmological simulation at $z=2$-$4$ by assigning absorption to catalogued dark-matter halos and by direct extraction of spectra from the simulation volume. To reproduce our results using the halo catalogs, we require circumgalactic (CGM) MgII envelopes within halos of progressively smaller mass at earlier times. This occurs naturally if we define the lower integration cutoff using SFR rather than mass. MgII profiles calculated directly from the Illustris volume yield far too few strong absorbers. We argue that this arises from unresolved phase space structure of CGM gas, particularly from turbulent velocities on sub-mesh scales. The presence of CGM MgII at $z>6$-- just $\sim 250$ Myr after the reionization redshift implied by Planck--suggests that enrichment of intra-halo gas may have begun before the presumed host galaxies' stellar populations were mature and dynamically relaxed. [abridged]

Quasars are galaxies hosting accreting supermassive black holes; due to their brightness, they are unique probes of the early universe. To date, only few quasars have been reported at $z > 6.5$ ($<$800 Myr after the Big Bang). In this work, we present six additional $z \gtrsim 6.5$ quasars discovered using the Pan-STARRS1 survey. We use a sample of 15 $z \gtrsim 6.5$ quasars to perform a homogeneous and comprehensive analysis of this highest-redshift quasar population. We report four main results: (1) the majority of $z\gtrsim$6.5 quasars show large blueshifts of the broad CIV 1549\AA$\,$emission line compared to the systemic redshift of the quasars, with a median value $\sim$3$\times$ higher than a quasar sample at $z\sim$1; (2) we estimate the quasars' black hole masses (M$\rm_{BH}\sim$0.3$-$5 $\times$ 10$^{9}$ M$_{\odot}$) via modeling of the MgII 2798\AA$\,$emission line and rest-frame UV continuum; we find that quasars at high redshift accrete their material (with $\langle (L_{\mathrm{bol}}/L_{\mathrm{Edd}}) \rangle = 0.39$) at a rate comparable to a luminosity-matched sample at lower$-$redshift, albeit with significant scatter ($0.4$ dex); (3) we recover no evolution of the FeII/MgII abundance ratio with cosmic time; (4) we derive near zone sizes; together with measurements for $z\sim6$ quasars from recent work, we confirm a shallow evolution of the decreasing quasar near zone sizes with redshift. Finally, we present new millimeter observations of the [CII] 158 $\mu$m emission line and underlying dust continuum from NOEMA for four quasars, and provide new accurate redshifts and [CII]/infrared luminosities estimates. The analysis presented here shows the large range of properties of the most distant quasars.

We utilize the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to search for extended Lyman-Alpha emission around the z~6.6 QSO J0305-3150. After carefully subtracting the point-spread-function, we reach a nominal 5-sigma surface brightness limit of SB = 1.9x10$^{-18}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, collapsing 5 wavelength slices centered at the expected location of the redshifted Lyman-Alpha emission (i.e. at 9256 Ang.). Current data suggest the presence (5-sigma, accounting for systematics) of a Lyman-Alpha nebula that extends for 9 kpc around the QSO. This emission is displaced and redshifted by 155 km/s with respect to the location of the QSO host galaxy traced by the [CII] emission line. The total luminosity is L = 3.0x10$^{42}$ erg/s. Our analysis suggests that this emission is unlikely to rise from optically thick clouds illuminated by the ionizing radiation of the QSO. It is more plausible that the Lyman-Alpha emission is due to fluorescence of the highly ionized optically thin gas. This scenario implies a high hydrogen volume density of n$_H$ ~ 6 cm$^{-3}$. In addition, we detect a Lyman-Alpha emitter (LAE) in the immediate vicinity of the QSO: i.e., with a projected separation of 12.5 kpc and a line-of-sight velocity difference of 560 km/s. The luminosity of the LAE is L = 2.1x10$^{42}$ erg/s and its inferred star-formation-rate is SFR ~ 1.3 M$_\odot$/yr. The probability of finding such a close LAE is one order of magnitude above the expectations based on the QSO-galaxy cross-correlation function. This discovery is in agreement with a scenario where dissipative interactions favour the rapid build-up of super-massive black holes at early Cosmic times.

The existence of massive ($10^{11}$ solar masses) elliptical galaxies by redshift z~4 (when the Universe was 1.5 billion years old) necessitates the presence of galaxies with star-formation rates exceeding 100 solar masses per year at z>6 (corresponding to an age of the Universe of less than 1 billion years). Surveys have discovered hundreds of galaxies at these early cosmic epochs, but their star-formation rates are more than an order of magnitude lower. The only known galaxies with very high star-formation rates at z>6 are, with only one exception, the host galaxies of quasars, but these galaxies also host accreting supermassive (more than $10^9$ solar masses) black holes, which probably affect the properties of the galaxies. Here we report observations of an emission line of singly ionized carbon ([CII] at a wavelength of 158 micrometres) in four galaxies at z>6 that are companions of quasars, with velocity offsets of less than 600 kilometers per second and linear offsets of less than 600 kiloparsecs. The discovery of these four galaxies was serendipitous; they are close to their companion quasars and appear bright in the far-infrared. On the basis of the [CII] measurements, we estimate star-formation rates in the companions of more than 100 solar masses per year. These sources are similar to the host galaxies of the quasars in [CII] brightness, linewidth and implied dynamical masses, but do not show evidence for accreting supermassive black holes. Similar systems have previously been found at lower redshift. We find such close companions in four out of twenty-five z>6 quasars surveyed, a fraction that needs to be accounted for in simulations. If they are representative of the bright end of the [CII] luminosity function, then they can account for the population of massive elliptical galaxies at z~4 in terms of cosmic space density.

Bright quasars, observed when the Universe was less than one billion years old (z>5.5), are known to host massive black holes (~10$^{9}$ M$_{\odot}$), and are thought to reside in the center of massive dark matter overdensities. In this picture, overdensities of galaxies are expected around high redshift quasars. However, observations based on the detection of Lyman Break Galaxies (LBGs) around these quasars do not offer a clear picture: this may be due to the uncertain redshift constraints of LBGs, which are selected through broad-band filters only. To circumvent such uncertainties, we here perform a search for Lyman Alpha Emitting galaxies (LAEs) in the field of the quasar PSO J215.1512-16.0417 at z~5.73, through narrow band, deep imaging with FORS2 at the VLT. We study an area of 37 arcmin$^{2}$, i.e. ~206 comoving Mpc$^{2}$ at the redshift of the quasar. We find no evidence for an overdensity of LAEs in the quasar field with respect to blank field studies. Possible explanations for these findings include that our survey volume is too small, or that the strong ionizing radiation from the quasar hinders galaxy formation in its immediate proximity. Another possibility is that these quasars are not situated in the dense environments predicted by some simulations.

We present ALMA observations of the inner 1' (1.2 kpc) of the Circinus galaxy, the nearest Seyfert. We target CO (1-0) in the region associated with a well-known multiphase outflow driven by the central active galactic nucleus (AGN). While the geometry of Circinus and its outflow make disentangling the latter difficult, we see indications of outflowing molecular gas at velocities consistent with the ionized outflow. We constrain the mass of the outflowing molecular gas to be 1.5e5 -5.1e6 solar masses, yielding a molecular outflow rate of 0.35-12.3 solar masses per year. The values within this range are comparable to the star formation rate in Circinus, indicating that the outflow indeed regulates star formation to some degree. The molecular outflow in Circinus is considerably lower in mass and energetics than previously-studied AGN-driven outflows, especially given its high ratio of AGN luminosity to bolometric luminosity. The molecular outflow in Circinus is, however, consistent with some trends put forth in Cicone et al. (2014), including a linear relation between kinetic power and AGN luminosity, as well as its momentum rate vs. bolometric luminosity (although the latter places Circinus among the starburst galaxies in that sample). We detect additional molecular species including CN and C17O.

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

We investigate the properties of the circumgalactic gas in the halo of quasar host galaxies from CIV absorption line systems. Optical spectroscopy of closely aligned pairs of quasars (projected distance \leq 200 kpc) obtained at the Gran Telescopio Canarias is used to investigate the distribution of the absorbing gas for a sample of 18 quasars at z \sim 2. We found that the detected absorption systems of EW \geq 0.3Ang associated with the foreground QSO are revealed up to 200 kpc from the center of the host galaxy. The structure of the absorbing gas is rather patchy with a covering fraction of the gas that quickly decreases beyond 100 kpc. These results are in qualitative agreement with those found for the lower ionisation metal Mg II 2800 Ang.

PKS 2155-304 is one of the brightest extragalactic source in the X-ray and EUV bands, and is a prototype for the BL Lac class of objects. In this paper we investigate the large-scale environment of this source using new multi-object as well as long-slit spectroscopy, together with archival spectra and optical images. We find clear evidence of a modest overdensity of galaxies at z=0.11610, consistent with previous determinations of the BL Lac redshift. The galaxy group has a radial velocity dispersion of 250km/s and a virial radius of 0.22Mpc, yielding a role-of-thumb estimate of the virial mass of M(vir)~1.5x10$^{13}$Msun, i.e., one order of magnitude less than what observed in other similar objects. This result hints toward a relatively wide diversity in the environmental properties of BL Lac objects.

The [CII] 158$\mu$m fine-structure line is known to trace regions of active star formation and is the main coolant of the cold, neutral atomic medium. In this \textit{Letter}, we report a strong detection of the [CII] line in the host galaxy of the brightest quasar known at $z>6.5$, the Pan-STARRS1 selected quasar PSO J036.5078+03.0498 (hereafter P036+03), using the IRAM NOEMA millimeter interferometer. Its [CII] and total far-infrared luminosities are $(5.8 \pm 0.7) \times 10^9 \,L_\odot$ and $(7.6\pm1.5) \times 10^{12}\,L_\odot$, respectively. This results in a $L_{[CII]} /L_{TIR}$ ratio of $\sim 0.8\times 10^{-3}$, which is at the high end for those found for active galaxies, though it is lower than the average found in typical main sequence galaxies at $z\sim 0$. We also report a tentative additional line which we identify as a blended emission from the $3_{22} - 3_{13}$ and $5_{23} - 4_{32}$ H$_2$O transitions. If confirmed, this would be the most distant detection of water emission to date. P036+03 rivals the current prototypical luminous J1148+5251 quasar at $z=6.42$, in both rest-frame UV and [CII] luminosities. Given its brightness and because it is visible from both hemispheres (unlike J1148+5251), P036+03 has the potential of becoming an important laboratory for the study of star formation and of the interstellar medium only $\sim 800\,$Myr after the Big Bang.

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{\nu,5{\rm GHz}}/f_{\nu,4400\AA}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

The number of known very high energy (VHE) blazars is $\sim\,50$, which is very small in comparison to the number of blazars detected in other frequencies. This situation is a handicap for population studies of blazars, which emit about half of their luminosity in the $\gamma$-ray domain. Moreover, VHE blazars, if distant, allow for the study of the environment that the high-energy $\gamma$-rays traverse in their path towards the Earth, like the extragalactic background light (EBL) and the intergalactic magnetic field (IGMF), and hence they have a special interest for the astrophysics community. We present the first VHE detection of 1ES\,0033+595 with a statistical significance of 5.5\,$\sigma$. The VHE emission of this object is constant throughout the MAGIC observations (2009 August and October), and can be parameterized with a power law with an integral flux above 150 GeV of $(7.1\pm1.3)\times 10^{-12} {\mathrm{ph\,cm^{-2}\,s^{-1}}}$ and a photon index of ($3.8\pm0.7$). We model its spectral energy distribution (SED) as the result of inverse Compton scattering of synchrotron photons. For the study of the SED we used simultaneous optical R-band data from the KVA telescope, archival X-ray data by \textit{Swift} as well as \textit{INTEGRAL}, and simultaneous high energy (HE, $300$\,MeV~--~$10$\,GeV) $\gamma$-ray data from the \textit{Fermi} LAT observatory. Using the empirical approach of Prandini et al. (2010) and the \textit{Fermi}-LAT and MAGIC spectra for this object, we estimate the redshift of this source to be $0.34\pm0.08\pm0.05$. This is a relevant result because this source is possibly one of the ten most distant VHE blazars known to date, and with further (simultaneous) observations could play an important role in blazar population studies, as well as future constraints on the EBL and IGMF.

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it was identified as a potential very high energy (VHE; E > 100 GeV) gamma-ray emitter. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 hours of good quality stereoscopic data. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3 {\sigma} (E > 70 GeV) during a 1.3-hour long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. We determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 +/- 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 +/- 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.

HESS J1857+026 is an extended TeV gamma-ray source that was discovered by H.E.S.S. as part of its Galactic plane survey. Given its broadband spectral energy distribution and its spatial coincidence with the young energetic pulsar PSR J1856+0245, the source has been put forward as a pulsar wind nebula (PWN) candidate. MAGIC has performed follow-up observations aimed at mapping the source down to energies approaching 100 GeV in order to better understand its complex morphology. HESS J1857+026 was observed by MAGIC in 2010, yielding 29 hours of good quality stereoscopic data that allowed us to map the source region in two separate ranges of energy. We detected very-high-energy gamma-ray emission from HESS J1857+026 with a significance of $12 \sigma$ above $150$ GeV. The differential energy spectrum between $100$ GeV and $13$ TeV is well described by a power law function $dN/dE = N_0(E/1\textrm{TeV})^{-\Gamma}$ with $N_0 = (5.37 \pm0.44_{stat} \pm1.5_{sys}) \times 10^{-12} (\textrm{TeV}^{-1} \textrm{cm}^{-2}$ $\textrm{ s}^{-1})$ and $\Gamma = 2.16\pm0.07_{stat} \pm0.15_{sys}$, which bridges the gap between the GeV emission measured by Fermi-LAT and the multi-TeV emission measured by H.E.S.S.. In addition, we present a detailed analysis of the energy-dependent morphology of this region. We couple these results with archival multi-wavelength data and outline evidence in favor of a two-source scenario, whereby one source is associated with a PWN, while the other could be linked with a molecular cloud complex containing an HII region and a possible gas cavity.

We investigate the properties of the galaxy environment of a sample of 14 low redshift (z $<$ 0.85) quasar physical pairs extracted from SDSS DR10 archives. The pairs have a systemic radial velocity difference $\Delta V_\parallel \leqslant$ 600 $km \ s^{-1}$ (based on [OIII]5007 \AA \ line) and projected distance $ R_\bot \leqslant$ 600 kpc. The physical association of the pairs is statistically confirmed at a level of $\sim$ 90 %. For most of the images of these quasars we are able to resolve their host galaxies that turn out to be on average similar to those of quasars not in pairs. We also found that quasars in a pair are on average in region of modest galaxy overdensity extending up 0.5 Mpc from the QSO. This galaxy overdensity is indistinguishable from that of a homogeneous sample of isolated quasars at the same redshift and with similar host galaxy luminosity. These results, albeit derived from a small (but homogeneous) sample of objects, suggest that the rare activation of two quasars with small physical separation does not require any extraordinary environment.

We study the multifrequency emission and spectral properties of the quasar 3C 279. We observed 3C 279 in very high energy (VHE, E>100GeV) gamma rays, with the MAGIC telescopes during 2011, for the first time in stereoscopic mode. We combine these measurements with observations at other energy bands: in high energy (HE, E>100MeV) gamma rays from Fermi-LAT, in X-rays from RXTE, in the optical from the KVA telescope and in the radio at 43GHz, 37GHz and 15GHz from the VLBA, Mets\"ahovi and OVRO radio telescopes and optical polarisation measurements from the KVA and Liverpool telescopes. During the MAGIC observations (February to April 2011) 3C 279 was in a low state in optical, X-ray and gamma rays. The MAGIC observations did not yield a significant detection. These upper limits are in agreement with the extrapolation of the HE gamma-ray spectrum, corrected for extragalactic background light absorption, from Fermi-LAT. The second part of the MAGIC observations in 2011 was triggered by a high activity state in the optical and gamma-ray bands. During the optical outburst the optical electric vector position angle rotatated of about 180 degrees. There was no simultaneous rotation of the 43GHz radio polarisation angle. No VHE gamma rays were detected by MAGIC, and the derived upper limits suggest the presence of a spectral break or curvature between the Fermi-LAT and MAGIC bands. The combined upper limits are the strongest derived to date for the source at VHE and below the level of the previously detected flux by a factor 2. Radiation models that include synchrotron and inverse Compton emissions match the optical to gamma-ray data, assuming an emission component inside the broad line region (BLR) responsible for the high-energy emission and one outside the BLR and the infrared torus causing optical and low-energy emission. We interpreted the optical polarisation with a bent trajectory model.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

High-redshift quasars are currently the only probes of the growth of supermassive black holes and potential tracers of structure evolution at early cosmic time. Here we present our candidate selection criteria from the Panoramic Survey Telescope & Rapid Response System 1 and follow-up strategy to discover quasars in the redshift range 5.7<z<6.2. With this strategy we discovered eight new 5.7<z<6.0 quasars, increasing the number of known quasars at z>5.7 by more than 10%. We additionally recovered 18 previously known quasars. The eight quasars presented here span a large range of luminosities (-27.3 < M_{1450} < -25.4; 19.6 < z_ps1 < 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show a weak or no Ly$\alpha$ emission line (25% with rest-frame equivalent width of the Ly$\alpha$ + Nv line lower than 15{\AA}). We find a larger fraction of weak-line emission quasars than in lower redshift studies. This may imply that the weak-line quasar population at the highest redshifts could be more abundant than previously thought. However, larger samples of quasars are needed to increase the statistical significance of this finding.

The Crab pulsar is the only astronomical pulsed source detected at very high energy (VHE, E>100GeV) gamma-rays. The emission mechanism of VHE pulsation is not yet fully understood, although several theoretical models have been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 4.6 years of {\it Fermi}-LAT data were also analyzed. The known two pulses per period were detected with a significance of $8.0 \sigma$ and $12.6 \sigma$. In addition, significant emission was found between the two pulses with $6.2 \sigma$. We discovered the bridge emission above 50 GeV between the two main pulses. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models.

We investigate the extent and the properties of the MgII cool, low-density absorbing gas located in the halo and in the circum-galactic environment of quasars, using a sample of 31 projected quasar pairs with impact parameter pd<200kpc in the redshift range 0.5<z<1.6. In the transverse direction, we detect 18 MgII absorbers associated with the foreground quasars, while no absorption system originated by the gas surrounding the quasar itself is found along the line-of-sight. This suggests that the quasar emission induces an anisotropy in the absorbing gas distribution. Our observations indicate that the covering fraction (fc) of MgII absorption systems with rest frame equivalent width Wr(2796)>0.3Ang. ranges from fc~1.0 at pd<65kpc to fc~0.2 at pd>150kpc, and appears to be higher than for galaxies. Our findings support a scenario where the luminosity/mass of the host galaxies affect the extent and the richness of the absorbing MgII circum-galactic medium.

The radio galaxy NGC 1275, recently identified as a very high energy (VHE, >100 GeV) gamma-ray emitter by MAGIC, is one of the few non-blazar AGN detected in the VHE regime. In order to better understand the origin of the gamma-ray emission and locate it within the galaxy, we studied contemporaneous multi-frequency observations of NGC 1275 and modeled the overall spectral energy distribution (SED). We analyzed unpublished MAGIC observations carried out between Oct. 2009 and Feb. 2010, and the previously published ones taken between Aug. 2010 and Feb. 2011. We studied the multi-band variability and correlations by analyzing data of Fermi-LAT (0.1-100 GeV), as well as Chandra (X-ray), KVA (optical) and MOJAVE (radio) data taken during the same period. Using customized Monte Carlo simulations corresponding to early MAGIC stereo data, we detect NGC 1275 also in the earlier campaign. The flux level and energy spectra are similar to the results of the second campaign. The monthly light curve >100 GeV shows a hint of variability at the 3.6 sigma level. In the Fermi-LAT band, both flux and spectral shape variabilities are reported. The optical light curve is variable and shows a clear correlation with the gamma-ray flux >100 MeV. In radio, 3 compact components are resolved in the innermost part of the jet. One of them shows a similar trend as the LAT and KVA light curves. The 0.1-650 GeV spectra measured simultaneously with MAGIC and Fermi-LAT can be well fitted either by a log-parabola or by a power-law with a sub-exponential cutoff for both campaigns. A single-zone synchrotron-self-Compton model, with an electron spectrum following a power-law with an exponential cutoff, can explain the broadband SED and the multi-band behavior of the source. However, this model suggests an untypical low bulk-Lorentz factor or a velocity alignment closer to the line of sight than the pc-scale radio jet.

We present the results of stereoscopic observations of the satellite galaxy Segue 1 with the MAGIC Telescopes, carried out between 2011 and 2013. With almost 160 hours of good-quality data, this is the deepest observational campaign on any dwarf galaxy performed so far in the very high energy range of the electromagnetic spectrum. We search this large data sample for signals of dark matter particles in the mass range between 100 GeV and 20 TeV. For this we use the full likelihood analysis method, which provides optimal sensitivity to characteristic gamma-ray spectral features, like those expected from dark matter annihilation or decay. In particular, we focus our search on gamma-rays produced from different final state Standard Model particles, annihilation with internal bremsstrahlung, monochromatic lines and box-shaped signals. Our results represent the most stringent constraints to the annihilation cross-section or decay lifetime obtained from observations of satellite galaxies, for masses above few hundred GeV. In particular, our strongest limit (95% confidence level) corresponds to a ~500 GeV dark matter particle annihilating into tau+tau-, and is of order <sigma v> ~ 1.2x10^{-24} cm^3 s^{-1} - a factor ~40 above the <sigma v> thermal value.

The bright gamma-ray quasar 4C +55.17 is a distant source ($z = 0.896$) with a hard spectrum at GeV energies as observed by the Large Area Telescope (LAT) on board the {{\it Fermi}} satellite. This source is identified as a good source candidate for very-high-energy (VHE; $> 30$ GeV) gamma rays. In general VHE gamma rays from distant sources provide an unique opportunity to study the extragalactic background light (EBL) and underlying astrophysics. The flux intensity of this source in the VHE range is investigated. Then, constraints on the EBL are derived from the attenuation of gamma-ray photons coming from the distant blazar. We searched for a gamma-ray signal from this object using the 35-hour observations taken by the MAGIC telescopes between November 2010 and January 2011. No significant VHE gamma-ray signal was detected. We computed the upper limits of the integrated gamma-ray flux at $95\%$ confidence level of $9.4 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ and $2.5 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ above $100$ GeV and $200$ GeV, respectively. The differential upper limits in four energy bins in the range from $80$ GeV to $500$ GeV are also derived. The upper limits are consistent with the attenuation predicted by low-flux EBL models on the assumption of a simple power-law spectrum extrapolated from LAT data.

Context. The radio galaxy IC 310 has recently been identified as a gamma-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E>100GeV) with the MAGIC telescopes. Originally classified as a head-tail radio galaxy, the nature of this object is subject of controversy since its nucleus shows blazar-like behavior. Aims. In order to understand the nature of IC 310 and the origin of the VHE emission we studied the spectral and flux variability of IC 310 from the X-ray band to the VHE gamma-ray regime. Methods. The light curve of IC 310 above 300GeV has been measured with the MAGIC telescopes from Oct. 2009 to Feb. 2010. Fermi-LAT data (2008-2011) in the 10-500GeV energy range were also analyzed. In X-ray, archival observations from 2003 to 2007 with XMM, Chandra, and Swift-XRT in the 0.5-10keV band were studied. Results. The VHE light curve reveals several high-amplitude and short-duration flares. Day-to-day flux variability is clearly present. The photon index between 120GeV and 8TeV remains at the value $\Gamma\sim2.0$ during both low and high flux states. The VHE spectral shape does not show significant variability, whereas the flux at 1TeV changes by a factor of $\sim7$. Fermi-LAT detected only eight gamma-ray events in the energy range 10GeV-500GeV in three years of observation. The measured photon index of $\Gamma=1.3\pm0.5$ in the Fermi-LAT range is very hard. The X-ray measurements show strong variability in flux and photon index. The latter varied from $1.76\pm0.07$ to $2.55\pm0.07$. Conclusion. The rapid variability measured confirms the blazar-like behavior of IC 310. The TeV emission seems to originate from scales of less than 80 Schwarzschild radii within the compact core of its FRI radio jet with orientation angle 10deg-38deg. The SED resembles that of an extreme blazar, albeit the luminosity is more than two orders of magnitude lower.

Indications of a GeV component in the emission from GRBs are known since the EGRET observations during the 1990's and they have been confirmed by the data of the Fermi satellite. These results have, however, shown that our understanding of GRB physics is still unsatisfactory. The new generation of Cherenkov observatories and in particular the MAGIC telescope, allow for the first time the possibility to extend the measurement of GRBs from several tens up to hundreds of GeV energy range. Both leptonic and hadronic processes have been suggested to explain the possible GeV/TeV counterpart of GRBs. Observations with ground-based telescopes of very high energy photons (E>30 GeV) from these sources are going to play a key role in discriminating among the different proposed emission mechanisms, which are barely distinguishable at lower energies. MAGIC telescope observations of the GRB 090102 (z=1.547) field and Fermi Large Area Telescope (LAT) data in the same time interval are analysed to derive upper limits of the GeV/TeV emission. We compare these results to the expected emissions evaluated for different processes in the framework of a relativistic blast wave model for the afterglow. Simultaneous upper limits with Fermi and a Cherenkov telescope have been derived for this GRB observation. The results we obtained are compatible with the expected emission although the difficulties in predicting the HE and VHE emission for the afterglow of this event makes it difficult to draw firmer conclusions. Nonetheless, MAGIC sensitivity in the energy range of overlap with space-based instruments (above about 40 GeV) is about one order of magnitude better with respect to Fermi. This makes evident the constraining power of ground-based observations and shows that the MAGIC telescope has reached the required performance to make possible GRB multiwavelength studies in the very high energy range.

We present Very Large Telescope optical spectroscopy of nine BL Lac objects of unknown redshift belonging to the list of optically selected radio loud BL Lacs candidates. We explore their spectroscopic properties and the possible link with gamma ray emission. From the new observations we determine the redshift of four objects from faint emission lines or from absorption features of the host galaxy. In three cases we find narrow intervening absorptions from which a lower limit to the redshift is inferred. For the remaining two featureless sources, lower limits to the redshift are deduced from the very absence of spectral lines. A search for gamma counterpart emission shows that six out of nine are Fermi gamma-ray emitters with two new detections. Our analysis suggests that most of the BL Lac still lacking of redshift information are most probably located at high redshift.

Motivated by the Costamante & Ghisellini (2002) predictions we investigated if the blazar 1ES 1727+502 (z=0.055) is emitting very high energy (VHE, E>100 GeV) gamma rays. We observed the BL Lac object 1ES 1727+502 in stereoscopic mode with the two MAGIC telescopes during 14 nights between May 6th and June 10th 2011, for a total effective observing time of 12.6 hours. For the study of the multiwavelength spectral energy distribution (SED) we use simultaneous optical R-band data from the KVA telescope, archival UV/optical and X-ray observations by instruments UVOT and XRT on board of the Swift satellite and high energy (HE, 0.1 GeV - 100 GeV) gamma-ray data from the Fermi-LAT instrument. We detect, for the first time, VHE gamma-ray emission from 1ES 1727+502 at a statistical significance of 5.5 sigma. The integral flux above 150 GeV is estimated to be (2.1\pm0.4)% of the Crab Nebula flux and the de-absorbed VHE spectrum has a photon index of (2.7\pm0.5). No significant short-term variability was found in any of the wavebands presented here. We model the SED using a one-zone synchrotron self-Compton model obtaining parameters typical for this class of sources.

We present optical spectroscopy of the nearest quasar pair listed in the 13th edition of the Veron-Cetty & Veron catalogue, i.e. the two quasars SDSS J15244+3032 and RXS J15244+3032 (redshift z~0.27, angular separation ~7 arcsec, and line-of-sight velocity difference ~1900 km/s). This system would be an optimal candidate to investigate the mutual interaction of the host galaxies with ground based optical imaging and spectroscopy. However, new optical data demonstrate that RXS J15244+3032 is indeed a star of spectral type G. This paper includes data gathered with the Asiago 1.82m telescope (Cima Ekar Observatory, Asiago, Italy).

We present the discovery of a triplet of quasars at z~1.51. The whole system is well accommodated within 25 arcsec (i.e., 200 kpc in projected distance). The velocity differences among the three objects (as measured through the broad MgII emission line) are less than 1000 km/s, suggesting that the quasars belong to the same physical structure. Broad band NIR images of the field do not reveal evidence of galaxies or galaxy clusters that could act as a gravitational lens, ruling out the possibility that two or all the three quasars are multiple images of a single, strongly lensed source. QQQ J1519+0627 is the second triplet of quasars known up to date. We estimate that these systems are extremely rare in terms of simple accidental superposition. The lack of strong galaxy overdensity suggests that this peculiar system is harboured in the seeds of a yet-to-be-formed massive structure. Based on observations collected at the La Silla Observatory with the New Technology Telescope (NTT) of the European Southern Observatory (ESO) and at the Calar Alto Observatory with the 3.5m telescope of the Centro Astron\'omico Hispano Alem\'an (CAHA).

We present optical spectroscopy of projected QSO pairs to investigate the MgII and the CIV absorption features imprinted on the spectrum of the background object by the gaseous halo surrounding the foreground QSO. We observed 13 projected pairs in the redshift range 0.7<z<2.2 spanning projected separations between 60 kpc and 120 kpc. In the spectra of the background QSOs, we identify MgII intervening absorption systems associated to the foreground QSOs in 7 out of 10 pairs, and 1 absorption system out of 3 is found for CIV. The distribution of the equivalent width as a function of the impact parameter shows that, unlike the case of normal galaxies, some strong absorption systems (EWr > 1 Ang) are present also beyond a projected radius of ~70 kpc. If we take into account the mass of the galaxies as an additional parameter that influence the extent of the gaseous haloes, the distribution of the absorptions connected to the QSOs is consistent to that of galaxies. In the spectra of the foreground QSOs we do not detect any MgII absorption lines originated by the gas surrounding the QSO itself, but in 2 cases these features are present for CIV. The comparison between the absorption features observed in the transverse direction and those along the line of sight allows us to comment on the distribution of the absorbing gas and on the emission properties of the QSOs. Based on observations undertaken at the European Southern Observatory (ESO) Very Large Telescope (VLT) under Programmes 085.B-0210(A) and 086.B-0028(A).

PKS0447-439 is one of the BL Lacertae objects that have been detected at very high energy. There has been a recent report of a lower limit of z>1.246 for the redshift of this blazar, challenging the current paradigm in which very high-energy gamma-rays cannot freely propagate in the z>1 universe. In this research note, we present a new MagE/Magellan spectrum of PKS0447-439 with exquisite signal-to-noise (S/N>150 at 6500A). Our analysis confirms the presence of the previously-reported absorption line at 6280A, which we identify, however, with a known telluric absorption, invalidating the claim that this blazar lies at z>1. Since no other extragalactic spectral features are detected, we cannot establish a redshift based on our spectrum.

We investigate the optical spectral properties of the blazar PKS 1222+216 during a period of 3 years. While the continuum is highly variable the broad line emission is practically constant. This supports a scenario in which the broad line region is not affected by jet continuum variations. We thus infer the thermal component of the continuum from the line luminosity and we show that it is comparable with the continuum level observed during the phases of minimum optical activity. The mass of the black hole is estimated through the virial method from the FWHM of MgII, Hbeta, and Halpha broad lines and from the thermal continuum luminosity. This yields a consistent black hole mass value of 6x10^8 solar masses.

The dynamical properties of six SDSS quasar pairs at z < 0.8 are investigated. The pairs have proper transverse separation < 500 kpc, and velocity difference along the line of sight < 500 km/s. If they are bound systems their dynamical mass can be evaluated and compared with that of host galaxies. Evidence is found of an excess of the former mass with respect to the latter. This suggests that these quasar pairs are hosted by galaxies with massive dark halos or that they reside in a group/cluster of galaxies.