CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

The main limiting factor of cosmological analyses based on thermal Sunyaev-Zel'dovich (SZ) cluster statistics comes from the bias and systematic uncertainties that affect the estimates of the mass of galaxy clusters. High-angular resolution SZ observations at high redshift are needed to study a potential redshift or morphology dependence of both the mean pressure profile and of the mass-observable scaling relation used in SZ cosmological analyses. The NIKA2 camera is a new generation continuum instrument installed at the IRAM 30-m telescope. With a large field of view, a high angular resolution and a high-sensitivity, the NIKA2 camera has unique SZ mapping capabilities. In this paper, we present the NIKA2 SZ large program, aiming at observing a large sample of clusters at redshifts between 0.5 and 0.9, and the characterization of the first cluster oberved with NIKA2.

NIKA2 is a dual-band millimetric continuum camera of 2900 Kinetic Inductance Detectors (KID), operating at $150$ and $260\,\rm{GHz}$, installed at the IRAM 30-meter telescope. We present the performance assessment of NIKA2 after one year of observation using a dedicated point-source calibration method, referred to as the \emph{baseline} method. Using a large data set acquired between January 2017 and February 2018 that span the whole range of observing elevations and atmospheric conditions encountered at the IRAM 30-m telescope, we test the stability of the performance parameters. We report an instantaneous field of view (FOV) of 6.5' in diameter, filled with an average fraction of $84\%$ and $90\%$ of valid detectors at $150$ and $260\,\rm{GHz}$, respectively. The beam pattern is characterized by a FWHM of $17.6'' \pm 0.1''$ and $11.1''\pm 0.2''$, and a beam efficiency of $77\% \pm 2\%$ and $55\% \pm 3\%$ at $150$ and $260\,\rm{GHz}$, respectively. The rms calibration uncertainties are about $3\%$ at $150\,\rm{GHz}$ and $6\%$ at $260\,\rm{GHz}$. The absolute calibration uncertainties are of $5\%$ and the systematic calibration uncertainties evaluated at the IRAM 30-m reference Winter observing conditions are below $1\%$ in both channels. The noise equivalent flux density (NEFD) at $150$ and $260\,\rm{GHz}$ are of $9 \pm 1\, \rm{mJy}\cdot s^{1/2}$ and $30 \pm 3\, \rm{mJy}\cdot s^{1/2}$. This state-of-the-art performance confers NIKA2 with mapping speeds of $1388 \pm 174$ and $111 \pm 11 \,\rm{arcmin}^2\cdot \rm{mJy}^{-2}\cdot \rm{h}^{-1}$ at $150$ and $260\,\rm{GHz}$. With these unique capabilities of fast dual-band mapping at high (better that 18'') angular resolution, NIKA2 is providing an unprecedented view of the millimetre Universe.

The NIKA2 polarization channel at 260 GHz (1.15 mm) has been proposed primarily to observe galactic star-forming regions and probe the critical scales between 0.01-0.05 pc at which magnetic field lines may channel the matter of interstellar filaments into growing dense cores. The NIKA2 polarimeter consists of a room temperature continuously rotating multi-mesh HWP and a cold polarizer that separates the two orthogonal polarizations onto two 260 GHz KIDs arrays. We describe in this paper the preliminary results obtained during the most recent commissioning campaign performed in December 2018. We concentrate here on the analysis of the extended sources, while the observation of compact sources is presented in a companion paper [12]. We present preliminary NIKA2 polarization maps of the Crab nebula. We find that the integrated polarization intensity flux measured by NIKA2 is consistent with expectations.In terms of polarization angle, we are still limited by systematic uncertainties that will be further investigated in the forthcoming commissioning campaigns.

We present here the polarization system of the NIKA camera and give a summary of the main results obtained and performed studies on Orion and the Crab nebula. The polarization system was equipped with a room temperature continuously rotating multi-mesh half wave plate and a grid polarizer facing the NIKA cryostat window. NIKA even though less sensitive than NIKA2 had polarization capability in both 1 and 2 millimiter bands. NIKA polarization observations demonstrated the ability of such a technology in detecting the polarization of different targets, compact and extended sources like the Crab nebula and Orion Molecular Cloud region OMC-1. These measurements together with the developed techniques to deal with systematics, opened the way to the current observations of NIKA2 in polarization that will provide important advances in the studies of galactic and extra-galactic emission and magnetic fields.

The new NIKA2 camera at the IRAM 30m radiotelescope was used to observe three known debris disks in order to constrain the SED of their dust emission in the millimeter wavelength domain. We have found that the spectral index between the two NIKA2 bands (1mm and 2mm) is consistent with the Rayleigh-Jeans regime (lambda^{-2}), unlike the steeper spectra (lambda^{-3}) measured in the submillimeter-wavelength domain for two of the three disks $-$ around the stars Vega and HD107146. We provide a succesful proof of concept to model this spectral inversion in using two populations of dust grains, those smaller and those larger than a grain radius a0 of 0.5mm. This is obtained in breaking the slope of the size distribution and the functional form of the absorption coefficient of the standard model at a0. The third disk - around the star HR8799 - does not exhibit this spectral inversion but is also the youngest.

The steps taken to tailor to NIKA2 observations the Scanamorphos algorithm (initially developed to subtract low-frequency noise from Herschel on-the-fly observations) are described, focussing on the consequences of the different instrument architecture and observation strategy. The method, making the most extensive use of the redundancy built in the multi-scan coverage with large arrays of a given region of the sky, is applicable to extended sources, while the pipeline is so far optimized for compact sources. An example of application is given. A related tool to build consistent broadband SEDs from 60 microns to 2 mm, combining Herschel and NIKA2 data, has also been developed. Its main task is to process the data least affected by low-frequency noise and coverage limitations (i.e. the Herschel data) through the same transfer function as the NIKA2 data, simulating the same scan geometry and applying the same noise and atmospheric signal as extracted from the 1 mm and 2 mm data.

We have conducted a multiwavelength survey of 42 radio loud narrow-1ine Seyfert 1 galaxies (RLNLS1s), selected by searching among all the known sources of this type and omitting those with steep radio spectra. We analyse data from radio frequencies to X-rays, and supplement these with information available from online catalogs and the literature in order to cover the full electromagnetic spectrum. This is the largest known multiwavelength survey for this type of source. We detected 90% of the sources in X-rays and found 17% at gamma rays. Extreme variability at high energies was also found, down to timescales as short as hours. In some sources, dramatic spectral and flux changes suggest interplay between a relativistic jet and the accretion disk. The estimated masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than those of blazars, while the accretion luminosities span a range from $\sim 0.01$ to $\sim 0.49$ times the Eddington limit, similar to those of quasars. The distribution of the calculated jet power spans a range from $\sim 10^{42.6}$ to $\sim 10^{45.6}$ erg s$^{-1}$, generally lower than quasars and BL Lac objects, but partially overlapping with the latter. Once normalised by the mass of the central black holes, the jet power of the three types of active galactic nuclei are consistent with each other, indicating the scalability of the jet. Despite the observational differences, the central engine of RLNLS1s is apparently quite similar to that of blazars. The historical difficulties in finding radio-loud narrow-line Seyfert 1 galaxies might be due to their low power and to intermittent jet activity.

Understanding where and when the mass of stars is determined is one of the fundamental, mostly unsolved, questions in astronomy. Here, we present the first results of GASTON, the Galactic Star Formation with NIKA2 large programme on the IRAM 30m telescope, that aims to identify new populations of low-brightness sources to tackle the question of stellar mass determination across all masses. In this paper, we focus on the high-mass star formation part of the project, for which we map a $\sim2$ deg$^2$ region of the Galactic plane around $l=24^\circ$ in both 1.2 mm and 2.0 mm continuum. Half-way through the project, we reach a sensitivity of 3.7 mJy/beam at 1.2mm. Even though larger than our target sensitivity of 2 mJy, the current sensitivity already allows the identification of a new population of cold, compact sources that remained undetected in any (sub-)mm Galactic plane survey so far. In fact, about 25% of the $\sim 1600$ compact sources identified in the 1.2 mm GASTON image are new detections. We present a preliminary analysis of the physical properties of the GASTON sources as a function of their evolutionary stage, arguing for a potential evolution of the mass distribution of these sources with time.

Clusters of galaxies, the largest bound objects in the Universe, constitute a cosmological probe of choice, which is sensitive to both dark matter and dark energy. Within this framework, the Sunyaev-Zel'dovich (SZ) effect has opened a new window for the detection of clusters of galaxies and for the characterization of their physical properties such as mass, pressure and temperature. NIKA, a KID-based dual band camera installed at the IRAM 30-m telescope, was particularly well adapted in terms of frequency, angular resolution, field-of-view and sensitivity, for the mapping of the thermal and kinetic SZ effect in high-redshift clusters. In this paper, we present the NIKA cluster sample and a review of the main results obtained via the measurement of the SZ effect on those clusters: reconstruction of the cluster radial pressure profile, mass, temperature and velocity.

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 45 galaxy clusters covering a wide mass range at high redshift in order to re-calibrate some of the tools needed for the cosmological exploitation of SZ surveys. We present the second cluster analysis of this program, targeting one of the faintest sources of the sample in order to tackle the difficulties in data reduction for such faint, low-SNR clusters. In this study, the main challenge is the precise estimation of the contamination by sub-millimetric point sources, which greatly affects the tSZ map of the cluster. We account for this contamination by performing a joint fit of the SZ signal and of the flux density of the compact sources. A prior knowledge of these fluxes is given by the adjustment of the SED of each source using data from both NIKA2 and the \textit{Herschel} satellite. The first results are very promising and demonstrate the possibility to estimate thermodynamic properties with NIKA2, even in a compact cluster heavily contaminated by point sources.

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

Filaments are a key step on the path that leads from molecular clouds to star formation. However, their characteristics are heavily debated, and the exact processes that lead to their formation and fragmentation into dense cores still remain to be fully understood. We aim at characterising the mass, kinematics, and stability against gravitational collapse of a statistically significant sample of filaments in the Orion B molecular cloud, which is renown for its very low star formation efficiency. We characterise the gas column densities and kinematics over a field of 1.9 deg$^2$, using C$^{18}$O(J=1-0) data from the IRAM-30m large programme ORION-B. Using two different Hessian-based filters, we extract and compare two filamentary networks, each containing over 100 filaments. Independent of the extraction method, the filaments have widths of 0.12$\pm$0.04 pc, and show a wide range of linear (1 - 100 $M_{\odot}$pc$^{-1}$) and volume densities (2.10$^3$ - 2.10$^5$ cm$^{-3}$). Compared to previous studies, the filament population is dominated by low-density, thermally sub-critical structures, suggesting that most of the identified filaments are not collapsing to form stars. In fact, only ~1% of the Orion B cloud mass covered by our observations can be found in super-critical, star-forming filaments, explaining the low star formation efficiency of the region. The velocity profiles observed across the filaments show quiescence in the centre, and coherency in the plane of the sky, despite being mostly supersonic. The filaments in Orion B apparently belong to a continuum which contains a few elements comparable to already studied star-forming filaments as well as many lower-density, gravitationally unbound structures. This comprehensive study of the Orion B filaments shows that the mass fraction in super-critical filaments is a key factor in determining star formation efficiency.

