Source author record

Cullan Howlett

Cullan Howlett 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

17works
5topics
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

17 published item(s)

preprint2026arXiv

Effective Theories of Redshift-Space Galaxy Peculiar Velocities

We present predictions for redshift-space peculiar velocity statistics in the Lagrangian and Eulerian formulations of the effective field theory (EFT) of large-scale structure. We compute 2-point pairwise velocity statistics up to the second moment at next-to-leading (1-loop) order, showing that they can be modeled together with redshift-space galaxy densities with a consistent set of EFT coefficients. We show that peculiar velocity statistics have a distinct dependence on long-wavelength bulk flows that necessitates a variation on the usual infrared (IR) resummation procedure used to model baryon acoustic oscillations (BAO) in galaxy clustering. This can be implemented recursively in powers of the velocity in both the Lagrangian and Eulerian frameworks. We validate our analytic calculations against fully nonlinear N-body simulations, demonstrating that they can be used to recover the growth rate at better than percent level precision, well beyond the statistical requirements of upcoming peculiar velocity surveys and measurements of the kinetic Sunyaev-Zeldovich (kSZ) effect. As part of this work, we release $\href{https://github.com/sfschen/velocisaurus}{\texttt{velocisaurus}}$, a fast $\texttt{Python}$ code for computing EFT predictions of peculiar velocity statistics.

preprint2026arXiv

Limits of self-interacting neutrinos from the BAO and CMB phase shift

Neutrinos with Standard Model interactions free-stream in the early Universe, leaving a distinct phase shift in the pattern of baryon acoustic oscillations (BAO). When isolated, this phase shift allows one to robustly infer the presence of the cosmic neutrino background in BAO and cosmic microwave background (CMB) data independently of other cosmological parameters. While in the context of the Standard Model, this phase shift follows a known scale-dependent relation, new physics in the cosmic neutrino background could alter the overall shape of this feature. In this paper, we discuss how changes in the neutrino phase shift could be used to constrain self-interactions among neutrinos. We produce simple models for this phase-shift assuming universal self-interactions, and use these in order to understand what constraining power is available for the strength of such interactions in BAO and CMB data. We find that, although challenging, it may be possible to use a detection of the phase to put a more robust limit on the strength of the self-interaction, $G_{\mathrm{eff}}$, which at present suffers from bimodality in cosmological constraints. Our forecast analysis reveals that BAO data alone will not provide the precision needed to tightly constrain self-interactions; however, the combined analysis of the phase shift signature in both CMB and BAO can potentially provide a way to detect the impact of new neutrino interactions. Our results could be extended upon for models with non-universal interactions.

preprint2022arXiv

The Sloan Digital Sky Survey Peculiar Velocity Catalogue

We present a new catalogue of distances and peculiar velocities (PVs) of $34,059$ early-type galaxies derived from Fundamental Plane (FP) measurements using data from the Sloan Digital Sky Survey (SDSS). This $7016\,\mathrm{deg}^{2}$ homogeneous sample comprises the largest set of peculiar velocities produced to date and extends the reach of PV surveys up to a redshift limit of $z=0.1$. Our SDSS-based FP distance measurements have a mean uncertainty of 23%. Alongside the data, we produce an ensemble of 2,048 mock galaxy catalogues that reproduce the data selection function, and are used to validate our fitting pipelines and check for systematic errors. We uncover a significant trend between group richness and mean surface brightness within the sample, which may hint at an environmental dependence within the FP or the presence of unresolved systematics, and can result in biased peculiar velocities. This is removed using multiple FP fits as function of group richness, a procedure made tractable through a new analytic derivation for the integral of a 3D Gaussian over non-trivial limits. Our catalogue is calibrated to the zero-point of the CosmicFlows-III sample with an uncertainty of $0.004$ dex (not including cosmic variance or the error within CosmicFlows-III itself), which is validated using independent cross-checks with the predicted zero-point from the 2M++ reconstruction of our local velocity field. Finally, as an example of what is possible with our new catalogue, we obtain preliminary bulk flow measurements up to a depth of $135\,h^{-1}\mathrm{Mpc}$. We find a slightly larger-than-expected bulk flow at high redshift, although this could be caused by the presence of the Shapley supercluster which lies outside the SDSS PV footprint.

preprint2021arXiv

CosmoBit: A GAMBIT module for computing cosmological observables and likelihoods

We introduce $\sf{CosmoBit}$, a module within the open-source $\sf{GAMBIT}$ software framework for exploring connections between cosmology and particle physics with joint global fits. $\sf{CosmoBit}$ provides a flexible framework for studying various scenarios beyond $Λ$CDM, such as models of inflation, modifications of the effective number of relativistic degrees of freedom, exotic energy injection from annihilating or decaying dark matter, and variations of the properties of elementary particles such as neutrino masses and the lifetime of the neutron. Many observables and likelihoods in $\sf{CosmoBit}$ are computed via interfaces to $\sf{AlterBBN}$, $\sf{CLASS}$, $\sf{DarkAges}$, $\sf{MontePython}$, $\sf{MultiModeCode}$, and $\sf{plc}$. This makes it possible to apply a wide range of constraints from large-scale structure, Type Ia supernovae, Big Bang Nucleosynthesis and the cosmic microwave background. Parameter scans can be performed using the many different statistical sampling algorithms available within the $\sf{GAMBIT}$ framework, and results can be combined with calculations from other $\sf{GAMBIT}$ modules focused on particle physics and dark matter. We include extensive validation plots and a first application to scenarios with non-standard relativistic degrees of freedom and neutrino temperature, showing that the corresponding constraint on the sum of neutrino masses is much weaker than in the standard scenario.

preprint2020arXiv

Barry and the BAO Model Comparison

We compare the performance of four state-of-the-art models for extracting isotropic measurements of the Baryon Acoustic Oscillation (BAO) scale. To do this, we created a new, public, modular code Barry, which contains datasets, model fitting tools, and model implementations incorporating different descriptions of non-linear physics and algorithms for isolating the BAO feature. These are then evaluated for bias, correlation, and fitting strength using mock power spectra and correlation functions developed for the Sloan Digital Sky Survey Data Release 12. Our main findings are as follows: 1) All of the models can recover unbiased constraints when fit to the pre- and post-reconstruction simulations. 2) Models that provide physical descriptions of the damping of the BAO feature (using e.g., standard perturbation or effective-field theory arguments) report smaller errors on average, although the distribution of mock $χ^{2}$ values indicates these are underestimated. 3) Allowing the BAO damping scale to vary can provide tighter constraints for some mocks, but is an artificial improvement that only arises when noise randomly sharpens the BAO peak. 4) Unlike recent claims in the literature when utilising a BAO Extractor technique, we find no improvement in the accuracy of the recovered BAO scale. 5) We implement a procedure for combining all models into a single consensus result that improves over the standard method without obviously underestimating the uncertainties. Overall, Barry provides an excellent framework for performing the cosmological analyses for upcoming surveys, and for rapidly testing and validating new models.

