Search

The time-variable emission from the accretion flow of Sgr A*, the supermassive black hole at the Galactic Center, has long been examined in the radio-to-mm, near-infrared (NIR), and X-ray regimes of the electromagnetic spectrum. However, until now, sensitivity and angular resolution have been insufficient in the crucial mid-infrared (MIR) regime. The MIRI instrument on JWST has changed that, and we report the first MIR detection of Sgr A*. The detection was during a flare that lasted about 40 minutes, a duration similar to NIR and X-ray flares, and the source's spectral index steepened as the flare ended. The steepening suggests synchrotron cooling is an important process for Sgr A*'s variability and implies magnetic field strengths $\sim$40--70 Gauss in the emission zone. Observations at $1.3~\mathrm{mm}$ with the Submillimeter Array revealed a counterpart flare lagging the MIR flare by $\approx$10 minutes. The observations can be self-consistently explained as synchrotron radiation from a single population of gradually cooling high-energy electrons accelerated through (a combination of) magnetic reconnection and/or magnetized turbulence., Comment: Accepted for publication ApJL

We describe the new design and current performance of the Mid-InfraRed Spectrometer and Imager (MIRSI) on the NASA Infrared Telescope Facility (IRTF). The system has been converted from a liquid nitrogen/liquid helium cryogen system to one that uses a closed-cycle cooler, which allows it to be kept on the telescope at operating temperature and available for observing on short notice, requiring less effort by the telescope operators and day crew to maintain operating temperature. Several other enhancements have been completed, including new detector readout electronics, an IRTF-style standard instrument user interface, new stepper motor driver electronics, and an optical camera that views the same field as the mid-IR instrument using a cold dichroic mirror, allowing for guiding and/or simultaneous optical imaging. The instrument performance is presented, both with an engineering-grade array used from 2021-2023, and a science-grade array installed in the fall of 2023. Some sample astronomical results are also shown. The upgraded MIRSI is a facility instrument at the IRTF available to all users., Comment: 15 pages, 17 figures, accepted to PASP

A systematic study, based on the third-moment structure function, of Sgr A*'s variability finds an exponential rise time $\tau_{1,\rm{obs}}=14.8^{+0.4}_{-1.5}~\mathrm{minutes}$ and decay time $\tau_{2,\rm{obs}}=13.1^{+1.3}_{-1.4}~\mathrm{minutes}$. This symmetry of the flux-density variability is consistent with earlier work, and we interpret it as caused by the dominance of Doppler boosting, as opposed to gravitational lensing, in Sgr~A*'s light curve. A relativistic, semi-physical model of Sgr~A* confirms an inclination angle $i<45$ degrees. The model also shows that the emission of the intrinsic radiative process can have some asymmetry even though the observed emission does not. The third-moment structure function, which is a measure of the skewness of the light-curve increments, may be a useful summary statistic in other contexts of astronomy because it senses only temporal asymmetry, i.e., it averages to zero for any temporally symmetric signal., Comment: Accepted for publication in A&A letters

We report on the final two days of a multiwavelength campaign of Sgr A* observing in the radio, submillimeter, infrared, and X-ray bands in July 2019. Sgr A* was remarkably active, showing multiple flaring events across the electromagnetic spectrum. We detect a transient $\sim35$-minute periodicity feature in Spitzer Space Telescope light curves on 21 July 2019. Time-delayed emission was detected in ALMA light curves, suggesting a hotspot within the accretion flow on a stable orbit. On the same night, we observe a decreased flux in the submillimeter light curve following an X-ray flare detected by the Chandra X-ray Observatory and model the feature with an adiabatically expanding synchrotron hotspot occulting the accretion flow. The event is produced by a plasma $0.55~R_{\text{S}}$ in radius with an electron spectrum $p=2.84$. It is threaded by a $\sim130$ Gauss magnetic field and expands at $0.6\%$ the speed of light. Finally, we reveal an unambiguous flare in the infrared, submillimeter, and radio, demonstrating that the variable emission is intrinsically linked. We jointly fit the radio and submillimeter light curves using an adiabatically expanding synchrotron hotspot and find it is produced by a plasma with an electron spectrum $p=0.59$, $187$ Gauss magnetic field, and radius $0.47~R_{\text{S}}$ that expands at $0.029c$. In both cases, the uncertainty in the appropriate lower and upper electron energy bounds may inflate the derived equipartition field strengths by a factor of 2 or more. Our results confirm that both synchrotron- and adiabatic-cooling processes are involved in the variable emission's evolution at submillimeter and infrared wavelengths., Comment: 24 pages, 13 figures, accepted to The Astrophysical Journal. Comments welcome! Paper I can be found here: arXiv:2107.09681