NIKA2 is a dual-band millimetric camera of thousands of Kinetic Inductance Detectors (KID) installed at the IRAM 30-meter telescope in the Spanish Sierra Nevada. The instrument commissioning was completed in September 2017, and NIKA2 is now open to the scientific community and will operate for the next decade. NIKA2 has well-adapted instrumental design and performance to produce high-resolution maps of the thermal Sunyaev-Zel'dovich (SZ) effect toward intermediate and high redshift galaxy clusters. Moreover, it benefits from a guaranteed time large program dedicated to mapping a representative sample of galaxy clusters via SZ and that includes X-ray follow-ups. The main expected outputs of the SZ large program are the constraints on the redshift evolution of the pressure profile and the mass-observable relation. The first SZ mapping of a galaxy cluster with NIKA2 was produced, as part of the SZ large program. We found a sizable impact of the intracluster medium dynamics on the integrated SZ observables. This shows NIKA2 capabilities for the precise characterisation of the mass-observable relation that is required for accurate cosmology with galaxy clusters.

The Crab nebula is a supernova remnant exhibiting a highly polarized synchrotron radiation at radio and millimeter wavelengths. It is the brightest source in the microwave sky with an extension of 7 by 5 arcminutes and commonly used as a standard candle for any experiment which aims at measuring the polarization of the sky. Though its spectral energy distribution has been well characterized in total intensity, polarization data are still lacking at millimetre wavelengths. We report in this paper high resolution (18 arcsec FWHM) observations of the Crab nebula in total intensity and linear polarization at 150 GHz with the NIKA camera. NIKA, operated at the IRAM 30 m telescope from 2012 to 2015, is a camera made of Lumped Element Kinetic Inductance Detectors (LEKIDs) observing the sky at 150 and 260 GHz. From these observations we are able to reconstruct the spatial distribution of the polarization degree and angle of the Crab nebula, which is found to be compatible with previous observations at lower and higher frequencies. Averaging across the source and using other existing data sets we find that the Crab nebula polarization angle is consistent with being constant over a wide range of frequencies with a value of -87.7$^\circ$ +- 0.3 in Galactic coordinates. We also present the first estimation of the Crab nebula spectral energy distribution polarized flux in a wide frequency range: 30-353 GHz. Assuming a single power law emission model we find that the polarization spectral index $\beta_{pol}$ = - 0.347 +- 0.026 is compatible with the intensity spectral index $\beta$ = - 0.323 +- 0.001.

The complete characterization of the pressure profile of high-redshift galaxy clusters, from their core to their outskirts, is a major issue for the study of the formation of large-scale structures. It is essential to constrain a potential redshift evolution of both the slope and scatter of the mass-observable scaling relations used in cosmology studies based on cluster statistics. In this paper, we present the first thermal Sunyaev-Zel'dovich (tSZ) mapping of a cluster from the sample of the NIKA2 SZ large program that aims at constraining the redshift evolution of cluster pressure profiles and the tSZ-mass scaling relation. We have observed the galaxy cluster PSZ2 G144.83+25.11 at redshift $z=0.58$ with the NIKA2 camera, a dual-band (150 and 260 GHz) instrument operated at the IRAM 30-meter telescope. We identify a thermal pressure excess in the south-west region of PSZ2 G144.83+25.11 and a high redshift sub-millimeter point source that affect the intracluster medium (ICM) morphology of the cluster. The NIKA2 data are used jointly with tSZ data acquired by the MUSTANG, Bolocam and $Planck$ experiments in order to non-parametrically set the best constraints on the electronic pressure distribution from the cluster core ($\rm{R} \sim 0.02 \rm{R_{500}}$) to its outskirts ($\rm{R} \sim 3 \rm{R_{500}} $). We investigate the impact of the over-pressure region on the shape of the pressure profile and on the constraints on the integrated Compton parameter $\rm{Y_{500}}$. A hydrostatic mass analysis is also performed by combining the tSZ-constrained pressure profile with the deprojected electronic density profile from XMM-$Newton$. This allows us to conclude that the estimates of $\rm{Y_{500}}$ and $\rm{M_{500}}$ obtained from the analysis with and without masking the disturbed ICM region differ by 65 and 79% respectively. (abridged)

The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds.

The New IRAM KID Arrays 2 (NIKA2) consortium has just finished installing and commissioning a millimetre camera on the IRAM 30 m telescope. It is a dual-band camera operating with three frequency multiplexed kilo-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK, designed to observe the intensity and polarisation of the sky at 260 and 150 GHz (1.15 and 2 mm). NIKA2 is today an IRAM resident instrument for millimetre astronomy, such as Intra Cluster Medium from intermediate to distant clusters and so for the follow-up of Planck satellite detected clusters, high redshift sources and quasars, early stages of star formation and nearby galaxies emission. We present an overview of the instrument performance as it has been evaluated at the end of the commissioning phase.

Sub-structures in the hot gas of galaxy clusters are related to their formation history and to the astrophysical processes at play in the intracluster medium (ICM). The thermal Sunyaev-Zel'dovich (tSZ) effect is directly sensitive to the line-of-sight integrated ICM pressure, and is thus particularly adapted to study ICM sub-structures. We apply structure-enhancement filtering algorithms to high resolution tSZ observations of distant clusters, in order to search for pressure discontinuities, compressions, as well as secondary peaks in the ICM. The same filters are applied to synthetic tSZ images extracted from RHAPSODY-G hydrodynamic simulations, in order to better interpret the extracted features. We also study the noise propagation trough the filters and quantify the impact of systematic effects, point source residuals being identified as the dominant potential contaminant. In 3 of our 6 NIKA clusters we identify features at high S/N that show clear evidence for merger events. In MACSJ0717 (z=0.55), three strong pressure gradients are observed on the E, SE and W sectors, and two main peaks in the pressure distribution are identified. We observe a lack of tSZ compact structure in the cool-core cluster PSZ1G045.85 (z=0.61), and a tSZ gradient ridge dominates in the SE. In the highest z cluster, CLJ1227 (z=0.89), we detect a ~45" (360 kpc) long ridge pressure gradient associated with a secondary pressure peak in the W region. Our results show that current tSZ facilities have now reached the angular resolution and sensitivity to allow an exploration of the details of pressure sub-structures in clusters, even at high z. This opens the possibility to quantify the impact of the dynamical state on the relation between the tSZ signal and the mass of clusters, which is important when using tSZ clusters to test cosmological models. This work also marks the first NIKA cluster sample data release.

Millimeter-wave continuum astronomy is today an indispensable tool for both general Astrophysics studies and Cosmology. General purpose, large field-of-view instruments are needed to map the sky at intermediate angular scales not accessible by the high-resolution interferometers and by the coarse angular resolution space-borne or ground-based surveys. These instruments have to be installed at the focal plane of the largest single-dish telescopes. In this context, we have constructed and deployed a multi-thousands pixels dual-band (150 and 260 GHz, respectively 2mm and 1.15mm wavelengths) camera to image an instantaneous field-of-view of 6.5arc-min and configurable to map the linear polarization at 260GHz. We are providing a detailed description of this instrument, named NIKA2 (New IRAM KID Arrays 2), in particular focusing on the cryogenics, the optics, the focal plane arrays based on Kinetic Inductance Detectors (KID) and the readout electronics. We are presenting the performance measured on the sky during the commissioning runs that took place between October 2015 and April 2017 at the 30-meter IRAM (Institut of Millimetric Radio Astronomy) telescope at Pico Veleta. NIKA2 has been successfully deployed and commissioned, performing in-line with the ambitious expectations. In particular, NIKA2 exhibits FWHM angular resolutions of around 11 and 17.5 arc-seconds at respectively 260 and 150GHz. The NEFD (Noise Equivalent Flux Densities) demonstrated on the maps are, at these two respective frequencies, 33 and 8 mJy*sqrt(s). A first successful science verification run has been achieved in April 2017. The instrument is currently offered to the astronomical community during the coming winter and will remain available for at least the next ten years.

Galaxy clusters constitute a major cosmological probe. However, Planck 2015 results have shown a weak tension between CMB-derived and cluster-derived cosmological parameters. This tension might be due to poor knowledge of the cluster mass and observable relationship. As for now, arcmin resolution Sunyaev-Zeldovich (SZ) observations ({\it e.g.} SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters ($z<0.2$). For high redshift clusters ($z>0.5$) high resolution and high sensitivity SZ observations are needed. With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain) is particularly well adapted for these observations. The NIKA2 SZ observation program will map a large sample of clusters (50) at redshifts between 0.5 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.}

Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases that constitute giant molecular clouds. We aim at using multiple tracers (sensitive to different physical processes) to segment a molecular cloud into physically/chemically similar regions (rather than spatially connected components). We use a machine learning clustering method (the Meanshift algorithm) to cluster pixels with similar molecular emission, ignoring spatial information. Simple radiative transfer models are used to interpret the astrophysical information uncovered by the clustering. A clustering analysis based only on the J=1-0 lines of 12CO, 13CO and C18O reveals distinct density/column density regimes (nH~100, 500, and >1000 cm-3), closely related to the usual definitions of diffuse, translucent and high-column-density regions. Adding two UV-sensitive tracers, the (1-0) lines of HCO+ and CN, allows us to distinguish two clearly distinct chemical regimes, characteristic of UV-illuminated and UV-shielded gas. The UV-illuminated regime shows overbright HCO+ and CN emission, which we relate to photochemical enrichment. We also find a tail of high CN/HCO+ intensity ratio in UV-illuminated regions. Finer distinctions in density classes (nH~7E3, and 4E4 cm-3) for the densest regions are also identified, likely related to the higher critical density of the CN and HCO+ (1-0) lines. The association of simultaneous multi-line, wide-field mapping and powerful machine learning methods such as the Meanshift algorithm reveals how to decode the complex information available in molecular tracers.

