We provide an accurate comparison, against large cosmological $N$-body simulations, of different prescriptions for modelling nonlinear matter power spectra in the presence of massive neutrinos and dynamical dark energy. We test the current most widely used approaches: fitting functions (HALOFIT and HMcode), the halo-model reaction (ReACT) and emulators (baccoemu and EuclidEmulator2). Focussing on redshifts $z\leq2$ and scales $k\lesssim 1 \ h/$Mpc (where the simulation mass resolution provides $\sim 1\%$ accuracy), we find that HMcode and ReACT considerably improve over the HALOFIT prescriptions of Smith and Takahashi (both combined with the Bird correction), with an overall agreement of 2\% for all the cosmological scenarios considered. Concerning emulators, we find that, especially at low redshifts, EuclidEmulator2 remarkably agrees with the simulated spectra at $\lesssim 1\%$ level in scenarios with dynamical dark energy and massless neutrinos, reaching a maximum difference of $\sim 2\%$ at $z=2$. baccoemu has a similar behaviour as EuclidEmulator2, except for a couple of dark energy models. In cosmologies with massive neutrinos, at $z=0$ all the nonlinear prescriptions improve their agreement with respect to the massless neutrino case, except for the Bird and TakaBird models which, however, are not tailored to $w_0$--$w_a$ models. At $z>0$ we do not find a similar improvement when including massive neutrinos, probably due to the lower impact of neutrino free-streaming at higher redshifts; rather at $z=2$ EuclidEmulator2 exceeds $2\%$ agreement for some dark energy equation of state. When considering ratios between the matter power spectrum computed in a given cosmological model and its $\Lambda$CDM counterpart, all the tested prescriptions agree with simulated data, at sub-percent or percent level, depending on $z$. [ABRIDGED]

We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from the combination of weak gravitational lensing and galaxy clustering data. The breath of Euclid's data will also foster a wide variety of scientific analyses. The Flagship simulation was developed to provide a realistic approximation to the galaxies that will be observed by Euclid and used in its scientific analyses. We ran a state-of-the-art N-body simulation with four trillion particles, producing a lightcone on the fly. From the dark matter particles, we produced a catalogue of 16 billion haloes in one octant of the sky in the lightcone up to redshift z=3. We then populated these haloes with mock galaxies using a halo occupation distribution and abundance matching approach, calibrating the free parameters of the galaxy mock against observed correlations and other basic galaxy properties. Modelled galaxy properties include luminosity and flux in several bands, redshifts, positions and velocities, spectral energy distributions, shapes and sizes, stellar masses, star formation rates, metallicities, emission line fluxes, and lensing properties. We selected a final sample of 3.4 billion galaxies with a magnitude cut of H_E<26, where we are complete. We have performed a comprehensive set of validation tests to check the similarity to observational data and theoretical models. In particular, our catalogue is able to closely reproduce the main characteristics of the weak lensing and galaxy clustering samples to be used in the mission's main cosmological analysis. (abridged)

Euclid will cover over 14000 $deg^{2}$ with two optical and near-infrared spectro-photometric instruments, and is expected to detect around ten million active galactic nuclei (AGN). This unique data set will make a considerable impact on our understanding of galaxy evolution and AGN. In this work we identify the best colour selection criteria for AGN, based only on Euclid photometry or including ancillary photometric observations, such as the data that will be available with the Rubin legacy survey of space and time (LSST) and observations already available from Spitzer/IRAC. The analysis is performed for unobscured AGN, obscured AGN, and composite (AGN and star-forming) objects. We make use of the spectro-photometric realisations of infrared-selected targets at all-z (SPRITZ) to create mock catalogues mimicking both the Euclid Wide Survey (EWS) and the Euclid Deep Survey (EDS). Using these catalogues we estimate the best colour selection, maximising the harmonic mean (F1) of completeness and purity. The selection of unobscured AGN in both Euclid surveys is possible with Euclid photometry alone with F1=0.22-0.23, which can increase to F1=0.43-0.38 if we limit at z>0.7. Such selection is improved once the Rubin/LSST filters (a combination of the u, g, r, or z filters) are considered, reaching F1=0.84 and 0.86 for the EDS and EWS, respectively. The combination of a Euclid colour with the [3.6]-[4.5] colour, which is possible only in the EDS, results in an F1-score of 0.59, improving the results using only Euclid filters, but worse than the selection combining Euclid and LSST. The selection of composite ($f_{\rm AGN}$=0.05-0.65 at 8-40 $\mu m$) and obscured AGN is challenging, with F1<0.3 even when including ancillary data. This is driven by the similarities between the broad-band spectral energy distribution of these AGN and star-forming galaxies in the wavelength range 0.3-5 $\mu m$.

