SEARCH · Engineering Papers

Results for “LSST”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Testing cosmology with double source lensing

Double source lensing provides a dimensionless ratio of distance ratios, a "remote viewing" of cosmology through distances relative to the gravitational lens, beyond the observer. We use this to test the cosmological framework, particularly with respect to spatial curvature and the distance duality relation. We derive a consistency equation for constant spatial curvature, allowing not only the investigation of flat vs curved but of the Friedmann-Lemaître-Robertson-Walker framework itself. For distance duality, we demonstrate that the evolution of the lens mass profile slope must be controlled to ≳ 5 times tighter fractional precision than a claimed distance duality violation. Furthermore, using LensPop forecasts of double source lensing systems in Euclid and LSST surveys we also explore constraints on dark energy equation of state parameters and any evolution of the lens mass profile slope.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Primordial non-Gaussianities with weak lensing: information on non-linear scales in the Ulagam full-sky simulations

Abstract Primordial non-Gaussianities (PNGs) are signatures in the density field that encode particle physics processes from the inflationary epoch. Such signatures have been extensively studied using the Cosmic Microwave Background, through constraining their amplitudes,f X NL , with future improvements expected from large-scale structure surveys; specifically, the galaxy correlation functions. We show that weak lensing fields can be used to achieve competitive and complementary constraints. This is shown via theUlagamsuite of N-body simulations, a subset of which evolves primordial fields with four types of PNGs. We create full-sky lensing maps and estimate the Fisher information from three summary statistics measured on the maps: the moments, the cumulative distribution function, and the 3-point correlation function. We find that the year 10 sample from the Rubin Observatory Legacy Survey of Space and Time (LSST) can constrain PNGs toσ(f NL eq ) ≈ 110,σ(f NL or, lss ) ≈ 120,σ(f NL loc ) ≈ 40. For the former two, this is better than or comparable to expected galaxy clustering-based constraints from the Dark Energy Spectroscopic Instrument (DESI). The PNG information in lensing fields is on non-linear scales and at low redshifts (z≲ 1.25), with a clear origin in the evolution history of massive halos. The constraining power degrades by ∼60% under scale cuts of ≳ 20 Mpc, showing there is still significant information on scales mostly insensitive to small-scale systematic effects (e.g., baryons). We publicly release theUlagamsuite to enable more survey-focused analyses.

Astronomy & Astrophysics↗

DOE OSTI2024

Light fields during inflation from BOSS and future galaxy surveys

Abstract Primordial non-Gaussianity generated by additional fields present during inflation offers a compelling observational target for galaxy surveys. These fields are of significant theoretical interest since they offer a window into particle physics in the inflaton sector. They also violate the single-field consistency conditions and induce a scale-dependent bias in the galaxy power spectrum. In this paper, we explore this particular signal for light scalar fields and study the prospects for measuring it with galaxy surveys. We find that the sensitivities of current and future surveys are remarkably stable for different configurations, including between spectroscopic and photometric redshift measurements. This is even the case at non-zero masses where the signal is not obviously localized on large scales. For realistic galaxy number densities, we demonstrate that the redshift range and galaxy bias of the sample have the largest impact on the sensitivity in the power spectrum. These results additionally motivated us to explore the potentially enhanced sensitivity of Vera Rubin Observatory's LSST through multi-tracer analyses. Finally, we apply this understanding to current data from the last data release of the Baryon Oscillation Spectroscopic Survey (BOSS DR12) and place new constraints on light fields coupled to the inflaton.

Astronomy & Astrophysics↗

DOE OSTI2024

On the detectability of the moving lens signal in CMB experiments