We analyse the circular polarisation data accumulated in the first 7 years of the POLAMI project introduced in an accompanying paper (Agudo et al.). In the 3mm wavelength band, we acquired more than 2600 observations, and all but one of our 37 sample sources were detected, most of them several times. For most sources, the observed distribution of the degree of circular polarisation is broader than that of unpolarised calibrators, indicating that weak (<0.5%) circular polarisation is present most of the time. Our detection rate and the maximum degree of polarisation found, 2.0%, are comparable to previous surveys, all made at much longer wavelengths. We argue that the process generating circular polarisation must not be strongly wavelength dependent, and we propose that the widespread presence of circular polarisation in our short wavelength sample dominated by blazars is mostly due to Faraday conversion of the linearly polarised synchrotron radiation in the helical magnetic field of the jet. Circular polarisation is variable, most notably on time scales comparable to or shorter than our median sampling interval <1 month. Longer time scales of about one year are occasionally detected, but severely limited by the weakness of the signal. At variance with some longer wavelength investigations we find that the sign of circular polarisation changes in most sources, while only 7 sources, including 3 already known, have a strong preference for one sign. The degrees of circular and linear polarisation do not show any systematic correlation. We do find however one particular event where the two polarisation degrees vary in synchronism during a time span of 0.9 years. The paper also describes a novel method for calibrating the sign of circular polarisation observations.

The development of precision cosmology with clusters of galaxies requires high-angular resolution Sunyaev-Zel'dovich (SZ) observations. As for now, arcmin resolution SZ observations (e.g. SPT, ACT and Planck) only allowed detailed studies of the intra cluster medium for low redshift clusters (z<0.2). With both a wide field of view (6.5 arcmin) and a high angular resolution (17.7 and 11.2 arcsec at 150 and 260 GHz), the NIKA2 camera installed at the IRAM 30-m telescope (Pico Veleta, Spain), will bring valuable information in the field of SZ imaging of clusters of galaxies. The NIKA2 SZ observation program will allow us to observe a large sample of clusters (50) at redshifts between 0.4 and 0.9. As a pilot study for NIKA2, several clusters of galaxies have been observed with the pathfinder, NIKA, at the IRAM 30-m telescope to cover the various configurations and observation conditions expected for NIKA2.

We propose a method to map the temperature distribution of the hot gas in galaxy clusters that uses resolved images of the thermal Sunyaev-Zel'dovich (tSZ) effect in combination with X-ray data. Application to images from the New IRAM KIDs Array (NIKA) and XMM-Newton allows us to measure and determine the spatial distribution of the gas temperature in the merging cluster MACS J0717.5+3745, at $z=0.55$. Despite the complexity of the target object, we find a good morphological agreement between the temperature maps derived from X-ray spectroscopy only -- using XMM-Newton ($T_{\rm XMM}$) and Chandra ($T_{\rm CXO}$) -- and the new gas-mass-weighted tSZ+X-ray imaging method ($T_{\rm SZ}$). We correlate the temperatures from tSZ+X-ray imaging and those from X-ray spectroscopy alone and find that $T_{\rm SZ}$ is higher than $T_{\rm XMM}$ and lower than $T_{\rm CXO}$ by $\sim 10\%$ in both cases. Our results are limited by uncertainties in the geometry of the cluster gas, contamination from kinetic SZ ($\sim 10\%$), and the absolute calibration of the tSZ map ($7\%$). Investigation using a larger sample of clusters would help minimise these effects.

Context: In the past decade, sensitive, resolved Sunyaev-Zel'dovich (SZ) studies of galaxy clusters have become common. Whereas many previous SZ studies have parameterized the pressure profiles of galaxy clusters, non-parametric reconstructions will provide insights into the thermodynamic state of the intracluster medium (ICM). Aims: We seek to recover the non-parametric pressure profiles of the high redshift ($z=0.89$) galaxy cluster CLJ 1226.9+3332 as inferred from SZ data from the MUSTANG, NIKA, Bolocam, and Planck instruments, which all probe different angular scales. Methods: Our non-parametric algorithm makes use of logarithmic interpolation, which under the assumption of ellipsoidal symmetry is analytically integrable. For MUSTANG, NIKA, and Bolocam we derive a non-parametric pressure profile independently and find good agreement among the instruments. In particular, we find that the non-parametric profiles are consistent with a fitted gNFW profile. Given the ability of Planck to constrain the total signal, we include a prior on the integrated Compton Y parameter as determined by Planck. Results: For a given instrument, constraints on the pressure profile diminish rapidly beyond the field of view. The overlap in spatial scales probed by these four datasets is therefore critical in checking for consistency between instruments. By using multiple instruments, our analysis of CLJ 1226.9+3332 covers a large radial range, from the central regions to the cluster outskirts: $0.05 R_{500} < r < 1.1 R_{500}$. This is a wider range of spatial scales than is typical recovered by SZ instruments. Similar analyses will be possible with the new generation of SZ instruments such as NIKA2 and MUSTANG2.

The characterization of dust properties in the interstellar medium (ISM) is key for star formation. Mass estimates are crucial to determine gravitational collapse conditions for the birth of new stellar objects in molecular clouds. However, most of these estimates rely on dust models that need further observational constraints from clouds to prestellar and protostellar cores. We present results of a study of dust emissivity changes based on mm-continuum data obtained with the NIKA camera at the IRAM-30m telescope. Observing dust emission at 1.15 mm and 2 mm allows us to constrain the dust emissivity index ($\beta$) in the Rayleigh-Jeans tail of the dust spectral energy distribution (SED) far from its peak emission, where the contribution of other parameters (i.e. dust temperature) is important. Focusing on the Taurus molecular cloud, a low-mass star-forming regions in the Gould Belt, we analyze the emission properties of several distinct objects in the B213 filament: three prestellar cores, two Class-0/I protostellar cores and one Class-II object. By means of the ratio of the two NIKA channel-maps, we show that in the Rayleigh-Jeans approximation the dust emissivity index varies among the objects. For one prestellar and two protostellar cores, we produce a robust study using Herschel data to constrain the dust temperature of the sources. By using the Abel transform inversion technique we get accurate radial $\beta$ profiles. We find systematic spatial variations of $\beta$ in the protostellar cores that is not observed in the prestellar core. While in the former case $\beta$ decreases toward the center, in the latter it remains constant. Moreover, $\beta$ appears anticorrelated with the dust temperature. We discuss the implication of these results in terms of dust grain evolution between pre- and protostellar cores.

We present \ci\,(2--1) and multi-transition $^{12}$CO observations of a dusty star-forming galaxy, ACT\,J2029+0120, which we spectroscopically confirm to lie at $z$\,=\,2.64. We detect CO(3--2), CO(5--4), CO(7--6), CO(8--7), and \ci\,(2--1) at high significance, tentatively detect HCO$^{+}$(4--3), and place strong upper limits on the integrated strength of dense gas tracers (HCN(4--3) and CS(7--6)). Multi-transition CO observations and dense gas tracers can provide valuable constraints on the molecular gas content and excitation conditions in high-redshift galaxies. We therefore use this unique data set to construct a CO spectral line energy distribution (SLED) of the source, which is most consistent with that of a ULIRG/Seyfert or QSO host object in the taxonomy of the \textit{Herschel} Comprehensive ULIRG Emission Survey. We employ RADEX models to fit the peak of the CO SLED, inferring a temperature of T$\sim$117 K and $n_{\text{H}_2}\sim10^5$ cm$^{-3}$, most consistent with a ULIRG/QSO object and the presence of high density tracers. We also find that the velocity width of the \ci\ line is potentially larger than seen in all CO transitions for this object, and that the $L'_{\rm C\,I(2-1)}/L'_{\rm CO(3-2)}$ ratio is also larger than seen in other lensed and unlensed submillimeter galaxies and QSO hosts; if confirmed, this anomaly could be an effect of differential lensing of a shocked molecular outflow.

Magnetic fields, which play a major role in a large number of astrophysical processes from galactic to cosmological scales, can be traced via observations of dust polarization as demonstrated by the Planck satellite results. In particular, low-resolution observations of dust polarization have demonstrated that Galactic filamentary structures, where star formation takes place, are associated to well organized magnetic fields. A better understanding of this process requires detailed observations of galactic dust polarization on scales of 0.01 to 0.1 pc. Such high-resolution polarization observations can be carried out at the IRAM 30 m telescope using the recently installed NIKA2 camera, which features two frequency bands at 260 and 150 GHz (respectively 1.15 and 2.05 mm), the 260 GHz band being polarization sensitive. NIKA2 so far in commissioning phase, has its focal plane filled with ~3300 detectors to cover a Field of View (FoV) of 6.5 arcminutes diameter. The NIKA camera, which consisted of two arrays of 132 and 224 Lumped Element Kinetic Inductance Detectors (LEKIDs) and a FWHM (Full-Width-Half-Maximum) of 12 and 18.2 arcsecond at 1.15 and 2.05 mm respectively, has been operated at the IRAM 30 m telescope from 2012 to 2015 as a test-bench for NIKA2. NIKA was equipped of a room temperature polarization system (a half wave plate (HWP) and a grid polarizer facing the NIKA cryostat window). The fast and continuous rotation of the HWP permits the quasi simultaneous reconstruction of the three Stokes parameters, I, Q and U at 150 and 260 GHz. This paper presents the first polarization measurements with KIDs and reports the polarization performance of the NIKA camera and the pertinence of the choice of the polarization setup in the perspective of NIKA2. (abridged)

Context. The combination of wideband receivers and spectrometers currently available in (sub-)millimeter observatories deliver wide- field hyperspectral imaging of the interstellar medium. Tens of spectral lines can be observed over degree wide fields in about fifty hours. This wealth of data calls for restating the physical questions about the interstellar medium in statistical terms. Aims. We aim at gaining information on the physical structure of the interstellar medium from a statistical analysis of many lines from different species over a large field of view, without requiring detailed radiative transfer or astrochemical modeling. Methods. We coupled a nonlinear rescaling of the data with one of the simplest multivariate analysis methods, namely the Principal Component Analysis, to decompose the observed signal into components that we interpret first qualitatively and then quantitatively based on our deep knowledge of the observed region and of the astrochemistry at play. Results. We identify 3 principal components, linear compositions of line brightness temperatures, that are correlated at various levels with the column density, the volume density and the UV radiation field. Conclusions. When sampling a sufficiently diverse mixture of physical parameters, it is possible to decompose the molecular emission in order to gain physical insight on the observed interstellar medium. This opens a new avenue for future studies of the interstellar medium.

