A closer look at the deep radio sky: Multi-component radio sources at 3 GHz VLA-COSMOS

Context. Given the unprecedented depth achieved in current large radio surveys, we are starting to probe populations of radio sources that have not been studied in the past. However, identifying and categorising these objects, differing in size, shape and physical properties, is becoming a more difficult task.Aims. 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(-1) rms), i.e. the radio sources which are composed of two or more radio blobs.Methods. 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 Ultra Deep Survey with the VISTA telescope (UltraVISTA) stacked mosaic available for COSMOS.Results. We have identified 67 multi-component radio sources at 3 GHz: 58 sources with active galactic nucleus (AGN) powered radio emission and nine star-forming galaxies. We report eight 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. We find that multi-component objects at 3 GHz VLA-COSMOS inhabit mainly massive galaxies (>10(10.5)M(circle dot)). 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. Furthermore, we provide detailed descriptions of the objects and find that amongst the AGN there are two head-tail, ten corelobe, nine wide-angle-tail (WAT), eight double-double or Z-/X-shaped, three 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. Additionally, we report a large number (32 out of 58) of disturbed/bent multi-component AGN, 18 of which do not lie within X-ray groups in COSMOS (redshift range 0.08 Conclusion. The high angular resolution and sensitivity of the 3 GHz VLA-COSMOS data set give us the opportunity to identify peculiar radio structures and sub-structures of multi-component objects, and relate them to physical phenomena such as AGN or star-forming galaxies. This study illustrates the complexity of the mu Jy radio-source population; at the sensitivity and resolution of 3 GHz VLA-COSMOS, the radio structures of AGN and SFG both emitting radio continuum emission, become comparable in the absence of clear, symmetrical jets. Thus, disentangling the AGN and SFG contributions using solely radio observations can be misleading in a number of cases. This has implications for future surveys, such as those done by square kilometre array (SKA) and precursors, which will identify hundreds of thousands of multi-component objects.</div

Gaia Focused Product Release: Radial velocity time series of long-period variables

The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged)Comment: 36 pages, 38 figure

Gaia Focused Product Release: Asteroid orbital solution: Properties and assessment

Context. We report the exploitation of a sample of Solar System observations based on data from the third Gaia Data Release (Gaia DR3) of nearly 157 000 asteroids. It extends the epoch astrometric solution over the time coverage planned for the Gaia DR4, which is not expected before the end of 2025. This data set covers more than one full orbital period for the vast majority of these asteroids. The orbital solutions are derived from the Gaia data alone over a relatively short arc compared to the observation history of many of these asteroids. Aims. The work aims to produce orbital elements for a large set of asteroids based on 66 months of accurate astrometry provided by Gaia and to assess the accuracy of these orbital solutions with a comparison to the best available orbits derived from independent observations. A second validation is performed with accurate occultation timings. Methods. We processed the raw astrometric measurements of Gaia to obtain astrometric positions of moving objects with 1D sub-mas accuracy at the bright end. For each asteroid that we matched to the data, an orbit fitting was attempted in the form of the best fit of the initial conditions at the median epoch. The force model included Newtonian and relativistic accelerations to derive the observation equations, which were solved with a linear least-squares fit. Results. Orbits are provided in the form of state vectors in the International Celestial Reference Frame for 156 764 asteroids, including near-Earth objects, main-belt asteroids, and Trojans. For the asteroids with the best observations, the (formal) relative uncertainty σa/a is better than 10−10. Results are compared to orbits available from the Jet Propulsion Laboratory and MPC. Their orbits are based on much longer data arcs, but from positions of lower quality. The relative differences in semi-major axes have a mean of 5 × 10−10 and a scatter of 5 × 10−9

Gaia Focused Product Release: Radial velocity time series of long-period variables

Context. The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity data planned with Data Release 4, this Focused Product Release (FPR) provides radial-velocity time series for a selection of LPV candidates with high-quality observations. Aims. We describe the production and content of the Gaia catalog of LPV radial-velocity time series, and the methods used to compute the variability parameters published as part of the Gaia FPR. Methods. Starting from the DR3 catalog of LPV candidates, we applied several filters to construct a sample of sources with high-quality radial-velocity measurements. We modeled their radial-velocity and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the radial-velocity period and at least one of the G, GBP, or GRP photometric periods. Results. The catalog includes radial-velocity time series and variability parameters for 9614 sources in the magnitude range 6 ≲ G/mag ≲ 14, including a flagged top-quality subsample of 6093 stars whose radial-velocity periods are fully compatible with the values derived from the G, GBP, and GRP photometric time series. The radial-velocity time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great majority of the sources (88%) as genuine LPV candidates, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% of the catalog consists of candidate ellipsoidal binaries. Quality checks against radial velocities available in the literature show excellent agreement. We provide some illustrative examples and cautionary remarks. Conclusions. The publication of radial-velocity time series for almost ten thousand LPV candidates constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog

Gaia Focused Product Release: Spatial distribution of two diffuse interstellar bands

Diffuse interstellar bands (DIBs) are absorption features seen in optical and infrared spectra of stars and extragalactic objects that are probably caused by large and complex molecules in the galactic interstellar medium (ISM). Here we investigate the Galactic distribution and properties of two DIBs identified in almost six million stellar spectra collected by the Gaia Radial Velocity Spectrometer. These measurements constitute a part of the Gaia Focused Product Release to be made public between the Gaia DR3 and DR4 data releases. In order to isolate the DIB signal from the stellar features in each individual spectrum, we identified a set of 160 000 spectra at high Galactic latitudes (|b| ⩾ 65°) covering a range of stellar parameters which we consider to be the DIB-free reference sample. Matching each target spectrum to its closest reference spectra in stellar parameter space allowed us to remove the stellar spectrum empirically, without reference to stellar models, leaving a set of six million ISM spectra. Using the star’s parallax and sky coordinates, we then allocated each ISM spectrum to a voxel (VOlume piXEL) on a contiguous three-dimensional grid with an angular size of 1.8° (level 5 HEALPix) and 29 unequally sized distance bins. Identifying the two DIBs at 862.1 nm (λ862.1) and 864.8 nm (λ864.8) in the stacked spectra, we modelled their shapes and report the depth, central wavelength, width, and equivalent width (EW) for each, along with confidence bounds on these measurements. We then explored the properties and distributions of these quantities and compared them with similar measurements from other surveys. Our main results are as follows: (1) the strength and spatial distribution of the DIB λ862.1 are very consistent with what was found in Gaia DR3, but for this work we attained a higher signal-to-noise ratio in the stacked spectra to larger distances, which allowed us to trace DIBs in the outer spiral arm and beyond the Scutum–Centaurus spiral arm; (2) we produced an all-sky map below ±65° of Galactic latitude to ~4000 pc of both DIB features and their correlations; (3) we detected the signals of DIB λ862.1 inside the Local Bubble (≲200 pc); and (4) there is a reasonable correlation with the dust reddening found from stellar absorption and EWs of both DIBs with a correlation coefficient of 0.90 for λ862.1 and 0.77 for λ864.8

Gaia Focused Product Release: Sources from Service Interface Function image analysis

Context. Gaia’s readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri (ω Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Aims. Our aim is to improve the completeness of the Gaia source inventory in a very dense region in the sky, ω Cen. Methods. An adapted version of Gaia’s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. These source detections were clustered and assigned to new SIF CF or existing Gaia sources by Gaia s cross-match software. For the new sources, astrometry was calculated using the Astrometric Global Iterative Solution software, and photometry was obtained in the Gaia DR3 reference system. We validated the results by comparing them to the public Gaia DR3 catalogue and external Hubble Space Telescope data. Results. With this Focused Product Release, 526 587 new sources have been added to the Gaia catalogue in ω Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal Gaia sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. Conclusions. This first SIF CF data publication already adds great value to the Gaia catalogue. It demonstrates what to expect for the fourth Gaia catalogue, which will contain additional sources for all nine SIF CF regions

Infrared-radio relation in the local Universe

Context. The Square Kilometer Array (SKA) is expected to detect high-redshift galaxies with star formation rates (SFRs) up to two orders of magnitude lower than Herschel surveys and will thus boost the ability of radio astronomy to study extragalactic sources. The tight infrared-radio correlation offers the possibility of using radio emission as a dust-unobscured star formation diagnostic. However, the physics governing the link between radio emission and star formation is poorly understood, and recent studies have pointed to differences in the exact calibration required when radio is to be used as a star formation tracer. Aims. We improve the calibration of the relation of the local radio luminosity–SFR and to test whether there are nonlinearities in it. Methods. We used a sample of Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) sources and investigated their radio luminosity, which was derived using the NRAO VLA Sky Survey (NVSS) and Faint Images of the Radio Sky at Twenty-cm (FIRST) maps. We stacked the bins of infrared luminosity and SFR and accounted for bins with no detections in the stacked images using survival analysis fitting. This approach was tested using Monte Carlo simulations. Results. After removing sources from the sample that have excess radio emission, which is indicative of nuclear radio activity, we found no deviations from linearity of the mean relations between radio luminosity and either SFR or infrared luminosity. Concluisions. We analyzed the link between radio emission and SFR or infrared luminosity using a local sample of star-forming galaxies without evidence of nuclear radio activity and found no deviations from linearity, although our data are also consistent with the small nonlinearity reported by some recent analyses. The normalizations of these relations are intermediate between those reported by earlier works

