BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ue-Li Pen

Ue-Li Pen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.HEastro-ph.IMastro-phastro-ph.GAhep-phhep-thgr-qcastro-ph.SRphysics.comp-phastro-ph.EPPerformancephysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

87works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

87 published item(s)

preprint2026arXiv

The Second CHIME/FRB Catalog of Fast Radio Bursts

We present a catalog of 4539 fast radio bursts (FRBs) observed with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope between 25 July 2018 and 15 September 2023. These bursts originate from 3641 unique sources, including 981 bursts from 83 known repeating sources. For each FRB, the catalog provides a $O(10')$ estimate of sky location along with corresponding measurements of cumulative exposure time and survey sensitivity over the observing period. It includes a total-intensity dynamic spectrum between 400 and 800 MHz at 0.983 ms resolution. From this spectrum, we constrain a model of the burst morphology and measure key parameters such as arrival time, intrinsic temporal width, dispersion measure, scattering time, and flux density. This second catalog includes all FRBs from the first catalog, with every event reprocessed using a uniform and improved analysis framework. We show that previously published inferences remain valid under the updated measurements. We assess consistency of the detection rate across observational parameters, present initial distributions of burst properties, and outline ongoing and future studies that will use this catalog to investigate the nature of FRBs and their utility as astrophysical and cosmological probes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A High-Time Resolution Search for Compact Objects using Fast Radio Burst Gravitational Lens Interferometry with CHIME/FRB

The gravitational field of compact objects, such as primordial black holes, can create multiple images of background sources. For transients such as fast radio bursts (FRBs), these multiple images can be resolved in the time domain. Under certain circumstances, these images not only have similar burst morphologies but are also phase-coherent at the electric field level. With a novel dechannelization algorithm and a matched filtering technique, we search for repeated copies of the same electric field waveform in observations of FRBs detected by the FRB backend of the Canadian Hydrogen Mapping Intensity Experiment (CHIME). An interference fringe from a coherent gravitational lensing signal will appear in the time-lag domain as a statistically-significant peak in the time-lag autocorrelation function. We calibrate our statistical significance using telescope data containing no FRB signal. Our dataset consists of $\sim$100-ms long recordings of voltage data from 172 FRB events, dechannelized to 1.25-ns time resolution. This coherent search algorithm allows us to search for gravitational lensing signatures from compact objects in the mass range of $10^{-4}-10^{4} ~\mathrm{M_{\odot}}$. After ruling out an anomalous candidate due to diffractive scintillation, we find no significant detections of gravitational lensing in the 172 FRB events that have been analyzed. In a companion work [Leung, Kader+2022], we interpret the constraints on dark matter from this search.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Overview of CHIME, the Canadian Hydrogen Intensity Mapping Experiment

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a drift scan radio telescope operating across the 400-800 MHz band. CHIME is located at the Dominion Radio Astrophysical Observatory near Penticton, BC Canada. The instrument is designed to map neutral hydrogen over the redshift range 0.8 to 2.5 to constrain the expansion history of the Universe. This goal drives the design features of the instrument. CHIME consists of four parallel cylindrical reflectors, oriented north-south, each 100 m $\times$ 20 m and outfitted with a 256 element dual-polarization linear feed array. CHIME observes a two degree wide stripe covering the entire meridian at any given moment, observing 3/4 of the sky every day due to Earth rotation. An FX correlator utilizes FPGAs and GPUs to digitize and correlate the signals, with different correlation products generated for cosmological, fast radio burst, pulsar, VLBI, and 21 cm absorber backends. For the cosmology backend, the $N_\mathrm{feed}^2$ correlation matrix is formed for 1024 frequency channels across the band every 31 ms. A data receiver system applies calibration and flagging and, for our primary cosmological data product, stacks redundant baselines and integrates for 10 s. We present an overview of the instrument, its performance metrics based on the first three years of science data, and we describe the current progress in characterizing CHIME's primary beam response. We also present maps of the sky derived from CHIME data; we are using versions of these maps for a cosmological stacking analysis as well as for investigation of Galactic foregrounds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraining Primordial Black Holes using Fast Radio Burst Gravitational-Lens Interferometry with CHIME/FRB

Fast radio bursts (FRBs) represent an exciting frontier in the study of gravitational lensing, due to their brightness, extragalactic nature, and the compact, coherent characteristics of their emission. In a companion work [Kader, Leung+2022], we use a novel interferometric method to search for gravitationally lensed FRBs in the time domain using bursts detected by CHIME/FRB. There, we dechannelize and autocorrelate electric field data at a time resolution of 1.25 ns. This enables a search for FRBs whose emission is coherently deflected by gravitational lensing around a foreground compact object such as a primordial black hole (PBH). Here, we use our non-detection of lensed FRBs to place novel constraints on the PBH abundance outside the Local Group. We use a novel two-screen model to take into account decoherence from scattering screens in our constraints. Our constraints are subject to a single astrophysical model parameter -- the effective distance between an FRB source and the scattering screen, for which we adopt a fiducial distance of 1 parsec. We find that coherent FRB lensing is a sensitive probe of sub-solar mass compact objects. Having observed no lenses in $172$ bursts from $114$ independent sightlines through the cosmic web, we constrain the fraction of dark matter made of compact objects, such as PBHs, to be $f \lesssim 0.8$, if their masses are $\sim 10^{-3} M_{\odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Correlating galaxy shapes and initial conditions: an observational study

Using data from the Sloan Digital Sky Survey we study correlations between directions of galaxy angular momenta determined from images of spiral galaxies and various observables derived from the reconstructed initial conditions. We find an apparent systematic effect consistent with galaxy-orientation-dependent selection function. After restricting our attention to the brightest half of the galaxies where this systematic effect is presumed to be absent, we find hints of excess/deficit correlation for two observables. Interestingly, tidal torque theory predicts excess/deficit correlation in exactly these two observables. After correcting for the redshift space distortions, the significance of these correlations drops below 3$σ$ threshold. We do not find any other systematic issues, but a thorough systematic analysis goes beyond the scope of this work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Localizing FRBs through VLBI with the Algonquin Radio Observatory 10-m Telescope

The CHIME/FRB experiment has detected thousands of Fast Radio Bursts (FRBs) due to its sensitivity and wide field of view; however, its low angular resolution prevents it from localizing events to their host galaxies. Very Long Baseline Interferometry (VLBI), triggered by FRB detections from CHIME/FRB will solve the challenge of localization for non-repeating events. Using a refurbished 10-m radio dish at the Algonquin Radio Observatory located in Ontario Canada, we developed a testbed for a VLBI experiment with a theoretical ~<30 masec precision. We provide an overview of the 10-m system and describe its refurbishment, the data acquisition, and a procedure for fringe fitting that simultaneously estimates the geometric delay used for localization and the dispersive delay from the ionosphere. Using single pulses from the Crab pulsar, we validate the system and localization procedure, and analyze the clock stability between sites, which is critical for phase-referencing an FRB event. We find a localization of 50 masec is possible with the performance of the current system. Furthermore, for sources with insufficient signal or restricted wideband to simultaneously measure both geometric and ionospheric delays, we show that the differential ionospheric contribution between the two sites must be measured to a precision of 1e-8 pc/cc to provide a reasonable localization from a detection in the 400--800 MHz band. Finally we show detection of an FRB observed simultaneously in the CHIME and the Algonquin 10-m telescope, the first FRB cross-correlated in this very long baseline. This project serves as a testbed for the forthcoming CHIME/FRB Outriggers project.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Measuring lens dimensionality in extreme scattering events through wave optics

Compact radio sources have been observed to undergo large, frequency dependent changes in intensity due to lensing by structures in the interstellar medium, in so-called "extreme scattering events" (ESEs). While the study of astrophysical plasma lensing has primarily focused on the geometric limit of optics, coherent radio sources such as pulsars exhibit wave effects when lensed. The additional phase information provided by interference effects in the wave regime may yield more information about the lens than could be obtained in the geometric regime. In this paper, we show that, using wave effects, one can potentially distinguish a one-dimensional lens (where "one-dimensional" includes both highly elongated lenses, as well as perfectly axisymmetric lenses) from a fully two-dimensional lens, with minimal assumptions on the form of the lensing potential.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observational search for primordial chirality violations using galaxy angular momenta

We search for evidence of primordial chirality violation in the galaxy data from the Sloan Digital Sky Survey by comparing how strongly directions of galaxy angular momenta correlate with left and right helical components of a spin vector field constructed from the initial density perturbations. Within uncertainties, galaxy spins correlate with these two helical components identically, which is consistent with Universe without primordial chirality violation. Given current data, it is not yet possible to rule out maximal chiral violation, although the case of vanishing correlation with the right helical component is ruled out at about 3.8$σ$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Sub-second periodicity in a fast radio burst

Fast radio bursts (FRBs) are millisecond-duration flashes of radio waves that are visible at distances of billions of light-years. The nature of their progenitors and their emission mechanism remain open astrophysical questions. Here we report the detection of the multi-component FRB 20191221A and the identification of a periodic separation of 216.8(1) ms between its components with a significance of 6.5 sigmas. The long (~3 s) duration and nine or more components forming the pulse profile make this source an outlier in the FRB population. Such short periodicity provides strong evidence for a neutron-star origin of the event. Moreover, our detection favours emission arising from the neutron-star magnetosphere, as opposed to emission regions located further away from the star, as predicted by some models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Moving Mirror model for Fast Radio Bursts

Recent observations of coherent radiation from the Crab pulsar (Bij et al 2021) suggest the emission is driven by an ultra - relativistic ($γ\sim 10^4$), cold plasma flow. A relativistically expanding plasma shell can compress the ambient magnetic field, like a moving mirror, and thus produce coherent radiation whose wavelength is shorter than that of the ambient medium by $γ^2$. This mechanism has been studied in the past by Colgate and Noerdelinger (1971), in the context of radio loud supernova explosions. In this work we propose that a similar mechanism drives the coherent emission in fast radio bursts. The high Lorenz factors dramatically lower the implied energy and magnetic field requirements, allowing the spin down energy of regular (or even recycled), fast spinning pulsars, rather than slow spinning magnetars, to explain FRBs. We show that this model can explain the frequency and the time evolution of observed FRBs, as well as their duration, energetics and absence of panchromatic counterparts. We also predict that the peak frequency of sub pulses decline with observation time as $ω_{\rm obs} \propto t_{\rm obs}^{-1/2}$. Unfortunately, with current capabilities it is not possible to constrain the shape of the curve $ω_{\rm obs} \left(t_{\rm obs} \right)$. Finally, we find that a variation of this model can explain weaker radio transients, such as the one observed from a galactic magnetar. In this variant, the shock wave produces low frequency photons which are then Compton scattered to the GHz range.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Cosmic Tidal Reconstruction with Halo Fields

The gravitational coupling between large-scale perturbations and small-scale perturbations leads to anisotropic distortions of the small-scale matter distribution. The measured local small-scale power spectrum can thus be used to infer the large-scale matter distribution. In this paper, we present a new tidal reconstruction algorithm for reconstructing large-scale modes using the full three-dimensional tidal shear information. We apply it to simulated dark matter halo fields and the reconstructed large-scale density field correlates well with the original matter density field on large scales, improving upon the previous tidal reconstruction method which only uses two transverse shear fields. This has profound implications for recovering lost 21~cm radial modes due to foreground subtraction and constraining primordial non-Gaussianity using the multi-tracer method with future cosmological surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Overview and status of EXCLAIM, the experiment for cryogenic large-aperture intensity mapping

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne far-infrared telescope that will survey star formation history over cosmological time scales to improve our understanding of why the star formation rate declined at redshift z < 2, despite continued clustering of dark matter. Specifically,EXCLAIM will map the emission of redshifted carbon monoxide and singly-ionized carbon lines in windows over a redshift range 0 < z < 3.5, following an innovative approach known as intensity mapping. Intensity mapping measures the statistics of brightness fluctuations of cumulative line emissions instead of detecting individual galaxies, thus enabling a blind, complete census of the emitting gas. To detect this emission unambiguously, EXCLAIM will cross-correlate with a spectroscopic galaxy catalog. The EXCLAIM mission uses a cryogenic design to cool the telescope optics to approximately 1.7 K. The telescope features a 90-cm primary mirror to probe spatial scales on the sky from the linear regime up to shot noise-dominated scales. The telescope optical elements couple to six μ-Spec spectrometer modules, operating over a 420-540 GHz frequency band with a spectral resolution of 512 and featuring microwave kinetic inductance detectors. A Radio Frequency System-on-Chip (RFSoC) reads out the detectors in the baseline design. The cryogenic telescope and the sensitive detectors allow EXCLAIM to reach high sensitivity in spectral windows of low emission in the upper atmosphere. Here, an overview of the mission design and development status since the start of the EXCLAIM project in early 2019 is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exploring the dispersion measure of the Milky Way halo

Fast radio bursts offer the opportunity to place new constraints on the mass and density profile of hot and ionized gas in galactic haloes. We test here the X-ray emission and dispersion measure predicted by different gas profiles for the halo of the Milky Way. We examine a range of models, including entropy stability conditions and external pressure continuity. We find that incorporating constraints from X-ray observations leads to favouring dispersion measures on the lower end of the range given by these models. We show that the dispersion measure of the Milky Way halo could be less than 10 cm$^{-3}$ pc in the most extreme model we consider, which is based on constraints from O VII absorption lines. However, the models allowed by the soft X-ray constraints span more than an order of magnitude in dispersion measures. Additional information on the distribution of gas in the Milky Way halo could be obtained from the signature of a dipole in the dispersion measure of fast radio bursts across the sky, but this will be a small effect for most cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing Primordial Chirality with Galaxy Spins

Chiral symmetry is maximally violated in weak interactions, and such microscopic asymmetries in the early Universe might leave observable imprints on astrophysical scales without violating the cosmological principle. In this Letter, we propose a helicity measurement to detect primordial chiral violation. We point out that observations of halo-galaxy angular momentum directions (spins), which are frozen in during the galaxy formation process, provide a fossil chiral observable. From the clustering mode of large scale structure of the Universe, we construct a spin mode in Lagrangian space and show in simulations that it is a good probe of halo-galaxy spins. In standard model, a strong symmetric correlation between the left and right helical components of this spin mode and galaxy spins is expected. Measurements of these correlations will be sensitive to chiral breaking, providing a direct test of chiral symmetry breaking in the early Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wave effects in the microlensing of pulsars and FRBs by point masses

Wave effects are often neglected in microlensing studies; however, for coherent point-like sources, such as pulsars and fast radio bursts (FRBs), wave effects will become important in their gravitational lensing. In this paper, we describe the wave optics formalism, its various limits, and the conditions for which these limits hold. Using the simple point lens as an example, we will show that the frequency dependence of wave effects breaks degeneracies that are present in the usual geometric optics limit, and constructive interference results in larger magnifications further from the lens. This latter fact leads to a generic increase in cross section for microlensing events in the wave-optics regime compared to the geometric optics regime. For realistic percent-level spectral sensitivities, this leads to a relative boost in lensing cross section of more than an order of magnitude. We apply the point-lens model to the lensing of FRBs and pulsars and find that these radio sources will be lensed in the full wave-optics regime by isolated masses in the range of $0.1-100\,M_\oplus$, which includes free-floating planets (FFPs), whose Einstein radius is smaller than the Fresnel scale. More generally, the interference pattern allows an instantaneous determination of lens masses, unlike traditional microlensing techniques which only yield a mass inference from the event timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Canada and the SKA from 2020-2030

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Accurate Polarization Calibration at 800 MHz with the Green Bank Telescope

Polarization leakage of foreground synchrotron emission is a critical issue in HI intensity mapping experiments. While the sought-after HI emission is unpolarized, polarized foregrounds such as Galactic and extragalactic synchrotron radiation, if coupled with instrumental impurity, can mimic or overwhelm the HI signals. In this paper we present the methodology for polarization calibration at 700-900 MHz, applied on data obtained from the Green Bank Telescope (GBT). We use astrophysical sources, both polarized and unpolarized sources including quasars and pulsars, as calibrators to characterize the polarization leakage and control systematic effects in our GBT HI intensity mapping project. The resulting fractional errors on polarization measurements on boresight are well controlled to within 0.6%-0.8% of their total intensity. The polarized beam patterns are measured by performing spider scans across both polarized quasars and pulsars. A dominant Stokes I to V leakage feature and secondary features of Stokes I to Q and I to U leakages in the 700-900 MHz frequency range are identified. These characterizations are important for separating foreground polarization leakage from the HI 21 cm signal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on the FRB rate at 700-900 MHz

Estimating the all-sky rate of fast radio bursts (FRBs) has been difficult due to small-number statistics and the fact that they are seen by disparate surveys in different regions of the sky. In this paper we provide limits for the FRB rate at 800 MHz based on the only burst detected at frequencies below 1.4 GHz, FRB 110523. We discuss the difficulties in rate estimation, particularly in providing an all-sky rate above a single fluence threshold. We find an implied rate between 700-900 MHz that is consistent with the rate at 1.4 GHz, scaling to $6.4^{+29.5}_{-5.0} \times 10^3$\,sky$^{-1}$\,day$^{-1}$ for an HTRU-like survey. This is promising for upcoming experiments below a GHz like CHIME and UTMOST, for which we forecast detection rates. Given 110523's discovery at 32$σ$ with nothing weaker detected, down to the threshold of 8$σ$, we find consistency with a Euclidean flux distribution but disfavour steep distributions, ruling out $γ> 2.2$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

FRB repetition and non-Poissonian statistics

We discuss some of the claims that have been made regarding the statistics of fast radio bursts (FRBs). In an earlier paper \citep{2015arXiv150505535C} we conjectured that flicker noise associated with FRB repetition could show up in non-cataclysmic neutron star emission models, like supergiant pulses. We show how the current limits of repetition would be significantly weakened if their repeat rate really were non-Poissonian and had a pink or red spectrum. Repetition and its statistics have implications for observing strategy, generally favouring shallow wide-field surveys, since in the non-repeating scenario survey depth is unimportant. We also discuss the statistics of the apparent latitudinal dependence of FRBs, and offer a simple method for calculating the significance of this effect. We provide a generalized Bayesian framework for addressing this problem, which allows for direct model comparison. It is shown how the evidence for a steep latitudinal gradient of the FRB rate is less strong than initially suggested and simple explanations like increased scattering and sky temperature in the plane are sufficient to decrease the low-latitude burst rate, given current data. The reported dearth of bursts near the plane is further complicated if FRBs have non-Poissonian repetition, since in that case the event rate inferred from observation depends on observing strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Holographic Beam Mapping of the CHIME Pathfinder Array

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Interpreting the unresolved intensity of cosmologically redshifted line radiation

Intensity mapping experiments survey the spectrum of diffuse line radiation rather than detect individual objects at high signal-to-noise. Spectral maps of unresolved atomic and molecular line radiation contain three-dimensional information about the density and environments of emitting gas, and efficiently probe cosmological volumes out to high redshift. Intensity mapping survey volumes also contain all other sources of radiation at the frequencies of interest. Continuum foregrounds are typically ~10^2-10^3 times brighter than the cosmological signal. The instrumental response to bright foregrounds will produce new spectral degrees of freedom that are not known in advance, nor necessarily spectrally smooth. The intrinsic spectra of foregrounds may also not be well-known in advance. We describe a general class of quadratic estimators to analyze data from single-dish intensity mapping experiments, and determine contaminated spectral modes from the data itself. The key attribute of foregrounds is not that they are spectrally smooth, but instead that they have fewer bright spectral degrees of freedom than the cosmological signal. Spurious correlations between the signal and foregrounds produce additional bias. Compensation for signal attenuation must estimate and correct this bias. A successful intensity mapping experiment will control instrumental systematics that spread variance into new modes, and it must observe a large enough volume that contaminant modes can be determined independently from the signal on scales of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measurement of Neutrino Masses from Relative Velocities

We present a new technique to measure neutrino masses using their flow field relative to dark matter. Present day streaming motions of neutrinos relative to dark matter and baryons are several hundred km/s, comparable with their thermal velocity dispersion. This results in a unique dipole anisotropic distortion of the matter-neutrino cross power spectrum, which is observable through the dipole distortion in the cross correlation of different galaxy populations. Such a dipole vanishes if not for this relative velocity and so it is a clean signature for neutrino mass. We estimate the size of this effect and find that current and future galaxy surveys may be sensitive to these signature distortions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Primordial density and BAO reconstruction

We present a new method to reconstruct the primordial (linear) density field using the estimated nonlinear displacement field. The divergence of the displacement field gives the reconstructed density field. We solve the nonlinear displacement field in the 1D cosmology and show the reconstruction results. The new reconstruction algorithm recovers a lot of linear modes and reduces the nonlinear damping scale significantly. The successful 1D reconstruction results imply the new algorithm should also be a promising technique in the 3D case.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Probing Neutrino Hierarchy and Chirality via Wakes

The relic neutrinos are expected to acquire a bulk relative velocity with respect to the dark matter at low redshifts, and neutrino wakes are expected to develop downstream of the dark matter halos. We propose a method of measuring the neutrino mass based on this mechanism. This neutrino wake will cause a dipole distortion of the galaxy-galaxy lensing pattern. This effect could be detected by combining upcoming lensing surveys with a low redshift galaxy survey or a 21 cm intensity mapping survey, which can map the neutrino flow field. The data obtained with LSST and Euclid should enable us to make a positive detection if the three neutrino masses are quasidegenerate with each neutrino mass of $\sim$0.1 eV, and a future high precision 21 cm lensing survey would allow the normal hierarchy and inverted hierarchy cases to be distinguished, and even the right-handed Dirac neutrinos may be detectable.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pulsar Acceleration Shifts from nearby Supernova Explosion

We show that when a supernova explodes, a nearby pulsar signal goes through a very specific change. The observed period first changes smoothly, then is followed by a sudden change in the time derivative. A stable millisecond pulsar can allow us to measure such an effect. This may improve our measurement of the total energy released in neutrinos and also the orientation of the supernova-pulsar system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shocks in the Early Universe

We point out a surprising consequence of the usually assumed initial conditions for cosmological perturbations. Namely, a spectrum of Gaussian, linear, adiabatic, scalar, growing mode perturbations not only creates acoustic oscillations of the kind observed on very large scales today, it also leads to the production of shocks in the radiation fluid of the very early universe. Shocks cause departures from local thermal equilibrium as well as creating vorticity and gravitational waves. For a scale-invariant spectrum and standard model physics, shocks form for temperatures $1$ GeV$<T<10^{7}$ GeV. For more general power spectra, such as have been invoked to form primordial black holes, shock formation and the consequent gravitational wave emission provides a signal detectable by current and planned gravitational wave experiments, allowing them to strongly constrain conditions present in the primordial universe as early as $10^{-30}$ seconds after the big bang.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Euclidean distribution of Fast Radio Bursts

We investigate whether current data on the distribution of observed flux densities of Fast Radio Bursts (FRBs) are consistent with a constant source density in Euclidean space. We use the number of FRBs detected in two surveys with different characteristics along with the observed signal-to-noise ratios of the detected FRBs in a formalism similar to a V/V_max-test to constrain the distribution of flux densities. We find consistency between the data and a Euclidean distribution. Any extension of this model is therefore not data-driven and needs to be motivated separately. As a byproduct we also obtain new improved limits for the FRB rate at 1.4 GHz, which had not been constrained in this way before.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Efficient Real-time Data Pipeline for the CHIME Pathfinder Radio Telescope X-Engine

The CHIME Pathfinder is a new interferometric radio telescope that uses a hybrid FPGA/GPU FX correlator. The GPU-based X-engine of this correlator processes over 819 Gb/s of 4+4-bit complex astronomical data from N=256 inputs across a 400 MHz radio band. A software framework is presented to manage this real-time data flow, which allows each of 16 processing servers to handle 51.2 Gb/s of astronomical data, plus 8 Gb/s of ancillary data. Each server receives data in the form of UDP packets from an FPGA F-engine over the eight 10 GbE links, combines data from these packets into large (32MB-256MB) buffered frames, and transfers them to multiple GPU co-processors for correlation. The results from the GPUs are combined and normalized, then transmitted to a collection server, where they are merged into a single file. Aggressive optimizations enable each server to handle this high rate of data; allowing the efficient correlation of 25 MHz of radio bandwidth per server. The solution scales well to larger values of N by adding additional servers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dense magnetized plasma associated with a fast radio burst

Fast Radio Bursts are bright, unresolved, non-repeating, broadband, millisecond flashes, found primarily at high Galactic latitudes, with dispersion measures much larger than expected for a Galactic source. The inferred all-sky burst rate is comparable to the core-collapse supernova rate out to redshift 0.5. If the observed dispersion measures are assumed to be dominated by the intergalactic medium, the sources are at cosmological distances with redshifts of 0.2 to 1. These parameters are consistent with a wide range of source models. One fast radio burst showed circular polarization [21(7)%] of the radio emission, but no linear polarization was detected, and hence no Faraday rotation measure could be determined. Here we report the examination of archival data revealing Faraday rotation in a newly detected burst - FRB 110523. It has radio flux at least 0.6 Jy and dispersion measure 623.30(5) pc cm$^{-3}$. Using Galactic contribution 45 pc cm$^{-3}$ and a model of intergalactic electron density, we place the source at a maximum redshift of 0.5. The burst has rotation measure -186.1(1.4) rad m$^{-2}$, much higher than expected for this line of sight through the Milky Way and the intergalactic medium, indicating magnetization in the vicinity of the source itself or within a host galaxy. The pulse was scattered by two distinct plasma screens during propagation, which requires either a dense nebula associated with the source or a location within the central region of its host galaxy. Keeping in mind that there may be more than one type of fast radio burst source, the detection in this instance of source-local magnetization and scattering favours models involving young stellar populations such as magnetars over models involving the mergers of older neutron stars, which are more likely to be located in low density regions of the host galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local Circumnuclear Magnetar Solution to Extragalactic Fast Radio Bursts

We synthesize the known information about Fast Radio Bursts and radio magnetars, and describe an allowed origin near nuclei of external, but non-cosmological, galaxies. This places them at $z\ll1$, within a few hundred megaparsecs. In this scenario, the high DM is dominated by the environment of the FRB, modelled on the known properties of the Milky Way Center, whose innermost 100pc provides 1000 pc/cm$^3$. A radio loud magnetar is known to exist in our galactic center, within $\sim$2 arc seconds of Sgr A*. Based on the polarization, DM, and scattering properties of this known magnetar, we extrapolate its properties to those of Crab-like giant pulses and SGR flares and point out their consistency with observed Fast Radio Bursts. We conclude galactic center magnetars could be the source of FRB's. This scenario is readily testable with VLBI measurements as well as with flux count statistics from large surveys such as CHIME or UTMOST.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-Cosmological FRBs from Young Supernova Remnant Pulsars

We propose a new extragalactic but non-cosmological explanation for fast radio bursts (FRBs) based on very young pulsars in supernova remnants. Within a few hundred years of a core-collapse supernova the ejecta is confined within $\sim$1 pc, providing a high enough column density of free electrons for the observed 375-1600 pc cm$^{-3}$ of dispersion measure (DM). By extrapolating a Crab-like pulsar to its infancy in an environment like that of SN 1987A, we hypothesize such an object could emit supergiant pulses sporadically which would be bright enough to be seen at a few hundred megaparsecs. We hypothesize that such supergiant pulses would preferentially occur early in the pulsar's life when the free electron density is still high, which is why we do not see large numbers of moderate DM FRBs ($\lesssim300$ pc cm$^{-3}$). In this scenario Faraday rotation at the source gives rotation measures (RMs) much larger than the expected cosmological contribution. If the emission were pulsar-like, then the polarization vector could swing over the duration of the burst, which is not expected from non-rotating objects. In this model, the scattering, large DM, and commensurate RM all come from one place which is not the case for the cosmological interpretation. The model also provides testable predictions of the flux distribution and repeat rate of FRBs, and could be further verified by spatial coincidence with optical supernovae of the past several decades and cross-correlation with nearby galaxy maps.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Precision reconstruction of the dark matter-neutrino relative velocity from N-body simulations

Discovering the mass of neutrinos is a principle goal in high energy physics and cosmology. In addition to cosmological measurements based on two-point statistics, the neutrino mass can also be estimated by observations of neutrino wakes resulting from the relative motion between dark matter and neutrinos. Such a detection relies on an accurate reconstruction of the dark matter-neutrino relative velocity which is affected by non-linear structure growth and galaxy bias. We investigate our ability to reconstruct this relative velocity using large N-body simulations where we evolve neutrinos as distinct particles alongside the dark matter. We find that the dark matter velocity power spectrum is overpredicted by linear theory whereas the neutrino velocity power spectrum is underpredicted. The magnitude of the relative velocity observed in the simulations is found to be lower than what is predicted in linear theory. Since neither the dark matter nor the neutrino velocity fields are directly observable from galaxy or 21 cm surveys, we test the accuracy of a reconstruction algorithm based on halo density fields and linear theory. Assuming prior knowledge of the halo bias, we find that the reconstructed relative velocities are highly correlated with the simulated ones with correlation coefficients of 0.94, 0.93, 0.91 and 0.88 for neutrinos of mass 0.05, 0.1, 0.2 and 0.4 eV. We confirm that the relative velocity field reconstructed from large scale structure observations such as galaxy or 21 cm surveys can be accurate in direction and, with appropriate scaling, magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pulsar lensing geometry

Our analysis of archival VLBI data of PSR 0834+06 revealed that its scintillation properties can be precisely modelled using the inclined sheet model (Pen & Levin 2014), resulting in two distinct lens planes. These data strongly favour the grazing sheet model over turbulence as the primary source of pulsar scattering. This model can reproduce the parameters of the observed diffractive scintillation with an accuracy at the percent level. Comparison with new VLBI proper motion results in a direct measure of the ionized ISM screen transverse velocity. The results are consistent with ISM velocities local to the PSR 0834+06 sight-line (through the Galaxy). The simple 1D structure of the lenses opens up the possibility of using interstellar lenses as precision probes for pulsar lens mapping, precision transverse motions in the ISM, and new opportunities for removing scattering to improve pulsar timing. We describe the parameters and observables of this double screen system. While relative screen distances can in principle be accurately determined, a global conformal distance degeneracy exists that allows a rescaling of the absolute distance scale. For PSR B0834+06, we present VLBI astrometry results that provide (for the fist time) a direct measurement of the distance of the pulsar. For targets where independent distance measurements are not available, which are the cases for most of the recycled millisecond pulsars that are the targets of precision timing observations, the degeneracy presented in the lens modelling could be broken if the pulsar resides in a binary system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Calibrating CHIME, A New Radio Interferometer to Probe Dark Energy

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coaxing Cosmic 21cm Fluctuations from the Polarized Sky using m-mode Analysis

In this paper we continue to develop the m-mode formalism, a technique for efficient and optimal analysis of wide-field transit radio telescopes, targeted at 21 cm cosmology. We extend this formalism to give an accurate treatment of the polarised sky, fully accounting for the effects of polarisation leakage and cross-polarisation. We use the geometry of the measured set of visibilities to project down to pure temperature modes on the sky, serving as a significant compression, and an effective first filter of polarised contaminants. We use the m-mode formalism with the Karhunen-Loeve transform to give a highly efficient method for foreground cleaning, and demonstrate its success in cleaning realistic polarised skies observed with an instrument suffering from substantial off axis polarisation leakage. We develop an optimal quadratic estimator in the m-mode formalism, which can be efficiently calculated using a Monte-Carlo technique. This is used to assess the implications of foreground removal for power spectrum constraints where we find that our method can clean foregrounds well below the foreground wedge, rendering only scales $k_\parallel < 0.02 h \,\mathrm{Mpc}^{-1}$ inaccessible. As this approach assumes perfect knowledge of the telescope, we perform a conservative test of how essential this is by simulating and analysing datasets with deviations about our assumed telescope. Assuming no other techniques to mitigate bias are applied, we recover unbiased power spectra when the per-feed beam width to be measured to 0.1%, and amplifier gains to be known to 1% within each minute. Finally, as an example application, we extend our forecasts to a wideband 400-800 MHz cosmological observation and consider the implications for probing dark energy, finding a medium-sized cylinder telescope improves the DETF Figure of Merit by around 70% over Planck and Stage II experiments alone.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Method for Direct Measurement of Cosmic Acceleration by 21-cm Absorption Systems

So far there is only indirect evidence that the Universe is undergoing an accelerated expansion. The evidence for cosmic acceleration is based on the observation of different objects at different distances, and requires invoking the Copernican cosmological principle, and Einstein's equations of motion. We examine the direct observability using recession velocity drifts (Sandage-Loeb effect) of 21cm hydrogen absorption systems in upcoming radio surveys. This measures the change in velocity of the {\it same} objects separate by a time interval and is a model-independent measure of acceleration. We forecast that for a CHIME-like survey with a decade time span, we can detect the acceleration of a $Λ$CDM Universe with $5σ$ confidence. This acceleration test requires modest data analysis and storage changes from the normal processing, and cannot be recovered retroactively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Possible Astrophysical Observables of Quantum Gravity Effects near Black Holes

Recent implications of results from quantum information theory applied to black holes has led to the confusing conclusions that requires either abandoning the equivalence principle (e.g. the firewall picture), or the no-hair theorem (e.g. the fuzzball picture), or even more unpalatable options. The recent discovery of a pulsar orbiting a black hole opens up new possibilities for tests of theories of gravity. We examine possible observational effects of semiclassical quantum gravity in the vicinity of black holes, as probed by pulsars and event horizon telescope imaging of flares. Pulsar radiation is observable at wavelengths only two orders of magnitude shorter than the Hawking radiation, so precision interferometry of lensed pulsar images may shed light on the quantum gravitational processes and interaction of Hawking radiation with the spacetime near the black hole. This paper discusses the impact on the pulsar radiation interference pattern, which is observable through the modulation index in the foreseeable future, and discusses a possible classical limit of BHC.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pulsar scintillations from corrugated reconnection sheets in the ISM

We show that surface waves along interstellar current sheets closely aligned with the line of sight lead to pulsar scintillation properties consistent with those observed. This mechanism naturally produces the length and density scales of the ISM scattering lenses that are required to explain the magnitude and dynamical spectrum of the scintillations. In this picture, the parts of warm ionized interstellar medium that are responsible for the scintillations are relatively quiescent, with scintillation and scattering resulting from weak waves propagating along magnetic domain boundary current sheets, which are both expected from helicity conservation and have been observed in numerical simulations. The model statistically predicts the spacing and amplitudes of inverted parabolic arcs seen in Fourier-transformed dynamical spectra of strongly scintillating pulsars with only 3 parameters. Multi-frequency, multi-epoch low frequency VLBI observations can quantitatively test this picture. If successful, in addition to mapping the ISM, this may open the door to precise nanoarcsecond pulsar astrometry, distance measurements, and emission studies using these 10AU interferometers in the sky.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong lensing interferometry for compact binaries

We propose a possibility to improve the current precision measurements on compact binaries. When the orbital axis is almost perpendicular to our line of sight, a pulsar behind its companion can form two strong-lensing images. These images cannot be resolved, but we can use multi-wavelength interferometry to accurately determine the passage through superior conjunction. This method does not depend strongly on the stability of the pulse profile, and applies equally well to both slow and fast pulsars. We discuss the possible improvement this can bring to the bound on stochastic gravitational wave background and to determine black hole spin. We also discuss the possibility of discovering a suitable binary system by the Square Kilometer Array that our method can apply to.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

50 pico arcsecond astrometry of pulsar emission

We use VLBI imaging of the interstellar scattering speckle pattern associated with the pulsar PSR 0834+06 to measure the astrometric motion of its emission. The ~ 5AU interstellar baselines, provided by interference between speckles spanning the scattering disk, enable us to detect motions with sub nanoarcsecond accuracy. We measure a small pulse deflection of ~18+/-2 km (not including geometric uncertainties), which is 100 times smaller than the native resolution of this interstellar interferometer. This implies that the emission region is small, and at an altitude of a few hundred km, with the exact value depending on field geometry. This is substantially closer to the star than to the light cylinder. Future VLBI measurements can improve on this finding. This new regime of ultra-precise astrometry may enable precision parallax distance determination of pulsar binary displacements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A simulation calibrated limit on the HI power spectrum from the GMRT Epoch of Reionization experiment

The GMRT Epoch of Reionization (EoR) experiment is an ongoing effort to measure the power spectrum from neutral hydrogen at high redshift. We have previously reported an upper limit of (70 mK)^2 at wavenumbers of k=0.65 h/Mpc using a basic piecewise-linear foreground subtraction. In this paper we explore the use of a singular value decomposition to remove foregrounds with fewer assumptions about the foreground structure. Using this method we also quantify, for the first time, the signal loss due to the foreground filter and present new power spectra adjusted for this loss, providing a revised measurement of a 2-sigma upper limit at (248 mK)^2 for k=0.50 h/Mpc. While this revised limit is larger than previously reported, we believe it to be more robust and still represents the best current constraints on reionization at z=8.6.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

All-Sky Interferometry with Spherical Harmonic Transit Telescopes

In this paper we describe the spherical harmonic transit telescope, a novel formalism for the analysis of transit radio telescopes. This all-sky approach bypasses the curved sky complications of traditional interferometry and so is particularly well suited to the analysis of wide-field radio interferometers. It enables compact and computationally efficient representations of the data and its statistics that allow new ways of approaching important problems like map-making and foreground removal. In particular, we show how it enables the use of the Karhunen-Loeve transform as a highly effective foreground filter, suppressing realistic foreground residuals for our fiducial example by at least a factor twenty below the 21cm signal even in highly contaminated regions of the sky. This is despite the presence of the angle-frequency mode mixing inherent in real-world instruments with frequency-dependent beams. We show, using Fisher forecasting, that foreground cleaning has little effect on power spectrum constraints compared to hypothetical foreground-free measurements. Beyond providing a natural real-world data analysis framework for 21cm telescopes now under construction and future experiments, this formalism allows accurate power spectrum forecasts to be made that include the interplay of design constraints and realistic experimental systematics with twenty-first century 21cm science.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High Performance P3M N-body code: CUBEP3M

This paper presents CUBEP3M, a publicly-available high performance cosmological N-body code and describes many utilities and extensions that have been added to the standard package. These include a memory-light runtime SO halo finder, a non-Gaussian initial conditions generator, and a system of unique particle identification. CUBEP3M is fast, its accuracy is tuneable to optimize speed or memory, and has been run on more than 27,000 cores, achieving within a factor of two of ideal weak scaling even at this problem size. The code can be run in an extra-lean mode where the peak memory imprint for large runs is as low as 37 bytes per particles, which is almost two times leaner than other widely used N-body codes. However, load imbalances can increase this requirement by a factor of two, such that fast configurations with all the utilities enabled and load imbalances factored in require between 70 and 120 bytes per particles. CUBEP3M is well designed to study large scales cosmological systems, where imbalances are not too large and adaptive time-stepping not essential. It has already been used for a broad number of science applications that require either large samples of non-linear realizations or very large dark matter N-body simulations, including cosmological reionization, halo formation, baryonic acoustic oscillations, weak lensing or non-Gaussian statistics. We discuss the structure, the accuracy, known systematic effects and the scaling performance of the code and its utilities, when applicable.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-Gaussian error bars in galaxy surveys -- II

(Abridged) Estimating the uncertainty on the matter power spectrum internally (i.e. directly from the data) is made challenging by the simple fact that galaxy surveys offer at most a few independent samples. In addition, surveys have non-trivial geometries, which make the interpretation of the observations even trickier, but the uncertainty can nevertheless be worked out within the Gaussian approximation. With the recent realization that Gaussian treatments of the power spectrum lead to biased error bars about the dilation of the baryonic acoustic oscillation scale, efforts are being directed towards developing non-Gaussian analyses, mainly from N-body simulations so far. Unfortunately, there is currently no way to tell how the non-Gaussian features observed in the simulations compare to those of the real Universe, and it is generally hard to tell at what level of accuracy the N-body simulations can model complicated non-linear effects such as mode coupling and galaxy bias. We propose in this paper a novel method that aims at measuring non-Gaussian error bars on the matter power spectrum directly from galaxy survey data. We utilize known symmetries of the 4-point function, Wiener filtering and principal component analysis to estimate the full covariance matrix from only four independent fields with minimal prior assumptions. With the noise filtering techniques and only four fields, we are able to recover the Fisher information obtained from a large N=200 sample to within 20 per cent, for k < 1.0 h/Mpc. Finally, we provide a prescription to extract a noise-filtered, non-Gaussian, covariance matrix from a handful of fields in the presence of a survey selection function.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pulsar scintillation patterns and strangelets

We propose that interstellar extreme scattering events, usually observed as pulsar scintillations, may be caused by a coherent agent rather than the usually assumed turbulence of $\rm H_2$ clouds. We find that the penetration of a flux of ionizing, positively charged strangelets or quark nuggets into a dense interstellar hydrogen cloud may produce ionization trails. Depending on the specific nature and energy of the incoming droplets, diffusive propagation or even capture in the cloud are possible. As a result, enhanced electron densities may form and constitute a lens-like scattering screen for radio pulsars and possibly for quasars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pulsar timing arrays as imaging gravitational wave telescopes: angular resolution and source (de)confusion

Pulsar timing arrays (PTAs) will be sensitive to a finite number of gravitational wave (GW) "point" sources (e.g. supermassive black hole binaries). N quiet pulsars with accurately known distances d_{pulsar} can characterize up to 2N/7 distant chirping sources per frequency bin Δf_{gw}=1/T, and localize them with "diffraction limited" precision δθ\gtrsim (1/SNR)(λ_{gw}/d_{pulsar}). Even if the pulsar distances are poorly known, a PTA with F frequency bins can still characterize up to (2N/7)[1-(1/2F)] sources per bin, and the quasi-singular pattern of timing residuals in the vicinity of a GW source still allows the source to be localized quasi-topologically within roughly the smallest quadrilateral of quiet pulsars that encircles it on the sky, down to a limiting resolution δθ\gtrsim (1/SNR) \sqrt{λ_{gw}/d_{pulsar}}. PTAs may be unconfused, even at the lowest frequencies, with matched filtering always appropriate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulating cosmic reionization: How large a volume is large enough?

We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $Δz$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmic Tides

We apply CMB lensing techniques to large scale structure and solve for the 3-D cosmic tidal field. We use small scale filamentary structures to solve for the large scale tidal shear and gravitational potential. By comparing this to the redshift space density field, one can measure the gravitational growth factor on large scales without cosmic variance. This potentially enables accurate measurements of neutrino masses and reconstruction of radial modes lost in 21 cm intensity mapping, which are essential for CMB and other cross correlations. We relate the tidal fields to the squeezed limit bispectrum, and present initial results from simulations and data from the SDSS.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-Gaussian errors of baryonic acoustic oscillations

We revisit the uncertainty in baryon acoustic oscillation (BAO) forecasts and data analyses. In particular, we study how much the uncertainties on both the measured mean dilation scale and the associated error bar are affected by the non-Gaussianity of the non-linear density field. We examine two possible impacts of non-Gaussian analysis: (1) we derive the distance estimators from Gaussian theory, but use 1000 N-Body simulations to measure the actual errors, and compare this to the Gaussian prediction, and (2) we compute new optimal estimators, which requires the inverse of the non-Gaussian covariance matrix of the matter power spectrum. Obtaining an accurate and precise inversion is challenging, and we opted for a noise reduction technique applied on the covariance matrices. By measuring the bootstrap error on the inverted matrix, this work quantifies for the first time the significance of the non-Gaussian error corrections on the BAO dilation scale. We find that the variance (error squared) on distance measurements can deviate by up to 12% between both estimators, an effect that requires a large number of simulations to be resolved. We next apply a reconstruction algorithm to recover some of the BAO signal that had been smeared by non-linear evolution, and we rerun the analysis. We find that after reconstruction, the rms error on the distance measurement improves by a factor of ~1.7 at low redshift (consistent with previous results), and the variance (σ^2) shows a change of up to 18% between optimal and sub-optimal cases (note, however, that these discrepancies may depend in detail on the procedure used to isolate the BAO signal). We finally discuss the impact of this work on current data analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observational Consequences of Dark Energy Decay

We consider the generic scenario of dark energy which arises through the latent heat of a hidden sector first order cosmological phase transition. This field could account for the extra radiation degree of freedom suggested by the CMB. We present the bubble nucleation solution for the viscous limit. The decay rate of the field is constrained by published KSZ data, and may be an explanation of current excess ISW correlations. Cross correlation of current and future surveys can further constrain or test the parameter space. The decay model is plausibly in the observable range, and avoids anthropic problems. This class of models is not well constrained by the popular dark energy figure of merit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimizing the Recovery of Fisher Information in the Dark Matter Power Spectrum

We combine two Gaussianization techniques - Wavelet Non-Linear Wiener Filter (WNLWF) and density reconstruction - to quantify the recovery of Fisher information that is lost in the gravitational collapse. We compute a displacement fields, in analogy with the Zel'dovich approximation, and apply a Wavelet Non-Linear Wiener Filter that decomposes the reconstructed density fields into a Gaussian and a non-Gaussian component. From a series of 200 realizations of N-body simulations, we compute the recovery performance for density fields obtained with both dark matter particles and haloes. We find that the height of the Fisher information trans-linear plateau is increased by more than an order of magnitude at k > 1.0h/Mpc for particles, whereas either technique alone offers an individual recovery boost of only a factor of three to five. We conclude that these two techniques work in a symbiosis, as their combined performance is stronger than the sum of their individual contribution. When applied to the halo catalogues, we find that the reconstruction has only a weak effect on the recovery of Fisher Information, while the non-linear wavelet filter boosts the information by about a factor offive. We also observe that non-Gaussian Poisson noise saturates the Fisher information, and that shot noise subtracted measurements exhibit a milder information recovery.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulating Cosmic Reionization and the Radiation Backgrounds from the Epoch of Reionization

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Can 21-cm observations discriminate between high-mass and low-mass galaxies as reionization sources?

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Information Content in the Angular Power Spectrum of Weak Lensing: Wavelet Method

We quantify the performance of a non-linear Wiener filter, constructed in wavelet space, at recovering some of the Fisher information that was lost in the weak lensing convergence field. The proposed method consists in a separation of the original field into the sum of a Gaussian and a non-Gaussian contribution. After filtering an ensemble of such fields, which are obtained from $N$-body simulations, we find that we can recapture about four times more Fisher information, an effect that can potentially improve by a significant amount the constraining power of weak lensing surveys on cosmological parameters, including the dark energy equation of state $ω$. We compare this performance with that of the logarithmic mapping and find that the wavelet method can recover up to three times more information.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Gaussian Error in Galaxy Survey (Part 1)

We propose a method to estimate non-Gaussian error bars on the matter power spectrum from galaxy surveys in the presence of non-trivial survey selection functions. The estimators are often obtained from formalisms like FKP and PKL, which rely on the assumption that the underlying field is Gaussian. The Monte Carlo method is more accurate but involves the tedious process of running and cross-correlating a large number of N-body simulations, in which the survey volume is embedded. From 200 N-body simulations, we extract a non-linear covariance matrix as a function of two scales and of the angle between two Fourier modes. All the non-Gaussian features of that matrix are then simply parametrized in terms of a few fitting functions and Eigenvectors. We furthermore develop a fast and accurate strategy that combines our parameterization with a general galaxy survey selection function, and incor- porate non-Gaussian Poisson uncertainty. We describe how to incorporate these two distinct non-Gaussian contributions into a typical analysis pipeline, and apply our method with the se- lection function from the 2dFGRS. We find that the observed Fourier modes correlate at much larger scales than that predicted by both FKP formalism or by pure N-body simulation in a "top hat" selection function. In particular, the observed Fourier modes are already 50 per cent correlated at k\sim0.1hMpc^(-1), and the non-Gaussian fractional variance on the power spectrum (σ /P(k))^2 is about a factor of 3.0 larger than the FKP prescription. At k~0.4hMpc^(-1), the deviations are an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Refractive Convergent Plasma Lenses explain ESE and pulsar scintillation

We propose convergent plasma lenses, possibly from current sheets, as a generic solution to strong interstellar scattering. These lenses resolve the overpressure problem by geometric alignment as noted by Goldreich and Shridhar (2006). They further quantitatively explain properties of extreme scattering events, and pulsar parabolic arcs. This model makes quantitative predictions testable by VLBI on scattering events. It differs conceptually from previous models by occurring through rare, localized underdense sheets. Such sheets are thermally and kinematically stable, and could be consequences of reconnection. The apparent diffractive effects are a result of coherent interference of refractive images. We propose that these lenses can be used for precision distance determination to pulsars, enabling accurate gravity source localization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A real-time software backend for the GMRT

The new era of software signal processing has a large impact on radio astronomy instrumentation. Our design and implementation of a 32 antennae, 33 MHz, dual polarization, fully real-time software backend for the GMRT, using only off-the-shelf components, is an example of this. We have built a correlator and a beamformer, using PCI-based ADC cards and a Linux cluster of 48 nodes with dual gigabit inter-node connectivity for real-time data transfer requirements. The highly optimized compute pipeline uses cache efficient, multi-threaded parallel code, with the aid of vectorized processing. This backend allows flexibility in final time and frequency resolutions, and the ability to implement algorithms for radio frequency interference rejection. Our approach has allowed relatively rapid development of a fairly sophisticated and flexible backend receiver system for the GMRT, which will greatly enhance the productivity of the telescope. In this paper we describe some of the first lights using this software processing pipeline. We believe this is the first instance of such a real-time observatory backend for an intermediate sized array like the GMRT.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cosmic Structure Formation at High Redshift

We present some preliminary results from a series of extremely large, high-resolution N-body simulations of the formation of early nonlinear structures. We find that the high-z halo mass function is inconsistent with the Sheth-Tormen mass function, which tends to over-estimate the abundance of rare halos. This discrepancy is in rough agreement with previous results based on smaller simulations. We also show that the number density of minihaloes is correlated with local matter density, albeit with a significant scatter that increases with redshift, as minihaloes become increasingly rare. The average correlation is in rough agreement with a simple analytical extended Press-Schechter model, but can differ by up to factor of 2 in some regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detection of Large-Scale Cosmic Magnetic Fields

Observational evidence for the existence of cosmic magnetic fields in intergalactic space coherent over Mpc scales is presented. Using an unprecedentedly large sample of Faraday rotation measures of radio sources from the NRAO VLA Sky Survey data and the photometric redshift galaxy catalog from the sixth Data Release of the Sloan Digital Sky Survey (SDSS), we measure the cross-correlations between the rotation measures and the galaxy density field distributed along the source sightlines. It is shown that the rotation measures for sightlines passing through high density regions at separation r >= 1 Mpc/h from the locations of background radio sources are significantly enhanced. We discuss possible generators of this enhancement and interpret it to be intergalactic magnetic fields coherent over 1 Mpc/h with mean field strength B~30 nG.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Enhanced Detectability of Pre-reionization 21-cm Structure

Before the universe was reionized, it was likely that the spin temperature of intergalactic hydrogen was decoupled from the CMB by UV radiation from the first stars through the Wouthuysen-Field effect. If the IGM had not yet been heated above the CMB temperature by that time, then the gas would appear in absorption relative to the CMB. Large, rare sources of X-rays could inject sufficient heat into the neutral IGM, so that the differential brightness temperature was greater than zero at comoving distances of tens to hundreds of Mpc, resulting in large 21-cm fluctuations with amplitudes of about 250 mK on arcminute to degree angular scales, an order of magnitude larger in amplitude than that caused by ionized bubbles during reionization, about 25 mK. This signal could therefore be easier to detect and probe higher redshifts than that due to patchy reionization. For the case in which the first objects to heat the IGM are QSOs hosting 10^7-solar mass black holes with an abundance exceeding about 1 per Gpc^3 at z~15, observations with either the Arecibo Observatory or the Five Hundred Meter Aperture Spherical Telescope (FAST) could detect and image their fluctuations at greater than 5-sigma significance in about a month of dedicated survey time. Additionally, existing facilities such as MWA and LOFAR could detect the statistical fluctuations arising from a population of 10^5-solar mass black holes with an abundance of about 10^4 per Gpc^3 at z~10-12.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast magnetic reconnection in three dimensional MHD simulations

We present a constructive numerical example of fast magnetic reconnection in a three dimensional periodic box. Reconnection is initiated by a strong, localized perturbation to the field lines. The solution is intrinsically three dimensional, and its gross properties do not depend on the details of the simulations. $\sim 50%$ of the magnetic energy is released in an event which lasts about one Alfven time, but only after a delay during which the field lines evolve into a critical configuration. We present a physical picture of the process. The reconnection regions are dynamical and mutually interacting. In the comoving frame of these regions, reconnection occurs through an X-like point, analogous to Petschek reconnection. The dynamics appear to be driven by global flows, not local processes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hydrogen 21-cm Intensity Mapping at redshift 0.8

Observations of 21-cm radio emission by neutral hydrogen at redshifts z ~ 0.5 to ~ 2.5 are expected to provide a sensitive probe of cosmic dark energy. This is particularly true around the onset of acceleration at z ~ 1, where traditional optical cosmology becomes very difficult because of the infrared opacity of the atmosphere. Hitherto, 21-cm emission has been detected only to z=0.24. More distant galaxies generally are too faint for individual detections but it is possible to measure the aggregate emission from many unresolved galaxies in the 'cosmic web'. Here we report a three dimensional 21-cm intensity field at z=0.53 to 1.12. We then co-add HI emission from the volumes surrounding about ten thousand galaxies (from the DEEP2 optical galaxy redshift survey. We detect the aggregate 21-cm glow at a significance of ~ 4 sigma.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Improved Dark Energy Detection through the Polarization-assisted WMAP-NVSS ISW Correlation

Integrated Sachs-Wolfe (ISW) effect can be estimated by cross-correlating Cosmic Microwave Background (CMB) sky with tracers of the local matter distribution. At late cosmic time, the dark energy induced decay of gravitation potential generates a cross-correlation signal on large angular scales. The dominant noise are the intrinsic CMB anisotropies from the inflationary epoch. In this Letter we use CMB polarization to reduce this intrinsic noise. We cross-correlate the microwave sky observed by Wilkinson Microwave Anisotropy Probe (WMAP) with the radio source catalog compiled by NRAO VLA Sky Survey (NVSS) to study the efficiency of the noise suppression . We find that the error bars are reduced about 5-12 %, improving the statistical power.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Increasing the Fisher Information Content in the Matter Power Spectrum by Non-linear Wavelet Weiner Filtering

We develop a purely mathematical tool to recover some of the information lost in the non-linear collapse of large-scale structure. From a set of 141 simulations of dark matter density fields, we construct a non-linear Weiner filter in order to separate Gaussian and non-Gaussian structure in wavelet space. We find that the non-Gaussian power is dominant at smaller scales, as expected from the theory of structure formation, while the Gaussian counterpart is damped by an order of magnitude on small scales. We find that it is possible to increase the Fisher information by a factor of three before reaching the translinear plateau, an effect comparable to other techniques like the linear reconstruction of the density field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large-scale BAO signatures of the smallest galaxies

Recent work has shown that at high redshift, the relative velocity between dark matter and baryonic gas is typically supersonic. This relative velocity suppresses the formation of the earliest baryonic structures like minihalos, and the suppression is modulated on large scales. This effect imprints a characteristic shape in the clustering power spectrum of the earliest structures, with significant power on 100 Mpc scales featuring highly pronounced baryon acoustic oscillations. The amplitude of these oscillations is orders of magnitude larger at z=20 than previously expected. This characteristic signature can allow us to distinguish the effects of minihalos on intergalactic gas at times preceding and during reionization. We illustrate this effect with the example of 21 cm emission and absorption from redshifts during and before reionization. This effect can potentially allow us to probe physics on kpc scales using observations on 100 Mpc scales. We present sensitivity forecasts for FAST and Arecibo. Depending on parameters, this enhanced structure may be detectable by Arecibo at redshifts near z=15-20, and with appropriate instrumentation FAST could measure the BAO power spectrum with high precision. In principle, this effect could also pose a serious challenge for efforts to constrain dark energy using observations of the BAO feature at low redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetohydrodynamics on Heterogeneous architectures: a performance comparison

We present magneto-hydrodynamic simulation results for heterogeneous systems. Heterogeneous architectures combine high floating point performance many-core units hosted in conventional server nodes. Examples include Graphics Processing Units (GPU's) and Cell. They have potentially large gains in performance, at modest power and monetary cost. We implemented a magneto-hydrodynamic (MHD) simulation code on a variety of heterogeneous and multi-core architectures --- multi-core x86, Cell, Nvidia and ATI GPU --- in different languages, FORTRAN, C, Cell, CUDA and OpenCL. We present initial performance results for these systems. To our knowledge, this is the widest comparison of heterogeneous systems for MHD simulations. We review the different challenges faced in each architecture, and potential bottlenecks. We conclude that substantial gains in performance over traditional systems are possible, and in particular that is possible to extract a greater percentage of peak theoretical performance from some systems when compared to x86 architectures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Near term measurements with 21 cm intensity mapping: neutral hydrogen fraction and BAO at z<2

It is shown that 21 cm intensity mapping could be used in the near term to make cosmologically useful measurements. Large scale structure could be detected using existing radio telescopes, or using prototypes for dedicated redshift survey telescopes. This would provide a measure of the mean neutral hydrogen density, using redshift space distortions to break the degeneracy with the linear bias. We find that with only 200 hours of observing time on the Green Bank Telescope, the neutral hydrogen density could be measured to 25% precision at redshift 0.54<z<1.09. This compares favourably to current measurements, uses independent techniques, and would settle the controversy over an important parameter which impacts galaxy formation studies. In addition, a 4000 hour survey would allow for the detection of baryon acoustic oscillations, giving a cosmological distance measure at 3.5% precision. These observation time requirements could be greatly reduced with the construction of multiple pixel receivers. Similar results are possible using prototypes for dedicated cylindrical telescopes on month time scales, or SKA pathfinder aperture arrays on day time scales. Such measurements promise to improve our understanding of these quantities while beating a path for future generations of hydrogen surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Parameter Survey of Nonradiative Black Hole Accretion -- Flow Structure and Variability of the Rotation Measure

We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars. Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow. We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Primordial gravity wave fossils and their use in testing inflation

A new effect is described by which primordial gravity waves leave a permanent signature in the large scale structure of the Universe. The effect occurs at second order in perturbation theory and is sensitive to the order in which perturbations on different scales are generated. We derive general forecasts for the detectability of the effect with future experiments, and consider observations of the pre-reionization gas through the 21 cm line. It is found that the Square Kilometre Array will not be competitive with current cosmic microwave background constraints on primordial gravity waves from inflation. However, a more futuristic experiment could, through this effect, provide the highest ultimate sensitivity to tensor modes and possibly even measure the tensor spectral index. It is thus a potentially quantitative probe of the inflationary paradigm.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Projected Constraints on Modified Gravity Cosmologies from 21 cm Intensity Mapping

We present projected constraints on modified gravity models from the observational technique known as 21 cm intensity mapping, where cosmic structure is detected without resolving individual galaxies. The resulting map is sensitive to both BAO and weak lensing, two of the most powerful cosmological probes. It is found that a 200 m x 200 m cylindrical telescope, sensitive out to z=2.5, would be able to distinguish DGP from most dark energy models, and constrain the Hu & Sawicki f(R) model to |f_{R0}| < 9*10^(-6) at 95% confidence. The latter constraint makes extensive use of the lensing spectrum in the nonlinear regime. These results show that 21 cm intensity mapping is not only sensitive to modifications of the standard model's expansion history, but also to structure growth. This makes intensity mapping a powerful and economical technique, achievable on much shorter time scales than optical experiments that would probe the same era.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self calibration of photometric redshift scatter in weak lensing surveys

Photo-z errors, especially catastrophic errors, are a major uncertainty for precision weak lensing cosmology. We find that the shear-(galaxy number) density and density-density cross correlation measurements between photo-z bins, available from the same lensing surveys, contain valuable information for self-calibration of the scattering probabilities between the true-z and photo-z bins. The self-calibration technique we propose does not rely on cosmological priors nor parameterization of the photo-z probability distribution function, and preserves all of the cosmological information available from shear-shear measurement. We estimate the calibration accuracy through the Fisher matrix formalism. We find that, for advanced lensing surveys such as the planned stage IV surveys, the rate of photo-z outliers can be determined with statistical uncertainties of 0.01-1% for $z<2$ galaxies. Among the several sources of calibration error that we identify and investigate, the {\it galaxy distribution bias} is likely the most dominant systematic error, whereby photo-z outliers have different redshift distributions and/or bias than non-outliers from the same bin. This bias affects all photo-z calibration techniques based on correlation measurements. Galaxy bias variations of $O(0.1)$ produce biases in photo-z outlier rates similar to the statistical errors of our method, so this galaxy distribution bias may bias the reconstructed scatters at several-$σ$ level, but is unlikely to completely invalidate the self-calibration technique.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The GMRT-EoR Experiment: A new upper limit on the neutral hydrogen power spectrum at z~8.6

We present a new upper limit to the 21cm power spectrum during the Epoch of Reionization (EoR) which constrains reionization models with an unheated IGM. The GMRT-EoR experiment is an ongoing effort to make a statistical detection of the power spectrum of 21cm neutral hydrogen emission at redshift z~9. Data from this redshift constrain models of the EoR, the end of the Dark Ages arising from the formation of the first bright UV sources, probably stars or mini-quasars. We present results from approximately 50 hours of observations at the Giant Metrewave Radio Telescope in India from December 2007. We describe radio frequency interference (RFI) localisation schemes which allow bright sources on the ground to be identified and physically removed in addition to automated flagging. Singular-value decomposition is used to remove remaining broadband RFI by identifying ground sources with large eigenvalues. Foregrounds are modelled using a piecewise linear filter and the power spectrum is measured using cross-correlations of foreground subtracted images.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dark Energy from Large-Scale Structure Lensing Information

Wide area Large-Scale Structure (LSS) surveys are planning to map a substantial fraction of the visible universe to quantify dark energy through Baryon Acoustic Oscillations (BAO). At increasing redshift, for example that probed by proposed 21-cm intensity mapping surveys, gravitational lensing potentially limits the fidelity (Hui et al., 2007) because it distorts the apparent matter distribution. In this paper we show that these distortions can be reconstructed, and actually used to map the distribution of intervening dark matter. The lensing information for sources at z=1-3 allows accurate reconstruction of the gravitational potential on large scales, l <~ 100, which is well matched for Integrated Sachs-Wolfe (ISW) effect measurements of dark energy and its sound speed, and a strong constraint for modified gravity models of dark energy. We built an optimal quadratic lensing estimator for non-Gaussian sources, which is necessary for LSS. The phenomenon of "information saturation" (Rimes & Hamilton, 2005) saturates reconstruction at mildly non-linear scales, where the linear source power spectrum Delta^2 ~ 0.2-0.5. We find that steeper power spectra, saturate more quickly. We compute the effective number densities of independent lensing sources for LSS lensing, and find that they increase rapidly with redshifts. For LSS/21-cm sources at z ~ 2-4, the lensing reconstruction is limited by cosmic variance at l<~ 100.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The GMRT EoR Experiment: Limits on Polarized Sky Brightness at 150 MHz

The GMRT reionization effort aims to map out the large scale structure of the Universe during the epoch of reionization (EoR). Removal of polarized Galactic emission is a difficult part of any 21 cm EoR program, and we present new upper limits to diffuse polarized foregrounds at 150 MHz. We find no high significance evidence of polarized emission in our observed field at mid galactic latitude (J2000 08h26m+26). We find an upper limit on the 2-dimensional angular power spectrum of diffuse polarized foregrounds of [l^2 C_l/(2 PI)]^{1/2}< 3K in frequency bins of width 1 MHz at 300<l<1000. The 3-dimensional power spectrum of polarized emission, which is most directly relevant to EoR observations, is [k^3 P_p(k)/(2 PI^2)]^{1/2}< 2K at k_perp > 0.03 h/Mpc, k < 0.1 h/Mpc. This can be compared to the expected EoR signal in total intensity of [k^3 P(k)/ (2 PI^2) ]^{1/2} ~ 10 mK. We find polarized structure is substantially weaker than suggested by extrapolation from higher frequency observations, so the new low upper limits reported here reduce the anticipated impact of these foregrounds on EoR experiments. We discuss Faraday beam and depth depolarization models and compare predictions of these models to our data. We report on a new technique for polarization calibration using pulsars, as well as a new technique to remove broadband radio frequency interference. Our data indicate that, on the edges of the main beam at GMRT, polarization squint creates ~ 3% leakage of unpolarized power into polarized maps at zero rotation measure. Ionospheric rotation was largely stable during these solar minimum night time observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Inhomogeneous Background of Hydrogen-Molecule Dissociating Radiation during Cosmic Reionization

The first, self-consistent calculations are presented of the cosmological, H2-dissociating UV background produced during the epoch of reionization by the sources of reionization. Large-scale radiative transfer simulations of reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down to dwarf galaxies of mass ~10^8 Msun, identified by very high-resolution N-body simulations, including the self-regulating effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the inhomogeneous radiation background in the Lyman-Werner bands of H2 between 11 and 13.6 eV. On average, the intensity of this Lyman-Werner background is found to rise to the threshold level at which dissociation suppresses H2 cooling and star formation inside minihalos, long before reionization is complete. Spatial variations in the Lyman-Werner background are found which result from the clustering of sources associated with large-scale structure formation, such that intensity fluctuations correlate with matter density fluctuations. As a result, the Lyman-Werner background rises to the threshold level for H2 suppression earlier in the vicinity of the reionization sources and their H II regions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Yuan-Tseh Lee Array for Microwave Background Anisotropy

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Halo stochasticity in global clustering analysis

In the present work we study the statistics of haloes, which in the halo model determines the distribution of galaxies. Haloes are known to be biased tracer of dark matter, and at large scales it is usually assumed there is no intrinsic stochasticity between the two fields. Following the work of Seljak & Warren (2004), we explore how correct this assumption is and, moving a step further, we try to qualify the nature of stochasticity. We use Principal Component Analysis applied to the outputs of a cosmological N-body simulation to: (1) explore the behaviour of stochasticity in the correlation between haloes of different masses; (2) explore the behaviour of stochasticity in the correlation between haloes and dark matter. We show results obtained using a catalogue with 2.1 million haloes, from a PMFAST simulation with box size of 1000h^{-1}Mpc. In the relation between different populations of haloes we find that stochasticity is not-negligible even at large scales. In agreement with the conclusions of Tegmark & Bromley (1999) who studied the correlations of different galaxy populations, we found that the shot-noise subtracted stochasticity is qualitatively different from `enhanced' shot noise and, specifically, it is dominated by a single stochastic eigenvalue. We call this the `minimally stochastic' scenario, as opposed to shot noise which is `maximally stochastic'. In the correlation between haloes and dark matter, we find that stochasticity is minimized, as expected, near the dark matter peak (k ~ 0.02 h Mpc^{-1} for a LambdaCDM cosmology) and, even at large scales, it is of the order of 15 per cent above the shot noise. Moreover, we find that the reconstruction of the dark matter distribution is improved when we use eigenvectors as tracers of the bias. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Quantifying the Cosmic Web I: The large-scale halo ellipticity-ellipticity and ellipticity-direction correlations

The formation of dark matter halos tends to occur anisotropically along the filaments of the Cosmic Web, which induces both ellipticity-ellipticity (EE) correlations between the shapes of halos, as well as ellipticity-direction (ED) cross-correlations between halo shapes and the directions to neighboring halos. We analyze the halo catalogue and the semi-analytic galaxy catalogue of the recent Millennium Run Simulation to measure the EE and ED correlations numerically at four different redshifts (z=0, 0.5, 1 and 2). For the EE correlations, we find that (i) the major-axis correlation is strongest while the intermediate-axis correlation is weakest; (ii) the signal is significant at distances out to 10 Mpc/h; (iii) the signal decreases as z decreases; (iv) and its behavior depends strongly on the halo mass scale, with larger masses showing stronger correlations at large distances. For the ED correlations, we find that (i) the correlations are much stronger than the EE correlations, and are significant even out to distances of 50 Mpc/h; (ii) the signal also decreases as z decreases; (iii) and it increases with halo mass at all distances. We also provide empirical fitting functions for the EE and ED correlations. The EE correlations are found to scale linearly with the linear density correlation function, xi(r). While the ED cross-correlation is found to scale as xi^{1/2}(r) at large distances beyond 10 Mpc/h. The best-fit values of the fitting parameters for the EE and the ED correlations are all determined through chi^{2}-statistics. Our results may be useful for quantifying the filamentary distribution of dark matter halos over a wide range of scales.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Kinetic Sunyaev-Zel'dovich Effect from Patchy Reionization: the View from the Simulations

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneus reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l(l+1)C_l/(2pi)_{max}~4-7x10^{-13}. The scale roughly corresponds to the typical ionized bubble sizes observed in our simulations, of ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. At large scales the patchy kSZ signal depends only on the source efficiencies. It is higher when sources are more efficient at producing ionizing photons, since such sources produce larger ionized regions, on average, than less efficient sources. The introduction of sub-grid gas clumping in the radiative transfer simulations produce significantly more power at small scales, but has little effect at large scales. The patchy reionization kSZ signal is dominated by the post-reionization signal from fully-ionized gas, but the two contributions are of similar order at scales l~3000-10^4, indicating that the kSZ anisotropies from reionization are an important component of the total kSZ signal at these scales.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Kinetic Sunyaev-Zel'dovich Effect from Radiative Transfer Simulations of Patchy Reionization

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneous reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l^2C_l~1\times10^{-12}. The peak scale is determined by the typical size of the ionized regions and roughly corresponds to the ionized bubble sizes observed in our simulations, ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. This predicted kSZ signal at arcminute scales is sufficiently strong to be detectable by upcoming experiments, like the Atacama Cosmology Telescope and South Pole Telescope which are expected to have ~1' resolution and ~muK sensitivity. The extended and patchy nature of the reionization process results in a boost of the peak signal in power by approximately one order of magnitude compared to a uniform reionization scenario, while roughly tripling the signal compared with that based upon the assumption of gradual but spatially uniform reionization. At large scales the patchy kSZ signal depends largely on the ionizing source efficiencies and the large-scale velocity fields: sources which produce photons more efficiently yield correspondingly higher signals. The introduction of sub-grid gas clumping in the radiative transfer simulations produces significantly more power at small scales, and more non-Gaussian features, but has little effect at large scales. The patchy nature of the reionization process roughly doubles the total observed kSZ signal for l~3000-10^4 compared to non-patchy scenarios with the same total electron-scattering optical depth.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Rotation in Gravitational Lenses

Gravitational lensing deflects light. A single lens deflector can only shear images, but cannot induce rotations. Multiple lens planes can induce rotations. Such rotations can be observed in quadruply imaged sources, and can be used to distinguish between two proposed solutions of the flux anomaly problem: substructures in lensing galaxies vs large scale structure. We predict the expected amount of rotation due to large scale structure in strong lensing systems, and show how this effect can be measured using ~ mas VLBI astrometry of quadruple lenses with extended source structures. The magnitude of rotation is around one degree. The biggest theoretical uncertainty is the power spectrum of dark matter on very small scales. This procedure can potentially be turned around to measure the dark matter power spectrum on very small scales. We list the predicted RMS rotation angles for several quadruple lenses with known lens and source redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

A Free, Fast, Simple and Efficient TVD MHD Code

We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin. The magnetic field is updated separately in two-dimensional advection-constraint steps. The electromotive force (EMF) is computed in the advection step using the TVD scheme, and this same EMF is used immediately in the constraint step in order to preserve \grad.B=0 without the need to store intermediate fluxes. Operator splitting is used to extend the code to three dimensions, and Runge-Kutta is used to get second order accuracy in time. The advantages of this code are high resolution per grid cell, second order accuracy in space and time, enforcement of the \grad.B=0 constraint to machine precision, no memory overhead, speed, and simplicity. A 3-D Fortran implementation less than 400 lines long is made freely available. We also implemented a fully scalable message-passing parallel MPI version. We present tests of the code on MHD waves and shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Hydrodynamical Simulations of the IGM at High Mach Numbers

We present a new approach to doing Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in hydrodynamical simulations of the IGM. In conventional Eulerian CFD, the thermal energy is poorly tracked in supersonic bulk flows where local fluid variables cannot be accurately separated from the much larger bulk flow components. We described a method in which local fluid quantities can be directly tracked and the Eulerian fluid equations solved in a local frame moving with the flow. The new algorithm has been used to run large hydrodynamical simulations on a 1024^3 grid to study the kinetic SZ effect. The KSZ power spectrum is broadly peaked at l~10^4 with temperature fluctuations on micro Kelvin levels.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Discriminating weak lensing from intrinsic spin correlations using the curl-gradient decomposition

The distortion field defined by the ellipticities of galaxy shapes projected on the sky can be uniquely decomposed into a gradient and a curl component. If the observed ellipticities are induced by weak gravitational lensing, then the distortion field is curl free. Here we show that, in contrast, the distortion field resulting from intrinsic spin alignments is not curl free. This provides a powerful discriminant between lensing and intrinsic contributions to observed ellipticity correlations. We also show how these contributions can be disentangled statistically from the ellipticity correlations or computed locally from circular integrals of the ellipticity field. This allows for an unambiguous detection of intrinsic galaxy alignments in the data. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors. These techniques can be applied equally well to the polarization of the microwave background, where it can be used to separate curl-free scalar perturbations from those produced by gravity waves or defects.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Brief Note: Analytical Fit to the Luminosity Distance for Flat Cosmologies with a Cosmological Constant

We present a fitting formula for the luminosity and angular diameter distances in cosmological models with pressureless matter, a cosmological constant and zero spatial curvature. The formula has a relative error of less than 0.4% for 0.2<Ω_0<1 for any redshift, and a global relative error of less than 4% for any choice of parameters.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Cosmic Microwave Anisotropies from Topological Defects in an Open Universe

We present a general formalism for computing Cosmic Background Radiation (CBR) and density fluctuations in open models with stiff sources. We find analytic Green's functions for the linearized Einstein equations in the presence of stiff sources and use this formalism to estimate the amplitude and harmonic spectrum of microwave background fluctuations produced by topological defects in an open universe. Unlike inflationary models that predict a flat universe and a spectrum of CBR fluctuations that is enhanced at large angular scales, defect models predict that CBR fluctuations are suppressed on angular scales larger than that subtended by the curvature scale. In an $Ω= 0.2 - 0.4$ universe, these models, when normalized to the amplitude of CBR fluctuations observed by COBE, require a moderate bias factor, $2-3$, to be compatible with the observed fluctuations in galaxy counts. In these models, accurate predictions can be made which are testable through CBR experiments in the near future. A CBR measurement of $Ω$ would then be possible, up to the limit imposed by cosmic variance. We discuss some of the philosophical implications of an open model and propose a solution to the flatness problem.

0 citations0 reviewsNo code linkedOpen access