KICC papers

arXiv:2505.05450v1 Announce Type: cross Abstract: We study the cosmology of multi-field Dark Energy, using a well-motivated axio-dilaton model that contains the minimal number of fields to have the 2-derivative sigma-model interactions that power-counting arguments show naturally compete with General Relativity at low energies. Our analysis differs from earlier, related, studies by treating the case where the dilaton's couplings to matter are large enough to require screening to avoid unacceptable dilaton-mediated forces in the solar system. We use a recently proposed screening mechanism that exploits the interplay between stronger-than-gravitational axion-matter couplings with the 2-derivative axion-dilaton interactions to suppress the couplings of the dilaton to bulk matter. The required axion-matter couplings also modify cosmology, with the axion's background energy density turning out to resemble early dark energy. We compute the properties of the axion fluid describing the rapid oscillations of the axion field around the time-dependent minimum of its matter-dependent effective potential, extending the usual formalism to include nontrivial kinetic sigma-model interactions. We explore the implications of these models for the Cosmic Microwave Background and the growth of structure and find that for dilaton potentials of the Albrecht-Skordis form (itself well-motivated by UV physics), successful screening can be consistent with the early dark energy temporarily comprising as much as 10% of the total density in the past. We find that increasing the dilaton-matter coupling decreases the growth of structure due to enhanced Hubble friction, an effect that dominates the usual fifth-force effects that amplify structure growth.

arXiv:2412.08809v2 Announce Type: replace Abstract: We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards ($16.5 < V < 19.5$) alongside three CALSPEC standards, from 912 \r{A} to 32 $\mu$m. The framework is the first of its kind to jointly infer photometric zeropoints and WD parameters (surface gravity $\log g$, effective temperature $T_{\text{eff}}$, extinction $A_V$, dust relation parameter $R_V$) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy and ground-based optical spectroscopy to sub-percent precision. Photometric residuals for the sample are the lowest yet yielding $<0.004$ mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory's Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

arXiv:2505.04610v1 Announce Type: new Abstract: We present ultraviolet (UV) through near-infrared (NIR) photometric and spectroscopic observations of the nearby SN 2024pxl, the third Type Ia supernova (SN Ia) in NGC 6384. SN 2024pxl is a Type Iax supernova (SN Iax) with an intermediate luminosity ($M_r = -16.99\pm0.32$ mag) and an average SN Iax light curve decline rate. SN 2024pxl was discovered $\sim$3 days after first light, and the rising light curve follows a single power law that is inconsistent with significant interaction with a companion star or circumstellar material. Our extensive NIR photometric coverage is comparable to that of the well-observed SNe Iax 2005hk and 2012Z, and we demonstrate that the $J-H$ colors of SNe Iax differ from normal SNe Ia and appear to be more homogeneous as a class. Spectroscopically, we report the earliest-ever NIR spectrum of a SN Iax as measured from maximum light ($t\approx-9$ days): a featureless continuum with similarities to a $\sim$9,000 K blackbody, and the line velocities are consistent with a mixed-ejecta structure, with C, Si, and Fe having similar velocities and velocity evolutions. We find a tentative correlation between the $H$-band break Co II velocity $\sim$20 days post-peak and absolute magnitude, with more luminous SNe Iax showing faster Co II velocities. Our observations suggest that SN 2024pxl resulted from the thermonuclear disruption of a CO white dwarf star that undergoes deflagration burning.

