Source author record

Michele Liguori

Michele Liguori appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Clustering of Gravitational Wave and Supernovae events: a multitracer analysis in Luminosity Distance Space

We study the clustering of Gravitational Wave (GW) merger events and Supernovae IA (SN), as cosmic tracers in Luminosity Distance Space. We modify the publicly available CAMB code to numerically evaluate auto- and cross- power spectra for the different sources, including Luminosity Distance Space distortion effects generated by peculiar velocities and lensing convergence. We perform a multitracer Fisher analysis to forecast expected constraints on cosmological and GW bias coefficients, using outputs from hydrodinamical N-body simulations to determine the bias fiducial model and considering future observations from the Vera Rubin Observatory and Einstein Telescope (ET), both single and in a 3 detector network configuration. We find that adding SN to the GW merger dataset considerably improves the forecast, mostly by breaking significant parameter degeneracies, with final constraints comparable to those obtainable from a Euclid-like survey. GW merger bias is forecasted to be detectable with good significance even in the single ET case.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

preprint2021arXiv

Cross-correlating Astrophysical and Cosmological Gravitational Wave Backgrounds with the Cosmic Microwave Background

General Relativity provides us with an extremely powerful tool to extract at the same time astrophysical and cosmological information from the Stochastic Gravitational Wave Backgrounds (SGWBs): the cross-correlation with other cosmological tracers, since their anisotropies share a common origin and the same perturbed geodesics. In this letter we explore the cross-correlation of the cosmological and astrophysical SGWBs with Cosmic Microwave Background (CMB) anisotropies, showing that future GW detectors, such as LISA or BBO, have the ability to measure such cross-correlation signals. We also present, as a new tool in this context, constrained realization maps of the SGWBs extracted from the high-resolution CMB {\it Planck} maps. This technique allows, in the low-noise regime, to faithfully reconstruct the expected SGWB map by starting from CMB measurements.

preprint2016arXiv

Angular dependence of primordial trispectra and CMB spectral distortions

Under the presence of anisotropic sources in the inflationary era, the trispectrum of the primordial curvature perturbation has a very specific angular dependence between each wavevector that is distinguishable from the one encountered when only scalar fields are present, characterized by an angular dependence described by Legendre polynomials. We examine the imprints left by curvature trispectra on the $TTμ$ bispectrum, generated by the correlation between temperature anisotropies (T) and chemical potential spectral distortions ($μ$) of the Cosmic Microwave Background (CMB). Due to the angular dependence of the primordial signal, the corresponding $TTμ$ bispectrum strongly differs in shape from $TTμ$ sourced by the usual $g_{\rm NL}$ or $τ_{\rm NL}$ local trispectra, enabling us to obtain an unbiased estimation. From a Fisher matrix analysis, we find that, in a cosmic-variance-limited (CVL) survey of $TTμ$, a minimum detectable value of the quadrupolar Legendre coefficient is $d_2 \sim 0.01$, which is 4 orders of magnitude better than the best value attainable from the $TTTT$ CMB trispectrum. In the case of an anisotropic inflationary model with a $f(ϕ)F^2$ interaction (coupling the inflaton field $ϕ$ with a vector kinetic term $F^2$), the size of the curvature trispectrum is related to that of quadrupolar power spectrum asymmetry, $g_*$. In this case, a CVL measurement of $TTμ$ makes it possible to measure $g_*$ down to $10^{-3}$.

preprint2016arXiv

Gravitational waves from inflation

The production of a stochastic background of gravitational waves is a fundamental prediction of any cosmological inflationary model. The features of such a signal encode unique information about the physics of the Early Universe and beyond, thus representing an exciting, powerful window on the origin and evolution of the Universe. We review the main mechanisms of gravitational-wave production, ranging from quantum fluctuations of the gravitational field to other mechanisms that can take place during or after inflation. These include e.g. gravitational waves generated as a consequence of extra particle production during inflation, or during the (p)reheating phase. Gravitational waves produced in inflation scenarios based on modified gravity theories and second-order gravitational waves are also considered. For each analyzed case, the expected power-spectrum is given. We discuss the discriminating power among different models, associated with the validity/violation of the standard consistency relation between tensor-to-scalar ratio $r$ and tensor spectral index $n_{\rm T}$. In light of the prospects for (directly/indirectly) detecting primordial gravitational waves, we give the expected present-day gravitational radiation spectral energy-density, highlighting the main characteristics imprinted by the cosmic thermal history, and we outline the signatures left by gravitational waves on the Cosmic Microwave Background and some imprints in the Large-Scale Structure of the Universe. Finally, current bounds and prospects of detection for inflationary gravitational waves are summarized.

