X-ray binaries exhibit a wide range of properties but there are few accepted methods to determine the nature of the compact object. Color-Color-Intensity diagrams have been suggested as a means of distinguishing between systems containing black holes from those containing neutron stars. However, this technique has been verified with data from only one instrument (RXTE/ASM) with a single set of X-ray colors defined using data available only in pre-determined energy bands. We test a selection of X-ray colors with a more sensitive instrument to determine the reliability of this method. We use data from the MAXI Gas Slit Camera, which allows users to specify energy-bands. We test X-ray colors that have been previously defined in the literature as well as ones that we define specifically in this paper. A representative set of systems are used to construct Color-Color-Intensity diagrams in each set of colors to determine which are best for separating different classes. For studying individual sources certain bands are more effective than others. For a specified energy range, the separation of soft states in black hole binaries was possible only where both soft and hard colors included information from the lowest energy band. We confirm that Color-Color-Intensity diagrams can distinguish between systems containing black holes or neutron stars in all X-ray colors tested; this suggests an universality in the accretion processes governing these different classes. We suggest possible physical processes driving different classes of X-ray binaries to different locations in Color-Color-Intensity diagrams.

Nearly 50% of all sources detected by the Fermi Large Area Telescope are classified as blazars or blazar candidates, one of the most elusive classes of active galaxies. Additional blazars can also be hidden within the sample of unidentified or unassociated $\gamma$-ray sources (UGSs) that constitute about one-third of all gamma-ray sources detected to date. We recently confirmed that the large majority of Fermi blazars of the BL Lac subclass have an X-ray counterpart. Using the X-ray properties of a BL Lac training set and combining these with archival multifrequency information, we aim to search for UGSs that could have a BL Lac source within their $\gamma$-ray positional uncertainty regions. We reduced and analyzed the Swift X-ray observations of a selected sample of 327 UGSs. We then compared the X-ray fluxes and hardness ratios of all sources detected in the pointed fields with those of known Fermi BL Lacs. We find at least one X-ray source, lying within the $\gamma$-ray positional uncertainty at 95% confidence level, for 223 UGSs and a total of 464 X-ray sources in all fields analyzed. The X-ray properties of a large fraction of them, eventually combined with radio, infrared, and optical information, exhibit BL Lac multi-frequency behavior, thus allowing us to select high-confidence BL Lac candidates; some of them were recently observed during our optical spectroscopic campaign which confirmed their nature. We find that out of 50 X-ray sources that were confirmed as BL Lacs through optical spectroscopy, 12 do not show canonical mid-infrared or radio BL Lac properties. This indicates that the selection of X-ray BL Lac candidates is a strong method to find new counterparts within Fermi UGSs. Finally, we pinpoint a sample of 32 Swift/XRT candidate counterparts to Fermi UGSs that are most likely BL Lac objects.

This paper continues the series of the Fornax Deep Survey (FDS). Following the previous studies dedicated to extended Fornax cluster members, we present the catalogs of compact stellar systems in the Fornax cluster as well as extended background sources and point-like sources. We derive ugri photometry of ~1.7 million sources over the $\sim$21 sq. degree area of FDS centered on NGC1399. For a wider area, of $\sim$27 sq. degs extending in the direction of NGC1316, we provide gri data for ~3.1 million sources. To improve the morphological characterization of sources we generate multi-band image stacks by coadding the best seeing gri-band single exposures with a cut at FWHM<=0.9 arcsec. We use the multi-band stacks as detection frames. The identification of compact sources is obtained from a combination of photometric and morphometric selection criteria taking as reference the properties of sources with well-defined classification from the literature. We present a preliminary analysis of globular cluster (GC) distributions in the Fornax area. The study confirms and extends further previous results. We observe the inter-galactic population of GCs, a population of mainly blue GCs centered on NGC1399, extends over $\sim$0.9Mpc, with an ellipticity $\sim$0.65. Several sub-structures extend over $\sim$0.5Mpc along various directions. Two of these structures do not cross any bright galaxy; one of them appears to be connected to NGC1404, a bright galaxy close to the cluster core and particularly poor of GCs. Using the gri catalogs we analyze the GC distribution over the extended FDS area, and do not find any obvious GC sub-structure bridging the two brightest cluster galaxies, NGC1316 and NGC1399. Although NGC1316 is twice brighter of NGC1399 in optical bands we estimate a factor of 3-4 richer GC population around NGC1399 compared to NGC1316, out to galactocentric distances of 40 arcmin

Associating gamma-ray sources to their low-energy counterparts is one of the major challenges of modern gamma-ray astronomy. In the context of the Fourth Fermi Large Area Telescope Source Catalog (4FGL), the associations rely mainly on parameters as apparent magnitude, integrated flux, and angular separation between the gamma-ray source and its low-energy candidate counterpart. In this work we propose a new use of likelihood ratio and a complementary supervised learning technique to associate gamma-ray blazars in 4FGL, based only on spectral parameters as gamma-ray photon index, mid-infrared colors and radio-loudness. In the likelihood ratio approach, we crossmatch the WISE Blazar-Like Radio-Loud Sources catalog with 4FGL and compare the resulting candidate counterparts with the sources listed in the gamma-ray blazar locus to compute an association probability for 1138 counterparts. In the supervised learning approach, we train a random forest algorithm with 869 high confidence blazar associations and 711 fake associations, and then compute an association probability for 1311 candidate counterparts. A list with all 4FGL blazar candidates of uncertain type associated by our method is provided to guide future optical spectroscopic follow up observations.

BL Lac objects are an extreme type of active galactic nuclei (AGNs) that belong to the largest population of $\gamma$-ray sources: blazars. This class of AGNs shows a double-bumped spectral energy distribution that is commonly described in terms of a synchrotron self-Compton (SSC) emission process, whereas the low-energy component that dominates their emission between the infrared and the X-ray band is tightly connected to the high-energy component that peaks in the $\gamma$-rays. Two strong connections that link radio and mid-infrared emission of blazars to the emission in the $\gamma$-ray band are well established. They constitute the basis for associating $\gamma$-ray sources with their low-energy counterparts. We searched for a possible link between X-ray and $\gamma$-ray emissions for the subclass of BL Lacs using all archival Swift/XRT observations combined with Fermi data for a selected sample of 351 sources. Analyzing $\sim$2400 ks of Swift/XRT observations that were carried out until December 2018, we discovered that above the $\gamma$-ray flux threshold $F_{\gamma}\approx3\times10^{-12}\,\rm{erg}\,\rm{cm}^{-2}\,\rm{s}^{-1}$, 96\% of all \emph{Fermi} BL Lacs have an X-ray counterpart that is detected with signal-to-noise ratio higher than 3. We did not find any correlation or clear trend between X-ray and $\gamma$-ray fluxes and/or spectral shapes, but we discovered a correlation between the X-ray flux and the mid-infrared color. Finally, we discuss on a possible interpretation of our results in the SSC framework.

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.

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.

The importance of archival science increases significantly for astrophysical observatories as they mature and their archive holdings grow in size and complexity. Further enhancing the science return of archival data requires engaging a larger audience than the mission reference community, mostly because of the growth of interest in multi-wavelength and transient/time variability research. Such a goal, though, can be difficult to achieve. In this paper I will describe a different approach to this question that, while minimizing technological friction and leveraging existing services, makes archival observations more accessible and increases our capability to proactively engage astronomers on potentially interesting archival records. Inspired by this strategy, the Chandra Data Archive team is working on two specific experimental projects that will hopefully demonstrate their potential while contributing to the maximization of the scientific return of the Chandra mission.

Using the data products of the Chandra Galaxy Atlas (Kim et al. 2019a), we have investigated the radial profiles of the hot gas temperature in 60 early type galaxies. Considering the characteristic temperature and radius of the peak, dip, and break (when scaled by the gas temperature and virial radius of each galaxy), we propose a universal temperature profile of the hot halo in ETGs. In this scheme, the hot gas temperature peaks at RMAX = 35 +/- 25 kpc (or ~0.04 RVIR) and declines both inward and outward. The temperature dips (or breaks) at RMIN (or RBREAK) = 3 - 5 kpc (or ~0.006 RVIR). The mean slope between RMIN (RBREAK) and RMAX is 0.3 +/- 0.1. Allowing for selection effects and observational limits, we find that the universal temperature profile can describe the temperature profiles of 72% (possibly up to 82%) of our ETG sample. The remaining ETGs (18%) with irregular or monotonically declining profiles do not fit the universal profile and require another explanation. The temperature gradient inside RMIN (RBREAK) varies widely, indicating different degrees of additional heating at small radii. Investigating the nature of the hot core (HC with a negative gradient inside RMIN), we find that HC is most clearly visible in small galaxies. Searching for potential clues associated with stellar, AGN feedback, and gravitational heating, we find that HC may be related to recent star formation. But we see no clear evidence that AGN feedback and gravitational heating play any significant role for HC.

