Source author record

Geraint Pratten

Geraint Pratten 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

21works
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

21 published item(s)

preprint2026arXiv

Improved parametrized test of general relativity using the IMRPhenomX waveform family: Including higher harmonics and precession

When testing general relativity (GR) with gravitational wave observations, parametrized tests of deviations from the expected strong-field source dynamics are one of the most widely used techniques. We present an updated version of the parametrized framework with the state-of-art IMRPhenomX waveform family. Our new framework incorporates deviations in the dominant mode as well as in the higher-order modes of the waveform. We demonstrate that the missing physics of either higher-order modes or precession in the parametrized model can lead to a biased conclusion of false deviation from GR. Our new implementation mitigates this issue and enables us to perform the tests for highly asymmetric and precessing binaries without being subject to systematic biases due to missing physics. Finally, we apply the improved test to analyze events observed during the second half of the third observing run of LIGO and Virgo (O3b). We provide constraints on GR deviations by combining O3b results with those from previous observation runs. Our findings show no evidence for violations of GR.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

Prospects for distinguishing dynamical tides in inspiralling binary neutron stars with third generation gravitational-wave detectors

Tidal effects in gravitational-wave (GW) observations from binary neutron star mergers have the potential to probe ultra-dense matter and shed light on the unknown nuclear equation of state of neutron stars. Tidal effects in inspiralling neutron star binaries become relevant at GW frequencies of a few hundred Hz and require detectors with exquisite high-frequency sensitivity. Third generation GW detectors such as the Einstein Telescope or Cosmic Explorer will be particularly sensitive in this high-frequency regime, allowing us to probe neutron star tides beyond the adiabatic approximation. Here we assess whether dynamical tides can be measured from a neutron star inspiral. We find that the measurability of dynamical tides depends strongly on the neutron star mass and equation of state. For a semi-realistic population of 10,000 inspiralling binary neutron stars, we conservatively estimate that on average $\mathcal{O}(50)$ binaries will have measurable dynamical tides. As dynamical tides are characterised not only by the star's tidal deformability but also by its fundamental ($f$-) mode frequency, they present a possibility of probing higher-order tidal effects and test consistency with quasi-universal relations. For a GW170817-like signal in a third generation detector network, we find that the stars' $f$-mode frequencies can be measured to within a few hundred Hz.

preprint2022arXiv

The last three years: multiband gravitational-wave observations of stellar-mass binary black holes

Understanding the formation and evolution of the stellar-mass binary black holes discovered by LIGO and Virgo is a challenge that spans many areas of astrophysics, from stellar evolution, dynamics and accretion disks, to possible exotic early universe processes. Over the final years of their lives, stellar-mass binaries radiate gravitational waves that are first observable by space-based detectors (such as LISA) and then ground-based instruments (such as LIGO, Virgo and the next generation observatories Cosmic Explorer and the Einstein Telescope). Using state-of-the-art waveform models and parameter-estimation pipelines for both ground- and space-based observations, we show that (the expected handful of) these multiband observations will allow at least percent-level measurements of all 17 parameters that describe the binary, the possible identification of a likely host galaxy, and the forewarning of the merger days in advance allowing telescopes at multiple wavelengths to search for any electromagnetic signature associated to it. Multiband sources will therefore be a gold mine for astrophysics, but we also show that they could be less useful as laboratories for fundamental tests of general relativity than has been previously suggested.

preprint2021arXiv

Mode-by-mode Relative Binning: Fast Likelihood Estimation for Gravitational Waveforms with Spin-Orbit Precession and Multiple Harmonics

