\section{Introduction}

Galactic structure in the Milky Way is difficult to determine because
of our location within the disk and the difficulties in deriving
accurate distances.  There is strong evidence, however, that the Milky
Way is a barred spiral galaxy \citep[e.g.,][]{churchwell2009}.  Most
barred spiral galaxies are Grand Design galaxies, with two prominent
symmetric spiral arms \citep{elmegreen1982}.  Understanding the Milky
Way structure is important since bars and spiral arms have dynamical
effects (e.g., radial migration) and influence star formation (e.g.,
shock gas).

To our knowledge, \citet{dame2011} were the first to connect the
Scutum-Centaurus (SC) spiral arm from its beginning at the end of the
bar in the inner Galaxy to the outer Galaxy in the first quadrant.
They called this outermost section the Outer Scutum-Centaurus (OSC)
arm. In the first quadrant, the OSC is at a distance between
\(15\ensuremath{\,{\rm kpc}}\xspace\) and \(19\ensuremath{\,{\rm kpc}}\xspace\) kpc from the Galactic Center (GC) and
between \(20\ensuremath{\,{\rm kpc}}\xspace\) and \(25\ensuremath{\,{\rm kpc}}\xspace\) from the Sun \citep{armentrout2017}.
The connection between the SC and OSC arm segments, however, remains
hypothetical. Toward the GC, velocity crowding and complex structure
prohibits clear identification of spiral structure.  The H\,{\sc i}\ emission
that has now been associated with the OSC arm was detected long
ago. The arm appears in the early 21 cm maps of both \citet{kerr1969}
and \citet{weaver1974}.  Yet it took decades to trace the SC spiral
arm from the bar to the outer Galaxy.

\citet{dame2011} carefully traced the OSC in the Leiden/Argentine/Bonn
(LAB) H\,{\sc i}\ 21\ensuremath{\,{\rm cm}}\xspace\ line all-sky survey data \citep{hartmann1997,
  arnal2000, bajaja2005, kalberla2005} and then, using the Center for
Astrophysics 1.2 m telescope, detected molecular gas in the arm for
the first time, at 10 locations coincident with HI emission peaks. A
CO map was made of one location revealing a molecular cloud with mass
and radius of \nexpo{5}{4}\msun\ and 47\ensuremath{\,{\rm pc}}\xspace, respectively.
\citet{koo2017} recently re-analyzed the LAB survey data using a
peak-finding algorithm and identified the OSC arm as a 20\ensuremath{\,{\rm kpc}}\xspace\ long
H\,{\sc i}\ structure coincident with several H\,{\sc ii}\ regions and molecular
clouds.

\citet{sun2015} suggested a possible extension of the OSC from the
first Galactic quadrant into the second quadrant.  Their results are
part of the Milky Way Imaging Scroll Painting (MWISP) project to map
\co\ and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ from Galactic longitudes $-10^{\ensuremath{\,^\circ}\xspace} < \ell <
250^{\ensuremath{\,^\circ}\xspace}$ and Galactic latitudes $-5^{\ensuremath{\,^\circ}\xspace} < b < 5^{\ensuremath{\,^\circ}\xspace}$
with the Purple Mountain Observatory Delingha 13.7\ensuremath{\,{\rm m}}\xspace\ telescope.  They
detected 72 molecular clouds with masses between
$10^{2}-10^{4}$\msun\ that lie between $120^{\ensuremath{\,^\circ}\xspace} < \ell <
150^{\ensuremath{\,^\circ}\xspace}$ and are roughly consistent with an extension of the OSC
assuming a log-periodic spiral arm model.  First quadrant data from
the MWISP for the longitude range of $35^{\ensuremath{\,^\circ}\xspace} < \ell <
45^{\ensuremath{\,^\circ}\xspace}$ reveal 168 molecular clouds consistent with the OSC
($\ell, v$) locus, defined by $V_{\rm LSR} = -1.6\, \ell \pm
13.2$\ensuremath{\,{\rm km\,s^{-1}}}\xspace. These molecular clouds have typical masses and
sizes of \nexpo{3}{3}\msun\ and 5\ensuremath{\,{\rm pc}}\xspace, respectively \citep{sun2017}.

Detection of CO does not necessarily imply star
formation. Observations of high-mass star formation tracers are
therefore critical for characterizing the OSC. H\,{\sc ii}\ regions are the
archetypical tracers of star formation and spiral structure. As part
of the H\,{\sc ii}\ Region Discovery Survey, \citet{anderson2015} discovered
six H\,{\sc ii}\ regions whose radio recombination line (RRL) emission
velocities are consistent with the OSC ($\ell, v$) locus defined by
\(V_{\rm LSR} = -1.6\, \ell \pm 15\ensuremath{\,{\rm km\,s^{-1}}}\xspace\).  There are four additional
H\,{\sc ii}\ regions in the {\it WISE} Catalog of Galactic H\,{\sc ii}\ Regions
consistent with this definition of the OSC \citep{anderson2012}.

\begin{figure*}[!htb]
\centering
\includegraphics[width=0.49\linewidth]{G026.418+1.683_12co.eps}
\includegraphics[width=0.49\linewidth]{G021.541+1.676_12co.eps}
\includegraphics[width=0.49\linewidth]{G054.491+1.579_12co.eps}
\includegraphics[width=0.49\linewidth]{G042.210+1.081_12co.eps}
\caption{Representative \co\ spectra for different quality factors
  (QF).  Plotted is the antenna temperature as a function of the LSR
  velocity.  The black curves are the data and the red curves are
  Gaussian fits to the data.  The LSR velocity of each Gaussian
  profile, together with the QF, is shown in the right-hand corner of the
  plot.}
\label{fig:qf}
\end{figure*}

\begin{figure*}[!htb]
\centering
\includegraphics[width=0.45\linewidth]{G050.901+1.056.eps}
\includegraphics[width=0.45\linewidth]{G061.154+2.170.eps} \\
\includegraphics[width=0.45\linewidth]{G062.194+1.715.eps}
\includegraphics[width=0.45\linewidth]{G066.608+2.061.eps}
\caption{\co\ (black) and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ (blue) spectra of targets where the
  association of \co\ emission and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission was not
  straightforward.  In G050.901+1.056 and G062.194+1.715 there are two
  \ensuremath{^{\rm 13}{\rm CO}}\xspace\ components that correspond to one \co\ component. In
  G066.608+2.061 there is one \ensuremath{^{\rm 13}{\rm CO}}\xspace\ component that corresponds to two
  \co\ components. Finally, in G061.154+2.170 there are
  \co\ components at velocities where the \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission is corrupted
  by emission in the Off position.}
\label{fig:special}
\end{figure*}

RRL emission is the best tracer of high-mass star formation, but it is
faint at the large distances of the OSC. \citet{armentrout2017} took a
different approach and used the Green Bank Telescope (GBT) to observe
the dense molecular gas tracers NH$_{3}$ (J,K) = (1,1), (2,2), (3,3)
and H$_{2}$O 6(1,6) $\rightarrow$ 5(2,3) toward 75 {\it WISE}
H\,{\sc ii}\ regions and H\,{\sc ii}\ region candidates located within the OSC
($\ell, b$) locus as defined by H\,{\sc i}\ \citep{dame2011}. Because these
targets were identified as H\,{\sc ii}\ region candidates based on their
mid-infrared morphology, these molecular lines trace dense gas that is
more likely to be associated with star formation. Any detected
spectral line {\it probably} indicates an active star formation
region.  About 20\% of the targets were detected in either ammonia or
water maser emission, but only two have velocities consistent with the
OSC.