In the past years, the results obtained by the WISSH quasar project provided a novel general picture on the distinctive multi-band properties of hyper-luminous ($L_{bol}>10^{47}$ erg/s) quasars at high redshift (z$\sim$2-4), unveiling interesting relations among active galactic nuclei, winds and interstellar medium, in these powerful sources at cosmic noon. Since 2022, we are performing a systematic and statistically-significant VLA study of the radio properties of WISSH. We carried out high-resolution VLA observations aiming at: 1) identifying young radio source from the broad-band spectral shape of these objects; 2) sample an unexplored high redshift/high luminosity regime, tracking possible evolutionary effects on the radio-loud/radio-quiet dichotomy; 3) quantifying orientation effects on the observed winds/outflows properties.

The Euclid Space Telescope will provide deep imaging at optical and near-infrared wavelengths, along with slitless near-infrared spectroscopy, across ~15,000 sq deg of the sky. Euclid is expected to detect ~12 billion astronomical sources, facilitating new insights into cosmology, galaxy evolution, and various other topics. To optimally exploit the expected very large data set, there is the need to develop appropriate methods and software. Here we present a novel machine-learning based methodology for selection of quiescent galaxies using broad-band Euclid I_E, Y_E, J_E, H_E photometry, in combination with multiwavelength photometry from other surveys. The ARIADNE pipeline uses meta-learning to fuse decision-tree ensembles, nearest-neighbours, and deep-learning methods into a single classifier that yields significantly higher accuracy than any of the individual learning methods separately. The pipeline has `sparsity-awareness', so that missing photometry values are still informative for the classification. Our pipeline derives photometric redshifts for galaxies selected as quiescent, aided by the `pseudo-labelling' semi-supervised method. After application of the outlier filter, our pipeline achieves a normalized mean absolute deviation of ~< 0.03 and a fraction of catastrophic outliers of ~< 0.02 when measured against the COSMOS2015 photometric redshifts. We apply our classification pipeline to mock galaxy photometry catalogues corresponding to three main scenarios: (i) Euclid Deep Survey with ancillary ugriz, WISE, and radio data; (ii) Euclid Wide Survey with ancillary ugriz, WISE, and radio data; (iii) Euclid Wide Survey only. Our classification pipeline outperforms UVJ selection, in addition to the Euclid I_E-Y_E, J_E-H_E and u-I_E,I_E-J_E colour-colour methods, with improvements in completeness and the F1-score of up to a factor of 2. (Abridged)

Lorentz invariance violation (LIV) has long been recognized as an observable low-energy signature of quantum gravity. In spite of a great effort to detect LIV effects, so far only lower bounds have been derived. The high energy photons from the gamma ray burst GRB 221009A have been detected by the LHAASO collaboration and one at ${\cal E} \simeq 251 \, \rm TeV$ by the Carpet collaboration using a partial data set. Very recently, the Carpet collaboration has completed the full data analysis, reporting further support for their previously detected photon now at ${\cal E} = 300^{+ 43}_{- 38} \, {\rm TeV}$, which manifestly clashes with conventional physics. Taking this result at face value, we derive the first evidence for LIV and we show that such a detection cannot be explained by axion-like particles (ALPs), which allow for the observation of the highest energy photons detected by LHAASO. We also outline a scenario in which ALPs and LIV naturally coexist. If confirmed by future observations our finding would represent the first positive result in quantum gravity phenomenology.