Recent population studies have shown that the variability Doppler factors can adequately describe blazars as a population. We use the flux density variations found within the extensive radio multi-wavelength datasets of the F-GAMMA program, a total of 10 frequencies from 2.64 up to 142.33 GHz, in order to estimate the variability Doppler factors for 58 $\gamma$-ray bright sources, for 20 of which no variability Doppler factor has been estimated before. We employ specifically designed algorithms in order to obtain a model for each flare at each frequency. We then identify each event and track its evolution through all the available frequencies for each source. This approach allows us to distinguish significant events producing flares from stochastic variability in blazar jets. It also allows us to effectively constrain the variability brightness temperature and hence the variability Doppler factor as well as provide error estimates. Our method can produce the most accurate (16\% error on average) estimates in the literature to date.

The nature of turbulence in molecular clouds is one of the key parameters that control star formation efficiency: compressive motions, as opposed to solenoidal motions, can trigger the collapse of cores, or mark the expansion of Hii regions. We try to observationally derive the fractions of momentum density ($\rho v$) contained in the solenoidal and compressive modes of turbulence in the Orion B molecular cloud and relate these fractions to the star formation efficiency in the cloud. The implementation of a statistical method developed by Brunt & Federrath (2014), applied to a $^{13}$CO(J=1-0) datacube obtained with the IRAM-30m telescope, allows us to retrieve 3-dimensional quantities from the projected quantities provided by the observations, yielding an estimate of the compressive versus solenoidal ratio in various regions of the cloud. Despite the Orion B molecular cloud being highly supersonic (mean Mach number $\sim$ 6), the fractions of motion in each mode diverge significantly from equipartition. The cloud's motions are on average mostly solenoidal (excess > 8 % with respect to equipartition), which is consistent with its low star formation rate. On the other hand, the motions around the main star-forming regions (NGC 2023 and NGC 2024) prove to be strongly compressive. We have successfully applied to observational data a method that was so far only tested on simulations, and have shown that there can be a strong intra-cloud variability of the compressive and solenoidal fractions, these fractions being in turn related to the star formation efficiency. This opens a new possibility for star-formation diagnostics in galactic molecular clouds.

Measurement of the gas velocity distribution in galaxy clusters provides insight into the physics of mergers, through which large scale structures form in the Universe. Velocity estimates within the intracluster medium (ICM) can be obtained via the Sunyaev-Zel'dovich (SZ) effect, but its observation is challenging both in term of sensitivity requirement and control of systematic effects, including the removal of contaminants. In this paper we report resolved observations, at 150 and 260 GHz, of the SZ effect toward the triple merger MACS J0717.5+3745 (z=0.55), using data obtained with the NIKA camera at the IRAM 30m telescope. Assuming that the SZ signal is the sum of a thermal (tSZ) and a kinetic (kSZ) component and by combining the two NIKA bands, we extract for the first time a resolved map of the kSZ signal in a cluster. The kSZ signal is dominated by a dipolar structure that peaks at -5.1 and +3.4 sigma, corresponding to two subclusters moving respectively away and toward us and coincident with the cold dense X-ray core and a hot region undergoing a major merging event. We model the gas electron density and line-of-sight velocity of MACS J0717.5+3745 as four subclusters. Combining NIKA data with X-ray observations from XMM-Newton and Chandra, we fit this model to constrain the gas line-of-sight velocity of each component, and we also derive, for the first time, a velocity map from kSZ data (i.e. that is model-dependent). Our results are consistent with previous constraints on the merger velocities, and thanks to the high angular resolution of our data, we are able to resolve the structure of the gas velocity. Finally, we investigate possible contamination and systematic effects with a special care given to radio and submillimeter galaxies. Among the sources that we detect with NIKA, we find one which is likely to be a high redshift lensed submillimeter galaxy.

We aim to develop the Orion B Giant Molecular Cloud (GMC) as a local template for interpreting extra-galactic molecular line observations. We use the wide-band receiver at the IRAM-30m to spatially and spectrally resolve the Orion B GMC. The observations cover almost 1 square degree at 26" resolution with a bandwidth of 32 GHz from 84 to 116 GHz in only two tunings. Among the mapped spectral lines are the 12CO, 13CO, C18O, C17O, HCN, HNC, 12CN, CCH, HCO+, N2H+ (1-0), and 12CS, 32SO, SiO, c-C3H2, CH3OH (2-1) transitions. We introduce the molecular anatomy of the Orion B GMC, including relations between line intensities and gas column density or far-UV radiation fields, and correlations between selected line and line ratios. We also obtain a dust-traced gas mass that is less than about one third the CO-traced mass, using the standard Xco conversion factor. The presence of overluminous CO can be traced back to the dependence of the CO intensity on UV illumination. In fact, while most lines show some dependence on the UV radiation field, CN and CCH are the most sensitive. Moreover dense cloud cores are almost exclusively traced by N2H+. Other traditional high density tracers, such as HCN (1-0), are also easily detected in extended translucent regions at a typical density of about 500 H2 cm-3. In general, we find no straightforward relation between line critical density and the fraction of the line luminosity coming from dense gas regions. Our initial findings demonstrate that the relations between line (ratio) intensities and environment in GMCs are more complicated than often assumed. Sensitivity (i.e., the molecular column density), excitation, and above all chemistry contribute to the observed line intensity distributions. They must be considered together when developing the next generation of extra-galactic molecular line diagnostics of mass, density, temperature and radiation field.

The determination of the thermodynamic properties of clusters of galaxies at intermediate and high redshift can bring new insights into the formation of large-scale structures. It is essential for a robust calibration of the mass-observable scaling relations and their scatter, which are key ingredients for precise cosmology using cluster statistics. Here we illustrate an application of high resolution $(< 20$ arcsec) thermal Sunyaev-Zel'dovich (tSZ) observations by probing the intracluster medium (ICM) of the \planck-discovered galaxy cluster \psz\ at redshift $z = 0.61$, using tSZ data obtained with the NIKA camera, which is a dual-band (150 and 260~GHz) instrument operated at the IRAM 30-meter telescope. We deproject jointly NIKA and \planck\ data to extract the electronic pressure distribution from the cluster core ($R \sim 0.02\, R_{500}$) to its outskirts ($R \sim 3\, R_{500}$) non-parametrically for the first time at intermediate redshift. The constraints on the resulting pressure profile allow us to reduce the relative uncertainty on the integrated Compton parameter by a factor of two compared to the \planck\ value. Combining the tSZ data and the deprojected electronic density profile from \xmm\ allows us to undertake a hydrostatic mass analysis, for which we study the impact of a spherical model assumption on the total mass estimate. We also investigate the radial temperature and entropy distributions. These data indicate that \psz\ is a massive ($M_{500} \sim 5.5 \times 10^{14}$ M$_{\odot}$) cool-core cluster. This work is part of a pilot study aiming at optimizing the treatment of the NIKA2 tSZ large program dedicated to the follow-up of SZ-discovered clusters at intermediate and high redshifts. (abridged)

To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region.

The Gamma-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features including an apparent (~100 deg) position angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps together with cm/mm total intensity and Gamma-ray observations from Fermi/LAT from 2008-2014, we aimed to determine the location of Gamma-ray emission and to explain the inner-mas structural changes. Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and allow us to observe above the synchrotron self-absorption peak frequency. The jet was spectrally decomposed at multiple locations along the jet. From this we derived estimates of the magnetic field. How the field decreases down the jet allowed an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics we attempt to locate the site of Gamma-ray activity, radio flares and spectral changes. Strong Gamma-ray flares appeared to originate from either the "core" region, a downstream stationary feature, or both, with Gamma-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be >1.6 G in the "core" and >0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the "core" as >6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.

The New IRAM KID Array 2 (NIKA 2) is a dual-band camera operating with three frequency-multiplexed kilopixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) cooled at 150 mK. NIKA 2 is designed to observe the intensity and polarisation of the sky at 1.15 and 2.0 mm wavelength from the IRAM 30 m telescope. The NIKA 2 instrument represents a huge step in performance as compared to the NIKA pathfinder instrument, which has already shown state-of-the-art detector and photometric performance. After the commissioning planned to be accomplished at the end of 2016, NIKA 2 will be an IRAM resident instrument for the next ten years or more. NIKA 2 should allow the astrophysical community to tackle a large number of open questions reaching from the role of the Galactic magnetic field in star formation to the discrepancy between cluster-based and CMB-based cosmology possibly induced by the unknown cluster physics. We present an overview of the commissioning phase together with some first results.

NIKA2 (New IRAM KID Arrays) is a dual band (150 and 260 GHz) imaging camera based on Kinetic Inductance Detectors (KIDs) and designed to work at the IRAM 30 m telescope (Pico Veleta, Spain). Built on the experience of the NIKA prototype, NIKA2 has been installed at the 30 m focal plane in October 2015 and the commissioning phase is now ongoing. Through the thermal Sunyaev-Zeldovich (tSZ) effect, NIKA2 will image the ionized gas residing in clusters of galaxies with a resolution of 12 and 18 arcsec FWHM (at 150 and 260 GHz, respectively). We report on the recent tSZ measurements with the NIKA camera and discuss the future objectives for the NIKA2 SZ large Program, 300h of observation dedicated to SZ science. With this program we intend to perform a high angular resolution follow-up of a cosmologically-representative sample of clusters belonging to SZ catalogues, with redshift greater than 0.5. The main output of the program will be the study of the redshift evolution of the cluster pressure profile as well as that of the scaling laws relating the cluster global properties.

Detailed investigation of broadband flux variability in the blazar 3C 273 allows us to probe the location and size of emission regions and their physical conditions. We report the results on correlation studies of the flaring activity observed between 2008 and 2012. The observed broadband variations were investigated using the structure function and the discrete correlation function, and power spectral density analysis (PSD) methods. The PSD analysis showed that the optical/IR light curve slopes are consistent with the slope of white noise processes, while, the PSD slopes at radio, X-ray and gamma-ray energies are consistent with red-noise processes. The flux variations at gamma-ray and mm-radio bands are found to be significantly correlated. Using the estimated time lag of (110\pm27) days between gamma-ray and radio light curves, we constrained the location of the gamma-ray emission region at a de-projected distance of 1.2\pm0.9 pc from the jet apex. Flux variations at X-ray bands were found to have a significant correlation with variations at both radio and \gamma-rays energies. The correlation between X-rays and gamma-rays light curves suggests presence of two components responsible for the X-ray emission. A negative time lag of -(50\pm20) days, where the X-rays are leading the emission, suggests X-rays are emitted closer to the jet apex from a compact region at a distance of ~(0.5\pm0.4) pc from the jet apex. A positive time lag of (110\pm20) days suggests jet-base origin of the other X-ray component at ~(4--5)~pc from the jet apex. The flux variations at radio frequencies were found to be well correlated with each other such that the variations at higher frequencies are leading the lower frequencies, which could be expected in the standard shock-in-jet model.

