We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $\mu$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

The relative positions of the high and low surface brightness regions of radio-loud active galaxies in the 3CR sample were found by Fanaroff and Riley to be correlated with their luminosity. We revisit this canonical relationship with a sample of 5805 extended radio-loud AGN from the LOFAR Two-Metre Sky Survey (LoTSS), compiling the most complete dataset of radio-galaxy morphological information obtained to date. We demonstrate that, for this sample, radio luminosity does *not* reliably predict whether a source is edge-brightened (FRII) or centre-brightened (FRI). We highlight a large population of low-luminosity FRIIs, extending three orders of magnitude below the traditional FR break, and demonstrate that their host galaxies are on average systematically fainter than those of high-luminosity FRIIs and of FRIs matched in luminosity. This result supports the jet power/environment paradigm for the FR break: low-power jets may remain undisrupted and form hotspots in lower mass hosts. We also find substantial populations that appear physically distinct from the traditional FR classes, including candidate restarting sources and ``hybrids''. We identify 459 bent-tailed sources, which we find to have a significantly higher SDSS cluster association fraction (at $z<0.4$) than the general radio-galaxy population, similar to the results of previous work. The complexity of the LoTSS faint, extended radio sources demonstrates the need for caution in the automated classification and interpretation of extended sources in modern radio surveys, but also reveals the wealth of morphological information such surveys will provide and its value for advancing our physical understanding of radio-loud AGN.

We present the first systematic study of the radio continuum size evolution of star-forming galaxies (SFGs) over the redshift range $0.35<z<2.25$. We use the VLA COSMOS 3GHz map (noise $\rm rms=2.3\,\mu Jy \,beam^{-1}$, $\theta_{\rm beam}=0.75\,\rm arcsec$) to construct a mass-complete sample of 3184 radio-selected SFGs that reside on and above the main-sequence (MS) of SFGs. We find no clear dependence between the radio size and stellar mass, $M_{\star}$, of SFGs with $10.5\lesssim\log(M_\star/\rm M_\odot)\lesssim11.5$. Our analysis suggests that MS galaxies are preferentially extended, while SFGs above the MS are always compact. The median effective radius of SFGs on (above) the MS of $R_{\rm eff}=1.5\pm0.2$ ($1.0\pm0.2$) kpc remains nearly constant with cosmic time; a parametrization of the form $R_{\rm eff}\propto(1+z)^\alpha$ yields a shallow slope of only $\alpha=-0.26\pm0.08\,(0.12\pm0.14)$ for SFGs on (above) the MS. The size of the stellar component of galaxies is larger than the extent of the radio continuum emission by a factor $\sim$2 (1.3) at $z=0.5\,(2)$, indicating star formation is enhanced at small radii. The galactic-averaged star formation rate surface density $(\Sigma_{\rm SFR})$ scales with the distance to the MS, except for a fraction of MS galaxies ($\lesssim10\%$) that harbor starburst-like $\Sigma_{\rm SFR}$. These "hidden" starbursts might have experienced a compaction phase due to disk instability and/or merger-driven burst of star formation, which may or may not significantly offset a galaxy from the MS. We thus propose to jointly use $\Sigma_{\rm SFR}$ and distance to the MS to better identify the galaxy population undergoing a starbursting phase.

In this data paper we present and characterise the multi-component radio sources identified in the VLA-COSMOS Large Project at 3 GHz (0.75 arcsec resolution, 2.3 {\mu}Jy/beam rms), i.e. the radio sources which are composed of two or more radio blobs.The classification of objects into multi-components was done by visual inspection of 351 of the brightest and most extended blobs from a sample of 10,899 blobs identified by the automatic code blobcat. For that purpose we used multi-wavelength information of the field, such as the 1.4-GHz VLA-COSMOS data and the UltraVISTA stacked mosaic available for COSMOS. We have identified 67 multi-component radio sources at 3 GHz: 58 sources with AGN powered radio emission and 9 star-forming galaxies. We report 8 new detections that were not observed by the VLA-COSMOS Large Project at 1.4 GHz, due to the slightly larger area coverage at 3 GHz. The increased spatial resolution of 0.75 arcsec has allowed us to resolve (and isolate) multiple emission peaks of 28 extended radio sources not identified in the 1.4-GHz VLA-COSMOS map. We report the multi-frequency flux densities (324 MHz, 325 MHz, 1.4 GHz & 3 GHz), star-formation-rates, and stellar masses of these objects. Multi-component objects at 3-GHz VLA-COSMOS inhabit mainly massive galaxies (>10^10.5 Msun). The majority of the multi-component AGN lie below the main-sequence of star-forming galaxies (SFGs), in the green valley and the quiescent region. We provide detailed description of the objects: amongst the AGN there are 2 head-tail, 10 core-lobe, 9 wide-angle-tail (WAT), 8 double-double or Z-/X-shaped, 3 bent-tail radio sources, and 26 symmetric sources, while amongst the SFGs we find the only star-forming ring seen in radio emission in COSMOS. We report a large number (32/58) of disturbed/bent multi-component AGN, 18 of which do not lie within X-ray groups in COSMOS (0.08 < z < 1.53). [abridged]

We construct the average radio spectral energy distribution (SED) of highly star-forming galaxies (HSFGs) up to z~4. Infrared and radio luminosities are bound by a tight correlation that is defined by the so-called q parameter. This infrared-radio correlation provides the basis for the use of radio luminosity as a star-formation tracer. Recent stacking and survival analysis studies find q to be decreasing with increasing redshift. It was pointed out that a possible cause of the redshift trend could be the computation of rest-frame radio luminosity via a single power-law assumption of the star-forming galaxies' (SFGs) SED.To test this, we constrained the shape of the radio SED of a sample of HSFGs. To achieve a broad rest-frame frequency range, we combined previously published VLA observations of the COSMOS field at 1.4 GHz and 3 GHz with unpublished GMRT observations at 325 MHz and 610 MHz by employing survival analysis to account for non-detections in the GMRT maps. We selected a sample of HSFGs in a broad redshift range (0.3<z<4,SFR>100M0/yr) and constructed the average radio SED. By fitting a broken power-law, we find that the spectral index changes from $\alpha_1=0.42\pm0.06$ below a rest-frame frequency of 4.3 GHz to $\alpha_2=0.94\pm0.06$ above 4.3 GHz. Our results are in line with previous low-redshift studies of HSFGs (SFR>10M0/yr) that show the SED of HSFGs to differ from the SED found for normal SFGs (SFR<10M0/yr). The difference is mainly in a steeper spectrum around 10 GHz, which could indicate a smaller fraction of thermal free-free emission. Finally, we also discuss the impact of applying this broken power-law SED in place of a simple power-law in K-corrections of HSFGs and a typical radio SED for normal SFGs drawn from the literature. We find that the shape of the radio SED is unlikely to be the root cause of the q-z trend in SFGs.

We investigate the linear radio size properties of the $\mu$Jy populations of radio-selected active galactic nuclei (AGN) and star-forming galaxies (SFGs) using a multi-resolution catalog based on the original VLA-COSMOS 3\,GHz 0\farcs75 resolution mosaic and its convolved images (up to a resolution of 2\farcs2). The final catalog contains 6\,399 radio sources above a 3\,GHz total flux density of $S_T>20$ $\mu$Jy (median $<S_T>=37$ $\mu$Jy), with redshift information (median $<z>=1.0$), and multi-wavelength classification as SFGs, radio-excess AGN (RX-AGN), or non-radio-excess AGN (NRX-AGN). RX-AGN are those whose radio emission exceeds the star formation rate derived by fitting the global spectral energy distribution. We derive the evolution with redshift and luminosity of the median linear sizes of each class of objects. We find that RX-AGN are compact, with median sizes of $\sim$ 1-2 kpc and increasing with redshift, corresponding to an almost constant angular size of 0\farcs25. NRX-AGN typically have radio sizes a factor of 2 larger than the RX-AGN. The median radio size of SFGs is about 5 kpc up to $z\sim 0.7$, and it decreases beyond this redshift. Using luminosity-complete subsamples of objects, we separately investigate the effect of redshift and luminosity dependance. We compare the radio sizes of SFGs with those derived in the rest-frame far-infrared (FIR) and UV bands. We find that SFGs have comparable sizes (within 15\%) in the radio and rest-frame FIR, while the sizes measured in the UV-band are systematically larger than the radio sizes.

Exploring how radio-emitting active galactic nuclei (AGN) behave and evolve with time is critical for understanding how AGN feedback impacts galaxy evolution. In this work, we investigate the relationship between 1.4 GHz radio continuum AGN luminosity ($L^{\rm AGN}_{\rm 1.4}$), specific black hole accretion rate (s-BHAR, defined as the accretion luminosity relative to the galaxy stellar mass) and redshift, for a luminosity-complete sample of radio-selected AGN in the VLA COSMOS 3 GHz Large Project. The sample was originally selected from radio-continuum observations at 3 GHz, and includes about 1800 radio AGN identified via ($>2\sigma$) radio-excess relative to the infrared-radio correlation of star-forming galaxies. We further select a subsample of over 1200 radio AGN that is complete in $L^{\rm AGN}_{\rm 1.4}$ over different redshift ranges, out to z~4, and use X-ray stacking to calculate the average s-BHAR in each $L^{\rm AGN}_{\rm 1.4}$-$z$ bin. We find that the average s-BHAR is independent of $L^{\rm AGN}_{\rm 1.4}$, at all redshifts. However, we see a strong increase of s-BHAR with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. This trend resembles the strong increase in the fraction of star-forming host galaxies (based on the $(NUV-r)$ / $(r-J)$ colours) with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. A possible explanation for this similarity might imply a link between average AGN radiative power and availability of cold gas supply within the host galaxy. This study corroborates the idea that radio-selected AGN become more radiatively efficient towards earlier epochs, independently of their radio power.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

We present the 25 square-degree GMRT-XXL-N 610 MHz radio continuum survey, conducted at 50~cm wavelength with the Giant Metrewave Radio Telescope (GMRT) towards the XXL Northern field (XXL-N). We combined previously published observations of the XMM-Large Scale Structure (XMM-LSS) field, located in the central part of XXL-N, with newly conducted observations towards the remaining XXL-N area, and imaged the combined data-set using the Source Peeling and Atmospheric Modeling ({\sc SPAM}) pipeline. The final mosaic encompasses a total area of $30.4$ square degrees, with ${\rm rms}<150~\mu$Jy/beam over 60\% of the area. The ${\rm rms}$ achieved in the inner 9.6 square degree area, enclosing the XMM-LSS field, is about $200~\mu$Jy/beam, while that over the outer 12.66 square degree area (which excludes the noisy edges) is about $45~\mu$Jy/beam. The resolution of the final mosaic is 6.5 arcsec. We present a catalogue of 5\,434 sources detected at $\geq7\times {\rm rms}$. We verify, and correct the reliability of, the catalog in terms of astrometry, flux, and false detection rate. Making use of the (to date) deepest radio continuum survey over a relatively large (2 square degree) field, complete at the flux levels probed by the GMRT-XXL-N survey, we also assess the survey's incompleteness as a function of flux density. The radio continuum sensitivity reached over a large field with a wealth of multi-wavelength data available makes the GMRT-XXL-N 610 MHz survey an important asset for studying the physical properties, environments and cosmic evolution of radio sources, in particular radio-selected active galactic nuclei (AGN).

We present an investigation of radio luminosity functions (LFs) and number counts based on the Karl G. Jansky Very Large Array-COSMOS 3 GHz Large Project. The radio-selected sample of 7826 galaxies with robust optical/near-infrared counterparts with excellent photometric coverage allows us to construct the total radio LF since z~5.7. Using the Markov chain Monte Carlo algorithm, we fit the redshift dependent pure luminosity evolution model to the data and compare it with previously published VLA-COSMOS LFs obtained on individual populations of radio-selected star-forming galaxies and galaxies hosting active galactic nuclei classified on the basis of presence or absence of a radio excess with respect to the star-formation rates derived from the infrared emission. We find they are in excellent agreement, thus showing the reliability of the radio excess method in selecting these two galaxy populations at radio wavelengths. We study radio number counts down to submicrojansky levels drawn from different models of evolving LFs. We show that our evolving LFs are able to reproduce the observed radio sky brightness, even though we rely on extrapolations toward the faint end. Our results also imply that no new radio-emitting galaxy population is present below 1 ujy. Our work suggests that selecting galaxies with radio flux densities between 0.1 and 10 ujy will yield a star-forming galaxy in 90-95 % of the cases with a high percentage of these galaxies existing around a redshift of z~2, thus providing useful constraints for planned surveys with the Square Kilometer Array and its precursors.

