The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

In our previous analysis we investigated the large-scale environment of two samples of radio galaxies (RGs) in the local Universe (i.e. with redshifts z<0.15), classified as FR I and FR II on the basis of their radio morphology. The analysis was carried out using i) extremely homogeneous catalogs and ii) a new method, known as cosmological overdensity, to investigate their large-scale environments. We concluded that, independently by the shape of their radio extended structure, RGs inhabit galaxy-rich large-scale environments with similar characteristics and richness. In the present work, we first highlight additional advantages of our procedure, that does not suffer cosmological biases and/or artifacts, and then we carry out an additional statistical test to strengthen our previous results. We also investigate properties of RG environments using those of the cosmological neighbors. We find that large-scale environments of both FRIs and FRIIs are remarkably similar and independent on the properties of central RG. Finally, we highlight the importance of comparing radio sources in the same redshift bins to obtain a complete overview of their large-scale environments.

In more than 17 years of operation in space INTEGRAL telescope has accumulated large data set that contains records of hard X-ray and soft gamma-ray astronomical sources. These data can be re-used in the context of multi-wavelength or multi-messenger studies of astronomical sources and have to be preserved on long time scales. We present a scientific validation of an interactive online INTEGRAL data analysis system, which generates publication-quality high-level data products: sky images, spectra and light-curves in response to user queries that define analysis parameters, such as source position, time and energy interval and binning. The data products can be requested via a web browser interface or via Application Programming Interface (API). The products for the ISGRI instrument of INTEGRAL are generated using the Offline Science Analysis (OSA) software which is provided by the instrument teams, but reorganized to preserve and re-use intermediate analysis products, ensuring that frequently requested results are available without delay. The platform is portable and can be deployed in any compatible infrastructure. We report the functionalities and performance of the online data analysis system by reproducing the benchmark INTEGRAL results on different types of sources, including bright steady and transient Galactic sources, and bright and weak variable extra-galactic sources. We compare the results obtained with the online data analysis system with previously published results on these sources. We consider the INTEGRAL online data analysis as a demonstrator of more general web-based "data analysis as a service" approach that provides a promising solution for preservation and maintenance of data analysis tools of astronomical telescopes on (multi)decade long time scales and facilitates combination of data in multi-wavelength and multi-messenger studies of astronomical sources.

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

The MAGIC telescopes are one of the three major IACTs (Imaging Atmospheric Cherenkov Telescopes) for observation of gamma rays in the TeV regime currently operative. MAGIC functions since 2003, and has published data from more than 60 sources, mostly blazars. MAGIC already provides astronomical \texttt{.fits} files with basic final scientific products such as spectral energy distributions, light curves and skymaps from published results. In future, the format of the files can be complemented with further relevant information to the community: a) by including the full multi-wavelength dataset enclosed in a publication, b) providing data in alternative easy-to-use formats such as ASCII or ECSV, which are accessible with other commonly used packages such as \texttt{astropy} or \texttt{gammapy}. Finally, besides high level products, activities have started to provide photon event lists and instrument response functions in a format such that scientists within and outside the community are allowed to perform higher level analysis. A second aim is to provide a full legacy of MAGIC data. This contribution will illustrate the achievements and plans of this activity.

The new generation of deep photometric surveys requires unprecedentedly precise shape and photometry measurements of billions of galaxies to achieve their main science goals. At such depths, one major limiting factor is the blending of galaxies due to line-of-sight projection, with an expected fraction of blended galaxies of up to 50%. Current deblending approaches are in most cases either too slow or not accurate enough to reach the level of requirements. This work explores the use of deep neural networks to estimate the photometry of blended pairs of galaxies in monochrome space images, similar to the ones that will be delivered by the Euclid space telescope. Using a clean sample of isolated galaxies from the CANDELS survey, we artificially blend them and train two different network models to recover the photometry of the two galaxies. We show that our approach can recover the original photometry of the galaxies before being blended with $\sim$7% accuracy without any human intervention and without any assumption on the galaxy shape. This represents an improvement of at least a factor of 4 compared to the classical SExtractor approach. We also show that forcing the network to simultaneously estimate a binary segmentation map results in a slightly improved photometry. All data products and codes will be made public to ease the comparison with other approaches on a common data set.

PKS 0438$-$436 at a redshift of $z=2.856$ has been previously recognized as possessing perhaps the most luminous known synchrotron jet. Little is known about this source since the maximum elevation above the horizon is low for the Very Large Array (VLA). We present the first VLA radio image that detects the radio lobes. We use both the 151 MHz luminosity, as a surrogate for the isotropic radio lobe luminosity, and the lobe flux density from the radio image to estimate a long term, time averaged, jet power, $\overline{Q} =1.5\pm 0.7 \times 10^{47} \rm{ergs~s^{-1}}$. We analyze two deep optical spectra with strong broad emission lines and estimate the thermal bolometric luminosity of the accretion flow, $L_{\rm{bol}} = 6.7 \pm 3.0 \times 10^{46} \rm{ergs~s^{-1}}$. The ratio, $\overline{Q}/L_{\rm{bol}} = 3.3 \pm 2.6 $, is at the limit of this empirical metric of jet dominance seen in radio loud quasars and this is the most luminous accretion flow to have this limiting behavior. Despite being a very luminous blazar, it previously had no $\gamma$-ray detections (EGRET, AGILE or FERMI) until December 11 - 13 2016 (54 hours) when FERMI detected a flare that we analyze here. The isotropic apparent luminosity from 100 MeV - 100 GeV rivals the most luminous detected blazar flares (averaged over 18 hours), $\sim 5-6 \times 10^{49} \rm{ergs~s^{-1}}$. The $\gamma$-ray luminosity varies over time by two orders of magnitude, highlighting the extreme role of Doppler abberation and geometric alignment in producing the inverse Compton emission.

