Verifying the fully kinematic nature of the cosmic microwave background (CMB) dipole is of fundamental importance in cosmology. In the standard cosmological model with the Friedman-Lemaitre-Robertson-Walker (FLRW) metric from the inflationary expansion the CMB dipole should be entirely kinematic. Any non-kinematic CMB dipole component would thus reflect the preinflationary structure of spacetime probing the extent of the FLRW applicability. Cosmic backgrounds from galaxies after the matter-radiation decoupling, should have kinematic dipole component identical in velocity with the CMB kinematic dipole. Comparing the two can lead to isolating the CMB non-kinematic dipole. It was recently proposed that such measurement can be done using the near-IR cosmic infrared background (CIB) measured with the currently operating Euclid telescope, and later with Roman. The proposed method reconstructs the resolved CIB, the Integrated Galaxy Light (IGL), from Euclid's Wide Survey and probes its dipole, with a kinematic component amplified over that of the CMB by the Compton-Getting effect. The amplification coupled with the extensive galaxy samples forming the IGL would determine the CIB dipole with an overwhelming signal/noise, isolating its direction to sub-degree accuracy. We develop details of the method for Euclid's Wide Survey in 4 bands spanning 0.6 to 2 mic. We isolate the systematic and other uncertainties and present methodologies to minimize them, after confining the sample to the magnitude range with negligible IGL/CIB dipole from galaxy clustering. These include the required star-galaxy separation, accounting for the extinction correction dipole using the method newly developed here achieving total separation, accounting for the Earth's orbital motion and other systematic effects. (Abridged)

Structure on sub-galactic scales provides important tests of galaxy formation models and the nature of dark matter. However, such objects are typically too faint to provide robust mass constraints. Here, we report the discovery of an extremely low-mass object detected via its gravitational perturbation to a thin lensed arc observed with milli-arcsecond-resolution very long baseline interferometry (VLBI). The object was identified using a non-parametric gravitational imaging technique and confirmed using independent parametric modelling. It contains a mass of $m_{\rm 80}=(1.13 \pm 0.04)\times 10^6{M_\odot}$ within a projected radius of 80 parsecs at an assumed redshift of 0.881. This detection is extremely robust and precise, with a statistical significance of 26$\sigma$, a 3.3 per cent fractional uncertainty on $m_{\rm 80}$, and an astrometric uncertainty of 194 $\mu$as. This is the lowest-mass object known to us, by two orders of magnitude, to be detected at a cosmological distance by its gravitational effect. This work demonstrates the observational feasibility of using gravitational imaging to probe the million-solar-mass regime far beyond our local Universe.

AGN jet structure and stability remain uncertain; recollimation shocks are linked to morphology and variability, but the role of downstream instabilities is still unclear. We aim to investigate how jet magnetization and other physical parameters influence the development of instabilities beyond the first recollimation shock. In particular, we focus on identifying the conditions under which the centrifugal instability (CFI) is effective. We perform high-resolution 2D and 3D simulations using the relativistic magnetohydrodynamics code PLUTO. The jets are initialized with a conical geometry and propagate into an ambient medium, and we follow by axisymmetric simulations how they evolve towards a steady-state. In 2D we explore a range of magnetizations (from 0 to 1), pressure contrasts, and inertia ratios to characterize the formation and evolution of recollimation shocks. The results are further evaluated using linear stability analysis to assess the growth and suppression of CFI. Finally, we perform 3D simulations of unstable and stable jets. We discuss how the different parameters of the axisymmetric steady solutions influence the location and strength of recollimation. We find that, even in moderately magnetized jets, $\sigma$=0.1, the CFI can still develop under suitable local conditions and disrupt the jet structure. This instability is governed by the jet radius, curvature, Lorentz factor, and magnetization, and is not always predictable from injection conditions. While magnetization can delay or locally suppress instability growth, it does not guarantee long-term jet stability. Our 3D results highlight the limitations of 2D models in capturing non-axisymmetric and nonlinear effects, and underline the complex interplay between magnetic confinement and destabilizing mechanisms. These findings have implications for interpreting variability, and polarization structure in AGN jets.