Using infrared data from the Herschel Space Observatory and Karl G. Jansky Very Large Array (VLA) 3 GHz observations in the COSMOS field, we investigate the redshift evolution of the infrared-radio correlation (IRRC) for star-forming galaxies (SFGs) we classify as either spheroid- or disc-dominated based on their morphology. The sample predominantly consists of disc galaxies with stellar mass ${\gtrsim}10^{10}\,M_{\odot}$, and residing on the star-forming main sequence (MS). After the removal of AGN using standard approaches, we observe a significant difference between the redshift-evolution of the median IR/radio ratio $\overline{q}_{\mathrm{TIR}}$ of (i) a sample of ellipticals, plus discs with a substantial bulge component (`spheroid-dominated' SFGs) and, (ii) virtually pure discs and irregular systems (`disc-dominated' SFGs). The spheroid-dominated population follows a declining $\overline{q}_{\mathrm{TIR}}$ vs. $z$ trend similar to that measured in recent evolutionary studies of the IRRC. However, for disc-dominated galaxies, where radio and IR emission should be linked to star formation in the most straightforward way, we measure very little change in $\overline{q}_{\mathrm{TIR}}$. This suggests that low-redshift calibrations of radio emission as an SFR-tracer may remain valid out to at least $z\,{\simeq}\,1\,{-}\,1.5$ for pure star-forming systems. We find that the different redshift-evolution of $q_{\rm TIR}$ for the spheroid- and disc-dominated sample is mainly due to an increasing radio excess for spheroid-dominated galaxies at $z\,{\gtrsim}\,$0.8, hinting at some residual AGN activity in these systems. This finding demonstrates that in the absence of AGN the IRRC is independent of redshift, and that radio observations can therefore be used to estimate SFRs at all redshifts for genuinely star-forming galaxies.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

Based on a sample of over 1,800 radio AGN at redshifts out to z~5, which have typical stellar masses within ~3x(10^{10}-10^{11}) Msol, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz radio luminosity functions for radio AGN (L_1.4GHz ~ 10^{22}-10^{27} W/Hz) out to z~5. We constrained the evolution of this population via continuous models of pure density and pure luminosity evolutions, and we found best-fit parametrizations of Phi*~(1+z)^{(2.00+/-0.18)-(0.60+/-0.14)z}, and L*~(1+z)^{(2.88+/-0.82)-(0.84+/-0.34)z}, respectively, with a turnover in number and luminosity densities of the population at z~1.5. We converted 1.4 GHz luminosity to kinetic luminosity taking uncertainties of the scaling relation used into account. We thereby derived the cosmic evolution of the kinetic luminosity density provided by the AGN and compared this luminosity density to the radio-mode AGN feedback assumed in the Semi-Analytic Galaxy Evolution (SAGE) model, i.e., to the redshift evolution of the central supermassive black hole accretion luminosity taken in the model as the source of heating that offsets the energy losses of the cooling, hot halo gas, and thereby limits further stellar mass growth of massive galaxies. We find that the kinetic luminosity exerted by our radio AGN may be high enough to balance the radiative cooling of the hot gas at each cosmic epoch since z~5. However, although our findings support the idea of radio-mode AGN feedback as a cosmologically relevant process in massive galaxy formation, many simplifications in both the observational and semi-analytic approaches still remain and need to be resolved before robust conclusions can be reached.

We present new observations of the GOODS-N field obtained at 5.5 GHz with the Karl G. Jansky Very Large Array (VLA). The central region of the field was imaged to a median r.m.s. of 3 microJy/beam with a resolution of 0.5 arcsec. From a 14-arcmin diameter region we extracted a sample of 94 radio sources with signal-to-noise ratio greater than 5. Near-IR identifications are available for about 88 percent of the radio sources. We used different multi-band diagnostics to separate active galactic nuclei (AGN), both radiatively efficient and inefficient, from star-forming galaxies. From our analysis, we find that about 80 percent of our radio-selected sample is AGN-dominated, with the fraction raising to 92 percent when considering only the radio sources with redshift >1.5. This large fraction of AGN-dominated radio sources at very low flux densities (the median flux density at 5.5 GHz is 42 microJy), where star-forming galaxies are expected to dominate, is somewhat surprising and at odds with other results. Our interpretation is that both the frequency and angular resolution of our radio observations strongly select against radio sources whose brightness distribution is diffuse on scale of several kpc. Indeed, we find that the median angular sizes of the AGN-dominated sources is around 0.2-0.3 arcsec against 0.8 arcsec for star-forming galaxies. This highlights the key role that high frequency radio observations can play in pinpointing AGN-driven radio emission at microJy levels. This work is part of the eMERGE legacy project.

We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.

(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.

We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies.

The sources SDSS\,J113126.08-020459.2 and SDSS\,J110851.04+065901.4 are two double-peaked [O III] emitting AGN, identified as candidate binary AGNs by optical and near infrared (NIR) observations. We observed the two sources with high resolution Very Long Baseline Interferometry (VLBI) using the European VLBI Network at 5\,GHz, reduced VLA observations at three frequencies available for one of the sources, and used archival HST observations. For the source SDSS\,J113126.08-020459.2, the VLBI observations detected only one single compact component associated with the eastern NIR nucleus. In SDSS\,J110851.04+065901.4, the VLBI observations did not detect any compact components, but the VLA observations allowed us to identify a possible compact core in the region of the north-western optical/NIR nucleus. In this source we find kpc-scale extended radio emission that is spatially coincident to the ultraviolet continuum and to the extended emission narrow line region. The UV continuum is significantly obscured since the amount of extended radio emission yields a star formation rate of about 110 \msunyr, which is an order of magnitude larger than implied by the observed ultraviolet emission. Our analysis confirms the presence of only one AGN in the two candidate binary AGNs.

We present 2.1 GHz imaging with the Australia Telescope Compact Array (ATCA) of a 6.5 deg^2 region within the XXM-Newton XXL South field using a band of 1.1-3.1 GHz. We achieve an angular resolution of 4.7" x 4.2" in the final radio continuum map with a median rms noise level of 50 uJy/beam. We identify 1389 radio sources in the field with peak S/N >=5 and present the catalogue of observed parameters. We find that 305 sources are resolved, of which 77 consist of multiple radio components. These number counts are in agreement with those found for the COSMOS-VLA 1.4 GHz survey. We derive spectral indices by a comparison with the Sydney University Molongolo Sky Survey (SUMSS) 843MHz data. We find an average spectral index of -0.78 and a scatter of 0.28, in line with expectations. This pilot survey was conducted in preparation for a larger ATCA program to observe the full 25 deg^2 southern XXL field. When complete, the survey will provide a unique resource of sensitive, wide-field radio continuum imaging with complementary X-ray data in the field. This will facilitate studies of the physical mechanisms of radio-loud and radio-quiet AGNs and galaxy clusters, and the role they play in galaxy evolution. The source catalogue is publicly available online via the XXL Master Catalogue browser and the Centre de Donn\'ees astronomiques de Strasbourg (CDS).

Radio emission at cm wavelengths from highly star-forming galaxies, such as SMGs, is dominated by synchrotron radiation arising from supernova activity. Using deep, high-resolution ($1\sigma=2.3$ $\mu$Jy beam$^{-1}$; $0.75^{"}$) cm radio-continuum observations taken by the VLA-COSMOS 3 GHz Large Project, we studied the radio-emitting sizes of a flux-limited sample of SMGs in the COSMOS field. Of the 39 SMGs studied here, 3 GHz emission was detected towards 18 of them ($\sim46\pm11\%$) with S/N ratios in the range of ${\rm S/N=4.2-37.4}$. Using 2D elliptical Gaussian fits, we derived a median deconvolved major axis FWHM size of $0.54^{"}\pm 0.11^{"}$ for our 18 SMGs detected at 3 GHz. For the 15 SMGs with known redshift we derived a median linear major axis FWHM of $4.2\pm0.9$ kpc. No clear correlation was found between the radio-emitting size and the 3 GHz or submm flux density, or the redshift of the SMG. However, there is a hint of larger radio sizes at $z\sim2.5-5$ compared to lower redshifts. The sizes we derived are consistent with previous SMG sizes measured at 1.4 GHz and in mid-$J$ CO emission, but significantly larger than those seen in the (sub)mm continuum emission. One possible scenario is that SMGs have i) an extended gas component with a low dust temperature, and which can be traced by low- to mid-$J$ CO line emission and radio continuum emission, and ii) a warmer, compact starburst region giving rise to the high-$J$ line emission of CO, which could dominate the dust continuum size measurements. Because of the rapid cooling of CR electrons in dense starburst galaxies ($\sim10^4-10^5$ yr), the more extended synchrotron radio-emitting size being a result of CR diffusion seems unlikely. Instead, if SMGs are driven by galaxy mergers the radio synchrotron emission might arise from more extended magnetised ISM around the starburst region.

We started a follow-up investigation of the Deep X-ray Radio Blazar Survey objects with declination >-10 deg. We undertook a survey with the EVN at 5GHz to make the first images of a complete sample of weak blazars, aiming at a comparison between high- and low-power samples of blazars. All of the 87 sources observed were detected. Point-like sources are found in 39 cases, and 48 show core-jet structure. According to the spectral indices previously obtained, 58 sources show a flat spectral index, and 29 sources show a steep spectrum or a spectrum peaking at a frequency around 1-2 GHz. Adding to the DXRBS objects we observed those already observed with ATCA in the southern sky, we found that 14 blazars and a SSRQ, are associated to gamma-ray emitters. We found that 56 sources can be considered blazars. We also detected 2 flat spectrum NLRGs. About 50% of the blazars associated to a gamma-ray object are BL Lacs, confirming that they are more likely detected among blazars gamma-emitters. We confirm the correlation found between the source core flux density and the gamma-ray photon fluxes down to fainter flux densities. We also found that weak blazars are also weaker gamma-ray emitters compared to bright blazars. Twenty-two sources are SSRQs or CSSs, and 7 are GPSs. The available X-ray ROSAT observations allow us to suggest that CSS and GPS quasars are not obscured by large column of cold gas surrounding the nuclei. We did not find any significant difference in X-ray luminosity between CSS and GPS quasars.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies

We present deep 3 GHz VLA observations of the potentially recoiling black hole CID-42 in the COSMOS field. This galaxy shows two optical nuclei in the HST/ACS image and a large velocity offset of ~ 1300 km/s between the broad and narrow H beta emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data has a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multi-scale multi-frequency synthesis and spectral windows stacking. The final resolutions and sensitivities of these maps are 0.7" with rms = 4.6 muJy/beam and 0.9" with rms = 4.8 muJy/beam respectively. With a 7 sigma detection we find that the entire observed 3 GHz radio emission can be associated with the South-Eastern component of CID-42, coincident with the detected X-ray emission. We use our 3 GHz data combined with other radio data from the literature ranging from 320 MHz to 9 GHz, which include the VLA, VLBA and GMRT data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum AGN in the South-Eastern component which may be surrounded by a more extended region of old synchrotron electron population or shocks generated by the outflow from the supermassive black hole. Our data are consistent with the recoiling black hole picture but cannot rule out the presence of an obscured and radio-quiet SMBH in the North-Western component.

We present 90 cm VLA imaging of the COSMOS field, comprising a circular area of 3.14 square degrees at 8.0"x6.0" angular resolution with an average rms of 0.5 mJy/beam. The extracted catalog contains 182 sources (down to 5.5sigma), 30 of which are multi-component sources. Using Monte Carlo artificial source simulations we derive the completeness of the catalog, and we show that our 90 cm source counts agree very well with those from previous studies. Using X-ray, NUV-NIR and radio COSMOS data to investigate the population mix of our 90 cm radio sample, we find that our sample is dominated by active galactic nuclei (AGN). The average 90-20 cm spectral index (S_nu~nu**alpha, where S_nu is the flux density at frequency nu, and alpha the spectral index) of our 90 cm selected sources is -0.70, with an interquartile range of -0.90 to -0.53. Only a few ultra-steep-spectrum sources are present in our sample, consistent with results in the literature for similar fields. Our data do not show clear steepening of the spectral index with redshift. Nevertheless, our sample suggests that sources with spectral indices steeper than -1 all lie at z>1, in agreement with the idea that ultra-steep-spectrum radio sources may trace intermediate-redshift galaxies (z>1).

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

We describe the completed VIMOS VLT Deep Survey, and the final data release of 35016 galaxies and type-I AGN with measured spectroscopic redshifts up to redshift z~6.7, in areas 0.142 to 8.7 square degrees, and volumes from 0.5x10^6 to 2x10^7h^-3Mpc^3. We have selected samples of galaxies based solely on their i-band magnitude reaching i_{AB}=24.75. Spectra have been obtained with VIMOS on the ESO-VLT, integrating 0.75h, 4.5h and 18h for the Wide, Deep, and Ultra-Deep nested surveys. A total of 1263 galaxies have been re-observed independently within the VVDS, and from the VIPERS and MASSIV surveys. They are used to establish the redshift measurements reliability, to assess completeness, and to provide a weighting scheme taking into account the survey selection function. We describe the main properties of the VVDS samples, and the VVDS is compared to other spectroscopic surveys. In total we have obtained spectroscopic redshifts for 34594 galaxies, 422 type-I AGN, and 12430 Galactic stars. The survey has enabled to identify galaxies up to very high redshifts with 4669 redshifts in 1<=z_{spec}<=2, 561 in 2<=z_{spec}<=3 and 468 with z_{spec}>3, and specific populations like LAE have been identified out to z=6.62. We show that the VVDS occupies a unique place in the parameter space defined by area, depth, redshift coverage, and number of spectra. The VVDS provides a comprehensive survey of the distant universe, covering all epochs since z, or more than 12 Gyr of cosmic time, with a uniform selection, the largest such sample to date. A wealth of science results derived from the VVDS have shed new light on the evolution of galaxies and AGN, and their distribution in space, over this large cosmic time. A final public release of the complete VVDS spectroscopic redshift sample is available at http://cesam.lam.fr/vvds.

We present 5 GHz e-MERLIN observations of the GOODS-N region at sub-arcsec resolution (0.2--0.5 arcsec). These data form part of the early commissioning observations for the e-MERLIN interferometer and a pilot for the e-MERLIN legacy program eMERGE. A total of 17 sources were detected with S/N>3. These observations provide unique information on the radio source morphology at sub-arcsec scales. For twelve of these sources, deeper 1.4 GHz MERLIN+VLA observations at the same spatial resolution are available, allowing radio spectral indices to be derived for ten sources on sub-arcsec angular scales. Via analysis of the spectral indices and radio morphologies, these sources have been identified as AGN cores in moderate-to-high redshift (1<z<4) galaxies. These results have provided AGN (or AGN candidate) classification for six previously unclassified sources and confirmed the AGN nature of the rest of the sample. Ultimately the eMERGE project will image the GOODS-N region at 1.4 and 5 GHz with higher resolution (about 50 mas at 5 GHz) and down to sub-microJy sensitivities. The unique combination of sensitivity and spatial resolution will be exploited to study star formation and AGN activity in distant galaxies.

The nuclear region of the Luminous Infra-red Galaxy Arp 299-A hosts a recent ($\simeq 10$ Myr), intense burst of massive star formation which is expected to lead to numerous core-collapse supernovae (CCSNe). Previous VLBI observations, carried out with the EVN at 5.0 GHz and with the VLBA at 2.3 and 8.4 GHz, resulted in the detection of a large number of compact, bright, non-thermal sources in a region $\lsim$150 pc in size. We aim at establishing the nature of all non-thermal, compact components in Arp 299-A, as well as estimating its core-collapse supernova rate. We use multi-epoch European VLBI Network (EVN) observations taken at 5.0 GHz to image with milliarcsecond resolution the compact radio sources in the nuclear region of Arp 299-A. We also use one single-epoch 5.0 GHz Multi-Element Radio Linked Interferometer Network (MERLIN) observation to image the extended emission in which the compact radio sources --traced by our EVN observations-- are embedded. Twenty-six compact sources are detected, 8 of them are new objects not previously detected. The properties of all detected objects are consistent with them being a mixed population of CCSNe and SNRs. We find clear evidence for at least two new CCSNe, implying a lower limit to the CCSN rate of $\nu_{\rm SN}\gsim$0.80 SN/yr indicating that the bulk of the current star formation in Arp 299-A is taking place in the innermost $\sim 150$ pc. Our MERLIN observations trace a region of diffuse, extended emission which is co-spatial to the region where all compact sources are found. From this diffuse, non-thermal radio emission we obtain an independent estimate for the core-collapse supernova rate, which is in the range $\nu_{\rm SN}=0.40$ - 0.65 SN/yr, roughly in agreement with previous estimates and our direct estimate of the CCSN rate from the compact radio emission.

We performed observations with the Effelsberg 100-m radio telescope to measure flux densities and polarised emission of sources selected from the "Deep X-ray Radio Blazar Survey" (DXRBS) to better define their spectral index behaviour in the radio band, with the aim to construct a homogeneous sample of blazars. Sources were observed at four different frequencies with the Effelsberg 100-m telescope. We complemented these measurements with flux density data at 1.4GHz derived from the NRAO VLA Sky Survey.The spectral indices of a sample of faint blazars were computed making use of almost simultaneous measurements. Sixty-six percent of the sources can be classified as "bona fide" blazars. Seven objects show a clearly inverted spectral index. Seventeen sources previously classified as flat spectrum radio quasars (FSRQs) are actually steep spectrum radio quasars (SSRQs). The flux densities obtained with the Effelsberg 100-m telescope at 5GHz are compared with the flux densities listed in the Green Bank GB6 survey and in the Parkes-MIT-NRAO PMN catalogue. About 43% of the sources in our sample exhibit flux density variations on temporal scales of 19 or 22 years. We confirm that 75 out of 103 sources of the DXRBS are indeed FSRQs. Twenty-seven sources show a spectral index steeper than -0.5 and should be classified as SSRQs. Polarised emission was detected for 36 sources at 4.85GHz. The median value of the percentage of polarised emission is (5.8+-0.9$)%. Five sources show rotation measure (RM) values >200 rad m^-2.

We present deep, high-resolution imaging of the nearby Fanaroff-Riley Class I (FR I) radio galaxies NGC 193, B2 0206+35, B2 0755+37 and M 84 at frequencies of 4.9 and 1.4 GHz using new and archival multi-configuration observations from the Very Large Array. In addition, we describe lower-resolution observations of B2 0326+39 and a reanalysis of our published images of 3C 296. All of these radio galaxies show twin jets and well-defined lobes or bridges of emission, and we examine the common properties of this class of source. We show detailed images of total intensity, brightness gradient, spectral index, degree of polarization and projected magnetic-field direction. The jet bases are very similar to those in tailed twin-jet sources and show the characteristics of decelerating, relativistic flows. Except on one side of M 84, we find that the jets can be traced at least as far as the ends of the lobes, where they often form structures which we call "caps" with sharp outer brightness gradients. Continuing, but less well collimated flows back into the lobes from the caps can often be identified by their relatively flat spectral indices. The lobes in these radio galaxies are similar in morphology, spectral-index distribution and magnetic-field structure to those in more powerful (FR II) sources, but lack hot-spots or other evidence for strong shocks at the ends of the jets. M 84 may be an intermediate case between lobed and tailed sources, in which one jet does not reach the end of its lobe, but disrupts to form a "bubble".

We present partial results from our monitoring of the nuclear region of the starburst galaxy IC 694 (=Arp 299-A) at radio wavelengths, aimed at discovering recently exploded CCSNe, as well as to determine their rate of explosion, which carries crucial information on star formation rates and starburst scenarios at work. Two epochs of eEVN observations at 5.0 GHz, taken in 2008, revealed the presence of a rich cluster of compact radio emitting sources in the central 150 pc of the nuclear starburst in Arp 299A. The large brightness temperatures observed for the compact sources indicate a non-thermal origin for the observed radio emission, implying that most, if not all, of those sources were young radio supernovae (RSNe) and supernova remnants (SNRs). More recently, contemporaneous EVN observations at 1.7 and 5.0 GHz taken in 2009 have allowed us to shed light on the compact radio emission of the parsec-scale structure in the nucleus of Arp 299-A. Namely, our EVN observations have shown that one of the compact VLBI sources, A1, previously detected at 5.0 GHz, has a flat spectrum between 1.7 and 5.0 GHz and is the brightest source at both frequencies. The morphology, radio luminosity, spectral index and ratio of radio-to-X-ray emission of the A1-A5 region allowed us to identify A1-A5 with long-sought AGN in Arp 299-A. This finding may suggest that both starburst and AGN are frequently associated phenomena in mergers. Finally, we also note that component A0, identified as a young RSN, exploded at the mere distance of two parsecs from the putative AGN in Arp 299-A, which makes this supernova one of the closest to a central supermassive black hole ever detected.

Context: The dusty nuclear regions of luminous infra-red galaxies (LIRGs) are heated by either an intense burst of massive star formation, an active galactic nucleus (AGN), or a combination of both. Disentangling the contribution of each of those putative dust-heating agents is a challenging task, and direct imaging of the innermost few pc can only be accomplished at radio wavelengths, using very high-angular resolution observations. Aims: We observed the nucleus A of the interacting starburst galaxy Arp 299, using very long baseline interferometry (VLBI) radio observations at 1.7 and 5.0 GHz. Our aim was to characterize the compact sources in the innermost few pc region of Arp 299-A, as well as to detect recently exploded core-collapse supernovae. Methods: We used the European VLBI Network (EVN) to image the 1.7 and 5.0 GHz compact radio emission of the parsec-scale structure in the nucleus of Arp 299-A with milliarcsecond resolution. Results: Our EVN observations show that one of the compact VLBI sources, A1, previously detected at 5.0 GHz, has a flat spectrum between 1.7 and 5.0 GHz and is the brightest source at both frequencies. Our 1.7 GHz EVN image shows also diffuse, low-surface brightness emission extending westwards from A1 and displays a prominent core-jet structure. Conclusions: The morphology, radio luminosity, spectral index and ratio of radio-to-X-ray emission of the A1-A5 region is consistent with a low-luminosity AGN (LLAGN), and rules out the possibility that it is a chain of young radio supernovae (RSNe) and supernova remnants (SNRs). We therefore conclude that A1-A5 is the long-sought AGN in Arp 299-A. This finding may suggest that both starburst and AGN are frequently associated phenomena in mergers.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

[Abridged]We exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent 24, 15, 12 and 8micron monochromatic and bolometric IR luminosity functions (LF) over the 0<z<2.5 redshift range. Our analysis is based on the combination of data from deep Spitzer surveys in the VVDS-SWIRE and GOODS areas. To our limiting flux of S(24)=400microJy our derived sample in VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S(24)=80microJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degs. We obtain reliable optical identifications and redshifts, providing us a rich and robust dataset for our luminosity function determination. Based on the multi-wavelength information available, we constrain the LFs at 8, 12, 15 and 24micron. We also extrapolate total IR luminosities from our best-fit to the observed SEDs of each source, and use this to derive the bolometric LF and comoving volume emissivity up to z~2.5. In the 0<z<1 interval, the bolometric IR luminosity density evolves as (1+z)^3.8+/-0.4. Although more uncertain at higher-z, our results show a flattening of the IR luminosity density at z>1. The mean redshift of the peak in the source number density shifts with luminosity: the brighest IR galaxies appear to be forming stars earlier in cosmic time (z>1.5), while the less luminous ones keep doing it at more recent epochs (z~1 for L(IR)<10^11L_sun). Our results suggest a rapid increase of the galaxy IR comoving volume emissivity back to z~1 and a constant average emissivity at z>1. We also seem to find a difference in the evolution rate of the source number densities as a function of luminosity, a downsizing evolutionary pattern similar to that reported from other samples of cosmic sources.

In the context of the VLA-COSMOS Deep project additional VLA A array observations at 1.4 GHz were obtained for the central degree of the COSMOS field and combined with the existing data from the VLA-COSMOS Large project. A newly constructed Deep mosaic with a resolution of 2.5" was used to search for sources down to 4 sigma with 1 sigma ~ 12 microJy/beam in the central 50'x50'. This new catalog is combined with the catalog from the Large project (obtained at 1.5"x1.4" resolution) to construct a new Joint catalog. All sources listed in the new Joint catalog have peak flux densities of >5 sigma at 1.5" and/or 2.5" resolution to account for the fact that a significant fraction of sources at these low flux levels are expected to be slighty resolved at 1.5" resolution. All properties listed in the Joint catalog such as peak flux density, integrated flux density and source size are determined in the 2.5" resolution Deep image. In addition, the Joint catalog contains 43 newly identified multi-component sources.

We present first pc-scale radio imaging of the radio-quiet candidate binary black hole system SDSS J1536+0441. The observations were carried out by the European VLBI Network at the frequency of 5 GHz and allowed to image SDSS J1536+0441 with a resolution of about 10 mas (50 pc). Two compact radio cores are detected at the position of the kpc-scale components VLA-A and VLA-B, proving the presence of two compact active nuclei with radio luminosity about 10^{40} erg/s, thus ruling out the possibility that the two radio sources are both powered by one 0.1 pc binary black hole. From a comparison with published 8.5 GHz flux densities we derived an estimate of the radio spectral index of the two pc-scale cores. Both cores have flat or inverted spectral index and, at least for the case of VLA-A, we can rule out the possibility that synchrotron self-absorption is responsible for the inverted radio spectrum. We suggest that thermal free-free emission from an X-ray heated disk wind may be powering the radio emission in VLA-A.

We model the evolution of the mean galaxy occupation of dark-matter halos over the range $0.1<z<1.3$, using the data from the VIMOS-VLT Deep Survey (VVDS). The galaxy projected correlation function $w_p(r_p)$ was computed for a set of luminosity-limited subsamples and fits to its shape were obtained using two variants of Halo Occupation Distribution models. These provide us with a set of best-fitting parameters, from which we obtain the average mass of a halo and average number of galaxies per halo. We find that after accounting for the evolution in luminosity and assuming that we are largely following the same population, the underlying dark matter halo shows a growth in mass with decreasing redshift as expected in a hierarchical structure formation scenario. Using two different HOD models, we see that the halo mass grows by 90% over the redshift interval z=[0.5,1.0]. This is the first time the evolution in halo mass at high redshifts has been obtained from a single data survey and it follows the simple form seen in N-body simulations with $M(z) = M_0 e^{-\beta z}$, and $\beta = 1.3 \pm 0.30$. This provides evidence for a rapid accretion phase of massive halos having a present-day mass $M_0 \sim 10^{13.5} h^{-1} M_\odot$, with a $m > 0.1 M_0$ merger event occuring between redshifts of 0.5 and 1.0. Futhermore, we find that more luminous galaxies are found to occupy more massive halos irrespectively of the redshift. Finally, the average number of galaxies per halo shows little increase from redshift z$\sim$ 1.0 to z$\sim$ 0.5, with a sharp increase by a factor $\sim$3 from z$\sim$ 0.5 to z$\sim$ 0.1, likely due to the dynamical friction of subhalos within their host halos.

The identities of the main processes triggering and quenching star-formation in galaxies remain unclear. A key stage in evolution, however, appears to be represented by post-starburst galaxies. To investigate their impact on galaxy evolution, we initiated a multiwavelength study of galaxies with k+a spectral features in the COSMOS field. We examine a mass-selected sample of k+a galaxies at z=0.48-1.2 using the spectroscopic zCOSMOS sample. K+a galaxies occupy the brightest tail of the luminosity distribution. They are as massive as quiescent galaxies and populate the green valley in the colour versus luminosity (or stellar mass) distribution. A small percentage (<8%) of these galaxies have radio and/or X-ray counterparts (implying an upper limit to the SFR of ~8Msun/yr). Over the entire redshift range explored, the class of k+a galaxies is morphologically a heterogeneous population with a similar incidence of bulge-dominated and disky galaxies. This distribution does not vary with the strength of the Hdelta absorption line but instead with stellar mass in a way reminiscent of the well-known mass-morphology relation. Although k+a galaxies are also found in underdense regions, they appear to reside typically in a similarly rich environment as quiescent galaxies on a physical scale of ~2-8Mpc, and in groups they show a morphological early-to-late type ratio similar to the quiescent galaxy class. With the current data set, we do not find evidence of statistical significant evolution in either the number/mass density of k+a galaxies at intermediate redshift with respect to the local values, or the spectral properties. Those galaxies, which are affected by a sudden quenching of their star-formation activity, may increase the stellar mass of the red-sequence by up to a non-negligible level of ~10%.

We investigate the properties and the environment of radio sources with optical counterpart from the combined VLA-COSMOS and zCOSMOS samples. The advantage of this sample is the availability of optical spectroscopic information, high quality redshifts, and accurate density determination. By comparing the star formation rates estimated from the optical spectral energy distribution with those based on the radio luminosity, we divide the radio sources in three families, passive AGN, non-passive AGN and star forming galaxies. These families occupy specific regions of the 8.0-4.5 $\mu$m infrared color--specific star formation plane, from which we extract the corresponding control samples. Only the passive AGN have a significantly different environment distribution from their control sample. The fraction of radio-loud passive AGN increases from ~2% in underdense regions to ~15% for overdensities (1+delta) greater than 10. This trend is also present as a function of richness of the groups hosting the radio sources. Passive AGN in overdensities tend to have higher radio luminosities than those in lower density environments. Since the black hole mass distribution is similar in both environments, we speculate that, for low radio luminosities, the radio emission is controlled (through fuel disponibility or confinement of radio jet by local gas pressure) by the interstellar medium of the host galaxy, while in other cases it is determined by the structure (group or cluster) in which the galaxy resides.

Hierarchical models of galaxy formation predict that the properties of a dark matter halo depend on the large-scale environment surrounding the halo. As a result of this correlation, we expect massive haloes to be present in larger number in overdense regions than in underdense ones. Given that a correlation exists between a galaxy stellar mass and the hosting dark matter halo mass, the segregation in dark matter halo mass should then result in a segregation in the distribution of stellar mass in the galaxy population. In this work we study the distribution of galaxy stellar mass and rest-frame optical color as a function of the large-scale galaxy distribution using the VLT VIMOS Deep Survey sample, in order to verify the presence of segregation in the properties of the galaxy population. We use the VVDS redshift measurements and multi-band photometric data to derive estimates of the stellar mass, rest-frame optical color, and of the large-scale galaxy density, on a scale of approximately 8 Mpc, for a sample of 5619 galaxies in the redshift range 0.2<z<1.4. We observe a significant mass and optical color segregation over the whole redshift interval covered by our sample, such that the median value of the mass distribution is larger and the rest-frame optical color is redder in regions of high galaxy density. The amplitude of the mass segregation changes little with redshift, at least in the high stellar mass regime that we can uniformely sample over the 0.2<z<1.4 redshift interval. The color segregation, instead, decreases significantly for z>0.7. However, when we consider only galaxies in narrow bins of stellar mass, in order to exclude the effects of the stellar mass segregation on the galaxy properties, we do not observe any more any significant color segregation.

From the VIMOS VLT Deep Survey we use a sample of 6447 galaxies with I_{AB} < 24 to identify 251 pairs of galaxies, each member with a secure spectroscopic redshift, which are close in both projected separation and in velocity. We find that at z ~ 0.9, 10.9 +/- 3.2 % of galaxies with M_B(z) < -18-Qz are in pairs with separations dr < 20 kpc/h, dv < 500 km/s, and with dM_B < 1.5, significantly larger than 3.76 +/- 1.71 % at z ~ 0.5; we find that the pair fraction evolves as (1+z)^m with m = 2.49 +/- 0.56. For brighter galaxies with M_B(z=0) < -18.77, the pair fraction is higher and its evolution with redshift is somewhat flatter with m=1.88 \pm 0.40, a property also observed for galaxies with increasing stellar masses. Early type, dry mergers, pairs increase their relative fraction from 3 % at z ~ 0.9 to 12 % at z ~ 0.5. We find that the merger rate evolves as N_{mg}=(9.05 +/- 3.76) * 10^{-4}) * (1+z)^{2.43 +/- 0.76}. We find that the merger rate of galaxies with M_B(z) < -18-Qz has significantly evolved since z ~ 1. The merger rate is increasing more rapidly with redshift for galaxies with decreasing luminosities, indicating that the flat evolution found for bright samples is not universal. The merger rate is also strongly dependent on the spectral type of galaxies involved, late type mergers being more frequent in the past, while early type mergers are more frequent today, contributing to the rise in the local density of early type galaxies. About 20 % of the stellar mass in present day galaxies with log(M/M_{sun}) > 9.5 has been accreted through major merging events since z ~ 1, indicating that major mergers have contributed significantly to the growth in stellar mass density of bright galaxies over the last half of the life of the Universe.

We explore the cosmic evolution of radio AGN with low radio powers (L_1.4GHz < 5\times10^25 W/Hz) out to z=1.3 using to-date the largest sample of ~600 low luminosity radio AGN at intermediate redshift drawn from the VLA-COSMOS survey. We derive the radio luminosity function for these AGN, and its evolution with cosmic time assuming two extreme cases: i) pure luminosity and ii) pure density evolution. The former and latter yield L_*\propto(1+z)^(0.8+/-0.1), and Phi*\propto (1+z)^(1.1+/-0.1), respectively, both implying a fairly modest change in properties of low radio-power AGN since z=1.3. We show that this is in stark contrast with the evolution of powerful (L_1.4GHz > 5\times10^25 W/Hz) radio AGN over the same cosmic time interval, constrained using the 3CRR, 6CE, and 7CRS radio surveys by Willott et al. (2001). We demonstrate that this can be explained through differences in black hole fueling and triggering mechanisms, and a dichotomy in host galaxy properties of weak and powerful AGN. Our findings suggest that high and low radio-power AGN activity is triggered in different stages during the formation of massive red galaxies. We show that weak radio AGN occur in the most massive galaxies already at z~1, and they may significantly contribute to the heating of their surrounding medium and thus inhibit gas accretion onto their host galaxies, as recently suggested for the `radio mode' in cosmological models.

We show preliminary results from a sample of Luminous and Ultra-Luminous Infrared Galaxies (LIRGs and ULIRGs, respectively) in the local universe, obtained from observations using the Very Large Array (VLA), the Multi-Element Radio Link Interferometer Network (MERLIN), and the European VLBI Network (EVN). The main goal of our high-resolution, high-sensitivity radio observations is to unveil the dominant gas heating mechanism in the central regions of local (U)LIRGs. The main tracer of recent star-formation in (U)LIRGs is the explosion of core-collapse supernovae (CCSNe), which are the endproducts of the explosion of massive stars and yield bright radio events. Therefore, our observations will not only allow us to answer the question of the dominant heating mechanism in (U)LIRGs, but will yield also the CCSN rate and the star-formation rate (SFR) for the galaxies of the sample.

We study the evolution of the low-order moments of the galaxy overdensity distribution over the redshift interval 0.7<z<1.5. We find that the variance and the normalized skewness evolve over this redshift interval in a way that is remarkably consistent with predictions of first- and second-order perturbation theory. This finding confirms the standard gravitational instability paradigm over nearly 9 Gyrs of cosmic time and demonstrates the importance of accounting for the non-linear component of galaxy biasing to avoid disagreement between theory and observations.

(abridged) We use the 1.4 GHz VIMOS-VLA Deep Survey and the optical VVDS and the CFHT-LS to compare the properties of radio loud galaxies with respect to the whole population of optical galaxies. The availability of multiband photometry and high quality photometric redshifts allows to derive rest frame colors and radio luminosity functions down to a limit of a B rest-frame magnitude of M=-20. Galaxy properties and luminosity functions (LFs) are estimated up to z~1 for radio loud and radio quiet early and late type galaxies. Radio loud late type galaxies are redder than radio quiet objects of the same class and this is an effect related to the presence of more dust in stronger star forming galaxies. Moreover, we estimate optical LFs, stellar masses and star formation rate distributions for radio sources and compare them with those derived for a well defined control sample, finding that the probability for a galaxy to be a radio emitter significantly increases at high values of these parameters. Radio loud early type galaxies show luminosity evolution of their bivariate radio-optical LF, due to an evolution in the radio-optical ratio. The lack of evolution of the mass function of radio loud early type galaxies means that no new AGN are formed at z<1. On the contrary, radio loud late type objects show a strong evolution, both in luminosity and in density, of the radio LF for z>0.7. This evolution is the direct effect of the strong optical evolution of this class and no significant change with redshift of the radio-optical ratio is required. With the knowledge of the radio-optical ratio and the optical and radio LFs for late type galaxies, we estimated the star formation history of the Universe up to z~1.5, using optical galaxies as tracers of the global radio emission.

We use one of the deepest spectroscopic samples of broad line active galactic nuclei (AGN) currently available, extracted from the VIMOS VLT Deep Survey (VVDS), to compute MgII and CIV virial masses estimate of 120 super-massive black holes in the redshift range 1.0<z<1.9 and 2.6<z<4.3. We find that the mass-luminosity relation shows considerably enhanced dispersion towards low AGN luminosities (log L_bol ~ 45). At these luminosities, there is a substantial fraction of black holes accreting far below their Eddington limit (L_bol/L_Edd < 0.1), in marked contrast to what is generally found for AGN of higher luminosities. We speculate that these may be AGN on the decaying branch of their light-curves, well past their peak activity. This would agree with recent theoretical predictions of AGN evolution. In the electronic Appendix of this paper we publish an update of the VVDS type-1 AGN sample, including the first and most of the second epoch observations. This sample contains 298 objects of which 168 are new.

Images made with the VLBA have resolved the region in a nearby radio galaxy, Pictor A, where the relativistic jet that originates at the nucleus terminates in an interaction with the intergalactic medium, a so-called radio galaxy hot spot. This image provides the highest spatial resolution view of such an object to date (16 pc), more than three times better than previous VLBI observations of similar objects. The north-west Pictor A hot spot is resolved into a complex set of compact components, seen to coincide with the bright part of the hot spot imaged at arcsecond-scale resolution with the VLA. In addition to a comparison with VLA data, we compare our VLBA results with data from the HST and Chandra telescopes, as well as new Spitzer data. The presence of pc-scale components in the hot spot, identifying regions containing strong shocks in the fluid flow, leads us to explore the suggestion that they represent sites of synchrotron X-ray production, contributing to the integrated X-ray flux of the hot spot, along with X-rays from synchrotron self-Compton scattering. This scenario provides a natural explanation for the radio morphology of the hot spot and its integrated X-ray emission, leading to very different predictions for the higher energy X-ray spectrum compared to previous studies. From the sizes of the individual pc-scale components and their angular spread, we estimate that the jet width at the hot spot is in the range 70 - 700 pc, which is comparable to similar estimates in PKS 2153-69, 3C 205, and 4C 41.17. The lower limit in this range arises from the suggestion that the jet may dither in its direction as it passes through hot spot backflow material close to the jet termination point, creating a "dentist drill" effect on the inside of a cavity 700 pc in diameter.

We derive the cosmic star formation history (CSFH) out to z=1.3 using a sample of ~350 radio-selected star-forming galaxies, a far larger sample than in previous, similar studies. We attempt to differentiate between radio emission from AGN and star-forming galaxies, and determine an evolving 1.4 GHz luminosity function based on these VLA-COSMOS star forming galaxies. We precisely measure the high-luminosity end of the star forming galaxy luminosity function (SFR>100 M_Sol/yr; equivalent to ULIRGs) out to z=1.3, finding a somewhat slower evolution than previously derived from mid-infrared data. We find that more stars are forming in luminous starbursts at high redshift. We use extrapolations based on the local radio galaxy luminosity function; assuming pure luminosity evolution, we derive $L_* \propto (1+z)^{2.1 \pm 0.2}$ or $L_* \propto (1+z)^{2.5 \pm 0.1}$, depending on the choice of the local radio galaxy luminosity function. Thus, our radio-derived results independently confirm the ~1 order of magnitude decline in the CSFH since z~1.

We selected a mass-limited sample of 4048 objects from the VIMOS VLT Deep Survey in the redshift interval 0.5<z<1.3. We used the amplitude of the 4000 Balmer break (Dn4000) to separate the galaxy population and the EW[OII]3727 line as proxy for the star formation activity. We discuss to what extent stellar mass drives galaxy evolution, showing for the first time the interplay between stellar ages and stellar masses over the past 8Gyr. Low-mass galaxies have small Dn4000 and at increasing stellar mass, the galaxy distribution moves to higher Dn4000 values as observed in the local Universe. As cosmic time goes by, we witness an increasing abundance of massive spectroscopically ET systems at the expense of the LT systems. This spectral transformation is a process started at early epochs and continuing efficiently down to the local Universe. This is confirmed by the evolution of our type-dependent stellar mass function. The underlying stellar ages of LT galaxies apparently do not show evolution, likely as a result of a continuous formation of new stars. All star formation activity indicators consistently point towards a star formation history peaked in the past for massive galaxies, with little or no residual star formation taking place in the most recent epochs. The activity and efficiency of forming stars are mechanisms that depend on stellar mass, and the mass assembly becomes progressively less efficient in massive systems as time elapses. The concepts of star formation downsizing and mass assembly downsizing describe a single scenario that has a top-down evolutionary pattern. The role of (dry) merging events seems to be only marginal at z<1.3, as our estimated efficiency in stellar mass assembly can possibly account for the progressive accumulation of passively evolving galaxies.

The VLA-COSMOS Large Project has imaged the 2 sq.deg. COSMOS field with a resolution of 1.5 arcsec and a sensitivity of about 11 microJy (1 sigma) yielding to a catalog of about 3600 radio sources. In this paper we present a further analysis of the VLA-COSMOS Large Project catalog of radio sources aimed to: 1) quantify and correct for the effect of bandwidth smearing in the catalog, 2) determine the incompleteness produced by the noise bias and the resolution bias in the new catalog and 3) derive the radio source counts at 1.4 GHz. The effect of bandwidth smearing on the radio sources in the catalog was quantified comparing the peak and total flux densities in the final mosaic and in each of the individual pointings where the source was closest to the center of the field. We find that the peak flux densities in the original VLA-COSMOS Large Project catalog have to be divided by a factor about 0.8 or 0.9, depending on the distance from the mosaic center. The completeness of the radio catalog has been tested using samples of simulated radio sources with different angular size distributions. These simulated sources have been added to the radio image and recovered using the same techniques used to produce the radio catalog. The fraction of missed sources as a function of the total flux density is a direct measure of the incompleteness. Finally, we derived the radio source counts down to 60 microJy with unprecedented good statistics. Comparison to the findings of other surveys shows good agreement in the flux density range 0.06-1 mJy confirming the upturn at about 0.5 mJy and a possible decline of the source counts below about 0.1 mJy.

We explore the properties of the submillijansky radio population at 20 cm by applying a newly developed optical color-based method to separate star forming (SF) from AGN galaxies at intermediate redshifts (z<1.3). Although optical rest-frame colors are used, our separation method is shown to be efficient, and not biased against dusty starburst galaxies. This classification method has been calibrated and tested on a local radio selected optical sample. Given accurate multi-band photometry and redshifts, it carries the potential to be generally applicable to any galaxy sample where SF and AGN galaxies are the two dominant populations. In order to quantify the properties of the submillijansky radio population, we have analyzed ~2,400 radio sources, detected at 20 cm in the VLA-COSMOS survey. 90% of these have submillijansky flux densities. We classify the objects into 1) star candidates, 2) quasi stellar objects, 3) AGN, 4) SF, and 5) high redshift (z>1.3) galaxies. We find, for the composition of the submillijansky radio population, that SF galaxies are not the dominant population at submillijansky flux levels, as previously often assumed, but that they make up an approximately constant fraction of 30-40% in the flux density range of ~50 microJy to 0.7 mJy. In summary, based on the entire VLA-COSMOS radio population at 20 cm, we find that the radio population at these flux densities is a mixture of roughly 30-40% of SF and 50-60% of AGN galaxies, with a minor contribution (~10%) of QSOs.

AIMS. We present a new Ks-band survey that represents a significant extension to the previous wide-field Ks-band imaging survey within the 0226-04 field of the VIMOS-VLT deep survey (VVDS). The new data add ~ 458 arcmin^2 to the previous imaging program, thus allowing us to cover a total contiguous area of ~ 600 arcmin^2 within this field. METHODS. Sources are identified both directly on the final K-band mosaic image and on the corresponding, deep chi^2-g'r'i' image from the CFHT Legacy Survey in order to reduce contamination while ensuring us the compilation of a truly K-selected catalogue down to the completeness limit of the Ks-band. The newly determined Ks-band magnitudes are used in combination with the ancillary multiwavelength data for the determination of accurate photometric redshifts. RESULTS. The final catalogue totals ~ 52000 sources, out of which ~ 4400 have a spectroscopic redshift from the VVDS first epoch survey. The catalogue is 90% complete down to K_Vega = 20.5 mag. We present K_s-band galaxy counts and angular correlation function measurements down to such magnitude limit. Our results are in good agreement with previously published work. We show that the use of K magnitudes in the determination of photometric redshifts significantly lowers the incidence of catastrophic errors. The data presented in this paper are publicly available through the CENCOS database.

We present a detailed analysis of the Galaxy Stellar Mass Function of galaxies up to z=2.5 as obtained from the VVDS. We estimate the stellar mass from broad-band photometry using 2 different assumptions on the galaxy star formation history and show that the addition of secondary bursts to a continuous star formation history produces systematically higher (up to 40%) stellar masses. At low redshift (z=0.2) we find a substantial population of low-mass galaxies (<10^9 Msun) composed by faint blue galaxies (M_I-M_K=0.3). In general the stellar mass function evolves slowly up to z=0.9 and more significantly above this redshift. Conversely, a massive tail is present up to z=2.5 and have extremely red colours (M_I-M_K=0.7-0.8). We find a decline with redshift of the overall number density of galaxies for all masses (59+-5% for M>10^8 Msun at z=1), and a mild mass-dependent average evolution (`mass-downsizing'). In particular our data are consistent with mild/negligible (<30%) evolution up to z=0.7 for massive galaxies (>6x10^10 Msun). For less massive systems the no-evolution scenario is excluded. A large fraction (>=50%) of massive galaxies have been already assembled and converted most of their gas into stars at z=1, ruling out the `dry mergers' as the major mechanism of their assembly history below z=1. This fraction decreases to 33% at z=2. Low-mass systems have decreased continuously in number and mass density (by a factor up to 4) from the present age to z=2, consistently with a prolonged mass assembly also at z<1.

By combining the VVDS with the SWIRE data, we have built the currently largest spectroscopic sample of galaxies selected in the rest-frame near-infrared. These allow us to investigate, for the first time using spectroscopic redshifts, the clustering evolution of galaxies selected from their rest-frame near-infrared luminosity in the redshift range 0.2<z<2.1. Therefore we use the projected two-point correlation function wp(rp) to study the three dimensional clustering properties of galaxies detected at 3.6 and 4.5 microns with IRAC and for which we have spectroscopic redshifts from first epoch VVDS. We find that in the flux limited samples at 3.6 and 4.5 microns, the clustering length does not change from redshift ~2 to the present. In addition, we find that the clustering lengths are systematically higher when galaxy samples are selected from increasingly redder wavelengths. These values are larger than the typical values found for I-band selected galaxies in the same redshift range, but we find that the difference in clustering length between I-band and 3.6/4.5 microns selected samples is decreasing with increasing redshift to become comparable at z~1.5. We interpret this as evidence that galaxies with older stellar populations and galaxies actively forming stars reside in comparably overdense environments at epochs earlier than z~1.5. The increasing difference in clustering length observed between rest-frame UV-optical and infrared selected samples could then be an indication that star formation is gradually shifting to lower density regions as cosmic time increases.

We present the evolution of the comoving SFR density in the redshift range 0 < z < 5 using the first epoch data release of the VVDS, that is 11564 spectra selected at I_AB=24 over 2200 arcmin^2 in two fields of view, the VVDS-0226-04 and the VVDS-CDFS-0332-27, and the cosmological parameters (Omega_M, Omega_L, h)=(0.3, 0.7, 0.7). We study the multi-wavelength non dust-corrected luminosity densities at 0 < z < 2 from the rest-frame FUV to the optical passbands, and the rest-frame 1500A luminosity functions and densities at 2.7 < z < 5. They evolve from z=1.2 to 0.05 according to (1+z)^{x} with x = 2.05, 1.94, 1.92, 1.14, 0.73, 0.42, 0.30 in the FUV-1500, NUV-2800, U-3600, B-4400, V-5500, R-6500, and I-7900 passbands, respectively. From z = 1.2 to 0.2 the B-band density for the irregular-like galaxies decreases markedly by a factor 3.5 while it increases by a factor 1.7 for the elliptical-like galaxies. We identify several SFR periods; from z = 5 to 3.4 the FUV-band density increases by at most 0.5dex, from z=3.4 to 1.2 it decreases by 0.08dex, from z=1.2 to 0.05 it declines steadily by 0.6dex. For the most luminous M_AB(1500) < -21 galaxies the FUV-band density drops by 2dex from z = 3.9 to 1.2, and for the intermediate -21 < M_AB(1500) < -20 galaxies it drops by 2dex from z = 0.2 to 0. Comparing with dust corrected surveys, at 0.4 < z < 2 the FUV seems obscured by a constant factor of ~1.8-2 mag, while at z < 0.5 it seems progressively less obscured by up to ~0.9-1 mag when the dust-deficient early-type population is increasingly dominating the B-band density. The VVDS results agree with a downsizing picture where the most luminous sources cease to efficiently produce new stars 12 Gyrs ago (at z~4), while intermediate luminosity sources keep producing stars until 2.5 Gyrs ago (at z~0.2).(abridged)

In a previous paper (Gavignaud et al. 2006), we presented the type-1 Active Galactic Nuclei (AGN) sample obtained from the first epoch data of the VIMOS-VLT Deep Survey (VVDS). The sample consists of 130 faint, broad-line AGN with redshift up to z=5 and 17.5< I <24.0, selected on the basis of their spectra. In this paper we present the measurement of the Optical Luminosity Function up to z=3.6 derived from this sample, we compare our results with previous results from brighter samples both at low and at high redshift. Our data, more than one magnitude fainter than previous optical surveys, allow us to constrain the faint part of the luminosity function up to high redshift. By combining our faint VVDS sample with the large sample of bright AGN extracted from the SDSS DR3 (Richards et al., 2006b) and testing a number of different evolutionary models, we find that the model which better represents the combined luminosity functions, over a wide range of redshift and luminosity, is a luminosity dependent density evolution (LDDE) model, similar to those derived from the major X-surveys. Such a parameterization allows the redshift of the AGN space density peak to change as a function of luminosity and explains the excess of faint AGN that we find at 1.0< z <1.5. On the basis of this model we find, for the first time from the analysis of optically selected samples, that the peak of the AGN space density shifts significantly towards lower redshift going to lower luminosity objects. This result, already found in a number of X-ray selected samples of AGN, is consistent with a scenario of "AGN cosmic downsizing", in which the density of more luminous AGN, possibly associated to more massive black holes, peaks earlier in the history of the Universe, than that of low luminosity ones.

We discuss the evolution of clustering of galaxies in the Universe from the present epoch back to z ~ 2, using the first-epoch data from the VIMOS-VLT Deep Survey (VVDS). We present the evolution of the projected two-point correlation function of galaxies for the global galaxy population, as well as its dependence on galaxy intrinsic luminosities and spectral types. While we do not find strong variations of the correlation function parameters with redshift for the global galaxy population, the clustering of objects with different intrinsic luminosities evolved significantly during last 8-10 billion years. Our findings indicate that bright galaxies in the past traced higher density peaks than they do now and that the shape of the correlation function of most luminous galaxies is different from observed for their local counterparts, which is a supporting evidence of a non-trivial evolution of the galaxy vs. dark matter bias.

We study the LF of the high-z galaxy population with 3<z<4 using a purely I-band magnitude-selected spectroscopic sample obtained in the framework of the VVDS. We determine the LF from the VVDS, taking care to add as few assumptions and as simple corrections as possible, and compare our results with those obtained from photometric studies, based on Lyman-break selections or photo-z measurements. We find that in the range 3<z<4, the VVDS LF is parameterized by phi*=1.24+-0.50 10-3 mag-1 Mpc-3 and M*=-21.49+-0.19, assuming a slope alpha=-1.4 consistent with most previous studies. While phi* is comparable to previously found values, M* is significantly brighter by about 0.5 mag at least. Using the conservative slope -1.4, we find a LD at 1700A rho(M<-18.5)=2.4 10 19 W Mpc-3 and rho Tot=3.1 10 19 W Mpc-3, comparable to that estimated in other studies. The unexpectedly large number of very bright galaxies found in the VVDS indicates that the color-selection and photo-z techniques that are generally used to build high-z galaxy samples may be affected by a significant fraction of color-measurement failures or by incomplete modelling of the mix of stellar emission, AGN contribution, dust absorption and intergalactic extinction assumed to identify high-z galaxies, making pure magnitude selection better able to trace the full population. Because of the difficulty to identify all low-luminosity galaxies in a spectroscopic survey, the LD could still be significantly underestimated. We also find that the relative contribution of the most luminous galaxies compared to the fainter ones is at least twice as large in the VVDS compared to former estimates. Therefore, the VVDS paints a quite different picture of the role of the most actively star-forming galaxies in the history of star formation.

The VLA-COSMOS large project is described and its scientific objective is discussed. We present a catalog of ~ 3,600 radio sources found in the 2deg^2 COSMOS field at 1.4 GHz. The observations in the VLA A and C configuration resulted in a resolution of 1.5''x1.4'' and a mean rms noise of ~ 10.5(15) uJy/beam in the central 1(2)deg^2. 80 radio sources are clearly extended consisting of multiple components, and most of them appear to be double-lobed radio galaxies. The astrometry of the catalog has been thoroughly tested and the uncertainty in the relative and absolute astrometry are 130mas and <55mas, respectively.

We investigate the redshift and luminosity evolution of the galaxy colour-density relation using the data from the First Epoch VIMOS-VLT Deep Survey (VVDS) on scales of R=5 h^(-1)Mpc up to redshift z ~ 1.5. While at lower redshift we confirm the existence of a steep colour-density relation, with the fraction of the reddest(/bluest) galaxies of the same luminosity increasing(/decreasing) as a function of density, this trend progressively disappears in the highest redshift bins investigated.Our results suggest the existence of an epoch (more remote for brighter galaxies) characterized by the absence of the colour-density relation on the R=5 h^(-1)Mpc scales investigated. The rest frame u*-g' colour-magnitude diagram shows a bimodal pattern in both low and high density environments up to redshift z ~ 1.5. We find that the bimodal distribution is not universal but strongly depends upon environment. Both the colour-density and colour-magnitude-density relations, on the R=5 h^(-1)Mpc scales, appear to be a transient, cumulative product of genetic and environmental factors that have been operating over at least a period of 9 Gyr. These findings support an evolutionary scenario in which star formation/gas depletion processes are accelerated in more luminous objects and in high density environments: star formation activity is progressively shifting with cosmic time towards lower luminosity galaxies (downsizing), and out of high density environments.

We present the low frequency (610 MHz) radio source counts of the VVDS-VLA field and investigate the radio spectral index properties of the sub-mJy population. We use new deep (r.m.s.$\simeq 50$ $\mu$Jy/beam) observations of the VVDS-VLA field obtained at 610 MHz with the GMRT and matched in resolution (6 arcsec) with already available VLA data at 1.4 GHz on the same field. We find evidence of a change of the dominant population of radio sources below 0.5 mJy (at 1.4 GHz): between 0.15 and 0.5 mJy the median spectral index is significantly flatter ($\alpha=-0.46\pm 0.03$) than that of brighter sources ($\alpha=-0.67\pm 0.05$). A relevant contribution below 0.5 mJy from a population of flat spectrum low luminosity compact AGNs and radio quiet QSOs could explain this effect. At even fainter flux density, between 0.10 and 0.15 mJy at 1.4 GHz, the median spectral index steepens again ($\alpha=-0.61\pm 0.04$) suggesting that the contribution of starburst galaxies becomes important below $\sim 0.2$ mJy. Finally we present a sample of 58 candidate ultra-steep sources with radio flux density from one to two orders of magnitude lower than any other sample of such objects.

We present a new K-band survey covering 623 arcmin$^2$ in the VVDS 0226-0430 deep field down to a limiting magnitude K$_{\rm{Vega}}$ $\leq$ 20.5. We use the spectroscopic sample extracted from this new K-band catalogue to assess the effectiveness of optical-near infrared color selections in identifying extreme classes of objects at high redshift.

We investigate the redshift and luminosity evolution of the galaxy colour-density relation using the data from the First Epoch VIMOS-VLT Deep Survey (VVDS). The size (6582 galaxies), depth (I_AB<=24) and redshift sampling rate of the survey enable us to reconstruct the 3D galaxy environment on relatively local scales (R=5 Mpc) up to z~1.5. Particular attention has been devoted to calibrate a density reconstruction scheme, which factors out survey selection effects and reproduces in an unbiased way the underlying `real' galaxy environment. While at lower redshift we confirm the existence of a steep colour-density relation, with the fraction of the reddest(/bluest) galaxies of the same luminosity increasing(/decreasing) as a function of density, this trend progressively disappears in the highest redshift bins investigated. The rest frame u*-g' colour-magnitude diagram shows a bimodal pattern in both low and high density environments up to z~1.5. We find that the bimodal distribution is not universal but strongly depends upon environment: at lower redshifts the colour-magnitude diagrams in low and high density regions are significantly different while the progressive weakening of the colour-density relation causes the two bimodal distributions to nearly mirror each other in the highest redshift bin investigated. Both the colour-density and the colour-magnitude-density relations appear to be a transient, cumulative product of genetic and environmental factors operating over at least a period of 9 Gyr. These findings support an evolutionary scenario in which star formation/gas depletion processes are accelerated in more luminous objects and in high density environments: star formation activity is shifting with cosmic time towards lower luminosity (downsizing), and out of high density environments.

We present the type-1 active galactic nuclei (AGN) sample extracted from the VIMOS VLT Deep Survey first observations of 21000 spectra in 1.75 square degree. This sample, which is purely magnitude limited, free of morphological or color selection biases, contains 130 broad line AGN (BLAGN) spectra with redshift up to 5. Our data are divided into a wide (Iab < 22.5) and a deep (Iab < 24) subsample containing 56 and 74 objects respectively. Because of its depth and selection criteria, this sample is uniquely suited to study the population of faint type-1 AGN. Our measured surface density (~ 472 +- 48 BLAGN per square degree with Iab < 24) is significantly higher than that of any other optically selected sample of BLAGN with spectroscopic confirmation. By applying a morphological and color analysis to our AGN sample we find that: (1)~23% of the AGN brighter than Iab=22.5 are classified as extended; this percentage increases to ~42% for those with z < 1.6; (2) a non-negligible fraction of our BLAGN are lying close to the color space area occupied by stars in u*-g' versus g'-r' color-color diagram. This leads us to the conclusion that classical optical ultraviolet preselection technique, if employed at such deep magnitudes (Iab=22.5) in conjuction with a preselection of point-like sources, can miss miss up to ~35% of the AGN population. Finally, we present a composite spectrum of our sample of objects. While the continuum shape is very similar to that of the SDSS composite at short wavelengths, it is much redder than it at lambda > 3000 A. We interpret this as due to significant contamination from emission of the host galaxies, as expected from the faint absolute magnitudes sampled by our survey.

From the first epoch observations of the VVDS up to z=1.5 we have derived luminosity functions (LF) of different spectral type galaxies. The VVDS data, covering ~70% of the life of the Universe, allow for the first time to study from the same sample and with good statistical accuracy the evolution of the LFs by galaxy type in several rest frame bands from a purely magnitude selected sample. The magnitude limit of the VVDS allows the determination of the faint end slope of the LF with unprecedented accuracy. Galaxies have been classified in four spectral classes, using their colours and redshift, and LFs have been derived in the U, B, V, R and I rest frame bands from z=0.05 to z=1.5. We find a significant steepening of the LF going from early to late types. The M* parameter is significantly fainter for late type galaxies and this difference increases in the redder bands. Within each of the galaxy spectral types we find a brightening of M* with increasing redshift, ranging from =< 0.5 mag for early type galaxies to ~1 mag for the latest type galaxies, while the slope of the LF of each spectral type is consistent with being constant with redshift. The LF of early type galaxies is consistent with passive evolution up to z~1.1, while the number of bright early type galaxies has decreased by ~40% from z~0.3 to z~1.1. We also find a strong evolution in the normalization of the LF of latest type galaxies, with an increase of more than a factor 2 from z~0.3 to z~1.3: the density of bright late type galaxies in the same redshift range increases of a factor ~6.6. These results indicate a strong type-dependent evolution and identifies the latest spectral types as responsible for most of the evolution of the UV-optical luminosity function out to z=1.5.

We have computed the evolution of the rest-frame B-band luminosity function (LF) for bulge and disk-dominated galaxies since z=1.2. We use a sample of 605 spectroscopic redshifts with I_{AB}<24 in the Chandra Deep Field South from the VIMOS-VLT Deep Survey, 3555 galaxies with photometric redshifts from the COMBO-17 multi-color data, coupled with multi-color HST/ACS images from the Great Observatories Origin Deep Survey. We split the sample in bulge- and disk-dominated populations on the basis of asymmetry and concentration parameters measured in the rest-frame B-band. We find that at z=0.4-0.8, the LF slope is significantly steeper for the disk-dominated population (\alpha=-1.19 \pm 0.07) compared to the bulge-dominated population (\alpha=-0.53 \pm 0.13). The LF of the bulge-dominated population is composed of two distinct populations separated in rest-frame color: 68% of red (B-I)_{AB}>0.9 and bright galaxies showing a strongly decreasing LF slope \alpha=+0.55 \pm 0.21, and 32% of blue (B-I)_{AB}<0.9 and more compact galaxies which populate the LF faint-end. We observe that red bulge-dominated galaxies are already well in place at z~1, but the volume density of this population is increasing by a factor 2.7 between z~1 and z~0.6. It may be related to the building-up of massive elliptical galaxies in the hierarchical scenario. In addition, we observe that the blue bulge-dominated population is dimming by 0.7 magnitude between z~1 and z~0.6. Galaxies in this faint and more compact population could possibly be the progenitors of the local dwarf spheroidal galaxies.

We present photometric redshifts for an uniquely large and deep sample of 522286 objects with i'_{AB}<25 in the Canada-France Legacy Survey ``Deep Survey'' fields, which cover a total effective area of 3.2 deg^2. We use 3241 spectroscopic redshifts with 0<z<5 from the VIMOS VLT Deep Survey as a calibration to derive these photometric redshifts. We devise a robust calibration method which removes systematic trends in the photometric redshifts and significantly reduces the fraction of catastrophic errors. We use our unique spectroscopic sample to present a detailed assessment of the robustness of the photometric redshift sample. For a sample selected at i'_{AB}<24, we reach a redshift accuracy of \sigma_{\Delta z/(1+z)}=0.037 with \eta=3.7% of catastrophic error. The reliability of our photometric redshifts is lower for fainter objects: we find \sigma_{\Delta z/(1+z)}=0.029, 0.043 and \eta=1.7%, 5.4% for samples selected at i'_{AB}=17.5-22.5 and 22.5-24 respectively. We find that the photometric redshifts of starburst galaxies in our sample are less reliable: although these galaxies represent only 18% of the spectroscopic sample they are responsible for 54% of the catastrophic errors. We find an excellent agreement between the photometric and the VVDS spectroscopic redshift distributions at i'_{AB}<24. Finally, we compare the redshift distributions of i' selected galaxies on the four CFHTLS deep fields, showing that cosmic variance is already present on fields of 0.8 deg^2.

We measure the evolution of clustering for galaxies with different spectral types from 6495 galaxies with 17.5<=I_AB<=24 and measured spectroscopic redshift in the first epoch VIMOS-VLT Deep Survey. We classify our sample into 4 classes, based on the fit of well-defined galaxy spectral energy distributions on observed multi-color data. We measure the projected function wp(rp) and estimate the best-fit parameters for a power-law real-space correlation function. We find the clustering of early-spectral-type galaxies to be markedly stronger than that of late-type galaxies at all redshifts up to z<=1.2. At z~0.8, early-type galaxies display a correlation length r_0=4.8+/-0.9h^{-1}Mpc, while late types have r_0=2.5+/-0.4h^{-1}Mpc. The clustering of these objects increases up to r_0=3.42+/-0.7h^{-1}Mpc for z~1.4. The relative bias between early- and late-type galaxies within our magnitude-limited survey remains approximately constant with b~1.7-1.8 from z~=0.2 up to z~=1, with indications for a decrease at z>1.2, due to the growth in clustering of the star-forming population. We find similar results when splitting the sample into `red' and `blue' galaxies using the observed color bi-modality. When compared to the expected linear growth of mass fluctuations, a natural interpretation of these observations is that: (a) the assembly of massive early type galaxies is already mostly complete in the densest dark matter halos at z~=1; (b) luminous late-type galaxies are located in higher-density, more clustered regions of the Universe at z~=1.5 than at present, indicating that star formation activity is progressively increasing, going back in time, in the higher-density peaks that today are mostly dominated by old galaxies.

We measure the evolution of the galaxy Luminosity Function as a function of large-scale environment up to z=1.5 from the VIMOS-VLT Deep Survey (VVDS) first epoch data. The 3D galaxy density field is reconstructed using a sample of 6582 galaxies with 17.5 < I_{AB} < 24 and measured spectroscopic redshifts. We split the sample in four redshift bins up to z=1.5 and in under-dense and over-dense environments according to the average density contrast \delta=0. There is a strong dependence of the Luminosity Function (LF) with large-scale environment up to z=1.2: the LF shape is observed to have a steeper slope in under-dense environments. We find a continuous brightening of \Delta M* ~0.6 mag from z=0.25 to z=1.5 both in under-dense and over-dense environments. The rest-frame B-band luminosity density continuously increases in under-dense environments from z=0.25 to z=1.5 whereas its evolution in over-dense environments presents a peak at z~0.9. We interpret the peak by a complex interplay between the decrease of the star formation rate and the increasing fraction of galaxies at \delta>0 due to hierarchical growth of structures. As the environmental dependency of the LF shape is already present at least up to z=1.2, we therefore conclude that either the shape of the LF is imprinted very early on in the life of the Universe, a `nature' process, or that `nurture' physical processes shaping up environment relation have already been efficient earlier than a look-back time corresponding to 30% of the current age of the Universe.

In this paper we present a new deep, wide-field near-infrared imaging survey. Our J- and K-band observations in four separate fields complement optical BVRI, ultraviolet and spectroscopic observations undertaken as part of the VIMOS-VLT deep survey (VVDS). In total, our survey spans ~400arcmis^2. Our catalogues are reliable in all fields to at least Kvega~20.75 and Jvega~21.50 (defined as the magnitude where object contamination is less than 10% and completeness greater than 90%). Taken together these four fields represents a unique combination of depth, wavelength coverage and area. We describe the complete data reduction process and outline a comprehensive series of tests carried out to characterise the reliability of the final catalogues. We compare the statistical properties of our catalogues with literature compilations. We find that our J- and K-selected galaxy counts are in good agreement with previously published works, as are our (J-K) versus K colour-magnitude diagrams. Stellar number counts extracted from our fields are consistent with a synthetic model of our galaxy. Using the location of the stellar locus in colour-magnitude space and the measured field-to-field variation in galaxy number counts we demonstrate that the absolute accuracy of our photometric calibration is at the 5% level or better. Finally, an investigation of the angular clustering of K- selected extended sources in our survey displays the expected scaling behaviour with limiting magnitude, with amplitudes in each magnitude bin in broad agreement with literature values.

This paper presents the ``First Epoch'' sample from the VIMOS VLT Deep Survey (VVDS). The VVDS goals, observations, data reduction with VIPGI, and redshift measurement with KBRED are discussed. Data have been obtained with the VIsible Multi Object Spectrograph (VIMOS) on the ESO-VLT UT3, allowing to observe ~600 slits simultaneously at R~230. A total of 11564 objects have been observed in the VVDS-02h and VVDS-CDFS Deep fields over a total area of 0.61deg^2, selected solely on the basis of apparent magnitude 17.5 <=I_{AB} <=24. The VVDS covers the redshift range 0 < z <= 5. It is successfully going through the ``redshift desert'' 1.5<z<2.2, while the range 2.2<z<2.7 remains of difficult access because of the VVDS wavelength coverage.A total of 9677 galaxies have a redshift measurement, 836 are stars, 90 are AGNs, and a redshift could not be measured for 961 objects. There are 1065 galaxies with a measured redshift z>1.4. The survey reaches a redshift measurement completeness of 78% overall (93% including less reliable objects), with a spatial sampling of the population of galaxies of 25% and ~30% in the VVDS-02h and VVDS-CDFS. The redshift accuracy measured from repeated observations with VIMOS and comparison to other surveys is ~276km/s. From this sample we present for the first time the redshift distribution of a magnitude limited spectroscopic sample down to IAB=24. The redshift distribution has a median of z=0.62, z=0.65, z=0.70, and z=0.76, for magnitude limited samples with IAB<=22.5, 23, 23.5, and 24. A high redshift tail above redshift 2 and up to redshift 5 becomes readily apparent for IAB>23.5, probing the bright star forming population of galaxies. This sample provides an unprecedented dataset to study galaxy evolution over 90% of the life of the universe

We are presenting in this paper a detailed account of the methods used to compute the three-dimensional two-point galaxy correlation function in the VIMOS-VLT deep survey (VVDS). We investigate how instrumental selection effects and observational biases affect the measurements and identify the methods to correct them. We quantify the accuracy of our correction method using an ensemble of fifty mock galaxy surveys generated with the GalICS semi-analytic model of galaxy formation which incorporate the same selection biases and tiling strategy as the real data does. We demonstrate that we are able to recover the real-space two-point correlation function xi(s) to an accuracy better than 10% on scales larger than 1 h^{-1} Mpc, and of about 30% on scales below 1 h^{-1} Mpc, with the sampling strategy used for the first epoch VVDS data. The projected correlation function w_p(r_p) is recovered with an accuracy better than 10% on all scales 0.1 <= r <= 10 h^{-1} Mpc. There is a tendency for a small but systematic under-estimate of the correlation length derived from w_p(r_p) of 6% on average, remaining after our correction process. The large number of simulated surveys allows us to provide a reliable estimate of the cosmic error on the measurements of the correlation length r_0, of about 15-20% for the first epoch VVDS observation (Le Fevre et al.2004, astro-ph/0409133). The error estimation and measurement techniques outlined in this paper are being used in several studies which investigate in detail the clustering properties of galaxies in the VVDS data.

We investigate the evolution of the galaxy luminosity function from the VIMOS-VLT Deep Survey (VVDS) from the present to z=2 in five (U, B, V, R and I) rest-frame band-passes. We use the first epoch VVDS deep sample of 11,034 spectra selected at 17.5 <= I_{AB} <= 24.0, on which we apply the Algorithm for Luminosity Function (ALF), described in this paper. We observe a substantial evolution with redshift of the global luminosity functions in all bands. From z=0.05 to z=2, we measure a brightening of the characteristic magnitude M* included in the magnitude range 1.8-2.5, 1.7-2.4, 1.2-1.9, 1.1-1.8 and 1.0-1.6 in the U, B, V, R and I rest-frame bands, respectively. We confirm this differential evolution of the luminosity function with rest-frame wavelength, from the measurement of the comoving density of bright galaxies (M < M*(z=0.1)). This density increases by a factor of around 2.6, 2.2, 1.8, 1.5, 1.5 between z=0.05 and z=1 in the U, B, V, R, I bands, respectively. We also measure a possible steepening of the faint-end slope of the luminosity functions, with \Delta\alpha ~ -0.3 between z=0.05 and z=1, similar in all bands.

This paper presents the evolution of the clustering of the main population of galaxies from z=2.1 to z=0.2, from the first epoch VIMOS VLT Deep Survey (VVDS), a magnitude limited sample with 17.5<=I_{AB}<=24. We have computed the correlation functions \xi(r_p,\pi) and w_p(r_p), and the correlation length r_0(z), for the VVDS-02h and VVDS-CDFS fields, for a total of 7155 galaxies in a 0.61deg^2 area. We find that the correlation length in this sample stays roughly constant from z=0.5 to z=1.1, with r_0(z)=2.5-2.8 h^{-1} Mpc (comoving), for galaxies comparable in luminosity to the local 2dFGRS and SDSS samples, indicating that the amplitude of the correlation function was ~2.5x lower at z~1 than observed locally. The correlation length in our lowest redshift bin z=[0.2,0.5] is r_0=2.4 h^{-1} Mpc, lower than for any other population at the same redshift, indicating the low clustering of very low luminosity galaxies, 1.5 magnitudes fainter than in the 2dFGRS or SDSS. The correlation length is increasing to r_0~3.0 h^{-1} Mpc at higher redshifts z=[1.3,2.1], as we are observing increasingly brighter galaxies, comparable to galaxies with MB_AB=-20.5 locally. We compare our measurement to the DEEP2 measurements in the range z=[0.7,1.35] \citep{coil} on the population selected applying the same magnitude and color selection criteria as in their survey, and find comparable results. The slowly varying clustering of VVDS galaxies as redshift increases is markedly different from the predicted evolution of the clustering of dark matter, indicating that bright galaxies are already tracing the large scale structures emerging from the dark matter distribution 9-10 billion years ago, a supporting evidence for a strong evolution of the galaxy vs. dark matter bias.

In this paper we present the optical and near-infrared identifications of the 1054 radio sources detected in the 20cm deep radio survey down to a 5sigma flux limit of about 80 microJy obtained with the VLA in the VIMOS VLT Deep Survey VVDS-02h deep field. Using U,B,V,R,I and K data, we identified 718 radio sources (~74% of the whole sample). The photometric redshift analysis shows that, in each magnitude bin, the radio sample has a higher median photometric redshift than the whole optical sample, while the median (V-I) color of the radio sources is redder than the median color of the whole optical sample. These results suggest that radio detection is preferentially selecting galaxies with higher intrinsic optical luminosity. From the analysis of the optical properties of the radio sources as function of the radio flux, we found that while about 35% of the radio sources are optically unidentified in the higher radio flux bin (S> 1.0 mJy), the percentage of unidentified sources decreases to about 25% in the faintest bins (S< 0.5 mJy). The median I magnitude for the total sample of radio sources,i.e. including also the unidentified ones, is brighter in the faintest radio bins than in the bin with higher radio flux. This suggests that most of the faintest radio sources are likely to be associated to relatively lower radio luminosity objects at relatively modest redshift, rather than radio-powerful, AGN type objects at high redshift.

We present the first measurements of the Probability Distribution Function (PDF) of galaxy fluctuations in the VIMOS-VLT Deep Survey (VVDS) cone, covering 0.4x0.4 deg between 0.4<z<1.5. The second moment of the PDF, i.e. the rms fluctuations of the galaxy density field, is with good approximation constant over the full redshift baseline investigated: we find that, in redshift space, sigma_8 for galaxies brighter than M=-20+5log h has a mean value of 0.94\pm0.07 in the redshift interval 0.7<z<1.5. The third moment, i.e. the skewness, increases with cosmic time: we find that the probability of having underdense regions is greater at z~0.7 than it was at z~1.5. By comparing the PDF of galaxy density contrasts with the theoretically predicted PDF of mass fluctuations we infer the redshift-, density-, and scale-dependence of the biasing function b(z, \delta, R) between galaxy and matter overdensities up to redshift z=1.5. Our results can be summarized as follows: i) the galaxy bias is an increasing function of redshift: evolution is marginal up to z~0.8 and more pronounced for z>0.8; ii) the formation of bright galaxies is inhibited below a characteristic mass-overdensity threshold whose amplitude increases with redshift and luminosity; iii) the biasing function is non linear in all the redshift bins investigated with non-linear effects of the order of a few to 10% on scales >5Mpc.