We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$\lambda$3727 $Å$ line luminosity. We find that 49% of the AGN host galaxies are on or above the main sequence (MS), 40% lie in the sub-MS locus, and 11% in the quiescent locus. Using the [OIII]$\lambda$5007 $Å$ line luminosity as a proxy of the AGN power, we find that at fixed AGN power Type II AGN host galaxies show a bimodal behaviour: systems with host galaxy stellar mass <10$^{10}$ M$_{\odot}$, reside along the MS or in the starbursts locus (high-SF Type II AGN), while systems residing in massive host-galaxies (>10$^{10}$ M$_{\odot}$) show a lower level of star formation (low-SF Type II AGN). At all stellar masses, the offset from the MS is positively correlated with the AGN power. We interpret this correlation as evidence of co-evolution between the AGN and the host, possibly due to the availability of cold gas. In the most powerful AGN with host galaxies below the MS we find a hint, though weak, of asymmetry in the [OIII] line profile, likely due to outflowing gas, consistent with a scenario in which AGN feedback removes the available gas and halts the star formation in the most massive hosts.

The massive, hot galaxy cluster PSZ2 G286.98+32.90 (hereafter PLCKG287, z=0.383) hosts a giant radio halo and two prominent radio relics which are signs of a disturbed dynamical state. However, despite optical and radio observations indicate a clear multiple merger, the X-ray emission of the cluster, derived from XMM-Newton observations, shows only moderate disturbance. We present new 200 ks Chandra observations of PLCKG287. We detect a shock front to the NW direction at a distance of ~390 kpc from the X-ray peak, characterized by a Mach number M~1.3, as well as a cold front at a distance of ~300 kpc from the X-ray peak, nested in the same direction of the shock in a typical configuration expected by a merger. We also find evidence for X-ray depressions to the E and W, that could be the signature of feedback from the active galactic nucleus (AGN). The radial profile of the thermodynamic quantities show a temperature and abundance peak in the cluster center, where also the pressure and entropy have a rapid increase. Based on these properties, we argue that PLCKG287 is what remains of a cool core after a heating event. We estimate that both the shock energy and the AGN feedback energy, implied by the analysis of the X-ray cavities, are sufficient to heat the core to the observed temperature of ~17 keV in the central ~160 kpc. We discuss the possible origin of the detected shock by investigating alternative scenarios of merger and AGN outburst, finding that they are both energetically viable. However, no single model seems able to explain all the X-ray features detected in this system. This suggests that the combined action of merger and central AGN feedback is likely necessary to explain the reheated cool core, the large-scale shock and the cold front. The synergy of these two processes may act in shaping the distribution of cool core and non cool core clusters. [Abridged]

We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$\lambda$3727 $Å$ line luminosity. We find that 49% of the AGN host galaxies are on or above the main sequence (MS), 40% lie in the sub-MS locus, and 11% in the quiescent locus. Using the [OIII]$\lambda$5007 $Å$ line luminosity as a proxy of the AGN power, we find that at fixed AGN power Type II AGN host galaxies show a bimodal behaviour: systems with host galaxy stellar mass <10$^{10}$ M$_{\odot}$, reside along the MS or in the starbursts locus (high-SF Type II AGN), while systems residing in massive host-galaxies (>10$^{10}$ M$_{\odot}$) show a lower level of star formation (low-SF Type II AGN). At all stellar masses, the offset from the MS is positively correlated with the AGN power. We interpret this correlation as evidence of co-evolution between the AGN and the host, possibly due to the availability of cold gas. In the most powerful AGN with host galaxies below the MS we find a hint, though weak, of asymmetry in the [OIII] line profile, likely due to outflowing gas, consistent with a scenario in which AGN feedback removes the available gas and halts the star formation in the most massive hosts.