We report the discovery and characterization of a wide binary population in the ultrafaint dwarf galaxy Boötes I using deep JWST/NIRCam imaging. Our sample consists of 52 candidate binaries with projected separations of 7,000 - 16,000 au and stellar masses from near the hydrogen-burning limit to the main-sequence turnoff ($\sim0.1$ - $0.8~{\rm M_\odot}$). By forward-modeling selection biases and chance alignments, we find that $1.25\pm0.25\%$ of Boötes I stars are members of wide binaries with separations beyond 5,000 au. This fraction, along with the distributions of separations and mass ratios, matches that in the Solar neighborhood, suggesting that wide binary formation is largely insensitive to metallicity, even down to [Fe/H] $\approx -2.5$. The observed truncation in the separation distribution near 16,000 au is well explained by stellar flyby disruptions. We also discuss how the binaries can be used to constrain the galaxy's dark matter properties. We show that our detection places new limits on primordial black hole dark matter, finding that compact objects with $M \gtrsim 5~{\rm M_\odot}$ cannot constitute more than $\sim1\%$ of the dark matter content. In contrast to previous work, we find that wide binaries are unlikely to provide robust constraints on the dark matter profile of ultrafaint galaxies given the uncertainties in the initial binary population, flyby disruptions, and contamination from chance alignments. These findings represent the most robust detection of wide binaries in an external galaxy to date, opening a new avenue for studying binary star formation and survival in extreme environments.

While the James Webb Space Telescope (JWST) now allows identifying quiescent galaxies (QGs) out to early epochs, the photometric selection of quiescent galaxy candidates (QGCs) and the derivation of key physical quantities are highly sensitive to the assumed star-formation histories (SFHs). We aim to quantify how the inclusion of JWST/MIRI data and different SFH models impacts the selection and characterisation of QGCs. We test the robustness of the physical properties inferred from the spectral energy distribution (SED) fitting, such as M*, age, star formation rate (SFR), and AV, and study how they impact the quiescence criteria of the galaxies across cosmic time. We perform SED fitting for ~13000 galaxies at z<6 from the CEERS/MIRI fields with up to 20 optical-mid infrared (MIR) broadband coverage. We implement three SFH prescriptions: flexible delayed, NonParametric, and extended Regulator. For each model, we compare results obtained with and without MIRI photometry and dust emission models. We evaluate the impact of these configurations on the number of candidate QGCs, selected based on rest UVJ colours, sSFR and main-sequence offset, and on their key physical properties such as M*, AV, and stellar ages. The number of QGCs selected varies significantly with the choice of SFH from 171 to 224 out of 13000 galaxies, depending on the model. This number increases to 222-327 when MIRI data are used (up to ~45% more QGCs). This enhancement is driven by improved constraints on dust attenuation and M*. We find a strong correlation between AV and M*, with massive galaxies (M*~10^11 M\odot) being 1.5-4.2 times more attenuated in magnitude than low-mass systems (M*~10^9 M\odot), depending on SFH. Regardless of the SFH assumption, ~13% of QGCs exhibit significant attenuation (AV > 0.5) in support of recent JWST studies challenging the notion that quiescent galaxies are uniformly dust-free.

We report the discovery of 164 compact (radius < 1 arcmin) radio rings using MeerKAT 1.3 GHz data from the SARAO MeerKAT Galactic Plane Survey (l=2-60deg, 252-358deg, |b|<1.5deg) and the Galactic Centre mosaic, from a search aimed at identifying previously uncatalogued radio sources. Within this sample, approximately 19 per cent of the rings contain a central point radio source. A multiwavelength analysis reveals a striking diversity: about 40 per cent of the rings enclose an isolated infrared point source, 50 per cent exhibit an extended counterpart in the mid- or far-infrared, and several are only detected in the radio band. We found that 17 per cent of the rings in the sample are positionally coincident (within 5 arcsec) with known entries in SIMBAD, including unclassified infrared sources, spiral galaxies, young stellar objects and long-period variable candidates. Based on these matches and exploiting ancillary multiwavelength data and catalogues, we explore several formation scenarios for the rings, such as HII regions, planetary nebulae, mass-loss relics from evolved massive stars, supernova remnants, nova shells, galaxies, galaxy cluster lenses and odd radio circles. Tentative classifications are proposed for nearly 60 per cent of the sample. These results highlight the potential of MeerKAT to uncover previously undetected compact radio structures and, particularly, recover missing Galactic radio-emitting objects.