preprint2016arXiv

Primordial trispectra and CMB spectral distortions

We study the $TTμ$ bispectrum, generated by correlations between Cosmic Microwave Background temperature (T) anisotropies and chemical potential ($μ$) distortions, and we analyze its dependence on primordial local trispectrum parameters $g_{\rm NL}$ and $τ_{\rm NL}$. We cross-check our results by comparing the full bispectrum calculation with the expectations from a general physical argument, based on predicting the shape of $μ$-T correlations from the couplings between short and long perturbation modes induced by primordial non-Gaussianity. We show that $both$ $g_{\rm NL}$ and $τ_{\rm NL}$-parts of the primordial trispectrum source a non-vanishing $TTμ$ signal, contrary to the $μμ$ auto-correlation function, which is sensitive only to the $τ_{\rm NL}$-component. A simple Fisher matrix-based forecast shows that a futuristic, cosmic-variance dominated experiment could in principle detect $g_{\rm NL} \sim 0.4$ and $τ_{\rm NL} \sim 40$ using $TTμ$.

preprint2016arXiv

Recent results and perspectives on cosmology and fundamental physics from microwave surveys

Recent cosmic microwave background data in temperature and polarization have reached high precision in estimating all the parameters that describe the current so-called standard cosmological model. Recent results about the integrated Sachs-Wolfe effect from cosmic microwave background anisotropies, galaxy surveys, and their cross-correlations are presented. Looking at fine signatures in the cosmic microwave background, such as the lack of power at low multipoles, the primordial power spectrum and the bounds on non-Gaussianities, complemented by galaxy surveys, we discuss inflationary physics and the generation of primordial perturbations in the early Universe. Three important topics in particle physics, the bounds on neutrinos masses and parameters, on thermal axion mass and on the neutron lifetime derived from cosmological data are reviewed, with attention to the comparison with laboratory experiment results. Recent results from cosmic polarization rotation analyses aimed at testing the Einstein equivalence principle are presented. Finally, we discuss the perspectives of next radio facilities for the improvement of the analysis of future cosmic microwave background spectral distortion experiments.

preprint2016arXiv

Science with the space-based interferometer LISA. IV: Probing inflation with gravitational waves

We investigate the potential for the LISA space-based interferometer to detect the stochastic gravitational wave background produced from different mechanisms during inflation. Focusing on well-motivated scenarios, we study the resulting contributions from particle production during inflation, inflationary spectator fields with varying speed of sound, effective field theories of inflation with specific patterns of symmetry breaking and models leading to the formation of primordial black holes. The projected sensitivities of LISA are used in a model-independent way for various detector designs and configurations. We demonstrate that LISA is able to probe these well-motivated inflationary scenarios beyond the irreducible vacuum tensor modes expected from any inflationary background.

preprint2015arXiv

Future Constraints on Angle-Dependent Non-Gaussianity from Large Radio Surveys

We investigate how well future large-scale radio surveys could measure different shapes of primordial non-Gaussianity; in particular we focus on angle-dependent non-Gaussianity arising from primordial anisotropic sources, whose bispectrum has an angle dependence between the three wavevectors that is characterized by Legendre polynomials $\mathcal{P}_L$ and expansion coefficients $c_L$. We provide forecasts for measurements of galaxy power spectrum, finding that Large-Scale Structure (LSS) data could allow measurements of primordial non-Gaussianity competitive or improving upon current constraints set by CMB experiments, for all the shapes considered. We argue that the best constraints will come from the possibility to assign redshift information to radio galaxy surveys, and investigate a few possible scenarios for the EMU and SKA surveys. A realistic (futuristic) modeling could provide constraints of $f_{\rm NL}^{\rm loc} \approx 1 (0.5)$ for the local shape, $f_{\rm NL}$ of $\mathcal{O}(10) (\mathcal{O}(1))$ for the orthogonal, equilateral and folded shapes, and $c_{L=1} \approx 80 (2)$, $c_{L=2} \approx 400 (10)$ for angle-dependent non-Gaussianity. The more futuristic forecasts show the potential of LSS analyses to considerably improve current constraints on non-Gaussianity, and so on models of the primordial Universe. Finally, we find the minimum requirements that would be needed to reach $σ(c_{L=1})=10$, which can be considered as a typical (lower) value predicted by some (inflationary) models.

