Search

Understanding telescope pointing (i.e., line of sight) is important for observing the cosmic microwave background (CMB) and astronomical objects. The Moon is a candidate astronomical source for pointing calibration. Although the visible size of the Moon ($\ang{;30}$) is larger than that of the planets, we can frequently observe the Moon once a month with a high signal-to-noise ratio. We developed a method for performing pointing calibration using observational data from the Moon. We considered the tilts of the telescope axes as well as the encoder and collimation offsets for pointing calibration. In addition, we evaluated the effects of the nonuniformity of the brightness temperature of the Moon, which is a dominant systematic error. As a result, we successfully achieved a pointing accuracy of $\ang{;3.3}$. This is one order of magnitude smaller than an angular resolution of $\ang{;36}$. This level of accuracy competes with past achievements in other ground-based CMB experiments using observational data from the planets., Comment: 18 pages, 17 figures, 3 tables

We use the new QUIJOTE-MFI wide survey (11, 13, 17 and 19 GHz) to produce spectral energy distributions (SEDs), on an angular scale of 1 deg, of the supernova remnants (SNRs) CTB 80, Cygnus Loop, HB 21, CTA 1, Tycho and HB 9. We provide new measurements of the polarized synchrotron radiation in the microwave range. For each SNR, the intensity and polarization SEDs are obtained and modelled by combining QUIJOTE-MFI maps with ancillary data. In intensity, we confirm the curved power law spectra of CTB 80 and HB 21 with a break frequency $\nu_{\rm b}$ at 2.0$^{+1.2}_{-0.5}$ GHz and 5.0$^{+1.2}_{-1.0}$ GHz respectively; and spectral indices respectively below and above the spectral break of $-0.34\pm0.04$ and $-0.86\pm0.5$ for CTB 80, and $-0.24\pm0.07$ and $-0.60\pm0.05$ for HB 21. In addition, we provide upper limits on the Anomalous Microwave Emission (AME), suggesting that the AME contribution is negligible towards these remnants. From a simultaneous intensity and polarization fit, we recover synchrotron spectral indices as flat as $-0.24$, and the whole sample has a mean and scatter of $-0.44\pm0.12$. The polarization fractions have a mean and scatter of $6.1\pm1.9$\%. When combining our results with the measurements from other QUIJOTE studies of SNRs, we find that radio spectral indices are flatter for mature SNRs, and particularly flatter for CTB 80 ($-0.24^{+0.07}_{-0.06}$) and HB 21 ($-0.34^{+0.04}_{-0.03}$). In addition, the evolution of the spectral indices against the SNRs age is modelled with a power-law function, providing an exponent $-0.07\pm0.03$ and amplitude $-0.49\pm0.02$ (normalised at 10 kyr), which are conservative with respect to previous studies of our Galaxy and the Large Magellanic Cloud., Comment: 33 pages, 15 figure, 15 tables. Submitted to MNRAS. QUIJOTE data maps available at https://research.iac.es/proyecto/quijote

The Andromeda Galaxy (M31) is the Local Group galaxy that is most similar to the Milky Way (MW). The similarities between the two galaxies make M31 useful for studying integrated properties common to spiral galaxies. We use the data from the recent QUIJOTE-MFI Wide Survey, together with new raster observations focused on M31, to study its integrated emission. The addition of raster data improves the sensitivity of QUIJOTE-MFI maps by almost a factor 3. Our main interest is to confirm if anomalous microwave emission (AME) is present in M31, as previous studies have suggested. To do so, we built the integrated spectral energy distribution of M31 between 0.408 and 3000 GHz. We then performed a component separation analysis taking into account synchrotron, free-free, AME and thermal dust components. AME in M31 is modelled as a log-normal distribution with maximum amplitude, $A_{\rm AME}$, equal to $1.03\pm0.32$ Jy. It peaks at $\nu_{\rm AME}=17.2\pm3.2$ GHz with a width of $W_{\rm AME}=0.58\pm0.16$. Both the Akaike and Bayesian Information Criteria find the model without AME to be less than 1 % as probable as the one taking AME into consideration. We find that the AME emissivity per 100 $\mu$m intensity in M31 is $\epsilon_{\rm AME}^{\rm 28.4\,GHz}=9.6\pm3.1$ $\mu$K/(MJy/sr), similar to that computed for the MW. We also provide the first upper limits for the AME polarization fraction in an extragalactic object. M31 remains the only galaxy where an AME measurement has been made of its integrated spectrum., Comment: 17 pages, 9 figures. Accepted for publication in MNRAS. QUIJOTE data maps available at https://research.iac.es/proyecto/quijote