We report on 1.6 and 5.0 GHz observations of the ultraluminous infrared galaxy (ULIRG) Mrk 273, using the European VLBI Network (EVN) and the Multi-Element Radio-Linked Interferometer Network (MERLIN). We also make use of published 1.4 GHz VLBA observations of Mrk 273 by Carilli & Taylor (2000). Our 5 GHz images have a maximum resolution of 5-10 mas, which corresponds to linear resolutions of 3.5-7 pc at the distance of Mrk 273, and are the most sensitive high-resolution radio observations yet made of this ULIRG. Component N1, often pinpointed as a possible AGN, displays a steep spectral index ($\alpha = 1.2 \pm 0.1; S_\nu \propto \nu^{-\alpha}$); hence it is very difficult to reconcile with N1 being an AGN, and rather suggests that the compact nonthermal radio emission is produced by an extremely high luminous individual radio supernova (RSN), or a combination of unresolved emission from nested supernova remnants (SNR), luminous RSNe, or both. Component N2 is partly resolved out into several compact radio sources --none of which clearly dominates-- and a region of extended emission about 30 pc in size. The integrated spectral index of this region is flat ($\alpha = 0.15 \pm 0.1$), which can be interpreted as due to a superposition of several unresolved components, e.g., RSNe or SNRs, whose radio emission peaks at different frequencies and is partially free-free absorbed. The overall extended radio emission from component N is typical of nonthermal, optically thin radio emission ($\alpha = 0.8 \pm 0.1$), and its 1.4 GHz luminosity ($L_{1.4 \rm GHz} = (2.2 \pm 0.1)\times 10^{23} $ WHz$^{-1}$) is consistent with being produced by relativistic electrons diffused away from supernova remnants in an outburst.