arXiv:2505.03873v1 Announce Type: new Abstract: Population III stars, the hypothetical first generation metal-free stars, have yet to be discovered. Even after three years of successful JWST operations, studies have shown that most galaxies identified to date at $z > 5$ exhibit a metallicity floor of $Z\gtrsim2\%\,Z_{\odot}$, possibly due to unknown selection biases toward bright galaxies or rapid metal enrichment. To address this question, we search for galaxies with low R3$=$[OIII]$\,\lambda$5008/H$\beta$ ratios as part of the JWST Cycle-3 large treasury program, the Slitless Areal Pure-Parallel HIgh-Redshift Emission Survey (SAPPHIRES). Using deep NIRCam Wide-Field Slitless Spectrscopy (WFSS) data, we report the discovery of seven extremely metal-poor galaxy candidates in the SAPPHIRES Early Data Release field, with estimated $12+{\rm log(O/H)}<7.0$ at $z\sim5-7$, including two sources with $Z<1\%\,Z_{\odot}$, significantly breaking the metallicity floor observed both locally and at high redshift. These candidates appear extremely faint ($\sim28-30\,$ F200W AB mag) and low-mass (${\rm log}(M_{*}/M_{\odot})\sim6.8-7.8$), as expected from the mass-metallicity relation. They also exhibit very blue UV slopes ($-2.6\lesssim\beta\lesssim-2.0$), likely due to low dust content $A_{V}\lesssim0.2\,{\rm mag}$ or young stellar ages $\sim5-20\,{\rm Myr}$. Compared to galaxies at similar redshift, they appear exceptionally bursty in their star formation activity. Our results highlight the power of NIRCam/WFSS in identifying extremely metal-poor galaxies, from just a single pointing, with more data to come in SAPPHIRES. This underscores the potential of pure-parallel programs towards achieving JWST's primary science goal: discovering the first pristine stars and galaxies. Deep JWST/NIRSpec follow-up observations will also be essential to confirm their nature and perform detailed chemical abundance analyses.

arXiv:2505.03884v1 Announce Type: new Abstract: We present ALMA CO observations of the molecular gas in a sample of 41 luminous unobscured quasars at z $\sim$ 2 from the Sloan Digital Sky Survey. 32 targets comprise the main sample observed in CO(3-2) and 9 targets have archival ALMA data of CO(3-2), CO(4-3) and CO(7-6). All quasars have rest-UV to optical spectra tracing ionised gas in the broad line region (e.g. CIV) and the narrow line region (e.g. [OIII]) and they cover the full range of outflow properties in the SDSS quasar population at these redshifts. 15 out of the 32 quasars in the main sample are detected in CO(3-2) and five out of the nine archival quasars are also detected in CO. The median gas mass for all 20 CO detected quasars is 8.0 $\pm$ 1.5 $\times$ 10$^9$ M$_{\odot}$ with a median M$_{dyn}$ of 1.4 $\pm$ 0.9 $\times$ 10$^{11}$ M$_{\odot}$. We find gas fractions in the range 0.02 - 0.32, which are generally lower than both inactive galaxies and obscured quasars at similar redshifts. We suggest an evolutionary trend in gas fractions of quasar host galaxies from obscured and gas rich to unobscured and gas poor. We note a tentative correlation between the gas fractions and the broad-line region properties with quasars showing high CIV blueshifts, indicating stronger broad-line region winds, having higher gas fractions. Six of the quasars corresponding to 15% of the sample also show evidence for at least one companion galaxy detected in CO at the same redshift.

arXiv:2501.08969v2 Announce Type: replace Abstract: In this work we estimate the explosion and progenitor properties of six Type II supernovae (SNe) at 0.675 <= z <= 3.61 discovered by the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES) transient survey by modeling their light curves. Two Type II SNe are found to have high explosion energies of 3e51 erg, while the other four Type II SNe are estimated to have typical explosion energies found in the local Universe [(0.5-2)e51 erg]. The fraction of Type II SNe with high explosion energies might be higher at high redshifts because of, e.g., lower metallicity, but it is still difficult to draw a firm conclusion because of the small sample size and potential observational biases. We found it difficult to constrain the progenitor masses for Type II SNe in our sample because of the sparse light-curve data. We found two Type II SN light curves can be better reproduced by introducing confined, dense circumstellar matter. Thus, the confined, dense circumstellar matter frequently observed in nearby Type II SNe is likely to exist in Type II SNe at high redshifts as well. Two Type II SNe are estimated to have high host galaxy extinctions, showing the ability of JWST to discover dust-obscured SNe at high redshifts. More high-redshift Type II SNe are required to investigate the differences in the properties of Type II SNe near and far, but here we show the first glimpse into the high-redshift population of Type II SNe.