Anomalous Microwave Emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquituous in our Galaxy and is observed at a broad range of angular scales. Here we use the new QUIJOTE-MFI wide survey data at 11, 13, 17 and 19 GHz to constrain the AME in the Galactic plane ($|b|<10^\circ$) on degree scales. We built the spectral energy distribution between 0.408 and 3000 GHz for each of the 5309 0.9$^\circ$ pixels in the Galactic plane, and fitted a parametric model by considering five emission components: synchrotron, free-free, AME, thermal dust and CMB anisotropies. We show that not including QUIJOTE-MFI data points leads to the underestimation (up to 50 %) of the AME signal in favour of free-free emission. The parameters describing these components are then intercompared, looking for relations that help to understand AME physical processes. We find median values for the AME width, $W_{\rm AME}$, and for its peak frequency, $\nu_{\rm AME}$, respectively of $0.560^{+0.059}_{-0.050}$ and $20.7^{+2.0}_{-1.9}$ GHz, slightly in tension with current theoretical models. We find spatial variations throughout the Galactic plane for $\nu_{\rm AME}$, but only with reduced statistical significance. We report correlations of AME parameters with certain ISM properties, such as that between the AME emissivity (which shows variations with the Galactic longitude) and the interstellar radiation field, and that between the AME peak frequency and dust temperature. Finally, we discuss the implications of our results on the possible molecules responsible for AME., Comment: 32 pages, 31 figures. Accepted for publication in MNRAS. Derived data products available at https://research.iac.es/proyecto/quijote

The Andromeda galaxy (M31) is our closest neighbouring spiral galaxy, making it an ideal target for studying the physics of the interstellar medium in a galaxy very similar to our own. Using new observations of M31 at 4.76GHz by the C-Band All-Sky Survey (C-BASS), and all available radio data at $1^\circ$ resolution, we produce the integrated spectrum and put new constraints on the synchrotron spectral index and anomalous microwave emission (AME) from M31. We use aperture photometry and spectral modelling to fit for the integrated spectrum of M31, and subtract a comprehensive model of nearby background radio sources. The AME in M31 is detected at $3\sigma$ significance with a peak near 30GHz and flux density $0.27\pm0.09$Jy. The synchrotron spectral index of M31 is flatter than our own Galaxy at $\alpha = -0.66 \pm 0.03$ with no strong evidence of spectral curvature. The emissivity of AME, averaged over the total emission from M31 is lower than typical AME sources in our Galaxy, implying that AME is not uniformly distributed throughout M31 and instead is likely confined to sub-regions -- this will need to be confirmed using future higher resolution observations around 20--30GHz., Comment: 16 pages, 6 figures. Accepted by MNRAS

The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets., Comment: Under review by MNRAS; 28 pages, 16 figures

We present jax-cosmo, a library for automatically differentiable cosmological theory calculations. It uses the JAX library, which has created a new coding ecosystem, especially in probabilistic programming. As well as batch acceleration, just-in-time compilation, and automatic optimization of code for different hardware modalities (CPU, GPU, TPU), JAX exposes an automatic differentiation (autodiff) mechanism. Thanks to autodiff, jax-cosmo gives access to the derivatives of cosmological likelihoods with respect to any of their parameters, and thus enables a range of powerful Bayesian inference algorithms, otherwise impractical in cosmology, such as Hamiltonian Monte Carlo and Variational Inference. In its initial release, jax-cosmo implements background evolution, linear and non-linear power spectra (using halofit or the Eisenstein and Hu transfer function), as well as angular power spectra with the Limber approximation for galaxy and weak lensing probes, all differentiable with respect to the cosmological parameters and their other inputs. We illustrate how autodiff can be a game-changer for common tasks involving Fisher matrix computations, or full posterior inference with gradient-based techniques. In particular, we show how Fisher matrices are now fast, exact, no longer require any fine tuning, and are themselves differentiable. Finally, using a Dark Energy Survey Year 1 3x2pt analysis as a benchmark, we demonstrate how jax-cosmo can be combined with Probabilistic Programming Languages to perform posterior inference with state-of-the-art algorithms including a No U-Turn Sampler, Automatic Differentiation Variational Inference,and Neural Transport HMC. We further demonstrate that Normalizing Flows using Neural Transport are a promising methodology for model validation in the early stages of analysis.