In the last years, Astroinformatics has become a well defined paradigm for many fields of Astronomy. In this work we demonstrate the potential of a multidisciplinary approach to identify globular clusters (GCs) in the Fornax cluster of galaxies taking advantage of multi-band photometry produced by the VLT Survey Telescope using automatic self-adaptive methodologies. The data analyzed in this work consist of deep, multi-band, partially overlapping images centered on the core of the Fornax cluster. In this work we use a Neural-Gas model, a pure clustering machine learning methodology, to approach the GC detection, while a novel feature selection method ($\Phi$LAB) is exploited to perform the parameter space analysis and optimization. We demonstrate that the use of an Astroinformatics based methodology is able to provide GC samples that are comparable, in terms of purity and completeness with those obtained using single band HST data (Brescia et al. 2012) and two approaches based respectively on a morpho-photometric (Cantiello et al. 2018b) and a PCA analysis (D'Abrusco et al. 2015) using the same data discussed in this work.

The formation of Low mass X-ray binaries (LMXB) is favored within dense stellar systems such as Globular Clusters (GCs). The connection between LMXB and Globular Clusters has been extensively studied in the literature, but these studies have always been restricted to the innermost regions of galaxies. We present a study of LMXB in GCs within the central 1.5 deg^2 of the Fornax cluster with the aim of confirming the existence of a population of LMXB in intra-cluster GCs and understand if their properties are related to the host GCs, to the environment or/and to different formation channels.

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.

We present a bibliographic analysis of Chandra, Hubble, and Spitzer publications. We find (a) archival data are used in >60% of the publication output and (b) archives for these missions enable a much broader set of institutions and countries to scientifically use data from these missions. Specifically, we find that authors from institutions that have published few papers from a given mission publish 2/3 archival publications, while those with many publications typically have 1/3 archival publications. We also show that countries with lower GDP per capita overwhelmingly produce archival publications, while countries with higher GDP per capital produce guest observer and archival publications in equal amounts. We argue that robust archives are thus not only critical for the scientific productivity of mission data, but also the scientific accessibility of mission data. We argue that the astronomical community should support archives to maximize the overall scientific societal impact of astronomy, and represent an excellent investment in astronomy's future.

The hot ISM in early type galaxies (ETGs) plays a crucial role in understanding their formation and evolution. The structural features of the hot gas identified by Chandra observations point to key evolutionary mechanisms, (e.g., AGN and stellar feedback, merging history). In our Chandra Galaxy Atlas (CGA) project, taking full advantage of the Chandra capabilities, we systematically analyzed the archival Chandra data of 70 ETGs and produced uniform data products for the hot gas properties. The primary data products are spatially resolved 2D spectral maps of the hot gas from individual galaxies. We emphasize that new features can be identified in the spectral maps which are not readily visible in the surface brightness maps. The high-level images can be viewed at the dedicated CGA website, and the CGA data products can be downloaded to compare with data at other wavelengths and to perform further analyses. Utilizing our data products, we address a few focused science topics.

The Fornax Deep Survey Dwarf galaxy Catalog (FDSDC) includes 564 dwarf galaxies in the Fornax cluster and the in-falling Fornax A subgroup. We use the FDSDC galaxies for statistical analysis of the structural and stellar population differences in the range of galactic environments within the Fornax cluster. We present the standard scaling relations for the dwarfs and analyze trends as a function of cluster-centric radius. We find a different behavior for the bright dwarfs (-18.5 mag < M$_r$ < -16 mag) as compared to the fainter ones (M$_r$ > -16 mag): While considering galaxies in the same magnitude-bins, we find that, while for fainter dwarfs the g'-r' color is redder for lower surface brightness objects (as expected from fading stellar populations), for brighter dwarfs the color is redder for the higher surface brightness and higher S\'ersic n objects. The trend of the bright dwarfs might be explained by those galaxies being affected by harassment and by slower quenching of star formation in their inner parts. As the fraction of early-type dwarfs with respect to late-types increases toward the central parts of the cluster, the color-surface brightness trends are also manifested in the cluster-centric trends, confirming that it is indeed the environment that changes the galaxies. We also estimate the strengths of the ram-pressure stripping, tidal disruption, and harassment in the Fornax cluster, and find that our observations are consistent with the theoretically expected ranges of galaxy properties where each of those mechanisms dominate. We furthermore find that the luminosity function, color-magnitude relation, and axis-ratio distribution of the dwarfs in the center of the Fornax cluster are similar to those in the center of the Virgo cluster.

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.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

This paper is based on the multi-band (ugri) Fornax Deep Survey (FDS) with the VLT Survey Telescope (VST). We study bright early-type galaxies (m_B<15 mag) in the 9 square degrees around the core of the Fornax cluster, which covers the virial radius. The main goal of the present work is to provide the analysis of the light distribution for all galaxies out to unprecedented limits (in radius and surface brightness) and to release the main products resulting from this analysis in all FDS bands. From the isophote fit, we derive the azimuthally averaged surface brightness profiles, the position angle and ellipticity profiles as a function of the semi-major axis. In each band, we derive the total magnitudes, effective radii, integrated colors and stellar mass-to-light ratios. The long integration times, the arcsec-level angular resolution of OmegaCam@VST and the large covered area of FDS allow us to map the light and color distributions out to large galactocentric distances (up to about 10-15 R_e) and surface brightness levels beyond mu_r = 27 mag/arcsec^2 (mu_B > 28 mag/arcsec^2). Therefore, the new FDS data allow us to explore in great detail the morphology and structure of cluster galaxies out to the region of the stellar halo. The observations suggest that the Fornax cluster is not completely relaxed inside the virial radius. The bulk of the gravitational interactions between galaxies happens in the W-NW core region of the cluster, where most of the bright early-type galaxies are located and where the intra-cluster baryons (diffuse light and GCs) are found. We suggest that the W-NW sub-clump of galaxies results from an infalling group onto the cluster, which has modified the structure of the galaxy outskirts (making asymmetric stellar halos) and has produced the intra-cluster baryons (ICL and GCs), concentrated in this region of the cluster.

The Fornax Deep Survey (FDS), an imaging survey in the u', g', r', and i'-bands, has a supreme resolution and image depth compared to the previous spatially complete Fornax Cluster Catalog (FCC). Our new data allows us to study the galaxies down to r'-band magnitude m$_{r'}\approx$21 mag (M$_{r'}\approx$-10.5 mag). These data provide an important legacy dataset to study the Fornax cluster. We aim to present the Fornax Deep Survey (FDS) dwarf galaxy catalog, focusing on explaining the data reduction and calibrations, assessing the quality of the data, and describing the methods used for defining the cluster memberships for the catalog objects. As a first step we used the SExtractor fine-tuned for dwarf galaxy detection, to find galaxies from the FDS data, covering a 26 deg$^2$ area of the main cluster, and the area around the Fornax A substructure. We made 2D-decompositions of the identified galaxies using GALFIT. We used color-magnitude, luminosity-radius and luminosity-concentration relations to separate the cluster galaxies from the background galaxies. We then divided the cluster galaxies into early- and late-type galaxies according to their morphology and gave first order morphological classifications. Our final catalog includes 14,095 galaxies. We classify 590 galaxies as being likely Fornax cluster galaxies, of which 564 are dwarfs (M$_{r'}$ > -18.5 mag) consisting our Fornax dwarf catalog. Of the cluster dwarfs we classify 470 as early-types, and 94 as late-type galaxies. Our final catalog reaches its 50% completeness limit at magnitude M$_{r'}$ = -10.5 mag and surface brightness $\bar{\mu}_{e,r'}$ = 26 mag arcsec-2, which is approximately three magnitudes deeper than the FCC. Based on previous works and comparison with a spectroscopically confirmed subsample, we estimate that our final Fornax dwarf galaxy catalog has < 10% contamination from the background objects.

We present the results of a wide spectroscopic survey aimed at detecting extragalactic globular clusters (GCs) in the core of the Fornax cluster. About 4500 low resolution spectra (from 4800 to 10000 \AA ) were observed in 25 VLT/VIMOS masks covering the central 1 deg$^{2}$ around the dominant galaxy NGC 1399 corresponding to $\sim$175 kpc galactocentric radius. We describe the methodology used for data reduction and data analysis. We found a total of 387 unique physical objects (372 GCs and 15 ultra compact dwarfs) in the field covered by our observations. Most of these objects are associated with NGC 1399, with only 10% likely belonging to other giant galaxies. The new VIMOS dataset is complementary to the many GC catalogues already present in the literature and it brings the total number of tracer particles around NGC 1399 to more than 1130 objects. With this comprehensive radial velocity sample we have found that the velocity dispersion of the GC population (equally for red and blue GC populations) shows a relatively sharp increase from low velocity dispersion ($\sim250$-$350$ kms$^{-1}$) to high velocity dispersion ($\sim300$-$400$ kms$^{-1}$) at projected radius of $\approx10$ arcmin ($\sim60$ kpc) from the galaxy centre. This suggests that at a projected radius of $\approx60$ kpc both blue and red GC populations begin to be governed by the dominating Fornax cluster potential, rather than by the central NGC 1399 galaxy potential. This kinematic evidence corroborates similar results found using surface brightness analysis and planetary nebulae kinematics.

We analyze the globular cluster (GC) systems in two very different galaxies, N3115 and N1399. With the papers of this series, we aim at highlighting common and different properties in the GC systems in galaxies covering a wide range of parameter space. We compare the GCs in N3115 and N1399 as derived from the analysis of u, g, and i band images taken with the VST telescope as part of the VEGAS and FDS surveys. We selected GC candidates using as reference the morpho-photometric and color properties of confirmed GCs. The surface density of GCs in N3115 reveal a morphology similar to the profile of field stars; the same is true when blue and red GCs are taken separately. The GC maps for N1399 are richer in structure and confirm the existence of an intracluster GC component. We confirm the presence of a spatial offset in the N1399 GC centroid and find that the centroid of the GCs for N3115 coincides well with the galaxy center. Both GC systems show unambiguous color bimodality in g-i and u-i. The azimuthal average of the radial density profiles in both galaxies reveals a larger spatial extent for the total GCs population with respect to the galaxy surface brightness profile. As for the specific frequency of GCs, Sn, we find it is a factor of two higher in N1399 than for N3115. By inspecting the radial behavior of the specific frequency, Sn(<r), for the total, blue, and red GCs, we find notable similarities between the trends for red GCs in the two targets. In spite of extremely different host environments, the red GCs in both cases appear closely linked to the light distribution of field stars. Blue GCs extend to larger galactocentric scales than red GCs, marking a significant difference between the two galaxies: the blue/red GCs and field stellar components of N3115 appear well thermalized with each other and the blue GCs in N1399 appear to fade into an unrelaxed intracluster GC population.