\citet{armentrout2017} also observed a similar sample of OSC targets
in radio continuum at 8--10\ensuremath{\,{\rm GHz}}\xspace\ with the Jansky Very Large Array
(JVLA). About 60\% of the targets were detected in radio continuum
with the JVLA. Five of these are associated with the locus of the OSC.
Together, RRL and continuum emission allowed various H\,{\sc ii}\ region
physical properties to be derived. Associating the radio continuum
with the molecular transition, however, is less secure since the
molecular cloud may not be associated with the H\,{\sc ii}\ region
\citep{anderson2009}. Nonetheless, assuming they are related yields a
distance. If the hydrogen-ionizing photon flux is produced by a single
star, then the observed radio continuum fluxes imply spectral types
ranging from O4 to O8.5. This suggests that even at such large
Galactocentric distances ($> 15$\ensuremath{\,{\rm kpc}}\xspace) high-mass star formation is
ongoing in the Milky Way.

Here we seek to explore the molecular content in the OSC using
\co\ and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ transitions, which are brighter than the molecular line
transitions observed by \citet{armentrout2017}. Since molecular clouds
need not be forming high-mass stars, carbon monoxide is not as good a
tracer of high-mass star formation as RRLs, NH$_{3}$, or H$_{2}$O, but
detection is much more likely at these large distances.  Following
\citet{armentrout2017} we only target H\,{\sc ii}\ regions and H\,{\sc ii}\ region
candidates from the {\it WISE} Catalog of Galactic
H\,{\sc ii}\ Regions. These targets all have the same characteristic infrared
morphology which increases the probability that the molecular gas is
associated with high-mass star formation.

\begin{figure*}[!htb]
\centering
\includegraphics[width=0.32\linewidth]{G028.319+1.243_12co.eps}
\includegraphics[width=0.32\linewidth]{G033.008+1.151_12co.eps}
\includegraphics[width=0.32\linewidth]{G038.627+0.813_12co.eps}
\includegraphics[width=0.32\linewidth]{G039.185-1.421_12co.eps}
\includegraphics[width=0.32\linewidth]{G040.955+2.473_12co.eps}
\includegraphics[width=0.32\linewidth]{G041.304+1.997_12co.eps}
\includegraphics[width=0.32\linewidth]{G041.810+1.503_12co.eps}
\includegraphics[width=0.32\linewidth]{G042.210+1.081_12co.eps}
\includegraphics[width=0.32\linewidth]{G046.368+0.802_12co.eps}
\includegraphics[width=0.32\linewidth]{G046.375+0.897_12co.eps}
\includegraphics[width=0.32\linewidth]{G054.094+1.749_12co.eps}
\includegraphics[width=0.32\linewidth]{G061.085+2.502_12co.eps}
\includegraphics[width=0.32\linewidth]{G061.180+2.447_12co.eps}
\includegraphics[width=0.32\linewidth]{G062.578+2.387_12co.eps}
\caption{\co\ spectra of emission lines originating within the OSC.
  Plotted is the antenna temperature as a function of the LSR
  velocity.  The black curves are the data and the red curves are
  Gaussian fits to the data.  The LSR velocity of each Gaussian
  profile, together with the QF, is shown in the right-hand corner of
  the plot. The vertical black line indicates the RRL velocity, when
  available.}
\label{fig:co12_osc}
\end{figure*}

\begin{figure*}[!htb]
\centering
\includegraphics[width=0.32\linewidth]{G033.008+1.151_13co.eps}
\includegraphics[width=0.32\linewidth]{G039.185-1.421_13co.eps}
\includegraphics[width=0.32\linewidth]{G040.955+2.473_13co.eps}
\includegraphics[width=0.32\linewidth]{G041.304+1.997_13co.eps}
\includegraphics[width=0.32\linewidth]{G046.368+0.802_13co.eps}
\includegraphics[width=0.32\linewidth]{G046.375+0.897_13co.eps}
\includegraphics[width=0.32\linewidth]{G061.180+2.447_13co.eps}
\caption{\ensuremath{^{\rm 13}{\rm CO}}\xspace\ spectra of emission lines originating within the OSC.
  Plotted is the antenna temperature as a function of the LSR
  velocity.  The black curves are the data and the red curves are
  Gaussian fits to the data.  The LSR velocity of each Gaussian
  profile, together with the QF, is shown in the right-hand corner of
  the plot. The vertical black line indicates the RRL velocity, when
  available.}
\label{fig:co13_osc}
\end{figure*}

\section{Observations and Data Reduction}

We used the Arizona Radio Observatory (ARO) 12\ensuremath{\,{\rm m}}\xspace\ telescope to observe
the \co\ J=1--0 transition at 115.27120\ensuremath{\,{\rm GHz}}\xspace\ toward 78 H\,{\sc ii}\ regions
and H\,{\sc ii}\ region candidates located in the first Galactic quadrant. We
selected all H\,{\sc ii}\ region candidates in the {\it WISE} Catalog of
Galactic H\,{\sc ii}\ Regions \citep{anderson2012} that lay near the OSC
($\ell, b$) locus defined as \(b = 0.375^\circ + 0.075\times\ell\)
within the range \(20^\circ < \ell < 70^\circ\). Many of these targets
have measured radio continuum emission and are most likely bona fide
H\,{\sc ii}\ regions \citep{armentrout2017}. In addition, we included 10
H\,{\sc ii}\ regions from \citep{armentrout2017} with RRL velocities within
the OSC (\(\ell, V\)) locus defined as $V_{\rm LSR} = -1.6\, \ell \pm
15\ensuremath{\,{\rm km\,s^{-1}}}\xspace$.  Since \co\ is often optically thick in Galactic molecular
clouds, we observed the optically thin \ensuremath{^{\rm 13}{\rm CO}}\xspace\ J=1--0 transition at
110.20132\ensuremath{\,{\rm GHz}}\xspace\ in a subset of targets with bright \co\ detections to
provide a more accurate measure of the molecular column density.