In this review, we present a self-contained introduction to axion-like particles (ALPs) with a particular focus on their effects on photon polarization: both theoretical and phenomenological aspects are discussed. We derive the photon survival probability in the presence of photon--ALP interaction, the corresponding final photon degree of linear polarization, and the polarization angle in a wide energy interval. The presented results can be tested by current and planned missions such as IXPE (already operative), eXTP, XL-Calibur, NGXP, XPP in the X-ray band and like COSI (approved to launch), e-ASTROGAM, and AMEGO in the high-energy range. Specifically, we describe ALP-induced polarization effects on several astrophysical sources, such as galaxy clusters, blazars, and gamma-ray bursts, and we discuss their real detectability. In particular, galaxy clusters appear as very good observational targets in this respect. Moreover, in the very-high-energy (VHE) band, we discuss a peculiar ALP signature in photon polarization, in principle capable of proving the ALP existence. Unfortunately, present technologies cannot detect photon polarization up to such high energies, but the observational capability of the latter ALP signature in the VHE band could represent an interesting challenge for the future. As a matter of fact, the aim of this review is to show new ways to make progress in the physics of ALPs, thanks to their effects on photon polarization, a topic that has aroused less interest in the past, but which is now timely with the advent of many new polarimetric missions.

The International Axion Observatory (IAXO) is a next-generation axion helioscope designed to search for solar axions with unprecedented sensitivity. IAXO holds a unique position in the global landscape of axion searches, as it will probe a region of the axion parameter space inaccessible to any other experiment. In particular, it will explore QCD axion models in the mass range from meV to eV, covering scenarios motivated by astrophysical observations and potentially extending to axion dark matter models. Several studies in recent years have demonstrated that IAXO has the potential to probe a wide range of new physics beyond solar axions, including dark photons, chameleons, gravitational waves, and axions from nearby supernovae. IAXO will build upon the two-decade experience gained with CAST, the detailed studies for BabyIAXO, which is currently under construction, as well as new technologies. If, in contrast to expectations, solar axion searches with IAXO ``only'' result in limits on new physics in presently uncharted parameter territory, these exclusions would be very robust and provide significant constraints on models, as they would not depend on untestable cosmological assumptions.

We report the results on the short gamma-ray burst GRB 241107A, obtained with the IBIS instrument on board the INTEGRAL satellite. The burst had a duration of about 0.2 s, a fluence of $8 \times 10^{-7}$ erg cm-2 in the 20 keV-10 MeV range and a hard spectrum, characterized by a peak energy of 680 keV. The position of GRB 241107A has been precisely determined because it fell inside the imaging field of view of the IBIS coded mask instrument. The presence of the nearby galaxy PGC 86046 in the 3 arcmin radius error region, suggests that GRB 241107A might be a giant flare from a magnetar rather than a canonical short GRB. For the 4.1 Mpc distance of PGC 86046, the isotropic energy of $1.6 \times 10^{45}$ erg is in agreement with this hypothesis, that is also supported by the time resolved spectral properties similar to those of the few other extragalactic magnetars giant flares detected so far.

Investigating the drivers of the quenching of star formation in galaxies is key to understanding their evolution. The Euclid mission will provide rich spatial and spectral data from optical to infrared wavelengths for millions of galaxies, enabling precise measurements of their star formation histories. Using the first Euclid Quick Data Release (Q1), we developed a probabilistic classification framework, that combines the average specific star-formation rate ($\rm sSFR_\tau$) inferred over two timescales ($\tau={10^8,10^9}$ yr), to categorize galaxies as `Ageing' (secularly evolving), `Quenched' (recently halted star formation), or `Retired' (dominated by old stars). We validated this methodology using synthetic observations from the IllustrisTNG simulation. Two classification methods were employed: a probabilistic approach, integrating posterior distributions, and a model-driven method optimizing sample purity and completeness using IllustrisTNG. At z&lt;0.1 and $M_\ast \gtrsim 3\times10^{8}\, M_\odot$, we obtain Euclid class fractions of 68-72%, 8-17%, and 14-19% for Ageing, Quenched, and Retired populations, respectively, consistent with previous studies. The evolution with redshift shows increasing/decreasing fraction of Ageing/Retired galaxies. The fraction of quenched systems shows a weaker dependence on stellar mass and redshift, varying between 5% and 15%. We analysed the mass-size-metallicity relation for each population. Ageing galaxies generally exhibit disc morphologies and low metallicities. Retired galaxies show compact structures and enhanced chemical enrichment, while Quenched galaxies form an intermediate population, more compact and chemically evolved than Ageing systems. This work demonstrates Euclid's great potential for elucidating the physical nature of the quenching mechanisms that govern galaxy evolution.