preprint2020arXiv

Multi-wavelength consensus of large-scale linear bias

We model the large-scale linear galaxy bias $b_g(x,z)$ as a function of redshift $z$ and observed absolute magnitude threshold $x$ for broadband continuum emission from the far infrared to ultra-violet, as well as for prominent emission lines, such as the H$α$, H$β$, Lya and [OII] lines. The modelling relies on the semi-analytic galaxy formation model GALFORM, run on the state-of-the-art $N$-body simulation SURFS with the Planck 2015 cosmology. We find that both the differential bias at observed absolute magnitude $x$ and the cumulative bias for magnitudes brighter than $x$ can be fitted with a five-parameter model: $b_g(x,z)=a + b(1+z)^e(1 + \exp{[(x-c)d]})$. We also find that the bias for the continuum bands follows a very similar form regardless of wavelength due to the mixing of star-forming and quiescent galaxies in a magnitude limited survey. Differences in bias only become apparent when an additional colour separation is included, which suggest extensions to this work could look at different colours at fixed magnitude limits. We test our fitting formula against observations, finding reasonable agreement with some measurements within $1σ$ statistical uncertainties, and highlighting areas of improvement. We provide the fitting parameters for various continuum bands, emission lines and intrinsic galaxy properties, enabling a quick estimation of the linear bias in any typical survey of large-scale structure.

