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.

Saturn's moon Iapetus, which is in synchronous rotation, is covered by an optically dark material mainly on its leading side, while its trailing side is significantly brighter. Because longer wavelengths probe deeper into the subsurface, observing both sides at a variety of wavelengths brings to light possible changes in thermal, compositional, and physical properties with depth. We have observed Iapetus's leading and trailing hemispheres at 1.2 and 2.0 mm, using the NIKA2 camera mounted on the IRAM 30-m telescope, and compared our observations to others performed at mm to cm wavelengths. We calibrate our observations on Titan, which is simultaneously observed within the field of view. Due to the proximity of Saturn, it is sometimes difficult to separate Iapetus's and Titan's flux from that of Saturn, detected in the telescope's side lobes. Preliminary results show that the trailing hemisphere brightness temperatures at the two wavelengths are equal within error bars, unlike the prediction made by Ries (2012). On the leading side, we report a steep spectral slope of increasing brightness temperature (by 10 K) from 1.2 to 2.0 mm, which may indicate rapidly varying emissivities within the top few centimeters of the surface. Comparison to a diffuse scattering model and a thermal model will be necessary to further constrain the thermophysical properties of the subsurface of Iapetus's two faces.

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)

Debris discs are second generation dusty discs formed by collisions of planetesimals. Many debris discs have been found and resolved around hot and solar-type stars. However, only a handful have been discovered around M-stars, and the reasons for their paucity remain unclear. Here we check whether the sensitivity and wavelength coverage of present-day telescopes are simply unfavourable for detection of these discs or if they are truly rare. We approach this question by looking at the Herschel/DEBRIS survey that has searched for debris discs including M-type stars. Assuming that these cool-star discs are "similar" to those of the hotter stars in some sense (i.e., in terms of dust location, temperature, fractional luminosity, or mass), we check whether this survey should have found them. With our procedure we can reproduce the $2.1^{+4.5}_{-1.7}$% detection rate of M-star debris discs of the DEBRIS survey, which implies that these discs can indeed be similar to discs around hotter stars and just avoid detection. We then apply this procedure to IRAM NIKA-2 and ALMA bands 3, 6 and 7 to predict possible detection rates and give recommendations for future observations. We do not favour observing with IRAM, since it leads to detection rates lower than for the DEBRIS survey, with 0.6%-4.5% for a 15 min observation. ALMA observations, with detection rates 0.9%-7.4%, do not offer a significant improvement either, and so we conclude that more sensitive far-infrared and single dish sub-millimetre telescopes are necessary to discover the missing population of M-star debris discs.

Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts.

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)

We present the results of a Herschel survey of 21 late-type stars that host planets discovered by the radial velocity technique. The aims were to discover new disks in these systems and to search for any correlation between planet presence and disk properties. In addition to the known disk around GJ 581, we report the discovery of two new disks, in the GJ 433 and GJ 649 systems. Our sample therefore yields a disk detection rate of 14%, higher than the detection rate of 1.2% among our control sample of DEBRIS M-type stars with 98% confidence. Further analysis however shows that the disk sensitivity in the control sample is about a factor of two lower in fractional luminosity than for our survey, lowering the significance of any correlation between planet presence and disk brightness below 98%. In terms of their specific architectures, the disk around GJ 433 lies at a radius somewhere between 1 and 30au. The disk around GJ 649 lies somewhere between 6 and 30au, but is marginally resolved and appears more consistent with an edge-on inclination. In both cases the disks probably lie well beyond where the known planets reside (0.06-1.1au), but the lack of radial velocity sensitivity at larger separations allows for unseen Saturn-mass planets to orbit out to $\sim$5au, and more massive planets beyond 5au. The layout of these M-type systems appears similar to Sun-like star + disk systems with low-mass planets.

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.}

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.

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.

We present the first ALMA observations of the closest known extrasolar debris disc. This disc orbits the star $\epsilon$ Eridani, a K-type star just 3.2pc away. Due to the proximity of the star, the entire disc cannot fit within the ALMA field of view. Therefore, the observations have been centred 18" North of the star, providing us with a clear detection of the northern arc of the ring, at a wavelength of 1.3mm. The observed disc emission is found to be narrow with a width of just 11-13AU. The fractional disc width we find is comparable to that of the Solar System's Kuiper Belt and makes this one of the narrowest debris discs known. If the inner and outer edges are due to resonances with a planet then this planet likely has a semi-major axis of 48AU. We find tentative evidence for clumps in the ring, although there is a strong chance that at least one is a background galaxy. We confirm, at much higher significance, the previous detection of an unresolved emission at the star that is above the level of the photosphere and attribute this excess to stellar chromospheric emission.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

While most of the known debris discs present cold dust at tens of AU, a few young systems exhibit hot dust analogous to the Zodiacal dust. $\eta$ Corvi is particularly interesting as it is old and it has both, with its hot dust significantly exceeding the maximum luminosity of an in-situ collisional cascade. Previous work suggested that this system could be undergoing an event similar to the Late Heavy Bombardment (LHB) soon after or during a dynamical instability. Here we present ALMA observations of $\eta$ Corvi with a resolution of 1."2 (~22au) to study its outer belt. The continuum emission is consistent with an axisymmetric belt, with a mean radius of 152au and radial FWHM of 46au, which is too narrow compared to models of inward scattering of an LHB-like scenario. Instead, the hot dust could be explained as material passed inwards in a rather stable planetary configuration. We also report a 4sigma detection of CO at ~ 20au. CO could be released in situ from icy planetesimals being passed in when crossing the H$_2$O or CO$_2$ ice lines. Finally, we place constraints on hidden planets in the disc. If a planet is sculpting the disc's inner edge, this should be orbiting at 75-100au, with a mass of 3-30 M$_\oplus$ and an eccentricity < 0.08. Such a planet would be able to clear its chaotic zone on a timescale shorter than the age of the system and scatter material inwards from the outer belt to the inner regions, thus feeding the hot dust.

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.

