Source author record

Christopher Hirata

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

14 published item(s)

preprint2022arXiv

Effect of dust in circumgalactic haloes on the cosmic shear power spectrum

Weak gravitational lensing is a powerful statistical tool for probing the growth of cosmic structure and measuring cosmological parameters. However, as shown by studies such as Ménard et al. (2010), dust in the circumgalactic region of haloes dims and reddens background sources. In a weak lensing analysis, this selects against sources behind overdense regions; since there is more structure in overdense regions, we will underestimate the amplitude of density perturbations $σ_8$ if we do not correct for the effects of circumgalactic dust. To model the dust distribution we employ the halo model. Assuming a fiducial dust mass profile based on measurements from Ménard et al. (2010), we compute the ratio $Z$ of the systematic error to the statistical error for a survey similar to the Nancy Grace Roman Space Telescope reference survey (2000 deg$^2$ area, single-filter effective source density 30 galaxies arcmin$^{-2}$). For a waveband centered at $1580$ nm ($H$-band), we find that $Z_{H} = 0.37$. For a similar survey with waveband centered at $620$ nm ($r$-band), we also computed $Z_{r} = 2.8$. Within our fiducial dust model, since $Z_{r} > 1$, the systematic effect of dust will be significant on weak lensing image surveys. We also computed the dust bias on the amplitude of the power spectrum, $σ_{8}$, and found it to be for each waveband $Δσ_8/σ_8 = -3.1\times 10^{-4}$ ($H$ band) or $-2.2\times 10^{-3}$ ($r$ band) if all other parameters are held fixed (the forecast Roman statistical-only error $σ(σ_8)/σ_8$ is $9\times 10^{-4}$).

preprint2021arXiv

Superresolution Reconstruction of Severely Undersampled Point-spread Functions Using Point-source Stacking and Deconvolution

Point-spread function (PSF) estimation in spatially undersampled images is challenging because large pixels average fine-scale spatial information. This is problematic when fine-resolution details are necessary, as in optimal photometry where knowledge of the illumination pattern beyond the native spatial resolution of the image may be required. Here, we introduce a method of PSF reconstruction where point sources are artificially sampled beyond the native resolution of an image and combined together via stacking to return a finely sampled estimate of the PSF. This estimate is then deconvolved from the pixel-gridding function to return a superresolution kernel that can be used for optimally weighted photometry. We benchmark against the < 1% photometric error requirement of the upcoming SPHEREx mission to assess performance in a concrete example. We find that standard methods like Richardson--Lucy deconvolution are not sufficient to achieve this stringent requirement. We investigate a more advanced method with significant heritage in image analysis called iterative back-projection (IBP) and demonstrate it using idealized Gaussian cases and simulated SPHEREx images. In testing this method on real images recorded by the LORRI instrument on New Horizons, we are able to identify systematic pointing drift. Our IBP-derived PSF kernels allow photometric accuracy significantly better than the requirement in individual SPHEREx exposures. This PSF reconstruction method is broadly applicable to a variety of problems and combines computationally simple techniques in a way that is robust to complicating factors such as severe undersampling, spatially complex PSFs, noise, crowded fields, or limited source numbers.

preprint2015arXiv

Sloan Digital Sky Survey III Photometric Quasar Clustering: Probing the Initial Conditions of the Universe using the Largest Volume

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky, and delivered over a trillion pixels of imaging data. We present the large-scale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that have been classified from this imaging, representing the highest density of quasars ever studied for clustering measurements. This data set spans ~11,000 square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large volume and medium density of tracers should facilitate high-precision cosmological constraints. We measure the angular clustering of photometrically classified quasars using an optimal quadratic estimator in four redshift slices with an accuracy of ~25% over a bin width of l ~10 - 15 on scales corresponding to matter-radiation equality and larger (l ~ 2 - 30). Observational systematics can strongly bias clustering measurements on large scales, which can mimic cosmologically relevant signals such as deviations from Gaussianity in the spectrum of primordial perturbations. We account for systematics by employing a new method recently proposed by Agarwal et al. (2014) to the clustering of photometrically classified quasars. We carefully apply our methodology to mitigate known observational systematics and further remove angular bins that are contaminated by unknown systematics. Combining quasar data with the photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et al. (2012), and marginalizing over all bias and shot noise-like parameters, we obtain a constraint on local primordial non-Gaussianity of fNL = -113+/-154 (1σerror). [Abridged]

preprint2015arXiv

Stability of small-scale baryon perturbations during cosmological recombination