Abstract Upcoming cosmic microwave background (CMB) experiments are expected to detect new signals probing interaction of CMB photons with intervening large-scale structure. Among these the moving-lens effect, the CMB temperature anisotropy induced by cosmological structures moving transverse to our line of sight, is anticipated to be measured to high significance in the near future. In this paper, we investigate two possible strategies for the detection of this signal: pairwise transverse-velocity estimation and oriented stacking. We expand on previous studies by including in the analysis realistic simulations of competing signals and foregrounds. We confirm that the moving lens effect can be detected at ≥ 10σlevel by a combination of CMB-S4 and LSST surveys. We show that the limiting factors in the detection depend on the strategy: for the stacking analysis, correlated extragalactic foregrounds, namely the cosmic infrared background and thermal Sunyaev Zel'dovich effect, play the most important role. The addition of foregrounds make the signal-to-noise ratio be most influenced by large and nearby objects. As for the pairwise detection, halo lensing and pair number counts are the main issues. In light of our findings, we elaborate on possible strategies to improve the analysis approach for the moving lens detection with upcoming experiments. We also deliver to the community all the simulations and tools we developed for this study.

Astronomy & Astrophysics↗

DOE OSTI2025

Aemulus ν: precision halo mass functions in wνCDM cosmologies

Precise and accurate predictions of the halo mass function for cluster mass scales in wνCDM cosmologies are crucial for extracting robust and unbiased cosmological information from upcoming galaxy cluster surveys. Here, we present a halo mass function emulator for cluster mass scales (≳ 1013 M ⊙/h) up to redshift z = 2 with comprehensive support for the parameter space of wνCDM cosmologies allowed by current data. Based on the Aemulus ν suite of simulations, the emulator marks a significant improvement in the precision of halo mass function predictions by incorporating both massive neutrinos and non-standard dark energy equation of state models. This allows for accurate modeling of the cosmology dependence in large-scale structure and galaxy cluster studies. We show that the emulator, designed using Gaussian Process Regression, has negligible theoretical uncertainties compared to dominant sources of error in future cluster abundance studies. Our emulator is publicly available (https://github.com/DelonShen/aemulusnu_hmf), providing the community with a crucial tool for upcoming cosmological surveys such as LSST and Euclid.

cluster counts↗

DOE OSTI2025

The HalfDome multi-survey cosmological simulations: N-body simulations

Upcoming cosmological surveys have the potential to reach groundbreaking discoveries on multiple fronts, including the neutrino mass, dark energy, and inflation. Most of the key science goals require the joint analysis of datasets from multiple surveys to break parameter degeneracies and calibrate systematics. To realize such analyses, a large set of mock simulations that realistically model correlated observables is required. In this paper we present the N-body component of the HalfDome cosmological simulations, designed for the joint analysis of Stage-IV cosmological surveys, such as Rubin LSST, Euclid, SPHEREx, Roman, DESI, PFS, Simons Observatory, CMB-S4, and LiteBIRD. Our 300TB initial data release includes full-sky lightcones and halo catalogs between z = 0–4 for 11 fixed cosmology realizations, as well as an additional run with local primordial non-Gaussianity (f NL = 20). The simulations evolve 61443 particles in a 3.75 h -1 Gpc box, reaching a minimum halo mass of ∼6 × 1012 h -1 M ⊙ and maximum scale of k ∼ h Mpc-1. Our data is publicly available: instructions to access the data and plans for future data releases can be found at https://halfdomesims.github.io.

Bayer, Adrian E↗

DOE OSTI2025

Selection of high-redshift Lyman-Break Galaxies from broadband and wide photometric surveys

Here, in this paper, we investigate the possibility of selecting high-redshift Lyman-Break Galaxies (LBG) using current and future broadband wide photometric surveys, such as the Ultraviolet Near Infrared Optical Northern Survey (UNIONS) or the Vera C. Rubin Legacy Survey of Space and Time (LSST), using a Random Forest algorithm. This work is conducted in the context of future large-scale structure spectroscopic surveys like DESI-II, the next phase of the Dark Energy Spectroscopic Instrument (DESI), which will start around 2029. We use deep imaging data from the Hyper Suprime Camera (HSC) and the Canada-France-Hawaii Telescope Large Area U-band Deep Survey (CLAUDS) on the COSMOS and XMM-LSS fields. To predict the selection performance of LBGs with image quality similar to UNIONS, we degrade the u,g,r,i and z bands to UNIONS depth. The Random Forest algorithm is trained with the u,g,r,i and z bands to classify LBGs in the 2.5 < z < 3.5 range. We find that fixing a target density budget of 1,100 deg -2 , the Random Forest approach gives a density of z > 2 targets of 873 deg -2 , and a density of 493 deg -2 of confirmed LBGs after spectroscopic confirmation with DESI. This UNIONS-like selection was tested in a dedicated spectroscopic observation campaign of 1,000 targets with DESI on the COSMOS field, providing a safe spectroscopic sample with a mean redshift of 3. This sample is used to derive forecasts for DESI-II, assuming a sky coverage of 5,000 deg 2 . We predict uncertainties on Alcock-Paczynski parameters α ⊥ and α ∥ to be 0.7% and 1% for 2.6 < z < 3.2, resulting in a potential 2% measurement of the dark energy fraction at high redshift. Additionally, we estimate the uncertainty in local non-Gaussianity and predict σ fNL ≈ 7, which would be comparable to the current best precision achieved by Planck. The latter forecast suggests that achieving the precision required to place stringent constraints on inflationary models (σ fNL ≈ 1) using spectroscopic galaxy surveys necessitates the development of a next-generation (Stage V) spectroscopic survey.

cosmological parameters from LSS↗

DOE OSTI2025

A direct detection method of galaxy intrinsic ellipticity-gravitational shear correlation in non-linear regimes using self-calibration

Intrinsic alignment (IA) of galaxies is a challenging source of contamination in the Cosmic shear (GG) signals. The galaxy intrinsic ellipticity-gravitational shear (IG) correlation is generally the most dominant component of such contamination for cross-correlating redshift bins. One of the most effective techniques to mitigate such contamination is the self-calibration (SC) method which extracts the IG correlation and allows for its removal from the GG signal. In a photometric survey, the SC method first extracts the galaxy number density-galaxy intrinsic ellipticity (gI) correlation from the observed galaxy-galaxy lensing correlation using the redshift dependence of lens-source pairs. The IG correlation is computed through a scaling relation using the gI correlation and other lensing observables. The applicability of the SC method has so far been focused on the linear IA scales and the linear galaxy bias. We extend the SC method beyond the linear regime by modifying its scaling relation which can account for the non-linear galaxy bias model and various IA models. In this study, we provide a framework to detect the IG correlation for the redshift bins for source galaxies for the proposed year 1 survey of the Rubin Legacy Survey of Space and Time (LSST Y1). We tested the method for the tidal alignment and tidal torquing (TATT) model of IA and we found that the scaling relation is accurate within 10% and 20% for cross-correlating and auto-correlating redshift bins, respectively. Hence the suppression of IG contamination in observed GG correlation can be accomplished with a factor of 10 and 5, for cross-correlating and auto-correlating redshift bins, respectively. We tested the method's robustness and found that the suppression of IG contamination by a factor of 5 is still achievable for all combinations of cross-correlating bins even with the inclusion of a moderate amount of uncertainties on IA and bias parameters, respectively. We also make available, a branch of the code FAST-PT to provide gI correlations up to 1-loop order term used by the new SC method.

gravitational lensing↗

DOE OSTI2025

Analyzing line-of-sight selection biases in galaxy-scale strong lensing with external convergence and shear

The upcoming Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will dramatically increase the number of strong gravitational lensing systems, requiring precise modeling of line-of-sight (LOS) effects to mitigate biases in lensing observations and cosmological inferences. We develop a method to construct joint distributions of external convergence (κ ext ) and shear (γ ext ) for strong lensing LOS by aggregating large-scale structure simulations with high-resolution halo renderings and non-linear correction. Our approach captures both smooth background matter and perturbations from halos, enabling accurate modeling of LOS effects. Here, we apply non-linear LOS corrections to κ ext and γ ext that address the non-additive lensing effects caused by objects along the LOS in strong lensing. We find that, with a minimum image separation of 1.0'', non-linear LOS correction due to the presence of a dominant deflector slightly increases the ratio of quadruple to double lenses; this non-linear LOS correction also introduces systematic biases of ∼ 0.1% for galaxy-AGN lens in the inferred Hubble constant (H 0 ) if not accounted for. We also observe a 0.66% bias for galaxy-galaxy lenses on H 0 , and even larger biases 1.02% for galaxy-AGN systems if LOS effects are not accounted for. These results highlight the importance of LOS for precision cosmology. The publicly available code and datasets provide tools for incorporating LOS effects in future analyses.

Hubble constant↗

DOE OSTI2022

The Sensitivity of GPz Estimates of Photo-z Posterior PDFs to Realistically Complex Training Set Imperfections

The accurate estimation of photometric redshifts is crucial to many upcoming galaxy surveys, for example, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Almost all Rubin extragalactic and cosmological science requires accurate and precise calculation of photometric redshifts; many diverse approaches to this problem are currently in the process of being developed, validated, and tested. In this work, we use the photometric redshift code GPz to examine two realistically complex training set imperfections scenarios for machine learning based photometric redshift calculation: (i) where the spectroscopic training set has a very different distribution in color–magnitude space to the test set, and (ii) where the effect of emission line confusion causes a fraction of the training spectroscopic sample to not have the true redshift. By evaluating the sensitivity of GPz to a range of increasingly severe imperfections, with a range of metrics (both of photo-z point estimates as well as posterior probability distribution functions, PDFs), we quantify the degree to which predictions get worse with higher degrees of degradation. In particular, we find that there is a substantial drop-off in photo-z quality when line-confusion goes above ~1%, and sample incompleteness below a redshift of 1.5, for an experimental setup using data from the Buzzard Flock synthetic sky catalogs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Modelling the Milky Way – I. Method and first results fitting the thick disc and halo with DES-Y3 data

ABSTRACT We present a technique to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) generated from trilegal models to real data. We apply this technique, which we call mwfitting, to photometric data from the first 3 yr of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud, and the Sagittarius Stream, our main sample spans a total area of ∼2300 deg2. We further explore a smaller subset (∼1300 deg2) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3 per cent. We fit the DES data with an exponential thick disc model and an oblate double power-law halo model. We find that the best-fitting thick disc model has radial and vertical scale heights of 2.67 ± 0.09 kpc and 925 ± 40 pc, respectively. The stellar halo is fit with a broken power-law density profile with an oblateness of 0.75 ± 0.01, an inner index of 1.82 ± 0.08, an outer index of 4.14 ± 0.05, and a break at 18.52 ± 0.27 kpc from the Galactic centre. Several previously discovered stellar overdensities are recovered in the residual stellar density map, showing the reliability of mwfitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict Milky Way star counts for surveys such as the LSST and Euclid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Reliable photometric membership (RPM) of galaxies in clusters – I. A machine learning method and its performance in the local universe