We present observations of Epsilon Eridani from the Submillimeter Array (SMA) at 1.3 millimeters and from the Australia Telescope Compact Array (ATCA) at 7 millimeters that reach an angular resolution of ~4" (13 AU). These first millimeter interferometer observations of Epsilon Eridani, which hosts the closest debris disk to the Sun, reveal two distinct emission components: (1) the well-known outer dust belt, which, although patchy, is clearly resolved in the radial direction, and (2) an unresolved source coincident with the position of the star. We use direct model-fitting of the millimeter visibilities to constrain the basic properties of these two components. A simple Gaussian shape for the outer belt fit to the SMA data results in a radial location of $64.4^{+2.4}_{-3.0}$ AU and FWHM of $20.2^{+6.0}_{-8.2}$ AU (fractional width $\Delta R/R = 0.3$. Similar results are obtained taking a power law radial emission profile for the belt, though the power law index cannot be usefully constrained. Within the noise obtained (0.2 mJy/beam), these data are consistent with an axisymmetric belt model and show no significant azimuthal structure that might be introduced by unseen planets in the system. These data also limit any stellocentric offset of the belt to $<9$ AU, which disfavors the presence of giant planets on highly eccentric ($>0.1$) and wide (10's of AU) orbits. The flux density of the unresolved central component exceeds predictions for the stellar photosphere at these long wavelengths, by a marginally significant amount at 1.3 millimeters but by a factor of a few at 7 millimeters (with brightness temperature $13000 \pm 1600$ K for a source size of the optical stellar radius). We attribute this excess emission to ionized plasma from a stellar corona or chromosphere.

Context. ABDoradus is the main system of the ABDoradus moving group. It is a quadruple system formed by two widely separated binaries of pre-main-sequence (PMS) stars: ABDor A/C and ABDor Ba/Bb. The pair ABDor A/C has been extensively studied and its dynamical masses have been determined with high precision, thus making of ABDor C a benchmark for calibrating PMS stellar models. If the orbit and dynamical masses of the pair ABDor Ba/Bb can be determined, they could not only play a similar role to that of ABDor C in calibrating PMS models, but would also help to better understand the dynamics of the whole ABDoradus system. Aims. We aim to determine the individual masses of the pair ABDor Ba/Bb using VLBI observations and archive infrared data, as part of a larger program directed to monitor binary systems in the ABDoradus moving group. Methods. We observed the system ABDor B between 2007 and 2013 with the Australian Long Baseline Array (LBA), at a frequency of 8.4 GHz in phase-reference mode. Results. We detected, for the first time, compact radio emission from both stars in the binary, ABDor Ba and ABDor Bb. This result allowed us to determine the orbital parameters of both the relative and absolute orbits and, consequently, their individual dynamical masses: 0.28+/-0.05 Msun and 0.25+/-0.05 Msun, respectively. Conclusions. Comparisons of the dynamical masses with the prediction of PMS evolutionary models show that the models underpredict the dynamical masses of the binary components Ba and Bb by 10-30% and 10-40%, respectively, although they all still agree at the 2-sigma level. Some of stellar models considered favour an age between 50 and 100 Myr for this system, meanwhile others predict older ages. We also discuss the evolutionary status of ABDor Ba/Bb in terms of an earlier double-double star scenario that might explain the strong radio emission detected in both components.

The debris disk closest to Earth is the one around the star epsilon Eridani at a distance of 3.2 pc. It is the prime target for detailed studies of a belt of planetesimals left from the early phase of planet formation other than the Kuiper Belt. The non-uniform ring-like structure around epsilon Eridani, originally discovered at lambda=850 microns with the bolometer camera SCUBA, could be the signpost of unseen long-period planets interior to the disk that gravitationally interact with it through mean-motion resonances. However, the reliability of the structure at 850 microns, which has been debated, has not been verified with independent observations until now. We present a high signal-to-noise ratio image of this structure at lambda=1.2 mm made with the bolometer camera MAMBO and compare this with the SCUBA image. We have found that three of the four emission clumps (NE, NW, SW) and the two deep hollows to the east and west are at the same positions in the MAMBO and SCUBA images within astrometric uncertainty. The SE clump is at odds, significantly brighter and more extended in the SCUBA than in the MAMBO images, but it is possible that this mismatch is an artifact. We conclude that this degree of positional coincidence provides tentative evidence that the observed structure is robust. In addition, we present the radial brightness profile of our MAMBO image and show that the width of the planetesimal belt around epsilon Eridani is narrower than 22 AU, a more stringent upper limit than determined from previous observations. The corresponding relative width is $0.1 < \Delta R / R < 0.4$, which is lower than for the Kuiper Belt.