We present results of a joint \textit{Chandra}/\textit{XMM-Newton} analysis of the early-type galaxies NGC 4649 and NGC 5846 aimed at investigating differences between mass profiles derived from X-ray data and those from optical data, to probe the state of the hot ISM in these galaxies. If the hot ISM is at a given radius in hydrostatic equilibrium (HE) the X-ray data can be used to measure the total enclosed mass of the galaxy. Differences from optically-derived mass distributions therefore yield information about departures from HE in the hot halos. The X-ray mass profiles in different angular sectors of NGC 4649 are generally smooth with no significant azimuthal asymmetries within \(12\) kpc. Extrapolation of these profiles beyond this scale yields results consistent with the optical estimate. However, in the central region (\(r < 3\) kpc) the X-ray data underpredict the enclosed mass, when compared with the optical mass profiles. Consistent with previous results we estimate a non-thermal pressure component accounting for \(30\%\) of the gas pressure, likely linked to nuclear activity. In NGC 5846 the X-ray mass profiles show significant azimuthal asymmetries, especially in the NE direction. Comparison with optical mass profiles in this direction suggests significant departures from HE, consistent with bulk gas compression and decompression due to sloshing on \(\sim 15\) kpc scales; this effect disappears in the NW direction where the emission is smooth and extended. In this sector we find consistent X-ray and optical mass profiles, suggesting that the hot halo is not responding to strong non-gravitational forces.

As part of the Fornax Deep Survey with the ESO VLT Survey Telescope, we present new $g$ and $r$ bands mosaics of the SW group of the Fornax cluster. It covers an area of $3 \times 2$ square degrees around the central galaxy NGC1316. The deep photometry, the high spatial resolution of OmegaCam and the large covered area allow us to study the galaxy structure, to trace stellar halo formation and look at the galaxy environment. We map the surface brightness profile out to 33arcmin ($\sim 200$kpc $\sim15R_e$) from the galaxy centre, down to $\mu_g \sim 31$ mag arcsec$^{-2}$ and $\mu_r \sim 29$ mag arcsec$^{-2}$. This allow us to estimate the scales of the main components dominating the light distribution, which are the central spheroid, inside 5.5 arcmin ($\sim33$ kpc), and the outer stellar envelope. Data analysis suggests that we are catching in act the second phase of the mass assembly in this galaxy, since the accretion of smaller satellites is going on in both components. The outer envelope of NGC1316 still hosts the remnants of the accreted satellite galaxies that are forming the stellar halo. We discuss the possible formation scenarios for NGC1316, by comparing the observed properties (morphology, colors, gas content, kinematics and dynamics) with predictions from cosmological simulations of galaxy formation. We find that {\it i)} the central spheroid could result from at least one merging event, it could be a pre-existing early-type disk galaxy with a lower mass companion, and {\it ii)} the stellar envelope comes from the gradual accretion of small satellites.

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.

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.

Using data from the WISE all-sky survey we discovered that the non-thermal infrared (IR) emission of blazars, the largest known population of extragalactic gamma-ray sources, has peculiar spectral properties. In this work, we confirm and strengthen our previous analyses using the latest available releases of both the WISE and the Fermi source catalogs. We also show that there is a tight correlation between the mid-IR colors and the gamma-ray spectral index of Fermi blazars. We name this correlation "the infrared--gamma-ray connection". We discuss how this connection links both the emitted powers and the spectral shapes of particles accelerated in jets arising from blazars over ten decades in energy. Based on this evidence, we argue that the infrared--gamma-ray connection is stronger than the well known radio--gamma-ray connection.

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.

Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle \alpha_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population.

We report the discovery of a complex extended density enhancement in the Globular Clusters (GCs) in the central $\sim 0.5(^{\circ})^2$ ($\sim 0.06$ Mpc$^2$) of the Fornax cluster, corresponding to $\sim 50\%$ of the area within 1 core radius. This overdensity connects the GC system of NGC1399 to most of those of neighboring galaxies within $\sim 0.6^{\circ}$ ($\sim 210$ kpc) along the W-E direction. The asymmetric density structure suggests that the galaxies in the core of the Fornax cluster experienced a lively history of interactions that have left a clear imprint on the spatial distribution of GCs. The extended central dominant structure is more prominent in the distribution of blue GCs, while red GCs show density enhancements that are more centrally concentrated on the host galaxies. We propose that the relatively small-scale density structures in the red GCs are caused by galaxy-galaxy interactions, while the extensive spatial distribution of blue GCs is due to stripping of GCs from the halos of core massive galaxies by the Fornax gravitational potential. Our investigations is based on density maps of candidate GCs extracted from the multi-band VLT Survey Telescope (VST) survey of Fornax (FDS), identified in a three-dimensional color space and further selected based on their $g$-band magnitude and morphology.

[Abrigded] We have started a new deep, multi-imaging survey of the Fornax cluster, dubbed Fornax Deep Survey (FDS), at the VLT Survey Telescope. In this paper we present the deep photometry inside two square degrees around the bright galaxy NGC1399 in the core of the cluster. We found a very extended and diffuse envelope surrounding the luminous galaxy NGC1399: we map the surface brightness out to 33 arcmin (~ 192 kpc) from the galaxy center and down to about 31 mag/arcsec^2 in the g band. The deep photometry allows us to detect a faint stellar bridge in the intracluster region between NGC1399 and NGC1387. By analyzing the integrated colors of this feature, we argue that it could be due to the ongoing interaction between the two galaxies, where the outer envelope of NGC1387 on its east side is stripped away. By fitting the light profile, we found that it exists a physical break radius in the total light distribution at R=10 arcmin (~58 kpc) that sets the transition region between the bright central galaxy and the outer exponential stellar halo. We discuss the main implications of this work on the build-up of the stellar halo at the center of the Fornax cluster. By comparing with the numerical simulations of the stellar halo formation for the most massive BCGs, we find that the observed stellar halo mass fraction is consistent with a halo formed through the multiple accretion of progenitors with a stellar mass in the range 10^8 - 10^11 M_sun. This might suggest that the halo of NGC1399 has also gone through a major merging event. The absence of a significant number of luminous stellar streams and tidal tails out to 192 kpc suggests that the epoch of this strong interaction goes back to an early formation epoch. Therefore, differently from the Virgo cluster, the extended stellar halo around NGC1399 is characterised by a more diffuse and well-mixed component, including the ICL.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

Blazars are radio-loud active galactic nuclei well known for their non thermal emission spanning a wide range of frequencies. The Roma-BZCAT is, to date, the most comprehensive list of these sources. We performed the cross-match of several catalogs obtained from recent surveys at different frequencies to search for new blazars. We cross-matched the 1$^{st}$ Swift-XRT Point Source catalog with the spectroscopic sample of the 9$^{th}$ Data Release of the Sloan Digital Sky Survey. Then, we performed further cross-matches with the catalogs corresponding to the Faint Images of the Radio Sky at Twenty cm survey and to the AllWISE Data release, focusing on sources with infrared colors similar to those of confirmed $\gamma$-ray blazars included in the Second Fermi-LAT catalog. As a result, we obtained a preliminary list of objects with all the elements needed for a proper blazar classification according to the prescriptions of the Roma-BZCAT. We carefully investigated additional properties such as their morphology and the slope of their spectral energy distribution in the radio domain, the features shown in their optical spectrum, and the luminosity in the soft X rays to exclude generic active galactic nuclei and focus on authentic blazar-like sources. At the end of our screening we obtained a list of 15 objects with firmly established blazar properties.

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.

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.

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.

We report the discovery of significant localized structures in the projected two-dimensional (2D) spatial distributions of the Globular Cluster (GC) systems of the ten brightest galaxies in the Virgo Cluster. We use catalogs of GCs extracted from the HST ACS Virgo Cluster Survey (ACSVCS) imaging data, complemented, when available, by additional archival ACS data. These structures have projected sizes ranging from $\sim\!5$ arcsec to few arc-minutes ($\sim\!1$ to $\sim\!25$ kpc). Their morphologies range from localized, circular, to coherent, complex shapes resembling arcs and streams. The largest structures are preferentially aligned with the major axis of the host galaxy. A few relatively smaller structures follow the minor axis. Differences in the shape and significance of the GC structures can be noticed by investigating the spatial distribution of GCs grouped by color and luminosity. The largest coherent GC structures are located in low-density regions within the Virgo cluster. This trend is more evident in the red GC population, believed to form in mergers involving late-type galaxies. We suggest that GC over-densities may be driven by either accretion of satellite galaxies, major dissipationless mergers or wet dissipation mergers. We discuss caveats to these scenarios, and estimate the masses of the potential progenitors galaxies. These masses range in the interval $10^{8.5}\!-\!10^{9.5}$ solar masses, larger than those of the Local Group dwarf galaxies.

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.