preprint2020arXiv

Standard Siren Speeds: Improving velocities in gravitational-wave measurements of $H_{0}$

We re-analyze data from the gravitational-wave event GW170817 and its host galaxy NGC4993 to demonstrate the importance of accurate total and peculiar velocities when measuring the Hubble constant using this nearby Standard Siren. We show that a number of reasonable choices can be made to estimate the velocities for this event, but that systematic differences remain between these measurements depending on the data used. This leads to significant changes in the Hubble constant inferred from GW170817. We present Bayesian Model Averaging as one way to account for these differences, and obtain $H_{0}=66.8^{+13.4}_{-9.2}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$. Adding additional information on the viewing angle from high resolution imaging of the radio counterpart refines this to $H_{0}=64.8^{+7.3}_{-7.2}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$. During this analysis we also present an alternative Bayesian model for the posterior on $H_{0}$ from Standard Sirens that works more closely with observed quantities from redshift and peculiar velocity surveys. Our results more accurately capture the true uncertainty on the total and peculiar velocities of NGC4993 and show that exploring how well different datasets characterize galaxy groups and the velocity field in the local Universe could improve this measurement further. These considerations impact any low-redshift distance measurement, and the improvements we suggest here can also be applied to standard candles like type Ia supernovae. GW170817 is particularly sensitive to peculiar velocity uncertainties because it is so close. For future standard siren measurements the importance of this error will decrease as (a) we will measure more distant standard sirens, and (b) the random direction of peculiar velocities will average out with more detections.

preprint2016arXiv

Cosmological Forecasts for Combined and Next Generation Peculiar Velocity Surveys

Peculiar velocity surveys present a very promising route to measuring the growth rate of large-scale structure and its scale dependence. However, individual peculiar velocity surveys suffer from large statistical errors due to the intrinsic scatter in the relations used to infer a galaxy's true distance. In this context we use a Fisher Matrix formalism to investigate the statistical benefits of combining multiple peculiar velocity surveys. We find that for all cases we consider there is a marked improvement on constraints on the linear growth rate $fσ_{8}$. For example, the constraining power of only a few peculiar velocity measurements is such that the addition of the 2MASS Tully-Fisher survey (containing only $\sim2,000$ galaxies) to the full redshift and peculiar velocity samples of the 6-degree Field Galaxy Survey (containing $\sim 110,000$ redshifts and $\sim 9,000$ velocities) can improve growth rate constraints by $\sim20\%$. Furthermore, the combination of the future TAIPAN and WALLABY+WNSHS surveys has the potential to reach a $\sim3\%$ error on $fσ_{8}$, which will place tight limits on possible extensions to General Relativity. We then turn to look at potential systematics in growth rate measurements that can arise due to incorrect calibration of the peculiar velocity zero-point and from scale-dependent spatial and velocity bias. For next generation surveys, we find that neglecting velocity bias in particular has the potential to bias constraints on the growth rate by over $5σ$, but that an offset in the zero-point has negligible impact on the velocity power spectrum.

preprint2016arXiv

Cosmology with Peculiar Velocities: Observational Effects

In this paper we investigate how observational effects could possibly bias cosmological inferences from peculiar velocity measurements. Specifically, we look at how bulk flow measurements are compared with theoretical predictions. Usually bulk flow calculations try to approximate the flow that would occur in a sphere around the observer. Using the Horizon Run 2 simulation we show that the traditional methods for bulk flow estimation can overestimate the magnitude of the bulk flow for two reasons: when the survey geometry is not spherical (the data do not cover the whole sky), and when the observations undersample the velocity distributions. Our results may explain why several bulk flow measurements found bulk flow velocities that seem larger than those expected in standard ΛCDM cosmologies. We recommend a different approach when comparing bulk flows to cosmological models, in which the theoretical prediction for each bulk flow measurement is calculated specifically for the geometry and sampling rate of that survey. This means that bulk flow values will not be comparable between surveys, but instead they are comparable with cosmological models, which is the more important measure.