It is often understood that supernova (SN) feedback in galaxies is responsible for regulating star formation and generating gaseous outflows. However, a detailed look at their effect on the local interstellar medium (ISM) on small mass scales in simulations shows that these processes proceed in clearly distinct channels. We demonstrate this finding in two independent simulations with solar-mass resolution, LYRA and RIGEL, of an isolated dwarf galaxy. Focusing on the immediate environment surrounding SNe, our findings suggest that the large-scale effect of a given SN on the galaxy is best predicted by its immediate local density. Outflows are driven by SNe in diffuse regions expanding to their cooling radii on large (\sim kpc) scales, while dense star-forming regions are disrupted in a localized (\sim pc) manner. However, these separate feedback channels are only distinguishable at very high numerical resolutions capable of following scales <<10^3 \msun. On larger scales, ISM densities are greatly mis-estimated, and differences between local environments of SNe become severely washed out. We demonstrate the practical implications of this effect by comparing with a mid-resolution simulation (Mptcl. \sim 200 \msun) of the same dwarf using the SMUGGLE model. The coarse-resolution simulation cannot self-consistently determine whether a given SN is responsible for generating outflows or suppressing star formation, suggesting that emergent galaxy physics such as star formation regulation through hot-phase outflows is fundamentally unresolvable by subgrid stellar feedback models, without appealing directly to simulations with highly resolved ISM.

We report the discovery of a dense molecular ring-like structure in a dense (10$^5$ cm$^{-3}$), cold (pc-scale CO depletion at a factor of 5), and young (10$^4$ year) star-forming region G34.74-0.12, revealed by C$^{18}$O (2-1), HNC (1-0), and N$_2$H$^+$ (1-0) observations with the Atacama Large Millimeter/submillimeter Array (ALMA). The ring-like structure is redshifted with respect to the clump, spanning from $V_{\rm sys,lsr} + 0.9$ to $V_{\rm sys,lsr} + 2.9$ km s$^{-1}$, with a total mass of 109 $M_{\odot}$. It is spatially coincident with 1.3 mm and 3.0 mm dust continuum emission from cores, and several protostellar outflows. However, no free-free emission or H\textsc{ii} region is detected in association with this structure. With a slow expansion speed indicated by the position-velocity diagram, this ring structure differs from rings previously identified in more evolved star-forming regions. Possible explanations for the ring-like structure include a relic wind-blown bubble produced by a deeply embedded young stellar object, a hollow cavity formed by cloud-cloud interactions, a gas ring resulting from a temperature gradient, or a line-of-sight superposition of multiple outflows or dense clouds. This discovery offers a rare observational glimpse into the earliest dynamical processes involved in massive star formation.

Red giants (RGs) provide a promising astrophysical environment for capturing dark matter (DM) via elastic scattering with stellar nuclei. Captured DM particles migrate toward the helium-rich core and accumulate into a compact configuration. As the DM population grows, it can become self-gravitating and undergo gravitational collapse, leading to adiabatic contraction through interactions with the ambient medium. The resulting energy release, through elastic scattering and, where relevant, DM annihilation during collapse, locally heats the stellar core and can trigger helium ignition earlier than that predicted by standard stellar evolution. We analyze the conditions under which DM-induced heating leads to runaway helium burning and identify the critical DM mass required for ignition. Imposing the observational constraint that helium ignition must not occur before the observed luminosity at the tip of the RG branch, we translate these conditions into bounds on DM properties. Remarkably, we find that RGs are sensitive to DM, particularly with masses around $10^{11} \,{\rm GeV}$ and spin-independent scattering cross sections near $10^{-37}\,{\rm cm}^2$, which is comparable to the reach of current terrestrial direct detection experiments.