The study of the planet-debris disk connection can shed light on the formation and evolution of planetary systems, and may help predict the presence of planets around stars with certain disk characteristics. In preliminary analyses of the Herschel DEBRIS and DUNES surveys, Wyatt et al. (2012) and Marshall et al. (2014) identified a tentative correlation between debris and low-mass planets. Here we use the cleanest possible sample out these surveys to assess the presence of such a correlation, discarding stars without known ages, with ages < 1 Gyr and with binary companions <100 AU, to rule out possible correlations due to effects other than planet presence. In our sample of 204 FGK stars, we do not find evidence that debris disks are more common or more dusty around stars harboring high-mass or low-mass planets compared to a control sample without identified planets, nor that debris disks are more or less common (or more or less dusty) around stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may account for the lack of correlations. The data show the correlation between the presence of high-mass planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to those expected from a Gaussian distribution, we find that a distribution centered on the Solar system's value fits well the data, while one centered at 10 times this value can be rejected. This is of interest in the context of future terrestrial planet characterization because it indicates that there are good prospects for finding a large number of debris disk systems (i.e. with evidence of harboring the building blocks of planets) with exozodiacal emission low enough to be appropriate targets for an ATLAST-type mission to search for biosignatures.

Debris disks have been found primarily around intermediate and solar mass stars (spectral types A-K), but rarely around low-mass M-type stars. This scarcity of detections in M star surveys can be confronted with the predictions of the steady state collisional evolution model. First, we determine the parameters of the disk population evolved with this model and fit to the distribution of the fractional dust luminosities measured in the surveys of A- and FGK-type stars observed by the infrared satellite Spitzer. Thus, in our approach, we stipulate that the initial disk mass distribution is bimodal and that only high-mass collisionally-dominated disks are detected. The best determined parameter is the diameter $D_c$ of the largest planetesimals in the collisional cascade of the model, which ranges between 2 and 60 km, consistently for disks around A- and FGK-type stars. Second, we assume that the same disk population surrounds the M dwarfs that have been the subjects of debris disk searches in the far-infrared with Spitzer and at submillimeter wavelengths with radiotelescopes. We find, in the framework of our study, that this disk population, which has been fit to the AFGK data, is still consistent with the observed lack of disks around M dwarfs with Spitzer.

The $\textit{Herschel}$ DEBRIS, DUNES and GT programmes observed 37 exoplanet host stars within 25 pc at 70, 100 and 160 $\mu$m with the sensitivity to detect far-infrared excess emission at flux density levels only an order of magnitude greater than that of the Solar system's Edgeworth-Kuiper belt. Here we present an analysis of that sample, using it to more accurately determine the (possible) level of dust emission from these exoplanet host stars and thereafter determine the links between the various components of these exoplanetary systems through statistical analysis. We have fitted the flux densities measured from recent \textit{Herschel} observations with a simple two parameter ($T_{d}$, $L_{\rm IR}/L_{\star}$) black body model (or to the 3-$\sigma$ upper limits at 100 $\mu$m). From this uniform approach we calculate the fractional luminosity, radial extent, dust temperature and disc mass. We then plotted the calculated dust luminosity or upper limits against the stellar properties, e.g. effective temperature, metallicity, age, and identified correlations between these parameters. A total of eleven debris discs are identified around the 37 stars in the sample. An incidence of ten cool debris discs around the Sun-like exoplanet host stars (29 $\pm$ 9 %) is consistent with the detection rate found by DUNES (20.2 $\pm$ 2.0 %). For the debris disc systems, the dust temperatures range from 20 to 80 K, and fractional luminosities ($L_{\rm IR}/L_{\star}$) between 2.4 $\times$10$^{-6}$ and 4.1 $\times$10$^{-4}$. In the case of non-detections, we calculated typical 3-$\sigma$ upper limits to the dust fractional luminosities of a few $\times10^{-6}$. We recover the previously identified correlation between stellar metallicity and hot Jupiter planets in our data set. We find a correlation between the increased presence of dust, lower planet masses and lower stellar metallicities. (abridged)

Fomalhaut is one of the most interesting and well studied nearby stars, hosting at least one planet, a spectacular debris ring, and two distant low-mass stellar companions (TW PsA and LP 876-10, a.k.a. Fomalhaut B & C). We observed both companions with Herschel, and while no disc was detected around the secondary, TW PsA, we have discovered the second debris disc in the Fomalhaut system, around LP 876-10. This detection is only the second case of two debris discs seen in a multiple system, both of which are relatively wide ($\gtrsim$3000 AU for HD 223352/40 and 158 kAU [0.77 pc] for Fomalhaut/LP 876-10). The disc is cool (24K) and relatively bright, with a fractional luminosity $L_{\rm disc}/L_\star = 1.2 \times 10^{-4}$, and represents the rare observation of a debris disc around an M dwarf. Further work should attempt to find if the presence of two discs in the Fomalhaut system is coincidental, perhaps simply due to the relatively young system age of 440 Myr, or if the stellar components have dynamically interacted and the system is even more complex than it currently appears.

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts.