ABSTRACT We introduce a new method to determine galaxy cluster membership based solely on photometric properties. We adopt a machine learning approach to recover a cluster membership probability from galaxy photometric parameters and finally derive a membership classification. After testing several machine learning techniques (such as stochastic gradient boosting, model averaged neural network and k-nearest neighbours), we found the support vector machine algorithm to perform better when applied to our data. Our training and validation data are from the Sloan Digital Sky Survey main sample. Hence, to be complete to $M_r^* + 3$, we limit our work to 30 clusters with $z$phot-cl ≤ 0.045. Masses (M200) are larger than $\sim 0.6\times 10^{14} \, \mathrm{M}_{\odot }$ (most above $3\times 10^{14} \, \mathrm{M}_{\odot }$). Our results are derived taking in account all galaxies in the line of sight of each cluster, with no photometric redshift cuts or background corrections. Our method is non-parametric, making no assumptions on the number density or luminosity profiles of galaxies in clusters. Our approach delivers extremely accurate results (completeness, C $\sim 92{\rm{ per\ cent}}$ and purity, P $\sim 87{\rm{ per\ cent}}$) within R200, so that we named our code reliable photometric membership. We discuss possible dependencies on magnitude, colour, and cluster mass. Finally, we present some applications of our method, stressing its impact to galaxy evolution and cosmological studies based on future large-scale surveys, such as eROSITA, EUCLID, and LSST.