We present an analysis of the optical spectra available in the Sloan Digital Sky survey data release nine (SDSS DR9) for the blazars listed in the ROMA-BZCAT and for the gamma-ray blazar candidates selected according to their IR colors. First, we adopt a statistical approach based on MonteCarlo simulations to find the optical counterparts of the blazarslisted in the ROMA-BZCAT catalog. Then we crossmatched the SDSS spectroscopic catalog with our selected samples of blazars and gamma-ray blazar candidates searching for those with optical spectra available to classify our blazar-like sources and, whenever possible, to confirm their redshifts. Our main objectives are determining the classification of uncertain blazars listed in the ROMA-BZCAT and discovering new gamma-ray blazars. For the ROMA-BZCAT sources we investigated a sample of 84 blazars confirming the classification for 20 of them and obtaining 18 new redshift estimates. For the gamma-ray blazars, indicated as potential counterparts of unassociated Fermi sources or with uncertain nature, we established the blazar-like nature of 8 out the 27 sources analyzed and confirmed 14 classifications.

A well known property of the gamma-ray sources detected by COS-B in the 1970s, by the Compton Gamma-ray Observatory in the 1990s and recently by the Fermi observations is the presence of radio counterparts, in particular for those associated to extragalactic objects. This observational evidence is the basis of the radio-gamma-ray connection established for the class of active galactic nuclei known as blazars. In particular, the main spectral property of the radio counterparts associated with gamma-ray blazars is that they show a flat spectrum in the GHz frequency range. Our recent analysis dedicated to search blazar-like candidates as potential counterparts for the unidentified gamma-ray sources (UGSs) allowed us to extend the radio-gamma-ray connection in the MHz regime. We also showed that below 1 GHz blazars maintain flat radio spectra. Thus on the basis of these new results, we assembled a low-frequency radio catalog of flat spectrum sources built by combining the radio observations of the Westerbork Northern Sky Survey (WENSS) and of the Westerbork in the southern hemisphere (WISH) catalog with those of the NRAO Very Large Array Sky survey (NVSS). This could be used in the future to search for new, unknown blazar-like counterparts of the gamma-ray sources. First we found NVSS counterparts of WSRT radio sources and then we selected flat spectrum radio sources according to a new spectral criterion specifically defined for radio observations performed below 1 GHz. We also described the main properties of the catalog listing 28358 radio sources and their logN-logS distributions. Finally a comparison with with the Green Bank 6-cm radio source catalog has been performed to investigate the spectral shape of the low-frequency flat spectrum radio sources at higher frequencies.

According to the second Fermi LAT Catalog (2FGL), about one third of the gamma-ray sources listed have no assigned counterparts at lower energies. Many statistical methods have been developed to find proper counterparts for these sources. We explore the sky area covered at low radio frequency by Westerbork in the Southern Hemisphere (WISH) survey to search for blazar-like associations among the unidentified gamma-ray sources listed in the 2FGL (UGSs). Searching the WISH and NRAO VLA Sky Survey (NVSS) radio surveys within the positional uncertainty regions of the 2FGL UGSs, we select as gamma-ray blazar candidates the radio sources characterized by flat radio spectra between 352 MHz and 1400 MHz. We propose new gamma-ray blazar associations for eight UGSs and we also discuss their spectral properties at low radio frequencies. We compare the radio flux density distribution of the low radio frequency gamma-ray blazar candidates with that of gamma-ray blazars associated with other methods. We find significant differences between these distributions. Finally, we discuss the results of this association method and its possible applicability to other regions of the sky and future radio surveys.

We present a catalog of radio-loud candidate gamma-ray emitting blazars with WISE mid-infrared colors similar to the colors of confirmed gamma-ray blazars. The catalog is assembled from WISE sources detected in all four WISE filters, with colors compatible with the three-dimensional locus of the WISE gamma-ray emitting blazars, and which can be spatially cross-matched with radio sources from either one of the three radio surveys: NVSS, FIRST and/or SUMSS. Our initial WISE selection uses a slightly modified version of previously successful algorithms. We then select only the radio-loud sources using a measure of the radio-to-IR flux, the q22 parameter, which is analogous to the q24 parameter known in the literature but which instead uses the WISE band-four flux at 22 micron. Our final catalog contains 7855 sources classified as BL Lacs, FSRQs or mixed candidate blazars; 1295 of these sources can be spatially re-associated with confirmed blazars. We describe the properties of the final catalog of WISE blazar-like radio-loud sources and consider possible contaminants. Finally, we discuss why this large catalog of candidate gamma-ray emitting blazars represents a new and useful resource to address the problem of finding low energy counterparts to currently unidentified high-energy sources.

Iris is an extensible application that provides astronomers with a user-friendly interface capable of ingesting broad-band data from many different sources in order to build, explore, and model spectral energy distributions (SEDs). Iris takes advantage of the standards defined by the International Virtual Observatory Alliance, but hides the technicalities of such standards by implementing different layers of abstraction on top of them. Such intermediate layers provide hooks that users and developers can exploit in order to extend the capabilities provided by Iris. For instance, custom Python models can be combined in arbitrary ways with the Iris built-in models or with other custom functions. As such, Iris offers a platform for the development and integration of SED data, services, and applications, either from the user's system or from the web. In this paper we describe the built-in features provided by Iris for building and analyzing SEDs. We also explore in some detail the Iris framework and software development kit, showing how astronomers and software developers can plug their code into an integrated SED analysis environment.

Astronomy is undergoing through a methodological revolution triggered by an unprecedented wealth of complex and accurate data. The new panchromatic, synoptic sky surveys require advanced tools for discovering patterns and trends hidden behind data which are both complex and of high dimensionality. We present DAMEWARE (DAta Mining & Exploration Web Application REsource): a general purpose, web-based, distributed data mining environment developed for the exploration of large datasets, and finely tuned for astronomical applications. By means of graphical user interfaces, it allows the user to perform classification, regression or clustering tasks with machine learning methods. Salient features of DAMEWARE include its capability to work on large datasets with minimal human intervention, and to deal with a wide variety of real problems such as the classification of globular clusters in the galaxy NGC1399, the evaluation of photometric redshifts and, finally, the identification of candidate Active Galactic Nuclei in multiband photometric surveys. In all these applications, DAMEWARE allowed to achieve better results than those attained with more traditional methods. With the aim of providing potential users with all needed information, in this paper we briefly describe the technological background of DAMEWARE, give a short introduction to some relevant aspects of data mining, followed by a summary of some science cases and, finally, we provide a detailed description of a template use case.

Nearly one-third of the gamma-ray sources detected by Fermi are still unidentified, despite significant recent progress in this effort. On the other hand, all the gamma-ray extragalactic sources associated in the second Fermi-LAT catalog have a radio counterpart. Motivated by this observational evidence we investigate all the radio sources of the major radio surveys that lie within the positional uncertainty region of the unidentified gamma-ray sources (UGSs) at 95% level of confidence. First we search for their infrared counterparts in the all-sky survey performed by the Wide-field Infrared Survey Explorer (WISE) and then we analyze their IR colors in comparison with those of the known gamma-ray blazars. We propose a new approach, based on a 2-dimensional kernel density estimation (KDE) technique in the single [3.4]-[4.6]-[12] micron WISE color-color plot, replacing the constraint imposed in our previous investigations on the detection at 22 micron of each potential IR counterpart of the UGSs with associated radio emission. The main goal of this analysis is to find distant gamma-ray blazar candidates that, being too faint at 22 micron, are not detected by WISE and thus are not selected by our purely IR based methods. We find fifty-five UGS's likely correspond to radio sources with blazar-like IR signatures. Additional eleven UGSs having, blazar-like IR colors, have been found within the sample of sources found with deep recent ATCA observations.

A significant fraction ($\sim 30$ %) of the gamma-ray sources listed in the second $\textit{Fermi}$ LAT (2FGL) catalog is still of unknown origin, being not yet associated with counterparts at lower energies. Using the available information at lower energies and optical spectroscopy on the selected counterparts of these gamma-ray objects we can pinpoint their exact nature. Here we present a pilot project pointing to assess the effectiveness of the several classification methods developed to select gamma-ray blazar candidates. To this end, we report optical spectroscopic observations of a sample of 5 gamma-ray blazar candidates selected on the basis of their infrared WISE colors or of their low-frequency radio properties. Blazars come in two main classes: BL Lacs and FSRQs, showing similar optical spectra except for the stronger emission lines of the latter. For three of our sources the almost featureless optical spectra obtained confirm their BL Lac nature, while for the source WISEJ022051.24+250927.6 we observe emission lines with equivalent width $EW\sim 31$ $\AA$, identifying it as a FSRQ with $z = 0.48$. The source WISEJ064459.38+603131.7, although not featuring a clear radio counterpart, shows a blazar-like spectrum with weak emission lines with $EW \sim 7$ $\AA$, yielding a redshift estimate of $z=0.36$. In addition we report optical spectroscopic observations of 4 WISE sources associated with known gamma-ray blazars without a firm classification or redshift estimate. For all of these latter sources we confirm a BL Lac classification, with a tentative redshift estimate for the source WISEJ100800.81+062121.2 of $z = 0.65$.

