General CMB and Primordial Trispectrum Estimation

We present trispectrum estimation methods which can be applied to general non-separable primordial and CMB trispectra. We present a general optimal estimator for the connected part of the trispectrum, for which we derive a quadratic term to incorporate the effects of inhomogeneous noise and masking. We describe a general algorithm for creating simulated maps with given arbitrary (and independent) power spectra, bispectra and trispectra. We propose a universal definition of the trispectrum parameter $T_{NL}$, so that the integrated bispectrum on the observational domain can be consistently compared between theoretical models. We define a shape function for the primordial trispectrum, together with a shape correlator and a useful parametrisation for visualizing the trispectrum. We derive separable analytic CMB solutions in the large-angle limit for constant and local models. We present separable mode decompositions which can be used to describe any primordial or CMB bispectra on their respective wavenumber or multipole domains. By extracting coefficients of these separable basis functions from an observational map, we are able to present an efficient estimator for any given theoretical model with a nonseparable trispectrum. The estimator has two manifestations, comparing the theoretical and observed coefficients at either primordial or late times. These mode decomposition methods are numerically tractable with order $l^5$ operations for the CMB estimator and approximately order $l^6$ for the general primordial estimator (reducing to order $l^3$ in both cases for a special class of models). We also demonstrate how the trispectrum can be reconstructed from observational maps using these methods.Comment: 38 pages, 9 figures. In v2 Figures 4-7 are altered slightly and some extra references are included in the bibliography. v3 matches version submitted to journal. Includes discussion of special case

Universal Non-Gaussian Initial Conditions for N-body Simulations

In this paper we present the implementation of an efficient formalism for the generation of arbitrary non-Gaussian initial conditions for use in N-body simulations. The methodology involves the use of a separable modal approach for decomposing a primordial bispectrum or trispectrum. This approach allows for the far more efficient generation of the non-Gaussian initial conditions already described in the literature, as well as the generation for the first time of non-separable bispectra and the special class of diagonal-free trispectra. The modal approach also allows for the reconstruction of the spectra from given realisations, a fact which is exploited to provide an accurate consistency check of the simulations.Comment: 7 pages, 3 figure

Rapid Separable Analysis of Higher Order Correlators in Large Scale Structure

We present an efficient separable approach to the estimation and reconstruction of the bispectrum and the trispectrum from observational (or simulated) large scale structure data. This is developed from general CMB (poly-)spectra methods which exploit the fact that the bispectrum and trispectrum in the literature can be represented by a separable mode expansion which converges rapidly (with $n_\textrm{max}={\cal{O}}(30)$ terms). With an effective grid resolution $l_\textrm{max}$ (number of particles/grid points $N=l_\textrm{max}^3$), we present a bispectrum estimator which requires only ${\cal O}(n_\textrm{max} \times l_\textrm{max}^3)$ operations, along with a corresponding method for direct bispectrum reconstruction. This method is extended to the trispectrum revealing an estimator which requires only ${\cal O}(n_\textrm{max}^{4/3} \times l_\textrm{max}^3)$ operations. The complexity in calculating the trispectrum in this method is now involved in the original decomposition and orthogonalisation process which need only be performed once for each model. However, for non-diagonal trispectra these processes present little extra difficulty and may be performed in ${\cal O}(l_\textrm{max}^4)$ operations. A discussion of how the methodology may be applied to the quadspectrum is also given. An efficient algorithm for the generation of arbitrary nonGaussian initial conditions for use in N-body codes using this separable approach is described. This prescription allows for the production of nonGaussian initial conditions for arbitrary bispectra and trispectra. A brief outline of the key issues involved in parameter estimation, particularly in the non-linear regime, is also given

The CMB Bispectrum