We present 5 GHz EVN+MERLIN observations of the nuclear region of the ultra luminous infrared galaxy Mrk 273. These observations confirms the detection, on the mas scale, of two resolved component labelled as N and SE in the literature. We use published VLBA observations at 1.4 GHz, at the same resolution, to derive spectral index information of component N and SE and discuss these findings in relation with different hypothesis (compact starburst or AGN) for the origin of the radio emission.

We present the results from a Chandra observation of the powerful radio galaxy 3C 265. We detect X-ray emission from the nucleus, the radio hotspots and lobes. In particular, the lobe X-ray emission is well explained as anisotropic inverse Compton scattering of the nuclear photons by the relativistic electrons in the radio lobes; the comparison between radio synchrotron and IC emission yields a magnetic field strength a factor about 2 lower than that calculated under minimum energy conditions. The X-ray spectrum of the nucleus is consistent with that of a powerful, strongly absorbed quasar and the X-ray emission of the south-eastern hotspot can be successfully reproduced by a combination of synchro-self Compton and inverse Compton emission assuming a magnetic field slightly lower than equipartition.

We present a catalog of 246 radio sources found in the central 1 degree of the COSMOS field at 1.4GHz. The VLA pilot project data have a resolution of 1.9"x1.6" and an rms noise limit of ~25-100uJy/beam covering 0.837 sqrdeg. About 20 radio sources are clearly extended and most of them appear to be double-lobed radio galaxies. We find evidence for a cluster of 7 radio sources with an extent of ~10' southeast of the COSMOS field center. This VLA pilot project was undertaken to demonstrate the feasibility of wide-field mosaicking at 2'' resolution at 1.4GHz using the VLA in its A array configuration. The 7-point mosaic data was used to develop the techniques necessary for reduction and analysis. These data will provide the initial astrometric frame for the optical (ground- and space-based) data of the COSMOS 2 sqrdeg survey. In addition, it demonstrates the feasibility of obtaining deep (rms ~ few uJy) radio imaging of this field at 1.4GHz, since the brightest radio emission peak detected within the area covered has a flux density of 13mJy/beam and no strong side-lobes from sources surrounding the COSMOS field were detected. Comparison of the number counts to other deep radio surveys shows that the COSMOS field is a representative deep field in the radio domain.