In this paper, we study small-scale fluctuations (baryon pressure sound waves) in the baryon fluid during recombination. In particular, we look at their evolution in the presence of relative velocities between baryons and photons on large scales ($k \sim 10^{-1} \ {\rm Mpc}^{-1}$), which are naturally present during the era of decoupling. Previous work concluded that the fluctuations grow due to an instability of sound waves in a recombining plasma, but that the growth factor is small for typical cosmological models. These analyses model recombination in an inhomogenous universe as a perturbation to the parameters of the homogenous solution. We show that for relevant wavenumbers $k\gtrsim 10^3\ {\rm Mpc}^{-1}$ the dynamics are significantly altered by the transport of both ionizing continuum ($hν>13.6$ eV) and Lyman-$α$ photons between crests and troughs of the density perturbations. We solve the radiative transfer of photons in both these frequency ranges and incorporate the results in a perturbed three-level atom model. We conclude that the instability persists at intermediate scales. We use the results to estimate a distribution of growth rates in $10^{7}$ random realizations of large-scale relative velocities. Our results indicate that there is no appreciable growth; out of these $10^7$ realizations, the maximum growth factor we find is less than $\approx 1.2$ at wavenumbers of $k \approx 10^{3} \ {\rm Mpc}^{-1}$. The instability's low growth factors are due to the relatively short duration of the recombination epoch during which the electrons and photons are coupled.

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.

preprint2013arXiv

Observational Probes of Cosmic Acceleration

The accelerating expansion of the universe is the most surprising cosmological discovery in many decades, implying that the universe is dominated by some form of "dark energy" with exotic physical properties, or that Einstein's theory of gravity breaks down on cosmological scales. The profound implications of cosmic acceleration have inspired ambitious experimental efforts to measure the history of expansion and growth of structure with percent-level precision or higher. We review in detail the four most well established methods for making such measurements: Type Ia supernovae, baryon acoustic oscillations (BAO), weak gravitational lensing, and galaxy clusters. We pay particular attention to the systematic uncertainties in these techniques and to strategies for controlling them at the level needed to exploit "Stage IV" dark energy facilities such as BigBOSS, LSST, Euclid, and WFIRST. We briefly review a number of other approaches including redshift-space distortions, the Alcock-Paczynski test, and direct measurements of H_0. We present extensive forecasts for constraints on the dark energy equation of state and parameterized deviations from GR, achievable with Stage III and Stage IV experimental programs that incorporate supernovae, BAO, weak lensing, and CMB data. We also show the level of precision required for other methods to provide constraints competitive with those of these fiducial programs. We emphasize the value of a balanced program that employs several of the most powerful methods in combination, both to cross-check systematic uncertainties and to take advantage of complementary information. Surveys to probe cosmic acceleration produce data sets with broad applications, and they continue the longstanding astronomical tradition of mapping the universe in ever greater detail over ever larger scales.

preprint2012arXiv

Exploring the NRO Opportunity for a Hubble-sized Wide-field Near-IR Space Telescope -- NEW WFIRST

We discuss scientific, technical and programmatic issues related to the use of an NRO 2.4m telescope for the WFIRST initiative of the 2010 Decadal Survey. We show that this implementation of WFIRST, which we call "NEW WFIRST," would achieve the goals of the NWNH Decadal Survey for the WFIRST core programs of Dark Energy and Microlensing Planet Finding, with the crucial benefit of deeper and/or wider near-IR surveys for GO science and a potentially Hubble-like Guest Observer program. NEW WFIRST could also include a coronagraphic imager for direct detection of dust disks and planets around neighboring stars, a high-priority science and technology precursor for future ambitious programs to image Earth-like planets around neighboring stars.

preprint2012arXiv

The signature of the first stars in atomic hydrogen at redshift 20

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.

preprint2011arXiv

Optimal Linear Image Combination

A simple, yet general, formalism for the optimized linear combination of astrophysical images is constructed and demonstrated. The formalism allows the user to combine multiple undersampled images to provide oversampled output at high precision. The proposed method is general and may be used for any configuration of input pixels and point spread function; it also provides the noise covariance in the output image along with a powerful metric for describing undesired distortion to the image convolution kernel. The method explicitly provides knowledge and control of the inevitable compromise between noise and fidelity in the output image. We present a first prototype implementation of the method, outlining steps taken to generate an efficient algorithm. This implementation is then put to practical use in reconstructing fully-sampled output images using simulated, undersampled input exposures that are designed to mimic the proposed \emph{Wide Field InfraRed Survey Telescope} (\emph{WFIRST}). We examine results using randomly rotated and dithered input images, while also assessing better-known "ideal" dither patterns: comparing results we illustrate the use of the method as a survey design tool. Finally, we use the method to test the robustness of linear image combination when subject to practical realities such as missing input pixels and focal plane plate scale variations.

preprint2010arXiv

Non-Gaussianity and large-scale structure in a two-field inflationary model

Single field inflationary models predict nearly Gaussian initial conditions and hence a detection of non-Gaussianity would be a signature of the more complex inflationary scenarios. In this paper we study the effect on the cosmic microwave background and on large scale structure from primordial non-Gaussianity in a two-field inflationary model in which both the inflaton and curvaton contribute to the density perturbations. We show that in addition to the previously described enhancement of the galaxy bias on large scales, this setup results in large-scale stochasticity. We provide joint constraints on the local non-Gaussianity parameter $\tilde f_{\rm NL}$ and the ratio $ξ$ of the amplitude of primordial perturbations due to the inflaton and curvaton using WMAP and SDSS data.