We use a separable mode expansion estimator with WMAP data to estimate the bispectrum for all the primary families of non-Gaussian models. We review the late-time mode expansion estimator methodology which can be applied to any non-separable primordial and CMB bispectrum model, and we demonstrate how the method can be used to reconstruct the CMB bispectrum from an observational map. We extend the previous validation of the general estimator using local map simulations. We apply the estimator to the coadded WMAP 5-year data, reconstructing the WMAP bispectrum using $l<500$ multipoles and $n=31$ orthonormal 3D eigenmodes. We constrain all popular nearly scale-invariant models, ensuring that the theoretical bispectrum is well-described by a convergent mode expansion. Constraints from the local model \fnl=54.4\pm 29.4 and the equilateral model \fnl=143.5\pm 151.2 (\Fnl = 25.1\pm 26.4) are consistent with previously published results. (Here, we use a nonlinearity parameter \Fnl normalised to the local case, to allow more direct comparison between different models.) Notable new constraints from our method include those for the constant model \Fnl = 35.1 \pm 27.4 , the flattened model \Fnl = 35.4\pm 29.2, and warm inflation \Fnl = 10.3\pm 27.2. We investigate feature models surveying a wide parameter range in both the scale and phase, and we find no significant evidence of non-Gaussianity in the models surveyed. We propose a measure \barFnl for the total integrated bispectrum and find that the measured value is consistent with the null hypothesis that CMB anisotropies obey Gaussian statistics. We argue that this general bispectrum survey with the WMAP data represents the best evidence for Gaussianity to date and we discuss future prospects, notably from the Planck satellite

Functional Imaging of the Outer Retinal Complex using High Fidelity Imaging Retinal Densitometry

We describe a new technique, high fidelity Imaging Retinal Densitometry (IRD), which probes the functional integrity of the outer retinal complex. We demonstrate the ability of the technique to map visual pigment optical density and synthesis rates in eyes with and without macular disease. A multispectral retinal imaging device obtained precise measurements of retinal reflectance over space and time. Data obtained from healthy controls and 5 patients with intermediate AMD, before and after photopigment bleaching, were used to quantify visual pigment metrics. Heat maps were plotted to summarise the topography of rod and cone pigment kinetics and descriptive statistics conducted to highlight differences between those with and without AMD. Rod and cone visual pigment synthesis rates in those with AMD (v = 0.043 SD 0.019 min-1 and v = 0.119 SD 0.046 min-1, respectively) were approximately half those observed in healthy controls (v = 0.079 SD 0.024 min-1 for rods and v = 0.206 SD 0.069 min-1 for cones). By mapping visual pigment kinetics across the central retina, high fidelity IRD provides a unique insight into outer retinal complex function. This new technique will improve the phenotypic characterisation, diagnosis and treatment monitoring of various ocular pathologies, including AMD

Oscillations in the bispectrum

There exist several models of inflation that produce primordial bispectra that contain a large number of oscillations. In this paper we discuss these models, and aim at finding a method of detecting such bispectra in the data. We explain how the recently proposed method of mode expansion of bispectra might be able to reconstruct these spectra from separable basis functions. Extracting these basis functions from the data might then lead to observational constraints on these models.Comment: 6 pages, 2 figures, submitted to JOP: Conference Series, PASCOS 201

The characteristics of sexual abuse in sport: A multidimensional scaling analysis of events described in media reports