Context. The origin of blazar variability, seen from radio up to gamma rays, is still a heavily debated matter and broadband flares offer a unique testbed towards a better understanding of these extreme objects. Such an energetic outburst was detected by Fermi/LAT in 2008 from the blazar PKS 1502+106. The outburst was observed from gamma rays down to radio frequencies. Aims. Through the delay between flare maxima at different radio frequencies, we study the frequency-dependent position of the unit-opacity surface and infer its absolute position with respect to the jet base. This nuclear opacity profile enables the magnetic field tomography of the jet. We also localize the gamma-ray emission region and explore the mechanism producing the flare. Methods. The radio flare of PKS 1502+106 is studied through single-dish flux density measurements at 12 frequencies in the range 2.64 to 226.5 GHz. To quantify it, we employ both a Gaussian process regression and a discrete cross-correlation function analysis. Results. We find that the light curve parameters (flare amplitude and cross-band delays) show a power-law dependence on frequency. Delays decrease with frequency, and the flare amplitudes increase up to about 43 GHz and then decay. This behavior is consistent with a shock propagating downstream the jet. The self-absorbed radio cores are located between about 10 and 4 pc from the jet base and their magnetic field strengths range between 14 and 176 mG, at the frequencies 2.64 to 86.24 GHz. Finally, the gamma-ray active region is located at (1.9 +/- 1.1) pc away from the jet base.

We present the outline and first results of a project using the synergies of the long term blazar radiomillimetre monitoring program F-GAMMA, the continued scanning of the millimetre-submillimetre sky by the Planck satellite, together with several dedicated observing programs at the Effelsberg 100m telescope, to obtain a data sample unprecedented in both time resolution and frequency span.

We report about recent updates of broad-band properties of radio-loud narrow-line Seyfert 1 galaxies.

The thermal Sunyaev-Zeldovich effect (tSZ) is a powerful probe to study clusters of galaxies and is complementary with respect to X-ray, lensing or optical observations. Previous arcmin resolution tSZ observations ({\it e.g.} SPT, ACT and Planck) only enabled detailed studies of the intra-cluster medium morphology for low redshift clusters ($z < 0.2$). Thus, the development of precision cosmology with clusters requires high angular resolution observations to extend the understanding of galaxy cluster towards high redshift. NIKA2 is a wide-field (6.5 arcmin field of view) dual-band camera, operated at $100 \ {\rm mK}$ and containing $\sim 3300$ KID (Kinetic Inductance Detectors), designed to observe the millimeter sky at 150 and 260 GHz, with an angular resolution of 18 and 12 arcsec respectively. The NIKA2 camera has been installed on the IRAM 30-m telescope (Pico Veleta, Spain) in September 2015. The NIKA2 tSZ observation program will allow us to observe a large sample of clusters (50) at redshift ranging between 0.5 and 1. As a pathfinder for NIKA2, several clusters of galaxies have been observed at the IRAM 30-m telescope with the NIKA prototype to cover the various configurations and observation conditions expected for NIKA2.

The prototype of the NIKA2 camera, NIKA, is an instrument operating at the IRAM 30-m telescope, which can observe simultaneously at 150 and 260GHz. One of the main goals of NIKA2 is to measure the pressure distribution in galaxy clusters at high resolution using the thermal SZ (tSZ) effect. Such observations have already proved to be an excellent probe of cluster pressure distributions even at high redshifts. However, an important fraction of clusters host submm and/or radio point sources, which can significantly affect the reconstructed signal. Here we report on <20" resolution observations at 150 and 260GHz of the cluster MACSJ1424, which hosts both radio and submm point sources. We examine the morphology of the tSZ signal and compare it to other datasets. The NIKA data are combined with Herschel satellite data to study the SED of the submm point source contaminants. We then perform a joint reconstruction of the intracluster medium (ICM) electronic pressure and density by combining NIKA, Planck, XMM-Newton, and Chandra data, focusing on the impact of the radio and submm sources on the reconstructed pressure profile. We find that large-scale pressure distribution is unaffected by the point sources because of the resolved nature of the NIKA observations. The reconstructed pressure in the inner region is slightly higher when the contribution of point sources are removed. We show that it is not possible to set strong constraints on the central pressure distribution without accurately removing these contaminants. The comparison with X-ray only data shows good agreement for the pressure, temperature, and entropy profiles, which all indicate that MACSJ1424 is a dynamically relaxed cool core system. The present observations illustrate the possibility of measuring these quantities with a relatively small integration time, even at high redshift and without X-ray spectroscopy.

NIKA 2 (New Instrument of Kids Array) is a next generation continuum and polarized instrument successfully installed in October 2015 at the IRAM 30 m telescope on Pico-Veleta (Granada, Spain). NIKA 2 is a high resolution dual-band camera, operating with frequency multiplexed LEKIDs (Lumped Element Kinetic Inductance Detectors) cooled at 100 mK. Dual color images are obtained thanks to the simultaneous readout of a 1020 pixels array at 2 mm and 1140 x 2 pixels arrays at 1.15 mm with a final resolution of 18 and 12 arcsec respectively, and 6.5 arcmin of Field of View (FoV). The two arrays at 1.15 mm allow us to measure the linear polarization of the incoming light. This will place NIKA 2 as an instrument of choice to study the role of magnetic fields in the star formation process. The NIKA experiment, a prototype for NIKA 2 with a reduced number of detectors (about 400 LEKIDs) and FoV (1.8 arcmin), has been successfully operated at the IRAM 30 telescope in several open observational campaigns. The performance of the NIKA 2 polarization setup has been successfully validated with the NIKA prototype.

NIKA2 (New IRAM KID Array 2) is a camera dedicated to millimeter wave astronomy based upon kilopixel arrays of Kinetic Inductance Detectors (KID). The pathfinder instrument, NIKA, has already shown state-of-the-art detector performance. NIKA2 builds upon this experience but goes one step further, increasing the total pixel count by a factor $\sim$10 while maintaining the same per pixel performance. For the next decade, this camera will be the resident photometric instrument of the Institut de Radio Astronomie Millimetrique (IRAM) 30m telescope in Sierra Nevada (Spain). In this paper we give an overview of the main components of NIKA2, and describe the achieved detector performance. The camera has been permanently installed at the IRAM 30m telescope in October 2015. It will be made accessible to the scientific community at the end of 2016, after a one-year commissioning period. When this happens, NIKA2 will become a fundamental tool for astronomers worldwide.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

We report on a multiwavelength survey of a sample of 42 flat-spectrum radio-loud narrow-line Seyfert 1 galaxies (RLNLS1s). This is the largest known sample of this type of active galactic nucleus (AGN) to date. We found that 17% of sources were detected at high-energy gamma rays (E>100 MeV), and 90% at X-rays (0.3-10 keV). The masses of the central black holes are in the range $\sim 10^{6-8}M_{\odot}$, smaller than the values of blazars. The disk luminosities are about 1-49% of the Eddington value, with one outlier at 0.3%, comparable with the luminosities observed in flat-spectrum radio quasars (FSRQs). The jet powers are $\sim 10^{42-46}$ erg s$^{-1}$, comparable with BL Lac Objects, yet relatively smaller than FSRQs. However, once renormalized by the mass of the central black hole, the jet powers of RLNLS1s, BL Lacs, and FSRQs are consistent each other, indicating the scalability of the jets. We found episodes of extreme variability at high energies on time scales of hours. In some cases, dramatic spectral and flux changes are interpreted as the interplay between the relativistic jet and the accretion disk. We conclude that, despite the distinct observational properties, the central engines of RLNLS1s are similar to those of blazars.

Context. Blazars are among the most energetic objects in the Universe. In 2008 August, Fermi/LAT detected the blazar PKS 1502+106 showing a rapid and strong gamma-ray outburst followed by high and variable flux over the next months. This activity at high energies triggered an intensive multi-wavelength campaign covering also the radio, optical, UV, and X-ray bands indicating that the flare was accompanied by a simultaneous outburst at optical/UV/X-rays and a delayed outburst at radio bands. Aims: In the current work we explore the phenomenology and physical conditions within the ultra-relativistic jet of the gamma-ray blazar PKS 1502+106. Additionally, we address the question of the spatial localization of the MeV/GeV-emitting region of the source. Methods: We utilize ultra-high angular resolution mm-VLBI observations at 43 and 86 GHz complemented by VLBI observations at 15 GHz. We also employ single-dish radio data from the F-GAMMA program at frequencies matching the VLBI monitoring. Results: PKS 1502+106 shows a compact core-jet morphology and fast superluminal motion with apparent speeds in the range 5--22 c. Estimation of Doppler factors along the jet yield values between ~7 up to ~50. This Doppler factor gradient implies an accelerating jet. The viewing angle towards the source differs between the inner and outer jet, with the former at ~3 degrees and the latter at ~1 degree, after the jet bends towards the observer beyond 1 mas. The de-projected opening angle of the ultra-fast, magnetically-dominated jet is found to be (3.8 +/- 0.5) degrees. A single jet component can be associated with the pronounced flare both at high-energies and in radio bands. Finally, the gamma-ray emission region is localized at less than 5.9 pc away from the jet base.

NIKA is a dual-band camera operating with 315 frequency multiplexed LEKIDs cooled at 100 mK. NIKA is designed to observe the sky in intensity and polarisation at 150 and 260 GHz from the IRAM 30-m telescope. It is a test-bench for the final NIKA2 camera. The incoming linear polarisation is modulated at four times the mechanical rotation frequency by a warm rotating multi-layer Half Wave Plate. Then, the signal is analysed by a wire grid and finally absorbed by the LEKIDs. The small time constant (< 1ms ) of the LEKID detectors combined with the modulation of the HWP enables the quasi-simultaneous measurement of the three Stokes parameters I, Q, U, representing linear polarisation. In this paper we present results of recent observational campaigns demonstrating the good performance of NIKA in detecting polarisation at mm wavelength.

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Mets\"ahovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