preprint2015arXiv

Cosmological implications of baryon acoustic oscillation (BAO) measurements

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $Ω_m=0.301 \pm 0.008$ and curvature $Ω_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $χ^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

preprint2015arXiv

L-PICOLA: A parallel code for fast dark matter simulation

Robust measurements based on current large-scale structure surveys require precise knowledge of statistical and systematic errors. This can be obtained from large numbers of realistic mock galaxy catalogues that mimic the observed distribution of galaxies within the survey volume. To this end we present a fast, distributed-memory, planar-parallel code, L-PICOLA, which can be used to generate and evolve a set of initial conditions into a dark matter field much faster than a full non-linear N-Body simulation. Additionally, L-PICOLA has the ability to include primordial non-Gaussianity in the simulation and simulate the past lightcone at run-time, with optional replication of the simulation volume. Through comparisons to fully non-linear N-Body simulations we find that our code can reproduce the $z=0$ power spectrum and reduced bispectrum of dark matter to within 2% and 5% respectively on all scales of interest to measurements of Baryon Acoustic Oscillations and Redshift Space Distortions, but 3 orders of magnitude faster. The accuracy, speed and scalability of this code, alongside the additional features we have implemented, make it extremely useful for both current and next generation large-scale structure surveys. L-PICOLA is publicly available at https://cullanhowlett.github.io/l-picola

preprint2015arXiv

Reconstruction in Fourier space

We present a fast iterative FFT-based reconstruction algorithm that allows for non- parallel redshift-space distortions (RSD). We test our algorithm on both N-body dark matter simulations and mock distributions of galaxies designed to replicate galaxy survey conditions. We compare solenoidal and irrotational components of the redshift distortion and show that an approximation of this distortion leads to a better estimate of the real-space potential (and therefore faster convergence) than ignoring the RSD when estimating the displacement field. Our iterative reconstruction scheme converges in two iterations for the mock samples corresponding to BOSS CMASS DR11 when we start with an approximation of the RSD. The scheme takes six iterations when the initial estimate, measured from the redshift-space overdensity, has no RSD correction. Slower convergence would be expected for surveys covering a larger angle on the sky. We show that this FFT based method provides a better estimate of the real space displacement field than a configuration space method that uses finite difference routines to compute the potential for the same grid resolution. Finally we show that a lognormal transform of the overdensity, used as a proxy for the linear overdensity, is beneficial in estimating the full displacement field from a dense sample of tracers. However the lognormal transform of the overdensity does not perform well when estimating the displacements from sparser simulations with a more realistic galaxy density.

preprint2015arXiv

The Clustering of the SDSS DR7 Main Galaxy Sample I: A 4 per cent Distance Measure at z=0.15

We create a sample of spectroscopically identified galaxies with $z < 0.2$ from the Sloan Digital Sky Survey (SDSS) Data Release 7, covering 6813 deg$^2$. Galaxies are chosen to sample the highest mass haloes, with an effective bias of 1.5, allowing us to construct 1000 mock galaxy catalogs (described in Paper II), which we use to estimate statistical errors and test our methods. We use an estimate of the gravitational potential to "reconstruct" the linear density fluctuations, enhancing the Baryon Acoustic Oscillation (BAO) signal in the measured correlation function and power spectrum. Fitting to these measurements, we determine $D_{V}(z_{\rm eff}=0.15) = (664\pm25)(r_d/r_{d,{\rm fid}})$ Mpc; this is a better than 4 per cent distance measurement. This "fills the gap" in BAO distance ladder between previously measured local and higher redshift measurements, and affords significant improvement in constraining the properties of dark energy. Combining our measurement with other BAO measurements from BOSS and 6dFGS galaxy samples provides a 15 per cent improvement in the determination of the equation of state of dark energy and the value of the Hubble parameter at $z=0$ ($H_0$). Our measurement is fully consistent with the Planck results and the $Λ$CDM concordance cosmology, but increases the tension between Planck$+$BAO $H_0$ determinations and direct $H_0$ measurements.

