The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

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 use the LOFAR Two-metre Sky Survey (LoTSS) Data Release I to identify the groups of galaxies (and individual galaxies) from the Hickson Compact Groups and Magnitude Limited Compact Groups samples that emit at the frequency of 150\,MHz, characterise their radio emission (extended or limited to the galaxies), and compare new results to earlier observations and theoretical predictions. The detection of 73 systems (and 7 more -- probably) out of 120, of which as many as 17 show the presence of extended radio structures, confirms the previous hypothesis of the common character of the magnetic field inside galaxy groups and its detectability. In order to investigate the future potential of low-frequency radio studies of galaxy groups, we also present a more detailed insight into four radio-emitting systems, for which the strength of the magnetic field inside their intergalactic medium (IGM) is calculated. The estimated values are comparable to that found inside star-forming galaxies, suggesting a dynamical and evolutionary importance of the magnetic field in galaxy groups.

We study a sample of 1,604 moderate-to-high radiative luminosity active galactic nuclei (HLAGN) selected at 3 GHz within the VLA-COSMOS 3 GHz Large Project. These were classified by combining multiple AGN diagnostics: X-ray data, mid-infrared data and broad-band spectral energy distribution fitting. We decompose the total radio 1.4 GHz luminosity ($\mathrm{L_{1.4\ GHz,TOT}}$) into the emission originating from star formation and AGN activity by measuring the excess in $\mathrm{L_{1.4\ GHz,TOT}}$ relative to the infrared-radio correlation of star-forming galaxies. To quantify the excess, for each source we calculate the AGN fraction ($\mathrm{f_{AGN}}$), the fractional contribution of AGN activity to $\mathrm{L_{1.4\ GHz,TOT}}$. The majority of the HLAGN, $(68.0\pm1.5)\%$, are dominated by star-forming processes ($f_{AGN}\leq0.5$), while $(32.0\pm1.5)\%$ are dominated by AGN-related radio emission ($0.5<f_{AGN}\leq1$). We use the AGN-related 1.4 GHz emission to derive the 1.4 GHz AGN luminosity functions of HLAGN. By assuming pure density and pure luminosity evolution models we constrain their cosmic evolution out to $z\sim6$, finding $\mathrm{\Phi^* (z) \propto (1+z)^{(2.64\pm0.10)+(-0.61\pm0.04) z}}$ and $\mathrm{L^* (z) \propto (1+z)^{(3.97\pm0.15) + (-0.92\pm0.06)z}}$. These evolutionary laws show that the number and luminosity density of HLAGN increased from higher redshifts ($z\sim6$) up to a maximum in the redshift range $ 1<z<2.5$, followed by a decline towards local values. By scaling the 1.4 GHz AGN luminosity to kinetic luminosity using the standard conversion, we estimate the kinetic luminosity density as a function of redshift. We compare our result to the semi-analytic models of radio mode feedback finding that this feedback could have played an important role in the context of AGN-host coevolution in HLAGN which show evidence of AGN-related radio emission ($f_{AGN}>0$).

We present very long baseline interferometry (VLBI) observations of 179 radio sources in the COSMOS field with extremely high sensitivity using the Green Bank Telescope (GBT) together with the Very Long Baseline Array (VLBA) (VLBA+GBT) at 1.4 GHz, to explore the faint radio population in the flux density regime of tens of $\mu$Jy. Here, the identification of active galactic nuclei (AGN) is based on the VLBI detection of the source, i.e., it is independent of X-ray or infrared properties. The milli-arcsecond resolution provided by the VLBI technique implies that the detected sources must be compact and have large brightness temperatures, and therefore they are most likely AGN (when the host galaxy is located at z$\geq$0.1). On the other hand, this technique allows us to only positively identify when a radio-active AGN is present, i.e., we cannot affirm that there is no AGN when the source is not detected. For this reason, the number of identified AGN using VLBI should be always treated as a lower limit. We present a catalogue containing the 35 radio sources detected with the VLBA+GBT, 10 of which were not previously detected using only the VLBA. We have constructed the radio source counts at 1.4 GHz using the samples of the VLBA and VLBA+GBT detected sources of the COSMOS field to determine a lower limit for the AGN contribution to the faint radio source population. We found an AGN contribution of >40-75% at flux density levels between 150 $\mu$Jy and 1 mJy. This flux density range is characterised by the upturn of the Euclidean-normalised radio source counts, which implies a contribution of a new population. This result supports the idea that the sub-mJy radio population is composed of a significant fraction of radio-emitting AGN, rather than solely by star-forming galaxies, in agreement with previous studies.