The gamma-ray emission in broad-line blazars is generally explained as inverse Compton (IC) radiation of relativistic electrons in the jet scattering optical-UV photons from the Broad Line Region (BLR), the so-called BLR External Compton scenario. We test this scenario on the Fermi gamma-ray spectra of 106 broad-line blazars detected with the highest significance or largest BLR, by looking for cut-off signatures at high energies compatible with gamma-gamma interactions with BLR photons. We do not find evidence for the expected BLR absorption. For 2/3 of the sources, we can exclude any significant absorption ($\tau_{max}<1$), while for the remaining 1/3 the possible absorption is constrained to be 1.5-2 orders of magnitude lower than expected. This result holds also dividing the spectra in high and low-flux states, and for powerful blazars with large BLR. Only 1 object out of 10 seems compatible with substantial attenuation ($\tau_{max}>5$). We conclude that for 9 out of 10 objects, the jet does not interact with BLR photons. Gamma-rays seem either produced outside the BLR most of the time, or the BLR is ~100x larger than given by reverberation mapping. This means that i) External Compton on BLR photons is disfavoured as the main gamma-ray mechanism, vs IC on IR photons from the torus or synchrotron self-Compton; ii) the Fermi gamma-ray spectrum is mostly intrinsic, determined by the interaction of the particle distribution with the seed-photons spectrum; iii) without suppression by the BLR, broad-line blazars can become copious emitters above 100 GeV, as demonstrated by 3C454.3. We expect the CTA sky to be much richer of broad-line blazars than previously thought.

Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected throughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{\rm radio}\propto\sigma_v^{3.3\pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.

Euclid is a Europe-led cosmology space mission dedicated to a visible and near infrared survey of the entire extra-galactic sky. Its purpose is to deepen our knowledge of the dark content of our Universe. After an overview of the Euclid mission and science, this contribution describes how the community is getting organized to face the data analysis challenges, both in software development and in operational data processing matters. It ends with a more specific account of some of the main contributions of the Swiss Science Data Center (SDC-CH).

We present a study on galaxy detection and shape classification using topometric clustering algorithms. We first use the DBSCAN algorithm to extract, from CCD frames, groups of adjacent pixels with significant fluxes and we then apply the DENCLUE algorithm to separate the contributions of overlapping sources. The DENCLUE separation is based on the localization of pattern of local maxima, through an iterative algorithm which associates each pixel to the closest local maximum. Our main classification goal is to take apart elliptical from spiral galaxies. We introduce new sets of features derived from the computation of geometrical invariant moments of the pixel group shape and from the statistics of the spatial distribution of the DENCLUE local maxima patterns. Ellipticals are characterized by a single group of local maxima, related to the galaxy core, while spiral galaxies have additional ones related to segments of spiral arms. We use two different supervised ensemble classification algorithms, Random Forest, and Gradient Boosting. Using a sample of ~ 24000 galaxies taken from the Galaxy Zoo 2 main sample with spectroscopic redshifts, and we test our classification against the Galaxy Zoo 2 catalog. We find that features extracted from our pipeline give on average an accuracy of ~ 93%, when testing on a test set with a size of 20% of our full data set, with features deriving from the angular distribution of density attractor ranking at the top of the discrimination power.

We present the results of a new search for galaxy-scale strong lensing systems in CFHTLS Wide. Our lens-finding technique involves a preselection of potential lens galaxies, applying simple cuts in size and magnitude. We then perform a Principal Component Analysis of the galaxy images, ensuring a clean removal of the light profile. Lensed features are searched for in the residual images using the clustering topometric algorithm DBSCAN. We find 1098 lens candidates that we inspect visually, leading to a cleaned sample of 109 new lens candidates. Using realistic image simulations we estimate the completeness of our sample and show that it is independent of source surface brightness, Einstein ring size (image separation) or lens redshift. We compare the properties of our sample to previous lens searches in CFHTLS. Including the present search, the total number of lenses found in CFHTLS amounts to 678, which corresponds to ~4 lenses per square degree down to i=24.8. This is equivalent to ~ 60.000 lenses in total in a survey as wide as Euclid, but at the CFHTLS resolution and depth.

We study the contribution of thermal and non-thermal processes to the inverse Compton emission of the radio galaxy M 87 by modelling its broad-band emission. Through this we aim to derive insight into where within the AGN the X-ray, gamma-ray, and VHE emission is produced. We have analysed all available INTEGRAL IBIS/ISGRI data on M 87, spanning almost 10 years, to set an upper limit to the average hard X-ray flux of $f(20 - 60 \rm \, keV) < 3\times 10^{-12}$ $\rm \, erg \, cm^{-2} \, s^{-1}$, using several techniques beyond the standard analysis which are also presented here. We also analysed hard X-ray data from Suzaku/PIN taken late November 2006, and we report the first hard X-ray detection of M 87 with a flux of $f(20 - 60 \rm \, keV) = 10^{-11}\rm \, erg \, cm^{-2} \, s^{-1} $. In addition we analyse data from Fermi/LAT, INTEGRAL/JEM-X, and Suzaku/XIS. We collected historical radio/IR/optical and VHE data and combined them with the X-ray and gamma-ray data, to create broad-band spectral energy distributions for the average low-flux state and the flaring state. The resulting spectral energy distributions are modelled by applying a single-zone SSC model with a jet angle of theta = 15 degrees. We also show that modelling the core emission of M 87 using a single-zone synchrotron self-Compton model does represent the SED, suggesting that the core emission is dominated by a BL Lac type AGN core. Using SED modelling we also show that the hard X-ray emission detected in 2006 is likely due to a flare of the jet knot HST-1, rather than being related to the core.

Aims. The high energy spectrum of 3C 273 is usually understood in terms of inverse-Compton emission in a relativistic leptonic jet. This model predicts variability patterns and delays that could be tested with simultaneous observations from the radio to the GeV range. Methods. The instruments IBIS, SPI, JEM-X on board INTEGRAL, PCA on board RXTE, and LAT on board Fermi have enough sensitivity to follow the spectral variability of 3C 273 from the keV to the GeV. We looked for correlations between the different energy bands, including radio data at 37 GHz collected at the Mets\"ahovi Radio Observatory and built quasi-simultaneous multiwavelength spectra in the high energy domain when the source is flaring either in the X-rays or in the {\gamma} rays. Results. Both temporal and spectral analysis suggest a two-component model to explain the complete high energy spectrum. X-ray emission is likely dominated by a Seyfert-like component while the {\gamma}-ray emission is dominated by a blazar-like component produced by the relativistic jet. The variability of the blazar-like component is discussed, comparing the spectral parameters in the two different spectral states. Changes of the electron Lorentz factor are found to be the most likely source of the observed variability.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of radio-loud Active Galactic Nuclei. For a summary, we refer to the paper.

The Large Area Telescope on board the Fermi Gamma-ray Space Telescope detected a strong gamma-ray flare on 2011 May 15 from a source identified as 4C 49.22, a flat spectrum radio quasar also known as S4 1150+49. This blazar, characterised by a prominent radio-optical-X-ray jet, was in a low gamma-ray activity state during the first years of Fermi observations. Simultaneous observations during the quiescent, outburst and post-flare gamma-ray states were obtained by Swift, Planck and optical-IR-radio telescopes (INAOE, Catalina CSS, VLBA, Metsahovi). The flare is observed from microwave to X-ray bands with correlated variability and the Fermi, Swift and Planck data for this FSRQ show some features more typical of BL Lac objects, like the synchrotron peak in the optical band that outshines the thermal blue-bump emission, and the X-ray spectral softening. Multi-epoch VLBA observations show the ejection of a new component close in time with the GeV gamma-ray flare. The radio-to-gamma-ray spectral energy distribution is modeled and fitted successfully for the outburst and the post-flare epochs using either a single flaring blob with two emission processes (synchrotron self Compton, and external-radiation Compton), and a two-zone model with SSC-only mechanism.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Blazars are highly variable sources over timescales that can be as low as minutes. This is the case of the High Energy Peaked BL Lac (HBL) objects showing strong variability in X-rays, which highly correlate with that of the TeV emission. The degree of this correlation is still debated, particularly when the flaring activity is followed down to very short time scales. This correlation could challenge the synchrotron-self-Compton scenario in which one relativistic electron population dominates the entire radiative output. We argue that the LOFT Large Area Detector (10 m$^2$, LAD), thanks to its unprecedented timing capability, will allow us to detect the X-ray counterpart (2-50 keV) of the very fast variability observed at TeV energies, sheding light on the nature of X-TeV connection. We will discuss the test case of PKS 2155-304, showing as it would be possible to look for any X-ray variability occurring at very short timescales, never explored so far. This will put strong constraints on the size and the location of any additional electron population in the multi-zone scenario. Under this perspective, LOFT and the CTA observatories, planned to operate in the same time frame, will allow us to investigate in depth the connection between X-ray and TeV emissions. We also discuss the potentialities of LOFT in measuring the change in spectral curvature of the synchrotron spectra in HBLs which will make possible to directly study the mechanism of acceleration of highly energetic electrons. LOFT timing capability will be also promising in the study of Flat Spectrum Radio Quasars (FSRQs) with flux $\ge 1$ mCrab. Constraints to the location of the high energy emission will be given by temporal investigation on second timescale and spectral trend analysis on minute timescales. This represents a further link with CTA because of the rapid (unexpected) TeV emission recently detected in some FSRQs.

The high energy spectral shape of 3C 273 is usually understood in terms of Inverse-Compton emission in a relativistic leptonic jet. This model predicts variability patterns and delays which could be tested if simultaneous observations are available from the infrared to the GeV range. The instruments IBIS, SPI, JEM-X on board INTEGRAL, PCA on board RXTE and LAT on board Fermi have enough sensitivity to follow the spectral variability from the keV to the GeV and to compare them with model predictions. We are presenting preliminary results on the high energy spectrum of 3C 273 and its variability and compare these results to predictions. We found that a single component is not able to adequately fit the multiwavelength spectrum (5 keV - 10 GeV) when the statistics in Fermi data is sufficient to constraint {\gamma}-ray emission. This suggests that the X-rays emission is not originated in the jet. A possible explanation could be that X-rays are dominated by a Seyfert-like thermal inverse Compton.

We report the detection of a new gamma-ray source in the Fermi-LAT sky using a source detection tool based on the Minimal Spanning Tree algorithm. The source, not reported in previous LAT catalogues but very recently observed in the X-rays and optical bands, is characterized by an increasing gamma-ray activity in 2012 June-September, reaching a weekly peak flux of (3.3+-0.6)*10^-7 photons cm^-2 s^-1. A search for a possible counterpart provides indication that it can be associated with the radio source NVSS J141828+354250 whose optical SDSS colours are typical of a blazar.

The Density Based Spatial Clustering of Applications with Noise (DBSCAN) is a topometric algorithm used to cluster spatial data that are affected by background noise. For the first time, we propose the use of this method for the detection of sources in gamma-ray astrophysical images obtained from the Fermi-LAT data, where each point corresponds to the arrival direction of a photon. We investigate the detection performance of the gamma-ray DBSCAN in terms of detection efficiency and rejection of spurious clusters, using a parametric approach, and exploring a large volume of the gamma-ray DBSCAN parameter space. By means of simulated data we statistically characterize the gamma-ray DBSCAN, finding signatures that differentiate purely random fields, from fields with sources. We define a significance level for the detected clusters, and we successfully test this significance with our simulated data. We apply the method to real data, and we find an excellent agreement with the results obtained with simulated data. We find that the gamma-ray DBSCAN can be successfully used in the detection of clusters in gamma-ray data. The significance returned by our algorithm is strongly correlated with that provided by the Maximum Likelihood analysis with standard Fermi-LAT software, and can be used to safely remove spurious clusters. The positional accuracy of the reconstructed cluster centroid compares to that returned by standard Maximum Likelihood analysis, allowing to look for astrophysical counterparts in narrow regions, minimizing the chance probability in the counterpart association. We find that gamma-ray DBSCAN is a powerful tool in the detection of clusters in gamma-ray data, this method can be used both to look for point-like sources, and extended sources, and can be potentially applied to any astrophysical field related with detection of clusters in data.

We present more than three years of observations at different frequencies, from radio to high-energy gamma-rays, of the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.585). This source is the first NLS1 detected at energies above 100 MeV and therefore can be considered the prototype of this emerging new class of gamma-ray emitting active galactic nuclei (AGN). The observations performed from 2008 August 1 to 2011 December 31 confirmed that PMN J0948+0022 generates a powerful relativistic jet, able to develop an isotropic luminosity at gamma-rays of the order of 10^48 erg s^-1, at the level of powerful quasars. The evolution of the radiation emission of this source in 2009 and 2010 followed the canonical expectations of relativistic jets, with correlated multiwavelength variability (gamma-rays followed by radio emission after a few months), but it was difficult to retrieve a similar pattern in the light curves of 2011. The comparison of gamma-ray spectra before and including 2011 data suggested that there was a softening of the high-energy spectral slope. We selected five specific epochs to be studied by modelling the broad-band spectrum, characterised by an outburst at gamma-rays or very low/high flux at other wavelengths. The observed variability can largely be explained either by changes in the injected power, the bulk Lorentz factor of the jet or the electron spectrum. The characteristic time scale of doubling/halving flux ranges from a few days to a few months, depending on the frequency and the sampling rate. The shortest doubling time scale at gamma-rays is 2.3+-0.5 days. These small values underline the need of highly-sampled multiwavelength campaigns to better understand the physics of these sources.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

