We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

We present the theoretical background and detailed equations for the synchrotron emission of a shock wave propagating in a relativistic jet. We then show how the evolution of an outburst in this shock-in-jet scenario can be analytically described and parameterized to be fitted to multi-frequency lightcurves of galactic and extragalactic sources. This is done here for the first time with a completely physical description of the jet and the shocked gas, while previous studies used a more phenomenological approach based on the observed properties of the outbursts. Another interesting addition to previous work is the introduction of a low-energy cut of the electron energy distribution that allows for much more diverse synchrotron spectral shapes. To demonstrate and illustrate the new methodology, we present results of infrared-to-radio lightcurve fitting of a succession of outbursts observed in 1994 in the microquasar Cyg X-3. We find that the diversity of outbursts in shape, amplitude, frequency range and timescale can be fairly described by varying only the strength of the shock and its build-up distance from the apex of the jet. A rapid build-up results in high-frequency outbursts evolving on short timescales, while slowly evolving, low-frequency outbursts form and evolve further out in the jet. We conclude by outlining future developments, in particular the inclusion of the associated synchrotron self-Compton emission at X-rays and gamma-rays.

We derive the spectra of the cosmic X-ray background (CXB) and of the Galactic ridge X-ray emission (GRXE) in the ~20-200 keV range from the data of the IBIS instrument aboard the INTEGRAL satellite obtained during the four dedicated Earth-occultation observations of early 2006. We analyse the modulation of the IBIS/ISGRI detector counts induced by the passage of the Earth through the field of view of the instrument. Unlike previous studies, we do not fix the spectral shape of the various contributions, but model instead their spatial distribution and derive for each of them the expected modulation of the detector counts. The spectra of the diffuse emission components are obtained by fitting the normalizations of the model lightcurves to the observed modulation in different energy bins. The obtained CXB spectrum is consistent with the historic HEAO-1 results and falls slightly below the spectrum derived with Swift/BAT. A 10% higher normalization of the CXB cannot be completely excluded, but it would imply an unrealistically high albedo of the Earth. The derived spectrum of the GRXE confirms the presence of a minimum around 80 keV with improved statistics and yields an estimate of ~0.6 M_Sun for the average mass of white dwarfs in the Galaxy. The analysis also provides updated normalizations for the spectra of the Earth's albedo and the cosmic-ray induced atmospheric emission.

The Planck Low Frequency Instrument (LFI) will observe the Cosmic Microwave Background (CMB) by covering the frequency range 30-70 GHz in three bands. The primary instrument data source are the temperature samples acquired by the 22 radiometers mounted on the Planck focal plane. Such samples represent the scientific data of LFI. In addition, the LFI instrument generates the so called housekeeping data by sampling regularly the on-board sensors and registers. The housekeeping data provides information on the overall health status of the instrument and on the scientific data quality. The scientific and housekeeping data are collected on-board into telemetry packets compliant with the ESA Packet Telemetry standards. They represent the primary input to the first processing level of the LFI Data Processing Centre. In this work we show the software systems which build the LFI Level 1. A real-time assessment system, based on the ESA SCOS 2000 generic mission control system, has the main purpose of monitoring the housekeeping parameters of LFI and detect possible anomalies. A telemetry handler system processes the housekeeping and scientific telemetry of LFI, generating timelines for each acquisition chain and each housekeeping parameter. Such timelines represent the main input to the subsequent processing levels of the LFI DPC. A telemetry quick-look system allows the real-time visualization of the LFI scientific and housekeeping data, by also calculating quick statistical functions and fast Fourier transforms. The LFI Level 1 has been designed to support all the mission phases, from the instrument ground tests and calibration to the flight operations, and developed according to the ESA engineering standards.