preprint2015arXiv

The Clustering of the SDSS Main Galaxy Sample II: Mock galaxy catalogues and a measurement of the growth of structure from Redshift Space Distortions at $z=0.15$

We measure Redshift-Space Distortions (RSD) in the two-point correlation function of a sample of $63,163$ spectroscopically identified galaxies with $z < 0.2$, an epoch where there are currently only limited measurements, from the Sloan Digital Sky Survey (SDSS) Data Release 7 Main Galaxy Sample. Our sample, which we denote MGS, covers 6,813 deg$^2$ with an effective redshift $z_{eff}=0.15$ and is described in our companion paper (Paper I), which concentrates on BAO measurements. In order to validate the fitting methods used in both papers, and derive errors, we create and analyse 1000 mock catalogues using a new algorithm called PICOLA to generate accurate dark matter fields. Haloes are then selected using a friends-of-friends algorithm, and populated with galaxies using a Halo-Occupation Distribution fitted to the data. Using errors derived from these mocks, we fit a model to the monopole and quadrupole moments of the MGS correlation function. If we assume no Alcock-Paczynski (AP) effect (valid at $z=0.15$ for any smooth model of the expansion history), we measure the amplitude of the velocity field, $fσ_{8}$, at $z=0.15$ to be $0.49_{-0.14}^{+0.15}$. We also measure $fσ_{8}$ including the AP effect. This latter measurement can be freely combined with recent Cosmic Microwave Background results to constrain the growth index of fluctuations, $γ$. Assuming a background $Λ$CDM cosmology and combining with current Baryon Acoustic Oscillation data we find $γ= 0.64 \pm 0.09$, which is consistent with the prediction of General Relativity ($γ\approx 0.55$), though with a slight preference for higher $γ$ and hence models with weaker gravitational interactions.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 10 and 11 galaxy samples

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $Λ$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,σ$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: mock galaxy catalogues for the low-redshift sample

We present one thousand mock galaxy catalogues for the analysis of the Low Redshift Sample (LOWZ, effective redshift z ~ 10.32) of the Baryon Oscillation Spectroscopic Survey Data Releases 10 and 11. These mocks have been created following the PTHalos method of Manera13 et al. (2013) revised to include new developments. The main improvement is the introduction of a redshift dependence in the Halo Occupation Distribution in order to account for the change of the galaxy number density with redshift. These mock catalogues are used in the analyses of the LOWZ galaxy clustering by the BOSS collaboration.

preprint2012arXiv

CMB power spectrum parameter degeneracies in the era of precision cosmology

Cosmological parameter constraints from the CMB power spectra alone suffer several well-known degeneracies. These degeneracies can be broken by numerical artefacts and also a variety of physical effects that become quantitatively important with high-accuracy data e.g. from the Planck satellite. We study degeneracies in models with flat and non-flat spatial sections, non-trivial dark energy and massive neutrinos, and investigate the importance of various physical degeneracy-breaking effects. We test the CAMB power spectrum code for numerical accuracy, and demonstrate that the numerical calculations are accurate enough for degeneracies to be broken mainly by true physical effects (the integrated Sachs-Wolfe effect, CMB lensing and geometrical and other effects through recombination) rather than numerical artefacts. We quantify the impact of CMB lensing on the power spectra, which inevitably provides degeneracy-breaking information even without using information in the non-Gaussianity. Finally we check the numerical accuracy of sample-based parameter constraints using CAMB and CosmoMC. In an appendix we document recent changes to CAMB's numerical treatment of massive neutrino perturbations, which are tested along with other recent improvements by our degeneracy exploration results.