A compact source, G0.02467-0.0727, was detected in ALMA \threemm observations in continuum and very broad line emission. The continuum emission has a spectral index $\alpha\approx3.3$, suggesting that the emission is from dust. The line emission is detected in several transitions of CS, SO, and SO$_2$ and exhibits a line width FWHM $\approx160$ \kms. The line profile appears Gaussian. The emission is weakly spatially resolved, coming from an area on the sky $\lesssim1"$ in diameter ($\lesssim10^4$ AU at the distance of the Galactic Center; GC). The centroid velocity is $v_{LSR}\approx40$-$50$ \kms, which is consistent with a location in the Galactic Center. With multiple SO lines detected, and assuming local thermodynamic equilibrium (LTE) conditions, $T_\mathrm{LTE} = 13$ K, which is colder than seen in typical GC clouds, though we cannot rule out low-density, subthermally excited, warmer gas. Despite the high velocity dispersion, no emission is observed from SiO, suggesting that there are no strong ($\gtrsim10~\mathrm{km~s}^{-1}$) shocks in the molecular gas. There are no detections at other wavelengths, including X-ray, infrared, and radio. We consider several explanations for the Millimeter Ultra-Broad Line Object (MUBLO), including protostellar outflow, explosive outflow, collapsing cloud, evolved star, stellar merger, high-velocity compact cloud, intermediate mass black hole, and background galaxy. Most of these conceptual models are either inconsistent with the data or do not fully explain it. The MUBLO is, at present, an observationally unique object., Comment: Accepted to ApJL

Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace spiral structure within galaxies. Over a dozen of these dense ($\sim$10$^4$\,cm$^{-3}$) and long ($>$10\,pc) filaments have been found within the Milky Way, and they are often referred to as "bones." Until now, none of these bones have had their magnetic field resolved and mapped in their entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping $\sim$10 of these Milky Way bones using the HAWC+ instrument at 214\,$\mu$m and 18$\farcs$2 resolution. Here we present a first result from this survey on the $\sim$60\,pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center-line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel or perpendicular to the Galactic plane nor the bone. The magnetic field strengths along the spine typically vary from $\sim$20 to $\sim$100\,$\mu$G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse., Comment: Accepted to ApJL; typo in author list metadata corrected

We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal' 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments., Comment: 28 pages, 12 figures, 1 table. Accepted for publication in ApJ

The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will have two instrument modules: a 10 - 400 micron hyperspectral imager with spectral resolution R = 8 (GEP-I) and a 24 - 193 micron, R = 200 grating spectrometer (GEP-S). GEP-I surveys will identify star-forming galaxies via their thermal dust emission and simultaneously measure redshifts using polycyclic aromatic hydrocarbon emission lines. Galaxy luminosities derived from star formation and nuclear supermassive black hole accretion will be measured for each source, enabling the cosmic star formation history to be measured to much greater precision than previously possible. Using optically thin far-infrared fine-structure lines, surveys with GEP-S will measure the growth of metallicity in the hearts of galaxies over cosmic time and extraplanar gas will be mapped in spiral galaxies in the local universe to investigate feedback processes. The science case and mission architecture designed to meet the science requirements are described, and the kinetic inductance detector and readout electronics state of the art and needed developments are described. This paper supersedes the GEP concept study report cited in it by providing new content, including: a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys. The reader interested in more technical details may want to consult the concept study report., Comment: 84 pages, 24 figures, published in JATIS

Pre-stellar cores represent the initial conditions of star formation. Although these initial conditions in nearby low-mass star-forming regions have been investigated in detail, such initial conditions remain vastly unexplored for massive star-forming regions. We report the detection of a cluster of low-mass starless and pre-stellar core candidates in a massive star protocluster forming cloud, NGC6334S. With the ALMA observations at a $\sim$0.02 pc spatial resolution, we identified 17 low-mass starless core candidates that do not show any evidence of protostellar activity. These candidates present small velocity dispersions, high fractional abundances of NH$_{2}$D, high NH$_{3}$ deuterium fractionations, and are completely dark in the infrared wavelengths from 3.6 up to 70~$\mu$m. Turbulence is significantly dissipated and the gas kinematics are dominated by thermal motions toward these candidates. Nine out of the 17 cores are gravitationally bound, and therefore are identified as pre-stellar core candidates. The embedded cores of NGC6334S show a wide diversity in masses and evolutionary stages., Comment: 13 pages, 4 figures, accepted for publication in ApJ Letters

It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $\lambda\ga100$ \micron. However, these models do not include AGN heating of host-galaxy-scale diffuse dust, which may have far-IR (FIR) colors comparable to cold diffuse dust heated by stars. To quantify the contribution of AGN heating to host-galaxy-scale cold dust emission at $\lambda\ga100$ \micron, we perform dust RT calculations on a simulated galaxy merger both including and excluding the bright AGN that it hosts. By differencing the SEDs yielded by RT calculations with and without AGN that are otherwise identical, we quantify the FIR cold dust emission arising solely from re-processed AGN photons. In extreme cases, AGN-heated host-galaxy-scale dust can increase galaxy-integrated FIR flux densities by factors of 2-4; star formation rates calculated from the FIR luminosity assuming no AGN contribution can overestimate the true value by comparable factors. Because the FIR colors of such systems are similar to those of purely star-forming galaxies and redder than torus models, broadband SED decomposition may be insufficient for disentangling the contributions of stars and heavily dust-enshrouded AGN in the most IR-luminous galaxies. We demonstrate how kpc-scale resolved observations can be used to identify deeply dust-enshrouded AGN with cool FIR colors when spectroscopic and/or X-ray detection methods are unavailable., Comment: 11 pages, 7 figures, Accepted to ApJ

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC) fit to the X-ray first-order structure function shows significantly less power at short timescales in the X-rays than in the NIR. Less X-ray variability at short timescales combined with the observed NIR-X-ray correlations means the variability can be described as the result of two strictly correlated stochastic processes, the X-ray process being the low-pass-filtered version of the NIR process. The NIR--X-ray linkage suggests a simple radiative model: a compact, self-absorbed synchrotron sphere with high-frequency cutoff close to NIR frequencies plus a synchrotron self-Compton scattering component at higher frequencies. This model, with parameters fit to the submillimeter, NIR, and X-ray structure functions, reproduces the observed flux densities at all wavelengths, the statistical properties of all light curves, and the time lags between bands. The fit also gives reasonable values for physical parameters such as magnetic flux density $B\approx13$ G, source size $L \approx2.2R_{S}$, and high-energy electron density $n_{e}\approx4\times10^{7}$ cm$^{-3}$. An animation illustrates typical light curves, and we make public the parameter chain of our Bayesian analysis, the model implementation, and the visualization code., Comment: Accepted by ApJS, June 3, 2021. The journal version of figure 21 is animated. The animation can also be found here: https://doi.org/10.17617/1.kctx3s25. This version (version 2) has been revised according to the referee's suggestions and includes optimized figures and some changes and corrections of the text; no scientific conclusions have changed

Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus (AGN). We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption/emission in 18 lines of H2O, which we combine with the ALMA H2O 423-330 448 GHz line (Eupper~400 K) and continuum images to study the nuclear region. Radiative transfer models indicate that 3 nuclear components are required to account for the H2O and continuum data. An envelope, with R~130-150 pc, T_dust~50 K, and N_H2~2x10^{23} cm^{-2}, surrounds a nuclear disk with R~40 pc and tau_100um~1.5-3 (N_H2~2x10^{24} cm^{-2}) and an extremely compact (R~12 pc), warm (~100 K), and buried (tau_100um>5, N_H2>~5x10^{24} cm^{-2}) core component. The three nuclear components account for 70% of the galaxy L_IR (SFR~16-18 Msun yr^{-1}). The nucleus is fed by a molecular inflow observed in CO 2-1 with ALMA, which is associated with the nuclear bar. With decreasing radius (r=450-225 pc), the mass inflow rate increases up to ~20 Msun yr^{-1}, which is similar to the nuclear SFR, indicating that the starburst is sustained by the inflow. At lower r, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated inflow rate of ~30 Msun yr^{-1}. The short timescale of ~20 Myr for nuclear gas replenishment indicates quick secular evolution, and indicates that we are witnessing an intermediate stage (<100 Myr) proto-pseudobulge fed by a massive inflow that is driven by a strong nuclear bar. We also apply the H2O model to the Herschel far-infrared spectroscopic observations of H2^{18}O, OH, $^{18}OH, OH+, H2O^+, H3O^+, NH, NH2, NH3, CH, CH^+, ^{13}CH^+, HF, SH, and C3, and estimate their abundances., Comment: Accepted for publication in Astronomy & Astrophysics