We present simultaneous Planck, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, hard X-ray, and gamma-ray bands. Our unique data set has allowed us to demonstrate that the selection method strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by Fermi-LAT, whereas ~40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are still below the gamma-ray detection limit even after integrating 27 months of Fermi-LAT data. The radio to sub-mm spectral slope of blazars is quite flat up to ~70GHz, above which it steepens to <\alpha>~-0.65. BL Lacs have significantly flatter spectra than FSRQs at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (\nupS) in the SED of FSRQs is the same in all the blazar samples with <\nupS>=10^13.1 Hz, while the mean inverse-Compton peak frequency, <\nupIC>, ranges from 10^21 to 10^22 Hz. The distributions of \nupS and of \nupIC of BL Lacs are much broader and are shifted to higher energies than those of FSRQs and strongly depend on the selection method. The Compton dominance of blazars ranges from ~0.2 to ~100, with only FSRQs reaching values >3. Its distribution is broad and depends strongly on the selection method, with gamma-ray selected blazars peaking at ~7 or more, and radio-selected blazars at values ~1, thus implying that the assumption that the blazar power is dominated by high-energy emission is a selection effect. Simple SSC models cannot explain the SEDs of most of the gamma-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-LAT may instead be consistent with SSC emission. Our data challenge the correlation between bolometric luminosity and \nupS predicted by the blazar sequence.