We study the faint radio population using wide-field very long baseline interferometry (VLBI) observations of 2865 known radio sources in the Cosmic Evolution Survey (COSMOS) field. The main objective of the project was to determine where active galactic nuclei (AGN) are present. The combination of number of sources, sensitivity, angular resolution and area covered by this project are unprecedented. We have detected 468 radio sources, expected to be AGNs, with the Very Long Baseline Array (VLBA) at 1.4 GHz. This is, to date, the largest sample assembled of VLBI detected sources in the sub-mJy regime. The input sample was taken from previous observations with the Very Large Array (VLA). We present the catalogue with additional multiwavelength information. We find a detection fraction of 20%, considering only those sources from the input catalogue which were in principle detectable with the VLBA (2361). As a function of redshift, we see no evolution of the detection fraction over the redshift range 0.5<z<3. In addition, we find that faint radio sources typically have a greater fraction of their radio luminosity in a compact core: ~70% of the sub-mJy sources detected with the VLBA have more than half of their total radio luminosity in a VLBI-scale component, whereas this is true for only ~30% of the sources that are brighter than 10 mJy. This suggests that fainter radio sources differ intrinsically from brighter ones. Across our entire sample, we find the predominant morphological classification of the host galaxies of the VLBA detected sources to be early type (57%), although this varies with redshift and at z>1.5 we find that spiral galaxies become the most prevalent (48%). We demonstrate that wide-field VLBI observations, together with new calibration methods such as multi-source self-calibration and mosaicing, result in information which is difficult or impossible to obtain otherwise.

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.

The origin of the radio emission in radio-quiet quasars (RQQs) has been a matter of debate for a long time. It is not well understood whether the emission is caused by star formation in the host galaxy or by black hole activity of the active galactic nuclei (AGN). We shed some light on these questions using the Very Long Baseline Interferometry (VLBI) technique to search for RQQs in the field of the Cosmological Evolution Survey (COSMOS). The extensive multi-wavelength coverage of the field (from radio to X-rays) was used to classify RQQs, and the milli-arcsecond resolution of VLBI provides a direct way to identify AGNs. In a sample of 18 RQQs we detected 3 using the Very Long Baseline Array (VLBA) at 1.4 GHz. In this letter we report for the first time on a sample of RQQs with a measured lower limit on the fraction of radio emission coming from the AGN, thus demonstrating that the radio emission of at least some RQQs is dominated by an AGN.

We are undertaking a search for high-redshift low luminosity Lyman Alpha sources in the SHARDS survey. Among the pre-selected Lyman Alpha sources 2 candidates were spotted, located 3.19 arcsec apart, and tentatively at the same redshift. Here we report on the spectroscopic confirmation with GTC of the Lyman Alpha emission from this pair of galaxies at a confirmed spectroscopic redshifts of z=5.07. Furthermore, one of the sources is interacting/merging with another close companion that looks distorted. Based on the analysis of the spectroscopy and additional photometric data, we infer that most of the stellar mass of these objects was assembled in a burst of star formation 100 Myr ago. A more recent burst (2 Myr old) is necessary to account for the measured Lyman Alpha flux. We claim that these two galaxies are good examples of Lyman Alpha sources undergoing episodic star formation. Besides, these sources very likely constitute a group of interacting Lyman Alpha emitters (LAEs).