preprint2015arXiv

Measuring primordial anisotropic correlators with CMB spectral distortions

We show that inflationary models with broken rotational invariance generate testable off-diagonal signatures in the correlation between the $μ$-type distortion and temperature fluctuations of the cosmic microwave background. More precisely, scenarios with a quadrupolar bispectrum asymmetry, usually generated by fluctuations of primordial vector fields, produce a nonvanishing $μ$-$T$ correlation when $|\ell_1-\ell_2|=2$. Since spectral distortions are sensitive to primordial fluctuations up to very small scales, a cosmic variance limited spectral distortion experiment can detect such effects with a high signal-to-noise ratio.

preprint2015arXiv

Observed parity-odd CMB temperature bispectrum

Parity-odd non-Gaussianities create a variety of temperature bispectra in the cosmic microwave background (CMB), defined in the domain: $\ell_1 + \ell_2 + \ell_3 = {\rm odd}$. These models are yet unconstrained in the literature, that so far focused exclusively on the more common parity-even scenarios. In this work, we provide the first experimental constraints on parity-odd bispectrum signals in WMAP 9-year temperature data, using a separable modal parity-odd estimator. Comparing theoretical bispectrum templates to the observed bispectrum, we place constraints on the so-called nonlineality parameters of parity-odd tensor non-Gaussianities predicted by several Early Universe models. Our technique also generates a model-independent, smoothed reconstruction of the bispectrum of the data for parity-odd configurations.

preprint2014arXiv

General parity-odd CMB bispectrum estimation

We develop a methodology for estimating parity-odd bispectra in the cosmic microwave background (CMB). This is achieved through the extension of the original separable modal methodology to parity-odd bispectrum domains ($\ell_1 + \ell_2 + \ell_3 = {\rm odd}$). Through numerical tests of the parity-odd modal decomposition with some theoretical bispectrum templates, we verify that the parity-odd modal methodology can successfully reproduce the CMB bispectrum, without numerical instabilities. We also present simulated non-Gaussian maps produced by modal-decomposed parity-odd bispectra, and show the consistency with the exact results. Our new methodology is applicable to all types of parity-odd temperature and polarization bispectra.

preprint2014arXiv

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): An Extended White Paper

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

preprint2013arXiv

Future CMB ISW-Lensing bispectrum constraints on modified gravity in the Parameterized Post-Friedmann formalism

We forecast the constraints on both Hu-Sawicki model and Bertschinger-Zukin model of modified gravity within the Parameterized Post-Friedmann (PPF) formalism for the Planck satellite experiment by performing the joint analysis of ISW-Lensing bispectrum and CMB power spectrum. We find that, even considering the temperature-temperature mode of CMB power spectrum only, Planck data are expected to reduce the error bars on the modified gravity parameter $B_0$ (related to the present value of Compton wavelength of the extra scalar degree of freedom) at least one order magnitude compared with WMAP. The spectrum-bispectrum joint analysis can further improve the results by a factor ranging from 1.14 to 5.32 depending on the specific modified gravity model. One of our main results is that the cross-correlation between ISW-Lensing bispectrum and power spectrum can be safely neglected when performing the joint analysis. For simplicity, we only investigate the likelihood of one parameter ({$B_0$}) and fix all other cosmological parameters to their best-fit values in WMAP7yr results.

preprint2013arXiv

Non-Gaussianity and CMB aberration and Doppler

The peculiar motion of an observer with respect to the CMB rest frame induces a deflection in the arrival direction of the observed photons (also known as CMB aberration) and a Doppler shift in the measured photon frequencies. As a consequence, aberration and Doppler effects induce non trivial correlations between the harmonic coefficients of the observed CMB temperature maps. In this paper we investigate whether these correlations generate a bias on Non-Gaussianity estimators $f_{NL}$. We perform this analysis simulating a large number of temperature maps with Planck-like resolution (lmax $= 2000$) as different realizations of the same cosmological fiducial model (WMAP7yr). We then add to these maps aberration and Doppler effects employing a modified version of the HEALPix code. We finally evaluate a generalization of the Komatsu, Spergel and Wandelt Non-Gaussianity estimator for all the simulated maps, both when peculiar velocity effects have been considered and when these phenomena have been neglected. Using the value $v/c=1.23 \times 10^{-3}$ for our peculiar velocity, we found that the aberration/Doppler induced Non-Gaussian signal is at most of about half of the cosmic variance $σ$ for $f_{NL}$ both in a full-sky and in a cut-sky experimental configuration, for local, equilateral and orthogonal estimators. We conclude therefore that when estimating $f_{NL}$ it is safe to ignore aberration and Doppler effects {\it if} the primordial map is already Gaussian. More work is necessary however to assess whether a map which contains Non-Gaussianity can be significantly distorted by a peculiar velocity.