The ESA Euclid mission will measure the photometric redshifts of billions of galaxies in order to provide an accurate 3D view of the Universe at optical and near-infrared wavelengths. Photometric redshifts are determined by the PHZ processing function on the basis of the multi-wavelength photometry of Euclid and ground-based observations. In this paper, we describe the PHZ processing used for the Euclid Quick Data Release, the output products, and their validation. The PHZ pipeline is responsible for the following main tasks: source classification into star, galaxy, and QSO classes based on photometric colours; determination of photometric redshifts and of physical properties of galaxies. The classification is able to provide a star sample with a high level of purity, a highly complete galaxy sample, and reliable probabilities of belonging to those classes. The identification of QSOs is more problematic: photometric information seems to be insufficient to accurately separate QSOs from galaxies. The performance of the pipeline in the determination of photometric redshifts has been tested using the COSMOS2020 catalogue and a large sample of spectroscopic redshifts. The results are in line with expectations: the precision of the estimates are compatible with Euclid requirements, while, as expected, a bias correction is needed to achieve the accuracy level required for the cosmological probes. Finally, the pipeline provides reliable estimates of the physical properties of galaxies, in good agreement with findings from the COSMOS2020 catalogue, except for an unrealistically large fraction of very young galaxies with very high specific star-formation rates. The application of appropriate priors is, however, sufficient to obtain reliable physical properties for those problematic objects. We present several areas for improvement for future Euclid data releases.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

