Aims. We study the coevolution between the black hole accretion rate (BHAR) and the star formation rate (SFR) in different phases of galaxy life: main-sequence star-forming galaxies, quiescent galaxies, and starburst galaxies at different cosmic epochs. Methods. We exploited the unique combination of depth and area in the COSMOS field and took advantage of the X-ray data from the Chandra COSMOS-Legacy survey and the extensive multiwavelength ancillary data presented in the COSMOS2015 catalog, including in particular the UVista Ultra-deep observations. These large datasets allowed us to perform an X-ray stacking analysis and combine it with detected sources in a broad redshift interval (0.1 <  z <  3.5) with unprecedented statistics for normal star-forming, quiescent, and starburst galaxies. The X-ray luminosity was used to predict the black holeAR, and a similar stacking analysis on far-infrared Herschel maps was used to measure the corresponding obscured SFR for statistical samples of sources in different redshifts and stellar mass bins. Results. We focus on the evolution of the average SFR-stellar mass (M*) relation and compare it with the BHAR-M* relation. This extends previous works that pointed toward the existence of almost linear correlations in both cases. We find that the ratio between BHAR and SFR does not evolve with redshift, although it depends on stellar mass. For the star-forming populations, this dependence on M* has a logarithmic slope of ∼0.6 and for the starburst sample, the slope is ∼0.4. These slopes are both at odds with quiescent sources, where the dependence remains constant (log(BHAR/SFR) ∼ −3.4). By studying the specific BHAR and specific SFR, we find signs of downsizing for M* and black hole mass (MBH) in galaxies in all evolutionary phases. The increase in black hole mass-doubling timescale was particularly fast for quiescents, whose super-massive black holes grew at very early times, while accretion in star-forming and starburst galaxies continued until more recent times. Conclusions. Our results support the idea that the same physical processes feed and sustain star formation and black hole accretion in star-forming galaxies while the starburst phase plays a lesser role in driving the growth of the supermassive black holes, especially at high redshift. Our integrated estimates of the M* − MBH relation at all redshifts are consistent with independent determinations of the local M* − MBH relation for samples of active galactic nuclei. This adds key evidence that the evolution in the BHAR/SFR is weak and its normalization is relatively lower than that of local dynamical M* − MBH relations.

Determining the frequency and duration of changing–look (CL) active galactic nuclei (AGNs) phenomena, where the optical broad emission lines appear or disappear, is crucial to understand the evolution of the accretion flow around supermassive black holes. We present a strategy to select new CL candidates starting from a spectroscopic type 2 AGNs sample and searching for current type 1 photometric variability. We use the publicly available Zwicky Transient Facility (ZTF) alert stream and the Automatic Learning for the Rapid Classification of Events (ALeRCE) light curve classifier to produce a list of CL candidates with a highly automated algorithm, resulting in 60 candidates. Visual inspection reduced the sample to 30. We performed new spectroscopic observations of six candidates of our clean sample, without further refinement, finding the appearance of clear broad Balmer lines in four of them and tentative evidence of type changes in the remaining two, which suggests a promising success rate of ≥66 per cent for this CL selection method.

The scarce optical variability studies in spectrally classified Type 2 active galactic nuclei (AGNs) have led to the discovery of anomalous objects that are incompatible with the simplest unified models (UM). This paper focuses on the exploration of different variability features that allows to separate between obscured, Type 2 AGNs, and the variable, unobscured Type 1s. We analyse systematically the Zwicky Transient Facility, 2.5 years long light curves of ∼ 15000 AGNs from the Sloan Digital Sky Survey Data Release 16, which are generally considered Type 2s due to the absence of strong broad emission lines (BELs). Consistently with the expectations from the UM, the variability features are distributed differently for distinct populations, with spectrally classified weak Type 1s showing 1 order of magnitude larger variances than the Type 2s. We find that the parameters given by the damped random walk model leads to broader Hα equivalent width for objects with τg &gt; 16 d and long term structure function SF∞, g&gt; 0.07 mag. By limiting the variability features, we find that ∼ 11 per cent of Type 2 sources show evidence for optical variations. A detailed spectral analysis of the most variable sources (∼1 per cent of the Type 2 sample) leads to the discovery of misclassified Type 1s with weak BELs and changing-state candidates. This work presents one of the largest systematic investigations of Type 2 AGN optical variability to date, in preparation for future large photometric surveys.

We present second epoch optical spectra for 30 changing-look (CL) candidates found by searching for Type-1 optical variability in a sample of active galactic nuclei (AGNs) spectroscopically classified as Type 2. We use a random-forest-based light curve classifier and spectroscopic follow-up, confirming 50 per cent of candidates as turning-on CLs. In order to improve this selection method and to better understand the nature of the not-confirmed CL candidates, we perform a multi-wavelength variability analysis including optical, mid-infrared (MIR) and X-ray data, and compare the results from the confirmed and not-confirmed CLs identified in this work. We find that most of the not-confirmed CLs are consistent with weak Type 1s dominated by host-galaxy contributions, showing weaker optical and MIR variability. On the contrary, the confirmed CLs present stronger optical fluctuations and experience a long (from five to ten years) increase in their MIR fluxes and the colour W1–W2 over time. In the 0.2–2.3 keV band, at least four out of 11 CLs with available SRG/eROSITA detections have increased their flux in comparison with archival upper limits. These common features allow us to select the most promising CLs from our list of candidates, leading to nine sources with similar multi-wavelength photometric properties to our CL sample. The use of machine learning algorithms with optical and MIR light curves will be very useful to identify CLs in future large-scale surveys.