(abridged) We analyse the spectral and pulse properties of the 11 Hz transient accreting pulsar, IGR J17480-2446, in the globular cluster Terzan 5, considering all the available RXTE, Swift and INTEGRAL observations performed between October and November, 2010. By measuring the pulse phase evolution we conclude that the NS spun up at an average rate of <nu_dot>=1.48(2)E-12 Hz/s, compatible with the accretion of the Keplerian angular momentum of matter at the inner disc boundary. Similar to other accreting pulsars, the stability of the pulse phases determined by using the second harmonic component is higher than that of the phases based on the fundamental frequency. Under the assumption that the second harmonic is a good tracer of the neutron star spin frequency, we successfully model its evolution in terms of a luminosity dependent accretion torque. If the NS accretes the specific Keplerian angular momentum of the in-flowing matter, we estimate the inner disc radius to lie between 47 and 93 km when the luminosity attains its peak value. Smaller values are obtained if the interaction between the magnetic field lines and the plasma in the disc is considered. The phase-averaged spectrum is described by thermal Comptonization of photons with energy of ~1 keV. A hard to soft state transition is observed during the outburst rise. The Comptonized spectrum evolves from a Comptonizing cloud at an electron temperature of ~20 keV towards an optically denser cloud at kT_e~3 keV. At the same time, the pulse amplitude decreases from 27% to few per cent and becomes strongly energy dependent. We discuss various possibilities to explain such a behaviour, proposing that at large accretion luminosities a significant fraction of the in-falling matter is not channelled towards the magnetic poles, but rather accretes more evenly onto the NS surface.

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an emphasis on proper, physical modelling of the synchrotron bump using multiplecomponents. Planck ERCSC data also suggest that the original accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data indicate that gamma-ray emission must originate in the same shocks that produce the radio emission.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

We report about the multiwavelength campaign on the Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z = 0.5846) performed in 2010 July-September and triggered by high activity as measured by Fermi/LAT. The peak luminosity in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^48 erg/s, a level comparable to the most powerful blazars. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