Lopes, Paulo A. A.↗

DOE OSTI2020

Cross-correlating Planck with VST ATLAS LRGs: a new test for the ISW effect in the Southern hemisphere

ABSTRACT The integrated Sachs–Wolfe (ISW) effect probes the late-time expansion history of the Universe, offering direct constraints on dark energy. Here, we present our measurements of the ISW signal at redshifts of $\bar{z}=0.35$, 0.55, and 0.68, using the cross-correlation of the Planck cosmic microwave background temperature map with ∼0.5 million luminous red galaxies (LRGs) selected from the VST ATLAS survey. We then combine these with previous measurements based on WMAP and similar SDSS LRG samples, providing a total sample of ∼2.1 million LRGs covering ∼12 000 deg2 of sky. At $\bar{z}=0.35$ and $\bar{z}=0.55$, we detect the ISW signal at 1.2σ and 2.3σ (or 2.6σ combined), in agreement with the predictions of lambda cold dark matter (ΛCDM). We verify these results by repeating the measurements using the BOSS LOWZ and CMASS, spectroscopically confirmed LRG samples. We also detect the ISW effect in three magnitude limited ATLAS + SDSS galaxy samples extending to z ≈ 0.4 at ∼2σ per sample. However, we do not detect the ISW signal at $\bar{z}=0.68$ when combining the ATLAS and SDSS results. Further tests using spectroscopically confirmed eBOSS LRGs at this redshift remain inconclusive due to the current low sky coverage of the survey. If the ISW signal is shown to be redshift dependent in a manner inconsistent with the predictions of ΛCDM, it could open the door to alternative theories such as modified gravity. It is therefore important to repeat the high-redshift ISW measurement using the completed eBOSS sample, as well as deeper upcoming surveys such as DESI and LSST.