Most research on sexual abuse has been conducted within family settings (Fergusson & Mullen, 1999). In recent years, following several high profile convictions and scandals, research into sexual abuse has also encompassed institutional and community settings such as sport and the church (Gallagher, 2000; Wolfe et al., 2003). Research into sexual abuse in sport, for example, began with both prevalence studies (Kirby & Greaves, 1996; Leahy, Pretty & Tenenbaum, 2002) and qualitative analyses of the processes and experiences of athlete sexual abuse (Brackenridge, 1997; Cense & Brackenridge, 2001, Toftegaard Nielsen, 2001). From such work, descriptions of the modus operandi of abusers in sport, and the experiences and consequences for athlete victims, have been provided, informing both abuse prevention work and coach education. To date, however, no study has provided empirical support for multiple associations or identified patterns of sex offending in sport in ways that might allow comparisons with research-generated models of offending outside sport. This paper reports on an analysis of 159 cases of criminally defined sexual abuse, reported in the print media over a period of 15 years. The main aim of the study was to identify the nature of sex offending in sport focusing on the methods and locations of offences. The data were analysed using multidimensional scaling (MDS), as a data reduction method, in order to identify the underlying themes within the abuse and explore the inter-relationships of behaviour, victim and context variables. The findings indicate that there are specific themes that can be identified within the perpetrator strategies that include ‘intimate’, ‘aggressive’, and ‘’dominant’ modes of interaction. The same patterns that are described here within the specific context of sport are consistent with themes that emerge from similar behavioural analyses of rapists (Canter & Heritage, 1990; Bishopp, 2003) and child molester groups (Canter, Hughes & Kirby, 1998). These patterns show a correspondence to a broader behavioural model – the interpersonal circumplex (e.g., Leary 1957). Implications for accreditation and continuing professional education of sport psychologists are noted

Optimal bispectrum constraints on single-field models of inflation

We use WMAP 9-year bispectrum data to constrain the free parameters of an 'effective field theory' describing fluctuations in single-field inflation. The Lagrangian of the theory contains a finite number of operators associated with unknown mass scales. Each operator produces a fixed bispectrum shape, which we decompose into partial waves in order to construct a likelihood function. Based on this likelihood we are able to constrain four linearly independent combinations of the mass scales. As an example of our framework we specialize our results to the case of 'Dirac-Born-Infeld' and 'ghost' inflation and obtain the posterior probability for each model, which in Bayesian schemes is a useful tool for model comparison. Our results suggest that DBI-like models with two or more free parameters are disfavoured by the data by comparison with single parameter models in the same class

Non-Gaussianities in N-flation

We compute non-Gaussianities in N-flation, a string motivated model of assisted inflation with quadratic, separable potentials and masses given by the Marcenko-Pastur distribution. After estimating parameters characterizing the bi- and trispectrum in the horizon crossing approximation, we focus on the non-linearity parameter $f_{NL}$, a measure of the bispectrum; we compute its magnitude for narrow and broad spreads of masses, including the evolution of modes after horizon crossing. We identify additional contributions due to said evolution and show that they are suppressed as long as the fields are evolving slowly. This renders $\mathcal{N}$-flation indistinguishable from simple single-field models in this regime. Larger non-Gaussianities are expected to arise for fields that start to evolve faster, and we suggest an analytic technique to estimate their contribution. However, such fast roll during inflation is not expected in N-flation, leaving (p)re-heating as the main additional candidate for generating non-Gaussianities.Comment: 27 pages, 4 figures, extended references to match version accepted in JCA

Effects and Detectability of Quasi-Single Field Inflation in the Large-Scale Structure and Cosmic Microwave Background

Quasi-single field inflation predicts a peculiar momentum dependence in the squeezed limit of the primordial bispectrum which smoothly interpolates between the local and equilateral models. This dependence is directly related to the mass of the isocurvatons in the theory which is determined by the supersymmetry. Therefore, in the event of detection of a non-zero primordial bispectrum, additional constraints on the parameter controlling the momentum-dependence in the squeezed limit becomes an important question. We explore the effects of these non-Gaussian initial conditions on large-scale structure and the cosmic microwave background, with particular attention to the galaxy power spectrum at large scales and scale-dependence corrections to galaxy bias. We determine the simultaneous constraints on the two parameters describing the QSF bispectrum that we can expect from upcoming large-scale structure and cosmic microwave background observations. We find that for relatively large values of the non-Gaussian amplitude parameters, but still well within current uncertainties, galaxy power spectrum measurements will be able to distinguish the QSF scenario from the predictions of the local model. A CMB likelihood analysis, as well as Fisher matrix analysis, shows that there is also a range of parameter values for which Planck data may be able distinguish between QSF models and the related local and equilateral shapes. Given the different observational weightings of the CMB and LSS results, degeneracies can be significantly reduced in a joint analysis.Comment: 27 pages, 14 figure