Because the 157.74 micron [C II] line is the dominant coolant of star-forming regions, it is often used to infer the global star-formation rates of galaxies. By characterizing the [C II] and far-infrared emission from nearby Galactic star-forming molecular clumps, it is possible to determine whether extragalactic [C II] emission arises from a large ensemble of such clumps, and whether [C II] is indeed a robust indicator of global star formation. We describe [C II] and far-infrared observations using the FIFI-LS instrument on the SOFIA airborne observatory toward four dense, high-mass, Milky Way clumps. Despite similar far-infrared luminosities, the [C II] to far-infrared luminosity ratio, L([C II])/L(FIR) varies by a factor of at least 140 among these four clumps. In particular, for AGAL313.576+0.324, no [C II] line emission is detected despite a FIR luminosity of 24,000 L_sun. AGAL313.576+0.324 lies a factor of more than 100 below the empirical correlation curve between L([C II])/L(FIR) and S_\nu (63 micron)/S_\nu (158 micron) found for galaxies. AGAL313.576+0.324 may be in an early evolutionary "protostellar" phase with insufficient ultraviolet flux to ionize carbon, or in a deeply embedded ``hypercompact' H II region phase where dust attenuation of UV flux limits the region of ionized carbon to undetectably small volumes. Alternatively, its apparent lack of \cii\, emission may arise from deep absorption of the \cii\, line against the 158 micron continuum, or self-absorption of brighter line emission by foreground material, which might cancel or diminish any emission within the FIFI-LS instrument's broad spectral resolution element (~250 km/s)

Emission from active galactic nuclei (AGNs) is known to play an important role in the evolution of many galaxies including luminous and ultraluminous systems (U/LIRGs), as well as merging systems. However, the extent, duration, and exact effects of its influence are still imperfectly understood. To assess the impact of AGNs on interacting systems, we present a Spectral Energy Distribution (SED) analysis of a sample of 189 nearby galaxies. We gather and systematically re-reduce archival broad-band imaging mosaics from the ultraviolet to the far-infrared using data from GALEX, SDSS, 2MASS, IRAS, WISE, Spitzer and Herschel. We use spectroscopy from Spitzer/IRS to obtain fluxes from fine-structure lines that trace star formation and AGN activity. Utilizing the SED modelling and fitting tool CIGALE, we derive the physical conditions of the ISM, both in star-forming regions and in nuclear regions dominated by the AGN in these galaxies. We investigate how the star formation rates (SFRs) and the fractional AGN contributions ($f_{\rm{AGN}}$) depend on stellar mass, galaxy type, and merger stage. We find that luminous galaxies more massive than about $10^{10} \rm{M}_{*}$ are likely to deviate significantly from the conventional galaxy main-sequence relation. Interestingly, infrared AGN luminosity and stellar mass in this set of objects are much tighter than SFR and stellar mass. We find that buried AGNs may occupy a locus between bright starbursts and pure AGNs in the $f_{\rm{AGN}}$-[Ne V]/[Ne II] plane. We identify a modest correlation between $f_{\rm{AGN}}$ and mergers in their later stages., Comment: Accepted for publication in MNRAS; 24 pages, 15 figures, 3 tables (plus appendix)