The broad-band spectral distributions of non-thermal sources, such as those of several known blazars, are well described by a log-parabolic fit. The second degree term in these fits measures the curvature in the spectrum. In this paper we investigate whether the curvature parameter observed in the spectra of the synchrotron emission can be used as a fingerprint of stochastic acceleration. As a first approach we use the multiplicative Central Limit theorem to show how fluctuations in the energy gain result in the broadening of the spectral shape, introducing a curvature into the energy distribution. Then, by means of a Monte-Carlo description, we investigate how the curvature produced in the electron distribution is linked to the diffusion in momentum space. To get a more generic description we turn to the diffusion equation in momentum space. We first study some "standard" scenarios, in order to understand the conditions that make the curvature in the spectra significant, and the relevance of the cooling over the acceleration process. We try to quantify the correlation between the curvature and the diffusive process in the pre-equilibrium stage, and investigate how the transition between the Klein-Nishina and the Thompson regime, in Inverse Compton cooling, can determine the curvature in the distribution at the equilibrium. We apply these results to some observed trends, such as the anticorrelation between the peak energy and the curvature term observed in the spectra of Mrk 421, and a sample of BL Lac objects whose synchrotron emission peaks at X-ray energies.

We present the results of our analysis devoted to the research of sources emitting in the energy bands surveyed by both the Swift-BAT and the Fermi-LAT telescopes. We cross-correlate the Fermi-LAT 1-year point source catalogue (1FGL) of {\gamma}-ray sources and the second Palermo BAT catalogue (2PBC) of hard X-ray sources, establishing a correspondence between sources when their error boxes overlap. We also extract the significance value in the BAT 15-150 keV map, obtained using a dedicated software for the reduction of BAT data, in the direction of the 1FGL sources and take into account those above the significance threshold {\sigma} = 3. We obtain a sample of common sources emitting in both the hard X- and the {\gamma}-ray energy bands and evaluate its content in galactic and extragalactic objects. We assess the fraction of unidentified sources and describe in greater detail the properties of two of them, 1FGL J0137.8+5814 and 1FGL J2056.7+4938, supporting their classification as blazars after the analysis of their broad-band spectral energy distribution. We discuss the blazar content of the collected 1FGL-2PBC sources: we build its redshift distibution and compare it with that of the whole blazar population as reported in the second edition of the BZCAT blazar catalogue.

We report on a multiwavelength campaign on the radio-loud Narrow-Line Seyfert 1 (NLS1) Galaxy PMN J0948+0022 (z=0.5846) performed in 2010 July-September and triggered by a high-energy gamma-ray outburst observed by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. The peak flux in the 0.1-100 GeV energy band exceeded, for the first time in this type of source, the value of 10^-6 ph cm^-2 s^-1, corresponding to an observed luminosity of 10^48 erg s^-1. Although the source was too close to the Sun position to organize a densely sampled follow-up, it was possible to gather some multiwavelength data that confirmed the state of high activity across the sampled electromagnetic spectrum. The comparison of the spectral energy distribution of the NLS1 PMN J0948+0022 with that of a typical blazar - like 3C 273 - shows that the power emitted at gamma rays is extreme.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

Narrow-Line Seyfert 1 (NLS1) class of active galactic nuclei (AGNs) is generally radio-quiet, but a small percent of them are radio-loud. The recent discovery by Fermi/LAT of high-energy gamma-ray emission from 4 NLS1s proved the existence of relativistic jets in these systems. It is therefore important to study this new class of gamma-ray emitting AGNs. Here we report preliminary results about the observations of the July 2010 gamma-ray outburst of PMN J0948+0022, when the source flux exceeded for the first time 10^-6 ph cm^-2 s^-1 (E > 100 MeV).

HESS J1632-478 is an extended and still unidentified TeV source in the galactic plane. In order to identify the source of the very high energy emission and to constrain its spectral energy distribution, we used a deep observation of the field obtained with XMM-Newton together with data from Molonglo, Spitzer and Fermi to detect counterparts at other wavelengths. The flux density emitted by HESS J1632-478 peaks at very high energies and is more than 20 times weaker at all other wavelengths probed. The source spectrum features two large prominent bumps with the synchrotron emission peaking in the ultraviolet and the external inverse Compton emission peaking in the TeV. HESS J1632-478 is an energetic pulsar wind nebula with an age of the order of 10^4 years. Its bolometric (mostly GeV-TeV) luminosity reaches 10% of the current pulsar spin down power. The synchrotron nebula has a size of 1 pc and contains an unresolved point-like X-ray source, probably the pulsar with its wind termination shock.

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

We present the results of a deep spectral analysis of all Swift observations of Mrk 421 between April 2006 and July 2006, when it reached its highest X-ray flux recorded until the end of 2006. We completed this data set with other historical X-ray observations. We used the full data set to investigate the correlation between the spectral parameters. We found a signature of stochastic acceleration in the anticorrelation between the peak energy (Ep) of the spectral energy distribution (SED) and the spectral curvature parameter (b). We found signature of energetic budget of the jet in the correlation between the peak flux of the SED (S p) and Ep. Moreover, using simultaneous Swift UVOT/XRT/BAT data, we demonstrated, that during the strongest flares, the UV-to-X-ray emission from Mrk 421 requires that the curved electron distribution develops a low energy power-law tail. The observed spectral curvature and its anticorrelation with Ep is consistent with both stochastic acceleration or energy-dependent acceleration probability mechanisms, whereas the power-law slope of XRT-UVOT data is close to that inferred from the GRBs X-ray afterglow and in agreement with the universal first-order relativistic shock acceleration models. This scenario implies that magnetic turbulence may play a twofold role: spatial diffusion relevant to the first order process and momentum diffusion relevant to the second order process.

We have studied the quasi-simultaneous Spectral Energy Distributions (SED) of 48 LBAS blazars, detected within the three months of the LAT Bright AGN Sample (LBAS) data taking period, combining Fermi and Swift data with radio NIR-Optical and hard-X/gamma-ray data. Using these quasi-simultaneous SEDs, sampling both the low and the high energy peak of the blazars broad band emission, we were able to apply a diagnostic tool based on the estimate of the peak frequencies of the synchrotron (S) and Inverse Compton (IC) components. Our analysis shows a Fermi blazars' divide based on the peak frequencies of the SED. The robust result is that the Synchrotron Self Compton (SSC) region divides in two the plane were we plot the peak frequency of the synchrotron SED vs the typical Lorentz factor of the electrons most contributing to the synchrotron emission and to the inverse Compton process. Objects within or below this region, radiating likely via the SSC process, are high-frequency-peaked BL Lac object (HBL), or low/intermediate-frequency peaked BL Lac object (LBL/IBL). All of the IBLs/LBLs within or below the SSC region are not Compton dominated. The objects lying above the SSC region, radiating likely via the External radiation Compton (ERC) process, are Flat Spectrum Radio Quasars and IBLs/LBLs. All of the IBLs/LBLs in the ERC region show a significant Compton dominance.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

