The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, in synergy with the SUDARE survey, is a deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). The observations for the CDFS field comprise about 4.38 deg$^2$ down to $r\sim26$ mag. The total on-sky time spans over four years in this field, distributed over four adjacent sub-fields. In this paper, we use the multi-epoch $r$-band imaging data to measure the variability of the detected objects and search for transients. We perform careful astrometric and photometric calibrations and point spread function (PSF) modeling. A new method, referring to as differential running-average photometry, is proposed to measure the light curves of the detected objects. With the method, the difference of PSFs between different epochs can be reduced, and the background fluctuations are also suppressed. Detailed uncertainty analysis and detrending corrections on the light curves are performed. We visually inspect the light curves to select variable objects, and present some objects with interesting light curves. Further investigation of these objects in combination with multi-band data will be presented in our forthcoming paper.

Variability has proven to be a powerful tool to detect active galactic nuclei (AGN) in multi-epoch surveys. The new-generation facilities expected to become operational in the next few years will mark a new era in time-domain astronomy and their wide-field multi-epoch campaigns will favor extensive variability studies. We present our analysis of AGN variability in the second half of the VST survey of the Wide Chandra Deep Field South (W-CDF-S), performed in the r band and covering a 2 sq. deg. area. The analysis complements a previous work, in which the first half of the area was investigated. We provide a reliable catalog of variable AGN candidates, which will be critical targets in future variability studies. We selected a sample of optically variable sources and made use of infrared data from the Spitzer mission to validate their nature by means of color-based diagnostics. We obtain a sample of 782 AGN candidates among which 12 are classified as supernovae, 54 as stars, and 232 as AGN. We estimate a contamination $\lesssim 20\%$ and a completeness $\sim 38\%$ with respect to mid-infrared selected samples.

The analysis of the variability of active galactic nuclei (AGNs) at different wavelengths and the study of possible correlations among different spectral windows are nowadays a major field of inquiry. Optical variability has been largely used to identify AGNs in multivisit surveys. The strength of a selection based on optical variability lies in the chance to analyze data from surveys of large sky areas by ground-based telescopes. However the effectiveness of optical variability selection, with respect to other multiwavelength techniques, has been poorly studied down to the depth expected from next generation surveys. Here we present the results of our r-band analysis of a sample of 299 optically variable AGN candidates in the VST survey of the COSMOS field, counting 54 visits spread over three observing seasons spanning > 3 yr. This dataset is > 3 times larger in size than the one presented in our previous analysis (De Cicco et al. 2015), and the observing baseline is ~8 times longer. We push towards deeper magnitudes (r(AB) ~23.5 mag) compared to past studies; we make wide use of ancillary multiwavelength catalogs in order to confirm the nature of our AGN candidates, and constrain the accuracy of the method based on spectroscopic and photometric diagnostics. We also perform tests aimed at assessing the relevance of dense sampling in view of future wide-field surveys. We demonstrate that the method allows the selection of high-purity (> 86%) samples. We take advantage of the longer observing baseline to achieve great improvement in the completeness of our sample with respect to X-ray and spectroscopically confirmed samples of AGNs (59%, vs. ~15% in our previous work), as well as in the completeness of unobscured and obscured AGNs. The effectiveness of the method confirms the importance to develop future, more refined techniques for the automated analysis of larger datasets.

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is a Guaranteed Time program carried out with the ESO/VST telescope to provide deep optical imaging over two 4 deg$^2$ patches of the sky centred on the CDFS and ES1 pointings. We present the cosmic shear measurement over the 4 deg$^2$ covering the CDFS region in the $r$-band using LensFit. Each of the four tiles of 1 deg$^2$ has more than one hundred exposures, of which more than 50 exposures passed a series of image quality selection criteria for weak lensing study. The $5\sigma$ limiting magnitude in $r$- band is 26.1 for point sources, which is $\sim$1 mag deeper than other weak lensing survey in the literature (e.g. the Kilo Degree Survey, KiDS, at VST). The photometric redshifts are estimated using the VOICE $u,g,r,i$ together with near-infrared VIDEO data $Y,J,H,K_s$. The mean redshift of the shear catalogue is 0.87, considering the shear weight. The effective galaxy number density is 16.35 gal/arcmin$^2$, which is nearly twice the one of KiDS. The performance of LensFit on such a deep dataset was calibrated using VOICE-like mock image simulations. Furthermore, we have analyzed the reliability of the shear catalogue by calculating the star-galaxy cross-correlations, the tomographic shear correlations of two redshift bins and the contaminations of the blended galaxies. As a further sanity check, we have constrained cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. For a flat $\Lambda$CDM model we have obtained $\Sigma_8$ = $\sigma_8(\Omega_m/0.3)^{0.5}$ = $0.68^{+0.11}_{-0.15}$.

The VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey is proposed to obtain deep optical $ugri$ imaging of the CDFS and ES1 fields using the VLT Survey Telescope (VST). At present, the observations for the CDFS field have been completed, and comprise in total about 4.9 deg$^2$ down to $r_\mathrm{AB}$$\sim$26 mag. In the companion paper by Fu et al. (2018), we present the weak lensing shear measurements for $r$-band images with seeing $\le$ 0.9 arcsec. In this paper, we perform image simulations to calibrate possible biases of the measured shear signals. Statistically, the properties of the simulated point spread function (PSF) and galaxies show good agreements with those of observations. The multiplicative bias is calibrated to reach an accuracy of $\sim$3.0%. We study the bias sensitivities to the undetected faint galaxies and to the neighboring galaxies. We find that undetected galaxies contribute to the multiplicative bias at the level of $\sim$0.3%. Further analysis shows that galaxies with lower signal-to-noise ratio (SNR) are impacted more significantly because the undetected galaxies skew the background noise distribution. For the neighboring galaxies, we find that although most have been rejected in the shape measurement procedure, about one third of them still remain in the final shear sample. They show a larger ellipticity dispersion and contribute to $\sim$0.2% of the multiplicative bias. Such a bias can be removed by further eliminating these neighboring galaxies. But the effective number density of the galaxies can be reduced considerably. Therefore efficient methods should be developed for future weak lensing deep surveys.

Broad absorption lines (BALs) in the spectra of quasi-stellar objects (QSOs) originate from outflowing winds along our line of sight; winds are thought to originate from the inner regions of the QSO accretion disk, close to the central supermassive black hole (SMBH). Winds likely play a role in galaxy evolution and aid the accretion mechanism onto the SMBH. BAL equivalent widths can change on typical timescales from months to years; such variability is generally attributed to changes in the covering factor and/or in the ionization level of the gas. We investigate BAL variability, focusing on BAL disappearance. We analyze multi-epoch spectra of more than 1500 QSOs -the largest sample ever used for such a study- observed by different programs from the Sloan Digital Sky Survey-I/II/III (SDSS), and search for disappearing C IV BALs. The spectra rest-frame time baseline ranges from 0.28 to 4.9 yr; the source redshifts range from 1.68 to 4.27. We detect 73 disappearing BALs in the spectra of 67 sources. This corresponds to 3.9% of disappearing BALs, and 5.1% of our BAL QSOs exhibit at least one disappearing BAL. We estimate the average lifetime of a BAL along our line of sight (~ 80-100 yr), which appears consistent with the accretion disk orbital time at distances where winds are thought to originate. We inspect properties of the disappearing BALs and compare them to the properties of our main sample. We also investigate the existence of a correlation in the variability of multiple troughs in the same spectrum, and find it persistent at large velocity offsets between BAL pairs, suggesting that a mechanism extending on a global scale is necessary to explain the phenomenon. We select a more reliable sample of disappearing BALs following Filiz Ak et al. (2012), where a subset of our sample was analyzed, and compare the findings from the two works, obtaining generally consistent results.

We describe the observing strategy, data reduction tools and early results of a supernova (SN) search project, named SUDARE, conducted with the ESO VST telescope aimed at measuring the rate of the different types of SNe in the redshift range 0.2<z<0.8. The search was performed in two of the best-studied extragalactic fields, CDFS and COSMOS, for which a wealth of ancillary data are available in the literature or public archives. (abridged) We obtained a final sample of 117 SNe, most of which are SNIa (57%) and the remaining core collapse events of which 44% type II, 22% type IIn and 34% type Ib/c. In order to link the transients, we built a catalog of ~1.3x10^5 galaxies in the redshift range 0<z<1 with a limiting magnitude K_AB=23.5 mag. We measured the SN rate per unit volume for SN Ia and core collapse SNe in different bin of redshifts. The values are consistent with other measurements from the literature. The dispersion of the rate measurements for SNe Ia is comparable with the scatter of the theoretical tracks for single (SD) and double degenerate (DD) binary systems models, therefore the data do not allow to disentangle among the two different progenitor scenarios. However, we may notice that among the three tested models, SD and two flavours of DD, either with a steep (DDC) or a wide (DDW) delay time distribution, the SD gives a better fit across the whole redshift range whereas the DDC better matches the steep rise up to redshift ~1.2. The DDW appears instead less favoured. The core collapse SN rate is fully consistent, unlike recent claims, with the prediction based on recent estimates of the star formation history, and standard progenitor mass range.

We present the VST Optical Imaging of the CDFS and ES1 Fields (VOICE) Survey, a VST INAF Guaranteed Time program designed to provide optical coverage of two 4 deg$^2$ cosmic windows in the Southern hemisphere. VOICE provides the first, multi-band deep optical imaging of these sky regions, thus complementing and enhancing the rich legacy of longer-wavelength surveys with VISTA, Spitzer, Herschel and ATCA available in these areas and paving the way for upcoming observations with facilities such as the LSST, MeerKAT and the SKA. VOICE exploits VST's OmegaCAM optical imaging capabilities and completes the reduction of WFI data available within the ES1 fields as part of the ESO-Spitzer Imaging Extragalactic Survey (ESIS) program providing $ugri$ and $uBVR$ coverage of 4 and 4 deg$^2$ areas within the CDFS and ES1 field respectively. We present the survey's science rationale and observing strategy, the data reduction and multi-wavelength data fusion pipeline. Survey data products and their future updates will be released at http://www.mattiavaccari.net/voice/ and on CDS/VizieR.