The ARO 12\ensuremath{\,{\rm m}}\xspace\ telescope is the European ALMA prototype antenna that
began operation on Kitt Peak in 2014.  The telescope's half-power
beam-width (HPBW) is 54\arcsec\ and 57\arcsec\ at 115\ensuremath{\,{\rm GHz}}\xspace\ and
110\ensuremath{\,{\rm GHz}}\xspace\ , respectively. The main beam efficiency is
\(\gsim90\%\) at these frequencies.  Our observations were performed
between February 11--16, 2016. The receiver consisted of the ALMA band
3, dual polarization, sideband-separating mixers with typical on-sky
system temperatures of $\sim 300$\ensuremath{\,{\rm K}}\xspace. We employed both the filter bank
spectrometer with 256 channels at 2\ensuremath{\,{\rm MHz}}\xspace\ spectral resolution
(5\ensuremath{\,{\rm km\,s^{-1}}}\xspace\ at 115\ensuremath{\,{\rm GHz}}\xspace), and the millimeter autocorrelator (MAC) with 4096
channels at 195\ensuremath{\,{\rm kHz}}\xspace\ spectral resolution (0.5\ensuremath{\,{\rm km\,s^{-1}}}\xspace\ at 115\ensuremath{\,{\rm GHz}}\xspace).  Both
spectrometers accept two intermediate frequency signals consisting of
the two orthogonal polarizations.  Here we only consider data from the
MAC spectrometer given the narrow line widths of CO ($\sim 1$\ensuremath{\,{\rm km\,s^{-1}}}\xspace).
The typical atmospheric optical depth at zenith was \(\tau_0 \sim
0.2\).

We made total power, position switched observations where the
reference position (Off) is offset by 20\arcmin\ in azimuth from the
source (On).  The On and Off positions were observed for 5 minutes
each with a switching rate of 30 seconds, for a total time of 10
minutes.  Typically 1--2 total power pairs were sufficient to detect
\co emission but in some cases we integrated longer.  The telescope
pointing and focus were corrected every 1--2 hours by peaking on
Jupiter, Venus, or Saturn.  The typical pointing accuracy was $\sim
2$\arcsec\ rms.  At the start of each session we checked the tuning of
the spectrometers by observing the test source M17SW.

