QML-FAST -- A Fast Code for low-$\ell$ Tomographic Maximum Likelihood Power Spectrum Estimation

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Kvasiuk, Yurii, Lai, Anderson, Münchmeyer, Moritz, Smith, Kendrick M.
Formato: Preprint
Publicado: 2025
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866910177712668672
author Kvasiuk, Yurii
Lai, Anderson
Münchmeyer, Moritz
Smith, Kendrick M.
author_facet Kvasiuk, Yurii
Lai, Anderson
Münchmeyer, Moritz
Smith, Kendrick M.
contents We present a novel implementation for the quadratic maximum likelihood (QML) power spectrum estimator for multiple correlated scalar fields on the sphere. Our estimator supports arbitrary binning in redshift and multipoles $\ell$ and includes cross-correlations of redshift bins. It implements a fully optimal analysis with a pixel-wise covariance model. We implement a number of optimizations which make the estimator and associated covariance matrix computationally tractable for a low-$\ell$ analysis, suitable for example for kSZ velocity reconstruction or primordial non-Gaussianity from scale-dependent bias analyses. We validate our estimator extensively on simulations and compare its features and precision with the common pseudo-$C_\ell$ method, showing significant gains at large scales. We make our code publicly available. In a companion paper, we apply the estimator to kSZ velocity reconstruction using data from ACT and DESI Legacy Survey and construct full set of QML estimators on 40 correlated fields up to $N_{\text{side}}= 32$ in timescale of an hour on a single 24-core CPU requiring $<256\ \mathrm{Gb}$ RAM, demonstrating the performance of the code.
format Preprint
id arxiv_https___arxiv_org_abs_2510_05215
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle QML-FAST -- A Fast Code for low-$\ell$ Tomographic Maximum Likelihood Power Spectrum Estimation
Kvasiuk, Yurii
Lai, Anderson
Münchmeyer, Moritz
Smith, Kendrick M.
Cosmology and Nongalactic Astrophysics
We present a novel implementation for the quadratic maximum likelihood (QML) power spectrum estimator for multiple correlated scalar fields on the sphere. Our estimator supports arbitrary binning in redshift and multipoles $\ell$ and includes cross-correlations of redshift bins. It implements a fully optimal analysis with a pixel-wise covariance model. We implement a number of optimizations which make the estimator and associated covariance matrix computationally tractable for a low-$\ell$ analysis, suitable for example for kSZ velocity reconstruction or primordial non-Gaussianity from scale-dependent bias analyses. We validate our estimator extensively on simulations and compare its features and precision with the common pseudo-$C_\ell$ method, showing significant gains at large scales. We make our code publicly available. In a companion paper, we apply the estimator to kSZ velocity reconstruction using data from ACT and DESI Legacy Survey and construct full set of QML estimators on 40 correlated fields up to $N_{\text{side}}= 32$ in timescale of an hour on a single 24-core CPU requiring $<256\ \mathrm{Gb}$ RAM, demonstrating the performance of the code.
title QML-FAST -- A Fast Code for low-$\ell$ Tomographic Maximum Likelihood Power Spectrum Estimation
topic Cosmology and Nongalactic Astrophysics
url https://arxiv.org/abs/2510.05215