This is the second paper of a series in which we present measurements of the Supernova (SN) rates from the SUDARE survey. In this paper, we study the trend of the SN rates with the intrinsic colours, the star formation activity and the mass of the parent galaxies. We have considered a sample of about 130000 galaxies and a SN sample of about 50 events. We found that the SN Ia rate per unit mass is higher by a factor of six in the star-forming galaxies with respect to the passive galaxies. The SN Ia rate per unit mass is also higher in the less massive galaxies that are also younger. These results suggest a distribution of the delay times (DTD) less populated at long delay times than at short delays. The CC SN rate per unit mass is proportional to both the sSFR and the galaxy mass. The trends of the Type Ia and CC SN rates as a function of the sSFR and the galaxy mass that we observed from SUDARE data are in agreement with literature results at different redshifts. The expected number of SNe Ia is in agreement with the observed one for all four DTD models considered both in passive and star-forming galaxies so we can not discriminate between different progenitor scenarios. The expected number of CC SNe is higher than the observed one, suggesting a higher limit for the minimum progenitor mass. We also compare the expected and observed trends of the SN Ia rate with the intrinsic U - J colour of the parent galaxy, assumed as a tracer of the age distribution. While the slope of the relation between the SN Ia rate and the U - J color in star-forming galaxies can be reproduced well by all four DTD models considered, only the steepest of them is able to account for the rates and colour in star-forming and passive galaxies with the same value of the SN Ia production efficiency.

This work makes use of the VST observations to select variable sources. We use also the IR photometry, SED fitting and X-ray information where available to confirm the nature of the AGN candidates. The IR data, available over the full survey area, allow to confirm the consistency of the variability selection with the IR color selection method, while the detection of variability may prove useful to detect the presence of an AGN in IR selected starburst galaxies.

One of the most peculiar characteristics of Active Galactic Nuclei (AGN) is their variability over all wavelengths. This property has been used in the past to select AGN samples and is foreseen to be one of the detection techniques applied in future multi-epoch surveys, complementing photometric and spectroscopic methods. In this paper, we aim to construct and characterise an AGN sample using a multi-epoch dataset in the r band from the SUDARE-VOICE survey. Our work makes use of the VST monitoring program of an area surrounding the Chandra Deep Field South to select variable sources. We use data spanning a six month period over an area of 2 square degrees, to identify AGN based on their photometric variability. The selected sample includes 175 AGN candidates with magnitude r < 23 mag. We distinguish different classes of variable sources through their lightcurves, as well as X-ray, spectroscopic, SED, optical and IR information overlapping with our survey. We find that 12% of the sample (21/175) is represented by SN. Of the remaining sources, 4% (6/154) are stars, while 66% (102/154) are likely AGNs based on the available diagnostics. We estimate an upper limit to the contamination of the variability selected AGN sample of about 34%, but we point out that restricting the analysis to the sources with available multi-wavelength ancillary information, the purity of our sample is close to 80% (102 AGN out of 128 non-SN sources with multi-wavelength diagnostics). Our work thus confirms the efficiency of the variability selection method in agreement with our previous work on the COSMOS field; in addition we show that the variability approach is roughly consistent with the infrared selection.

Optical variability has proven to be an effective way of detecting AGNs in imaging surveys, lasting from weeks to years. In the present work we test its use as a tool to identify AGNs in the VST multi-epoch survey of the COSMOS field, originally tailored to detect supernova events. We make use of the multi-wavelength data provided by other COSMOS surveys to discuss the reliability of the method and the nature of our AGN candidates. Our selection returns a sample of 83 AGN candidates; based on a number of diagnostics, we conclude that 67 of them are confirmed AGNs (81% purity), 12 are classified as supernovae, while the nature of the remaining 4 is unknown. For the subsample of AGNs with some spectroscopic classification, we find that Type 1 are prevalent (89%) compared to Type 2 AGNs (11%). Overall, our approach is able to retrieve on average 15% of all AGNs in the field identified by means of spectroscopic or X-ray classification, with a strong dependence on the source apparent magnitude. In particular, the completeness for Type 1 AGNs is 25%, while it drops to 6% for Type 2 AGNs. The rest of the X-ray selected AGN population presents on average a larger r.m.s. variability than the bulk of non variable sources, indicating that variability detection for at least some of these objects is prevented only by the photometric accuracy of the data. We show how a longer observing baseline would return a larger sample of AGN candidates. Our results allow us to assess the usefulness of this AGN selection technique in view of future wide-field surveys.