This work considers which higher-order effects in modelling the cosmic shear angular power spectra must be taken into account for Euclid. We identify which terms are of concern, and quantify their individual and cumulative impact on cosmological parameter inference from Euclid. We compute the values of these higher-order effects using analytic expressions, and calculate the impact on cosmological parameter estimation using the Fisher matrix formalism. We review 24 effects and find the following potentially need to be accounted for: the reduced shear approximation, magnification bias, source-lens clustering, source obscuration, local Universe effects, and the flat Universe assumption. Upon computing these explicitly, and calculating their cosmological parameter biases, using a maximum multipole of $\ell=5000$, we find that the magnification bias, source-lens clustering, source obscuration, and local Universe terms individually produce significant ($\,>0.25\sigma$) cosmological biases in one or more parameters, and accordingly must be accounted for. In total, over all effects, we find biases in $\Omega_{\rm m}$, $\Omega_{\rm b}$, $h$, and $\sigma_{8}$ of $0.73\sigma$, $0.28\sigma$, $0.25\sigma$, and $-0.79\sigma$, respectively, for flat $\Lambda$CDM. For the $w_0w_a$CDM case, we find biases in $\Omega_{\rm m}$, $\Omega_{\rm b}$, $h$, $n_{\rm s}$, $\sigma_{8}$, and $w_a$ of $1.49\sigma$, $0.35\sigma$, $-1.36\sigma$, $1.31\sigma$, $-0.84\sigma$, and $-0.35\sigma$, respectively; which are increased relative to the $\Lambda$CDM due to additional degeneracies as a function of redshift and scale., Comment: 20 pages, submitted to A&A

The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consistent with the implementation within the Euclid Consortium of the dedicated processing function and find that, when jointly modelling mass and the concentration parameter of the Navarro-Frenk-White halo profile, the weak lensing masses tend to be, on average, biased low by 5-10% with respect to the true mass, up to z=0.5. Using a fixed value for the concentration $c_{200} = 3$, the mass bias is diminished below 5%, up to z=0.7, along with its relative uncertainty. Simulating the weak lensing signal by projecting along the directions of the axes of the moment of inertia tensor ellipsoid, we find that orientation matters: when clusters are oriented along the major axis, the lensing signal is boosted, and the recovered weak lensing mass is correspondingly overestimated. Typically, the weak lensing mass bias of individual clusters is modulated by the weak lensing signal-to-noise ratio, related to the redshift evolution of the number of galaxies used for weak lensing measurements: the negative mass bias tends to be larger toward higher redshifts. However, when we use a fixed value of the concentration parameter, the redshift evolution trend is reduced. These results provide a solid basis for the weak-lensing mass calibration required by the cosmological application of future cluster surveys from Euclid and Rubin., Comment: Accepted for publication in Astronomy & Astrophysics

Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten different HOS on a common set of $Euclid$-like mocks, derived from N-body simulations. In this first paper of the HOWLS series, we computed the nontomographic ($\Omega_{\rm m}$, $\sigma_8$) Fisher information for the one-point probability distribution function, peak counts, Minkowski functionals, Betti numbers, persistent homology Betti numbers and heatmap, and scattering transform coefficients, and we compare them to the shear and convergence two-point correlation functions in the absence of any systematic bias. We also include forecasts for three implementations of higher-order moments, but these cannot be robustly interpreted as the Gaussian likelihood assumption breaks down for these statistics. Taken individually, we find that each HOS outperforms the two-point statistics by a factor of around two in the precision of the forecasts with some variations across statistics and cosmological parameters. When combining all the HOS, this increases to a $4.5$ times improvement, highlighting the immense potential of HOS for cosmic shear cosmological analyses with $Euclid$. The data used in this analysis are publicly released with the paper., Comment: 33 pages, 24 figures, main results in Fig. 19 & Table 5, version published in A&A