We address the question of why low-luminosity radio sources with similar flat radio spectra show a range of optical activity. The investigation is based on the spectral energy distributions (SEDs) of objects from the 200 mJy sample. We gathered new data from the VLA at 43 GHz, from SCUBA in the JCMT at 2000, 1350 and 850 micron and from the ISOPHOT instrument on ISO at 170, 90, 60 and 25 micron. There is considerable diversity amongst the SEDs of the objects: there are objects with steep broadband spectra between centimetre and millimetre bands (14% of the sample), there are those with flat broadband spectra over most of the spectral range (48% of the sample) and there are those which show pronounced sub-mm/infrared excesses (27% of the sample). Some objects of the first group have two sided radio morphology indicating that their pc-scale emission is not dominated by beamed jet emission. Amongst the objects that have smooth broadband spectra from the radio to the infrared there are passive elliptical galaxies as well as the expected BL Lacs objects. The most pronounced sub-mm/infrared excesses are shown by the broad emission line objects.

The VIMOS VLT Deep Survey (VVDS) is underway to study the evolution of galaxies, large scale structures and AGNs, from the measurement of more than 100000 spectra of faint objects. We present here the results from the first epoch observations of more than 20000 spectra. The main challenge of the program, the redshift measurements, is described, in particular entering the ``redshift desert'' in the range $1.5<z<3$ for which only very weak features are observed in the observed wavelength range. The redshift distribution of a magnitude limited sample brighter than $I_{AB}=24$ is presented for the first time, showing a peak at a low redshift $z\sim0.7$, and a tail extending all the way above $z=4$. The evolution of the luminosity function out to $z=1.5$ is presented, with the LF of blue star forming galaxies carrying most of the evolution, with L* changing by more than two magnitudes for this sub-sample.