arXiv:2505.02895v1 Announce Type: new Abstract: We present a spectroscopic search for broad H$\alpha$ emitters at z~3.7-6.5 in the GOODS-N field, utilizing JWST/NIRCam slitless spectroscopy from FRESCO and CONGRESS, complemented by JADES imaging. We identify 19 broad H$\alpha$ emitters with FWHM > 1000 km/s at z~4-5.5, including 9 new sources. The black hole masses and AGN bolometric luminosities, inferred from the broad H$\alpha$ components, indicate that most sources are accreting at ~10% of the Eddington limit. We derive their host stellar masses via SED fitting and find higher $M_{BH}/M_{*}$ ratios relative to the local $M_{BH}$-$M_{*}$ relations, consistent with previous studies. We find that 42% of the sample do not satisfy the widely-used color selection criteria for Little Red Dots (LRDs), with the majority of these sources lacking the characteristic steep red slope. A comparison of the average SEDs between our sample and LRDs selected in the same field reveals that the steep red slopes observed in some LRDs are likely due to line-boosting effects as previously suggested. Furthermore, we find that 68% of color-selected LRDs with H$\alpha$ detections in the NIRCam/Grism spectra do not exhibit broad-line features. While the limited sensitivity of the grism spectra may hinder the detection of broad-line components in faint sources, our findings still highlight the enigmatic nature of the LRD population.

arXiv:2505.02896v1 Announce Type: new Abstract: We study the large-scale environments and clustering properties of 28 low-luminosity AGNs at $z=3.9-6$ in the GOODS-N field. Our sample, identified from the JWST NIRCam Imaging and WFSS data in CONGRESS and FRESCO surveys with either broad H$\alpha$ emission lines or V-shape continua, are compared to 782 H$\alpha$ emitters (HAEs) selected from the same data. These AGNs are located in diverse large-scale environments and do not preferentially reside in denser environments compared to HAEs. Their overdensity field, $\delta$, averaged over (15 $h^{-1}$cMpc)$^3$, ranges from $-0.56$ to 10.56, and shows no clear correlation with broad-line luminosity, black hole (BH) masses, or the AGN fraction. It suggests that $> 10$ cMpc structures do not significantly influence BH growth. We measure the two-point cross-correlation function of AGNs with HAEs, finding a comparable amplitude to that of the HAE auto-correlation. This indicates similar bias parameters and host dark matter halo masses for AGNs and HAEs. The correlation length of field AGNs is 4.26 $h^{-1}$cMpc, and 7.66 $h^{-1}$cMpc at $3.9 < z < 5$ and $5 < z < 6$, respectively. We infer a median host dark matter halo mass of $\log (M_h/M_\odot)\approx 11.0-11.2$ and host stellar masses of $\log (M_\star/M_\odot) \approx 8.4-8.6$ by comparing with the UniverseMachine simulation. Our clustering analysis suggests that low-luminosity AGNs at high redshift reside in normal star-forming galaxies with overmassive BHs. They represent an intrinsically distinct population from luminous quasars and could be a common phase in galaxy evolution.

arXiv:2505.03737v1 Announce Type: new Abstract: Weak gravitational lensing of the CMB has been established as a robust and powerful observable for precision cosmology. However, the impact of Galactic foregrounds, which has been studied less extensively than many other potential systematics, could in principle pose a problem for CMB lensing measurements. These foregrounds are inherently non-Gaussian and hence might mimic the characteristic signal that lensing estimators are designed to measure. We present an analysis that quantifies the level of contamination from Galactic dust in lensing measurements, focusing particularly on measurements with the Atacama Cosmology Telescope and the Simons Observatory. We employ a whole suite of foreground models and study the contamination of lensing measurements with both individual frequency channels and multifrequency combinations. We test the sensitivity of different estimators to the level of foreground non-Gaussianity, and the dependence on sky fraction and multipole range used. We find that Galactic foregrounds do not present a problem for the Atacama Cosmology Telescope experiment (the bias in the inferred CMB lensing power spectrum amplitude remains below $0.3\sigma$). For Simons Observatory, not all foreground models remain below this threshold. Although our results are conservative upper limits, they suggest that further work on characterizing dust biases and determining the impact of mitigation methods is well motivated, especially for the largest sky fractions.