We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extended catalogue of 739 sources either selected from the \emph{Planck} Second Catalogue of Compact Sources or found through a blind search carried out with a Mexican Hat 2 wavelet. A significant fraction of the sources in our catalogue (38.7 per cent) are within the $|b| \leq 20^\circ$ region of the Galactic plane. We determine statistical properties for those sources that are likely to be extragalactic. We find that these statistical properties are compatible with currently available models, with a $\sim$1.8 Jy completeness limit at 11 GHz. We provide the polarimetric properties of (38, 33, 31, 23) sources with P detected above the $99.99\%$ significance level at (11, 13, 17, 19) GHz, respectively. Median polarization fractions are in the $2.8$-$4.7$\% range in the 11-19 GHz frequency interval. We do not distinguish between Galactic and extragalactic sources here. The results presented here are consistent with those reported in the literature for flat- and steep-spectrum radio sources., Comment: 21 pages, 13 figures, 7 tables

We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the WMAP K and Ka bands (23 and 33 GHz) and all Planck polarized channels (30-353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foregrounds, especially the estimation of the synchrotron spectral index, $\beta_s$. We present the first detailed $\beta_s$ map of the Northern Celestial Hemisphere at a smoothing scale of $2^{\circ}$. We find statistically significant spatial variability across the sky. We obtain an average value of $-3.08$ and a dispersion of $0.13$, considering only pixels with reliable goodness-of-fit. The power law model of the synchrotron emission provides a good fit to the data outside the Galactic plane but fails to track the complexity within this region. Moreover, when we assume a synchrotron model with uniform curvature, $c_s$, we find a value of $c_s = -0.0797 \pm 0.0012$. However, there is insufficient statistical significance to determine which model is favoured, either the power law or the power law with uniform curvature. Furthermore, we estimate the thermal dust spectral parameters in polarization. Our CMB, synchrotron, and thermal dust maps are highly correlated with the corresponding products of the PR4 Planck release, although some large-scale differences are observed in the synchrotron emission. Finally, we find that the $\beta_s$ estimation in the high signal-to-noise synchrotron emission areas is prior-independent while, outside these regions, the prior governs the $\beta_s$ estimation., Comment: 25 pages, 29 figures, Accepted for publication in MNRAS

The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ at 11, 13, 17 and 19$\,$GHz. These data are combined with ancillary data to produce Spectral Energy Distributions in intensity in the frequency range 0.4--3\,000$\,$GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation analysis at 1$^\circ$ scale on the full sample from which we identify 44 sources with high AME significance. We explore correlations between different fitted parameters on this last sample. QUIJOTE-MFI data contribute to notably improve the characterisation of the AME spectrum, and its separation from the other components. In particular, ignoring the 10--20\,GHz data produces on average an underestimation of the AME amplitude, and an overestimation of the free-free component. We find an average AME peak frequency of 23.6 $\pm$ 3.6$\,$GHz, about 4$\,$GHz lower than the value reported in previous studies. The strongest correlation is found between the peak flux density of the thermal dust and of the AME component. A mild correlation is found between the AME emissivity ($A_{\rm AME}/\tau_{250}$) and the interstellar radiation field. On the other hand no correlation is found between the AME emissivity and the free-free radiation Emission Measure. Our statistical results suggest that the interstellar radiation field could still be the main driver of the intensity of the AME as regards spinning dust excitation mechanisms. On the other hand, it is not clear whether spinning dust would be most likely associated with cold phases of the interstellar medium rather than with hot phases dominated by free-free radiation., Comment: 36 pages. Accepted for publication in MNRAS. Maps and derived data products available at https://research.iac.es/proyecto/quijote/

The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $\gamma$-ray Fermi bubbles, and with the S-PASS radio polarization plumes, suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic centre and dark matter decay in the Galactic halo, do not yet provide a clear physical interpretation. In this paper we present a re-analysis of the Haze including new observations from the Multi-Frequency Instrument (MFI) of the Q-U-I JOint TEnerife (QUIJOTE) experiment, at 11 and 13 GHz. We analyze the Haze in intensity and polarization, characterizing its spectrum. We detect an excess of diffuse intensity signal ascribed to the Haze. The spectrum at frequencies 11$\,\leq\nu\leq\,$70 GHz is a power-law with spectral index $\beta^{\rm H}=-2.79\pm0.08$, which is flatter than the Galactic synchrotron in the same region ($\beta^{\rm S}=-2.98\pm0.04$), but steeper than that obtained from previous works ($\beta^{\rm H}\sim-2.5$ at 23$\,\leq\,\nu\leq\,$70 GHz). We also observe an excess of polarized signal in the QUIJOTE-MFI maps in the Haze area. This is a first hint detection of polarized Haze, or a consequence of curvature of the synchrotron spectrum in that area. Finally, we show that the spectrum of polarized structures associated with Galactic centre activity is steep at low frequencies ($\beta \sim -3.2$ at 2.3 $\leq\nu\leq$ 23 GHz), and becomes flatter above 11 GHz., Comment: 31 pages. Accepted for publication in MNRAS

We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 ${\mu}$K/beam for total intensity and 13-23 ${\mu}$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7${\sigma}$ and W51 at 4.0${\sigma}$ with peak frequencies ${\nu}_{AME}$ = (20.0 $\pm$ 1.4) GHz and ${\nu}_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices ${\alpha}_s$ = -0.67 $\pm$ 0.10 in W49 and ${\alpha}_s$ = -0.51 $\pm$ 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency ${\nu}_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study., Comment: 29 pages, 11 figures. Accepted for publication in MNRAS

We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polarization within the range 35-40 $\mu$K per 1-degree beam, being a factor $\sim 2$-$4$ worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterise the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5%. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A and 3C274) which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11GHz and WMAP 23GHz map peaks at $\beta=-3.09$ with a standard deviation of 0.14. The measured BB/EE ratio at scales of $\ell=80$ is $0.26\pm 0.07$ for a Galactic cut $|b|>10^\circ$. We find a positive TE correlation for 11GHz at large angular scales ($\ell \lesssim 50$), while the EB and TB signals are consistent with zero in the multipole range $30 \lesssim \ell \lesssim 150$. The maps discussed in this paper are publicly available., Comment: 60 pages. Accepted for publication in MNRAS. Maps and derived data products available at https://research.iac.es/proyecto/quijote/