We present Atacama Large Millimeter/submillimeter Array (ALMA) and Karl G. Jansky Very Large Array (JVLA) observations of the massive infrared dark cloud NGC 6334S (also known as IRDC G350.56+0.44), located at the southwestern end of the NGC 6334 molecular cloud complex. The H$^{13}$CO$^{+}$ and the NH$_{2}$D lines covered by the ALMA observations at a $\sim$3$^{\prime\prime}$ angular resolution ($\sim$0.02 pc) reveal that the spatially unresolved non-thermal motions are predominantly subsonic and transonic, a condition analogous to that found in low-mass star-forming molecular clouds. The observed supersonic non-thermal velocity dispersions in massive star forming regions, often reported in the literature, might be significantly biased by poor spatial resolutions that broaden the observed line widths due to unresolved motions within the telescope beam. Our 3~mm continuum image resolves 49 dense cores, whose masses range from 0.17 to 14 $M_{\odot}$. The majority of them are resolved with multiple velocity components. Our analyses of these gas velocity components find an anti-correlation between the gas mass and the virial parameter. This implies that the more massive structures tend to be more gravitationally unstable. Finally, we find that the external pressure in the NGC 6334S cloud is important in confining these dense structures, and may play a role in the formation of dense cores, and subsequently, the embedded young stars., Comment: 30 pages, 10 figures, 2 tables, accepted for publication in The Astrophysical Journal

The fourth Fermi Large Area Telescope Source Catalog (4FGL) lists over 5000 gamma-ray sources with statistical significance above 4$\sigma$. About 23% of the sources listed in this catalog are unidentified/unassociated gamma-ray sources while ~26% of the sources are classified as blazar candidates of uncertain type (BCUs), lacking optical spectroscopic information. To probe the blazar nature of candidate counterparts of UGSs and BCUs, we started our optical spectroscopic follow up campaign in 2012, which up to date account for more than 350 observed sources. In this paper, the tenth of our campaign, we report on the spectroscopic observations of 37 sources, mostly BCUs, whose observations were carried out predominantly at the Observatorio Astron\'omico Nacional San Pedro M\'artir and the Southern Astrophysical Research Observatory between August 2018 and September 2019. We confirm the BL Lac nature of 27 sources and the flat spectrum radio quasar nature of three sources. The remaining ones are classified as six BL Lacs galaxy-dominated and one normal galaxy. We were also able to measure the redshifts for 20 sources, including 10 BL Lacs. As in previous analyses, the largest fraction of BCUs revealed to be BL Lac objects., Comment: 24 pages, 14 figures, 3 tables (Pre-proof version)

Nearly one-third of the sources in the $Fermi$-LAT catalogs lack a lower energy counterpart, hence being referred as unidentified/unassociated gamma-ray sources (UGSs). In order to firmly classify them, dedicated multifrequency follow-up campaigns are necessary. These will permit to unveil their nature and identify the fraction that could belong to the class of active galaxies known as blazars that is the largest population of extragalactic $\gamma$-ray sources. In $Fermi$-LAT catalogs there are also gamma-ray sources associated with multifrequency blazar-like objects known as Blazars Candidates of Uncertain type (i.e., BCUs) for which follow up spectroscopic campaigns are mandatory to confirm their blazar nature. Thus, in 2013 we started an optical spectroscopic campaign to identify blazar-like objects potential counterparts of UGSs and BCUs. Here we report the spectra of 31 additional targets observed as part of our follow up campaign. Thirteen of them are BCUs for which we acquired spectroscopic observations at Observatorio Astrof\'isico Guillermo Haro (OAGH) and at Southern Astrophysical Research Observatory (SOAR) telescopes, while the rest has been identified thanks to the archival observations available from the Sloan Digital Sky Survey (SDSS). We confirm the blazar nature of all BCUs: three of them are in blazar of quasar type (BZQs) while the remaining ones can be spectroscopically classified as BL Lac objects (BZBs). Then we also discovered 18 BL Lac objects lying within the positional uncertainty regions of UGSs that could be their potential counterparts., Comment: 30 pages and 65 figures

Over the last decade more than five thousand gamma-ray sources were detected by the Large Area Telescope (LAT) on board Fermi Gamma-ray Space Telescope. Given the positional uncertainty of the telescope, nearly 30% of these sources remain without an obvious counterpart in lower energies. This motivated the release of new catalogs of gamma-ray counterpart candidates and several follow up campaigns in the last decade. Recently, two new catalogs of blazar candidates were released, they are the improved and expanded version of the WISE Blazar-Like Radio-Loud Sources (WIBRaLS2) catalog and the Kernel Density Estimation selected candidate BL Lacs (KDEBLLACS) catalog, both selecting blazar-like sources based on their infrared colors from the Wide-field Infrared Survey Explorer (WISE). In this work we characterized these two catalogs, clarifying the true nature of their sources based on their optical spectra from SDSS data release 15, thus testing how efficient they are in selecting true blazars. We first selected all WIBRaLS2 and KDEBLLACS sources with available optical spectra in the footprint of Sloan Digital Sky Survey data release 15. Then we analyzed these spectra to verify the nature of each selected candidate and see which fraction of the catalogs is composed by spectroscopically confirmed blazars. Finally, we evaluated the impact of selection effects, specially those related to optical colors of WIBRaLS2/KDEBLLACS sources and their optical magnitude distributions. We found that at least ~ 30% of each catalog is composed by confirmed blazars, with quasars being the major contaminants in the case of WIBRaLS2 (~ 58%) and normal galaxies in the case of KDEBLLACS (~ 38.2%). The spectral analysis also allowed us to identify the nature of 11 blazar candidates of uncertain type (BCUs) from the Fermi-LAT 4th Point Source Catalog (4FGL) and to find 25 new BL Lac objects., Comment: 11 pages, 11 figures

Thermal infrared observations are the most effective way to measure asteroid diameter and albedo for a large number of near-Earth objects. Major surveys like NEOWISE, NEOSurvey, ExploreNEOs, and NEOLegacy find a small fraction of high albedo objects that do not have clear analogs in the current meteorite population. About 8% of Spitzer-observed near-Earth objects have nominal albedo solutions greater than 0.5. This may be a result of lightcurve variability leading to an incorrect estimate of diameter or inaccurate absolute visual magnitudes. For a sample of 23 high albedo NEOs we do not find that their shapes are significantly different from the McNeill et al. (2019) near-Earth object shape distribution. We performed a Monte Carlo analysis on 1505 near-Earth objects observed by Spitzer, sampling the visible and thermal fluxes of all targets to determine the likelihood of obtaining a high albedo erroneously. Implementing the McNeill shape distribution for near-Earth objects, we provide an upper-limit on the geometric albedo of 0.5+/-0.1 for the near-Earth population., Comment: 19 pages, accepted by the Astronomical Journal