arXiv:2505.02233v1 Announce Type: new Abstract: We investigate the degeneracy between the effects of ultra-light axion dark matter and baryonic feedback in suppressing the matter power spectrum. We forecast that galaxy shear data from the Rubin Observatory's Legacy Survey of Space and Time (LSST) could limit an axion of mass $m = 10^{-25}\,\mathrm{eV}$ to be $\lesssim 5\%$ of the dark matter, stronger than any current bound, if the interplay between axions and feedback is accurately modelled. Using a halo model emulator to construct power spectra for mixed cold and axion dark matter cosmologies, including baryonic effects, we find that galaxy shear is sensitive to axions from $10^{-27}\,\mathrm{eV}$ to $10^{-21}\,\mathrm{eV}$, with the capacity to set competitive bounds across much of this range. For axions with $m \sim 10^{-25}\,\mathrm{eV}$, the scales at which axions and feedback impact structure formation are similar, introducing a parameter degeneracy. We find that, with an external feedback constraint, we can break the degeneracy and constrain the axion transfer function, such that LSST could detect a $10^{-25}\,\mathrm{eV}$ axion comprising 10\% of the dark matter at $\sim 3 \sigma$ significance. Direct reconstruction of the non-linear matter power spectrum provides an alternative way of analysing weak lensing surveys, with the advantage of identifying the scale-dependent features in the data that the dark matter model imposes. We advocate for dedicated cosmological hydrodynamical simulations with an axion dark matter component so that upcoming galaxy and cosmic microwave background lensing surveys can disentangle the dark matter-baryon transfer function.

arXiv:2409.11464v3 Announce Type: replace Abstract: For the first time, we present a systematic search for galaxies with extended emission line and potential outflow features using \textit{JWST} medium-band images in the GOODS-S field. This is done by comparing the morphology in medium-band images to adjacent continuum and UV bands. We look for galaxies that have a maximum extent 50\% larger, an excess area 30\% greater, or an axis ratio difference of more than 0.3 in the medium band compared to the reference bands. After visual inspection, we find 326 candidate galaxies at $1.4 < z < 8.4$, with a peak in the population near cosmic noon, benefiting from the good coverage of the medium-band filters. By fitting their SEDs, we find that the candidate galaxies are at least 20\% more bursty in their star-forming activity and have 50\% more young stellar populations compared to a control sample selected based on the continuum band flux. Additionally, these candidates exhibit a significantly higher production rate of ionizing photons. We further find that candidates hosting known AGN produce extended emission that is more anisotropic compared to non-AGN candidates. A few of our candidates have been spectroscopically confirmed to have prominent outflow signatures through NIRSpec observations, showcasing the robustness of the photometric selection. Future spectroscopic follow-up will better help verify and characterize the kinematics and chemical properties of these systems.

arXiv:2505.01512v1 Announce Type: new Abstract: The discovery of inhabited exoplanets hinges on identifying biosignature gases. JWST is revealing potential biosignatures in exoplanet atmospheres, though their presence is yet to provide strong evidence for life. The central challenge is attribution: how to confidently identify biogenic sources while ruling out, or deeming unlikely, abiotic explanations? Attribution is particularly difficult for individual planets, especially regarding system-scale stochastic processes that could set atmospheric conditions. To address this, we here propose a comparative multi-planet approach: comparing atmospheric compositions across multiple planets within a system to empirically define the 'abiotic baseline'. This baseline serves as a reference point for biosignatures, and enables marginalisation over inaccessible, shared abiotic parameters. This is possible because planets within a system are linked by their birth in the same natal disk, having been irradiated by the same evolving star, and having a related dynamical history. Observations aligning with the abiotic baseline, where the locally informed abiotic models demonstrate high out-of-sample predictive accuracy, are likely non-biological. Deviations from the baseline -- potentially biotic anomalies -- suggest an alternative origin. We present the application of Bayesian leave-one-out cross-validation to evaluate the performance of geochemical- and biogeochemical-climate models in explaining these anomalies, using the expected log pointwise predictive density as a diagnostic. When biogeochemical models outperform their abiotic counterparts, the anomaly may be shaped by life, and constitutes a comparative biosignature. If both models perform poorly, the anomaly is flagged as an "unknown unknown" -- a signature of either unrecognised abiotic chemistry, or life as we don't yet know it.