Faster likelihood evaluation enhances the efficiency of gravitational wave signal analysis. We present Mode-by-mode Relative Binning (MRB), a new method designed for obtaining fast and accurate likelihoods for advanced waveform models that include spin-orbit precession effects and multiple radiation harmonics from compact binary coalescence. Leveraging the "twisting-up" procedure of constructing precessing waveform modes from non-precessing ones, the new method mitigates degrade of relative binning accuracy due to interference from superimposed modes. Additionally, we supplement algorithms for optimizing the choice of frequency bins specific to any given strain signal under analysis. Using the new method, we are able to evaluate the likelihood with up to an order of magnitude reduction in the number of waveform model calls per frequency compared to the previously used relative binning scheme, and achieve better likelihood accuracy than is sufficient for obtaining source parameter posterior distributions that are indistinguishable from the exact ones.

preprint2020arXiv

A first survey of spinning eccentric black hole mergers: Numerical relativity simulations, hybrid waveforms, and parameter estimation

We analyze a new numerical relativity dataset of spinning but nonprecessing binary black holes on eccentric orbits, with eccentricities from approximately $0.1$ to $0.5$, with dimensionless spins up to $0.75$ included at mass ratios $q=m_1/m_2 = (1, 2)$, and further nonspinning binaries at mass ratios $q = (1.5, 3, 4)$. A comparison of the final mass and spin of these simulations with noneccentric data extends previous results in the literature on circularization of eccentric binaries to the spinning case. For the $(l,m)=(2,2)$ spin-weighted spherical harmonic mode we construct eccentric hybrid waveforms that connect the numerical relativity data to a post-Newtonian description for the inspiral, and we discuss the limitations in the current knowledge about post-Newtonian theory which complicate the generation of eccentric hybrid waveforms. We also perform a Bayesian parameter estimation study, quantifying the parameter biases introduced when using three different quasicircular waveform models to estimate the parameters of highly eccentric binary systems. We find that the used aligned-spin quasicircular model including higher order modes produces lower bias in certain parameters than the nonprecessing quasicircular model without higher order modes and the quasicircular precessing model.

preprint2020arXiv

A multipolar effective one body waveform model for spin-aligned black hole binaries

We introduce \TEOBiResumSM{}, an improved version of the effective-one-body (EOB) waveform model \TEOBResumS{} for spin-aligned, coalescing black hole binaries, that includes subdominant gravitational waveform modes completed through merger and ringdown. Beyond the dominant $(\ell,|m|)=(2,2)$ one, the more robust multipoles all over the parameter space are: $(2,1)$, $(3,3)$, $(3,2)$, $(4,4)$ and $(5,5)$. The multipolar ringdown EOB waveform stems from suitably fitting many numerical relativity (NR) waveform data from the Simulating eXtreme Spacetimes (SXS) collaboration together with test-mass waveform data. Mode-mixing effects are not incorporated. The orbital (nonspinning) part of the multipolar waveform amplitudes includes test-mass results up to (relative) 6PN order and, for most modes,is Padé resummed. The $m$=odd waveform multipoles (up to $\ell=5$) incorporate most of the currently available spin-dependent analytical information. Improving on previous work, we confirm that certain $m=\text{odd}$ modes, e.g. the $(2,1)$, and even the $(3,1)$, may develop a zero (or a minimum) in the amplitude for nearly equal-mass binaries and for several combinations of the individual spins. A remarkable EOB/NR agreement around such zero is found for these modes. The new waveform, and radiation reaction, prompts a new NR-calibration of the spinning sector of the model, done with only $32$ datasets.The maximum $(2,2)$ EOB/NR unfaithfulness $\bar{F}$ with Advanced LIGO noise against the SXS catalog ($\sim 595$ datasets) is always below $0.5\%$ for binaries with total mass $M$ as $10M_\odot\leq M \leq 200M_\odot$, except for a single outlier with $\max{(\bar{F})}\sim 0.85\%$. When $(2,1)$, $(3,3)$ and $(4,4)$ modes are included, one finds an excellent EOB/NR agreement up to $M\sim 120M_\odot$, above which the performance degrades slightly and moves above $3\%$

preprint2020arXiv

Assessing gravitational-wave binary black hole candidates with Bayesian odds