Debris disks have been found primarily around intermediate and solar mass stars (spectral types A-K) but rarely around low mass M-type stars. We have spatially resolved a debris disk around the remarkable M3-type star GJ581 hosting multiple planets using deep PACS images at 70, 100 and 160 microns as part of the DEBRIS Program on the Herschel Space Observatory. This is the second spatially resolved debris disk found around an M-type star, after the one surrounding the young star AU Mic (12 Myr). However, GJ 581 is much older (2-8 Gyr), and is X-ray quiet in the ROSAT data. We fit an axisymmetric model of the disk to the three PACS images and found that the best fit model is for a disk extending radially from 25+/-12 AU to more than 60 AU. Such a cold disk is reminiscent of the Kuiper Belt but it surrounds a low mass star (0.3 M_sol) and its fractional dust luminosity L_dust/L_* of \sim 10^-4 is much higher. The inclination limits of the disk found in our analysis make the masses of the planets small enough to ensure the long-term stability of the system according to some dynamical simulations. The disk is collisionally dominated down to submicron-sized grains and the dust cannot be expelled from the system by radiation or wind pressures because of the low luminosity and low X-ray luminosity of GJ581. We suggest that the correlation between low-mass planets and debris disks recently found for G-type stars also applies to M-type stars.Finally, the known planets, of low masses and orbiting within 0.3 AU from the star, cannot dynamically perturb the disk over the age of the star, suggesting that an additional planet exists at larger distance that is stirring the disk to replenish the dust.

This paper describes Herschel observations of the nearby (8.5pc) G5V multi-exoplanet host star 61 Vir at 70-500micron carried out by the DEBRIS survey. These reveal emission that is extended out to >15arcsec with a morphology that can be fitted by a nearly edge-on (77deg inclination) radially broad (from 30AU to >100AU) debris disk of fractional luminosity 2.7x10^-5, with two unrelated sources nearby that are more prominent at longer wavelengths. Chance alignment with a background object seen at 1.4GHz provides potential for confusion, but the star's 1.4"/yr proper motion allows Spitzer 70micron images to confirm that what we are interpreting as disk emission really is circumstellar. Although the exact shape of the disk's inner edge is not well constrained, the region inside 30AU must be significantly depleted in planetesimals. This is readily explained if there are additional planets in the 0.5-30AU region, but is also consistent with collisional erosion. We also find tentative evidence that the presence of detectable debris around nearby stars correlates with the presence of the lowest mass planets that are detectable in current radial velocity surveys. Out of an unbiased sample of the nearest 60 G stars, 11 are known to have planets, of which 6 (including 61 Vir) have planets that are all less massive than Saturn, and 4 of these have evidence for debris. The debris toward one of these planet-hosts (HD20794) is reported here for the first time. This fraction (4/6) is higher than that expected for nearby field stars (15%), and implies that systems that form low-mass planets are also able to retain bright debris disks. We suggest that this correlation could arise because such planetary systems are dynamically stable and include regions that are populated with planetesimals in the formation process where the planetesimals can remain unperturbed over Gyr timescales.

We describe the NASA/Stanford gyroscope relativity mission, Gravity Probe B (GP-B), and provide an overview of the following series of six astrometric and astrophysical papers that report on our radio observations and analyses made in support of this mission. The main goal of this 8.5 year program of differential VLBI astrometry was to determine the proper motion of the guide star of the GP-B mission, the RS CVn binary IM Pegasi (IM Peg; HR 8703). This proper motion is determined with respect to compact, extragalactic reference sources. The results are: -20.833 +- 0.090 mas/yr and -27.267 +- 0.095 mas/yr for, respectively, the right ascension and declination, in local Cartesian coordinates, of IM Peg's proper motion, and 10.370 +- 0.074 mas (i.e., 96.43 +- 0.69 pc) for its parallax (and distance). Each quoted uncertainty is meant to represent an ~70% confidence interval that includes the estimated contribution from systematic error. These results are accurate enough not to discernibly degrade the GP-B estimates of its gyroscopes' relativistic precessions: the frame-dragging and geodetic effects.

We present a physical interpretation for the locations of the sources of radio emission in IM Pegasi (IM Peg, HR 8703), the guide star for the NASA/Stanford relativity mission Gravity Probe B. This emission is seen in each of our 35 epochs of 8.4-GHz VLBI observations taken from 1997 to 2005. We found that the mean position of the radio emission is at or near the projected center of the primary to within about 27% of its radius, identifying this active star as the radio emitter. The positions of the radio brightness peaks are scattered across the disk of the primary and slightly beyond, preferentially along an axis with position angle, p.a. = (-38 +- 8) deg, which is closely aligned with the sky projections of the orbit normal (p.a. = -49.5 +- 8.6 deg) and the expected spin axis of the primary. Comparison with simulations suggests that brightness peaks are 3.6 (+0.4,-0.7) times more likely to occur (per unit surface area) near the pole regions of the primary (|latitude| >= 70 deg) than near the equator (|latitude| <= 20 deg), and to also occur close to the surface with ~2/3 of them at altitudes not higher than 25% of the radius of the primary.