We study the quasi-simultaneous Spectral Energy Distributions (SED) of 48 LBAS blazars, collected within three months of the Fermi LAT Bright AGN Sample (LBAS) data taking period, combining Fermi and Swift data with radio NIR-Optical and hard-X/gamma-ray data. We have used these SEDs to characterize the peak position and intensity of both the low and the high-energy features of blazar spectra. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical, alpha_ro, and optical to X-ray, alpha_ox, spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency is positioned between 10^(12.5) and 10^(14.5) Hz in broad-lined FSRQs and between 10^(13) and 10^(17) Hz in featureless BL Lacertae objects. We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. More complex models involving External Compton Radiation or multiple SSC components are required to reproduce the overall SEDs and the observed spectral variability.

We present results from a deep spectral analysis of all the Swift observations of Mrk 421 from April 2006 to July 2006, when it reached its largest X-ray flux recorded until 2006. The peak flux was about 85 milli-Crab in the 2.0-10.0 keV band, with the peak energy (Ep) of the spectral energy distribution (SED) laying often at energies larger than 10 keV. We performed spectral analysis of the Swift observations investigating the trends of the spectral parameters in terms of acceleration and energetic features phenomenologically linked to the SSC model parameters, predicting their effects in the gamma-ray band, in particular the spectral shape expected in the Fermi Gamma-ray Space Telescope-LAT band. We confirm that the X-ray spectrum is well described by a log-parabolic distribution close to Ep, with the peak flux of the SED (Sp) being correlated with Ep, and Ep anti-correlated with the curvature parameter (b). During the most energetic flares the UV-to-soft-X-ray spectral shape requires an electron distribution spectral index s about 2.3. Present analysis shows that the UV-to-X-ray emission from Mrk 421 is likely to be originated by a population of electrons that is actually curved, with a low energy power-law tail. The observed spectral curvature is consistent both with stochastic acceleration or energy dependent acceleration probability mechanisms, whereas the power-law slope form XRT-UVOT data is very close to that inferred from the GRBs X-ray afterglow and in agreement with the universal first-order relativistic shock acceleration models. This scenario hints that the magnetic turbulence may play a twofold role: spatial diffusion relevant to the first order process and momentum diffusion relevant to the second order process.

Many of the extragalactic sources detected in $\gamma$ rays at TeV energies are BL Lac objects. In particular, they belong to the subclass of ``high frequency peaked BL Lacs" (HBLs), as their spectral energy distributions exhibit a first peak in the X-ray band. At a closer look, their X-ray spectra appear to be generally curved into a log-parabolic shape. In a previous investigation of Mrk 421, two correlations were found between the spectral parameters. One involves the height $S_p$ increasing with the position $E_p$ of the first peak; this was interpreted as a signature of synchrotron emission from relativistic electrons. The other involves the curvature parameter $b$ decreasing as $E_p$ increases; this points toward statistical/stochastic acceleration processes for the emitting electrons. We analyse X-ray spectra of several TeV HBLs to pinpoint their behaviours in the $E_p-S_p$ and $E_p-b$ planes and to compare them with Mrk 421. We perfom X-ray spectral analyses of a sample of 15 BL Lacs. We report the whole set of observations obtained with the \sax, \xmm and \swf satellites between 29/06/96 and 07/04/07. We focus on five sources (PKS 0548-322, 1H 1426+418, Mrk 501, 1ES 1959+650, PKS2155-304) whose X-ray observations warrant detailed searching of correlations or trends. Within our database, we find that four out of five sources, namely PKS 0548-322, 1H 1426+418, Mrk 501 and 1ES 1959+650, follow similar trends as Mrk 421 in the $E_p-S_p$ plane, while PKS 2155-304 differs. As for the $E_p-b$ plane, all TeV HBLs follow a similar behaviour. The trends exhibited by Mrk 421 appear to be shared by several TeV HBLs, such as to warrant discussing predictions from the X-ray spectral evolution to that of TeV emissions.

We developed a source detection algorithm based on the Minimal Spanning Tree (MST), that is a graph-theoretical method useful for finding clusters in a given set of points. This algorithm is applied to gamma-ray bidimensional images where the points correspond to the arrival direction of photons, and the possible sources are associated with the regions where they clusterize. Some filters to select these clusters and to reduce the spurious detections are introduced. An empirical study of the statistical properties of MST on random fields is carried in order to derive some criteria to estimate the best filter values. We introduce also two parameters useful to verify the goodness of candidate sources. To show how the MST algorithm works in the practice, we present an application to an EGRET observation of the Virgo field, at high galactic latitude and with a low and rather uniform background, in which several sources are detected.

We studied the temporal and spectral evolution of the synchrotron emission from the high energy peaked BL Lac object 1E 1207.9+3945. Two recent observations have been performed by the XMM-Newton and Swift satellites; we carried out X-ray spectral analysis for both of them, and photometry in optical-ultraviolet filters for the Swift one. Combining the results thus obtained with archival data we built the long-term X-ray light curve, spanning a time interval of 26 years, and the Spectral Energy Distribution (SED) of this source. The light curve shows a large flux increasing, about a factor of six, in a time interval of a few years. After reaching its maximum in coincidence with the XMM-Newton pointing in December 2000 the flux decreased in later years, as revealed by Swift. The very good statistics available in the 0.5-10 keV XMM-Newton X-ray spectrum points out a highly significant deviation from a single power law. A log-parabolic model with a best fit curvature parameter of 0.25 and a peak energy at ~1 keV describes well the spectral shape of the synchrotron emission. The simultaneous fit of Swift UVOT and XRT data provides a milder curvature (b~0.1) and a peak at higher energies (~15 keV), suggesting a different state of source activity. In both cases UVOT data support the scenario of a single synchrotron emission component extending from the optical/UV to the X-ray band. New X-ray observations are important to monitor the temporal and spectral evolution of the source; new generation gamma-ray telescopes like AGILE and GLAST could for the first time detect its inverse Compton emission.