We report the results from a Chandra observation of the powerful nearby (z=0.1744) radio galaxy 3C 219. We find evidence for non-thermal X-ray emission from the radio lobes which fits fairly well with a combination of inverse Compton scattering of Cosmic Microwave Background radiation and of nuclear photons with the relativistic electrons in the lobes. The comparison between radio synchrotron and IC emission yields a magnetic field strength significantly lower (about a factor 3) than that calculated under minimum energy conditions; the source energetics is then dominated by the relativistic particles.

We present HST observations of previously undetected optical jets in the low-luminosity radio galaxies B2 0755+37 and B2 1553+24. We show that there is accurate spatial coincidence between optical and radio emission, implying that the former is likely to be synchrotron radiation. The physical properties of the jets are similar to those known previously: their radio-optical spectral indices are ~0.7 and in B2 0755+37, the spectrum steepens between optical and X-ray wavelengths. Our results support the hypothesis that optical emission is detectable from jets orientated within ~20 degrees of the line of sight for the B2 sample.

We present MERLIN, global VLBI and VLBA observations of the high-luminosity, compact steep-spectrum quasar B1524$-$136 at cm wavelengths. These observations reveal well-defined radio jets on both sides of the active nucleus, a situation which is almost unique amongst high-luminosity radio quasars. However, the radio jets on opposite sides are very dissimilar, and the overall radio structure appears highly distorted. We discuss possible implications of these observations.