The data were reduced and analyzed using TMBIDL, an IDL single-dish
software package \citep{bania2016}.\footnote{V7.1, see
  \url{https://github.com/tvwenger/tmbidl}} The data reduction and
analysis were performed independently by three of the authors.  Each
spectrum was visually inspected. We discarded \({\sim}1\%\) of the
spectra due to poor baseline structure.  For each target the data were
averaged over all total power pairs and polarizations to produce a
single, averaged spectrum.  We modeled the spectral baselines with a
third-order polynomial function that was subtracted from the data to
remove any sky continuum emission or residual baseline structure in
the spectrum.  We fitted a Gaussian function to each profile using a
Levenberg-Markwardt least squares method \citep{markwardt2009} to
derive the peak intensity, the full width at half-maximum (FWHM) line
width, and the local standard of rest (LSR) velocity.

We compared the results of each independent analysis. The results were
similar for \({\sim}90\%\) of the spectra.  Any differences typically
involved how best to fit complex, blended profiles.  Multiple
\co\ emission components were detected in about 2/3 of the targets.
These directions thus contain several \co\ clouds along the
line-of-sight.  Furthermore, because CO is pervasive in the inner
Galactic plane, there was often \co\ emission detected in the Off
position. This creates apparent absorption lines in the processed
position-switched spectrum. Since any emission components near these
absorption features are contaminated by this Off emission, we
eliminated these components from any further analysis. The
\ensuremath{^{\rm 13}{\rm CO}}\xspace\ spectra were analyzed in the same manner.

\begin{figure*}[!htb]
\centering
\includegraphics[width=\linewidth]{ARO_LB_LV.eps}
\caption{The outer Scutum-Centaurus spiral arm as traced by integrated
  H\,{\sc i}\ emission, CO emission, and H\,{\sc ii}\ regions. Top:
  velocity-integrated H\,{\sc i}\ emission (grayscale image) tracing the OSC
  spiral arm, summed over a 14\ensuremath{\,{\rm km\,s^{-1}}}\xspace\ wide window following the center
  velocity given by $V_{\rm LSR} = -1.6\ensuremath{\,{\rm km\,s^{-1}}}\xspace\,{\rm deg}^{-1} \times
  \ell$.  CO emission detected here toward {\it WISE} H\,{\sc ii}\ regions is
  shown as red circles with sizes scaled by the \co\ peak line
  intensity.  Gray stars correspond to H\,{\sc ii}\ regions consistent with
  the location of the OSC \citep{armentrout2017}.  Bottom:
  longitude-velocity diagram of H\,{\sc i}\ emission (grayscale image),
  summed over a 3\ensuremath{\rlap.{^{\circ}}} 5 window following the arm in latitude with
  $b = $0\ensuremath{\rlap.{^{\circ}}} 75 + 0.75 $\times \ell$.  The symbols are the same as
  in the top panel.  The OSC locus is defined by $V_{\rm LSR} = -1.6\,
  \ell \pm 15$\ensuremath{\,{\rm km\,s^{-1}}}\xspace\ and shown by the dashed lines. (Reproduced from
  \citet{dame2011} and \citet{armentrout2017}.)}
\label{fig:lv_diagram}
\end{figure*}

We assigned a Quality Factor (A, B, or C) to every CO emission
line. This judgment was based on the signal-to-noise ratio and
expected line properties (e.g., the line widths should be a few
\(\ensuremath{\,{\rm km\,s^{-1}}}\xspace\) at most). A quality factor of ``A'' was given to profiles
with a signal-to-noise ratio larger than 10 with no blending and flat
spectral baselines.  Targets with either lower signal-to-noise ratio
or some blending of lines were given a quality factor ``B''.  Here the
blending should not be so severe that the individual peaks could not
be visually detected.  A quality factor of ``C'' was given to the
remaining targets which were often blended and more difficult to fit.
A quality factor of D was assigned to spectra with no CO detections.

\section{Results and Discussion}

We detect CO toward $\sim 80$\% of our targets. About 2/3 of our
targets have at least one CO emission line at negative velocity,
placing it beyond the Solar orbit.  The results are summarized in
Figures~\ref{fig:qf}--\ref{fig:co13_osc} and
Tables~\ref{tab:co12}--\ref{tab:co13}.  We measure 117 distinct
\co\ components toward 62 of 78 targets, and 40 \ensuremath{^{\rm 13}{\rm CO}}\xspace\ components
toward 27 of 30 targets. Representative \co\ spectra that span the
range of quality factors A--C are shown in Figure~\ref{fig:qf}.  All
line intensities are in units of the antenna temperature corrected for
atmospheric attenuation, radiative loss, and rearward scattering and
spillover, $T_{\rm A}^{*}$.  The Gaussian fit parameters are shown in
Table~\ref{tab:co12} and Table~\ref{tab:co13} for the \co\ and
\ensuremath{^{\rm 13}{\rm CO}}\xspace\ spectra, respectively.  Listed are the target name, the Gaussian
fit parameters, the root-mean-square noise in the line-free baseline,
RMS, the total integration time, $t_{\rm intg}$, and the quality
factor, QF.  The Gaussian fits include the peak line intensity,
$T_{\rm L}$, the LSR velocity, $V_{\rm LSR}$, and the FWHM line width,
$\Delta{V}$, and their associated 1-$\sigma$ uncertainties.  For
\co\ spectra with multiple emission components, in
Table~\ref{tab:co12} we append a letter at the end of the target name
(a, b, c, etc.) in order of decreasing line intensity.

In most cases we can associate each \ensuremath{^{\rm 13}{\rm CO}}\xspace\ component with a
corresponding \co\ component using the LSR velocity (i.e. the \co\ and
\ensuremath{^{\rm 13}{\rm CO}}\xspace\ components have the same LSR velocity, within \(2\ensuremath{\,{\rm km\,s^{-1}}}\xspace\)). There
are four exceptions.  In G050.901+1.056 and G062.194+1.715 there are
two \ensuremath{^{\rm 13}{\rm CO}}\xspace\ components that correspond to one \co\ component (identified
as G050.901+1.056a, G050.901+1.056b, G062.194+1.715a1, and
G062.194+1.715a2 in Table~\ref{tab:co13}). The reverse occurs in
G066.608+2.061 where there is one \ensuremath{^{\rm 13}{\rm CO}}\xspace\ component corresponding to two
\co\ components (identified as G066.608+2.061ab in
Table~\ref{tab:co13}). These differences may be caused by opacity
effects where the \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission is arising from a region that is
optically thick in \co. This can cause self-absorption in the
\co\ spectrum. The signal-to-noise ratio will also affect the ability
to resolve multiple components. Finally, for G061.154+2.170 there are
\co\ components at velocities where the \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission is corrupted by
emission in the Off position. The \co\ and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ spectra for these
four targets are shown in Figure~\ref{fig:special}.

Following \citet{armentrout2017}, we deem molecular emission
components to be located in the OSC whenever the molecular line LSR
velocity is within the ($\ell, v$) locus defined by $V_{\rm LSR} =
-1.6\, \ell \pm 15$\ensuremath{\,{\rm km\,s^{-1}}}\xspace.  We find that 17 \co\ emission lines toward
14 targets and 8 \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission lines toward 7 targets originate
within the OSC. Of these detections, 7 \co\ targets (G028.319+1.243,
G033.008+1.151, G039.185-1.421, G041.304+1.997, G041.810+1.503,
G054.094+1.749, G062.578+2.387) and 3 \ensuremath{^{\rm 13}{\rm CO}}\xspace\ targets (G033.008+1.151,
G039.185-1.421, G041.304+1.997) have a measured RRL velocity in the
\textit{WISE} Catalog of Galactic H\,{\sc ii}\ regions \citep{anderson2012}.
We show spectra for all of our OSC detections in
Figure~\ref{fig:co12_osc} and \ref{fig:co13_osc} for \co\ and \ensuremath{^{\rm 13}{\rm CO}}\xspace,
respectively. Two OSC targets have multiple \co\ or \ensuremath{^{\rm 13}{\rm CO}}\xspace\ spectral
line components that lie within the OSC: G038.627+0.813 and
G046.375+0.897.  Our OSC targets may not all be unique objects.  There
are two target pairs --- G046.368+0.802/G046.375+0.897 and
G061.085+2.502/G061.180+2.447 --- that lie near the same ($\ell, b,
v$) locus and therefore may be part of the same star formation
complex.  This may account for the multiple component spectrum seen
toward G046.375+0.897.  Finally, six of our OSC targets are probably
related to sources detected by \citet{sun2017} in CO emission:
G038.627+0.813ab (MWISP G38.533+0.892), G039.185$-$1.421 (MWISP
G39.175$-$1.425), G040.955+2.473 (MWISP G40.958+2.483), G041.304+1.997
(MWISP G41.308+2.000), G041.810+1.503 (MWISP G41.758+1.567, MWISP
G41.733+1.517, MWISP G41.742+1.458), and G042.210+1.081 (MWISP
G42.192+1.083).

\startlongtable
\input{sources12.tex}

\subsection{CO Distribution throughout the Milky Way}

How are the molecular clouds with detected CO emission distributed in
the Galaxy?  Our targets were selected from the {\it WISE}
H\,{\sc ii}\ region catalog and therefore their distribution should reflect
the Galactic H\,{\sc ii}\ region or high-mass star formation distribution.
We expect that in many cases the CO emission is associated with the
H\,{\sc ii}\ region along the line-of-sight, but this is not always true
\citep{russeil2004,anderson2009}.

The most direct view is given by the observed, model independent,
parameters: $\ell, b, v$.  Figure~\ref{fig:lv_diagram} shows the
($\ell, b$) (top) and the ($\ell, v$) diagrams (bottom) for the OSC
region of the Galaxy.  Similar plots are found in \citet{dame2011} and
\citet{armentrout2017}.  Integrated H\,{\sc i}\ emission is shown together
with our \co\ detections (red circles).  The molecular clouds trace
the warp in the disk as defined by H\,{\sc i}\ emission but are concentrated
at the lower latitude end of the H\,{\sc i}\ envelope.  The bottom panel of
Figure~\ref{fig:lv_diagram} shows the molecular clouds located in the
OSC as red circles between the dashed lines. Several of these are
clearly associated with H\,{\sc ii}\ regions indicated by the gray stars.

Since the majority of our targets have multiple emission components,
we cannot identify which CO component might be associated with the
H\,{\sc ii}\ region.  Nevertheless, most of the molecular clouds in our
sample have negative velocities and are thus located beyond the Solar
orbit. In fact, Figure~\ref{fig:lv_diagram} shows that the majority of
the negative velocity components CO seem to be located in the Outer
Arm and not the more-distant OSC.

\startlongtable
\input{sources13.tex}

\startlongtable
\input{properties.tex}

Face-on maps of the Milky Way disk are difficult to make since
distances to the sources are required.  Parallax distances are the
most accurate but can only be measured for a relatively small sample
\citep[e.g., see][]{reid2014}.  Kinematic distances are far easier to
derive since they only require a spectral line velocity and a model of
Galactic rotation.  The main limitations are that the rotation curves
suffer from non-circular (streaming) motions, and that within the
Solar orbit there is a kinematic distance ambiguity
\citep[e.g.,][]{anderson2012}.  Studies of the OSC arm avoid the
latter problem because the OSC is located beyond the Solar orbit and
so does not suffer from the distance ambiguity issue.

\begin{figure}[!htb]
\centering
\includegraphics[width=\linewidth]{ARO_FaceOn.eps}
\caption{Face-on map of H\,{\sc ii}\ regions, CO emission components, and
  H\,{\sc i}\ emission in the first Galactic quadrant.  Plotted are
  H\,{\sc ii}\ regions from the {\it WISE} H\,{\sc ii}\ region catalog (gray
  crosses), H\,{\sc ii}\ regions consistent with the location of the OSC
  (gray stars), and CO emission detected here in the direction of {\it
    WISE} H\,{\sc ii}\ regions (red circles).  The red circle sizes are
  scaled by the intensity of the \co\ peak line intensity.  The Solar
  circle and tangent points are indicated by the solid and dashed
  curves, respectively.  The dotted lines connect the near and far
  kinematic distances for sources within the Solar orbit since we did
  not resolve the kinematic distance ambiguity.  The gray scale shows
  OSC H\,{\sc i}\ emission from the LAB survey. Here the OSC is defined in
  latitude and velocity as: $V_{\rm LSR} = -1.6\,{\rm km}\,{\rm
    s}^{-1}\,{\rm deg}^{-1} \times \ell$ and $b = $0\ensuremath{\rlap.{^{\circ}}} 75 + 0.75
  $\times \ell$.  This is the same H\,{\sc i}\ emission shown between the
  dashed lines in Figure~\ref{fig:lv_diagram} (bottom panel), but
  transformed using the Brand rotation curve
  \citep{brand1993}. (Reproduced from \citet{armentrout2017}.)}
\label{fig:faceon_diagram}
\end{figure}

We use the \citet{brand1993} rotation curve to determine kinematic
distances to H\,{\sc i}\ and CO emission as well as H\,{\sc ii}\ regions.  The
resulting face-on map for the first quadrant OSC zone studied here is
shown in Figure~\ref{fig:faceon_diagram} where the gray crosses are
taken from the {\it WISE} H\,{\sc ii}\ region catalog and the red circles are
CO emission from this study. We also include previously detected OSC
H\,{\sc ii}\ regions (gray stars) and H\,{\sc i}\ data (gray contours). The
H\,{\sc ii}\ regions and molecular clouds associated with the OSC span a wide
range of longitudes in the first Galactic quadrant.  Their
distribution is well-matched to the OSC extent delineated by
H\,{\sc i}\ emission.

\begin{figure}[!htb]
\centering
\includegraphics[width=\linewidth]{fwhm_vs_R.eps}
\caption{\co\ (black) and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ (blue) emission line FWHM line width,
  \(\Delta V\), as a function of Galactocentric radius, \(R\). The
  black circles and black squares are the \co\ and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ FWHM line
  widths, respectively. The error bars are the \(1\sigma\)
  uncertainties in the measured FWHM line widths.}
\label{fig:fwhm_vs_R}
\end{figure}

\begin{figure}[!htb]
\centering
\includegraphics[width=\linewidth]{NH2_vs_R.eps}
\caption{Molecular column density, \(N({\rm H}_2)\), as a function of
  Galactocentric radius, \(R\), for molecular clouds with
  \ensuremath{^{\rm 13}{\rm CO}}\xspace\ detections.}
\label{fig:NH2_vs_R}
\end{figure}

\subsection{Comparison of OSC and non-OSC CO Clouds}

We search for differences between the emission line properties of
molecular clouds in the OSC and clouds not in the OSC.
Figure~\ref{fig:fwhm_vs_R} shows the \co\ and \ensuremath{^{\rm 13}{\rm CO}}\xspace\ FWHM line widths
as a function of Galactocentric radius, \(R\). The ratio of the mean
\co\ FWHM line width of OSC clouds to non--OSC clouds is \(1.13\) with
a standard deviation of \(0.52\). For \ensuremath{^{\rm 13}{\rm CO}}\xspace\, the ratio is \(1.23\)
with a standard deviation of \(0.74\). This evidence suggests
molecular cloud velocity dispersions are fairly uniform for all
molecular clouds in our sample. We note, however, that our
observations are a single pointing in the direction of the cloud
rather than a full map. Our derived FWHM line widths are therefore
likely lower limits of the true molecular cloud FWHM line widths,
especially if the clouds are much larger than our beam
(\({\sim}1\)\arcmin\ ).

The \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission lines should be optically thin and may be used to
derive the total molecular column densities of the molecular clouds.
We derive the total H\(_2\) column density of each of our
\ensuremath{^{\rm 13}{\rm CO}}\xspace\ detections using Equation (2) from \citet{simon2001}. This
equation assumes a \co/\ensuremath{^{\rm 13}{\rm CO}}\xspace\ abundance ratio of 45, an excitation
temperature of \(10\,\ensuremath{\,{\rm K}}\xspace\), and an X-factor (\co/H$_{2}$) of
\nexpo{8}{-5}.  We assume the main beam efficiency is near unity and
use the \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission line properties from Table~\ref{tab:co13}.  The
molecular cloud column densities are compiled in Table~\ref{tab:prop}
where we list the \co\ component name, the LSR velocity, the
Galactocentric radius, the distance from the Sun, the \co\ and
\ensuremath{^{\rm 13}{\rm CO}}\xspace\ line intensities, and the molecular column density.  Sources
within the Solar circle have two values for the distance since we did
not resolve the kinematic distance ambiguity.

The total molecular column density of these clouds also does not have
any trend with \(R\). Figure~\ref{fig:NH2_vs_R} shows the derived
H\(_2\) column densities as a function of \(R\) for each of our clouds
with \ensuremath{^{\rm 13}{\rm CO}}\xspace\ detections. The ratio of the mean molecular column
densities of OSC clouds to non--OSC clouds is \(0.97\) with a standard
deviation of \(1.11\). This evidence again suggests no difference
between OSC and non--OSC molecular clouds within our sample with the
caveat that our pointed observations likely underestimate the true
molecular cloud FWHM line widths, and therefore these derived column
densities may represent lower limits.

Since we expect molecular clouds to be extended in CO emission and
have not mapped these regions, we cannot accurately derive the total
molecular mass of these complexes.  The eight sources from the MWISP
survey in common with our survey have molecular masses and radii that
range from ${\sim}200-7000$\msun\ and ${\sim}3-13$\ensuremath{\,{\rm pc}}\xspace, respectively.
Since the molecular masses were calculated using \co\ intensities that
are likely optically thick, these should be interpreted as lower
limits.

\section{Conclusion}

The Outer Scutum-Centaurus (OSC) spiral arm is the most distant site
of massive star formation in the Milky Way. The OSC contains at least
12 known H\,{\sc ii}\ regions ionized by high-mass stars with spectral types
as early as O4.  Molecular clouds have been detected within the OSC
but most have not been associated with H\,{\sc ii}\ regions.  Here we use the
Arizona Radio Observatory (ARO) 12\ensuremath{\,{\rm m}}\xspace\ telescope to observe CO emission
from 78 {\it WISE} H\,{\sc ii}\ region candidates located within the Galactic
longitude-latitude ($\ell, b$) locus of the OSC spiral arm in the
first quadrant ($20^{\circ} < \ell < 90^{\circ}$).  We detect 117
\co\ spectral line components in 62 of 78 targets, and 40
\ensuremath{^{\rm 13}{\rm CO}}\xspace\ components in 27 of 30 targets.  About 2/3 of the molecular
clouds reside beyond the Solar orbit and are associated with the Outer
Arm. We discovered 17 \co\ emission lines and 8 \ensuremath{^{\rm 13}{\rm CO}}\xspace\ emission lines
consistent with the ($\ell, v$) OSC locus. These OSC molecular clouds
have the same physical properties (FWHM line widths and molecular
column densities) as non--OSC clouds within our sample.

\acknowledgments

We thank T. M. Dame and the anonymous referee for their useful
comments and suggestions which improved the quality of this
paper. T.V.W. is supported by the NSF through the Grote Reber
Fellowship Program administered by Associated Universities,
Inc./National Radio Astronomy Observatory, the D.N. Batten Foundation
Fellowship from the Jefferson Scholars Foundation, the Mars Foundation
Fellowship from the Achievement Rewards for College Scientists
Foundation, and the Virginia Space Grant Consortium.  L.D.A. is
supported by NSF grant AST1516021. The Kitt Peak 12 Meter is operated
by the Arizona Radio Observatory (ARO), Steward Observatory,
University of Arizona. \nraoblurb

\facilities{ARO$\,12\,$m}

\software{TMBIDL \citep{bania2016}}


\section{Keywords:} active galactic nuclei, $\gamma$-ray emission, super-massive black hole, narrow-line Seyfert 1 galaxy, blazar, relativistic jet, accretion process, cosmological evolution}
\end{abstract}

\section{Introduction}

Relativistic jets are one of the most powerful manifestations of the release
of energy related to the super-massive black hole (SMBH) at the center of
active galactic nuclei (AGN). In about 10\% of AGN, termed radio-loud AGN, the accretion disc is at
the base of a bipolar outflow of relativistic plasma, which may extend well
beyond the host galaxy, forming the spectacular lobes of plasma visible in the
radio band. The jet emission is observed across the entire electromagnetic
spectrum, from radio to $\gamma$ rays. When the jet axis is closely aligned
with our line of sight, the rest-frame radiation is strongly amplified due to
the Doppler boosting with a large fraction of the output observed at higher energies, and giving rise to the blazar phenomenon. Blazars are traditionally
divided into flat spectrum radio quasars (FSRQ) and BL Lac objects, based on
the presence or not, respectively, of broad emission lines (i.e. Equivalent Width $>$
5 \AA) in their optical and UV spectrum \citep[e.g.,][]{stickel91}. Recently,
a new classification based on the luminosity of the broad line
region (BLR) in Eddington luminosity was proposed by \citet{ghisellini11}: sources with
L$_{BLR}$/L$_{Edd}$ higher or lower than 5$\times$10$^{-4}$ being classified
as FSRQ or BL Lacs, respectively, in agreement with a transition of the
accretion regime from efficient to inefficient between these classes. The
spectral energy distribution (SED) of blazars are characterized by two bumps with
the lower energy peak occurring in the IR/optical band
in the FSRQ and at UV/X-rays in the BL Lacs. This first peak is univocally interpreted as synchrotron radiation from
highly relativistic electrons in a jet. The SED higher energy peak, observed at MeV--GeV energies in the FSRQ and at GeV--TeV energies in the BL Lacs, is commonly interpreted as inverse Compton scattering of seed photons, internal or external to the jet,
by the same relativistic electrons \citep[e.g.,][]{ulrich97}. However, other models involving hadronic processes have been proposed \citep[e.g.,][]{boettcher13}. \citet{fossati98} proposed that the SEDs of blazars form a spectral sequence, with the position of the two peaks governed by the observed bolometric luminosity: blazars with lower luminosities have the peaks at higher energies. This was theoretically interpreted by
\citet{ghisellini98} in terms of different radiative cooling suffered by the electrons emitting at the two peaks.

Radio galaxies are viewed at larger angles than blazars, with less severe
boosting effects and classified as Fanaroff-Riley type I (FR I) and type II (FR II) based on their radio power and morphology \citep{fanaroff74}. Following
the Unified model of AGN proposed by \citet{urry95}, FR I and FR II radio
galaxies are the non-aligned (to the observer viewing angle) parent
populations of BL Lac and FSRQ, respectively. The discovery by {\em Fermi}-LAT
of variable $\gamma$-ray emission from a few radio-loud narrow-line Seyfert 1
galaxies (NLSy1) suggested this as the third class of AGN with a relativistic jet \citep{abdo09}.

The {\em Fermi  Gamma-ray  Space  Telescope} with its main instrument on-board, the Large Area Telescope (LAT), opened a new era in the study of
high-energy emission from AGN. The $\gamma$-ray sky observed by {\em Fermi}-LAT is dominated by AGN that are $\sim$60\% of the sources
included in the Third {\em Fermi} LAT source catalog \citep[3FGL;][]{acero15}. Apart from a handful of starburst galaxies, for which the $\gamma$-ray emission originates from the interaction of cosmic rays with gas and interstellar radiation fields \citep{ackermann12}, almost all the extragalactic sources are associated with radio-loud AGN.

Important findings regarding the blazar population included in the 3FGL and the NLSy1 are discussed in Section \ref{3LAC} and \ref{NLSy1}, respectively. The high-redshift blazars detected in $\gamma$ rays are discussed in Section \ref{highredshift}, while the third catalog of hard LAT sources is presented in Section \ref{3FHL}. Concluding remarks are reported in Section \ref{concluding}.

\section{The Third LAT AGN Catalogue}\label{3LAC}

Several $\gamma$-ray source catalogs, both general and for specific class of
objects (i.e. AGN, pulsars, supernova remnants, pulsar wind nebulae, gamma-ray
bursts) have already been produced by the LAT Collaboration. The 3FGL was used
as a starting point for producing a catalog of AGN only: the Third LAT AGN
catalog \citep[3LAC;][]{ackermann15}. The 3FGL includes 3033 sources detected
after four years of operation with a Test Statistic\footnote{The Test Statistic is defined as TS = 2$\times$(log$L_1$ - log$L_0$), where $L$ is the likelihood of the data given the model with ($L_1$) or without ($L_0$) a point source at the position of the target \citep[e.g.][]{mattox96}} greater than 25,
corresponding to a significance $>$ 4$\sigma$; 2192 sources are detected at Galactic latitude $\mid$b$\mid$ $>$ 10$^{\circ}$. 1563 sources (71\% of the 3FGL objects at $\mid$b$\mid$ $>$ 10$^{\circ}$) are associated at high-confidence with 1591 AGN (28 objects have two possible
associations), which constitute, together with the low-latitude AGN, the 3LAC sample. Most of the high-latitude objects (98\%) are blazars. The 3LAC includes radio-loud AGN of different types: 467 FSRQ, 632 BL Lacs, 460 BCU (Blazar with Uncertain Classification) and 32 non-blazar AGN. Among the non-blazar AGN there are 12 FR I, 3 FR II, 8 steep spectrum radio sources and 5 NLSy1.

Associations for 182 low-latitude ($\mid$b$\mid$ $<$ 10$^{\circ}$) AGN are reported in the 3LAC: 24 FSRQ, 30 BL Lacs, 125 BCU and 3 non-blazar AGN. Extrapolating the number of high-latitude sources to low-latitude and assuming the same sensitivity, $\sim$340 sources would have been expected. The discrepancy between expected and associated source numbers is likely due to a higher Galactic diffuse emission background and a higher incompleteness of the catalogs of counterparts at low-latitude.

On average, sources with high $\gamma$-ray luminosity (mostly FSRQ) are found
to have softer spectra than sources with low $\gamma$-ray luminosity  (mostly BL
Lacs), in agreement with the ``blazar sequence'' \citep{ghisellini98} and the
``blazar divide'' \citep{ghisellini09}, with the exception of a few outliers
(high-synchrotron-peaked-FSRQ and high-luminosity high-synchrotron-peaked BL
Lacs). A strong anti-correlation between the synchrotron peak position
($\nu_{peak}$) and the spectral index ($\Gamma_{\gamma}$) is observed for FSRQ
and BL Lacs. A similar trend is noticed for BCU supporting the idea that BCU
with low $\nu_{peak}$ and high $\Gamma_{\gamma}$ are likely FSRQ, while BCU
with high $\nu_{peak}$ and low $\Gamma_{\gamma}$ are likely BL Lacs.
 Based on the 3LAC results, the blazar sequence was recently revised by \citet{ghisellini17}: FSRQ display approximately the same SED as the luminosity increases,
 following a sequence only in Compton dominance (i.e. the ratio of the Compton to synchrotron peak luminosities) and in the X–ray slope, while
 BL Lacs become `redder' (i.e. the peak frequencies becomes smaller) when more luminous. Moreover, a correlation between the
 jet power and the accretion power has been found in blazars, with the
 jet power dominating over the accretion disc luminosity by a factor of 10 and somewhat
 larger than the entire gravitational power \citep{ghisellini14}.

At the time of writing, 70 AGN have been detected at TeV energies and listed in the TeVCat\footnote{http://tevcat.uchicago.edu/}. All these sources are present in the 3LAC except for HESS J1943$+$213. A large fraction ($>$ 75\%) of blazars detected by {\em Swift}-BAT in hard X-rays have also been detected by {\em Fermi}-LAT.

The blazars detected by {\em Fermi}-LAT in $\gamma$ rays after 4 years of operation represent a sizeable fraction of the entire population of known blazars as listed in the BZCAT \citep{massaro09}. The overall LAT-detected fraction is 24\% (409/1707)
for FSRQ, 44\% (543/1221) for BL Lacs and 27\% (59/221) for BCU. No strong
differences in the radio, optical and X-ray flux distributions are observed
between $\gamma$-ray detected and non-detected blazars in the BZCAT, suggesting that
all blazars could eventually shine in $\gamma$ rays at LAT-detection
level. More BL Lacs than FSRQ have been detected by {\em Fermi}-LAT so
far. This may be related to a longer duty cycle of FSRQ with respect to BL
Lacs in $\gamma$ rays.

\section{Narrow-line Seyfert 1 galaxies}\label{NLSy1}

The discovery of variable $\gamma$-ray emission from a few NLSy1 confirmed the
presence of relativistic jets in these objects. In addition to the 5 objects
reported in the 3LAC, {\em Fermi}-LAT has recently detected $\gamma$ rays from
other 3 new NLSy1: FBQS J1644$+$2619 \citep{dammando15}, B3 1441$+$476 and
NVSS J124634$+$476 \citep{dammando16}. Luminosity, variability, and spectral
properties of these NLSy1 in $\gamma$ rays indicate a blazar-like behaviour
\citep[e.g.,][]{dammando16}. Apparent superluminal jet components were detected in SBS 0846+513 \citep{dammando12}, PMN J0948$+$0022, and 1H
0323$+$342 \citep{lister16}, supporting the presence of relativistic jets in
this class of objects.

The detection of relativistic jets in a class of AGN thought to be hosted in
spiral galaxies with a BH mass typically of 10$^{6}$-10$^{7}$ M$_{\odot}$
\citep[e.g.,][]{deo06}, challenges the theoretical scenarios of jet formation \citep[e.g.,][]{boettcher02}, suggesting two possible interpretations: either relativistic jets in NLSy1 are produced by a different mechanism or the BH mass in NLSy1 is largely underestimated.

In the last years it has been claimed that the BH mass of NLSy1 maybe underestimated due either to the effect of radiation pressure from ionizing photons on BLR \citep{marconi08} or to projection effects \citep{baldi16}. By considering these effects, NLSy1 have BH masses of 10$^{8}$--10$^{9}$ M$_\odot$, in agreement with the values estimated by modelling the optical-UV part of their spectra with a Shakura and Sunyaev disc spectrum \citep[e.g.,][]{calderone13}. This may solve the problem of the minimum BH mass predicted in different theoretical scenarios of relativistic jet formation, but leaves open the host galaxy issue.

Spiral galaxies are usually formed by minor mergers, with BH masses typically ranging between 10$^{6}$--10$^{7}$ M$_\odot$ \citep[e.g.,][]{woo02}, so it would not be clear how powerful relativistic jets could form in these galaxies. It is worth mentioning that the morphological classification has been done mainly for radio-quiet nearby NLSy1. Among the NLSy1 detected by {\em Fermi}-LAT up to now, the morphology of the host galaxy has been investigated only for 1H 0323$+$342, PKS 2004$-$447, and FBQS J1644$+$2619. Observations of 1H 0323$+$342 with the Hubble Space Telescope and the Nordic Optical Telescope revealed a structure that may be interpreted either as a one-armed spiral galaxy \citep{zhou07} or as a circumnuclear ring produced by a recent merger \citep{anton08, leon14}. A pseudo-bulge morphology of the host galaxy of the NLSy1 PKS 2004$-$447 and FBQS J1644$+$2619 have been claimed by \citet{kotilainen16} and \citet{olguin17}, respectively, but no conclusive results have been obtained so far. Hence, it is crucial to determine the type of galaxy hosting $\gamma$-ray emitting NLSy1 and their BH mass.

For this reason near-infrared observations in $J$ band of FBQS J1644$+$2619
were performed using the Canarias Infrared Camera Experiment (CIRCE) at the
Gran Telescopio Canaries. The 2D surface brightness profile of the source is
modelled up to 5 arcsec by the combination of a nuclear component, associated
with the AGN contribution, and a bulge component with a S\'ersic index $n$ =
3.7, indicative of an elliptical galaxy. The structural parameters of the host
are consistent with the correlations of effective radius and surface
brightness against absolute magnitude measured for elliptical galaxies. From
the infrared bulge luminosity, a BH mass of (2.1 $\pm$ 0.2)$\times$10$^{8}$
M$_{\odot}$ was estimated \citep{dammando17}. All these pieces of evidence
strongly indicate that the relativistic jet in the NLSy1 FBQS J1644$+$2619 is produced by a massive BH in an elliptical galaxy, as expected for radio-loud AGN.

\begin{figure}[h!]
\begin{center}
\includegraphics[width=10cm]{Highz.jpeg}
\end{center}
\caption{Photon index vs. $\gamma$-ray luminosity plane for 3LAC and high-redshift blazars at $z$ $>$ 3.1 detected in $\gamma$ rays. Adapted from \citet{ackermann17}}\label{fig1}
\end{figure}

\section{High-redshift blazars}\label{highredshift}

High-redshift blazars are the most powerful radio-loud AGN in the Universe and
are bright targets in hard X-rays, representing a significant fraction of the
extragalactic hard X-ray sky. However, they are not commonly detected in
$\gamma$ rays. In fact, high-redshift ($z > $ 2) blazars represent $<$10\% of the AGN population observed by {\em Fermi}-LAT so far. Flaring activity in the $\gamma$-ray energy range from these sources is even more uncommon, with only fourteen FSRQ at $z$ $>$ 2 detected during a $\gamma$-ray flare.

In the 3LAC there are 64 objects at $z$ $>$ 2 ($\sim$3.7 per cent of the $\gamma$-ray sources associated with AGN) at $z$ $>$ 2, and only 2 at $z$ $>$ 3: PKS 0537$-$286 ($z$ = 3.104) and TXS 0800$+$618 ($z$ = 3.033). On the other hand, 13 blazars at $z$ $>$ 3 have been detected in hard X-rays by {\em Swift}-BAT, {\em INTEGRAL}-IBIS, and {\em NuSTAR} so far. Hard X-ray observations are more suitable for detecting blazars at $z$ $>$ 3, and this is mainly due to a spectral bias. In fact, the inverse Compton peak of high-redshift blazars is shifted towards lower energies as the bolometric luminosity increases. Only 10 sources at $z$ $>$ 2 are in both the 3LAC and the {\em Swift}-BAT 70-month catalog \citep{baumgartner13}. All
high-redshift blazars listed in both 3LAC and {\em Swift}-BAT catalogues have an average $\gamma$-ray luminosity L$_{\gamma}$ $>$ 2$\times$10$^{48}$ erg s$^{-1}$, indicating that only the most luminous blazars have been detected by both instruments. Furthermore, only blazars with an X-ray photon index $\Gamma_{\rm\,X}$ $<$1.6 have been detected in $\gamma$ rays, while no dependence on the X-ray spectral luminosity seems to be present \citep{dammando15b}.

As said above, high-redshift blazars at $z >$ 3.1 are missing in the {\em Fermi}
catalogs. These objects typically have large bolometric luminosities
(L$_{bol}$ $>$ 10$^{48}$ erg s$^{-1}$) and harbor extremely massive BH
(M$_{BH}$ $\sim$ 10$^{9}$ M$_{\odot}$). The new Pass 8 data set, with an improved event-level analysis, substantially enhances the
sensitivity of the LAT, in particular at lower energies, increasing the
capability of the LAT to detect sources with soft spectra like the
high-redshift blazars. By analysing 92 months of Pass 8 data between 60 MeV
and 300 GeV of a large sample of radio-loud quasars, 5 new $\gamma$-ray emitting blazars at $z$ $>$ 3.1 have
been detected with high significance. Among them, NVSS J151002+570243 ($z$ =
4.31) is now the most distant $\gamma$-ray emitting blazar so far. All the blazars discovered show steep $\gamma$-ray spectra
($\Gamma_{\gamma}$ $>$ 2.5), indicating an IC peak at MeV energies. These five
sources lie in the region of high $\gamma$-ray luminosities (L$_{\gamma}$ $>$
10$^{47}$ erg s$^{-1}$) and soft photon indices (Fig. \ref{fig1}) typical of
powerful blazars \citep{ackermann17}. Among the 5 new high-redshift blazars
there are (at least) two with redshift between 3 and 4 with a M$_{BH}$ $>$
10$^{9}$ M$_{\odot}$, implying the presence of 2$\times$2$\Gamma^{2}$
(i.e. 675, adopting $\Gamma$ = 13) similar objects but with a misaligned jet
in the same range of redshift. This changes the estimate of the space density
of very massive BH hosted in jetted sources to 68$^{+36}_{-24}$ Gpc$^{-3}$. As a consequence at $z \sim$ 4 we should have a similar number of SMBH hosted in radio-loud and radio-quiet sources and, given their strong evolution, above that redshift most massive BH might be hosted in radio-loud AGN. This suggests that the radio-loud phase can be a key ingredient for a rapid BH growth in the early Universe.

\section{The Third Catalog of Hard Fermi-LAT Sources}\label{3FHL}

\begin{figure}[h!]
\begin{center}
\includegraphics[width=10cm]{CH.jpeg}
\end{center}
\caption{The highest-energy photons vs. redshift for 3FHL blazars, color coded by the optical depth calculated from the model presented by \citet{dominguez11}. The cosmic $\gamma$-ray horizon based on different EBL models is shown. Adapted from \citet{ajello17}}\label{fig2}
\end{figure}

In addition to the {\em Fermi}-LAT catalogs with the standard low-energy threshold of 100
MeV, three hard source catalogs have been released: the First {\em Fermi}-LAT
Catalog of Sources above 10 GeV \citep[1FHL;][]{ackermann13}, based on the first three
years of data analyzed in the 10--500 GeV energy range, the Second Catalog of Hard
{\em Fermi}-LAT Sources \citep[2FHL;][]{ackermann16}, based on 80 months of
data analyzed in the 50 GeV--2 TeV energy range, and the Third Catalog of Hard
{\em Fermi}-LAT Sources \citep[3FHL;][]{ajello17}, based on 7 years of data in the
10 GeV--2 TeV energy range. The 3FHL contains 1556 objects and takes advantage
of the improvement provided by Pass 8 by using the PSF-type event
classification. The 3FHL includes 214 new $\gamma$-ray sources never appeared
in previous {\em Fermi} catalogs. Three of these 214 have been detected with
the Imaging Atmospheric Cherenkov Telescopes (IACT). The vast majority of
detected sources (79\%) are associated with extragalactic counterparts at
other wavelengths, including 16 sources located at high-redshift ($z$ $>$ 2):
11 FSRQ, 3 BL Lac, and 2 BCU. BL Lacs are the most numerous extragalactic
population (61\%) followed by BCU (23\%) and FSRQ (14\%). Only 72 of the 3FHL
extragalactic sources have been already detected by current IACT. In this context, the 3FHL is a resource for planning observations of the current
(MAGIC, VERITAS, H.E.S.S.) and future (Cherenkov Telescope Array) IACT observatories.
Interestingly, a few highest-energy photons from distant blazars included in the 3FHL catalogue are
in the region around and beyond the cosmic $\gamma$-ray horizon (i.e. the energy at which the cosmic optical depth $\tau$ = 1, see e.g. \citet{dominguez13}) as shown in Fig. \ref{fig2}. These photons provide important constraints on extragalactic background light models as they may also help in the understanding of $\gamma$-ray propagation over cosmological distances.

\section{Concluding Remarks}\label{concluding}

{\em Fermi}-LAT has been performing the first all-sky survey in $\gamma$ rays, gathering well-sampled, continuous light curves for hundreds
of AGN and compiling source catalogs for different energy ranges and time
periods. These observations constitute important resources to the astronomical
community for a better understanding of the jet physics, cosmological
evolution, and accretion processes of SMBH. In fact, the {\em Fermi}-LAT
observations, in conjunction with the multi-frequency data collected from
radio to VHE, are key to revealing the nature of jet physics in different
classes of AGN, including particle acceleration, environmental effects, and
interaction processes. In addition, Pass 8 LAT data have increased the
sensitivity for hard-spectrum sources, which are important targets for
ground-based VHE telescopes including the planned Cherenkov Telescope
Array. Moreover, the extended energy range reached by the Pass 8 data opens
new opportunities for the study of blazars at high redshifts. Over the next
years the {\em Fermi} satellite will provide a fundamental contribution in time domain astronomy and multi-messenger/multi-wavelength studies.

\section*{Conflict of Interest Statement}

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

\section*{Author Contributions}

The author confirms being the sole contributor of this work and approved it for publication.

\section*{Acknowledgments}

The \textit{Fermi}-LAT Collaboration acknowledges support for LAT development, operation and data analysis from NASA and DOE (United States), CEA/Irfu and IN2P3/CNRS (France), ASI and INFN (Italy), MEXT, KEK, and JAXA (Japan), and the K.A.~Wallenberg Foundation, the Swedish Research Council and the National Space Board (Sweden). Science analysis support in the operations phase from INAF (Italy) and CNES (France) is also gratefully acknowledged. This work performed in part under DOE Contract DE-AC02-76SF00515.