Ansarinejad, Behzad↗

DOE OSTI2020

First cosmology results using type Ia supernovae from the Dark Energy Survey: the effect of host galaxy properties on supernova luminosity

ABSTRACT We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. For the DES-SN sample, when considering a 5D (z, x1, c, α, β) bias correction, we find evidence of a Hubble residual ‘mass step’, where SNe Ia in high-mass galaxies (>1010M⊙) are intrinsically more luminous (after correction) than their low-mass counterparts by $\gamma =0.040\pm 0.019$ mag. This value is larger by 0.031 mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $\gamma =0.066\pm 0.020$ mag. The 1D−5D γ difference for DES-SN is $0.026\pm 0.009$ mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the x1 component of the 5D distance-bias correction. Including an intrinsic correlation between the observed properties of SNe Ia, stretch and colour, and stellar mass in simulated SN Ia samples, we show that a 5D fit recovers γ with −9 mmag bias compared to a +2 mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modelling correlations between galaxy properties and SN is necessary to ensure unbiased precision estimates of the dark energy equation of state as we enter the era of LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Predicting bulge to total luminosity ratio of galaxies using deep learning

ABSTRACT We present a deep learning model to predict the r-band bulge-to-total luminosity ratio (B/T) of nearby galaxies using their multiband JPEG images alone. Our Convolutional Neural Network (CNN) based regression model is trained on a large sample of galaxies with reliable decomposition into the bulge and disc components. The existing approaches to estimate the B/T ratio use galaxy light-profile modelling to find the best fit. This method is computationally expensive, prohibitively so for large samples of galaxies, and requires a significant amount of human intervention. Machine learning models have the potential to overcome these shortcomings. In our CNN model, for a test set of 20 000 galaxies, 85.7 per cent of the predicted B/T values have absolute error (AE) less than 0.1. We see further improvement to 87.5 per cent if, while testing, we only consider brighter galaxies (with r-band apparent magnitude <17) with no bright neighbours. Our model estimates the B/T ratio for the 20 000 test galaxies in less than a minute. This is a significant improvement in inference time from the conventional fitting pipelines, which manage around 2–3 estimates per minute. Thus, the proposed machine learning approach could potentially save a tremendous amount of time, effort, and computational resources while predicting B/T reliably, particularly in the era of next-generation sky surveys such as the Legacy Survey of Space and Time (LSST) and the Euclid sky survey which will produce extremely large samples of galaxies.