Gamma Ray Bursts (GRBs) show evidence of different spectral shapes, light curves, duration, host galaxies and they explode within a wide redshift range. However, the most of them seems to follow very tight correlations among some observed quantities relating to their energetic. If true, these correlations have significant implications on burst physics, giving constraints on theoretical models. Moreover, several suggestions have been made to use these correlations in order to calibrate GRBs as standard candles and to constrain the cosmological parameters. We investigate the cosmological relation between low energy $\alpha$ index in GRBs prompt spectra and the redshift $z$. We present a statistical analysis of the relation between the total isotropic energy $E_{iso}$ and the peak energy $E_p$ (also known as Amati relation) in GRBs spectra searching for possible functional biases. Possible implications on the $E_{iso}$ vs $E_p$ relation of the $\alpha$ vs $(1+z)$ correlation are evaluated. We used MonteCarlo simulations and the boostrap method to evaluate how large are the effects of functional biases on the $E_{iso}$ vs $E_p$. We show that high values of the linear correlation coefficent, up to about 0.8, in the $E_{iso}$ vs $E_p$ relation are obtained for random generated samples of GRBs, confirming the relevance of functional biases. Astrophysical consequences from $E_{iso}$ vs $E_p$ relation are then to be revised after a more accurate and possibly bias free analysis.

Almost the totality of the bright foreground sources in the WMAP CMB maps are blazars, a class of sources that show usually also X-ray emission. However, 23 objects in a flux-limited sample of 140 blazars of the WMAP catalog (first year) were never reported before as X-ray sources. We present here the results of 41 Swift observations which led to the detection of all these 23 blazars in the 0.3-10 keV band. We conclude that all micro-wave selected blazars are X-ray emitters and that the distribution of the micro-wave to X-ray spectral slope $\alpha_{mu x}$ of LBL blazars is very narrow, confirming that the X-ray flux of most blazars is a very good estimator of their micro-wave emission. The X-ray spectral shape of all the objects that were observed long enough to allow spectral analysis is flat and consistent with inverse Compton emission within the commonly accepted view where the radiation from blazars is emitted in a Sychrotron-Inverse-Compton scenario. We predict that all blazars and most radio galaxies above the sensitivity limit of the WMAP and of the Planck CMB missions are X-ray sources detectable by the present generation of X-ray satellites. An hypothetical all-sky soft X-ray survey with sensitivity of approximately $10^{-15}$ erg/s would be crucial to locate and remove over 100,000 blazars from CMB temperature and polarization maps and therefore accurately clean the primordial CMB signal from the largest population of extragalactic foreground contaminants.

We present an empirical spectral modelling of the high energy emission of the anomalous X-ray pulsar 4U 0142+614, based on simultaneous Swift and INTEGRAL observations from X to gamma-ray energies. We adopted models contained in the XSPEC analysis package, as well as models based on recent theoretical studies, and restricted ourselves to those combinations of up to three components which produce a good fit while requiring the lowest number of free parameters. Only three models were found to fit satisfactorily the 0.5-250keV spectrum of 4U 0142+614: i) a ~0.4keV blackbody and two power-laws, ii) a resonant cyclotron scattering model plus a power-law and iii) two log-parabolic functions. We found that only the latter two models do not over-predict the infrared/optical emission observed simultaneously from this AXP, and only the log-parabolic functions can naturally account for the upper limits set by COMPTEL in the gamma-ray range. A possible interpretation of the two log-parabolae in terms of inverse Compton scattering of soft X-ray photons by very energetic particles is discussed.

BL Lac objects undergo strong flux variations involving considerable changes in their spectral shapes. We specifically investigate the X-ray spectral evolution of Mrk 421 over a time span of about nine years. We aim at statistically describing and physically understanding the large spectral changes in X rays observed in Mrk 421 over this time span. We perform a homogeneous spectral analysis of a wide data set including archived observations with ASCA, BeppoSax, RXTE, as well as published and unpublished XMM-Newton data. The presence of uncertainties is taken into account in our correlation analysis. The significance of the correlations found and possible spurious effects are studied with Monte Carlo simulations. We find that the Mrk421 spectral energy distribution (SED) has a lower peak at energies that vary in the range, 0.1-10 keV while its X-ray spectrum is definitely curved. Parameterizing the X-ray spectra with a log-parabolic model, we find a positive correlation between the position and the height of the SED peak. In addition, we find a negative trend of the spectral curvature parameter vs. the SED peak energy. We show that these relations between the spectral parameters are consistent with statistical or stochastic acceleration of the emitting particles, and provide insight into the physical processes occurring in BL Lac nuclei.

Context: We present the results of a set of observations of nine TeV detected BL Lac objects performed by the XRT and UVOT detectors on board the Swift satellite between March and December 2005. Aims: We are mainly interested in measuring the spectral parameters, and particularly the intrinsic curvature in the X-ray band. Methods: We perform X-ray spectral analysis of observed BL Lac TeV objects using either a log-parabolic or a simple power-law model . Results: We found that many of the objects in our sample do show significant spectral curvature, whereas those having the peak of the spectral energies distribution at energies lower than ~0.1 keV show power law spectra. In these cases, however, the statistics are generally low thus preventing a good estimate of the curvature. Simultaneous UVOT observations are important to verify how X-ray spectra can be extrapolated at lower frequencies and to search for multiple emission components. Conclusions: The results of our analysis are useful for the study of possible signatures of statistical acceleration processes predicting intrinsically curved spectra and for modelling the SED of BL Lacertae objects up to TeV energies where a corresponding curvature is likely to be present.