arXiv:2410.19905v2 Announce Type: replace Abstract: Energetic feedback processes associated with accreting supermassive black holes can expel gas from massive haloes and significantly alter various measures of clustering on ~Mpc scales, potentially biasing the values of cosmological parameters inferred from analyses of large-scale structure (LSS) if not modelled accurately. Here we use the state-of-the-art FLAMINGO suite of cosmological hydrodynamical simulations to gauge the impact of feedback on large-scale structure by comparing to Planck + ACT stacking measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect of SDSS BOSS galaxies. We make careful like-with-like comparisons to the observations, aided by high precision KiDS and DES galaxy-galaxy lensing measurements of the BOSS galaxies to inform the selection of the simulated galaxies. In qualitative agreement with several recent studies using dark matter only simulations corrected for baryonic effects, we find that the kSZ effect measurements prefer stronger feedback than predicted by simulations which have been calibrated to reproduce the gas fractions of low redshift X-ray-selected groups and clusters. We find that the increased feedback can help to reduce the so-called S8 tension between the observed and CMB-predicted clustering on small scales as probed by cosmic shear (although at the expense of agreement with the X-ray group measurements). However, the increased feedback is only marginally effective at reducing the reported offsets between the predicted and observed clustering as probed by the thermal SZ (tSZ) effect power spectrum and tSZ effect--weak lensing cross-spectrum, both of which are sensitive to higher halo masses than cosmic shear.