We present two catalogs of radio-loud candidate blazars whose WISE mid-infrared colors are selected to be consistent with the colors of confirmed gamma-ray emitting blazars. The first catalog is the improved and expanded release of the WIBRaLS catalog presented by D'Abrusco et al. (2014): it includes sources detected in all four WISE filters, spatially cross-matched with radio source in one of three radio surveys and radio-loud based on their q22 spectral parameter. WIBRaLS2 includes 9541 sources classified as BL Lacs, FSRQs or mixed candidates based on their WISE colors. The second catalog, called KDEBLLACS, based on a new selection technique, contains 5579 candidate BL Lacs extracted from the population of WISE sources detected in the first three WISE passbands ([3.4], [4.6] and [12]) only, whose mid-infrared colors are similar to those of confirmed, gamma-ray BL Lacs. KDBLLACS members area also required to have a radio counterpart and be radio-loud based on the parameter q12, defined similarly to q22 used for the WIBRaLS2. We describe the properties of these catalogs and compare them with the largest samples of confirmed and candidate blazars in the literature. We crossmatch the two new catalogs with the most recent catalogs of gamma-ray sources detected by Fermi LAT instrument. Since spectroscopic observations of candidate blazars from the first WIBRaLS catalog within the uncertainty regions of gamma-ray unassociated sources confirmed that ~90% of these candidates are blazars, we anticipate that these new catalogs will play again an important role in the identification of the gamma-ray sky., Comment: 20 pages, 7 figures. Accepted for publication in The Astrophysical Journal Supplement Series

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion., Comment: 9 pages, 4 figures, submitted for Astro2020 White Paper

We present lightcurves and derive periods and amplitudes for a subset of 38 near earth objects (NEOs) observed at 4.5 microns with the IRAC camera on the the Spitzer Space Telescope, many of them having no previously reported rotation periods. This subset was chosen from about 1800 IRAC NEO observations as having obvious periodicity and significant amplitude. For objects where the period observed did not sample the full rotational period, we derived lower limits to these parameters based on sinusoidal fits. Lightcurve durations ranged from 42 to 544 minutes, with derived periods from 16 to 400 minutes. We discuss the effects of lightcurve variations on the thermal modeling used to derive diameters and albedos from Spitzer photometry. We find that both diameters and albedos derived from the lightcurve maxima and minima agree with our previously published results, even for extreme objects, showing the conservative nature of the thermal model uncertainties. We also evaluate the NEO rotation rates, sizes, and their cohesive strengths., Comment: 16 pages, 4 figures, 3 tables, to appear in the Astrophysical Journal Supplement Series

Stars form in clustered environments, but how they form when the available resources are shared is still not well understood. A related question is whether the IMF is in fact universal across galactic environments, a galactic initial mass function (IGIMF), or whether it is an average of local IMFs. One of the long-standing problems in resolving this question and in the study of young clusters is observational: the emission from multiple sources is frequently seen as blended because at different wavelengths or with different telescopes the beam sizes are different. The confusion hinders our ability to fully characterize clustered star formation. Here we present a new method that uses a genetic algorithm and Bayesian inference to fit the blended SEDs and images of individual YSOs in confused clusters. We apply this method to the infrared photometry of a sample comprising 70 Spitzer-selected, low-mass ($M_{\rm{cl}}<100~\rm{M}_{\odot}$) young clusters in the galactic plane, and use the derived physical parameters to investigate the distributions of masses and evolutionary stages of clustered YSOs, and the implications of those distributions for studies of the IMF and the different models of star formation. We find that for low-mass clusters composed of class I and class II YSOs, there exists a non-trivial relationship between the total stellar mass of the cluster ($M_{\rm{cl}}$) and the mass of its most massive member ($m_{\rm{max}}$). The properties of the derived correlation are most compatible with the random sampling of a Kroupa IMF, with a fundamental high-mass limit of $150~\rm{M}_{\odot}$. Our results are also compatible with SPH models that predict a dynamical termination of the accretion in protostars, with massive stars undergoing this stopping at later times in their evolution., Comment: 31 pages, 11 figures. Submitted to ApJ. Comments are welcome

Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associated with any star formation indicators, hence are prestellar core candidates. Evolutionary phases of these cores and their linewidths, temperatures, NH$_3$ abundances, and virial parameters are found to be correlated. In a sub-sample of 4 morphologically well-defined filaments, we find that their fragmentation can not be solely explained by thermal or turbulence pressure support. We also investigate distributions of gas temperatures and non-thermal motions along the filaments, and find a spatial correlation between non-thermal linewidths and star formation activities. We find evidence of gas flows along these filaments, and derive an accretion rate along filaments of $\sim$10$^{-4}$ M$_\odot$ yr$^{-1}$. These results suggest a strong relationship between massive filaments and high-mass star formation, through i) filamentary fragmentation in very early evolutionary phases to form dense cores, ii) accretion flows along filaments that are important for the growth of dense cores and protostars, and iii) enhancement of non-thermal motion in the filaments by the feedback or accretion during star formation., Comment: 31 pages, 11 figures, 8 tables, ApJ Accepted. 3-D interactive figures to visualize outflows (see Figure 7) will be published online, and are available at https://xinglunju.github.io/outflows.html