Gravitational waves from the coalescence of binary black holes can be distinguished from noise transients in a detector network through Bayesian model selection by exploiting the coherence of the signal across the network. We present a Bayesian framework for calculating the posterior probability that a signal is of astrophysical origin, agnostic to the specific search strategy, pipeline or search domain with which a candidate is identified. We apply this framework under \textit{identical} assumptions to all events reported in the LIGO-Virgo GWTC-1 catalog, GW190412 and numerous event candidates reported by independent search pipelines by other authors. With the exception of GW170818, we find that all GWTC-1 candidates, and GW190412, have odds overwhelmingly in favour of the astrophysical hypothesis, including GW170729, which was assigned significantly different astrophysical probabilities by the different search pipelines used in GWTC-1. GW170818 is de-facto a single detector trigger, and is therefore of no surprise that it is disfavoured as being produced by an astrophysical source in our framework. We find \textit{three} additional event candidates, GW170121, GW170425 and GW170727, that have significant support for the astrophysical hypothesis, with a probability that the signal is of astrophysical origin of 0.53, 0.74 and 0.64 respectively. We carry out a hierarchical population study which includes these three events in addition to those reported in GWTC-1, finding that the main astrophysical results are unaffected.

preprint2020arXiv

Setting the cornerstone for the IMRPhenomX family of models for gravitational waves from compact binaries: The dominant harmonic for non-precessing quasi-circular black holes

In this paper we present IMRPhenomXAS, a thorough overhaul of the IMRPhenomD [1,2] waveform model, which describes the dominant $l=2, \:| m | = 2$ spherical harmonic mode of non-precessing coalescing black holes in terms of piecewise closed form expressions in the frequency domain. Improvements include in particular the accurate treatment of unequal spin effects, and the inclusion of extreme mass ratio waveforms. IMRPhenomD has previously been extended to approximately include spin precession [3] and subdominant spherical harmonics [4], and with its extensions it has become a standard tool in gravitational wave parameter estimation. Improved extensions of IMRPhenomXAS are discussed in companion papers [5,6].

preprint2020arXiv

Simple procedures to reduce eccentricity of binary black hole simulations

We present simple procedures to construct quasi-circular initial data for numerical evolutions of binary black hole spacetimes. Our method consists of using Post-Newtonian theory in three ways: first to provide an initial guess for the initial momenta at 3.5PN order that implies low residual eccentricity, second to measure the resulting eccentricity, and third to calculate corrections to the momenta or initial separation which further reduce the eccentricity. Regarding the initial guess, we compare numerical evolutions in post-Newtonian theory to the post-circular and post-post-circular analytical approximations to quasi-circular data. We discuss a robust fitting procedure to measure eccentricity from numerical simulations using the orbital frequency, and derive from the quasi-Keplerian parametrization at 1PN oder the correction factors for the tangential and radial momentum components required to achieve reduce the measured eccentricity to zero. We first test our procedure integrating PN equations of motion at $3.5$PN where low eccentric initial data is easily obtained, and then apply our method to sets of binary black hole numerical relativity simulations with different mass ratios (q=1,2,...,8), spin configurations and separations. Our set of simulations contains non-spinning, spin-aligned and precessing simulations. We observe that the iterative procedure produces low eccentric simulations with eccentricities of the order (10^{-4}) with only one iteration. The simplicity of the procedure allows to obtain low eccentric NR simulations easily and saving computational resources. Moreover, the analytical PN formulas derived in this paper will be useful to generate eccentric hybrid waveforms.

preprint2020arXiv

Towards the routine use of subdominant harmonics in gravitational-wave inference: re-analysis of GW190412 with generation X waveform models

We re-analyse the gravitational-wave event GW190412 with state-of-the-art phenomenological waveform models. This event, which has been associated with a black hole merger, is interesting due to the significant contribution from subdominant harmonics. We use both frequency-domain and time-domain waveform models. The PhenomX waveform models constitute the fourth generation of frequency-domain phenomenological waveforms for black hole binary coalescence; they have more recently been complemented by the time-domain PhenomT models, which open up new strategies to model precession and eccentricity, and to perform tests of general relativity with the phenomenological waveforms approach. Both PhenomX and PhenomT have been constructed with similar techniques and accuracy goals, and due to their computational efficiency this "generation X" model family allows the routine use of subdominant spherical harmonics in Bayesian inference. We show the good agreement between these and other state-of-the-art waveform models for GW190412, and discuss the improvements over the previous generation of phenomenological waveform models. We also discuss practical aspects of Bayesian inference such as run convergence, variations of sampling parameters, and computational cost.