The thermal Sunyaev-Zel'dovich (tSZ) effect is expected to provide a low scatter mass proxy for galaxy clusters since it is directly proportional to the cluster thermal energy. The tSZ observations have proven to be a powerful tool for detecting and studying them, but high angular resolution observations are now needed to push their investigation to a higher redshift. In this paper, we report high angular (< 20 arcsec) resolution tSZ observations of the high-redshift cluster CL J1226.9+3332 (z=0.89). It was imaged at 150 and 260 GHz using the NIKA camera at the IRAM 30-meter telescope. The 150 GHz map shows that CL J1226.9+3332 is morphologically relaxed on large scales with evidence of a disturbed core, while the 260 GHz channel is used mostly to identify point source contamination. NIKA data are combined with those of Planck and X-ray from Chandra to infer the cluster's radial pressure, density, temperature, and entropy distributions. The total mass profile of the cluster is derived, and we find $M_{500} = 5.96^{+1.02}_{-0.79} $ x $10^{14} M_{\odot}$ within the radius $R_{500} = 930^{+50}_{-43}$ kpc, at a 68% confidence level. ($R_{500}$ is the radius within which the average density is 500 times the critical density at the cluster's redshift.) NIKA is the prototype camera of NIKA2, a KIDs (kinetic inductance detectors) based instrument to be installed at the end of 2015. This work is, therefore, part of a pilot study aiming at optimizing tSZ NIKA2 large programs.

Blazars are among the most variable objects in the universe. They feature energetic jets of plasma that launch from the cores of these active galactic nuclei (AGN), triggering activity from radio up to gamma-ray energies. Spatial localization of the region of their MeV/GeV emission is a key question in understanding the blazar phenomenon. The flat spectrum radio quasar (FSRQ) PKS 1502+106 has exhibited extreme and correlated, radio and high-energy activity that triggered intense monitoring by the Fermi-GST AGN Multi-frequency Monitoring Alliance (F-GAMMA) program and the Global Millimeter VLBI Array (GMVA) down to $\lambda$3 mm (or 86 GHz), enabling the sharpest view to date towards this extreme object. Here, we report on preliminary results of our study of the gamma-ray loud blazar PKS 1502+106, combining VLBI and single dish data. We deduce the critical aspect angle towards the source to be $\theta_{\rm c} = 2.6^{\circ}$, calculate the apparent and intrinsic opening angles and constrain the distance of the 86 GHz core from the base of the conical jet, directly from mm-VLBI but also through a single dish relative timing analysis. Finally, we conclude that gamma rays from PKS 1502+106 originate from a region between ~1-16 pc away from the base of the hypothesized conical jet, well beyond the bulk of broad-line region (BLR) material of the source.

We studied the radio emission from four radio-loud and gamma-ray-loud narrow-line Seyfert 1 galaxies. The goal was to investigate whether a relativistic jet is operating at the source, and quantify its characteristics. We relied on the most systematic monitoring of such system in the cm and mm radio bands which is conducted with the Effelsberg 100 m and IRAM 30 m telescopes and covers the longest time-baselines and the most radio frequencies to date. We extract variability parameters and compute variability brightness temperatures and Doppler factors. The jet powers were computed from the light curves to estimate the energy output. The dynamics of radio spectral energy distributions were examined to understand the mechanism causing the variability. All the sources display intensive variability that occurs at a pace faster than what is commonly seen in blazars. The flaring events show intensive spectral evolution indicative of shock evolution. The brightness temperatures and Doppler factors are moderate, implying a mildly relativistic jet. The computed jet powers show very energetic flows. The radio polarisation in one case clearly implies a quiescent jet underlying the recursive flaring activity. Despite the generally lower flux densities, the sources appear to show all typical characteristics seen in blazars that are powered by relativistic jets.

NIKA (New IRAM KID Arrays) is a dual-band imaging instrument installed at the IRAM (Institut de RadioAstronomie Millimetrique) 30-meter telescope at Pico Veleta (Spain). Two distinct Kinetic Inductance Detectors (KID) focal planes allow the camera to simultaneously image a field-of-view of about 2 arc-min in the bands 125 to 175 GHz (150 GHz) and 200 to 280 GHz (240 GHz). The sensitivity and stability achieved during the last commissioning Run in June 2013 allows opening the instrument to general observers. We report here the latest results, in particular in terms of sensitivity, now comparable to the state-of-the-art Transition Edge Sensors (TES) bolometers, relative and absolute photometry. We describe briefly the next generation NIKA-2 instrument, selected by IRAM to occupy, from 2015, the continuum imager/polarimeter slot at the 30-m telescope.

Clusters of galaxies provide valuable information on the evolution of the Universe and large scale structures. Recent cluster observations via the thermal Sunyaev-Zel'dovich (tSZ) effect have proven to be a powerful tool to detect and study them. In this context, high resolution tSZ observations (~ tens of arcsec) are of particular interest to probe intermediate and high redshift clusters. Observations of the tSZ effect will be carried out with the millimeter dual-band NIKA2 camera, based on Kinetic Inductance Detectors (KIDs) to be installed at the IRAM 30-meter telescope in 2015. To demonstrate the potential of such an instrument, we present tSZ observations with the NIKA camera prototype, consisting of two arrays of 132 and 224 detectors that observe at 140 and 240 GHz with a 18.5 and 12.5 arcsec angular resolution, respectively. The cluster RX J1347.5-1145 was observed simultaneously at 140 and 240 GHz. We used a spectral decorrelation technique to remove the atmospheric noise and obtain a map of the cluster at 140 GHz. The efficiency of this procedure has been characterized through realistic simulations of the observations. The observed 140 GHz map presents a decrement at the cluster position consistent with the tSZ nature of the signal. We used this map to study the pressure distribution of the cluster by fitting a gNFW model to the data. Subtracting this model from the map, we confirm that RX J1347.5-1145 is an ongoing merger, which confirms and complements previous tSZ and X-ray observations. For the first time, we demonstrate the tSZ capability of KID based instruments. The NIKA2 camera with ~ 5000 detectors and a 6.5 arcmin field of view will be well-suited for in-depth studies of the intra cluster medium in intermediate to high redshifts, which enables the characterization of recently detected clusters by the Planck satellite.

High resolution observations of the thermal Sunyaev-Zel'dovich (tSZ) effect are necessary to allow the use of clusters of galaxies as a probe for large scale structures at high redshifts. With its high resolution and dual-band capability at millimeter wavelengths, the NIKA camera can play a significant role in this context. NIKA is based on newly developed Kinetic Inductance Detectors (KIDs) and operates at the IRAM 30m telescope, Pico Veleta, Spain. In this paper, we give the status of the NIKA camera, focussing on the KID technology. We then present observations of three galaxy clusters: RX J1347.5-1145 as a demonstrator of the NIKA capabilities and the recent observations of CL J1226.9+3332 (z = 0.89) and MACS J0717.5+3745 (z = 0.55). We also discuss prospects for the final NIKA2 camera, which will have a 6.5 arcminute field of view with about 5000 detectors in two bands at 150 and 260 GHz.

The New IRAM KID Array (NIKA) is a dual-band camera operating with frequency multiplexed arrays of Lumped Element Kinetic Inductance Detectors (LEKIDs) cooled to 100 mK. NIKA is designed to observe the intensity and polarisation of the sky at 1.25 and 2.14 mm from the IRAM 30 m telescope. We present the improvements on the control of systematic effects and astrophysical results made during the last observation campaigns between 2012 and 2014.

The New IRAM KID Array (NIKA) instrument is a dual-band imaging camera operating with Kinetic Inductance Detectors (KID) cooled at 100 mK. NIKA is designed to observe the sky at wavelengths of 1.25 and 2.14 mm from the IRAM 30 m telescope at Pico Veleta with an estimated resolution of 13\,arcsec and 18 arcsec, respectively. This work presents the performance of the NIKA camera prior to its opening to the astrophysical community as an IRAM common-user facility in early 2014. NIKA is a test bench for the final NIKA2 instrument to be installed at the end of 2015. The last NIKA observation campaigns on November 2012 and June 2013 have been used to evaluate this performance and to improve the control of systematic effects. We discuss here the dynamical tuning of the readout electronics to optimize the KID working point with respect to background changes and the new technique of atmospheric absorption correction. These modifications significantly improve the overall linearity, sensitivity, and absolute calibration performance of NIKA. This is proved on observations of point-like sources for which we obtain a best sensitivity (averaged over all valid detectors) of 40 and 14 mJy.s$^{1/2}$ for optimal weather conditions for the 1.25 and 2.14 mm arrays, respectively. NIKA observations of well known extended sources (DR21 complex and the Horsehead nebula) are presented. This performance makes the NIKA camera a competitive astrophysical instrument.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

The exact location of the gamma-ray emitting region in blazars is still controversial. In order to attack this problem we present first results of a cross-correlation analysis between radio (11 cm to 0.8 mm wavelength, F-GAMMA program) and gamma-ray (0.1-300 GeV) ~ 3.5 year light curves of 54 Fermi-bright blazars. We perform a source stacking analysis and estimate significances and chance correlations using mixed source correlations. Our results reveal: (i) the first highly significant multi-band radio and gamma-ray correlations (radio lagging gamma rays) when averaging over the whole sample, (ii) average time delays (source frame: 76+/-23 to 7+/-9 days), systematically decreasing from cm to mm/sub-mm bands with a frequency dependence tau_r,gamma (nu) ~ nu^-1, in good agreement with jet opacity dominated by synchrotron self-absorption, (iii) a bulk gamma-ray production region typically located within/upstream of the 3 mm core region (tau_3mm,gamma=12+/-8 days), (iv) mean distances between the region of gamma-ray peak emission and the radio "tau=1 photosphere" decreasing from 9.8+/-3.0 pc (11 cm) to 0.9+/-1.1 pc (2 mm) and 1.4+/-0.8 pc (0.8 mm), (v) 3 mm/gamma-ray correlations in 9 individual sources at a significance level where one is expected by chance (probability: 4 times 10^-6), (vi) opacity and "time lag core shift" estimates for quasar 3C 454.3 providing a lower limit for the distance of the bulk gamma-ray production region from the supermassive black hole (SMBH) of ~ 0.8-1.6 pc, i.e. at the outer edge of the Broad Line Region (BLR) or beyond. A 3 mm tau=1 surface at ~ 2-3 pc from the jet-base (i.e. well outside the "canonical BLR") finally suggests that BLR material extends to several pc distances from the SMBH.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

The New IRAM KIDs Array (NIKA) is a pathfinder instrument devoted to millimetric astronomy. In 2009 it was the first multiplexed KID camera on the sky; currently it is installed at the focal plane of the IRAM 30-meters telescope at Pico Veleta (Spain). We present preliminary data from the last observational run and the ongoing developments devoted to the next NIKA-2 kilopixels camera, to be commissioned in 2015. We also report on the latest laboratory measurements, and recent improvements in detector cosmetics and read-out electronics. Furthermore, we describe a new acquisition strategy allowing us to improve the photometric accuracy, and the related automatic tuning procedure.

We present the first detection of the thermal Sunyaev-Zel'dovich (tSZ) effect from a cluster of galaxies performed with a KIDs (Kinetic Inductance Detectors) based instrument. The tSZ effect is a distortion of the black body CMB (Cosmic Microwave Background) spectrum produced by the inverse Compton interaction of CMB photons with the hot electrons of the ionized intra-cluster medium. The massive, intermediate redshift cluster RX J1347.5-1145 has been observed using NIKA (New IRAM KIDs arrays), a dual-band (140 and 240 GHz) mm-wave imaging camera, which exploits two arrays of hundreds of KIDs: the resonant frequencies of the superconducting resonators are shifted by mm-wave photons absorption. This tSZ cluster observation demonstrates the potential of the next generation NIKA2 instrument, being developed for the 30m telescope of IRAM, at Pico Veleta (Spain). NIKA2 will have 1000 detectors at 140GHz and 2x2000 detectors at 240GHz, providing in that band also a measurement of the linear polarization. NIKA2 will be commissioned in 2015.

[Abridged] Context. To construct and interpret the spectral energy distribution (SED) of BL Lacertae objects, simultaneous broad-band observations are mandatory. Aims. We present the results of a dedicated multi-wavelength study of the high-frequency peaked BL Lacertae (HBL) object and known TeV emitter 1ES 2344+514 by means of a pre-organised campaign. Methods. The observations were conducted during simultaneous visibility windows of MAGIC and AGILE in late 2008. The measurements were complemented by Mets\"ahovi, RATAN-600, KVA+Tuorla, Swift and VLBA pointings. Additional coverage was provided by the ongoing long-term F-GAMMA and MOJAVE programs, the OVRO 40-m and CrAO telescopes as well as the Fermi satellite. The obtained SEDs are modelled using a one-zone as well as a self-consistent two-zone synchrotron self-Compton model. Results. 1ES 2344+514 was found at very low flux states in both X-rays and very high energy gamma rays. Variability was detected in the low frequency radio and X-ray bands only, where for the latter a small flare was observed. The X-ray flare was possibly caused by shock acceleration characterised by similar cooling and acceleration time scales. MOJAVE VLBA monitoring reveals a static jet whose components are stable over time scales of eleven years, contrary to previous findings. There appears to be no significant correlation between the 15 GHz and R-band monitoring light curves. The observations presented here constitute the first multi-wavelength campaign on 1ES 2344+514 from radio to VHE energies and one of the few simultaneous SEDs during low activity states. The quasi-simultaneous Fermi-LAT data poses some challenges for SED modelling, but in general the SEDs are described well by both applied models. The resulting parameters are typical for TeV emitting HBLs. Consequently it remains unclear whether a so-called quiescent state was found in this campaign.

The detection of gamma rays from a small number of Narrow Line Seyfert 1 galaxies by the LAT instrument onboard Fermi seriously challenged our understanding of AGN physics. Among the most important findings associated with their discovery has been the realisation that smaller-mass black holes seem to be hosted by these systems. Immediately after their discovery a radio multi- frequency monitoring campaign was initiated to understand their jet radio emission. Here the first results of the campaign are presented. The light curves and some first variability analyses are discussed, showing that the brightness temperatures and Doppler factors are moderate. The phenomenologies are typically blazar-like. The frequency domain on the other hand indicates intense spectral evolution and the variability patterns indicate mechanisms similar to those acting in the jets of BL Lacs and FSRQs. Finally, the linear polarisation also reveals the presence of a quiescent, optically thin jet in certain cases.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

We present the first detection of the l-C3H+ hydrocarbon in the interstellar medium. The Horsehead WHISPER project, a millimeter unbiased line survey at two positions, namely the photo-dissociation region (PDR) and the nearby shielded core, revealed a consistent set of eight unidentified lines toward the PDR position. Six of them are detected with a signal-to-noise ratio from 6 to 19, while the two last ones are tentatively detected. Mostly noise appears at the same frequency toward the dense core, located less than 40" away. We simultaneously fit 1) the rotational and centrifugal distortion constants of a linear rotor, and 2) the Gaussian line shapes located at the eight predicted frequencies. The observed lines can be accurately fitted with a linear rotor model, implying a 1Sigma ground electronic state. The deduced rotational constant value is Be= 11244.9512 +/- 0.0015 MHz, close to that of l-C3H. We thus associate the lines to the l-C3H+ hydrocarbon cation, which enables us to constrain the chemistry of small hydrocarbons. A rotational diagram is then used to infer the excitation temperature and the column density. We finally compare the abundance to the results of the Meudon PDR photochemical model.

We investigate the radio and gamma-ray variability of the flat spectrum radio quasar PKS 1510-089 in the time range between 2010 November and 2012 January. In this period the source showed an intense activity, with two major gamma-ray flares detected in 2011 July and October. During the latter episode both the gamma-ray and the radio flux density reached their historical peak. Multiwavelength analysis shows a rotation of about 380 deg of the optical polarization angle close in time with the rapid and strong gamma-ray flare in 2011 July. An enhancement of the optical emission and an increase of the fractional polarization both in the optical and in radio bands is observed about three weeks later, close in time with another gamma-ray outburst. On the other hand, after 2011 September a huge radio outburst has been detected, first in the millimeter regime followed with some time delay at centimeter down to decimeter wavelengths. This radio flare is characterized by a rising and a decaying stage, in agreement with the formation of a shock and its evolution, as a consequence of expansion and radiative cooling. If the gamma-ray flare observed in 2011 October is related to this radio outburst, then this strongly indicates that the region responsible for the gamma-ray variability is not within the broad line, but a few parsecs downstream along the jet.

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.

One of the least understood properties of comets is the compositional structure of their nuclei, which can either be homogeneous or heterogeneous. The nucleus structure can be conveniently studied at millimeter wavelengths, using velocity-resolved spectral time series of the emission lines, obtained simultaneously for multiple molecules as the body rotates. Using this technique, we investigated the sources of CH3OH and HCN in comet 103P/Hartley 2, the target of NASA's EPOXI mission, which had an exceptionally favorable apparition in late 2010. Our monitoring with the IRAM 30 m telescope shows short-term variability of the spectral lines caused by nucleus rotation. The varying production rates generate changes in brightness by a factor of 4 for HCN and by a factor of 2 for CH3OH, and they are remarkably well correlated in time. With the addition of the velocity information from the line profiles, we identify the main sources of outgassing: two jets, oppositely directed in a radial sense, and icy grains, injected into the coma primarily through one of the jets. The mixing ratio of CH3OH and HCN is dramatically different in the two jets, which evidently shows large-scale chemical heterogeneity of the nucleus. We propose a network of identities linking the two jets with morphological features reported elsewhere, and postulate that the chemical heterogeneity may result from thermal evolution. The model-dependent average production rates are 3.5x10**26 molec/s for CH3OH and 1.25x10**25 molec/s for HCN, and their ratio of 28 is rather high but not abnormal. The rotational temperature from CH3OH varied strongly, presumably due to nucleus rotation, with the average value being 47 K.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

The F-GAMMA program is among the most comprehensive programs that aim at understanding the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are presented. (a) We discuss that the variability patterns traced can be classified into two classes: (1) to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic behaviour. We show that a simple two-component model can very well reproduce the observed phenomenologies. (b) We present the cm-to-mm behaviour of three gamma-ray bright Narrow Line Seyfert 1 galaxies over time spans varying between ~1.5 and 3 years and compare their variability characteristics with typical blazars.

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.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

We discuss the parsec-scale structural variability of the extragalactic jet 3C 111 related to a major radio flux density outburst in 2007. The data analyzed were taken within the scope of the MOJAVE, UMRAO, and F-GAMMA programs, which monitor a large sample of the radio brightest compact extragalactic jets with the VLBA, the University of Michigan 26 m, the Effelsberg 100 m, and the IRAM 30 m radio telescopes. The analysis of the VLBA data is performed by fitting Gaussian model components in the visibility domain. We associate the ejection of bright features in the radio jet with a major flux-density outburst in 2007. The evolution of these features suggests the formation of a leading component and multiple trailing components.

Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity to follow the radio spectra of a large sample of extragalactic sources to frequencies 2-20 times higher than allowed by past, large area, ground-based surveys. We combine the results of the Planck Early Release Compact Source Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as archival data, at frequencies below or overlapping Planck frequency bands, to validate the astrometry and photometry of the ERCSC radio sources and study the spectral features shown in this new frequency window opened by Planck. The ERCSC source positions and flux density scales are found to be consistent with the ground-based observations. We present and discuss the spectral energy distributions (SEDs) of a sample of "extreme" radio sources to illustrate the richness of the ERCSC for the study of extragalactic radio sources. Variability is found to play a role in the unusual spectral features of some of these sources.

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 F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of {\em Fermi} blazars. The current study is concerned with the broad-band radio spectra composed of measurement at ten frequencies between 2.64 and 142 GHz. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. The first four types are dominated by spectral evolution and can be reproduced by a simple two-component system made of the quiescent spectrum of a large scale jet populated with a flaring event evolving according to Marscher & Gear (1985). The last type is characterized by an achromatic change of the broad-band spectrum which must be attributed to a completely different mechanism. Here are presented, the classification, the assumed physical system and the results of simulations that have been conducted.

We study the relation between molecular gas and star formation in a volume-limited sample of 222 galaxies from the COLD GASS survey, with measurements of the CO(1-0) line from the IRAM 30m telescope. The galaxies are at redshifts 0.025<z<0.05 and have stellar masses in the range 10.0<log(M*/Msun)<11.5. The IRAM measurements are complemented by deep Arecibo HI observations and homogeneous SDSS and GALEX photometry. A reference sample that includes both UV and far-IR data is used to calibrate our estimates of star formation rates from the seven optical/UV bands. The mean molecular gas depletion timescale, tdep(H2), for all the galaxies in our sample is 1 Gyr, however tdep(H2) increases by a factor of 6 from a value of ~0.5 Gyr for galaxies with stellar masses of 10^10 Msun to ~3 Gyr for galaxies with masses of a few times 10^11 Msun. In contrast, the atomic gas depletion timescale remains contant at a value of around 3 Gyr. This implies that in high mass galaxies, molecular and atomic gas depletion timescales are comparable, but in low mass galaxies, molecular gas is being consumed much more quickly than atomic gas. The strongest dependences of tdep(H2) are on the stellar mass of the galaxy (parameterized as log tdep(H2)= (0.36+/-0.07)(log M* - 10.70)+(9.03+/-0.99)), and on the specific star formation rate. A single tdep(H2) versus sSFR relation is able to fit both "normal" star-forming galaxies in our COLD GASS sample, as well as more extreme starburst galaxies (LIRGs and ULIRGs), which have tdep(H2) < 10^8 yr. Normal galaxies at z=1-2 are displaced with respect to the local galaxy population in the tdep(H2) versus sSFR plane and have molecular gas depletion times that are a factor of 3-5 times longer at a given value of sSFR due to their significantly larger gas fractions.

We present results of 2 mm observations of the Crab Nebula, obtained using the Goddard-IRAM Superconducting 2 Millimeter Observer (GISMO) bolometer camera on the IRAM 30 m telescope. Additional 3.3 mm observations with the MUSTANG bolometer array on the Green Bank Telescope are also presented. The integrated 2 mm flux density of the Crab Nebula provides no evidence for the emergence of a second synchrotron component that has been proposed. It is consistent with the radio power law spectrum, extrapolated up to a break frequency of log(nu_{b} [GHz]) = 2.84 +/- 0.29 or nu_{b} = 695^{+651}_{-336} GHz. The Crab Nebula is well-resolved by the ~16.7" beam (FWHM) of GISMO. Comparison to radio data at comparable spatial resolution enables us to confirm significant spatial variation of the spectral index between 21 cm and 2 mm. The main effect is a spectral flattening in the inner region of the Crab Nebula, correlated with the toroidal structure at the center of the nebula that is prominent in the near-IR through X-ray regime.

The nuclei of active comets emit molecules anisotropically from discrete vents. As the nucleus rotates, we expect to observe periodic variability in the molecular emission line profiles, which can be studied through mm/submm spectroscopy. Using this technique we investigated the HCN atmosphere of comet 103P/Hartley 2, the target of NASA's EPOXI mission, which had an exceptionally favorable apparition in late 2010. We detected short-term evolution of the spectral line profile, which was stimulated by the nucleus rotation, and which provides evidence for rapid deceleration and excitation of the rotation state. The measured rate of change in the rotation period is +1.00 \pm 0.15 min per day and the period itself is 18.32 \pm 0.03 hr, both applicable at the epoch of the EPOXI encounter. Surprisingly, the spin-down efficiency is lower by two orders of magnitude than the measurement in comet 9P/Tempel 1 and the best theoretical prediction. This secures rotational stability of the comet's nucleus during the next few returns, although we anticipate a catastrophic disruption from spin-up as its ultimate fate.

We are conducting COLD GASS, a legacy survey for molecular gas in nearby galaxies. Using the IRAM 30m telescope, we measure the CO(1-0) line in a sample of ~350 nearby (D=100-200 Mpc), massive galaxies (log(M*/Msun)>10.0). The sample is selected purely according to stellar mass, and therefore provides an unbiased view of molecular gas in these systems. By combining the IRAM data with SDSS photometry and spectroscopy, GALEX imaging and high-quality Arecibo HI data, we investigate the partition of condensed baryons between stars, atomic gas and molecular gas in 0.1-10L* galaxies. In this paper, we present CO luminosities and molecular hydrogen masses for the first 222 galaxies. The overall CO detection rate is 54%, but our survey also uncovers the existence of sharp thresholds in galaxy structural parameters such as stellar mass surface density and concentration index, below which all galaxies have a measurable cold gas component but above which the detection rate of the CO line drops suddenly. The mean molecular gas fraction MH2/M* of the CO detections is 0.066+/-0.039, and this fraction does not depend on stellar mass, but is a strong function of NUV-r colour. Through stacking, we set a firm upper limit of MH2/M*=0.0016+/-0.0005 for red galaxies with NUV-r>5.0. The average molecular-to-atomic hydrogen ratio in present-day galaxies is 0.3, with significant scatter from one galaxy to the next. The existence of strong detection thresholds in both the HI and CO lines suggests that "quenching" processes have occurred in these systems. Intriguingly, atomic gas strongly dominates in the minority of galaxies with significant cold gas that lie above these thresholds. This suggests that some re-accretion of gas may still be possible following the quenching event.

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

We have used the new 90GHz MUSTANG camera on the Robert C. Byrd Green Bank Telescope (GBT) to map the bright Huygens region of the star-forming region M42 with a resolution of 9" and a sensitivity of 2.8mJy/beam. 90GHz is an interesting transition frequency, as MUSTANG detects both the free-free emission characteristic of the HII region created by the Trapezium stars, normally seen at lower frequencies, and thermal dust emission from the background OMC1 molecular cloud, normally mapped at higher frequencies. We also present similar data from the 150GHz GISMO camera taken on the IRAM telescope. This map has 15" resolution. By combining the MUSTANG data with 1.4, 8, and 21GHz radio data from the VLA and GBT, we derive a new estimate of the emission measure (EM) averaged electron temperature of Te = 11376K by an original method relating free-free emission intensities at optically thin and optically thick frequencies. Combining ISO-LWS data with our data, we derive a new estimate of the dust temperature and spectral emissivity index within the 80" ISO-LWS beam toward OrionKL/BN, Td = 42K and Beta=1.3. We show that both Td and Beta decrease when going from the HII region and excited OMC1 interface to the denser UV shielded part of OMC1 (OrionKL/BN, Orion S). With a model consisting of only free-free and thermal dust emission we are able to fit data taken at frequencies from 1.5GHz to 854GHz.

Based on millimeter-wavelength continuum observations we suggest that the recent 'spectacle' of comet 17P/Holmes can be explained by a thick, air-tight dust cover and the effects of H2O sublimation, which started when the comet arrived at the heliocentric distance <= 2.5 AU. The porous structure inside the nucleus provided enough surface for additional sublimation, which eventually led to the break up of the dust cover and to the observed outburst. The magnitude of the particle burst can be explained by the energy provided by insolation, stored in the dust cover and the nucleus within the months before the outburst: the subliming surface within the nucleus is more than one order of magnitude larger than the geometric surface of the nucleus -- possibly an indication of the latter's porous structure. Another surprise is that the abundance ratios of several molecular species with respect to H2O are variable. During this apparition, comet Holmes lost about 3% of its mass, corresponding to a 'dirty ice' layer of 20m.

The Star Formation (SF) rate in galaxies is an important parameter at all redshifts and evolutionary stages of galaxies. In order to understand the increased SF rates in intermediate redshift galaxies one possibility is to study star formation in local galaxies with properties frequently found at this earlier epoch like low metallicity and small size. We present sensitive observations of the molecular gas in M 33, a small Local Group spiral at a distance of 840 kpc which shares many of the characteristics of the intermediate redshift galaxies. The observations were carried out in the CO(2--1) line with the HERA heterodyne array on the IRAM 30 m telescope. A 11\arcmin$\times$22\arcmin region in the northern part of M 33 was observed, reaching a detection threshold of a few 10$^{3}$ \msol. The correlation in this field between the CO emission and tracers of SF (8\mum, 24\mum, \Ha, FUV) is excellent and CO is detected very far North, showing that molecular gas forms far out in the disk even in a small spiral with a subsolar metallicity. One major molecular cloud was discovered in an interarm region with no HI peak and little if any signs of SF -- without a complete survey this cloud would never have been found. The radial dependence of the CO emission has a scale length similar to the dust emission, less extended than the \Ha or FUV. If, however, the $\ratioo$ ratio varies inversely with metallicity, then the scale length of the H$_2$ becomes similar to that of the \Ha or FUV. Comparing the SF rate to the H$_2$ mass shows that M 33, like the intermediate redshift galaxies it resembles, has a significantly higher SF efficiency than large local universe spirals.

We discuss the interpretation of the dust SED from the mid-infrared to the millimeter range of NGC 1569. The model developed by D\'esert et al. (1990) including three dust components (Polyaromatic Hydrocarbons, Very Small Grains and big grains) can explain the data using a realistic interstellar radiation field and adopting an enhanced abundance of VSGs. A simple three-temperature model is also able to reproduce the data but requires a very low dust temperature which is considered to be unlikely in this low-metallicity starburst galaxy. The high abundance of Very Small Grains might be due to large grain destruction in supernova shocks. This possibility is supported by ISO data showing that the emission at 14.3 $\mu$m, tracing VSGs, is enhanced with respect to the emission at 6.7 $\mu$m and 850 $\mu$m in regions of high star formation.

We present new data of the dwarf galaxy NGC 1569 at 450 $\mu$m, 850 $\mu$m and 1200$\mu$m taken with SCUBA at the JCMT and the bolometer array at the IRAM 30m telescope. After including data from IRAS at 12, 25, 60 and 100 $\mu$m, we have successfully fitted the dust grain population model of D\'esert et al. (1990) to the observed midinfrared-to-millimeter spectrum. The fit requires a combination of both large and very small grains exposed to a strong radiation field as well as an enhancement of the number of very small grains relative to the number of large grains. We interpret this as the consequence of large grain destruction due to shocks in the turbulent interstellar medium of NGC 1569. The contribution of polyaromatic hydrocarbons (PAH's) is found to be negligible. Comparison of the dust emission maps with an HI map of similar resolution shows that both dust and molecular gas distributions peak close to the radio continuum maximum and at a minimum in the HI distribution. From a comparison of these three maps and assuming that the gas-to-dust mass ratio is the same everywhere, we estimate the ratio of molecular hydrogen column density to integrated CO intensity to be about 25 -- 30 times the local Galactic value. The gas-to-dust ratio is 1500 -- 2900, about an order of magnitude higher than in the Solar Neighbourhood.