Galaxy mergers are key events in galaxy evolution, often causing massive starbursts and fueling active galactic nuclei (AGN). In these highly dynamic systems, it is not yet precisely known how much starbursts and AGN respectively contribute to the total luminosity, at what interaction stages they occur, and how long they persist. Here we estimate the fraction of the bolometric infrared (IR) luminosity that can be attributed to AGN by measuring and modeling the full ultraviolet to far-infrared spectral energy distributions (SEDs) in up to 33 broad bands for 24 merging galaxies with the Code for Investigating Galaxy Emission. In addition to a sample of 12 confirmed AGN in late-stage mergers, found in the $Infrared$ $Array$ $Satellite$ Revised Bright Galaxy Sample or Faint Source Catalog, our sample includes a comparison sample of 12 galaxy mergers from the $Spitzer$ Interacting Galaxies Survey, mostly early-stage. We perform identical SED modeling of simulated mergers to validate our methods, and we supplement the SED data with mid-IR spectra of diagnostic lines obtained with $Spitzer$ InfraRed Spectrograph. The estimated AGN contributions to the IR luminosities vary from system to system from 0% up to 91% but are significantly greater in the later-stage, more luminous mergers, consistent with what is known about galaxy evolution and AGN triggering., Comment: 26 pages, 10 figures

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars., Comment: 75 pages. See page 3 for Table of Contents and page 4 for Executive Summary

The existence of "radio weak BL Lac objects" (RWBLs) has been an open question, still unsolved, since the discovery that quasars could be radio-quiet or radio-loud. Recently several groups identified RWBL candidates, mostly found while searching for low energy counterparts of the unidentified/unassociated gamma-ray sources listed in the Fermi catalogs. Confirming RWBLs is a challenging task since they could be confused with white dwarfs (WDs) or weak emission line quasars (WELQs) when there are not sufficient data to precisely draw their broad band spectral energy distribution and their classification is mainly based on a featureless optical spectra. Motivated by the recent discovery that Fermi BL Lacs appear to have very peculiar mid-IR emission, we show that it is possible to distinguish between WDs, WELQs and BL Lacs using the [3.4]-[4.6]-[12]$\mu$m color-color plot built using the WISE magnitudes when the optical spectrum is available. On the basis of this analysis, we identify WISE J064459.38+603131 and WISE J141046.00+740511.2 as the first two genuine RWBLs, both potentially associated with Fermi sources. Finally, to strengthen our identification of these objects as true RWBLs, we present multifrequency observations for these two candidates to show that their spectral behavior is indeed consistent with those of the BL Lac population., Comment: 12 pages, 9 figures, ApJ accepted for publication, pre-proof version

We recently developed a procedure to recognize gamma-ray blazar candidates within the positional uncertainty regions of the unidentified/unassociated gamma-ray sources (UGSs). Such procedure was based on the discovery that Fermi blazars show peculiar infrared colors. However, to confirm the real nature of the selected candidates, optical spectroscopic data are necessary. Thus, we performed an extensive archival search for spectra available in the literature in parallel with an optical spectroscopic campaign aimed to reveal and confirm the nature of the selected gamma-ray blazar candidates. Here, we first search for optical spectra of a selected sample of gamma-ray blazar candidates that can be potential counterparts of UGSs using the Sloan Digital Sky Survey (SDSS DR12). This search enables us to update the archival search carried out to date. We also describe the state-of-art and the future perspectives of our campaign to discover previously unknown gamma-ray blazars., Comment: 11 pages, 2 figures, 5 tables, pre-proof version, accepted for publication of Astrophysics and Space Science

Despite the fact that blazars constitute the rarest class among active galactic nuclei (AGNs) they are the largest known population of associated $\gamma$-ray sources. Many of the $\gamma$-ray objects listed in the Fermi-Large Area Telescope Third Source catalog (3FGL) are classified as blazar candidates of uncertain type (BCUs), either because they show multifrequency behaviour similar to blazars but lacking optical spectra in the literature, or because the quality of such spectra is too low to confirm their nature. Here we select, out of 585 BCUs in the 3FGL, 42 BCUs which we identify as probable blazars by their WISE infrared colors and which also have optical spectra that are available in the Sloan Digital Sky Survey (SDSS) and/or Six-Degree Field Galaxy Survey Database (6dFGS). We confirm the blazar nature of all of the sources. We furthermore conclude that 28 of them are BL Lacs, 8 are radio-loud quasars with flat radio spectrum and 6 are BL Lac whose emission is dominated by their host galaxy.

Blazars, one of the most extreme class of active galaxies, constitute so far the largest known population of $\gamma$-ray sources and their number is continuously growing in the Fermi catalogs. However in the latest release of the Fermi catalog there is still a large fraction of sources that are classified as blazar candidate of uncertain type (BCUs) for which optical spectroscopic observations are necessary to confirm their nature and their associations. In addition about 1/3 of the $\gamma$-ray point sources listed in the Third Fermi-LAT Source Catalog (3FGL) are still unassociated and lacking an assigned lower energy counterpart. Since 2012 we have been carrying out an optical spectroscopic campaign to observe blazar candidates to confirm their nature. In this paper, the sixth of the series, we present optical spectroscopic observations for 30 $\gamma$-ray blazar candidates from different observing programs we carried out with the TNG, WHT, OAN, SOAR and Magellan telescopes. We found that 21 out of 30 sources investigated are BL Lac objects while the remaining targets are classified as flat spectrum radio quasars showing the typical broad emission lines of normal quasi stellar objects. We conclude that our selection of $\gamma$-ray blazars candidates based on their multifrequency properties continues to be a successful way to discover potential low-energy counterparts of the FermiUnidentified Gamma-ray Sources (UGSs) and to confirm the nature of BCUs.