Two Compact Steep-spectrum Sources (CSSs), 0548+165 and 1524-136, chosen from a list of CSSs with polarization percentages that decrease with decreasing frequency and high rotation measure values (RM > 450 rad/m**2) on arcsecond scales, were observed with the VLBA at 4.9 and 8.4 GHz. RM values up to ~10**4 rad/m**2 were found in several regions along the jets in both sources. We suggest that a thin screen of magneto-ionic material with about 1 kpc thickness is responsible for these high RMs. The observed depolarization may be due to beam depolarization and/or inhomogeneities in the magnetic field.

We report the results from a {\it Chandra} observation of the powerful nearby (z=0.1744) radio galaxy 3C 219. We find evidence for non--thermal X-ray emission from the radio lobes which fits fairly well with a combination of IC scattering of CMB and of nuclear photons. The comparison between radio synchrotron and IC fluxes yields a magnetic field strength significantly lower ($\sim 3$) than that calculated under minimum energy conditions; the majority of the energetics is then associated to the relativistic particles.

Several blazars, and BL Lac objects in particular, show a misalignment between the jet orientation on parsec and kiloparsec scales. Some authors (i.e. Conway & Murphy, 1993) have attempted to explain this behaviour invoking helical jets for misalignment angles around 90$\degr$, showing how in this case there are interesting implications for the understanding of the medium into which the jet is expanding. By comparing sensitive VLA observations (Cassaro et al., 1999) with images available in the literature for the BL Lac objects from the 1-Jy Sample (Stickel et al., 1991), it is clear that there is a wide range of misalignments between the initial jet direction and the kpc-scale jet, when detected. We have carried out VLBI observations of these BL Lac objects, in order to investigate the spatial evolution of the radio jets from few tens to hundreds of mas, and to search for helical jets in this class of sources. We present here the first dataset obtained from EVN+MERLIN observations at 5 GHz for seven objects. From these observations we never have a clear detection of helical jets, we only have a possible signature of their presence in 2 objects. In only one of the sources with a misalignment angle around 90$\degr$ the presence of helical jets can be ruled out. This implies that it is not possible to invoke helical jets to explain the morphology of all the sources showing a misalignment of about 90$\degr$ between the parsec and the kiloparsec scale jets.