We present the principal astrometric results of the very-long-baseline interferometry (VLBI) program undertaken in support of the Gravity Probe B (GP-B) relativity mission. VLBI observations of the GP-B guide star, the RS CVn binary IM Pegasi (HR 8703), yielded positions at 35 epochs between 1997 and 2005. We discuss the statistical assumptions behind these results and our methods for estimating the systematic errors. We find the proper motion of IM Peg in an extragalactic reference frame closely related to the International Celestial Reference Frame 2 (ICRF2) to be -20.83 +- 0.03 +- 0.09 mas/yr in right ascension and -27.27 +- 0.03 +- 0.09 mas/yr in declination. For each component the first uncertainty is the statistical standard error and the second is the total standard error (SE) including plausible systematic errors. We also obtain a parallax of 10.37 +- 0.07 mas (distance: 96.4 +- 0.7 pc), for which there is no evidence of any significant contribution of systematic error. Our parameter estimates for the ~25-day-period orbital motion of the stellar radio emission have SEs corresponding to ~0.10 mas on the sky in each coordinate. The total SE of our estimate of IM Peg's proper motion is ~30% smaller than the accuracy goal set by the GP-B project before launch: 0.14 mas/yr for each coordinate of IM Peg's proper motion. Our results ensure that the uncertainty in IM Peg's proper motion makes only a very small contribution to the uncertainty of the GP-B relativity tests.

When VLBI observations are used to determine the position or motion of a radio source relative to reference sources nearby on the sky, the astrometric information is usually obtained via: (i) phase-referenced maps; or (ii) parametric model fits to measured fringe phases or multiband delays. In this paper we describe a "merged" analysis technique which combines some of the most important advantages of these other two approaches. In particular, our merged technique combines the superior model-correction capabilities of parametric model fits with the ability of phase-referenced maps to yield astrometric measurements of sources that are too weak to be used in parametric model fits. We compare the results from this merged technique with the results from phase-referenced maps and from parametric model fits in the analysis of astrometric VLBI observations of the radio-bright star IM Pegasi (HR 8703) and the radio source B2252+172 nearby on the sky. In these studies we use central-core components of radio sources 3C 454.3 and B2250+194 as our positional references. We obtain astrometric results for IM Peg with our merged technique even when the source is too weak to be used in parametric model fits, and we find that our merged technique yields superior astrometric results to the phase-referenced mapping technique. We used our merged technique to estimate the proper motion and other astrometric parameters of IM Peg in support of the NASA/Stanford Gravity Probe B mission.