The extragalactic $\gamma$-ray sky is dominated by emission from blazars, a peculiar class of active galactic nuclei (AGNs). Many of the $\gamma$-ray sources included in Fermi -Large Area Telescope Third Source catalog (3FGL) are classified as a blazar candidate of uncertain type (BCU) because there is no optical spectra available in the literature to confirm their nature. In 2013 we started a spectroscopic campaign to look for the optical counterparts of the BCUs and of the Unidentified $\gamma$-ray Sources. The main goal of our investigation is to confirm the blazar nature of these sources having peculiar properties as compact radio emission and/or selected on the basis of their infrared (IR) colors. Whenever possible we also determine their redshifts. Here we present the results of the observations carried out in the Northern hemisphere in 2013 and 2014 at Telescopio Nazionale Galilleo (TNG), Kitt Peak National Observatory (KPNO) and Observatorio Astron\'omico Nacional (OAN) in San Pedro M\'artir. In this paper we describe the optical spectra of 25 sources. We confirmed that all the 15 BCUs observed in our campaign and included in our sample are blazars and we estimated the redshift for 3 of them. In addition, we present the spectra for 3 sources classified as BL Lacs in the literature but with no optical spectra available to date. We found that one of them is a quasar (QSO) at a redshift $z$ = 0.208 and the other 2 are BL Lacs. Moreover, we also present 7 new spectra for known blazars listed in the Roma-BZCAT having an uncertain redshift or being classified as BL Lac candidates. We found that one of them, 5BZB J0724+2621 is a `changing look' blazar. According to the spectrum available in the literature it was classified as a BL Lac but in our observation we clearly detected a broad emission line that lead to classify this source as a QSO at $z$=1.17., Comment: 21 pages, 26 figures, 3 tables, pre-proof version, published on the Astronomical Journal. arXiv admin note: text overlap with arXiv:1503.04202

Near Earth Objects (NEOs) are small Solar System bodies whose orbits bring them close to the Earth's orbit. We are carrying out a Warm Spitzer Cycle 11 Exploration Science program entitled NEOSurvey --- a fast and efficient flux-limited survey of 597 known NEOs in which we derive diameter and albedo for each target. The vast majority of our targets are too faint to be observed by NEOWISE, though a small sample has been or will be observed by both observatories, which allows for a cross-check of our mutual results. Our primary goal is to create a large and uniform catalog of NEO properties. We present here the first results from this new program: fluxes and derived diameters and albedos for 80 NEOs, together with a description of the overall program and approach, including several updates to our thermal model. The largest source of error in our diameter and albedo solutions, which derive from our single band thermal emission measurements, is uncertainty in eta, the beaming parameter used in our thermal modeling; for albedos, improvements in Solar System absolute magnitudes would also help significantly. All data and derived diameters and albedos from this entire program are being posted on a publicly accessible webpage at nearearthobjects.nau.edu ., Comment: AJ in press

Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion., Comment: 13 pages, 8 figures, ApJ accepted

We present dust column densities and dust temperatures for $\sim3000$ young high-mass molecular clumps from the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey, derived from adjusting single temperature dust emission models to the far-infrared intensity maps measured between 160 and 870 \micron\ from the Herschel/Hi-Gal and APEX/ATLASGAL surveys. We discuss the methodology employed in analyzing the data, calculating physical parameters, and estimating their uncertainties. The population average dust temperature of the clumps are: $16.8\pm0.2$ K for the clumps that do not exhibit mid-infrared signatures of star formation (Quiescent clumps), $18.6\pm0.2$ K for the clumps that display mid-infrared signatures of ongoing star formation but have not yet developed an HII region (Protostellar clumps), and $23.7\pm0.2$ and $28.1\pm0.3$ K for clumps associated with HII and photo-dissociation regions, respectively. These four groups exhibit large overlaps in their temperature distributions, with dispersions ranging between 4 and 6 K. The median of the peak column densities of the Protostellar clump population is $0.20\pm0.02$ gr cm$^{-2}$, which is about 50% higher compared to the median of the peak column densities associated with clumps in the other evolutionary stages. We compare the dust temperatures and column densities measured toward the center of the clumps with the mean values of each clump. We find that in the Quiescent clumps the dust temperature increases toward the outer regions and that they are associated with the shallowest column density profiles. In contrast, molecular clumps in the Protostellar or HII region phase have dust temperature gradients more consistent with internal heating and are associated with steeper column density profiles compared with the Quiescent clumps., Comment: 41 pages, 9 figures, accepted in ApJ

The Fermi-Large Area Telescope (LAT) First Source Catalog (1FGL) was released in February 2010 and the Fermi-LAT 2-Year Source Catalog (2FGL) appeared in April 2012, based on data from 24 months of operation. Since their releases, many follow up observations of unidentified gamma-ray sources (UGSs) were performed and new procedures to associate gamma-ray sources with potential counterparts at other wavelengths were developed. Here we review and characterize all the associations as published in the 1FGL and 2FGL catalog on the basis of multifrequency archival observations. In particular we located 177 spectra for the low-energy counterparts that were not listed in the previous Fermi catalogs, and in addition we present new spectroscopic observations of 8 gamma-ray blazar candidates. Based on our investigations, we introduce a new counterpart category of "candidate associations" and propose a refined classification for the candidate low-energy counterparts of the Fermi sources. We compare the 1FGL-assigned counterparts with those listed in the 2FGL to determine which unassociated sources became associated in later releases of the Fermi catalogs. We also search for potential counterparts to all the remaining unassociated Fermi sources. Finally, we prepare a refined and merged list of all the associations of the 1FGL plus 2FGL catalogs that includes 2219 unique Fermi objects. This is the most comprehensive and systematic study of all the associations collected for the gamma-ray sources available to date. We conclude that 80% of the Fermi sources have at least one known plausible gamma-ray emitter within their positional uncertainty regions., Comment: 26 pages, 24 figures, 7 tables, ApJS accepted for publication (pre-proof version uploaded)