A significant fraction ($\sim 30$ %) of the high-energy gamma-ray sources listed in the second $\textit{Fermi}$ LAT (2FGL) catalog are still of unknown origin, being not yet associated with counterparts at lower energies. In order to investigate the nature of these enigmatic sources, we present here an extensive search of X-ray sources lying in the positional uncertainty region of a selected sample of these Unidentified Gamma-ray Sources (UGSs) that makes use of all available observations performed by the $\textit{Swift}$ X-ray Telescope before March 31, 2013, available for 205 UGSs. To detect the fainter sources, we merged all the observations covering the $\textit{Fermi}$ LAT positional uncertainty region at 95 % level of confidence of each UGSs. This yields a catalog of 357 X-ray sources, finding {candidate} X-ray counterparts for $\sim 70$ % of the selected sample. In particular, 25 % of the UGSs feature a single X-ray source within their positional uncertainty region while 45 % have multiple X-ray sources. For each X-ray source we also looked in the corresponding $\textit{Swift}$ UVOT merged images for optical and ultraviolet counterparts, also performing source photometry. We found ultraviolet-optical correspondences for $\sim 70$ % of the X-ray sources. We searched several major radio, infrared, optical and ultraviolet surveys for possible counterparts within the positional error of the sources in the X-ray catalog to obtain additional information on their nature. Applying the kernel density estimator technique to infrared colors of WISE counterparts of our X-ray sources we select 6 $\gamma$-ray blazar candidates. In addition, comparing our results with previous analyses, we select 11 additional $\gamma$-ray blazar candidates.

We report significant inhomogeneities in the projected two-dimensional (2D) spatial distributions of Low-Mass X-ray Binaries (LMXBs) and Globular Clusters (GCs) of the intermediate mass elliptical galaxy NGC4278. In the inner region of NGC4278, a significant arc-like excess of LMXBs extending south of the center at ~50" in the western side of the galaxy can be associated to a similar over-density of the spatial distribution of red GCs from~Brassington et al. (2009). Using a recent catalog of GCs produced by Usher et al.(2013) and covering the whole field of the NGC4278 galaxy, we have discovered two other significant density structures outside the D25 isophote to the W and E of the center of NGC4278, associated to an over-density and an under-density respectively. We discuss the nature of these structures in the context of the similar spatial inhomogeneities discovered in the LMXBs and GCs populations of NGC4649 and NGC4261, respectively. These features suggest streamers from disrupted and accreted dwarf companions.

We report significant anisotropies in the projected two-dimensional (2D) spatial distributions of Globular Clusters (GCs) of the giant Virgo elliptical galaxy NGC4649 (M60). Similar features are found in the 2D distribution of low-mass X-ray binaries (LMXBs), both associated with GCs and in the stellar field. Deviations from azimuthal symmetry suggest an arc-like excess of GCs extending north at 4-15 kpc galactocentric radii in the eastern side of major axis of NGC4649. This feature is more prominent for red GCs, but still persists in the 2D distribution of blue GCs. High and low luminosity GCs also show some segregation along this arc, with high-luminosity GCs preferentially located in the southern end and low-luminosity GCs in the northern section of the arc. GC-LMXBs follow the anisotropy of red-GCs, where most of them reside; however, a significant overdensity of (high-luminosity) field LMXBs is present to the south of the GC arc. These results suggest that NGC4649 has experienced mergers and/or multiple accretions of less massive satellite galaxies during its evolution, of which the GCs in the arc may be the fossil remnant. We speculate that the observed anisotropy in the field LMXB spatial distribution indicates that these X-ray binaries may be the remnants of a star formation event connected with the merger, or maybe be ejected from the parent red GCs, if the bulk motion of these clusters is significantly affected by dynamical friction. We also detect a luminosity enhancement in the X-ray source population of the companion spiral galaxy NGC4647. We suggest that these may be younger high mass X-ray binaries formed as a result of the tidal interaction of this galaxy with NGC4649.

We investigate the radial distribution of the low-mass X-ray binary (LMXB) population in the elliptical galaxy NGC 4649, using Chandra and Hubble data to separate the field and globular cluster (GC) populations. GCs with LMXBs have the same radial distribution as the parent red and blue GCs. The radial profile of field LMXBs follows the V-band profile within the D25 of NGC 4649. Using the spatial information provided by our data, we find that the global galaxy-wide relations between cumulative number and luminosity of LMXBs and the integrated stellar mass hold on local scales within D25. An excess of field LMXBs with respect to the V-light is observed in the galaxy's outskirt, which may be partially due to unidentified GC sources or to a rejuvenated field LMXB population caused by past merging interactions.

In this paper we discuss an application of machine learning based methods to the identification of candidate AGN from optical survey data and to the automatic classification of AGNs in broad classes. We applied four different machine learning algorithms, namely the Multi Layer Perceptron (MLP), trained respectively with the Conjugate Gradient, Scaled Conjugate Gradient and Quasi Newton learning rules, and the Support Vector Machines (SVM), to tackle the problem of the classification of emission line galaxies in different classes, mainly AGNs vs non-AGNs, obtained using optical photometry in place of the diagnostics based on line intensity ratios which are classically used in the literature. Using the same photometric features we discuss also the behavior of the classifiers on finer AGN classification tasks, namely Seyfert I vs Seyfert II and Seyfert vs LINER. Furthermore we describe the algorithms employed, the samples of spectroscopically classified galaxies used to train the algorithms, the procedure followed to select the photometric parameters and the performances of our methods in terms of multiple statistical indicators. The results of the experiments show that the application of self adaptive data mining algorithms trained on spectroscopic data sets and applied to carefully chosen photometric parameters represents a viable alternative to the classical methods that employ time-consuming spectroscopic observations.

Blazars are a highly-variable, radio-loud subclass of active galactic nuclei (AGN). In order to better understand such objects we must be able to easily identify candidate blazars from the growing population of unidentified sources. Working towards this goal we attempt to identify new gamma-ray blazar candidates from a sample of 102 previously unidentified sources. These sources are selected from Astronomer's Telegrams and the literature on the basis of non-periodic variability and multi-wavelength behavior. We then attempt to associate these objects to an IR counterpart in the WISE all-sky survey. We are able to identify sixteen candidate sources whose IR colors are consistent with those of the blazar population. Of those sixteen, thirteen sources have IR colors indicative of being gamma-ray emitting blazar candidates. These sources all possess archival multi-wavelength observations that support their blazar-like nature.

Blazars are the most extreme class of active galactic nuclei (AGNs). Despite a previous investigation at 102 MHz for a small sample of BL Lacs and our recent analysis of blazars detected in the Westerbork Northern Sky Survey (WENSS), a systematic study of the blazar spectral properties at frequencies below 100 MHz has been never carried out. In this paper, we present the first analysis of the radio spectral behavior of blazars based on the recent Very Large Array Low-Frequency Sky Survey (VLSS) at 74 MHz. We search for blazar counterparts in the VLSS catalog confirming that they are detected at 74 MHz. We then show that blazars present radio flat spectra {(i.e., radio spectral indices ~0.5)} when evaluated also about an order of magnitude in frequency lower than previous analyses. Finally, we discuss the implications of our findings in the context of the blazars - radio galaxies connection since the low frequency radio data provide a new diagnostic tool to verify the expectations of the unification scenario for radio-loud active galaxies.

About one third of the gamma-ray sources listed in the second Fermi LAT catalog (2FGL) have no firmly established counterpart at lower energies so being classified as unidentified gamma-ray sources (UGSs). Here we propose a new approach to find candidate counterparts for the UGSs based on the 325 MHz radio survey performed with Westerbork Synthesis Radio Telescope (WSRT) in the northern hemisphere. First we investigate the low-frequency radio properties of blazars, the largest known population of gamma-ray sources; then we search for sources with similar radio properties combining the information derived from the Westerbork Northern Sky Survey (WENSS) with those of the NRAO VLA Sky survey (NVSS). We present a list of candidate counterparts for 32 UGSs with at least one counterpart in the WENSS. We also performed an extensive research in literature to look for infrared and optical counterparts of the gamma-ray blazar candidates selected with the low-frequency radio observations to confirm their nature. On the basis of our multifrequency research we identify 23 new gamma-ray blazar candidates out of 32 UGSs investigated. Comparison with previous results on the UGSs are also presented. Finally, we speculate on the advantages on the use of the low-frequency radio observations to associate UGSs and to search for gamma-ray pulsar candidates.

We searched for gamma-ray blazar candidates among the 382 unidentified hard X-ray sources of the 3rd Palermo BAT Catalog (3PBC) obtained from the analysis of 66 months of SWIFT-BAT survey data and listing 1586 sources. We adopted a recently developed association method based on the peculiar infrared colors which characterize the gamma-ray blazars included in the second catalog of active galactic nuclei detected by the Fermi Large Area Telescope (2LAC). We used this method exploiting the data of the all-sky survey performed by the Wide-Field Infrared Survey Explorer (WISE) to establish correspondences between unidentified 3PBC sources and WISE gamma-ray blazar candidates located within the BAT positional uncertainty region at 99% confidence level. We obtained a preliminary list of candidates for which we analysed all the available data in the SWIFT archive to complement the information in the literature and in the radio, infrared and optical catalogs with the information on their optical-UV and soft X-ray emission. Requiring the presence of radio and soft X-ray counterparts consistent with the infrared positions of the selected WISE sources, as well as a blazar-like radio morphology, we finally obtained a list of 24 gamma-ray blazar candidates.