The presence of O VI ions can be indicative of plasma temperatures of a few times 10^5 K that is expected in heat conduction layers between the hot shocked stellar wind gas at several 10^6 K and the cooler (~10,000 K) nebular gas of planetary nebulae (PNe). We have used FUSE observations of PNe to search for nebular O VI emission or absorption as a diagnostic of conduction layer to ensure the presence of hot interior gas. Three PNe showing nebular O VI, namely IC 418, NGC 2392, and NGC 6826, have been selected for Chandra observations and diffuse X-ray emission is indeed detected in each of these PNe. Among the three, NGC 2392 has peculiarly high diffuse X-ray luminosity and plasma temperature compared with those expected from its stellar wind's mechanical luminosity and terminal velocity. The limited effects of heat conduction on the plasma temperature of a hot bubble at the low terminal velocity of the stellar wind of NGC 2392 may partially account for its high plasma temperature, but the high X-ray luminosity needs to be powered by processes other than the observed stellar wind, probably caused by the presence of an unseen binary companion of the CSPN of NGC 2392. We have compiled relevant information on the X-ray, stellar, and nebular properties of PNe with a bubble morphology and found that the expectations of bubble models including heat conduction compare favorably with the present X-ray observations of hot bubbles around H-rich CSPNe, but have notable discrepancies for those around H-poor [WR] CSPNe. We note that PNe with more massive central stars can produce hotter plasma and higher X-ray surface brightness inside central hot bubbles.

The planetary nebula (PN) A30 is believed to have undergone a very late thermal pulse resulting in the ejection of knots of hydrogen-poor material. Using HST images we have detected the angular expansion of these knots and derived an age of 850+280-150 yr. To investigate the spectral and spatial properties of the soft X-ray emission detected by ROSAT, we have obtained Chandra and XMM-Newton observations of A30. The X-ray emission from A30 can be separated into two components: a point-source at the central star and diffuse emission associated with the hydrogen-poor knots and the cloverleaf structure inside the nebular shell. To help us assess the role of the current stellar wind in powering this X-ray emission, we have determined the stellar parameters of the central star of A 30 using a non-LTE model fit to its optical and UV spectrum. The spatial distribution and spectral properties of the diffuse X-ray emission is suggestive that it is generated by the post-born-again and present fast stellar winds interacting with the hydrogen-poor ejecta of the born-again event. This emission can be attributed to shock-heated plasma, as the hydrogen-poor knots are ablated by the stellar winds, under which circumstances the efficient mass-loading of the present fast stellar wind raises its density and damps its velocity to produce the observed diffuse soft X-rays. Charge transfer reactions between the ions of the stellar winds and material of the born-again ejecta has also been considered as a possible mechanism for the production of diffuse X-ray emission, and upper limits on the expected X-ray production by this mechanism have been derived. The origin of the X-ray emission from the central star of A 30 is puzzling: shocks in the present fast stellar wind and photospheric emission can be ruled out, while the development of a new, compact hot bubble confining the fast stellar wind seems implausible.

One key feature of the interacting stellar winds model of the formation of planetary nebulae (PNe) is the presence of shock-heated stellar wind confined in the central cavities of PNe. This so-called hot bubble should be detectable in X-rays. Here we present XMM-Newton observations of NGC 3242, a multiple-shell PN whose shell morphology is consistent with the interacting stellar winds model. Diffuse X-ray emission is detected within its inner shell with a plasma temperature ~2.35\times10^6 K and an intrinsic X-ray luminosity ~2\times10^30 ergs s^(-1) at the adopted distance of 0.55 kpc. The observed X-ray temperature and luminosity are in agreement with "ad-hoc" predictions of models including heat conduction. However, the chemical abundances of the X-ray-emitting plasma seem to imply little evaporation of cold material into the hot bubble, whereas the thermal pressure of the hot gas is unlikely to drive the nebular expansion as it is lower than that of the inner shell rim. These inconsistencies are compounded by the apparent large filling factor of the hot gas within the central cavity of NGC 3242. Subject headings: planetary nebulae: individual (NGC 3242)