The recently discovered population of faint FR0 radiogalaxies has been interpreted as the extension to low power of the classical FRI sources. Their radio emission appears to be concentrated in very compact (pc-scale) cores, any extended emission is very weak or absent and VLBI observations show that jets are already mildly or sub-relativistic at pc scales. Based on these observational properties we propose here that the jets of FR0s are strongly decelerated and disturbed at pc scale by hydrodynamical instabilities.}{With the above scenario in mind, we study the dynamics of a low-power relativistic jet propagating into a confining external medium, focusing on the effects of entrainment and mixing promoted by the instabilities developing at the jet-environment interface downstream of a recollimation shock. We perform a 3D relativistic hydrodynamical simulation of a recollimated jet by means of the state-of-the-art code PLUTO. The jet is initially conical, relativistic (with initial Lorentz Factor $\Gamma$=5), cold and light with respect to the confining medium, whose pressure is assumed to slowly decline with distance. The magnetic field is assumed to be dynamically unimportant. The 3D simulation shows that, after the first recollimation/reflection shock system, a rapidlygrowing instability develops, as a result of the interplay between recollimation-induced instabilities and Richtmyer-Meshkov modes. In turn, the instabilities promote strong mixing and entrainment that rapidly lead to the deceleration of the jet and spread its momentum to slowly moving, highly turbulent external gas. We argue that this mechanism could account for the peculiarities of the low-power FR0 jets. For outflows with higher power, Lorentz factor or magnetic field, we expect that the destabilizing effects are less effective, allowing the survival of the jet up to the kpc scale, as observed in FRIs.

This is the second paper in the HOWLS (higher-order weak lensing statistics) series exploring the usage of non-Gaussian statistics for cosmology inference within \textit{Euclid}. With respect to our first paper, we develop a full tomographic analysis based on realistic photometric redshifts which allows us to derive Fisher forecasts in the ($\sigma_8$, $w_0$) plane for a \textit{Euclid}-like data release 1 (DR1) setup. We find that the 5 higher-order statistics (HOSs) that satisfy the Gaussian likelihood assumption of the Fisher formalism (1-point probability distribution function, $\ell$1-norm, peak counts, Minkowski functionals, and Betti numbers) each outperform the shear 2-point correlation functions by a factor $2.5$ on the $w_0$ forecasts, with only marginal improvement when used in combination with 2-point estimators, suggesting that every HOS is able to retrieve both the non-Gaussian and Gaussian information of the matter density field. The similar performance of the different estimators\inlinecomment{, with a slight preference for Minkowski functionals and 1-point probability distribution function,} is explained by a homogeneous use of multi-scale and tomographic information, optimized to lower computational costs. These results hold for the $3$ mass mapping techniques of the \textit{Euclid} pipeline: aperture mass, Kaiser--Squires, and Kaiser--Squires plus, and are unaffected by the application of realistic star masks. Finally, we explore the use of HOSs with the Bernardeau--Nishimichi--Taruya (BNT) nulling scheme approach, finding promising results towards applying physical scale cuts to HOSs.

We present new ALMA [OIII]$_{88}$ observations of eight previously [CII]$_{158}$-detected galaxies at $6.8 \lesssim z \lesssim 7.7$. Six of our targets -- the primary sample -- are massive, UV-luminous galaxies drawn from the REBELS survey, while the remaining two are UV-fainter galaxies that were previously serendipitously detected through their luminous [CII] lines in the REBELS fields. We detect [OIII]$_{88}$ emission in all eight galaxies at $6.2 - 17.7\sigma$ significance, and find them to be consistent with the local dwarf galaxy relation between $L_\mathrm{[OIII]}$ and star formation rate. Our sample spans [OIII]/[CII] $\approx 1.9 - 9.6$, which is typical for the high-redshift galaxy population. Five of the primary targets benefit from JWST/NIRSpec observations, enabling a direct comparison of the [OIII]/[CII] ratio against rest-optical ISM diagnostics. We supplement our high-redshift sample with eleven $z\approx6-14$ galaxies in the literature for which similar ALMA and JWST observations are available, and furthermore compare to the [OIII]/[CII] ratios measured for local dwarf galaxies. We find that, at fixed metallicity and ionization parameter, $z>6$ galaxies show elevated [OIII]/[CII] ratios compared to local dwarfs. Instead, we find that a large [OIII]$_{4959,5007}$+H$\beta$ equivalent width -- a proxy for burstiness -- is the main driver of the high [OIII]/[CII] ratios seen in the early Universe, which is primarily due to [CII] being suppressed in bursty galaxies. Given the apparent validity of the [OIII]$_{88}$-SFR relation across most of cosmic time, as well as the abundance of young, bursty galaxies at high redshift, [OIII]$_{88}$ is set to remain a powerful ISM tracer at the cosmic dawn.

The existence of galaxies with no heavy elements is a key prediction of cosmological models. So far no "zero-metallicity", or Population III, galaxies have been identified. Here, we report the identification of an extremely metal-poor galaxy at redshift z = 5.725 ("AMORE6"), multiply imaged by a foreground galaxy cluster. JWST spectra consistently detect H$\beta$ at both positions, but [OIII]${\lambda}{\lambda}$4960,5008 remains undetected. This places a firm upper limit on its oxygen abundance, < 0.12 % of solar metallicity at 2${\sigma}$, establishing itself as the most pristine object by far. AMORE6 exhibits exceptional properties that indicate the presence of pristine massive star formation. The finding of such an example at a relatively late cosmic time is surprising, but it also validates the basic ideas behind the Big Bang model.

We use the first observation of Betelgeuse in hard $X$-rays to perform a novel search for axion-like particles (ALPs). Betelgeuse is not expected to be a standard source of $X$-rays, but light ALPs produced in the stellar core could be converted back into photons in the Galactic magnetic field, producing a detectable flux that peaks in the hard $X$-ray band (E_\gamma&gt;10\mathrm{\,keV}). Using a 50 ks observation of Betelgeuse by the $NuSTAR$ satellite telescope, we find no significant excess of events above the expected background. Using models of the regular Galactic magnetic field in the direction of Betelgeuse, we set a 95% C.L. upper limit on the ALP-photon coupling of {g_{a\gamma}&lt;(0.5-1.8)\times10^{-11}} GeV$^{-1}$ (depending on magnetic field model) for ALP masses {m_{a}&lt;(5.5-3.5) \times10^{-11}} eV.

The star-forming main sequence (SFMS) is a tight relation observed between stellar masses and star formation rates (SFR) in a population of galaxies. This relation is observed at different redshifts, in various morphological, and environmental domains, and is key to understanding the underlying relations between a galaxy budget of cold gas and its stellar content. Euclid Quick Data Release 1 (Q1) gives us the opportunity to investigate this fundamental relation in galaxy formation and evolution. We complement the Euclid release with public IRAC observations of the Euclid Deep Fields, improving the quality of recovered photometric redshifts, stellar masses, and SFRs, as is shown both with simulations and a comparison with available spectroscopic redshifts. From Q1 data alone, we recover more than $\sim 30\,\mathrm{k}$ galaxies with $\log_{10} (M_\ast/M_\odot) > 11$, giving a precise constraint of the SFMS at the high-mass end. We investigated the SFMS, in a redshift interval between $0.2$ and $3.0$, comparing our results with the existing literature and fitting them with a parameterisation taking into account the presence of a bending of the relation at the high-mass end, depending on the bending mass, $M_0$. We find good agreement with previous results in terms of $M_0$ values, and an increasing trend for the relation scatter at higher stellar masses. We also investigate the distribution of physical (e.g. dust absorption, $A_V$, and formation age) and morphological properties (e.g., Sérsic index and radius) in the SFR--stellar mass plane, and their relation with the SFMS. These results highlight the potential of Euclid in studying the fundamental scaling relations that regulate galaxy formation and evolution in anticipation of the forthcoming Data Release 1.

Turbulence in the interstellar medium (ISM) plays an important role in many physical processes, including forming stars and shaping complex ISM structures. In this work, we investigate the HI turbulent properties of the Small Magellanic Cloud (SMC) to reveal what physical mechanisms drive the turbulence and at what scales. Using the high-resolution HI data of the Galactic ASKAP (GASKAP) survey and multi-point structure functions (SF), we perform a statistical analysis of HI turbulence in 34 subregions of the SMC. Two-point SFs tend to show a linear trend, and their slope values are relatively uniform across the SMC, suggesting that large-scale structures exist and are dominant in the two-point SFs. On the other hand, seven-point SF enables us to probe small-scale turbulence by removing large-scale fluctuations, which is difficult to achieve with the two-point SFs. In the seven-point SFs, we find break features at scales of 34-84 pc, with a median scale of $\sim$50 pc. This result indicates the presence of small-scale turbulent fluctuations in the SMC and quantifies its scale. In addition, we find strong correlations between slope values of the seven-point SFs and the stellar feedback-related quantities (e.g., H$\alpha$ intensities, the number of young stellar objects, and the number of HI shells), suggesting that stellar feedback may affect the small-scale turbulent properties of the HI gas in the SMC. Lastly, estimated sonic Mach numbers across the SMC are subsonic, which is consistent with the fact that the HI gas of the SMC primarily consists of the warm neutral medium.

We assess how physically realistic the ''simulation hypothesis'' for this Universe is, based on physical constraints arising from the link between information and energy, and on known astrophysical constraints. We investigate three cases: the simulation of the entire visible Universe, the simulation of Earth only, or a low resolution simulation of Earth, compatible with high-energy neutrino observations. In all cases, the amounts of energy or power required by any version of the simulation hypothesis are entirely incompatible with physics, or (literally) astronomically large, even in the lowest resolution case. Only universes with very different physical properties can produce some version of this Universe as a simulation. On the other hand, our results show that it is just impossible that this Universe is simulated by a universe sharing the same properties, regardless of technological advancements of the far future.

One of the most puzzling properties of the high-redshift AGN population recently discovered by JWST, including both broad-line and narrow-line sources, is their X-ray weakness. With very few exceptions, and regardless of the optical classification, they are undetected at the limits of the deepest Chandra fields, even when stacking signals from tens of sources in standard observed-frame energy intervals (soft, hard, and full bands). It has been proposed that their elusive nature in the X-ray band is due to heavy absorption by dust-free gas or intrinsic weakness, possibly due to high, super-Eddington accretion. In this work, we perform X-ray stacking in three customized rest-frame energy ranges (1-4, 4-7.25, and 10-30 keV) of a sample of 50 Type 1 and 38 Type 2 AGN identified by JWST in the CDFS and CDFN fields. For the Type 2 sub-sample, we reach a total of about 210 Ms exposure, and we report a significant ($\sim 3\sigma$) detection in the hardest (10-30 keV rest frame) band, along with relatively tight upper limits in the rest frame softer energy bands. The most straightforward interpretation is in terms of heavy obscuration due to gas column densities well within the Compton thick regime ($> 2 \times 10^{24}$cm$^{-2}$) with a large covering factor, approaching 4$\pi$. The same procedure applied to the Type 1 sub-sample returns no evidence for a significant signal in about 140 Ms stacked data in any of the adopted bands, confirming their surprisingly elusive nature in the X-ray band obtained with previous stacking experiments. A brief comparison with the current observations and the implications for the evolution of AGN are discussed.