preprint2020arXiv

Validity of common modelling approximations for precessing binary black holes with higher-order modes

The current paradigm for constructing waveforms from precessing compact binaries is to first construct a waveform in a non-inertial, co-precessing binary source frame followed by a time-dependent rotation to map back to the physical, inertial frame. A key insight in the construction of these models is that the co-precessing waveform can be effectively mapped to some equivalent aligned spin waveform. Secondly, the time-dependent rotation implicitly introduces $m$-mode mixing, necessitating an accurate description of higher-order modes in the co-precessing frame. We assess the efficacy of this modelling strategy in the strong field regime using Numerical Relativity simulations. We find that this framework allows for the highly accurate construction of $(2,\pm 2)$ modes in our data set, while for higher order modes, especially the $(2,|1|), (3,|2|)$ and $(4,|3|)$ modes, we find rather large mismatches. We also investigate a variant of the approximate map between co-precessing and aligned spin waveforms, where we only identify the slowly varying part of the time domain co-precessing waveforms with the aligned-spin one, but find no significant improvement. Our results indicate that the simple paradigm to construct precessing waveforms does not provide an accurate description of higher order modes in the strong-field regime, and demonstrate the necessity for modelling mode asymmetries and mode-mixing to significantly improve the description of precessing higher order modes.

preprint2019arXiv

A Multipolar Effective One Body Model for Non-Spinning Black Hole Binaries

We introduce \TEOBiResumSM{}, a nonspinning inspiral-merger-ringdown waveform model built within the effective one body (EOB) framework that includes gravitational waveform modes beyond the dominant quadrupole $(\ell,|m|) = (2,2)$. The model incorporates: (i) an improved Padé resummation of the factorized waveform amplitudes $ρ_{\ell m}^{\rm orb}$ entering the EOB-resummed waveform where the 3PN, mass-ratio dependent, terms are hybridized with test-mass limit terms up to 6PN relative order for most of the multipoles up to $\ell=6$ included; (ii) an improved determination of the effective 5PN function $a_6^c(ν)$ entering the EOB interaction potential done using the most recent, error-controlled, nonspinning numerical relativity (NR) waveforms from the Simulating eXtreme Spacetimes (SXS) collaboration; and (iii) a NR-informed phenomenological description of the multipolar ringdown. Such representation stems from 19 NR waveforms with mass ratios up to $m_1/m_2=18$ as well as test-mass waveform data, although it does not incorporate mode-mixing effects. The NR-completed higher modes through merger and ringdown considered here are: $(\ell,|m|) = \lbrace (2,1), (3,3), (3,2),(3,1),(4,4), (4,3),(4,2), (4,1),(5,5)\rbrace$. For simplicity, the other subdominant modes, up to $\ell=8$, are approximated by the corresponding, purely analytical, factorized and resummed EOB waveform. To attempt an estimate of (some of) the underlying analytic uncertainties of the model, we also contrast the effect of the 6PN-hybrid Padé-resummed $ρ_{\ell m}$'s with the standard $3^{+2}$PN, Taylor-expanded, ones used in previous EOB works. The maximum unfaithfulness $\bar{F}$ against the SXS waveforms including all NR-completed modes up to $\ell=m=5$ is always $\lesssim 2\%$ for binaries with total mass $M$ as $50 M_{\odot} \leq M \lesssim 200 M_{\odot}$.

preprint2019arXiv

Gravitational-Wave Asteroseismology with Fundamental Modes from Compact Binary Inspirals