preprint2010arXiv

Relative velocity of dark matter and baryonic fluids and the formation of the first structures

At the time of recombination, baryons and photons decoupled and the sound speed in the baryonic fluid dropped from relativistic to the thermal velocities of the hydrogen atoms. This is less than the relative velocities of baryons and dark matter computed via linear perturbation theory, so we infer that there are supersonic coherent flows of the baryons relative to the underlying potential wells created by the dark matter. As a result, the advection of small-scale perturbations (near the baryonic Jeans scale) by large-scale velocity flows is important for the formation of the first baryonic structures. This effect involves a quadratic term in the cosmological perturbation theory equations and hence has not been included in studies based on linear perturbation theory. We show that the relative motion suppresses the abundance of the first bound objects, even if one only investigates dark matter haloes, and leads to qualitative changes in their spatial distribution, such as introducing scale-dependent bias and stochasticity. We discuss the possible observable implications for high-redshift galaxy clustering and reionization.

preprint2010arXiv

Suppression and Spatial Variation of Early Galaxies and Minihalos

We study the effect of the relative velocity of dark matter and baryonic fluids after the epoch of recombination on the evolution of the first bound objects in the early universe. Recent work has shown that, although relative motion of the two fluids is formally a second order effect in density, it has a dramatic impact on the formation and distribution of the first cosmic structures. Focusing on the gas content, we analyze the effect of relative velocity on the properties of halos over a wide range of halo masses and redshifts. We calculate accurately the linear evolution of the baryon and dark matter fluctuations, and quantify the resulting effect on halos based on an analytical formalism that has been carefully checked with simulations in the case with no relative velocity. We estimate the effect on the abundance of and gas fraction in early halos. We find that the relative velocity effect causes several changes: (i) the characteristic mass that divides gas-rich and gas-poor halos is increased by roughly an order of magnitude, from 2 10^4 Msun to about 2 10^5 Msun; (ii) this characteristic mass has a large scatter (full width at half maximum is ~ 1.5 10^5 Msun at z=20); (iii) the fraction of baryons in star-less gas minihalos is suppressed by a factor of 4 at z=20; (iv) the fraction of baryons in halos that can cool and form stars is suppressed by a factor of 1.5 at z=20; and (v) there are enhanced spatial variations of these various fractions.

preprint2010arXiv

Tidal alignments as a contaminant of the galaxy bispectrum

If the orientations of galaxies are correlated with large-scale structure, then anisotropic selection effects such as preferential selection of face-on disc galaxies can contaminate large scale structure observables. Here we consider the effect on the galaxy bispectrum, which has attracted interest as a way to break the degeneracy between galaxy bias and the amplitude of matter fluctuations sigma_8. We consider two models of intrinsic galaxy alignments: one where the probability distribution for the galaxy's orientation contains a term linear in the local tidal field, appropriate for elliptical galaxies; and one with a term quadratic in the local tidal field, which may be applicable to disc galaxies. We compute the correction to the redshift-space bispectrum in the quasilinear regime, and then focus on its effects on parameter constraints from the transverse bispectrum, i.e. using triangles in the plane of the sky. We show that in the linear alignment model, intrinsic alignments result in an error in the galaxy bias parameters, but do not affect the inferred value of sigma_8. In contrast, the quadratic alignment model results in a systematic error in both the bias parameters and sigma_8. However, the quadratic alignment effect has a unique configuration dependence that should enable it to be removed in upcoming surveys.

preprint2008arXiv

Constraints on local primordial non-Gaussianity from large scale structure

Recent work has shown that the local non-Gaussianity parameter f_NL induces a scale-dependent bias, whose amplitude is growing with scale. Here we first rederive this result within the context of peak-background split formalism and show that it only depends on the assumption of universality of mass function, assuming halo bias only depends on mass. We then use extended Press-Schechter formalism to argue that this assumption may be violated and the scale dependent bias will depend on other properties, such as merging history of halos. In particular, in the limit of recent mergers we find the effect is suppressed. Next we use these predictions in conjunction with a compendium of large scale data to put a limit on the value of f_NL. When combining all data assuming that halo occupation depends only on halo mass, we get a limit of -29 ~ (-65)< f_NL < +70 ~(+93) at 95% (99.7%) confidence. While we use a wide range of datasets, our combined result is dominated by the signal from the SDSS photometric quasar sample. If the latter are modeled as recent mergers then the limits weaken to -31 ~(-96) < f_NL < +70 ~ (+96) . These limits are comparable to the strongest current limits from the WMAP 5 year analysis, with no evidence of a positive signal in f_NL. While the method needs to be thoroughly tested against large scale structure simulations with realistic quasar and galaxy formation models, our results indicate that this is a competitive method relative to CMB and should be further pursued both observationally and theoretically.