We present the results of a spectral analysis of 5 Swift XRT and UVOT observations of the BL Lac object PKS 0548-322 carried out over the period April-June 2005. The X-ray flux of this high energy peaked BL Lac (HBL) source was found to be approximately constant at a level of F(2-10 keV) ~ 4x10^-11 erg cm^-2 s^-1, a factor of 2 brighter than when observed by BeppoSAX in 1999 and close to the maximum intensity reported in the Einstein Slew Survey. The very good statistics obtained in the 0.3-10 keV Swift X-ray spectrum allowed us to detect highly significant deviations from a simple power law spectral distribution. A log-parabolic model describes well the X-ray data and gives a best fit curvature parameter of 0.18 and peak energy in the Spectral Energy Distribution of about 2 keV. The UV spectral data from Swift UVOT join well with a power law extrapolation of the soft X-ray data points suggesting that the same component is responsible for the observed emission in the two bands. The combination of synchrotron peak in the X-ray band and high intensity state confirms PKS 0548-322 as a prime target for TeV observations. X-ray monitoring and coordinated TeV campaigns are highly advisable.

We present the results of a series of Swift and quasi simultaneous ground-based infra-red observations of the blazar 3C 454.3 carried out in April-May 2005 when the source was 10 to 30 times brighter than previously observed. We found 3C 454.3 to be very bright and variable at all frequencies covered by our instrumentation. The broad-band Spectral Energy Distribution (SED) shows the usual two-bump shape (in nu-nu f(nu) space) with the Infra-red, optical and UV data sampling the declining part of the synchrotron emission that, even during this extremely large outburst, had its maximum in the far-infrared. The X-ray spectral data from the XRT and BAT instruments are flat and due to inverse Compton emission. The remarkable SED observed implies that at the time of the Swift pointings 3C 454.3 was one of the brightest objects in the extragalactic sky with a gamma-ray emission similar or brighter than that of 3C 279 when observed in a high state by EGRET. Time variability in the optical-UV flux is very different from that in the X-ray data: while the first component varied by about a factor two within a single exposure, but remained approximately constant between different observations, the inverse Compton component did not vary on short time-scales but changed by more than a factor of 3 between observations separated by a few days. This different dynamical behaviour illustrates the need to collect simultaneous multi-frequency data over a wide range of time-scales to fully constrain physical parameters in blazars.

Curved broad-band spectral distributions of non-thermal sources like blazars are described well by a log-parabolic (LP) law where the second degree term measures the curvature. LP energy spectra can be obtained for relativistic electrons by means of a statistical acceleration mechanism whose probability of acceleration depends on energy. In this paper we compute the spectra radiated by an electron population via synchrotron (S) and Synchro-Self Compton(SSC) processes to derive the relations between the LP parameters. These spectra were obtained by means of an accurate numerical code. We found that the ratio between the curvature parameters of the S spectrum to that of the electrons is equal to about 0.2 instead of 0.25, the value foreseen in the delta approximation. Inverse Compton spectra are also intrinsically curved and can be approximated by a log-parabola only in limited ranges. The curvature parameter, estimated around the SED peak, may vary from a lower value than that of the S spectrum up to that of emitting electrons depending on whether the scattering is in the Thomson or in the Klein-Nishina regime. We applied this analysis to computing the SSC emission from the BL Lac object Mkn 501 during the large flare of April 1997. We fit simultaneous BeppoSAX and CAT data and reproduced intensities and spectral curvatures of both components with good accuracy. The large curvature observed in the TeV range was found to be mainly intrinsic, and therefore did not require a large pair production absorption against the extragalactic background. We regard this finding as an indication that the Universe is more transparent at these energies than previously assumed by several models found in the literature. This conclusion is supported by recent detection of two relatively high redshift blazars with H.E.S.S.

In this paper we present a broad-band study of the X-ray emission of the blazar PKS1830-211 based on Chandra and Integral observations. Notwithstanding the high redshift (z=2.507), it is a bright X-ray source (F(2-10 keV)~10^{-11} erg cm^{-2} s^{-1}), due to gravitational lensing by an intervening galaxy at z=0.89. Previous X-ray observations attribute the observed absorption at E<2 keV to the lensing galaxy. Our analysis, although not in contrast with this hypothesis, suggests also the possibility of an intrinsic (ionized) absorption, taking place at the front side of the jet. This scenario is also supported by some evidence, in the same data, of a feature observed at 2.15 keV which can be interpreted as a blueshifted iron line (v/c ~ 0.18). The SED of PKS1830-211 can be well modelled by combining a Synchrotron Self-Compton component and an external source of photons to be scattered up to \gamma-ray energies by relativistic electrons moving outward in the jet. The main source of low energy photons is a dust torus at the temperature of 10^3 K as expected in MeV blazars.

We present a numerical code, written in C, which can be used to simulate or to analyze the emission of Blazars over the entire electromagnetic spectrum. Our code can reproduce the following features: synchrotron emission, inverse Compton emission (Thomson Klein-Nishina regime) external Compton emission, accretion disk variability using a Cellular Automata algorithm, temporal evolution of the emitting plasma energy distribution, flaring phenomena, light curves in the rest and in the observer frame (taking account for time crossing effects). In this paper we will show mainly the accretion disk simulation, and the implications in the External Compton scenario.

We present preliminary results of a statistical analysis obtained with a sample of blazars observed at the Perugia Astronomical Observatory since 1992. We briefly show how these statistical results can be useful to discriminate faint variable sources against the background noise. This technique, together with the more traditional ones, may be used to discover and identify high-energy point sources.