BL Lac objects are the most numerous class of extragalactic TeV-detected sources. One of the biggest difficulties in investigating their TeV emission resides in their limited number, since only 47 BL Lacs are known as TeV emitters. In this paper, we propose new criteria to select TeV BL Lac candidates based on the infrared (IR) and X-ray observations. We apply our selection criteria to the BL Lac objects listed in the ROMA-BZCAT catalog so identifying 41 potential TeV emitters. We then consider a search over a more extended sample combining the ROSAT bright source catalog and the WISE all-sky survey revealing 54 additional candidates for TeV observations. Our investigation also led to a tentative classification of 16 unidentified X-ray sources as BL Lac candidates. This analysis provides new interesting BL Lac targets for future observations with ground based Cherenkov telescopes.

We present a new method for the determination of the two-dimensional (2D) projected spatial distribution of globular clusters (GCs) in external galaxies. This method is based on the K-Nearest Neighbor density estimator of Dressler (1980), complemented by MonteCarlo simulations to establish the statistical significance of the results. We apply this method to NGC4261, a "test galaxy" where significant 2D anisotropy in the GC distribution has been reported. We confirm that the 2D distribution of GC is not azimuthally isotropic. Moreover, we demonstrate that the 2D distribution departures from the average GC radial distribution results in highly significant spiral-like or broken shell features. Overall, the same perturbations are found in "red" and "blue" GCs, but with some differences. In particular, we observe a central feature, roughly aligned with the minor axis of NGC4261, composed of red and most luminous GCs. Blue and fainter GCs are more frequent at large radial distances and follow the spiral-like features of the overall density structure. These results suggest a complex merging history for NGC4261.

MLPQNA stands for Multi Layer Perceptron with Quasi Newton Algorithm and it is a machine learning method which can be used to cope with regression and classification problems on complex and massive data sets. In this paper we give the formal description of the method and present the results of its application to the evaluation of photometric redshifts for quasars. The data set used for the experiment was obtained by merging four different surveys (SDSS, GALEX, UKIDSS and WISE), thus covering a wide range of wavelengths from the UV to the mid-infrared. The method is able i) to achieve a very high accuracy; ii) to drastically reduce the number of outliers and catastrophic objects; iii) to discriminate among parameters (or features) on the basis of their significance, so that the number of features used for training and analysis can be optimized in order to reduce both the computational demands and the effects of degeneracy. The best experiment, which makes use of a selected combination of parameters drawn from the four surveys, leads, in terms of DeltaZnorm (i.e. (zspec-zphot)/(1+zspec)), to an average of DeltaZnorm = 0.004, a standard deviation sigma = 0.069 and a Median Absolute Deviation MAD = 0.02 over the whole redshift range (i.e. zspec <= 3.6), defined by the 4-survey cross-matched spectroscopic sample. The fraction of catastrophic outliers, i.e. of objects with photo-z deviating more than 2sigma from the spectroscopic value is < 3%, leading to a sigma = 0.035 after their removal, over the same redshift range. The method is made available to the community through the DAMEWARE web application.

A significant fraction (~30%) of the high-energy gamma-ray sources listed in the second Fermi LAT catalog (2FGL) are still of unknown origin, being not yet associated with counterparts at low energies. We recently developed a new association method to identify if there is a gamma-ray blazar candidate within the positional uncertainty region of a generic 2FGL source. This method is entirely based on the discovery that blazars have distinct infrared colors with respect to other extragalactic sources found thanks, to the Wide-field Infrared Survey Explorer (WISE) all-sky observations. Several improvements have been also performed to increase the efficiency of our method in recognizing gamma-ray blazar candidates. In this paper we applied our method to two different samples, the first constituted by the unidentified gamma-ray sources (UGSs) while the second by the active galaxies of uncertain type (AGUs), both listed in the 2FGL. We present a catalog of IR counterparts for ~20% of the UGSs investigated. Then, we also compare our results on the associated sources with those present in literature. In addition, we illustrate the extensive archival research carried out to identify the radio, infrared, optical and X-ray counterparts of the WISE selected, gamma-ray blazar candidates. Finally, we discuss the future developments of our method based on ground-based follow-up observations.

One of the main scientific objectives of the recent Fermi mission is investigating the origin of the unidentified gamma-ray sources (UGSs). Despite the large improvements of Fermi in the gamma-ray source localization with respect to the past gamma-ray missions, about 1/3 of the gamma-ray objects detected do not have yet an assigned counterpart a low energies. We recently developed a new association method to identify if there is a gamma-ray blazar candidate within the positional uncertainty region of a generic gamma-ray source. This method is entirely based on the discovery that blazars, the largest known class of gamma-ray sources, can be recognized and separated from other extragalactic sources on the basis of their infrared colors. Here we summarize the results obtained by applying our association procedure to the unidentified gamma-ray sources (UGSs) and to the active galaxies of uncertain type (AGUs) listed in the Second Fermi Large Area Telescope (LAT) catalog.

We present a new method for identifying blazar candidates by examining the locus, i.e. the region occupied by the Fermi gamma-ray blazars in the three-dimensional color space defined by the WISE infrared colors. This method is a refinement of our previous approach that made use of the two-dimensional projection of the distribution of WISE gamma-ray emitting blazars (the Strip) in the three WISE color-color planes (Massaro et al. 2012a). In this paper, we define the three-dimensional locus by means of a Principal Component (PCs) analysis of the colors distribution of a large sample of blazars composed by all the ROMA-BZCAT sources with counterparts in the WISE All-Sky Catalog and associated to gamma-ray source in the second Fermi LAT catalog (the WISE Fermi Blazars sample, WFB). Our new procedure yields a total completeness of c~81% and total efficiency of e~97%. We also obtain local estimates of the efficiency and completeness as functions of the WISE colors and galactic coordinates of the candidate blazars. The catalog of all WISE candidate blazars associated to the WFB sample is also presented, complemented by archival multi-frequency information for the alternative associations. Finally, we apply the new association procedure to all gamma-ray blazars in the 2FGL and provide a catalog containing all the gamma-ray candidates blazars selected according to our procedure.