The first detection of gravitational waves (GWs) from the binary neutron star (NS) inspiral GW170817 has opened a unique channel for probing the fundamental properties of matter at supra-nuclear densities inaccessible elsewhere in the Universe. This observation yielded the first constraints on the equation of state (EoS) of NS matter from the GW imprint of tidal interactions. Tidal signatures in the GW arise from the response of a matter object to the spacetime curvature sourced by its binary companion. They crucially depend on the EoS and are predominantly characterised by the tidal deformability parameters $Λ_{\ell}$, where $\ell=2,3$ denotes the quadrupole and octupole respectively. As the binary evolves towards merger, additional dynamical tidal effects become important when the orbital frequency approaches a resonance with the stars' internal oscillation modes. Among these modes, the fundamental ($f_\ell$-)modes have the strongest tidal coupling and can give rise to a cumulative imprint in the GW signal even if the resonance is not fully excited. Here we present the first direct constraints on fundamental oscillation mode frequencies for GW170817 using an inspiral GW phase model with an explicit dependence on the $f$-mode frequency and without assuming any relation between $f_\ell$ and $Λ_\ell$. We rule out anomalously small values of $f_\ell$ and, for the larger companion, determine a lower bound on the $f_2$-mode ($f_3$-mode) frequency of $\geq 1.39$ kHz ($\geq 1.86$ kHz) at the 90\% credible interval (CI). We then show that networks of future GW detectors will be able to measure $f$-mode frequencies to within tens of Hz from the inspiral alone. Such precision astroseismology will enable novel tests of fundamental physics and the nature of compact binaries.

preprint2016arXiv

3D Weak Lensing: Modified Theories of Gravity

Weak lensing (WL) promises to be a particularly sensitive probe of both the growth of large scale structure (LSS) as well as the fundamental relation between matter density perturbations and metric perturbations, thus providing a powerful tool with which we may constrain modified theories of gravity (MG) on cosmological scales. Future deep, wide-field WL surveys will provide an unprecedented opportunity to constrain deviations from General Relativity (GR). Employing a three-dimensional (3D) analysis based on the spherical Fourier-Bessel (sFB) expansion, we investigate the extent to which MG theories will be constrained by a typical 3D WL survey configuration including noise from the intrinsic ellipticity distribution $σ_ε$ of source galaxies. Here we focus on two classes of screened theories of gravity: i) $f(R)$ chameleon models and ii) environmentally dependent dilaton models. We use one-loop perturbation theory combined with halo models in order to accurately model the evolution of matter power-spectrum with redshift in these theories. Using a Fisher information matrix based approach, we show that for an all-sky spectroscopic survey, the parameter $f_{R_0}$ can be constrained in the range $f_{R_0}< 5\times 10^{-6}(9\times 10^{-6})$ for $n=1(2)$ with a 3$σ$ confidence level. This can be achieved by using relatively low order angular harmonics $\ell<100$. Including higher order harmonics $\ell>100$ can further tighten the constraints, making them comparable to current solar-system constraints. We also employ a Principal Component Analysis (PCA) in order to study the parameter degeneracies in the MG parameters. Our results can trivially be extended to other MG theories, such as the K-mouflage models. The confusion from intrinsic ellipticity correlation and modification of the matter power-spectrum at small scale due to feedback mechanisms is briefly discussed.

preprint2015arXiv

Covariant Perturbations of f(R) Black Holes: The Weyl Terms

In this paper we revisit non-spherical perturbations of the Schwarzschild black hole in the context of $f(R)$ gravity. Previous studies were able to demonstrate the stability of the $f(R)$ Schwarzschild black hole against gravitational perturbations in both the even and odd parity sectors. In particular, it was seen that the Regge-Wheeler and Zerilli equations in $f(R)$ gravity obey the same equations as their General Relativity counterparts. More recently, the 1+1+2 semi-tetrad formalism has been used to derive a set of two wave equations: one for transverse, trace-free (tensor) perturbations and one for the additional scalar modes that characterise fourth-order theories of gravitation. The master variable governing tensor perturbations was shown to be a modified Regge-Wheeler tensor obeying the same equation as in General Relativity. However, it is well known that there is a non-uniqueness in the definition of the master variable. In this paper we derive a set of two perturbation variables and their concomitant wave equations that describe gravitational perturbations in a covariant and gauge invariant manner. These variables can be related to the Newman-Penrose (NP) Weyl scalars as well as the master variables from the 2+2 formalism.