The extragalactic gamma-ray sky is dominated by the emission arising from blazars, one of the most peculiar classes of radio-loud active galaxies. Since the launch of Fermi several methods were developed to search for blazars as potential counterparts of unidentified gamma-ray sources (UGSs). To confirm the nature of the selected candidates, optical spectroscopic observations are necessary. In 2013 we started a spectroscopic campaign to investigate gamma-ray blazar candidates selected according to different procedures. The main goals of our campaign are: 1) to confirm the nature of these candidates, and 2) whenever possible determine their redshifts. Optical spectroscopic observations will also permit us to verify the robustness of the proposed associations and check for the presence of possible source class contaminants to our counterpart selection. This paper reports the results of observations carried out in 2014 in the Northern hemisphere with Kitt Peak National Observatory (KPNO) and in the Southern hemisphere with the Southern Astrophysical Research (SOAR) telescopes. We also report three sources observed with the Magellan and Palomar telescopes. Our selection of blazar-like sources that could be potential counterparts of UGSs is based on their peculiar IR colors and on their combination with radio observations both at high and low frequencies (i.e., above and below ~ 1 GHz) in publicly available large radio surveys. We present the optical spectra of 27 objects. [....] We also report the case for WISE J173052.85-035247.2, candidate counterpart of the source 2FGL J1730.6-0353, which has no radio counterpart in the major radio surveys. We confirm that our selection of gamma-ray blazars candidates can successfully indentify low-energy counterparts to Fermi unassociated sources and allow us to discover new blazars., Comment: 20 pages, 30 figures, 1 table, AJ accepted for publication (pre-proof version), Abstract cut

The classification of galaxy mergers and isolated disks is key for understanding the relative importance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic properties of galaxies as traced by emission lines have been used to suggest the existence of a significant population of high-z star-forming galaxies consistent with isolated rotating disks. However, recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. To further investigate the robustness of merger/disk classifications based on kinematic properties, we carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially redshift the WiFeS observations of these local (U)LIRGs to z=1.5 to make a realistic comparison with observations at high-z, and also to ensure that all galaxies have the same spatial sampling of ~900 pc. Using both kinemetry-based and visual classifications, we find that the reliability of kinematic classification shows a strong trend with the interaction stage of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties compared to isolated disks, whereas a significant population of the interacting disks and merger remnants are indistinguishable from isolated disks. The high fraction of late-stage mergers showing disk-like kinematics reflects the complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified disks and mergers depend on the definition of kinematic asymmetries and the classification threshold when using kinemetry-based classifications. Our results suggest that additional indicators such as morphologies traced by stars or molecular gas are required to further constrain the merger/disk classifications at high-z., Comment: 16 pages, 5 figures, ApJ accepted

We report the results of our exploratory program carried out with the Southern Astrophysical Research (SOAR) telescope aimed at associating counterparts and establishing the nature of the Fermi Unidentified gamma-ray Sources (UGS). We selected the optical counterparts of 6 UGSs from the Fermi catalog on the basis of our recently discovered tight connection between infrared and gamma-ray emission found for the gamma-ray blazars detected by the Wide-Field Infrared Survey Explorer (WISE) in its the all-sky survey. We perform for the first time a spectroscopic study of the low-energy counterparts of Fermi UGS, in the optical band, confirming the blazar-like nature for the whole sample. We also present new spectroscopic observations of 6 Active Galaxies of Uncertain type associated with Fermi sources (AGUs) that appear to be BL Lac objects. Finally, we report the spectra collected for 6 known gamma-ray blazars belonging to the Roma BZCAT that were obtained to establish their nature or better estimate their redshifts. Two interesting cases of high redshift and extremely luminous BL Lac objects (z>1.18 and z>1.02, based on the detection of Mg II intervening systems) are also discussed., Comment: 41 pages (single column), 22 figures, 1 table, AJ accepted for publication (pre-proof version)

We recently started a systematic search of low-energy counterparts of the unidentified gamma-ray sources (UGSs) listed in the Fermi-Large Area Telescope (LAT) First Source Catalog (1FGL) and the Fermi-LAT 2-Year Source Catalog (2FGL).} The main goal of our investigation is to find active galaxies belonging to the blazar class that lie within the positional uncertainty region of the UGSs and thus could be their potential low-energy counterparts. To achieve our aims, we first adopted several procedures based on the peculiar observational properties of blazars in the radio and in the IR. Then we carried out a follow-up spectroscopic campaign in the optical band to verify the nature of the candidates selected as potential counterparts of the UGSs. Here we present the results of the observations carried out in 2013 in the Northern Hemisphere at Kitt Peak National Observatory (KPNO). Optical spectroscopy is crucial to confirm the nature of the sources and can be used to estimate their redshifts; it will also allow us to test the robustness of our methods when the whole campaign is completed. Here we present the optical spectroscopic observations of 39 sources. Within our sample we found that 6 sources are blazars, candidates to be low-energy counterparts of the UGSs listed in the 2FGL. We confirm that an additional 8 sources, previously classified as active galaxies of uncertain type and associated in the 2FGL, are also all BL Lac objects. Moreover, we also present 20 new spectra for known blazars listed in the Multi-frequency Catalogue of Blazars as having an uncertain redshift and/or being classified as BL Lac candidates. We conclude that our methods for selecting gamma-ray blazar candidates allows us to discover new blazars and increase the list of potential low-energy counterparts for the Fermi UGSs., Comment: 27 pages, 39 figures, 1 table, A&A accepted for publication (pre-proof version)