We present secure [NII[ detections in two mm-bright, strongly lensed objects at high redshift, APM08279+5255 (z=3.911) and MM18423+5938 (z=3.930), using the IRAM Plateau de Bure Interferometer. Due to its ionization energy [NII] is a good tracer of the ionized gas phase in the interstellar medium. The measured fluxes are S([NII])=(4.8+/-0.8) Jy km/s and (7.4+/-0.5) Jy km/s respectively, yielding line luminosities of L([NII]) =(1.8+/-0.3) x 10^9 \mu^{-1} Lsun for APM08279+5255 and L([NII]) =(2.8+/-0.2) x 10^9 \mu^{-1} Lsun for MM18423+5938. Our high-resolution map of the [NII] and 1 mm continuum emission in MM18423+5938 clearly resolves an Einstein ring in this source, and reveals a velocity gradient in the dynamics of the ionized gas. A comparison of these maps with high-resolution EVLA CO observations enables us to perform the first spatially-resolved study of the dust continuum-to-molecular gas surface brightness (Sigma_{FIR} Sigma_{CO}^N, which can be interpreted as the star formation law) in a high-redshift object. We find a steep relation (N=1.4+/-0.2), consistent with a starbursting environment. We measure a [NII]/FIR luminosity ratio in APM0828+5255 and MM18423+5938 of 9.0 x 10^{-6} and 5.8 x 10^{-6}, respectively. This is in agreement with the decrease of the [NII]/FIR ratio at high FIR luminosities observed in local galaxies.

A debris disk is a constituent of any planetary system surrounding a main sequence star. We study whether close stellar encounters can disrupt and strip a debris disk of its planetesimals in the expanding open cluster of its birth with a decreasing star number density over 100 Myrs. Such stripping would affect the dust production and hence detectability of the disk. We tabulated the fractions of planetesimals stripped off during stellar flybys of miss distances between 100 and 1000 AU and for several mass ratios of the central to passing stars. We then estimated the numbers of close stellar encounters over the lifetime of several expanding open clusters characterized by their initial star densities. We found that a standard disk, with inner and outer radii of 40 and 100 AU, suffers no loss of planetesimals over 100 Myrs around a star born in a common embedded cluster with star density <1000 pc^-3. In contrast, we found that such a disk is severely depleted of its planetesimals over this timescale around a star born in an Orion-type cluster where the star density is >20 000 pc^-3. In this environment, a disk loses >97% of its planetesimals around an M-dwarf, >63% around a solar-type star, and >42% around an A-dwarf, over 100 Myrs. We roughly estimate that two-thirds of the stars may be born in such high star density clusters. This might explain in part why fewer debris disks are observed around lower mass stars.

We present high resolution imaging of the low order (J=1 and 2) CO line emission from the z = 3.93 submillimeter galaxy (SMG) MM18423+5938 using the Expanded Very Large Array, and optical and near-IR imaging using the Canada-France-Hawaii Telescope. This SMG with a spectroscopic redshift was thought to be gravitationally lensed given its enormous apparent brightness. We find that the CO emission is consistent with a complete Einstein ring with a major axis diameter of ~ 1.4", indicative of lensing. We have also identified the lens galaxy as a very red elliptical coincident with the geometric center of the ring and estimated its photometric redshift z~1.1. A first estimate of the lens magnification factor is m~12. The luminosity L'_CO(1-0) of the CO(1-0) emission is 2.71+/-0.38 x10^11/m K km s^-1 pc^2, and, adopting the commonly used conversion factor for ULIRGs, the molecular gas mass is M(H_2) = 2.2 x10^11/m M_sol, comparable to unlensed SMGs if corrected by m ~ 12. Our revised estimate of the far-IR luminosity of MM18423+5938 is 2 x10^13/m < L_FIR < 3 x10^14/m L_sol, comparable to that of ULIRGs. Further observations are required to quantify the star formation rate in MM18423+5938 and to constrain the mass model of the lens in more detail.

We have used multi-epoch long-baseline radio interferometry to determine the proper motion and orbital elements of Algol and UX Arietis, two radio-bright, close binary stellar systems with distant tertiary components. For Algol, we refine the proper motion and outer orbit solutions, confirming the recent result of Zavala et al. (2010) that the inner orbit is retrograde. The radio centroid closely tracks the motion of the KIV secondary. In addition, the radio morphology varies from double-lobed at low flux level to crescent-shaped during active periods. These results are most easily interpreted as synchrotron emission from a large, co-rotating meridional loop centered on the K-star. If this is correct, it provides a radio-optical frame tie candidate with an uncertainty {\pm}0.5 mas. For UX Arietis, we find a outer orbit solution that accounts for previous VLBI observations of an acceleration term in the proper motion fit. The outer orbit solution is also consistent with previously published radial velocity curves and speckle observations of a third body. The derived tertiary mass, 0.75 solar masses, is consistent with the K1 main-sequence star detected spectroscopically. The inner orbit solution favors radio emission from the active K0IV primary only. The radio morphology, consisting of a single, partially resolved emission region, may be associated with the persistent polar spot observed using Doppler imaging.

Pulsars PSR J0248+6021 (rotation period P=217 ms and spin-down power Edot = 2.13E35 erg/s) and PSR J2240+5832 (P=140 ms, Edot = 2.12E35 erg/s) were discovered in 1997 with the Nancay radio telescope during a northern Galactic plane survey, using the Navy-Berkeley Pulsar Processor (NBPP) filter bank. GeV gamma-ray pulsations from both were discovered using the Fermi Large Area Telescope. Twelve years of radio and polarization data allow detailed investigations. The two pulsars resemble each other both in radio and in gamma-ray data. Both are rare in having a single gamma-ray pulse offset far from the radio peak. The high dispersion measure for PSR J0248+6021 (DM = 370 pc cm^-3) is most likely due to its being within the dense, giant HII region W5 in the Perseus arm at a distance of 2 kpc, not beyond the edge of the Galaxy as obtained from models of average electron distributions. Its high transverse velocity and the low magnetic field along the line-of-sight favor this small distance. Neither gamma-ray, X-ray, nor optical data yield evidence for a pulsar wind nebula surrounding PSR J0248+6021. The gamma-ray luminosity for PSR J0248+6021 is L_ gamma = (1.4 \pm 0.3)\times 10^34 erg/s. For PSR J2240+5832, we find either L_gamma = (7.9 \pm 5.2) \times 10^34 erg/s if the pulsar is in the Outer arm, or L_gamma = (2.2 \pm 1.7) \times 10^34 erg/s for the Perseus arm. These luminosities are consistent with an L_gamma ~ sqrt(Edot) rule. Comparison of the gamma-ray pulse profiles with model predictions, including the constraints obtained from radio polarization data, favor emission in the far magnetosphere. These two pulsars differ mainly in their inclination angles and acceleration gap widths, which in turn explains the observed differences in the gamma-ray peak widths.

Serendipitously we have discovered a rare, bright submillimeter galaxy (SMG) with a flux density of 30 +/- 2 mJy at lambda=1.2mm, using MAMBO2 at the IRAM 30-meter millimeter telescope. Although no optical counterpart is known for MM18423+5938, we were able to measure the redshift z=3.92960 +/- 0.00013 from the detection of CO lines using the IRAM Eight MIxer Receiver (EMIR). In addition, by collecting all available photometric data in the far-infrared and radio to constrain its spectral energy distribution, we derive the FIR luminosity 4.8 10^14/m Lsol and mass 6.0 10^9/m Msol for its dust, allowing for a magnification factor m caused by a probable gravitational lens. The corresponding star-formation rate is 8.3 10^4/m Msol/yr. The detection of three lines of the CO rotational ladder, and a significant upper limit for a fourth CO line, allow us to estimate an H2 mass of between 1.9 10^11/m Msol and 1.1 10^12/m Msol. The two lines CI(3p1-3p0) and CI(3p2-3p1) were clearly detected and yield a [CI]/[H2] number abundance between 1.4 10^-5 and 8.0 10^-5. Upper limits are presented for emission lines of HCN, HCO^+, HNC, H_2O and other molecules observed. The moderate excitation of the CO lines is indicative of an extended starburst, and excludes the dominance of an AGN in heating this high-redshift SMG.

We present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

Although 70 % of the stars in the Galaxy are M-dwarfs, thermal emission searches for cold debris disks have been conducted mostly for A-type and solar-type stars. We report on new lambda=1.2 mm continuum observations of thirty M-dwarfs, using the MAMBO-2 bolometer array camera at the IRAM 30m telescope. For a statistical analysis, we combine these data with our prior SCUBA and MAMBO-2 observations of 20 other M-dwarfs. Our total sample divides in M-dwarfs in moving groups, with relatively young ages, and in nearby M-dwarfs with unknown ages. Only one cold debris disk (GJ842.2) was detected significantly. We compare the implied disk abundance constraints with those found in two comparable submillimeter surveys of 10 to 190 Myr old A- and FGK-type stars. For the 19 youngest (ages less than 200 Myr) M-dwarfs in our sample, we derive a cold disk fraction of 5.3^{+10.5}_{-5.0} %, compared to 15 +/-11.5 % for FGK-stars and 22^{+33}_{-20} % for A-stars. Hence, for this age group, there is an apparent trend of fewer cold disks for later stellar types, i.e., lower star masses. Although its statistical significance is marginal, this trend is strengthened by the deeper sensitivity of observations in the M-dwarf sample. We derive a cold disk fraction of < 10 % for the older (likely a few Gyr) M-dwarfs in our sample. Finally, although inconclusively related to a debris disk, we present the complex millimeter structure found around the position of the M1.5 dwarf GJ526 in our sample.

The rare and conspicuous flux density variations of some radio sources (extragalactic and pulsars) for periods of weeks to months have been denoted Extreme Scattering Events (ESE's) by Fiedler et al. (1987). Presently, there is no astrophysical mechanism that satisfactorily produces this phenomenon. In this paper, we conjecture that inhomogeneities of the electronic density in the turbulent interstellar medium might be the origin of this phenomenon. We have tested this conjecture by a simulation of the scintillation of the pulsar B1937+21 at 1.4 GHz and 1.7 GHz for a period of six months. To this end, we have constructed a large square Kolmogorov phase screen made of 131k x 131k pixels with electron inhomogeneity scales ranging from 6 x $10^6$ m to $10^{12}$ m and used the Kirchhoff-Fresnel integral to simulate dynamic spectra of a pulsar within the framework of Physical Optics. The simulated light curves exhibit a 10 day long variation simultaneously at 1.41 and 1.7 GHz that is alike the ``ESE'' observed with the Nancay radiotelescope toward the pulsar B1937+21 in October 1989. Consequently, we conclude that ``ESE'' toward pulsars can be caused naturally by the turbulence in the ionized interstellar medium instead of invoking the crossing of discrete over pressured ionized clouds on the line of sight as in the model of Fiedler et al. (1987). We suggest that longer events could occur in a simulation of scintillation, if larger electron inhomogeneities > $10^{12}$ m were included in the construction of the Kolmogorov phase screen. This next step requires a supercomputer.

Debris disks are believed to be related to planetesimals left over around stars after planet formation has ceased. The frequency of debris disks around M-dwarfs which account for 70% of the stars in the Galaxy is unknown while constrains have already been found for A- to K-type stars. We have searched for cold debris disks around 32 field M-dwarfs by conducting observations at lambda = 850 microns with the SCUBA bolometerarray camera at the JCMT and at lambda = 1.2mm with the MAMBO array at the IRAM 30-m telescopes. This is the first survey of a large sample of M-dwarfs conducted to provide statistical constraints on debris disks around this type of stars. We have detected a new debris disk around the M0.5 dwarf GJ842.2 at lambda = 850 microns, providing evidence for cold dust at large distance from this star (~ 300AU). By combining the results of our survey with the ones of Liu et al. (2004), we estimate for the first time the detection rate of cold debris disks around field M-dwarfs with ages between 20 and 200 Myr. This detection rate is 13^{+6}_{-8} % and is consistent with the detection rate of cold debris disks (9 - 23 %) around A- to K-type main sequence stars of the same age. This is an indication that cold disks may be equally prevalent across stellar spectral types.

Astrometric accuracy of complex modern VLBI arrays cannot be calculated analytically. We study the astrometric accuracy of phase-referenced VLBI observations for the VLBA, EVN and global VLBI array by simulating VLBI data for targets at declinations -25$^\circ$, 0$^\circ$, 25$^\circ$, 50$^\circ$, 75$^\circ$ and 85$^\circ$. The systematic error components considered in this study are calibrator position, station coordinate, Earth orientation and troposphere parameter uncertainties. We provide complete tables of the astrometric accuracies of these arrays for a source separation of 1$^\circ$ either along the right ascension axis or along the declination axis. Astrometric accuracy is 50microas at mid declination and is 300microas at low (-25$^\circ$) and high (85$^\circ$) declinations for the VLBA and EVN. In extending our simulations to source separations of 0.5$^\circ$ and 2$^\circ$, we establish the formula for the astrometric accuracy of the VLBA: Delta = (Delta_1$^\circ$-14)*d+ 14 (microas) where Delta_1$^\circ$ is the astrometric accuracy for a separation d=1$^\circ$ provided in our tables for various declinations and conditions of the wet troposphere. We argue that this formula is also valid for the astrometric accuracy of the EVN and global VLBI array.

We present an astrometric analysis of the binary systems ABDorA /ABDorC and ABDorBa / ABDorBb. These two systems of well-known late-type stars are gravitationally associated and they constitute the quadruple ABDoradus system. From the astrometric data available at different wavelengths, we report: (i) a determination of the orbit of ABDorC, the very low mass companion to ABDorA, which confirms the mass estimate of 0.090Msun reported in previous works; (ii) a measurement of the parallax of ABDorBa, which unambiguously confirms the long-suspected physical association between this star and ABDorA; and (iii) evidence of orbital motion of ABDorBa around ABDorA, which places an upper bound of 0.4Msun on the mass of the pair ABDorBa / ABDorBb (50% probability). Further astrometric monitoring of the system at all possible wavelengths would determine with extraordinary precision the dynamical mass of its four components.

The most complete catalogue of nearby stars (Gliese and Jarheiss, 1991) has been used to search for the stars whose photospheric thermal emission can be detected by the future millimeter array ALMA. We found that 446 nearby stars with spectral types ranging from A to M are detectable. A long-term astrometric observation programme (> 10 yr) of these stars with ALMA (theoretical astrometric precision = 0.1 milliarcsecond) would be sensitive to wobbles caused by the gravitational pull of possible unseen planets. Such a programme would probe minimum planetary masses as low as 0.1 Jupiter for an orbital period of 10 years. We provide histogrammes of these minimum masses and of ALMA integration times for these 446 ALMA stars.

Spatial phase-referencing in dual-field optical interferometry is reconsidered. Our analysis is based on the 2-sample variance of the differential phase between target and reference star. We show that averaging over time of the atmospheric effects depends on this 2-sample phase variance (Allan variance) rather than on the true variance. The proper expression for fringe smearing beyond the isoplanatic angle is derived. With simulations of atmospheric effects, based on a Paranal turbulence model, we show how the performances of a dual-field optical interferometer can be evaluated in a diagram 'separation angle' versus 'magnitude of faint object'. In this diagram, a domain with short exposure is found to be most useful for interferometry, with about the same magnitude limits in the H and K bands. With star counts from a Galaxy model, we evaluate the sky coverage for differential astrometry and detection of exoplanets, i.e. likelihood of faint reference stars in the vicinity of a bright target. With the 2mass survey, we evaluate sky coverage for phase-referencing, i.e. avaibility of a bright enough star for main delay tracking in the vicinity of any target direction.

We present VLBI images of the RS CVn binary star HR 5110 (=BH CVn; HD 118216), obtained from observations made at 8.4 GHz on 1994 May 29/30 in support of the NASA/Stanford Gravity Probe B project. Our images show an emission region with a core-halo morphology. The core was 0.39 +/- 0.09 mas (FWHM) in size, or 66% +/- 20% of the 0.6 +/- 0.1 mas diameter of the chromospherically active K subgiant star in the binary system. The halo was 1.95 +/- 0.22 mas (FWHM) in size, or 1.8 +/- 0.2 times the 1.1 +/- 0.1 mas separation of the centers of the K and F stars. The core increased significantly in brightness over the course of the observations and seems to have been the site of flare activity that generated an increase in the total flux density of ~200% in 12 hours. The fractional circular polarization simultaneously decreased from ~10% to 2.5%.

We present a spectral analysis of 35 GRBs detected with the HETE-2 gamma-ray detectors (the FREGATE instrument) in the energy range 7-400 keV. The GRB sample analyzed is made of GRBs localized with the Wide Field X-ray Monitor onboard HETE-2 or with the GRB Interplanetary Network. We derive the spectral parameters of the time-integrated spectra, and present the distribution of the low-energy photon index, alpha, and of the peak energy, e_peak . We then discuss the existence and nature of the recently discovered X-Ray Flashes and their relationship with classical GRBs.

We present 5 GHz space-VLBI observations of the quasar pair B1342+662/B1342+663 that demonstrate the feasibility of phase-reference techniques using an antenna in space. The space-based data were taken by the satellite HALCA, of the space-VLBI mission VSOP. From residual (referenced) phases we derive an upper bound of 10 meters to the uncertainty of the spacecraft orbit reconstruction. An analysis of the phase-reference maps of the sources additionally suggests that the above mentioned uncertaintyis likely not larger than 3 meters. With errors of this magnitude, HALCA is a useful tool for astrometric studies of close pairs of radio sources.

The Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (CGRO) has triggered on 1637 cosmic gamma-ray bursts between 1991 April 19 and 1996 August 29. These events constitute the Fourth BATSE burst catalog. The current version (4Br) has been revised from the version first circulated on CD-ROM in September 1997 (4B) to include improved locations for a subset of bursts that have been reprocssed using additional data. A significant difference from previous BATSE catalogs is the inclusion of bursts from periods when the trigger energy range differed from the nominal 50-300 keV. We present tables of the burst occurrence times, locations, peak fluxes, fluences, and durations. In general, results from previous BATSE catalogs are confirmed here with greater statistical significance.

We investigate spectral evolution in 37 bright, long gamma-ray bursts observed with the BATSE Spectroscopy Detectors. High resolution spectra are characterized by the energy of the peak of \nfn~and the evolution of this quantity is examined relative to the emission intensity. In most cases it is found that this peak energy either rises with or slightly precedes major intensity increases and softens for the remainder of the pulse. Inter-pulse emission is generally harder early in the burst. For bursts with multiple intensity pulses, later spikes tend to be softer than earlier ones indicating that the energy of the peak of \nfn~is bounded by an envelope which decays with time. Evidence is found that bursts in which the bulk of the flux comes well after the event which triggers the instrument tend to show less peak energy variability and are not as hard as several bursts in which the emission occurs promptly after the trigger. Several recently proposed burst models are examined in light of these results and no qualitative conflicts with the observations presented here are found.