preprint2015arXiv

Galaxy Clustering in 3D and Modified Gravity Theories

We study Modified Gravity (MG) theories by modelling the redshifted matter power spectrum in a spherical Fourier-Bessel (sFB) basis. We use a fully non-linear description of the real-space matter power-spectrum and include the lowest-order redshift-space correction (Kaiser effect), taking into account some additional non-linear contributions. Ignoring relativistic corrections, which are not expected to play an important role for a shallow survey, we analyse two different modified gravity scenarios, namely the generalised Dilaton scalar-tensor theories and the $f({R})$ models in the large curvature regime. We compute the 3D power spectrum ${\cal C}^s_{\ell}(k_1,k_2)$ for various such MG theories with and without redshift space distortions, assuming precise knowledge of background cosmological parameters. Using an all-sky spectroscopic survey with Gaussian selection function $φ(r)\propto \exp(-{r^2 / r^2_0})$, $r_0 = 150 \, h^{-1} \, {\textrm{Mpc}}$, and number density of galaxies $\bar {\textrm{N}} =10^{-4}\;{\textrm{Mpc}}^{-3}$, we use a $χ^2$ analysis, and find that the lower-order $(\ell \leq 25)$ multipoles of ${\cal C}^s_\ell(k,k')$ (with radial modes restricted to $k < 0.2 \, h \,{\textrm{Mpc}}^{-1}$) can constraint the parameter $f_{R_0}$ at a level of $2\times 10^{-5} (3\times 10^{-5})$ with $3 σ$ confidence for $n=1(2)$. Combining constraints from higher $\ell > 25$ modes can further reduce the error bars and thus in principle make cosmological gravity constraints competitive with solar system tests. However this will require an accurate modelling of non-linear redshift space distortions. Using a tomographic $β(a)$-$m(a)$ parameterization we also derive constraints on specific parameters describing the Dilaton models of modified gravity.

preprint2014arXiv

A simple model of complete precessing black-hole-binary gravitational waveforms

The construction of a model of the gravitational-wave (GW) signal from generic configurations of spinning-black-hole binaries, through inspiral, merger and ringdown, is one of the most pressing theoretical problems in the build-up to the era of GW astronomy. We present the first such model in the frequency domain, "PhenomP", which captures the basic phenomenology of the seven-dimensional parameter space of binary configurations with only three key physical parameters. Two of these (the binary's mass ratio and an effective total spin parallel to the orbital angular momentum, which determines the inspiral rate) define an underlying non-precessing-binary model. The non-precessing-binary waveforms are then "twisted up" with approximate expressions for the precessional motion, which require only one additional physical parameter, an effective precession spin, $χ_p$. All other parameters (total mass, sky location, orientation and polarisation, and initial phase) can be specified trivially. The model is constructed in the frequency domain, which will be essential for efficient GW searches and source measurements. We have tested the model's fidelity for GW applications by comparison against hybrid post-Newtonian-numerical-relativity waveforms at a variety of configurations --although we did not use these numerical simulations in the construction of the model. Our model can be used to develop GW searches, to study the implications for astrophysical measurements, and as a simple conceptual framework to form the basis of generic-binary waveform modelling in the advanced-detector era.

preprint2014arXiv

Reconstructing the Thermal Sunyaev-Zel'dovich Effect in 3D

The thermal Sunyaev-Zel'dovich (tSZ) effect measures the line-of-sight projection of the thermal pressure of free electrons and lacks any redshift information. By cross correlating the tSZ effect with an external cosmological tracer we can recover a good fraction of this lost information. Weak lensing (WL) is thought to provide an unbiased probe of the dark Universe, with many WL surveys having sky coverage that overlaps with tSZ surveys. Generalising the tomographic approach, we advocate the use of the spherical Fourier-Bessel (sFB) expansion to perform an analysis of the cross-correlation between the projected (2D) tSZ Compton $y$-parameter maps and 3D weak lensing convergence maps. We use redshift dependent linear biasing and the halo model as a tool to investigate the tSZ-WL cross-correlations in 3D. We use the Press-Schechter (PS) and the Sheth-Tormen (ST) mass-functions in our calculations, finding that the results are quite sensitive to detailed modelling. We provide detailed analysis of surveys with photometric and spectroscopic redshifts. The signal-to-noise (S/N) of the cross-spectra $\mathcal{C}_{\ell} (k)$ for individual 3D modes, defined by the radial and tangential wave numbers $(k;\ell)$, remains comparable to, but below, unity though optimal binning is expected to improve this. The results presented can be generalised to analyse other CMB secondaries, such as the kinetic Sunyaev-Zel'dovich (kSZ) effect.

preprint2013arXiv

Effects of Linear Redshift Space Distortions and Perturbation Theory on BAOs: A 3D Spherical Analysis

The Baryon Acoustic Oscillations (BAO) are features in the matter power spectrum on scales of order 100-150 h^{-1} Mpc that promise to be a powerful tool to constrain and test cosmological models. The BAO have attracted such attention that future upcoming surveys have been designed with the BAO at the forefront of the primary science goals. Recent studies have advocated the use of a spherical-Fourier Bessel (sFB) expansion for future wide field surveys that cover both wide and deep regions of the sky necessitating the simultaneous treatment of the spherical sky geometry as well as the extended radial coverage. Ignoring the possible effects of growth, which is not expected to be significant at low redshifts, we present an extended analysis of the BAOs using the sFB formalism by taking into account the role of non-linearities and linear redshift distortions in the oscillations observed in the galaxy power spectrum. The sFB power spectrum has both radial and tangential dependence and it has been shown that in the limit that we approach a deep survey the sFB power spectrum is purely radial and collapses to the Cartesian Fourier power spectrum. This radialisation of information is shown to hold even in the presence of redshift space distortions (RSD) and 1-loop corrections to the galaxy power spectrum albeit with modified tangential and radial dependence. As per previous studies we find that the introduction of non-linearities leads to a damping of the oscillations in the matter power spectrum.

preprint2012arXiv

Non-Gaussianity in Large Scale Structure and Minkowski Functionals

Minkowski Functionals (MFs) are topological statistics that have become one of many standard tools used for investigating the statistical properties of cosmological random fields. They have found regular use in studies of departures from Gaussianity in a number of important cosmological scenarios. Important examples include the Cosmic Microwave Background (CMB), weak lensing studies, 21cm surveys and large scale structure (LSS). To lowest order the MFs depend on three generalised skewness parameters that can be shown to probe the bispectrum with differing weights. Recent studies have advocated the use of a power spectrum associated with the bispectrum, called the skew-spectrum, that has more power to distinguish between various contributions to the bispectrum than the conventional formalism adopted when using the Minkowski Functionals. In this article we review the motivations for studying non-Gaussianity and emphasize the importance of the momentum dependence of higher order correlators in investigating both inflationary and early Universe models as well as analytical models for gravitational instability. We then introduce the skew-spectra, applied to galaxy surveys, as a tool for investigating various models for primordial and gravitationally induced non-Gaussianities. We present analytical expressions for the skew-spectra for the density field and divergence of the velocity field in 3D and for projected surveys as a function of redshift and a smoothing angular scale. A Gaussian window function is assumed throughout this paper. Analytical results are derived for the case of gravitationally induced non-Gaussianity. These results can be generalised to incorporate redshift space effects. This will be useful in probing primordial and gravitationally induced non-Gaussianity from ongoing and future galaxy surveys.