The brightest ever observed gamma ray burst GRB 221009A at redshift $z = 0.151$ was detected on October 9, 2022. Its highest energy photons have been recorded by the LHAASO collaboration up to above $12 \, \rm TeV$, and one of the at ${\cal E} = 251 \, \rm TeV$ by the Carpet-2 collaboration. Very recently, the Carpet-3 collaboration has completed the data analysis, showing that the evidence of the $251 \, {\rm TeV}$ photon is quite robust. Still, according to conventional physics photons with ${\cal E} \gtrsim 10 \, \rm TeV$ cannot be observed owing to the absorption by the extragalactic background light (EBL). Previously it has been demonstrated that an axion-like particle (ALP) with allowed parameters ensures the observability of the LHAASO photons. Here we show that the Lorentz invariance violation allows the ${\cal E} = 251 \, {\rm TeV}$ (now around 300 TeV) Carpet photon to be detected.

Heavy chemical elements such as iron in the intra-cluster medium (ICM) of galaxy clusters are a signpost of the interaction between the gas and stellar components. Observations of the ICM metallicity in present-day massive systems, however, pose a challenge to the underlying assumption that the cluster galaxies have produced the amount of iron that enriches the ICM. We evaluate the iron share between ICM and stars within simulated galaxy clusters with the twofold aim of investigating the origin of possible differences with respect to observational findings and of shedding light on the observed excess of iron on the ICM with respect to expectations based on the observed stellar population. We evaluated the iron mass in gas and stars in a sample of 448 simulated systems with masses M500 > 1e14 Msun at z=0.07. These were extracted from the high-resolution (352 cMpc/h)^3 volume of the Magneticum cosmological hydrodynamical simulations. We compared our results with observational data of low-redshift galaxy clusters. The iron share in simulated clusters features a shallow dependence on the total mass, and its value is close to unity on average. In the most massive simulated systems, the iron share is thus smaller than observational values by almost an order of magnitude. The dominant contribution to this difference is related to the stellar component, whereas the chemical properties of the ICM agree well overall with the observations. We find larger stellar mass fractions in simulated massive clusters, which in turn yield higher stellar iron masses, than in observational data. Consistently with the modelling, we confirm that the stellar content within simulated present-day massive systems causes the metal enrichment in the ICM. It will be crucial to alleviate the stellar mass discrepancy between simulations and observations to definitely assess the iron budget in galaxy clusters.

We present a multi-band study of the diffuse emission in the galaxy cluster Abell 3558, located in the core of the Shapley Supercluster. Using new MeerKAT UHF-Band and uGMRT Band-3 observations, and published MeerKAT L-band and ASKAP 887 MHz data, we perform a detailed analysis of the diffuse emission in the cluster centre. We complement with XMM-Newton X-ray information for a thorough study of the connection between the thermal and non-thermal properties of the cluster. We find that the diffuse radio emission in the cluster centre is more extended than published earlier, with a previously undetected extension spanning 100 kpc towards the north beyond the innermost cold front, increasing the total size of the emission to 550 kpc, and shows a clear spatial correlation with the X-ray features. The overall radio spectrum is steep ($\alpha_{\rm 400\,MHz}^{\rm 1569\,MHz}=1.18\pm0.10$), with local fluctuations which show several connections with the X-ray surface brightness, cold fronts, and residual emission. The point-to-point correlation between the radio and X-ray surface brightness is sub-linear, and steepens with increasing frequency. We discuss the classification of the diffuse emission considering its overall properties, those of the ICM, and the existing scaling laws between the radio and X-ray quantities in galaxy clusters. We conclude that it is a mini-halo, powered by turbulent (re)-acceleration induced by sloshing motions within the cluster region delimited by the cold fronts, and it supports the picture of a known minor merger between A\,3558 and the group SC\,1327--312 with mass ratio 5:1.