The XXL Survey: XLIII. The quasar radio loudness dichotomy exposed via radio luminosity functions obtained by combining results from COSMOS and XXL-S X-ray selected quasars

We studied a sample of 274 radio and X-ray selected quasars (XQSOs) detected in the COSMOS and XXL-S radio surveys at 3 GHz and 2.1 GHz, respectively. This sample was identified by adopting a conservative threshold in X-ray luminosity, LX [2-10 keV] ≥ 1044 erg s-1, selecting only the most powerful quasars. A number of previous studies on the origin of radio emission in type-1 quasars have focused on the radio loudness distributions, some claiming to have found evidence for bimodality, pointing toward the existence of two physically different mechanisms for the radio emission. Using available multiwavelength data, we examined various criteria for the selection of radio-loud (RL) and radio-quiet (RQ) XQSOs and found that the number of RL/RQ XQSOs changes significantly depending on the chosen criterion. This discrepancy arises due to the different criteria tracing different physical processes and due to the fact that our sample was selected from flux-limited radio and X-ray surveys. Another approach to study the origin of radio emission in XQSOs is via their radio luminosity functions (RLF). We constructed the XQSO 1.4 GHz RLFs in six redshift bins at 0:5 ≤ z ≤ 3:75. The lower-1.4 GHz luminosity end shows a higher normalization than expected only from AGN contribution in all studied redshift bins. We found that the so-called "bump"is mostly dominated by emission due to star-forming processes within the host galaxies of XQSOs. As expected, AGN-related radio emission is the dominant contribution at the higher-luminosity end of RLF. To study the evolution of the XQSO RLF, we used a combination of analytic forms from the literature to constrain the "bump"due to star formation and the higher-luminosity AGN part of the RLF. We defined two 1.4 GHz luminosity thresholds, Lth;SF and Lth;AGN, below and above which more than 80% of sources contributing to the RLF are dominated by star formation and AGN-related activity, respectively. The two thresholds evolve with redshift, which is most likely driven by the strong evolution of star formation rates of the XQSO host galaxies. We found that both the lower and higher luminosity ends evolve significantly in density, while their luminosity evolution parameters are consistent with being constant. We found that the lower-luminosity end evolves both in density and luminosity, while the higher-luminosity end evolves significantly only in density. Our results expose the dichotomy of the origin of radio emission: while the higher-luminosity end of the XQSO RLF is dominated by AGN activity, the lower-luminosity end is dominated by the star formation-related processes

The XXL Survey: XLI. Radio AGN luminosity functions based on the GMRT 610 MHz continuum observations

We study the space density evolution of active galactic nuclei (AGN) using the 610 MHz radio survey of the XXL-North field, performed with the Giant Metrewave Radio Telescope. The survey covers an area of 30.4 deg2, with a beamsize of 6.5 arcsec. The survey is divided into two parts, one covering an area of 11.9 deg2 with 1σ rms noise of 200 μJy beam−1 and the other spanning 18.5 deg2 with rms noise of 45 μJy beam−1. We extracted the catalog of radio components above 7σ. The catalog was cross-matched with a multi-wavelength catalog of the XXL-North field (covering about 80% of the radio XXL-North field) using a likelihood ratio method, which determines the counterparts based on their positions and their optical properties. The multi-component sources were matched visually with the aid of a computer code: Multi-Catalog Visual Cross-Matching. A flux density cut above 1 mJy selects AGN hosts with a high purity in terms of star formation contamination based on the available source counts. After cross-matching and elimination of observational biases arising from survey incompletenesses, the number of remaining sources was 1150. We constructed the rest-frame 1.4 GHz radio luminosity functions of these sources using the maximum volume method. This survey allows us to probe luminosities of 23 ≲ log(L1.4 GHz[W Hz−1]) ≲ 28 up to redshifts of z ≈ 2.1. Our results are consistent with the results from the literature in which AGN are comprised of two differently evolving populations, where the high luminosity end of the luminosity functions evolves more strongly than the low-luminosity end