This paper is part of the Prelaunch status LFI papers published on JINST (http://www.iop.org/EJ/journal/-page=extra.proc5/1748-0221). Planck's Low Frequency Instrument is an array of 22 pseudo-correlation radiometers at 30, 44, and 70 GHz. Before integrating the overall array assembly, a first set of tests has been performed for each radiometer chain assembly (RCA), consisting of two radiometers. In this paper, we describe Rachel, a software application which has been purposely developed and used during the RCA test campaign to carry out both near-realtime on-line data analysis and data storage (in FITS format) of the raw output from the radiometric chains.

The Level 1 of the Planck LFI Data Processing Centre (DPC) is devoted to the handling of the scientific and housekeeping telemetry. It is a critical component of the Planck ground segment which has to strictly commit to the project schedule to be ready for the launch and flight operations. In order to guarantee the quality necessary to achieve the objectives of the Planck mission, the design and development of the Level 1 software has followed the ESA Software Engineering Standards. A fundamental step in the software life cycle is the Verification and Validation of the software. The purpose of this work is to show an example of procedures, test development and analysis successfully applied to a key software project of an ESA mission. We present the end-to-end validation tests performed on the Level 1 of the LFI-DPC, by detailing the methods used and the results obtained. Different approaches have been used to test the scientific and housekeeping data processing. Scientific data processing has been tested by injecting signals with known properties directly into the acquisition electronics, in order to generate a test dataset of real telemetry data and reproduce as much as possible nominal conditions. For the HK telemetry processing, validation software have been developed to inject known parameter values into a set of real housekeeping packets and perform a comparison with the corresponding timelines generated by the Level 1. With the proposed validation and verification procedure, where the on-board and ground processing are viewed as a single pipeline, we demonstrated that the scientific and housekeeping processing of the Planck-LFI raw data is correct and meets the project requirements.

The Planck Low Frequency Instrument (LFI) is an array of 22 pseudo-correlation radiometers on-board the Planck satellite to measure temperature and polarization anisotropies in the Cosmic Microwave Background (CMB) in three frequency bands (30, 44 and 70 GHz). To calibrate and verify the performances of the LFI, a software suite named LIFE has been developed. Its aims are to provide a common platform to use for analyzing the results of the tests performed on the single components of the instrument (RCAs, Radiometric Chain Assemblies) and on the integrated Radiometric Array Assembly (RAA). Moreover, its analysis tools are designed to be used during the flight as well to produce periodic reports on the status of the instrument. The LIFE suite has been developed using a multi-layered, cross-platform approach. It implements a number of analysis modules written in RSI IDL, each accessing the data through a portable and heavily optimized library of functions written in C and C++. One of the most important features of LIFE is its ability to run the same data analysis codes both using ground test data and real flight data as input. The LIFE software suite has been successfully used during the RCA/RAA tests and the Planck Integrated System Tests. Moreover, the software has also passed the verification for its in-flight use during the System Operations Verification Tests, held in October 2008.

To asses stability against 1/f noise, the Low Frequency Instrument (LFI) onboard the Planck mission will acquire data at a rate much higher than the data rate allowed by its telemetry bandwith of 35.5 kbps. The data are processed by an onboard pipeline, followed onground by a reversing step. This paper illustrates the LFI scientific onboard processing to fit the allowed datarate. This is a lossy process tuned by using a set of 5 parameters Naver, r1, r2, q, O for each of the 44 LFI detectors. The paper quantifies the level of distortion introduced by the onboard processing, EpsilonQ, as a function of these parameters. It describes the method of optimizing the onboard processing chain. The tuning procedure is based on a optimization algorithm applied to unprocessed and uncompressed raw data provided either by simulations, prelaunch tests or data taken from LFI operating in diagnostic mode. All the needed optimization steps are performed by an automated tool, OCA2, which ends with optimized parameters and produces a set of statistical indicators, among them the compression rate Cr and EpsilonQ. For Planck/LFI the requirements are Cr = 2.4 and EpsilonQ <= 10% of the rms of the instrumental white noise. To speedup the process an analytical model is developed that is able to extract most of the relevant information on EpsilonQ and Cr as a function of the signal statistics and the processing parameters. This model will be of interest for the instrument data analysis. The method was applied during ground tests when the instrument was operating in conditions representative of flight. Optimized parameters were obtained and the performance has been verified, the required data rate of 35.5 Kbps has been achieved while keeping EpsilonQ at a level of 3.8% of white noise rms well within the requirements.

We present INTEGRAL light curves and spectra of the black-hole binary Cyg X-1 during a bright event that occurred in 2006 September, and which was simultaneous with a detection at 0.15-1 TeV energies by the MAGIC telescope. We analyse the hard X-ray emission from 18 to 700 keV with the INTEGRAL data taken on 2006 September 24-26. We present the light curves and fit the high energy spectrum with various spectral models. Despite variations in the flux by a factor of ~2 in the the 20-700 keV energy band, the shape of the energy spectrum remained remarkably stable. It is very well represented by an e-folded power law with the photon index of 1.4 and a high energy cut-off around 130-140 keV. The spectrum is also well described by thermal Comptonisation including a moderate reflection component, with a reflectionamplitude R around 0.4. The temperature of the hot Comptonising electrons is ~70 keV and their Thomson optical depth is ~2.5. These spectral properties are typical of those observed in the low/hard state. This shows that Cyg X-1 may stay in the low hard state at least up to the flux level of 2 Crab, which corresponds to ~2-3 percent of the Eddington luminosity. It is the first time a persistent high-mass black-hole binary is observed at a few percent of the Eddington luminosity with a stable low/hard state spectrum over a period of a few days. The TeV detection coincides with the peak of a small X-ray flare just after a very fast rise in hard X-ray flux. In contrast, the source remained undetected by MAGIC at the peak of a larger X-ray flare occurring one day later and corresponding to the maximum of the X-ray luminosity of the whole outburst. We do not find any obvious correlation between the X-ray and TeV emission.

We study the spectrum of the cosmic X-ray background (CXB) in energy range $\sim$5-100 keV. Early in 2006 the INTEGRAL observatory performed a series of four 30ksec observations with the Earth disk crossing the field of view of the instruments. The modulation of the aperture flux due to occultation of extragalactic objects by the Earth disk was used to obtain the spectrum of the Cosmic X-ray Background(CXB). Various sources of contamination were evaluated, including compact sources, Galactic Ridge emission, CXB reflection by the Earth atmosphere, cosmic ray induced emission by the Earth atmosphere and the Earth auroral emission. The spectrum of the cosmic X-ray background in the energy band 5-100 keV is obtained. The shape of the spectrum is consistent with that obtained previously by the HEAO-1 observatory, while the normalization is $\sim$10% higher. This difference in normalization can (at least partly) be traced to the different assumptions on the absolute flux from the Crab Nebulae. The increase relative to the earlier adopted value of the absolute flux of the CXB near the energy of maximum luminosity (20-50 keV) has direct implications for the energy release of supermassive black holes in the Universe and their growth at the epoch of the CXB origin.

The aim of this paper is to study the evolution of the broadband spectrum of one of the brightest and nearest quasars 3C 273. We analyze the data obtained during quasi-simultaneous INTEGRAL and XMM monitoring of the blazar 3C 273 in 2003--2005 in the UV, X-ray and soft gamma-ray bands and study the results in the context of the long-term evolution of the source. The 0.2-100 keV spectrum of the source is well fitted by a combination of a soft cut-off power law and a hard power law. No improvement of the fit is achieved if one replaces the soft cut-off power law by either a blackbody, or a disk reflection model. During the observation period the source has reached the historically softest state in the hard X-ray domain with a photon index $\Gamma=1.82\pm 0.01$. Comparing our data with available archived X-ray data from previous years, we find a secular evolution of the source toward softer X-ray emission (the photon index has increased by $\Delta\Gamma\simeq 0.3-0.4$ over the last thirty years). We argue that existing theoretical models have to be significantly modified to account for the observed spectral evolution of the source.

We discuss differences and similarities between jets powered by super-massive black holes in quasars and by stellar-mass black holes in microquasars. The comparison is based on multi-wavelength radio-to-infrared observations of the two active galactic nuclei 3C 273 and 3C 279, as well as the two galactic binaries GRS 1915+105 and Cyg X-3. The physical properties of the jet are derived by fitting the parameters of a shock-in-jet model simultaneously to all available observations. We show that the variable jet emission of galactic sources is, at least during some epochs, very similar to that of extra-galactic jets. As for quasars, their observed variability pattern can be well reproduced by the emission of a series of self-similar shock waves propagating down the jet and producing synchrotron outbursts. This suggests that the physical properties of relativistic jets is independent of the mass of the black hole.

Aims. The aim of this work is to identify and study spectral features in the quasar 3C 273 usually blended by its strong jet emission. Method. A historic minimum in the sub-millimetre emission of 3C 273 triggered coordinated multi-wavelength observations in June 2004. X-ray observations from the INTEGRAL, XMM-Newton and RXTE satellites are complemented by ground-based optical, infrared, millimetre and radio observations. The overall spectrum is used to model the infrared and X-ray spectral components. Results. Three thermal dust emission components are identified in the infrared. The dust emission on scales from 1 pc to several kpc is comparable to that of other quasars, as expected by AGN unification schemes. The observed weakness of the X-ray emission supports the hypothesis of a synchrotron self-Compton origin for the jet component. There is a clear soft-excess and we find evidence for a very broad iron line which could be emitted in a disk around a Kerr black hole. Other signatures of a Seyfert-like X-ray component are not detected.

We discuss accretion flows on massive black holes in which different elements of the flow (clumps) have velocities that may differ substantially. We estimate the consequence of collisions between these clumps as they come close to the central object and calculate the resulting radiation. We show that this radiation is similar to that observed in the optical to X-ray spectral domain in Seyfert galaxies and quasars. We also show that the large scale accretion is likely to be clumpy when arriving in the active region and that the clumps keep their identity between collisions.

In this paper we apply the reverberation method to determine the mass of the black hole in 3C273 from the Ly a and C iv emission lines using archival IUE observations. Following the standard assumptions of the method, we find a maximum-likelihood estimate of the mass of 6.59 10^9 Mo, with a 1 sigma confidence interval 5.69-8.27 10^9 Mo. This estimate is more than one order of magnitude larger than that obtained in a previous study using Balmer lines. We reanalyze the optical data and show that the method applied to the Ha, Hb, and Hg Balmer lines produce mass estimates lower by a factor 2.5, but already much larger than the previous estimate derived from the same lines. The finding of such a high mass in a face-on object is a strong indication that the gas motion is not confined to the accretion disk. The new mass estimate makes 3C273 accreting with an accretion rate about six times lower than the Eddington rate. We discuss the implications of our result for the broad-line-region size and black-hole mass vs luminosity relationships for the set of objects for which reverberation black-hole masses have been obtained. We find that, while objects with super-Eddington luminosities might theoretically be possible, their existence is not necessarily implied by this sample.

The INTEGRAL Science Data Centre (ISDC) routinely monitors the Near Real Time data (NRT) from the INTEGRAL satellite. A first scientific analysis is made in order to check for the detection of new, transient or highly variable sources in the data. Of primary importance for this work is the Interactive Quick Look Analysis (IQLA), which produces JEM-X and ISGRI images and monitors them for interesting astrophysical events

The true nature of the radio emitting material observed to be moving relativistically in quasars and microquasars is still unclear. The microquasar community usually interprets them as distinct clouds of plasma, while the extragalactic community prefers a shock wave model. Here we show that the synchrotron variability pattern of the microquasar GRS 1915+105 observed on 15 May 1997 can be reproduced by the standard shock model for extragalactic jets, which describes well the long-term behaviour of the quasar 3C 273. This strengthens the analogy between the two classes of objects and suggests that the physics of relativistic jets is independent of the mass of the black hole. The model parameters we derive for GRS 1915+105 correspond to a rather dissipative jet flow, which is only mildly relativistic with a speed of 0.60 c. We can also estimate that the shock waves form in the jet at a distance of about 1 AU from the black hole.

We present here a new determination of the ratio of energies radiated by active galactic nuclei and by stars and discuss the reasons for the apparently conflicting results found in previous studies. We conclude that the energy radiated by accretion processes onto super massive black holes is about 1 to 5% of the energy radiated by stars. We also estimate that the total mass accreted on average by a super massive black hole at the centre of a typical 10^11 Msol galaxy is of about 7 10^7 Msol.

In this paper we study the variability properties of the Lyman alpha and C IV emission lines in 3C273 using archival IUE observations. Our data show for the first time the existence of variability on time scales of several years. We study the spatial distribution and the velocity field of the emitting gas by performing detailed analyses on the line variability using correlations, 1D and 2D response functions, and principal component analysis. In both lines we find evidence for two components, one which has the dynamic properties of gas in Keplerian motion around a black hole with a mass of the order of 10^9 Mo, and one which is characterized by high, blue-shifted velocities at large lag. There is no indication of the presence of optically thick emission medium neither in the Lya, nor in the Civ response functions. The component characterized by blue-shifted velocities, which is comparatively much stronger in Civ than in Lya, is more or less compatible with being the result of gas falling towards the central black hole with free-fall acceleration. We propose however that the line emission at high, blue-shifted velocities is better explained in terms of entrainment of gas clouds by the jet. This gas is therefore probably collisionally excited as a result of heating due to the intense infrared radiation from the jet, which would explain the strength of this component in Civ relative to Lya. This phenomenon might be a signature of disk-jet interaction.

We present a phenomenological jet model which is able to reproduce well the observed variations of the submillimetre-to-radio emission of the bright quasar 3C 273 during the last 20 years. It is a generalization of the original shock model of Marscher & Gear (1985), which is now able to describe an accelerating or decelerating shock wave, in a curved, non-conical and non-adiabatic jet. The model defines the properties of a synchrotron outburst which is expected to be emitted by the jet material in a small region just behind the shock front. By a proper parameterization of the average outburst's evolution and of the peculiarities of individual outbursts, we are able to decompose simultaneously thirteen long-term light-curves of 3C 273 in a series of seventeen distinct outbursts. It is the first time that a model is so closely confronted to the long-term multi-wavelength variability properties of a quasar. The ability of the model to reproduce the very different shapes of the submillimetre-to-radio light curves of 3C 273 gives strong support to the shock model of Marscher & Gear (1985). Indirectly, it also reinforces the idea that the outbursts seen in the light-curves are physically linked to the distinct features observed to move along the jet with apparently superluminal velocities.

We present a new approach to derive the observed properties of synchrotron outbursts in relativistic jets. The idea is to use the very well sampled submillimetre-to-radio long-term light curves of 3C 273 to extract the spectral and temporal evolution of a typical outburst. The method consists in a decomposition of these light curves into a series of twelve self-similar flares. With a model-independent parameterization, we find that the obtained outburst's evolution is in good qualitative agreement with the expectations of shock models in relativistic jets. We then derive, by a second approach, the relevant parameters of three-stage shock models. We observe for the first time that the optically thin spectral index is steeper during the initial rising phase of the evolution than during the final declining phase as expected by the shock model of Marscher & Gear (1985). We obtain that this index flattens from alpha=-1.1 to alpha=-0.5, in good agreement with what is expected from a power law electron energy distribution with an index of -2.0. The observed flattening gives support to the idea that radiative (synchrotron and/or Compton) losses are the dominant cooling process of the electrons during the initial phase of the shock evolution. Two other results give us confidence in our decomposition: 1) the outbursts that we identify do well correspond to the VLBI components observed in the jet and 2) there is strong evidence that high-frequency peaking outbursts evolve faster than low-frequency peaking outbursts. We propose that this last correlation is related to the distance from the core of the jet at which the shock forms.

We present a wide multi-wavelength database of most observations of the quasar 3C 273 obtained during the last 30 years. This database is the most complete set of observations available for an active galactic nucleus (AGN). It contains nearly 20'000 observations grouped together into 70 light curves covering 16 orders of magnitude in frequency from the radio to the gamma-ray domain. The database is constituted of many previously unpublished observations and of most publicly available data gathered in the literature and on the World Wide Web (WWW). It is complete to the best of our knowledge, except in the optical (UBV) domain where we chose not to add all observations from the literature. In addition to the photometric data, we present the spectra of 3C 273 obtained by the International Ultraviolet Explorer (IUE) satellite. In the X-ray domain, we used the spectral fit parameters from the literature to construct the light curves. Apart from describing the data, we show the most representative light curves and the average spectrum of 3C 273. The database is available on the WWW in a homogeneous and clear form and we wish to update it regularly by adding new observations.

In 3C 273, ultraviolet flux and X-ray flux measured by BATSE are not well correlated, contrarily to predictions of several models, unless the X-ray flux lags the UV emission by 1.75 yr. The absence of observed correlation at small lag cannot be due to spectral variability. A Comptonizing corona model is however compatible with all UV and X-ray observations covering the BATSE period.

We apply a principal component analysis (PCA) to the spectra of each of the 18 Seyfert 1-like objects observed more than 15 times by the international ultraviolet explorer (IUE) from 1978 until the end of 1991. PCA allows us to decompose the Lyman alpha and the C IV (1549 A) emission-lines of each object into two components with uncorrelated variations. We find that the principal component describes correlated continuum and line variations, whereas the rest component reveals the line-parts that do not vary in tune with the continuum. A cross-correlation analysis in NGC 5548 reveals that variable line features in the rest component follow the continuum by about 25 days. The symmetry usually observed in the component's line profiles excludes that the velocity field of the line-emitting region is dominated by radial motion. In some objects, the principal component reveals a clearly double-peaked line profile. The results suggest that the variable line feature in the principal and the rest component is respectively emitted in the inner and outer parts of an extended broad-line region (BLR). The fraction of the line that varies together with the continuum seems to decrease significantly with increasing luminosity. We interpret this trend as a consequence of a larger BLR in more luminous objects. Finally, we propose a solution to the C IV/Ly_a ratio problem without referring to optically thin emission-line clouds.