X-ray and optical observations have shown that the supernova explosions and stellar winds of the starburst in NGC 2146 produce an outflow of hot material along the minor axis. This outflow is expected to have a more or less conical shape, on either side of the galactic plane, and cone walls of cooler material where the outflow is in shock contact with halo gas. We attempt to determine the geometry (diameter at the base and opening angle) and the physical parameters (velocity and density) of the material in the cone walls from the optical emission line and radio observations presented here, and from published X-ray, radio, and optical observations. We compare the outflow of NGC 2146 with the outflow of M 82.

We report on observations of CO(1-0) emission from the chemically young starburst galaxy Mrk 109. These observations were part of a study to constrain the O2/CO ratio in metal-deficient galaxies, which were motivated by theoretical work that suggests the possible enhancement of the O2/CO ratio in chemically young systems. Five low metallicity (Z < 0.5 Z(solar)) IRAS galaxies at redshifts z > 0.02 (required to shift the 118.75 GHz O2 line away from the atmospheric line) were searched for CO emission. We detected the CO line in only Mrk 109. From O2 observations of Mrk 109, we achieved an upper limit for the O2 column density abundance ratio of N(O2)/N(CO) < 0.31. These results provide useful constraints for the theoretical models of chemically young galaxies. We argue that either most of the molecular gas in Mrk 109 does not reside in dark clouds (Av >~ 5), or the standard equilibrium chemistry models are inadequate for metal-poor systems. The molecular gas mass implied by the CO observations of Mrk 109 is M(H2) ~= 4 x 10^{9} M(sun), and the CO data are consistent with a central starburst induced by the interaction with a nearby companion.