preprint2013arXiv

Optimal bispectrum estimator and simulations of the the CMB Lensing-ISW non-Gaussian signal

In this paper we present the tools to optimally extract the Lensing-Integrated Sachs Wolfe (L-ISW) bispectrum signal from future CMB data. We implement two different methods to simulate the non-Gaussian CMB maps with the L-ISW signal: a non-perturbative method based on the FLINTS lensing code and the separable mode expansion method. We implement the Komatsu, Spergel and Wandelt (KSW) optimal estimator analysis for the Lensing-ISW bispectrum and we test it on the non-Gaussian simulations in the case of a realistic CMB experimental settings with an inhomogeneous sky coverage. We show that the estimator approaches the Cramer-Rao bound and that Wiener filtering the L-ISW simulations gives a slight improvement on the estimate of $f_{NL}^{L-ISW}$ of $\leq 10%$. For a realistic CMB experimental setting accounting for anisotropic noise and masked sky, we show that the linear term of the estimator is highly correlated to the cubic term and it is necessary to recover the signal and the optimal error bars. We also show that the L-ISW bispectrum, if not correctly accounted for, yields an underestimation of the $f_{NL}^{local}$ error bars of $\simeq 4%$. A joint analysis of the non-Gaussian shapes and/or L-ISW template subtraction is needed in order to recover unbiased results of the primordial non-Gaussian signal from ongoing and future CMB experiments.

preprint2013arXiv

Parametrized modified gravity constraints after Planck

We constrain $f(R)$ and chameleon-type modified gravity in the framework of the Berstchinger-Zukin parametrization using the recent released Planck data, including both CMB temperature power spectrum and lensing potential power spectrum. Some other external data sets are included, such as BAO measurements from the 6dFGS, SDSS DR7 and BOSS DR9 surveys, HST $H_0$ measurement and supernovae from Union2.1 compilation. We also use WMAP9yr data for consistency check and comparison. For $f(R)$ gravity, WMAP9yr results can only give quite a loose constraint on the modified gravity parameter $B_0$, which is related to the present value of the Compton wavelength of the extra scalar degree of freedom, $B_0<3.37$ at $95\% {\rm C.L.}$ We demonstrate that this constraint mainly comes from the late ISW effect. With only Planck CMB temperature power-spectrum data, we can improve the WMAP9yr result by a factor $3.7$ ($B_0<0.91$ at $95\% {\rm C.L.}$). If the Planck lensing potential power-spectrum data are also taken into account, the constraint can be further strenghtened by a factor $5.1$ ($B_0<0.18$ at $95\% {\rm C.L.}$). This major improvement mainly comes from the small-scale lensing signal. Furthermore, BAO, HST and supernovae data could slightly improve the $B_0$ bound ($B_0<0.12$ at $95\% {\rm C.L.}$).For the chameleon-type model, we find that the data set which we used cannot constrain the Compton wavelength $B_0$ and the potential index $s$ of chameleon field, but can give a tight constraint on the parameter $β_1=1.043^{+0.163}_{-0.104}$ at $95\% {\rm C.L.}$ ($β_1=1$ in general relativity), which accounts for the non-minimal coupling between the chameleon field and the matter component. In addition, we find that both modified gravity models we considered favor a relatively higher Hubble parameter than the concordance LCDM model in general relativity.

preprint2013arXiv

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): A White Paper on the Ultimate Polarimetric Spectro-Imaging of the Microwave and Far-Infrared Sky

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

preprint2012arXiv

On the linear term correction for needlets/wavelets non-Gaussianity estimators

We derive the linear correction term for needlet and wavelet estimators of the bispectrum and the non-linearity parameter fNL on cosmic microwave background radiation data. We show that on masked WMAP-like data with anisotropic noise, the error bars improve by 10-20% and almost reach the optimal error bars obtained with the KSW estimator (Komatsu et al 2005). In the limit of full-sky and isotropic noise, this term vanishes. We apply needlet and wavelet estimators to the WMAP 7-year data and obtain our best estimate fNL=37.5 \pm 21.8.

preprint2010arXiv