We present a new method for identifying blazar candidates by examining the locus, i.e. the region occupied by the Fermi gamma-ray blazars in the three-dimensional color space defined by the WISE infrared colors. This method is a refinement of our previous approach that made used of the two-dimensional projection of the the distribution of WISE gamma-ray emitting blazars (the Strip) in the three WISE color-color planes. In this paper, we define the three-dimensional locus by means of a Principal Component analysis of the the colors distribution of a large sample of blazars composed by all the ROMA-BZCAT sources with counterparts in the WISE All-Sky Catalog and associated to a gamma-ray source in the second Fermi LAT catalog (2FGL) (the WISE Fermi Blazars sample, WFB). Our new procedure, as reported in [D'Abrusco et al. 2013], yields a total completeness of c=81% and total efficiency of e=97%.

The existence of correlations between nuclear properties of galaxies, such as the mass of their central black holes, and larger scale features, like the bulge mass and luminosity, represent a fundamental constraint on galaxy evolution. Although the actual reasons for these relations have not yet been identified, it is widely believed that they could stem from a connection between the processes that lead to black hole growth and stellar mass assembly. The problem of understanding how the processes of nuclear activity and star formation can affect each other became known to the literature as the Starburst-AGN connection. Despite years of investigation, the physical mechanisms which lie at the basis of this relation are known only in part. In this work, we analyze the problem of star formation and nuclear activity in a large sample of galaxies. We study the relations between the properties of the nuclear environments and of their host galaxies. We find that the mass of the stellar component within the galaxies of our sample is a critical parameter, that we have to consider in an evolutionary sequence, which provides further insight in the connection between AGN and star formation processes.

The U.S. Virtual Astronomical Observatory (VAO) is a product-driven organization that provides new scientific research capabilities to the astronomical community. Software development for the VAO follows a lightweight framework that guides development of science applications and infrastructure. Challenges to be overcome include distributed development teams, part-time efforts, and highly constrained schedules. We describe the process we followed to conquer these challenges while developing Iris, the VAO application for analysis of 1-D astronomical spectral energy distributions (SEDs). Iris was successfully built and released in less than a year with a team distributed across four institutions. The project followed existing International Virtual Observatory Alliance inter-operability standards for spectral data and contributed a SED library as a by-product of the project. We emphasize lessons learned that will be folded into future development efforts. In our experience, a well-defined process that provides guidelines to ensure the project is cohesive and stays on track is key to success. Internal product deliveries with a planned test and feedback loop are critical. Release candidates are measured against use cases established early in the process, and provide the opportunity to assess priorities and make course corrections during development. Also key is the participation of a stakeholder such as a lead scientist who manages the technical questions, advises on priorities, and is actively involved as a lead tester. Finally, frequent scheduled communications (for example a bi-weekly tele-conference) assure issues are resolved quickly and the team is working toward a common vision

In this paper we present the Clustering-Labels-Score Patterns Spotter (CLaSPS), a new methodology for the determination of correlations among astronomical observables in complex datasets, based on the application of distinct unsupervised clustering techniques. The novelty in CLaSPS is the criterion used for the selection of the optimal clusterings, based on a quantitative measure of the degree of correlation between the cluster memberships and the distribution of a set of observables, the labels, not employed for the clustering. In this paper we discuss the applications of CLaSPS to two simple astronomical datasets, both composed of extragalactic sources with photometric observations at different wavelengths from large area surveys. The first dataset, CSC+, is composed of optical quasars spectroscopically selected in the SDSS data, observed in the X-rays by Chandra and with multi-wavelength observations in the near-infrared, optical and ultraviolet spectral intervals. One of the results of the application of CLaSPS to the CSC+ is the re-identification of a well-known correlation between the alphaOX parameter and the near ultraviolet color, in a subset of CSC+ sources with relatively small values of the near-ultraviolet colors. The other dataset consists of a sample of blazars for which photometric observations in the optical, mid and near infrared are available, complemented for a subset of the sources, by Fermi gamma-ray data. The main results of the application of CLaSPS to such datasets have been the discovery of a strong correlation between the multi-wavelength color distribution of blazars and their optical spectral classification in BL Lacs and Flat Spectrum Radio Quasars and a peculiar pattern followed by blazars in the WISE mid-infrared colors space. This pattern and its physical interpretation have been discussed in details in other papers by one of the authors.

We present Iris, the VAO (Virtual Astronomical Observatory) application for analyzing SEDs (spectral energy distributions). Iris is the result of one of the major science initiatives of the VAO, and the first version was released in September 2011. Iris combines key features of several existing software applications to streamline and enhance SED analysis. With Iris, users may read and display SEDs, select data ranges for analysis, fit models to SEDs, and calculate confidence limits on best-fit parameters. SED data may be uploaded into the application from IVOA-compliant VOTable and FITS format files, or retrieved directly from NED. Data written in unsupported formats may be converted using SedImporter, a new application provided with Iris. The components of Iris have been contributed by members of the VAO. Specview, contributed by STScI, provides a GUI for reading, editing, and displaying SEDs, as well as defining models and parameter values. Sherpa, contributed by the Chandra project at SAO, provides a library of models, fit statistics, and optimization methods; the underlying I/O library, SEDLib, is a VAO product written by SAO to current IVOA (International Virtual Observatory Alliance) data model standards. NED is a service provided by IPAC for easy location of data for a given extragalactic source, including SEDs. SedImporter is a new tool for converting non-standard SED data files into a format supported by Iris. We demonstrate the use of SedImporter to retrieve SEDs from a variety of sources--from the NED SED service, from the user's own data, and from other VO applications using SAMP (Simple Application Messaging Protocol). We also demonstrate the use of Iris to read, display, select ranges from, and fit models to SEDs. Finally, we discuss the architecture of Iris, and the use of IVOA standards so that Specview, Sherpa, SEDLib and SedImporter work together seamlessly.

The identification of low-energy counterparts for gamma-ray sources is one of the biggest challenge in modern gamma-ray astronomy. Recently, we developed and successfully applied a new association method to recognize gamma-ray blazar candidates that could be possible counterparts for the unidentified gamma-ray sources above 100 MeV in the second Fermi Large Area Telescope (LAT) catalog (2FGL). This method is based on the Infrared (IR) colors of the recent Wide-Field Infrared Survey Explorer (WISE) all-sky survey. In this letter we applied our new association method to the case of unidentified INTEGRAL sources (UISs) listed in the fourth soft gamma-ray source catalog (4IC). Only 86 UISs out of the 113 can be analyzed, due to the sky coverage of the WISE Preliminary data release. Among these 86 UISs, we found that 18 appear to have a gamma-ray blazar candidate within their positional error region. Finally, we analyzed the SWIFT archival data available for 10 out these 18 gamma-ray blazar candidates, and we found that 7 out of 10 are clearly detected in soft X-rays and/or in the optical-ultraviolet band. We cannot confirm the associations between the UISs and the selected gamma-ray blazar candidates due to the discrepancies between the INTEGRAL and the soft X-ray spectra. However, the discovery of the soft X-ray counterparts for the selected gamma-ray blazar candidates adds an important clue to help understand their origin and to confirm their blazar nature.

One of the main scientific objectives of the ongoing Fermi mission is unveiling the nature of the unidentified gamma-ray sources (UGSs). Despite the large improvements of Fermi in the localization of gamma-ray sources with respect to the past gamma-ray missions, about one third of the Fermi-detected objects are still not associated to low energy counterparts. Recently, using the Wide-Field Infrared Survey Explorer (WISE) survey, we discovered that blazars, the rarest class of Active Galactic Nuclei and the largest population of gamma-ray sources, can be recognized and separated from other extragalactic sources on the basis of their infrared (IR) colors. Based on this result, we designed an association method for the gamma-ray sources to reognize if there is a blazar candidate within the positional uncertainty region of a generic gamma-ray source. With this new IR diagnostic tool, we searched for gamma-ray blazar candidates associated to the UGS sample of the second Fermi gamma-ray catalog (2FGL). We found that our method associates at least one gamma-ray blazar candidate as a counterpart each of 156 out of 313 UGSs analyzed. These new low-energy candidates have the same IR properties as the blazars associated to gamma-ray sources in the 2FGL catalog.

Despite the large number of discoveries made recently by Fermi, the origins of the so called unidentified gamma-ray sources remain unknown. The large number of these sources suggests that among them there could be a population that significantly contributes to the isotropic gamma-ray background and is therefore crucial to understand their nature. The first step toward a complete comprehension of the unidentified gamma-ray source population is to identify those that can be associated with blazars, the most numerous class of extragalactic sources in the gamma-ray sky. Recently, we discovered that blazars can be recognized and separated from other extragalactic sources using the infrared (IR) WISE satellite colors. The blazar population delineates a remarkable and distinctive region of the IR color-color space, the WISE blazar strip. In particular, the subregion delineated by the gamma-ray emitting blazars is even narrower and we named it as the WISE Gamma-ray Strip (WGS). In this paper we parametrize the WGS on the basis of a single parameter s that we then use to determine if gamma-ray Active Galactic Nuclei of the uncertain type (AGUs) detected by Fermi are consistent with the WGS and so can be considered blazar candidates. We find that 54 AGUs out of a set 60 analyzed have IR colors consistent with the WGS; only 6 AGUs are outliers. This result implies that a very high percentage (i.e., in this sample about 90%) of the AGUs detected by Fermi are indeed blazar candidates.

Blazars constitute the most enigmatic class of extragalactic gamma-ray sources, and their observational features have been ascribed to a relativistic jet closely aligned to the line of sight. They are generally divided in two main classes: the BL Lac objects (BL Lacs) and the Flat Spectrum Radio Quasars (FSRQs). In the case of BL Lacs the double bumped spectral energy distribution (SED) is generally described by the Synchrotron Self Compton (SSC) emission, while for the FSRQs it is interpreted as due to External Compton (EC) emission. Recently, we showed that in the [3.4]-[4.6]-[12] micron color- color diagram the blazar population covers a distinct region (i.e., the WISE blazar Strip, WBS), clearly separated from the other extragalactic sources that are dominated by thermal emission. In this paper we investigate the relation between the infrared and gamma-ray emission for a subset of confirmed blazars from the literature, associated with Fermi sources, for which WISE archival observations are available. This sample is a proper subset of the sample of sources used previously, and the availability of Fermi data is critical to constrain the models on the emission mechanisms for the blazars. We found that the selected blazars also lie on the WISE blazar Strip covering a narrower region of the infrared color-color planes than the overall blazars population. We then found an evident correlation between the IR and gamma-ray spectral indices expected in the SSC and EC frameworks. Finally, we determined the ratio between their gamma-ray and infrared fluxes, a surrogate of the ratio of powers between the inverse Compton and the synchrotron SED components, and used such parameter to test different emitting scenarios blazars.

Blazars constitute the most interesting and enigmatic class of extragalactic gamma-ray sources dominated by non-thermal emission. In this Letter, we show how the WISE infrared data make possible to identify a distinct region of the [3.4]-[4.6]-[12] micron color-color diagram where the sources dominated by the the thermal radiation are separated from those dominated by non-thermal emission, in particular the blazar population. This infrared non-thermal region delineated as the WISE Blazar Strip (WBS), it is a powerful new diagnostic tool when the full WISE survey data is released. The WBS can be used to extract new blazar candidates, to identify those of uncertain type and also to search for the counterparts of unidentified gamma-ray sources. We show one example of the value of the use of the WBS identifying the TeV source VER J 0648+152, recently discovered by VERITAS.

Astronomical data are gathered through a very large number of heterogeneous techniques and stored in very diversified and often incompatible data repositories. Moreover in the e-science environment, it is needed to integrate services across distributed, heterogeneous, dynamic "virtual organizations" formed by different resources within a single enterprise and/or external resource sharing and service provider relationships. The DAME/VONeural project, run jointly by the University Federico II, INAF (National Institute of Astrophysics) Astronomical Observatories of Napoli and the California Institute of Technology, aims at creating a single, sustainable, distributed e-infrastructure for data mining and exploration in massive data sets, to be offered to the astronomical (but not only) community as a web application. The framework makes use of distributed computing environments (e.g. S.Co.P.E.) and matches the international IVOA standards and requirements. The integration process is technically challenging due to the need of achieving a specific quality of service when running on top of different native platforms. In these terms, the result of the DAME/VO-Neural project effort will be a service-oriented architecture, obtained by using appropriate standards and incorporating Grid paradigms and restful Web services frameworks where needed, that will have as main target the integration of interdisciplinary distributed systems within and across organizational domains.

Nowadays, many scientific areas share the same broad requirements of being able to deal with massive and distributed datasets while, when possible, being integrated with services and applications. In order to solve the growing gap between the incremental generation of data and our understanding of it, it is required to know how to access, retrieve, analyze, mine and integrate data from disparate sources. One of the fundamental aspects of any new generation of data mining software tool or package which really wants to become a service for the community is the possibility to use it within complex workflows which each user can fine tune in order to match the specific demands of his scientific goal. These workflows need often to access different resources (data, providers, computing facilities and packages) and require a strict interoperability on (at least) the client side. The project DAME (DAta Mining & Exploration) arises from these requirements by providing a distributed WEB-based data mining infrastructure specialized on Massive Data Sets exploration with Soft Computing methods. Originally designed to deal with astrophysical use cases, where first scientific application examples have demonstrated its effectiveness, the DAME Suite results as a multi-disciplinary platform-independent tool perfectly compliant with modern KDD (Knowledge Discovery in Databases) requirements and Information & Communication Technology trends.

With the availability of the huge amounts of data produced by current and future large multi-band photometric surveys, photometric redshifts have become a crucial tool for extragalactic astronomy and cosmology. In this paper we present a novel method, called Weak Gated Experts (WGE), which allows to derive photometric redshifts through a combination of data mining techniques. \noindent The WGE, like many other machine learning techniques, is based on the exploitation of a spectroscopic knowledge base composed by sources for which a spectroscopic value of the redshift is available. This method achieves a variance \sigma^2(\Delta z)=2.3x10^{-4} (\sigma^2(\Delta z) =0.08), where \Delta z = z_{phot} - z_{spec}) for the reconstruction of the photometric redshifts for the optical galaxies from the SDSS and for the optical quasars respectively, while the Root Mean Square (RMS) of the \Delta z variable distributions for the two experiments is respectively equal to 0.021 and 0.35. The WGE provides also a mechanism for the estimation of the accuracy of each photometric redshift. We also present and discuss the catalogs obtained for the optical SDSS galaxies, for the optical candidate quasars extracted from the DR7 SDSS photometric dataset {The sample of SDSS sources on which the accuracy of the reconstruction has been assessed is composed of bright sources, for a subset of which spectroscopic redshifts have been measured.}, and for optical SDSS candidate quasars observed by GALEX in the UV range. The WGE method exploits the new technological paradigm provided by the Virtual Observatory and the emerging field of Astroinformatics.