We report on the {\it Chandra} discovery of large scale non--thermal emission features in the double lobed SSRL quasar 3C 207 (z=0.684). These are: a diffuse emission well correlated with the western radio lobe, a bright one sided jet whose structure coincides with that of the eastern radio jet and an X-ray source at the tip of the jet coincident with the hot spot of the eastern lobe. The diffuse X-ray structure is best interpreted as inverse Compton (IC) scattering of the IR photons from the nuclear source and provides direct observational support to an earlier conjecture (Brunetti et al., 1997) that the spectrum of the relativistic electrons in the lobes of radio galaxies extends to much lower energies than those involved in the synchrotron radio emission. The X-ray luminous and spatially resolved knot along the jet is of particular interest: by combining VLA and {\it Chandra} data we show that a SSC model is ruled out, while the X-ray spectrum and flux can be accounted for by the IC scattering of the CMB photons (EIC) under the assumptions of a relatively strong boosting and of an energy distribution of the relativistic electrons as that expected from shock acceleration mechanisms. The X-ray properties of the hot spot are consistent with a SSC model. In all cases we find that the inferred magnetic field strength are lower, but close to the equipartition values. The constraints on the energy distribution of the relativistic electrons, imposed by the X-ray spectra of the observed features, are discussed. To this aim we derive in the Appendices precise semi--analytic formulae for the emissivities due to the SSC and EIC processes.

We present the results of a multi-frequency study of the inversion symmetric radio galaxy NGC326 based on Very Large Array observations at 1.4, 1.6, 4.8, 8.5 and 14.9 GHz. The morphological, spectral and polarization properties of this peculiar object are studied at different levels of spatial resolutions. The interpretation of the data will be discussed in forthcoming papers.

In this letter we report a Chandra X-ray detection of the double northern hot spot of the radio quasar 3C 351. The hot spot has also been observed in the optical with the Hubble Space Telescope (R-band) and with the 3.5m. Telescopio Nazionale Galileo (B-band). The radio-to-optical and X-ray spectra are interpreted as the results of the synchrotron and synchrotron-self-Compton (SSC) mechanisms, respectively, with hot-spot magnetic field strengths ~3 times smaller than the equipartition values. In the framework of shock acceleration theory, we show that the requirement for such a relatively small field strength is in agreement with the fitted synchrotron spectral models and with the sizes of the hot spots. Finally, we show that the combination of a lower magnetic field strength with the high frequencies of the synchrotron cut-off in the fitted synchrotron spectra provides strong evidence for electron acceleration in the hot spots.

In this paper we present the results of a 4-year (1996 - 1999) radio flux density monitoring program for a sample of X- and $\gamma$-ray loud blazars. Our program started in January 1996 and was carried out on monthly basis at the frequencies of 5 GHz and 8.4 GHz with the 32-m antennas located in Medicina (Bologna, Italy) and Noto (Siracusa, Italy). 22 GHz data collected in Medicina from January 1996 to June 1997 will also be presented. The sample of selected sources comprises most radio loud blazars with $\delta \ge -10^{\circ}$ characterised by emission in the X- and $\gamma$-ray regimes, and target sources for the BeppoSAX X-ray mission. All sources in the sample, except J1653+397 (MKN 501), are variable during the four years of our monitoring program. We classified the type of variability in each source by means of a structure function analysis. We also computed th$\alpha_{5}^{8.4}$ for all epochs with nearly simultaneous observations ate spectral index and found that $\alpha_{5}^{8.4}$ starts flattening at the very beginning o f a radio flare, or flux density increase.

We have observed 28 sources selected from the 1Jy sample of BL Lac objects (Stickel et al. 1991) with the Very Large Array (VLA) in A, B and D configurations at 1.36, 1.66 and 4.85 GHz, and/or with the Westerbork Synthesis Radio Telescope (WSRT) at 1.40 GHz. In this paper we present high sensitivity images at arcsecond resolution of the 18 objects showing extended structure in our images, and of another source from the FIRST (Faint Images of the Radio Sky at Twenty-cm) survey (Becker et al. 1995). In general our high sensitivity images reveal an amount of extended emission larger than previously reported. In some objects the luminosity of the extended structure is comparable with that of FR~II radio sources. A future paper will be devoted to the interpretation of these results.

We report on centimeter VLA and VLBI observations of the giant, low power radio galaxy 1144+35. These observations are sensitive to structures on scales from less than 1 parsec to greater than 1 megaparsec. Diffuse steep spectrum lobes on the megaparsec scale are consistent with an age of $\sim$ 10$^8$ years. On the parsec scale, a complex jet component is seen to move away from the center of activity with an apparent velocity 2.7 h$^{-1}_{50}$ c. It shows a central spine -- shear layer morphology. A faint parsec scale counterjet is detected and an intrinsic jet velocity of 0.95 c and angle to the line of sight of 25$^\circ$ are derived, consistent with an intrinsically symmetric ejection. The central spine in the parsec scale jet is expected to move at a higher velocity and a Lorentz factor $\gamma$ $\sim$ 15 has been estimated near the core.The age of this inner VLBI structure is $\sim$ 300 years. Assuming a constant angle to the line-of-sight, the jet velocity is found to decrease from 0.95 c at 20 mas (32 pc on the plane of the sky) to 0.02 c at 15 arcsec (24 kpc on the plane of the sky). These findings lend credence to the claim that (1) even the jets of low power radio galaxies start out relativistic; and (2) these jets are decelerated to subrelativistic velocities by the time they reach kiloparsec scales.

The class of Compact Steep-spectrum Sources is dominated by double-lobed objects (70%). The remaining 30% are jet-dominated objects, with the jet brightened either by Doppler boosting or by interaction with the ambient media. We show that there is both observational and statistical evidence in favour of an interaction between jets and dense gas clouds. Such an interaction should happen in the Narrow Line Regions. The images of four CSSs observed by us with VLBI are also presented. These sources do show large bent jets in the first kpc from the nucleus.