CMB lensing and primordial non-Gaussianity

We study the effects of gravitational lensing on the estimation of non-Gaussianity from the bispectrum of the cosmic microwave background (CMB) temperature anisotropies. We find that the effect of lensing on the bispectrum may qualitatively be described as a smoothing of the acoustic features analogous to the temperature power spectrum. In contrast to previous results, for a Planck-like experiment which is cosmic-variance limited to L=2000, we find that lensing causes no significant degradation of our ability to constrain the non-Gaussianity amplitude fNL for both local and equilateral configurations, provided that the biases due to the cross correlation between the lensing potential and the integrated-Sachs-Wolfe (ISW) contribution to the CMB temperature are adequately understood. With numerical simulations, we also verify that low-order Taylor approximations to the lensed bispectrum and ISW-lensing biases are accurate.

preprint2009arXiv

Constraining Running Non-Gaussianity

The primordial non-Gaussian parameter fNL has been shown to be scale-dependent in several models of inflation with a variable speed of sound. Starting from a simple ansatz for a scale-dependent amplitude of the primordial curvature bispectrum for two common phenomenological models of primordial non-Gaussianity, we perform a Fisher matrix analysis of the bispectra of the temperature and polarization of the Cosmic Microwave Background (CMB) radiation and derive the expected constraints on the parameter nNG that quantifies the running of fNL(k) for current and future CMB missions such as WMAP, Planck and CMBPol. We find that CMB information alone, in the event of a significant detection of the non-Gaussian component, corresponding to fNL = 50 for the local model and fNL = 100 for the equilateral model of non-Gaussianity, is able to determine nNG with a 1-sigma uncertainty of Delta nNG = 0.1 and Delta nNG = 0.3, respectively, for the Planck mission. In addition, we consider a Fisher matrix analysis of the galaxy power spectrum to determine the expected constraints on the running parameter nNG for the local model and of the galaxy bispectrum for the equilateral model from future photometric and spectroscopic surveys. We find that, in both cases, large-scale structure observations should achieve results comparable to or even better than those from the CMB, while showing some complementarity due to the different distribution of the non-Gaussian signal over the relevant range of scales. Finally, we compare our findings to the predictions on the amplitude and running of non-Gaussianity of DBI inflation, showing how the constraints on a scale-dependent fNL(k) translate into constraints on the parameter space of the theory.

preprint2007arXiv

Fast Estimator of Primordial Non-Gaussianity from Temperature and Polarization Anisotropies in the Cosmic Microwave Background II: Partial Sky Coverage and Inhomogeneous Noise

In our recent paper (Yadav et al. 2007) we described a fast cubic (bispectrum) estimator of the amplitude of primordial non-Gaussianity of local type, f_{NL}, from a combined analysis of the Cosmic Microwave Background (CMB) temperature and E-polarization observations. In this paper we generalize the estimator to deal with a partial sky coverage as well as inhomogeneous noise. Our generalized estimator is still computationally efficient, scaling as O(N^3/2) compared to the O(N^5/2) scaling of the brute force bispectrum calculation for sky maps with N pixels. Upcoming CMB experiments are expected to yield high-sensitivity temperature and E-polarization data. Our generalized estimator will allow us to optimally utilize the combined CMB temperature and E-polarization information from these realistic experiments, and to constrain primordial non-Gaussianity.

preprint2007arXiv

Temperature and Polarization CMB Maps from Primordial non-Gaussianities of the Local Type

The forthcoming Planck experiment will provide high sensitivity polarization measurements that will allow us to further tighten the f_NL bounds from the temperature data. Monte Carlo simulations of non-Gaussian CMB maps have been used as a fundamental tool to characterize non-Gaussian signatures in the data, as they allow us to calibrate any statistical estimators and understand the effect of systematics, foregrounds and other contaminants. We describe an algorithm to generate high-angular resolution simulations of non-Gaussian CMB maps in temperature and polarization. We consider non-Gaussianities of the local type, for which the level of non-Gaussianity is defined by the dimensionless parameter, f_NL. We then apply the temperature and polarization fast cubic statistics recently developed by Yadav et al. to a set of non-Gaussian temperature and polarization simulations. We compare our results to theoretical expectations based on a Fisher matrix analysis, test the unbiasedness of the estimator, and study the dependence of the error bars on f_NL. All our results are in very good agreement with theoretical predictions, thus confirming the reliability of both the simulation algorithm and the fast cubic temperature and polarization estimator.