Broadband receiver data need colour corrections applying to correct for the different source spectra across their wide bandwidths. The full integration over a receiver bandpass may be computationally expensive and redundant when repeated many times. Colour corrections can be applied, however, using a simple quadratic fit based on the full integration instead. Here we describe fastcc and interpcc, quick Python and IDL codes that return, respectively, colour correction coefficients for different power-law spectral indices and modified black bodies for various Cosmic Microwave Background related experiments. The codes are publicly available, and can be easily extended to support additional telescopes., Comment: 3 pages, 1 figure. Published in RNAAS

Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential at any position throughout the image plane, allowing for a nearly model-independent description of the lensing mass. We apply our method on simulated Hubble Space Telescope imaging data containing different types of perturbations to a smooth mass distribution: a localized dark subhalo, a population of subhalos, and an external shear perturbation. Assuming knowledge of the source surface brightness, we use the continuous neural field to model either the perturbations alone or the full lensing potential. In both cases, the resulting model is able to fit the imaging data, and we are able to accurately recover the properties of both the smooth potential and of the perturbations. Unlike many other deep learning methods, ours explicitly retains lensing physics (i.e., the lens equation) and introduces high flexibility in the model only where required, namely, in the lens potential. Moreover, the neural network does not require pre-training on large sets of labelled data and predicts the potential from the single observed lensing image. Our model is implemented in the fully differentiable lens modeling code Herculens.

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), Comment: 37 pages (including appendices), 26 figures; accepted for publication in Astronomy & Astrophysics