This document describes the scientific requirements for the SED builder and analysis tool that will be designed and built by the US Virtual Astronomical Observatory (VAO). VAO is the VO effort based in the US, whose primary emphasis is to provide new scientific research capabilities to the astronomy community. The near-term goal of the VAO is to put useful and efficient tools in the hands of research astronomers as soon as possible. The VAO has identified eight major research initiatives that are of high priority, including the science obtained with the creation and analysis of the spectral energy distributions of astronomical sources. This document contains the high-level scientific requirements that will drive the design and implementation of the VAO SED tool.

Nowadays, many scientific areas share the same need of being able to deal with massive and distributed datasets and to perform on them complex knowledge extraction tasks. This simple consideration is behind the international efforts to build virtual organizations such as, for instance, the Virtual Observatory (VObs). DAME (DAta Mining & Exploration) is an innovative, general purpose, Web-based, VObs compliant, distributed data mining infrastructure specialized in Massive Data Sets exploration with machine learning methods. Initially fine tuned to deal with astronomical data only, DAME has evolved in a general purpose platform which has found applications also in other domains of human endeavor. We present the products and a short outline of a science case, together with a detailed description of main features available in the beta release of the web application now released.

Aim of this communication is to describe the first results of a work-in-progress regarding the chemical properties of gas and stars in the circumnuclear regions of nearby galaxies. Different techniques have been employed to estimate the abundances of chemical elements in the gaseous and stellar components of nuclear surroundings in different classes of galaxies according to the level of activity of the nucleus (normal or passive, star forming galaxies and AGNs).

In order to characterize the amount of recent or ongoing stellar formation in the circumnuclear region of active galaxies on a statistically sound basis, we have studied the stellar component of the nuclear spectra in three different samples of galaxies, namely Seyfert 2 galaxies (hereafter S2G), star-forming galaxies (SFG) and passive normal galaxies (NG), i.e., no emission lines observed, using Sloan Digital Sky Survey data (SDSS) (Adelman-McCarthy, 2008). The stellar component of the observed spectra has been extracted using STARLIGHT (Cid Fernandes et al., 2004), which fits an observed spectrum with a model (template) spectrum obtained by combining a library of pre-defined simple stellar populations spectra, with distinct ages and metallicities. The resulting template spectra for the different samples of galaxies have been compared to determine the features of the stellar emission component and to evaluate the presence and intensity of the star formation in the nuclear regions of different families of galaxies. From a first qualitative analysis it results that the shape of the Spectral Energy Distribution (SED) of S2G and NG is very similar, while that of SFG is characterized by a strong blue excess. The presence of the 4000 A break in the spectra of S2G and NG together with the lack of a strong blue continuum clearly indicate the absence of ongoing star formation in the circumnuclear regions of S2G and obviously of NG. Anyway traces of a recent star formation history are evident in the spectra of S2G galaxies, which show a 4000 A break systematically shallower than in NG.

We present a systematic study of the sub-sample of Shakhbazyan groups (SHKs) covered by the Sloan Digital Sky Survey Data Release--5 (SDSS-5). SHKs probe an environment with characteristics which are intermediate between those of loose and very compact groups. Surprisingly, we found that several groups identifying algorithms (e.g. Berlind et al. 2006, Tago et al. 2008) miss this type of structures. Using the SDSS-5 spectroscopic data and the photometric redshifts derived in D'Abrusco et al. 2007, we identified possible group members in photometric redshift space and derived, for each group, several individual properties. We also combined pointed and stacked Rosat All Sky Survey data to investigate the X-ray luminosities of these systems. Our study confirms that the majority of groups are physical entities with richness in the range 3--13 galaxies, and properties ranging between those of loose and compact groups. We confirm that SHK groups are richer in early-type galaxies than the surrounding environment and the field, as expected from the morphology-density relation and from the selection of groups of red galaxies. Furthermore, our work supports the existence of two sub-classes of structures, the first one being formed by compact and isolated groups and the second formed by extended structures. We suggest that while the first class of objects dwells in less dense regions like the outer parts of clusters or the field, possibly sharing the properties of Hickson Compact Groups, the more extended structures represent a mixture of [core+halo] configurations and cores of rich clusters. X-ray luminosities for SHKs are generally consistent with these results and with the expectations for the L_X-sigma_v relation, but also suggest the velocity dispersions reported in literature are underestimated for some of the richest systems.

We present a method for the photometric selection of candidate quasars in multiband surveys. The method makes use of a priori knowledge derived from a subsample of spectroscopic confirmed QSOs to map the parameter space. The disentanglement of QSOs candidates and stars is performed in the colour space through the combined use of two algorithms, the Probabilistic Principal Surfaces and the Negative Entropy clustering, which are for the first time used in an astronomical context. Both methods have been implemented in the VONeural package on the Astrogrid VO platform. Even though they belong to the class of the unsupervised clustering tools, the performances of the method are optimized by using the available sample of confirmed quasars and it is therefore possible to learn from any improvement in the available "base of knowledge". The method has been applied and tested on both optical and optical plus near infrared data extracted from the visible SDSS and infrared UKIDSS-LAS public databases. In all cases, the experiments lead to high values of both efficiency and completeness, comparable if not better than the methods already known in the literature. A catalogue of optical candidate QSOs extracted from the SDSS DR7 Legacy photometric dataset has been produced and is publicly available at the URL voneural.na.infn.it/qso.html.

S.Co.P.E. is one of the four projects funded by the Italian Government in order to provide Southern Italy with a distributed computing infrastructure for fundamental science. Beside being aimed at building the infrastructure, S.Co.P.E. is also actively pursuing research in several areas among which astrophysics and observational cosmology. We shortly summarize the most significant results obtained in the first two years of the project and related to the development of middleware and Data Mining tools for the Virtual Observatory.

VO-Neural is the natural evolution of the Astroneural project which was started in 1994 with the aim to implement a suite of neural tools for data mining in astronomical massive data sets. At a difference with its ancestor, which was implemented under Matlab, VO-Neural is written in C++, object oriented, and it is specifically tailored to work in distributed computing architectures. We discuss the current status of implementation of VO-Neural, present an application to the classification of Active Galactic Nuclei, and outline the ongoing work to improve the functionalities of the package.

Modern photometric multiband digital surveys produce large amounts of data that, in order to be effectively exploited, need automatic tools capable to extract from photometric data an objective classification. We present here a new method for classifying objects in large multi-parametric photometric data bases, consisting of a combination of a clustering algorithm and a cluster agglomeration tool. The generalization capabilities and the potentialities of this approach are tested against the complexity of the Sloan Digital Sky Survey archive, for which an example of application is reported.

We present a new analysis of the Sloan Digital Sky Survey data aimed at producing a detailed map of the nearby (z < 0.5) universe. Using neural networks trained on the available spectroscopic base of knowledge we derived distance estimates for about 30 million galaxies distributed over ca. 8,000 sq. deg. We also used unsupervised clustering tools developed in the framework of the VO-Tech project, to investigate the possibility to understand the nature of each object present in the field and, in particular, to produce a list of candidate AGNs and QSOs.

In the framework of the European VO-Tech project, we are implementing new machine learning methods specifically tailored to match the needs of astronomical data mining. In this paper, we shortly present the methods and discuss an application to the Sloan Digital Sky Survey public data set. In particular, we discuss some preliminary results on the 3-D taxonomy of the nearby (z < 0.5) universe. Using neural networks trained on the available spectroscopic base of knowledge we derived distance estimates for ca. 30 million galaxies distributed over 8,000 sq. deg. We also use unsupervised clustering tools to investigate whether it is possible to characterize in broad morphological bins the nature of each object and produce a reliable list of candidate AGNs and QSOs.