Grover, Harsh (ORCID:0000000321338142)↗

DOE OSTI2021

Optimizing a magnitude-limited spectroscopic training sample for photometric classification of supernovae

ABSTRACT In preparation for photometric classification of transients from the Legacy Survey of Space and Time (LSST) we run tests with different training data sets. Using estimates of the depth to which the 4-m Multi-Object Spectroscopic Telescope (4MOST) Time Domain Extragalactic Survey (TiDES) can classify transients, we simulate a magnitude-limited sample reaching rAB ≈ 22.5 mag. We run our simulations with the software snmachine, a photometric classification pipeline using machine learning. The machine-learning algorithms struggle to classify supernovae when the training sample is magnitude limited, in contrast to representative training samples. Classification performance noticeably improves when we combine the magnitude-limited training sample with a simulated realistic sample of faint high-redshift supernovae observed from larger spectroscopic facilities; the algorithms’ range of average area under receiver operator characteristic curve (AUC) scores over 10 runs increases from 0.547–0.628 to 0.946–0.969 and purity of the classified sample reaches 95 per cent in all runs for two of the four algorithms. By creating new, artificial light curves using the augmentation software avocado, we achieve a purity in our classified sample of 95 per cent in all 10 runs performed for all machine-learning algorithms considered. We also reach a highest average AUC score of 0.986 with the artificial neural network algorithm. Having ‘true’ faint supernovae to complement our magnitude-limited sample is a crucial requirement in optimization of a 4MOST spectroscopic sample. However, our results are a proof of concept that augmentation is also necessary to achieve the best classification results.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Modelling type 1 quasar colours in the era of Rubin and Euclid

We construct a parametric SED model which is able to reproduce the average observed SDSS–UKIDSS–WISE quasar colours to within one-tenth of a magnitude across a wide range of redshift (0 < z < 5) and luminosity (–22 > M i > –29). This model is shown to provide accurate predictions for the colours of known quasars which are less luminous than those used to calibrate the model parameters, and also those at higher redshifts z > 5. Using a single parameter, the model encapsulates an up-to-date understanding of the intra-population variance in the rest-frame ultraviolet and optical emission lines of luminous quasars. At fixed redshift, there are systematic changes in the average quasar colours with apparent i-band magnitude, which we find to be well explained by the contribution from the host galaxy and our parametrization of the emission-line properties. By including redshift as an additional free parameter, the model could be used to provide photometric redshifts for individual objects. For the population as a whole we find that the average emission line and host-galaxy contributions can be well described by simple functions of luminosity which account for the observed changes in the average quasar colours across 18.1 < i AB < 21.5. Finally, we use these trends to provide predictions for quasar colours at the luminosities and redshifts which will be probed by the Rubin Observatory LSST and ESA-Euclid wide survey. The model code is applicable to a wide range of upcoming photometric and spectroscopic surveys, and is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