arXiv:2505.00095v1 Announce Type: new Abstract: Metallicity ($Z$) estimates based on ultraviolet (UV) emission lines from the narrow-line regions (NLRs) of active galactic nuclei (AGNs) have been found to differ from those derived from optical lines. However, the origin of this discrepancy ($ZR$) remains poorly understood. To investigate the source of $ZR$, we compiled from the literature the fluxes of narrow near-UV ($1000 < \lambda(\angstrom) < 2000)$ and optical ($3000 < \lambda(\angstrom) < 7000)$ emission line measurements for a sample of 11 AGNs (9 at $z<0.4$ and 2 at $z\sim2.4$). Metallicity values for our sample were derived using a semi-empirical calibration based on the $C43$=log[(\ion{C{iv}$\lambda$1549+\ion{C{iii}]$\lambda$1909)/\ion{He}{ii}$\lambda$1640] emission-line ratio and compared with those obtained via direct measurement of the electron temperature ($T_{\rm e}$-method) and via calibrations based on optical emission-lines. The source of the discrepancy was investigated in terms of the ionization parameter ($U$), electron density ($N_{\rm e}$), and carbon abundance (C/H). We found a weak correlation between $ZR$, $U$ and $N_{\rm e}$. However, a moderate correlation was observed between $ZR$ and direct estimates of C/H, suggesting that the previously assumed (C/O)-$Z$ relations in photoionization models used to derive UV carbon-line calibrations may not be valid for AGNs. By combining a large set of abundance estimates for local star-forming regions with those of our AGN sample, we derived a new (C/O)-$Z$ relation. Comparisons between the results of photoionization models that assume this new abundance relation and the UV observational data of our sample produce $Z$ values derived from the $C43$ index that are consistent with those obtained using the $T_{\rm e}$-method.

arXiv:2409.20379v2 Announce Type: replace Abstract: A precise measurement of photometric redshifts (photo-z) is key for the success of modern photometric galaxy surveys. Machine learning (ML) methods show great promise in this context, but suffer from covariate shift (CS) in training sets due to selection bias where interesting sources are underrepresented, and the corresponding ML models show poor generalisation properties. We present an application of the StratLearn method to the estimation of photo-z, validating against simulations where we enforce the presence of CS to different degrees. StratLearn is a statistically principled approach that relies on splitting the source and target datasets into strata based on estimated propensity scores (i.e. the probability for an object to be in the source set given its observed covariates). After stratification, two conditional density estimators are fit separately to each stratum, then combined via a weighted average. We benchmark our results against the GPz algorithm, quantifying the performance of the two codes with a set of metrics. Our results show that the StratLearn-z metrics are only marginally affected by the presence of CS, while GPz shows a significant degradation of performance in the photo-z prediction for fainter objects. For the strongest CS scenario, StratLearn-z yields a reduced fraction of catastrophic errors, a factor of 2 improvement for the RMSE and one order of magnitude improvement on the bias. We also assess the quality of the conditional redshift estimates with the probability integral transform (PIT). The PIT distribution obtained from StratLearn-z features fat fewer outliers and is symmetric, i.e. the predictions appear to be centered around the true redshift value, despite showing a conservative estimation of the spread of the conditional redshift distributions. Our julia implementation of the method is available at https://github.com/chiaramoretti/StratLearn-z.

arXiv:2504.21087v1 Announce Type: new Abstract: The Einstein Telescope (ET), along with other third-generation gravitational wave (GW) detectors, will be a key instrument for detecting GWs in the coming decades. However, analyzing the data and estimating source parameters will be challenging, especially given the large number of expected detections - of order $10^5$ per year - which makes current methods based on stochastic sampling impractical. In this work, we use Dingo-IS to perform Neural Posterior Estimation (NPE) of high-redshift events detectable with ET in its triangular configuration. NPE is a likelihood-free inference technique that leverages normalizing flows to approximate posterior distributions. After training, inference is fast, requiring only a few minutes per source, and accurate, as corrected through importance sampling and validated against standard Bayesian inference methods. To confirm previous findings on the ability to estimate parameters for high-redshift sources with ET, we compare NPE results with predictions from the Fisher information matrix (FIM) approximation. We find that FIM underestimates sky localization errors substantially for most sources, as it does not capture the multimodalities in sky localization introduced by the geometry of the triangular detector. FIM also overestimates the uncertainty in luminosity distance by a factor of $\sim 3$ on average when the injected luminosity distance is $d^{\mathrm{inj}}_{\mathrm{L}} > 10^5~$Mpc, further confirming that ET will be particularly well suited for studying the early Universe.

arXiv:2504.20144v1 Announce Type: cross Abstract: We report the existence of two exotic compact objects in the leading relativistic model of modified Newtonian dynamics, namely aether-scalar-tensor theory. This model is consistent with precision cosmology and gravitational wave constraints on tensor speed. Black hole mimickers could subtly change observations: gravitational waves from their mergers might show unusual echoes or altered ringdown patterns, and images of their horizon-scale shadows might be slightly different from those of a true black hole. Shapiro-free lenses are massless objects that deflect light without any gravitational time delay, producing distinctive lensing events. These predictions connect to ongoing and future gravitational-wave searches, horizon-scale imaging, and time-domain lensing surveys.

arXiv:2503.17342v2 Announce Type: replace Abstract: We present an updated reconstruction of the dark energy equation of state, $w(a)$, using the newly released DESI DR2 Baryon Acoustic Oscillation (BAO) data in combination with Pantheon+ and DES5Y Type Ia supernovae measurements, respectively. Building on our previous analysis in arXiv:2503.08658, which employed a nonparametric flexknot reconstruction approach, we examine whether the evidence for dynamical dark energy persists with the improved precision of the DESI DR2 dataset. We find that while the overall qualitative structure of $w(a)$ remains consistent with our earlier findings, the statistical support for dynamical dark energy is reduced when considering DESI DR2 data alone, particularly for more complex flexknot models with higher numbers of knots. However, the evidence for simpler dynamical models, such as $w$CDM and CPL (which correspond to $n=1$ and $n=2$ knots respectively), increases relative to $\Lambda$CDM with DESI DR2 alone, with CPL being the preferred dynamical model, consistent with previous DESI analyses. When combined with Pantheon+ data, the conclusions remain broadly consistent with our earlier work, but when instead combined with DES5Y supernovae data, there is an increased preference for flexknot models for all values of $n$ considered. This results in all such models being preferred over $\Lambda$CDM, with the CPL model being the most favoured by a Bayes factor of $\sim 2.3$ relative to $\Lambda$CDM.

arXiv:2504.20992v1 Announce Type: new Abstract: We analyse the large-scale angular clustering of quasars in the Gaia-unWISE quasar catalog, Quaia, and their cross-correlation with maps of the lensing convergence of the Cosmic Microwave Background (CMB), to constrain the level of primordial non-Gaussianity (PNG). Specifically, we target the scale-dependent bias that would be induced by PNG on biased tracers of the matter inhomogeneities on large scales. The Quaia sample is particularly well suited for this analysis, given the large effective volume covered, and our ability to map out the main potential sources of systematic contamination and mitigate their impact. Using the universality relation to characterise the response of the quasar overdensity to PNG ($p_\phi=1$), we report constraints on the local-type PNG parameter $f_{\rm NL}$ of $f_{\rm NL}=-20.5^{+19.0}_{-18.1}$ (68\% C.L.) by combining the quasar auto-correlation and its cross-correlation with CMB lensing in two tomographic redshift bins (or $f_{\rm NL}=-28.7^{+26.1}_{-24.6}$ if assuming a lower response for quasars, $p_\phi=1.6$). Using the CMB lensing cross-correlations alone, we find $f_{\rm NL}=-13.8^{+26.7}_{-25.0}$. These are the tightest constraints on $f_{\rm NL}$ to date from angular clustering statistics and cross-correlations with CMB lensing.