Showing 1–50 of 142 results for author: Hirata, C M

Recommendations for Early Definition Science with the Nancy Grace Roman Space Telescope

Authors: Robyn E. Sanderson, Ryan Hickox, Christopher M. Hirata, Matthew J. Holman, Jessica R. Lu, Ashley Villar

Abstract: The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recom… ▽ More The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recommend that a survey of the Galactic plane be defined early, based on the broad range of stakeholders for such a survey, the added scientific value of a first pass to obtain a baseline for proper motions complementary to Gaia's, and the significant potential synergies with ground-based surveys, notably the Legacy Survey of Space and Time (LSST) on Rubin. We also found strong motivation to follow a community definition process for ultra-deep observations with Roman. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Final report of the committee convened by NASA to assess community proposals for early definition science with Roman

Compton scattering of electrons in the intergalactic medium

Authors: Yuanyuan Yang, Heyang Long, Christopher M. Hirata

Abstract: This paper investigates the distribution and implications of cosmic ray electrons within the intergalactic medium (IGM). Utilizing a synthesis model of the extragalactic background, we evolve the spectrum of Compton-included cosmic rays. The energy density distribution of cosmic ray electrons peaks at redshift $z \approx2$, and peaks in the $\sim$MeV range. The fractional contribution of cosmic ra… ▽ More This paper investigates the distribution and implications of cosmic ray electrons within the intergalactic medium (IGM). Utilizing a synthesis model of the extragalactic background, we evolve the spectrum of Compton-included cosmic rays. The energy density distribution of cosmic ray electrons peaks at redshift $z \approx2$, and peaks in the $\sim$MeV range. The fractional contribution of cosmic ray pressure to the general IGM pressure progressively increases toward lower redshift. At mean density, the ratio of cosmic ray electron to thermal pressure in the IGM $ P_{\rm CRe} / P_{\rm th}$ is 0.3% at $z=2$, rising to 1.0% at $z=1$, and 1.8% at $z=0.1$ (considering only the cosmic rays produced locally by Compton scattering). We compute the linear Landau damping rate of plasma oscillations in the IGM caused by the $\sim$MeV cosmic ray electrons, and find it to be of order $\sim 10^{-6}\,\rm s^{-1}$ for wavenumbers $1.2\lesssim ck/ω_{\rm p}\lesssim 5$ at $z=2$ and mean density (where $ω_{\rm p}$ is the plasma frequency). This strongly affects the fate of TeV $e^+e^-$ pair beams produced by blazars, which are potentially unstable to oblique instabilities involving plasma oscillations with wavenumber $ck/ω_{\rm p}\approx\secθ$ ($θ$ being the angle between the beam and wave vector). Linear Landau damping is at least thousands of times faster than either pair beam instability growth or collisional effects; it thus turns off the pair beam instability except for modes with very small $θ$ ($ck/ω_{\rm p}\rightarrow 1$, where linear Landau damping is kinematically suppressed). This leaves open the question of whether the pair beam instability is turned off entirely, or can still proceed via the small-$θ$ modes. △ Less

Submitted 12 March, 2024; v1 submitted 30 November, 2023; originally announced November 2023.

Comments: 22 pages, 9 figures

Simulating image coaddition with the Nancy Grace Roman Space Telescope: II. Analysis of the simulated images and implications for weak lensing

Authors: Masaya Yamamoto, Katherine Laliotis, Emily Macbeth, Tianqing Zhang, Christopher M. Hirata, M. A. Troxel, Kaili Cao, Ami Choi, Jahmour Givans, Katrin Heitmann, Mustapha Ishak, Mike Jarvis, Eve Kovacs, Heyang Long, Rachel Mandelbaum, Andy Park, Anna Porredon, Christopher W. Walter, W. Michael Wood-Vasey

Abstract: One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the ou… ▽ More One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their 1- and 2-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1 - 6) x 10-4 per component depending on the band; the correlation functions are >= 2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman 2nd moment error budget, although the 4th moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all 4 filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman. △ Less

Submitted 12 January, 2024; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 25 pages, 20 figures. Submitted to MNRAS

Simulating image coaddition with the Nancy Grace Roman Space Telescope: I. Simulation methodology and general results

Authors: Christopher M. Hirata, Masaya Yamamoto, Katherine Laliotis, Emily Macbeth, M. A. Troxel, Tianqing Zhang, Kaili Cao, Ami Choi, Jahmour Givans, Katrin Heitmann, Mustapha Ishak, Mike Jarvis, Eve Kovacs, Heyang Long, Rachel Mandelbaum, Andy Park, Anna Porredon, Christopher W. Walter, W. Michael Wood-Vasey

Abstract: The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that… ▽ More The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired "target" point spread function (PSF). We present here an initial application of Imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that Imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the "leakage"), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ("Paper II") explores the implications for weak lensing analyses. △ Less

Submitted 12 January, 2024; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 28 pages, 19 figures, matches version accepted by MNRAS

Lyman-$α$ polarization from cosmological ionization fronts: II. Implications for intensity mapping

Authors: Emily Koivu, Heyang Long, Yuanyuan Yang, Christopher M. Hirata

Abstract: This is the second paper in a series whose aim is to predict the power spectrum of intensity and polarized intensity from cosmic reionization fronts. After building the analytic models for intensity and polarized intensity calculations in paper I, here we apply these models to simulations of reionization. We construct a geometric model for identifying front boundaries, calculate the intensity and… ▽ More This is the second paper in a series whose aim is to predict the power spectrum of intensity and polarized intensity from cosmic reionization fronts. After building the analytic models for intensity and polarized intensity calculations in paper I, here we apply these models to simulations of reionization. We construct a geometric model for identifying front boundaries, calculate the intensity and polarized intensity for each front, and compute a power spectrum of these results. This method was applied to different simulation sizes and resolutions, so we ensure that our results are convergent. We find that the power spectrum of fluctuations at $z=8$ in a bin of width $Δz=0.5$ ($λ/Δλ=18$) is $Δ_\ell \equiv [\ell(\ell+1)C_\ell/2π]^{1/2}$ is $3.2\times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ sr$^{-1}$ for the intensity $I$, $7.6\times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ sr$^{-1}$ for the $E$-mode polarization, and $5.8\times10^{-13}$ erg s$^{-1}$ cm$^{-2}$ sr$^{-1}$ for the $B$-mode polarization at $\ell=1.5\times10^4$. After computing the power spectrum, we compare results to detectable scales and discuss implications for observing this signal based on a proposed experiment. We find that, while fundamental physics does not exclude this kind of mapping from being attainable, an experiment would need to be highly ambitious and require significant advances to make mapping Lyman-$α$ polarization from cosmic reionization fronts a feasible goal. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 18 pages, 9 figures, to be submitted to JCAP

Lyman-α polarization from cosmological ionization fronts: I. Radiative transfer simulations

Authors: Yuanyuan Yang, Emily Koivu, Chenxiao Zeng, Heyang Long, Christopher M. Hirata

Abstract: In this paper, we present the formalism of simulating Lyman-$α$ emission and polarization around reionization ($z$ = 8) from a plane-parallel ionization front. We accomplish this by using a Monte Carlo method to simulate the production of a Lyman-$α$ photon, its propagation through an ionization front, and the eventual escape of this photon. This paper focuses on the relation of the input paramete… ▽ More In this paper, we present the formalism of simulating Lyman-$α$ emission and polarization around reionization ($z$ = 8) from a plane-parallel ionization front. We accomplish this by using a Monte Carlo method to simulate the production of a Lyman-$α$ photon, its propagation through an ionization front, and the eventual escape of this photon. This paper focuses on the relation of the input parameters of ionization front speed $U$, blackbody temperature $T_{\rm bb}$, and neutral hydrogen density $n_{\rm HI}$, on intensity $I$ and polarized intensity $P$ as seen by a distant observer. The resulting values of intensity range from $3.18\times 10^{-14}$ erg/cm$^{2}$/s/sr to $1.96 \times 10^{-9}$ erg/cm$^{2}$/s/sr , and the polarized intensity ranges from $5.73\times 10^{-17}$ erg/cm$^{2}$/s/sr to $5.31 \times 10^{-12}$ erg/cm$^{2}$/s/sr. We found that higher $T_{\rm bb}$, higher $U$, and higher $n_{\rm HI}$ contribute to higher intensity, as well as polarized intensity, though the strongest dependence was on the hydrogen density. The dependence of viewing angle of the front is also explored. We present tests to support the validity model, which makes the model suitable for further use in a following paper where we will calculate the intensity and polarized intensity power spectrum on a full reionization simulation. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 29 pages, 13 figures, to be submitted to JCAP

LISA Galactic Binaries in the Roman Galactic Bulge Time-Domain Survey

Authors: Matthew C. Digman, Chris M. Hirata

Abstract: Short-period Galactic white dwarf binaries detectable by LISA are the only guaranteed persistent sources for multi-messenger gravitational-wave astronomy. Large-scale surveys in the 2020s present an opportunity to conduct preparatory science campaigns to maximize the science yield from future multi-messenger targets. The Nancy Grace Roman Space Telescope Galactic Bulge Time Domain Survey will (in… ▽ More Short-period Galactic white dwarf binaries detectable by LISA are the only guaranteed persistent sources for multi-messenger gravitational-wave astronomy. Large-scale surveys in the 2020s present an opportunity to conduct preparatory science campaigns to maximize the science yield from future multi-messenger targets. The Nancy Grace Roman Space Telescope Galactic Bulge Time Domain Survey will (in its Reference Survey design) image seven fields in the Galactic Bulge approximately 40,000 times each. Although the Reference Survey cadence is optimized for detecting exoplanets via microlensing, it is also capable of detecting short-period white dwarf binaries. In this paper, we present forecasts for the number of detached short-period binaries the Roman Galactic Bulge Time Domain Survey will discover and the implications for the design of electromagnetic surveys. Although population models are highly uncertain, we find a high probability that the baseline survey will detect of order ~5 detached white dwarf binaries. The Reference Survey would also have a $\gtrsim20\%$ chance of detecting several known benchmark white dwarf binaries at the distance of the Galactic Bulge. △ Less

Submitted 30 December, 2022; originally announced December 2022.

Comments: 9 pages, 4 figure, 1 table

Self-calibrating optical galaxy cluster selection bias using cluster, galaxy, and shear cross-correlations

Authors: Chenxiao Zeng, Andrés N. Salcedo, Hao-Yi Wu, Christopher M. Hirata

Abstract: The clustering signals of galaxy clusters are known to be powerful tools for self-calibrating the mass-observable relation and are complementary to cluster abundance and lensing. In this work, we explore the possibility of combining three correlation functions -- cluster lensing, the cluster-galaxy cross-correlation function, and the galaxy auto-correlation function -- to self-calibrate optical cl… ▽ More The clustering signals of galaxy clusters are known to be powerful tools for self-calibrating the mass-observable relation and are complementary to cluster abundance and lensing. In this work, we explore the possibility of combining three correlation functions -- cluster lensing, the cluster-galaxy cross-correlation function, and the galaxy auto-correlation function -- to self-calibrate optical cluster selection bias, the boosted clustering and lensing signals in a richness-selected sample mainly caused by projection effects. We develop mock catalogues of redMaGiC-like galaxies and redMaPPer-like clusters by applying Halo Occupation Distribution (HOD) models to N-body simulations and using counts-in-cylinders around massive haloes as a richness proxy. In addition to the previously known small-scale boost in projected correlation functions, we find that the projection effects also significantly boost 3D correlation functions out to scales 100 $h^{-1} \mathrm{Mpc}$. We perform a likelihood analysis assuming survey conditions similar to that of the Dark Energy Survey (DES) and show that the selection bias can be self-consistently constrained at the 10% level. We discuss strategies for applying this approach to real data. We expect that expanding the analysis to smaller scales and using deeper lensing data would further improve the constraints on cluster selection bias. △ Less

Submitted 3 July, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 12 pages, 10 figures; replaced to match published version

Impact of inhomogeneous reionization on post-reionization 21 cm intensity mapping measurement of cosmological parameters

Authors: Heyang Long, Catalina Morales-Gutiérrez, Paulo Montero-Camacho, Christopher M. Hirata

Abstract: 21 cm intensity mapping (IM) has the potential to be a strong and unique probe of cosmology from redshift of order unity to redshift potentially as high as 30. For post-reionization 21 cm observations, the signal is modulated by the thermal and dynamical reaction of gas in the galaxies to the passage of ionization fronts during the Epoch of Reionization. In this work, we investigate the impact of… ▽ More 21 cm intensity mapping (IM) has the potential to be a strong and unique probe of cosmology from redshift of order unity to redshift potentially as high as 30. For post-reionization 21 cm observations, the signal is modulated by the thermal and dynamical reaction of gas in the galaxies to the passage of ionization fronts during the Epoch of Reionization. In this work, we investigate the impact of inhomogeneous reionization on the post-reionization 21 cm power spectrum and the induced shifts of cosmological parameters at redshifts $3.5 \lesssim z \lesssim 5.5$. We make use of hydrodynamics simulations that could resolve small-scale baryonic structure evolution to quantify HI abundance fluctuation, while semi-numerical large box 21cmFAST simulations capable of displaying inhomogeneous reionization process are deployed to track the inhomogeneous evolution of reionization bubbles. We discussed the prospects of capturing this effect in two post-reionization 21 cm intensity mapping experiments: SKA1-LOW and PUMA. We find the inhomogeneous reionization effect could impact the HI power spectrum up to tens of percent level and shift cosmological parameters estimation from sub-percent to tens percent in the observation of future post-reionization 21 cm intensity mapping experiments such as PUMA, while SKA1-LOW is likely to miss this effect at the redshifts of interest given the considered configuration. In particular, the shift is up to 0.0206 in the spectral index $n_s$ and 0.0192 eV in the sum of the neutrino masses $\sum m_ν$ depending on the reionization model and the observational parameters. We discuss strategies to mitigate and separate these biases. △ Less

Submitted 2 October, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 12 pages, 9 figures. Published in MNRAS

Probing Large Scale Ionizing Background Fluctuation with Lyman $α$ Forest and Galaxy Cross-correlation at z=2.4

Authors: Heyang Long, Christopher M. Hirata

Abstract: The amplitude of the metagalactic ultraviolet background (UVB) at large-scales is impacted by two factors. First, it naturally attenuates at scales larger than mean-free-path of UVB photons due to the absorption by neutral intergalactic medium. Second, there are discrete and rare ionizing sources distributing in the Universe, emitting the UVB photons, and thus enhancing the local UVB amplitude. Th… ▽ More The amplitude of the metagalactic ultraviolet background (UVB) at large-scales is impacted by two factors. First, it naturally attenuates at scales larger than mean-free-path of UVB photons due to the absorption by neutral intergalactic medium. Second, there are discrete and rare ionizing sources distributing in the Universe, emitting the UVB photons, and thus enhancing the local UVB amplitude. Therefore, for cosmological probe that is sensitive to the UVB amplitude and capable of detecting the large scale like Lyman-$α$ forest spectrum, the fluctuation due to the clustering of ionizing sources becomes a significant factor for Lyman-$α$ flux transmission and leave imprints on Lyman-$α$ flux power spectrum at these large scales. In this work, we make use of a radiative transfer model that parametrizes the UVB source distribution by its bias $b_{\rm j}$ and shot noise $\overline{n}_{\rm j}$. We estimate the constraints on this model through the cross-correlation between Lyman-$α$ forest survey and galaxy survey, using the DESI Lyman-$α$ forest survey and the Roman Space Telescope emission line galaxy survey as an example. We show the detection sensitivity improvement for UVB parameters from disjoint to maximal overlap of DESI+Roman survey strategy. We also show that the degeneracy of two ionizing source parameters can be broken by increasing the overlapping survey area. Our results motivate survey strategies more dedicated to probe the UVB large-scale fluctuations. △ Less

Submitted 18 January, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 11 pages, 5 figures. Added quasar lifetime modeling, consistent with published version on MNRAS

Pixel centroid characterization with laser speckle and application to the Nancy Grace Roman Space Telescope detector arrays

Authors: Christopher M. Hirata, Christopher Merchant

Abstract: The Nancy Grace Roman Space Telescope will use its wide-field instrument to carry out a suite of sky surveys in the near infrared. Several of the science objectives of these surveys, such as the measurement of the growth of cosmic structure using weak gravitational lensing, require exquisite control of instrument-related distortions of the images of astronomical objects. Roman will fly 4kx4k Teled… ▽ More The Nancy Grace Roman Space Telescope will use its wide-field instrument to carry out a suite of sky surveys in the near infrared. Several of the science objectives of these surveys, such as the measurement of the growth of cosmic structure using weak gravitational lensing, require exquisite control of instrument-related distortions of the images of astronomical objects. Roman will fly 4kx4k Teledyne H4RG-10 infrared detector arrays. This paper investigates whether the pixel centroids are located on a regular grid by projecting laser speckle patterns through a double slit aperture onto a non-flight detector array. We develop a method to reconstruct the pixel centroid offsets from the stochastic speckle pattern. Due to the orientation of the test setup, only x-offsets are measured here. We test the method both on simulations, and by injecting artificial offsets into the real images. We use cross-correlations of the reconstructions from different speckle realizations to determine how much of the variance in the pixel offset maps is signal (fixed to the detector) and how much is noise. After performing this reconstruction on 64x64 pixel patches, and fitting out the best-fit linear mapping from pixel index to position, we find that there are residual centroid offsets in the x (column) direction from a regular grid of 0.0107 pixels RMS (excluding shifts of an entire row relative to another, which our speckle patterns cannot constrain). This decreases to 0.0097 pix RMS if we consider residuals from a quadratic rather than linear mapping. These RMS offsets include both the physical pixel offsets, as well as any apparent offsets due to cross-talk and remaining systematic errors in the reconstruction. We comment on the advantages and disadvantages of speckle scene measurements as a tool for characterizing the pixel-level behavior in astronomical detectors. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 24 pages, 9 figures, to be submitted to PASP

Journal ref: Publications of the Astronomical Society of the Pacific, 134:115001 (2022)

Towards Powerful Probes of Neutrino Self-Interactions in Supernovae

Authors: Po-Wen Chang, Ivan Esteban, John F. Beacom, Todd A. Thompson, Christopher M. Hirata

Abstract: Neutrinos remain mysterious. As an example, enhanced self-interactions ($ν$SI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, $ν$SI should be important, but robust observables have been lacking. We show that $ν$SI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either… ▽ More Neutrinos remain mysterious. As an example, enhanced self-interactions ($ν$SI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, $ν$SI should be important, but robust observables have been lacking. We show that $ν$SI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either a burst or a steady-state wind; which occurs here is uncertain. Though the diffusive environment where neutrinos are produced may make a wind more likely, further work is needed to determine when each case is realized. In the burst-outflow case, $ν$SI increase the duration of the neutrino signal, and even a simple analysis of SN 1987A data has powerful sensitivity. For the wind-outflow case, we outline several promising ideas that may lead to new observables. Combined, these results are important steps towards solving the 35-year-old puzzle of how $ν$SI affect supernovae. △ Less

Submitted 28 June, 2023; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 5+19 pages, 4+14 figures. Accepted for publication in Phys.Rev.Lett. Match the version to be published: minor revisions for clarity; moved discussion about supernova neutrino properties to Supplemental Material C; updated references

Journal ref: Phys. Rev. Lett. 131, 071002 (2023)

Quantum yield and charge diffusion in the Nancy Grace Roman Space Telescope infrared detectors

Authors: Jahmour J. Givans, Ami Choi, Anna Porredon, Jenna K. C. Freudenburg, Christopher M. Hirata, Robert J. Hill, Christopher Bennett, Roger Foltz, Lane Meier

Abstract: The shear signal required for weak lensing analyses is small, so any detector-level effects which distort astronomical images can contaminate the inferred shear. The Nancy Grace Roman Space Telescope (Roman) will fly a focal plane with 18 Teledyne H4RG-10 near infrared (IR) detector arrays; these have never been used for weak lensing and they present unique instrument calibration challenges. A pai… ▽ More The shear signal required for weak lensing analyses is small, so any detector-level effects which distort astronomical images can contaminate the inferred shear. The Nancy Grace Roman Space Telescope (Roman) will fly a focal plane with 18 Teledyne H4RG-10 near infrared (IR) detector arrays; these have never been used for weak lensing and they present unique instrument calibration challenges. A pair of previous investigations (Hirata & Choi 2020; Choi & Hirata 2020) demonstrated that spatiotemporal correlations of flat fields can effectively separate the brighter-fatter effect (BFE) and interpixel capacitance (IPC). Later work (Freudenburg et al. 2020) introduced a Fourier-space treatment of these correlations which allowed the authors to expand to higher orders in BFE, IPC, and classical nonlinearity (CNL). This work expands the previous formalism to include quantum yield and charge diffusion. We test the updated formalism on simulations and show that we can recover input characterization values. We then apply the formalism to visible and IR flat field data from three Roman flight candidate detectors. We fit a 2D Gaussian model to the charge diffusion at 0.5 $μ$m wavelength, and find variances of $C_{11} = 0.1066\pm 0.0011$ pix$^2$ in the horizontal direction, $C_{22} = 0.1136\pm 0.0012$ pix$^2$ in the vertical direction, and a covariance of $C_{12} = 0.0001\pm 0.0007$ pix$^2$ (stat) for SCA 20829. Last, we convert the asymmetry of the charge diffusion into an equivalent shear signal, and find a contamination of the shear correlation function to be $ξ_+ \sim 10^{-6}$ for each detector. This exceeds Roman's allotted error budget for the measurement by a factor of $\mathcal{O}(10)$ in power (amplitude squared) but can likely be mitigated through standard methods for fitting the point spread function (PSF) since charge diffusion can be treated as a contribution to the PSF. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 26 pages, 9 figures

The High Latitude Spectroscopic Survey on the Nancy Grace Roman Space Telescope

Authors: Yun Wang, Zhongxu Zhai, Anahita Alavi, Elena Massara, Alice Pisani, Andrew Benson, Christopher M. Hirata, Lado Samushia, David H. Weinberg, James Colbert, Olivier Doré, Tim Eifler, Chen Heinrich, Shirley Ho, Elisabeth Krause, Nikhil Padmanabhan, David Spergel, Harry I. Teplitz

Abstract: The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $μ$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$σ$, requiring $\sim$0.6 yrs of observing time.… ▽ More The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $μ$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$σ$, requiring $\sim$0.6 yrs of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semi-analytic galaxy formation model to a cosmological N-body simulation, and use this mock survey to create pixel-level simulations of 4 deg$^2$ of HLSS grism spectroscopy. We find that the reference HLSS would measure $\sim$ 10 million H$α$ galaxy redshifts that densely map large scale structure at $z=1-2$ and 2 million [OIII] galaxy redshifts that sparsely map structures at $z=2-3$. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the growth of large scale structure with redshift space distortions. We also study possible deviations from the reference design, and find that a deep HLSS at $f_{\rm line}>7\times10^{-17}$erg/s/cm$^2$ over 4000 deg$^2$ (requiring $\sim$1.5 yrs of observing time) provides the most compelling stand-alone constraints on dark energy from Roman alone. This provides a useful reference for future optimizations. The reference survey, simulated data sets, and forecasts presented here will inform community decisions on the final scope and design of the Roman HLSS. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 29 pages, 8 figures, ApJ submitted

Streaming Velocity Effects on the Post-reionization 21 cm Baryon Acoustic Oscillation Signal

Authors: Heyang Long, Jahmour J. Givans, Christopher M. Hirata

Abstract: The relative velocity between baryons and dark matter in the early Universe can suppress the formation of small-scale baryonic structure and leave an imprint on the baryon acoustic oscillation (BAO) scale at low redshifts after reionization. This "streaming velocity" affects the post-reionization gas distribution by directly reducing the abundance of pre-existing mini-halos (… ▽ More The relative velocity between baryons and dark matter in the early Universe can suppress the formation of small-scale baryonic structure and leave an imprint on the baryon acoustic oscillation (BAO) scale at low redshifts after reionization. This "streaming velocity" affects the post-reionization gas distribution by directly reducing the abundance of pre-existing mini-halos ($\lesssim 10^7 M_{\bigodot}$) that could be destroyed by reionization and indirectly modulating reionization history via photoionization within these mini-halos. In this work, we investigate the effect of streaming velocity on the BAO feature in HI 21 cm intensity mapping after reionization, with a focus on redshifts $3.5\lesssim z\lesssim5.5$. We build a spatially modulated halo model that includes the dependence of the filtering mass on the local reionization redshift and thermal history of the intergalactic gas. In our fiducial model, we find isotropic streaming velocity bias coefficients $b_v$ ranging from $-0.0043$ at $z=3.5$ to $-0.0273$ at $z=5.5$, which indicates that the BAO scale is stretched (i.e., the peaks shift to lower $k$). In particular, streaming velocity shifts the transverse BAO scale between 0.121% ($z=3.5$) and 0.35% ($z=5.5$) and shifts the radial BAO scale between 0.167% ($z=3.5$) and 0.505% ($z=5.5$). These shifts exceed the projected error bars from the more ambitious proposed hemispherical-scale surveys in HI (0.13% at $1σ$ per $Δz = 0.5$ bin). △ Less

Submitted 7 March, 2022; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 12 pages, 5 figures, accepted for publication in MNRAS

Impact of Image Persistence in the Roman Space Telescope High-Latitude Survey

Authors: Chien-Hao Lin, Rachel Mandelbaum, M. A. Troxel, Christopher M. Hirata, Mike Jarvis

Abstract: The High Latitude Survey of the Nancy Grace Roman Space Telescope is expected to measure the positions and shapes of hundreds of millions of galaxies in an area of 2220 deg$^2$. This survey will provide high-quality weak lensing data with unprecedented systematics control. The Roman Space Telescope will survey the sky in near infrared (NIR) bands using Teledyne H4RG HgCdTe photodiode arrays. These… ▽ More The High Latitude Survey of the Nancy Grace Roman Space Telescope is expected to measure the positions and shapes of hundreds of millions of galaxies in an area of 2220 deg$^2$. This survey will provide high-quality weak lensing data with unprecedented systematics control. The Roman Space Telescope will survey the sky in near infrared (NIR) bands using Teledyne H4RG HgCdTe photodiode arrays. These NIR arrays exhibit an effect called persistence: charges that are trapped in the photodiodes during earlier exposures are gradually released into later exposures, leading to contamination of the images and potentially to errors in measured galaxy properties such as fluxes and shapes. In this work, we use image simulations that incorporate the persistence effect to study its impact on galaxy shape measurements and weak lensing signals. No significant spatial correlations are found between the galaxy shape changes induced by persistence. On the scales of interest for weak lensing cosmology, the effect of persistence on the weak lensing correlation function is about two orders of magnitude lower than the Roman Space Telescope additive shear error budget, indicating that the persistence effect is expected to be a subdominant contributor to the systematic error budget for weak lensing with the Roman Space Telescope given its current design. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 7 pages, 7 figures

Line confusion in spectroscopic surveys and its possible effects: Shifts in Baryon Acoustic Oscillations position

Authors: Elena Massara, Shirley Ho, Christopher M. Hirata, Joseph DeRose, Risa H. Wechsler, Xiao Fang

Abstract: Roman Space Telescope will survey about 17 million emission-line galaxies over a range of redshifts. Its main targets are H$α$ emission-line galaxies at low redshifts and [O III] emission-line galaxies at high redshifts. The Roman Space Telescope will estimate the redshift these galaxies with single line identification. This suggests that other emission-line galaxies may be misidentified as the ma… ▽ More Roman Space Telescope will survey about 17 million emission-line galaxies over a range of redshifts. Its main targets are H$α$ emission-line galaxies at low redshifts and [O III] emission-line galaxies at high redshifts. The Roman Space Telescope will estimate the redshift these galaxies with single line identification. This suggests that other emission-line galaxies may be misidentified as the main targets. In particular, it is hard to distinguish between the H$β$ and [O III] lines as the two lines are close in wavelength and hence the photometric information may not be sufficient to separate them reliably. Misidentifying H$β$ emitter as [O III] emitter will cause a shift in the inferred radial position of the galaxy by approximately 90 Mpc/h. This length scale is similar to the Baryon Acoustic Oscillation (BAO) scale and could shift and broaden the BAO peak, possibly introduce errors in determining the BAO peak position. We qualitatively describe the effect of this new systematic and further quantify it with a lightcone simulation with emission-line galaxies. △ Less

Submitted 30 September, 2020; originally announced October 2020.

Non-equilibrium temperature evolution of ionization fronts during the Epoch of Reionization

Authors: Chenxiao Zeng, Christopher M. Hirata

Abstract: The epoch of reionization (EoR) marks the end of the Cosmic Dawn and the beginning of large-scale structure formation in the universe. The impulsive ionization fronts (I-fronts) heat and ionize the gas within the reionization bubbles in the intergalactic medium (IGM). The temperature during this process is a key yet uncertain ingredient in current models. Typically, reionization simulations assume… ▽ More The epoch of reionization (EoR) marks the end of the Cosmic Dawn and the beginning of large-scale structure formation in the universe. The impulsive ionization fronts (I-fronts) heat and ionize the gas within the reionization bubbles in the intergalactic medium (IGM). The temperature during this process is a key yet uncertain ingredient in current models. Typically, reionization simulations assume that all baryonic species are in instantaneous thermal equilibrium with each other during the passage of an I-front. Here we present a new model of the temperature evolution for the ionization front by studying non-equilibrium effects. In particular, we include the energy transfer between major baryon species ($e^{-}$, \HI, \HII, \HeI, and \HeII) and investigate their impacts on the post-ionization front temperature $T_{\mathrm{re}}$. For a better step-size control when solving the stiff equations, we implement an implicit method and construct an energy transfer rate matrix. We find that the assumption of equilibration is valid for a low-speed ionization front ($\lessapprox\ 10^9~\mathrm{cm}/\mathrm{s}$), but deviations from equilibrium occur for faster fronts. The post-front temperature $T_{\mathrm{re}}$ is lower by up to 19.7\% (at $3\times 10^9$ cm/s) or 30.8\% (at $10^{10}$ cm/s) relative to the equilibrium case. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 9 pages, 6 figures

Brighter-fatter effect in near-infrared detectors -- III. Fourier-domain treatment of flat field correlations and application to WFIRST

Authors: Jenna K. C. Freudenburg, Jahmour J. Givans, Ami Choi, Christopher M. Hirata, Chris Bennett, Stephanie Cheung, Analia Cillis, Dave Cottingham, Robert J. Hill, Jon Mah, Lane Meier

Abstract: Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. WFIRST will fly a focal plane containing 18 Teledyne H4RG-10 near infrared detector arrays, which present different instrument cali… ▽ More Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. WFIRST will fly a focal plane containing 18 Teledyne H4RG-10 near infrared detector arrays, which present different instrument calibration challenges from previous weak lensing observations. Previous work has shown that correlation functions of flat field images are effective tools for disentangling linear and non-linear inter-pixel capacitance (IPC) and the brighter-fatter effect (BFE). Here we present a Fourier-domain treatment of the flat field correlations, which allows us to expand the previous formalism to all orders in IPC, BFE, and classical non-linearity. We show that biases in simulated flat field analyses in Paper I are greatly reduced through the use of this formalism. We then apply this updated formalism to flat field data from three WFIRST flight candidate detectors, and explore the robustness to variations in the analysis. We find that the BFE is present in all three detectors, and that its contribution to the flat field correlations dominates over the non-linear IPC. The magnitude of the BFE is such that the effective area of a pixel is increased by $(3.54\pm0.03)\times 10^{-7}$ for every electron deposited in a neighboring pixel. We compare IPC maps from flat field autocorrelation measurements to those obtained from the single pixel reset method and find a median difference of 0.113%. After further diagnosis of this difference, we ascribe it largely to an additional source of cross-talk, the vertical trailing pixel effect, and recommend further work to develop a model for this effect. These results represent a significant step toward calibration of the non-ideal effects in WFIRST detectors. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Redshift-space streaming velocity effects on the Lyman-$α$ forest baryon acoustic oscillation scale

Authors: Jahmour J. Givans, Christopher M. Hirata

Abstract: The baryon acoustic oscillation (BAO) scale acts as a standard ruler for measuring cosmological distances and has therefore emerged as a leading probe of cosmic expansion history. However, any physical effect that alters the length of the ruler can lead to a bias in our determination of distance and expansion rate. One of these physical effects is the streaming velocity, the relative velocity betw… ▽ More The baryon acoustic oscillation (BAO) scale acts as a standard ruler for measuring cosmological distances and has therefore emerged as a leading probe of cosmic expansion history. However, any physical effect that alters the length of the ruler can lead to a bias in our determination of distance and expansion rate. One of these physical effects is the streaming velocity, the relative velocity between baryons and dark matter in the early Universe, which couples to the BAO scale due to their common origin in acoustic waves at recombination. In this work, we investigate the impact of streaming velocity on the BAO feature of the Lyman-$α$ forest auto-power spectrum, one of the main tracers being used by the recently commissioned DESI spectrograph. To do this, we develop a new perturbative model for Lyman-$α$ flux density contrast which is complete to second order for a certain set of fields, and applicable to any redshift-space tracer of structure since it is based only on symmetry considerations. We find that there are 8 biasing coefficients through second order. We find streaming velocity-induced shifts in the BAO scale of 0.081--0.149% (transverse direction) and 0.053--0.058% (radial direction), depending on the model for the biasing coefficients used. These are smaller than, but not negligible compared to, the DESI Lyman-$α$ BAO error budget, which is 0.46% on the overall scale. The sensitivity of these results to our choice of bias parameters underscores the need for future work to measure the higher-order biasing coefficients from simulations, especially for future experiments beyond DESI. △ Less

Submitted 9 August, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: 25 pages, 5 figures. This version is the closest match to the PRD final print

Journal ref: Phys. Rev. D 102, 023515 (2020)

A Synthetic Roman Space Telescope High-Latitude Imaging Survey: Simulation Suite and the Impact of Wavefront Errors on Weak Gravitational Lensing

Authors: M. A. Troxel, H. Long, C. M. Hirata, A. Choi, M. Jarvis, R. Mandelbaum, K. Wang, M. Yamamoto, S. Hemmati, P. Capak

Abstract: The Nancy Grace Roman Space Telescope (Roman) mission is expected to launch in the mid-2020s. Its weak lensing program is designed to enable unprecedented systematics control in photometric measurements, including shear recovery, point-spread function (PSF) correction, and photometric calibration. This will enable exquisite weak lensing science and allow us to adjust to and reliably contribute to… ▽ More The Nancy Grace Roman Space Telescope (Roman) mission is expected to launch in the mid-2020s. Its weak lensing program is designed to enable unprecedented systematics control in photometric measurements, including shear recovery, point-spread function (PSF) correction, and photometric calibration. This will enable exquisite weak lensing science and allow us to adjust to and reliably contribute to the cosmological landscape after the initial years of observations from other concurrent Stage IV dark energy experiments. This potential requires equally careful planning and requirements validation as the mission prepares to enter its construction phase. We present a suite of image simulations based on GalSim that are used to construct a complex, synthetic Roman weak lensing survey that incorporates realistic input galaxies and stars, relevant detector non-idealities, and the current reference five-year Roman survey strategy. We present a first study to empirically validate the existing Roman weak lensing requirements flowdown using a suite of 12 matched image simulations, each representing a different perturbation to the wavefront or image motion model. These are chosen to induce a range of potential static and low- and high-frequency time-dependent PSF model errors. We analyze the measured shapes of galaxies from each of these simulations and compare them to a reference, fiducial simulation to infer the response of the shape measurement to each of these modes in the wavefront model. We then compare this to existing analytic flowdown requirements, and find general agreement between the empirically derived response and that predicted by the analytic model. △ Less

Submitted 4 December, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: Matches accepted version; substantial additions to appendices

The Impact of Light Polarization Effects on Weak Lensing Systematics

Authors: Chien-Hao Lin, Brent Tan, Rachel Mandelbaum, Christopher M. Hirata

Abstract: A fraction of the light observed from edge-on disk galaxies is polarized due to two physical effects: selective extinction by dust grains aligned with the magnetic field, and scattering of the anisotropic starlight field. Since the reflection and transmission coefficients of the reflecting and refracting surfaces in an optical system depend on the polarization of incoming rays, this optical polari… ▽ More A fraction of the light observed from edge-on disk galaxies is polarized due to two physical effects: selective extinction by dust grains aligned with the magnetic field, and scattering of the anisotropic starlight field. Since the reflection and transmission coefficients of the reflecting and refracting surfaces in an optical system depend on the polarization of incoming rays, this optical polarization produces both (a) a selection bias in favor of galaxies with specific orientations and (b) a polarization-dependent PSF. In this work we build toy models to obtain for the first time an estimate for the impact of polarization on PSF shapes and the impact of the selection bias due to the polarization effect on the measurement of the ellipticity used in shear measurements. In particular, we are interested in determining if this effect will be significant for WFIRST. We show that the systematic uncertainties in the ellipticity components are $8\times 10^{-5}$ and $1.1 \times 10^{-4}$ due to the selection bias and PSF errors respectively. Compared to the overall requirements on knowledge of the WFIRST PSF ellipticity ($4.7\times 10^{-4}$ per component), both of these systematic uncertainties are sufficiently close to the WFIRST tolerance level that more detailed studies of the polarization effects or more stringent requirements on polarization-sensitive instrumentation for WFIRST are required. △ Less

Submitted 19 January, 2021; v1 submitted 11 October, 2019; originally announced October 2019.

Comments: 8 pages, 3 figures, updated to published version

A Framework for Measuring Weak-Lensing Magnification Using the Fundamental Plane

Authors: Jenna K. C. Freudenburg, Eric M. Huff, Christopher M. Hirata

Abstract: Galaxy-galaxy lensing is an essential tool for probing dark matter halos and constraining cosmological parameters. While galaxy-galaxy lensing measurements usually rely on shear, weak-lensing magnification contains additional constraining information. Using the fundamental plane (FP) of elliptical galaxies to anchor the size distribution of a background population is one method that has been propo… ▽ More Galaxy-galaxy lensing is an essential tool for probing dark matter halos and constraining cosmological parameters. While galaxy-galaxy lensing measurements usually rely on shear, weak-lensing magnification contains additional constraining information. Using the fundamental plane (FP) of elliptical galaxies to anchor the size distribution of a background population is one method that has been proposed for performing a magnification measurement. We present a formalism for using the FP residuals of elliptical galaxies to jointly estimate the foreground mass and background redshift errors for a stacked lens scenario. The FP residuals include information about weak-lensing magnification $κ$, and therefore foreground mass, since to first order, nonzero $κ$ affects galaxy size but not other FP properties. We also present a modular, extensible code that implements the formalism using emulated galaxy catalogs of a photometric galaxy survey. We find that combining FP information with observed number counts of the source galaxies constrains mass and photo-z error parameters significantly better than an estimator that includes number counts only. In particular, the constraint on the mass is 17.0\% if FP residuals are included, as opposed to 27.7\% when only number counts are included. The effective size noise for a foreground lens of mass $M_H=10^{14}M_\odot$, with a conservative selection function in size and surface brightness applied to the source population, is $σ_{κ,\mathrm{eff}}=0.250$. We discuss the improvements to our FP model necessary to make this formalism a practical companion to shear analyses in weak lensing surveys. △ Less

Submitted 3 November, 2019; v1 submitted 7 October, 2019; originally announced October 2019.

Detecting Magnetic Fields in Exoplanets with Spectropolarimetry of the Helium Line at 1083 nm

Authors: Antonija Oklopčić, Makana Silva, Paulo Montero-Camacho, Christopher M. Hirata

Abstract: The magnetic fields of the solar system planets provide valuable insights into the planets' interiors and can have dramatic consequences for the evolution of their atmospheres and interaction with the solar wind. However, we have little direct knowledge of magnetic fields in exoplanets. Here we present a method for detecting magnetic fields in the atmospheres of close-in exoplanets based on spectr… ▽ More The magnetic fields of the solar system planets provide valuable insights into the planets' interiors and can have dramatic consequences for the evolution of their atmospheres and interaction with the solar wind. However, we have little direct knowledge of magnetic fields in exoplanets. Here we present a method for detecting magnetic fields in the atmospheres of close-in exoplanets based on spectropolarimetric transit observations at the wavelength of the helium line at 1083 nm. This methodology has been successfully applied before for exploring magnetic fields in solar coronal filaments. Strong absorption signatures (transit depths on the order of a few percent) in the 1083 nm line have recently been observed for several close-in exoplanets. We show that in the conditions in these escaping atmospheres, metastable helium atoms should be optically pumped by the starlight and, for field strengths more than a few $\times 10^{-4}$ G, should align with the magnetic field. This results in linearly polarized absorption at 1083 nm that traces the field direction (the Hanle effect), which we explore by both analytic computation and with the Hazel numerical code. The linear polarization $\sqrt{Q^2+U^2}/I$ ranges from $\sim 10^{-3}$ in optimistic cases down to a few $\times 10^{-5}$ for particularly unfavorable cases, with very weak dependence on field strength. The line-of-sight component of the field results in a slight circular polarization (the Zeeman effect), also reaching $V/I\sim {\rm few}\times 10^{-5}(B_\parallel/10\,{\rm G})$. We discuss the detectability of these signals with current (SPIRou) and future (extremely large telescope) high-resolution infrared spectropolarimeters, and we briefly comment on possible sources of astrophysical contamination. △ Less

Submitted 10 April, 2020; v1 submitted 6 October, 2019; originally announced October 2019.

Comments: published in The Astrophysical Journal on February 14, 2020; updated to match the published version

Journal ref: ApJ 890, 88 (2020)

ATLAS Probe: Breakthrough Science of Galaxy Evolution, Cosmology, Milky Way, and the Solar System

Authors: Yun Wang, Mark Dickinson, Lynne Hillenbrand, Massimo Robberto, Lee Armus, Mario Ballardini, Robert Barkhouser, James Bartlett, Peter Behroozi, Robert A. Benjamin, Jarle Brinchmann, Ranga-Ram Chary, Chia-Hsun Chuang, Andrea Cimatti, Charlie Conroy, Robert Content, Emanuele Daddi, Megan Donahue, Olivier Dore, Peter Eisenhardt, Henry C. Ferguson, Andreas Faisst, Wesley C. Fraser, Karl Glazebrook, Varoujan Gorjian , et al. (23 additional authors not shown)

Abstract: ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and… ▽ More ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope. △ Less

Submitted 30 August, 2019; originally announced September 2019.

Comments: APC white paper submitted to Astro2020. arXiv admin note: substantial text overlap with arXiv:1802.01539

(Not as) Big as a Barn: Upper Bounds on Dark Matter-Nucleus Cross Sections

Authors: Matthew C. Digman, Christopher V. Cappiello, John F. Beacom, Christopher M. Hirata, Annika H. G. Peter

Abstract: Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is… ▽ More Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $σ_{χA} \propto A^2 μ_A^2 σ_{χN}\propto A^4 σ_{χN}$, where the $A^2$ comes from coherence and the $μ_A^2\simeq A^2 m_N^2$ from kinematics for $m_χ\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $σ_{χA} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $σ_{χN}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $σ_{χN}>10^{-32}\;\text{cm}^2$ with $m_χ\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted. △ Less

Submitted 30 December, 2022; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: 18 pages, 7 figures, revised to match journal version including erratum, comments are welcome

Journal ref: Phys. Rev. D 100, 063013 (2019)

Revisiting constraints on asteroid-mass primordial black holes as dark matter candidates

Authors: Paulo Montero-Camacho, Xiao Fang, Gabriel Vasquez, Makana Silva, Christopher M. Hirata

Abstract: As the only dark matter candidate that does not invoke a new particle that survives to the present day, primordial black holes (PBHs) have drawn increasing attention recently. Up to now, various observations have strongly constrained most of the mass range for PBHs, leaving only small windows where PBHs could make up a substantial fraction of the dark matter. Here we revisit the PBH constraints fo… ▽ More As the only dark matter candidate that does not invoke a new particle that survives to the present day, primordial black holes (PBHs) have drawn increasing attention recently. Up to now, various observations have strongly constrained most of the mass range for PBHs, leaving only small windows where PBHs could make up a substantial fraction of the dark matter. Here we revisit the PBH constraints for the asteroid-mass window, i.e., the mass range $3.5\times 10^{-17}M_\odot < m_{\mathrm{PBH}} < 4\times 10^{-12}M_\odot$. We revisit 3 categories of constraints. (1) For optical microlensing, we analyze the finite source size and diffractive effects and discuss the scaling relations between the event rate, $m_{\mathrm{PBH}}$ and the event duration. We argue that it will be difficult to push the existing optical microlensing constraints to much lower m$_{\mathrm{PBH}}$. (2) For dynamical capture of PBHs in stars, we derive a general result on the capture rate based on phase space arguments. We argue that survival of stars does not constrain PBHs, but that disruption of stars by captured PBHs should occur and that the asteroid-mass PBH hypothesis could be constrained if we can work out the observational signature of this process. (3) For destruction of white dwarfs by PBHs that pass through the white dwarf without getting gravitationally captured, but which produce a shock that ignites carbon fusion, we perform a 1+1D hydrodynamic simulation to explore the post-shock temperature and relevant timescales, and again we find this constraint to be ineffective. In summary, we find that the asteroid-mass window remains open for PBHs to account for all the dark matter. △ Less

Submitted 26 August, 2019; v1 submitted 13 June, 2019; originally announced June 2019.

Comments: Updated to match accepted version for publication. 43 pages

Journal ref: JCAP08(2019)031

Brighter-fatter effect in near-infrared detectors -- II. Auto-correlation analysis of H4RG-10 flats

Authors: Ami Choi, Christopher M. Hirata

Abstract: The Wide Field Infrared Survey Telescope (WFIRST) will investigate the origins of cosmic acceleration using weak gravitational lensing at near infrared wavelengths. Lensing analyses place strict constraints on the precision of size and ellipticity measurements of the point spread function. WFIRST will use infrared detector arrays, which must be fully characterized to inform data reduction and cali… ▽ More The Wide Field Infrared Survey Telescope (WFIRST) will investigate the origins of cosmic acceleration using weak gravitational lensing at near infrared wavelengths. Lensing analyses place strict constraints on the precision of size and ellipticity measurements of the point spread function. WFIRST will use infrared detector arrays, which must be fully characterized to inform data reduction and calibration procedures such that unbiased cosmological results can be achieved. Hirata & Choi 2019 introduces formalism to connect the cross-correlation signal of different flat field time samples to non-linear detector behaviors such as the brighter fatter effect (BFE) and non-linear inter-pixel capacitance (NL-IPC), and this paper applies that framework to a WFIRST development detector, SCA 18237. We find a residual correlation signal after accounting for classical non-linearity. This residual correlation contains a combination of the BFE and NL-IPC; however, further tests suggest that the BFE is the dominant mechanism. If interpreted as a pure BFE, it suggests that the effective area of a pixel is increased by $(2.87\pm0.03)\times 10^{-7}$ (stat.) for every electron in the 4 nearest neighbors, with a rapid $\sim r^{-5.6\pm0.2}$ fall-off of the effect for more distant neighbors. We show that the IPC inferred from hot pixels contains the same large-scale spatial variations as the IPC inferred from auto-correlations, albeit with an overall offset of $\sim 0.06\%$. The NL-IPC inferred from hot pixels is too small to explain the cross-correlation measurement, further supporting the BFE hypothesis. This work presents the first evidence for the BFE in an H4RG-10 detector, demonstrates some of the useful insights that can be gleaned from flat field statistics, and represents a significant step towards calibration of WFIRST data. △ Less

Submitted 23 January, 2020; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: 28 pages, 8 figures, 3 tables, version accepted by PASP

Journal ref: Publications of the Astronomical Society of the Pacific, Volume 132, Issue 1007, pp. 014502 (2020)

Brighter-fatter effect in near-infrared detectors -- I. Theory of flat auto-correlations

Authors: Christopher M. Hirata, Ami Choi

Abstract: Weak gravitational lensing studies aim to measure small distortions in the shapes of distant galaxies, and thus place very tight demands on the understanding of detector-induced systematic effects in astronomical images. The Wide-Field Infrared Survey Telescope (WFIRST) will carry out weak lensing measurements in the near infrared using the new Teledyne H4RG-10 detector arrays, which makes the ran… ▽ More Weak gravitational lensing studies aim to measure small distortions in the shapes of distant galaxies, and thus place very tight demands on the understanding of detector-induced systematic effects in astronomical images. The Wide-Field Infrared Survey Telescope (WFIRST) will carry out weak lensing measurements in the near infrared using the new Teledyne H4RG-10 detector arrays, which makes the range of possible detector systematics very different from traditional weak lensing measurements using optical CCDs. One of the non-linear detector effects observed in CCDs is the brighter-fatter effect (BFE), in which charge already accumulated in a pixel alters the electric field geometry and causes new charge to be deflected away from brighter pixels. Here we describe the formalism for measuring the BFE using flat field correlation functions in infrared detector arrays. The auto-correlation of CCD flat fields is often used to measure the BFE, but because the infrared detector arrays are read out with the charge "in place," the flat field correlations are dominated by capacitive cross-talk between neighboring pixels (the inter-pixel capacitance, or IPC). Conversely, if the BFE is present and one does not account for it, it can bias correlation measurements of the IPC and photon transfer measurements of the gain. We show that because the infrared detector arrays can be read out non-destructively, one can compute numerous cross-correlation functions between different time slices of the same flat exposures, and that correlations due to IPC and BFE leave distinct imprints. We generate a suite of simulated flat fields and show that the underlying IPC and BFE parameters can be extracted, even when both are present in the simulation. There are some biases in the BFE coefficients up to 12%, which are likely caused by higher order terms that are dropped from this analysis. △ Less

Submitted 23 January, 2020; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: 45 pages, 4 figures, 1 table, version accepted by PASP, abstract abridged

Journal ref: Publications of the Astronomical Society of the Pacific, Volume 132, Issue 1007, pp. 014501 (2020)

Forecasting Super-Sample Covariance in Future Weak Lensing Surveys with SuperSCRAM

Authors: Matthew C. Digman, Joseph E. McEwen, Christopher M. Hirata

Abstract: The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constra… ▽ More The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to super-sample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0. △ Less

Submitted 26 April, 2019; originally announced April 2019.

Comments: 32 pages, 7 figures, 1 table, comments/suggestions welcome!

An Ultra Deep Field survey with WFIRST

Authors: Anton M. Koekemoer, R. J. Foley, D. N. Spergel, M. Bagley, R. Bezanson, F. B. Bianco, R. Bouwens, L. Bradley, G. Brammer, P. Capak, I. Davidzon, G. De Rosa, M. E. Dickinson, O. Doré, J. S. Dunlop, R. S. Ellis, X. Fan, G. G. Fazio, H. C. Ferguson, A. V. Filippenko, S. Finkelstein, B. Frye, E. Gawiser, N. A. Grogin, N. P. Hathi , et al. (47 additional authors not shown)

Abstract: Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their siz… ▽ More Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.) △ Less

Submitted 19 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Cosmic voids: a novel probe to shed light on our Universe

Authors: Alice Pisani, Elena Massara, David N. Spergel, David Alonso, Tessa Baker, Yan-Chuan Cai, Marius Cautun, Christopher Davies, Vasiliy Demchenko, Olivier Doré, Andy Goulding, Mélanie Habouzit, Nico Hamaus, Adam Hawken, Christopher M. Hirata, Shirley Ho, Bhuvnesh Jain, Christina D. Kreisch, Federico Marulli, Nelson Padilla, Giorgia Pollina, Martin Sahlén, Ravi K. Sheth, Rachel Somerville, Istvan Szapudi , et al. (4 additional authors not shown)

Abstract: Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust… ▽ More Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper for the Astro2020 Decadal Survey, 5 pages text, 14 pages with references and affiliations

Messengers from the Early Universe: Cosmic Neutrinos and Other Light Relics

Authors: Daniel Green, Mustafa A. Amin, Joel Meyers, Benjamin Wallisch, Kevork N. Abazajian, Muntazir Abidi, Peter Adshead, Zeeshan Ahmed, Behzad Ansarinejad, Robert Armstrong, Carlo Baccigalupi, Kevin Bandura, Darcy Barron, Nicholas Battaglia, Daniel Baumann, Keith Bechtol, Charles Bennett, Bradford Benson, Florian Beutler, Colin Bischoff, Lindsey Bleem, J. Richard Bond, Julian Borrill, Elizabeth Buckley-Geer, Cliff Burgess , et al. (114 additional authors not shown)

Abstract: The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic… ▽ More The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 5 pages + references; 1 figure; science white paper submitted to the Astro2020 decadal survey

WFIRST: Enhancing Transient Science and Multi-Messenger Astronomy

Authors: Ryan J. Foley, Joshua S. Bloom, S. Bradley Cenko, Ryan Chornock, Georgios Dimitriadis, Olivier Dore, Alexei V. Filippenko, Ori D. Fox, Christopher M. Hirata, Saurabh W. Jha, David O. Jones, Mansi Kasliwal, Patrick L. Kelly, Charles D. Kilpatrick, Robert P. Kirshner, Anton M. Koekemoer, Jeffrey W. Kruk, Kaisey S. Mandel, Raffaella Margutti, Vivian Miranda, Samaya Nissanke, Armin Rest, Jason Rhodes, Steven A. Rodney, Benjamin M. Rose , et al. (9 additional authors not shown)

Abstract: Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first obj… ▽ More Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020

Astro2020: Empirically Constraining Galaxy Evolution

Authors: Peter Behroozi, Matthew R. Becker, Frank C. van den Bosch, Jarle Brinchmann, Charlie Conroy, Mark Dickinson, Christopher M. Hirata, Andrew Hearin, Alexie Leauthaud, Chun Ly, Yao-Yuan Mao, Benjamin P. Moster, Christine O'Donnell, Casey Papovich, Aldo Rodríguez-Puebla, Rachel Somerville, Erik Tollerud, Jeremy Tinker, Yun Wang, Risa H. Wechsler, Charity Woodrum, Ann Zabludoff, Dennis Zaritsky, Andrew R. Zentner, Huanian Zhang

Abstract: Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive bl… ▽ More Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Primordial Non-Gaussianity

Authors: P. Daniel Meerburg, Daniel Green, Muntazir Abidi, Mustafa A. Amin, Peter Adshead, Zeeshan Ahmed, David Alonso, Behzad Ansarinejad, Robert Armstrong, Santiago Avila, Carlo Baccigalupi, Tobias Baldauf, Mario Ballardini, Kevin Bandura, Nicola Bartolo, Nicholas Battaglia, Daniel Baumann, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Colin Bischoff, Jonathan Blazek, J. Richard Bond, Julian Borrill , et al. (153 additional authors not shown)

Abstract: Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with… ▽ More Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale. △ Less

Submitted 14 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

Science from an Ultra-Deep, High-Resolution Millimeter-Wave Survey

Authors: Neelima Sehgal, Ho Nam Nguyen, Joel Meyers, Moritz Munchmeyer, Tony Mroczkowski, Luca Di Mascolo, Eric Baxter, Francis-Yan Cyr-Racine, Mathew Madhavacheril, Benjamin Beringue, Gil Holder, Daisuke Nagai, Simon Dicker, Cora Dvorkin, Simone Ferraro, George M. Fuller, Vera Gluscevic, Dongwon Han, Bhuvnesh Jain, Bradley Johnson, Pamela Klaassen, Daan Meerburg, Pavel Motloch, David N. Spergel, Alexander van Engelen , et al. (44 additional authors not shown)

Abstract: Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey woul… ▽ More Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: 5 pages + references; Submitted to the Astro2020 call for science white papers

Impact of inhomogeneous reionization on the Lyman-α forest

Authors: Paulo Montero-Camacho, Christopher M. Hirata, Paul Martini, Klaus Honscheid

Abstract: The \lya forest at high redshifts is a powerful probe of reionization. Modeling and observing this imprint comes with significant technical challenges: inhomogeneous reionization must be taken into account while simultaneously being able to resolve the web-like small-scale structure prior to reionization. In this work we quantify the impact of inhomogeneous reionization on the \lya forest at lower… ▽ More The \lya forest at high redshifts is a powerful probe of reionization. Modeling and observing this imprint comes with significant technical challenges: inhomogeneous reionization must be taken into account while simultaneously being able to resolve the web-like small-scale structure prior to reionization. In this work we quantify the impact of inhomogeneous reionization on the \lya forest at lower redshifts ($2 < z < 4$), where upcoming surveys such as DESI will enable precision measurements of the flux power spectrum. We use both small box simulations capable of handling the small-scale structure of the \lya forest and semi-numerical large box simulations capable of representing the effects of inhomogeneous reionization. We find that inhomogeneous reionization could produce a measurable effect on the \lya forest power spectrum. The deviation in the 3D power spectrum at $z_{\rm obs} = 4$ and $k = 0.14 \ \rm{Mpc}^{-1}$ ranges from $19 - 36\%$, with a larger effect for later reionization. The corrections decrease to $2.0 - 4.1\%$ by $z_{\rm obs} = 2$. The impact on the 1D power spectrum is smaller, and ranges from $3.3 - 6.5\%$ at $z_{\rm obs}=4$ to $0.35 - 0.75\%$ at $z_{\rm obs}=2$, values which are comparable to the statistical uncertainties in current and upcoming surveys. Furthermore, we study how can this systematic be constrained with the help of the quadrupole of the 21 cm power spectrum. △ Less

Submitted 20 May, 2019; v1 submitted 7 February, 2019; originally announced February 2019.

Comments: New version matches the accepted version for publication. 11 pages, 5 figures

The Population of Eccentric Binary Black Holes: Implications for mHz Gravitational Wave Experiments

Authors: Xiao Fang, Todd A. Thompson, Christopher M. Hirata

Abstract: The observed binary black hole (BBH) mergers indicate a large Galactic progenitor population continuously evolving from large orbital separations and low gravitational wave (GW) frequencies to the final merger phase. We investigate the equilibrium distribution of binary black holes in the Galaxy. Given the observed BBH merger rate, we contrast the expected number of systems radiating in the low-fr… ▽ More The observed binary black hole (BBH) mergers indicate a large Galactic progenitor population continuously evolving from large orbital separations and low gravitational wave (GW) frequencies to the final merger phase. We investigate the equilibrium distribution of binary black holes in the Galaxy. Given the observed BBH merger rate, we contrast the expected number of systems radiating in the low-frequency $0.1-10\,$mHz GW band under two assumptions: (1) that all merging systems originate from near-circular orbits, as may be indicative of isolated binary evolution, and (2) that all merging systems originate at very high eccentricity, as predicted by models of dynamically-formed BBHs and triple and quadruple systems undergoing Lidov-Kozai eccentricity oscillations. We show that the equilibrium number of systems expected at every frequency is higher in the eccentric case (2) than in the circular case (1) by a factor of $\simeq 2-15$. This follows from the fact that eccentric systems spend more time than circular systems radiating in the low-frequency GW bands. The GW emission comes in pulses at periastron separated by the orbital period, which may be days to years. For a LISA-like sensitivity curve, we show that if eccentric systems contribute significantly to the observed merger rate, then $\simeq 10$ eccentric systems should be seen in the Galaxy. △ Less

Submitted 7 March, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Comments: 12 pages, 6 figures, matches ApJ accepted version

Journal ref: ApJ, 875, 75 (2019)

Time evolution of intrinsic alignments of galaxies

Authors: Denise M. Schmitz, Christopher M. Hirata, Jonathan Blazek, Elisabeth Krause

Abstract: Intrinsic alignments (IA), correlations between the intrinsic shapes and orientations of galaxies on the sky, are both a significant systematic in weak lensing and a probe of the effect of large-scale structure on galactic structure and angular momentum. In the era of precision cosmology, it is thus especially important to model IA with high accuracy. Efforts to use cosmological perturbation theor… ▽ More Intrinsic alignments (IA), correlations between the intrinsic shapes and orientations of galaxies on the sky, are both a significant systematic in weak lensing and a probe of the effect of large-scale structure on galactic structure and angular momentum. In the era of precision cosmology, it is thus especially important to model IA with high accuracy. Efforts to use cosmological perturbation theory to model the dependence of IA on the large-scale structure have thus far been relatively successful; however, extant models do not consistently account for time evolution. In particular, advection of galaxies due to peculiar velocities alters the impact of IA, because galaxy positions when observed are generally different from their positions at the epoch when IA is believed to be set. In this work, we evolve the galaxy IA from the time of galaxy formation to the time at which they are observed, including the effects of this advection, and show how this process naturally leads to a dependence of IA on the velocity shear. We calculate the galaxy-galaxy-IA bispectrum to tree level (in the linear matter density) in terms of the evolved IA coefficients. We then discuss the implications for weak lensing systematics as well as for studies of galaxy formation and evolution. We find that considering advection introduces nonlocality into the bispectrum, and that the degree of nonlocality represents the memory of a galaxy's path from the time of its formation to the time of observation. We discuss how this result can be used to constrain the redshift at which IA is determined and provide Fisher estimation for the relevant measurements using the example of SDSS-BOSS. △ Less

Submitted 16 July, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 30 pages, 5 figures, 2 tables

Journal ref: JCAP07(2018)030

The ultimate frontier of 21cm cosmology

Authors: Patrick C. Breysse, Yacine Ali-Haïmoud, Christopher M. Hirata

Abstract: We present the most detailed computation to date of the 21-cm global signal and fluctuations at $z\gtrsim 500$. Our calculations include a highly precise estimate of the Wouthuysen-Field (WF) effect and the first explicit calculation of the impact of free-free processes, the two dominant components of the signal at $z\gtrsim 800$. We implement a new high-resolution Ly$α$ radiative transfer calcula… ▽ More We present the most detailed computation to date of the 21-cm global signal and fluctuations at $z\gtrsim 500$. Our calculations include a highly precise estimate of the Wouthuysen-Field (WF) effect and the first explicit calculation of the impact of free-free processes, the two dominant components of the signal at $z\gtrsim 800$. We implement a new high-resolution Ly$α$ radiative transfer calculation, coupled to a state-of-the-art primordial recombination code. Using these tools, we find a global signal from 21-cm processes alone of roughly 0.01mK at $z\sim1000$, slightly larger than it would be without the WF effect, but much weaker than previous estimates including this effect. We also find that this signal is swamped by a smooth $1-2$ mK signal due to free-free absorption at high redshift by the partially ionized gas along the line of sight. In addition, we estimate the amplitude of 21-cm fluctuations, of order $\sim 10^{-7}$ mK at $z\sim1000$. Unfortunately, we find that due to the brightness of the low-frequency sky, these fluctuations will not be observable beyond $z\sim$ a few hundred by even extremely futuristic observations. The 21 cm fluctuations are exponentially suppressed at higher redshifts by the large free-free optical depth, making this the ultimate upper redshift limit for 21-cm surveys. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Comments: 14 pages, 9 figures, submitted to PRD

Journal ref: Phys. Rev. D 98, 043520 (2018)

Effects of NII and H$α$ Line Blending on the WFIRST Galaxy Redshift Survey

Authors: Daniel Martens, Xiao Fang, M. A. Troxel, Joe DeRose, Christopher M. Hirata, Risa H. Wechsler, Yun Wang

Abstract: The Wide Field Infrared Survey Telescope (WFIRST) will conduct a galaxy redshift survey using the H$α$ emission line primarily for spectroscopic redshift determination. Due to the modest spectroscopic resolution of the grism, the H$α$ and the neighboring [NII] lines are blended, leading to a redshift bias that depends on the [NII]/H$α$ ratio, which is correlated with a galaxy's metallicity, hence… ▽ More The Wide Field Infrared Survey Telescope (WFIRST) will conduct a galaxy redshift survey using the H$α$ emission line primarily for spectroscopic redshift determination. Due to the modest spectroscopic resolution of the grism, the H$α$ and the neighboring [NII] lines are blended, leading to a redshift bias that depends on the [NII]/H$α$ ratio, which is correlated with a galaxy's metallicity, hence mass and ultimately environment. We investigate how this bias propagates into the galaxy clustering and cosmological parameters obtained from the WFIRST. Using simulation, we explore the effect of line blending on redshift-space distortion and baryon acoustic oscillation (BAO) measurements. We measure the BAO parameters $α_{\parallel}$, $α_{\perp}$, the logarithmic growth factor $f_{v}$, and calculate their errors based on the correlations between the line ratio and large-scale structure. We find $Δα_{\parallel} = 0.31 \pm 0.23 \%$ ($0.26\pm0.17\%$), $Δα_{\perp} = -0.10\pm0.10\%$ ($-0.12 \pm 0.11 \%$), and $Δf_{v} = 0.17\pm0.33\%$ ($-0.20 \pm 0.30\%$) for redshift 1.355--1.994 (0.700--1.345), which use approximately 18$\%$, 9$\%$, and 7$\%$ of the systematic error budget in a root-sum-square sense. These errors may already be tolerable but further mitigations are discussed. Biases due to the environment-independent redshift error can be mitigated by measuring the redshift error probability distribution function. High-spectral-resolution re-observation of a few thousand galaxies would be required (if by direct approach) to reduce them to below 25$\%$ of the error budget. Finally, we outline the next steps to improve the modeling of [NII]-induced blending biases and their interaction with other redshift error sources. △ Less

Submitted 7 February, 2019; v1 submitted 21 April, 2018; originally announced April 2018.

Comments: 19 pages, 11 figures, minor changes, matched with MNRAS accepted version

Exploring circular polarization in the CMB due to conventional sources of cosmic birefringence

Authors: Paulo Montero-Camacho, Christopher M. Hirata

Abstract: The circular polarization of the cosmic microwave background (CMB) is usually taken to be zero since it is not generated by Thomson scattering. Here we explore the actual level of circular polarization in the CMB generated by conventional cosmological sources of birefringence. We consider two classes of mechanisms for birefringence. One is alignment of the matter to produce an anisotropic suscepti… ▽ More The circular polarization of the cosmic microwave background (CMB) is usually taken to be zero since it is not generated by Thomson scattering. Here we explore the actual level of circular polarization in the CMB generated by conventional cosmological sources of birefringence. We consider two classes of mechanisms for birefringence. One is alignment of the matter to produce an anisotropic susceptibility tensor: the hydrogen spins can be aligned either by density perturbations or CMB anisotropies themselves. The other is anisotropy of the radiation field coupled to the non-linear response of the medium to electromagnetic fields: this can occur either via photon-photon scattering (non-linear response of the vacuum); atomic hyperpolarizability (non-linear response of neutral atoms); or plasma delay (non-linear response of free electrons). The strongest effect comes from photon-photon scattering from recombination at a level of $\sim 10^{-14}$ K. Our results are consistent with a negligible circular polarization of the CMB in comparison with the linear polarization or the sensitivity of current and near-term experiments. △ Less

Submitted 21 August, 2018; v1 submitted 12 March, 2018; originally announced March 2018.

Comments: New version matches the accepted version for publication. 45 pages, 7 figures

Journal ref: JCAP08(2018)040

A Radial Measurement of the Galaxy Tidal Alignment Magnitude with BOSS Data

Authors: Daniel Martens, Christopher M. Hirata, Ashley J. Ross, Xiao Fang

Abstract: The anisotropy of galaxy clustering in redshift space has long been used to probe the rate of growth of cosmological perturbations. However, if galaxies are aligned by large-scale tidal fields, then a sample with an orientation-dependent selection effect has an additional anisotropy imprinted onto its correlation function. We use the LOWZ and CMASS catalogs of SDSS-III BOSS Data Release 12 to divi… ▽ More The anisotropy of galaxy clustering in redshift space has long been used to probe the rate of growth of cosmological perturbations. However, if galaxies are aligned by large-scale tidal fields, then a sample with an orientation-dependent selection effect has an additional anisotropy imprinted onto its correlation function. We use the LOWZ and CMASS catalogs of SDSS-III BOSS Data Release 12 to divide galaxies into two sub-samples based on their offset from the Fundamental Plane, which should be correlated with orientation. These sub-samples must trace the same underlying cosmology, but have opposite orientation-dependent selection effects. We measure the clustering parameters of each sub-sample and compare them in order to calculate the dimensionless parameter $B$, a measure of how strongly galaxies are aligned by gravitational tidal fields. We found that for CMASS (LOWZ), the measured $B$ was $-0.024 \pm 0.015$ ($-0.030 \pm 0.016$). This result can be compared to the theoretical predictions of Hirata 2009, who argued that since galaxy formation physics does not depend on the direction of the observer, the same intrinsic alignment parameters that describe galaxy-ellipticity correlations should also describe intrinsic alignments in the radial direction. We find that the ratio of observed to theoretical values is $0.51\pm 0.32$ ($0.77\pm0.41$) for CMASS (LOWZ). We combine the results to obtain a total ${\rm {Obs}/{Theory}} = 0.61\pm 0.26$. This measurement constitutes evidence (between 2 and 3$σ$) for radial intrinsic alignments, and is consistent with theoretical expectations ($<2σ$ difference). △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 23 pages, 18 figures, to be submitted to MNRAS

A New Window into Escaping Exoplanet Atmospheres: 10830 Å Line of Helium

Authors: Antonija Oklopčić, Christopher M. Hirata

Abstract: Observational evidence for escaping exoplanet atmospheres has been obtained for a few exoplanets to date. It comes from strong transit signals detected in the ultraviolet, most notably in the wings of the hydrogen Lyman-$α$ (Ly$α$) line. However, the core of the Ly$α$ line is often heavily affected by interstellar absorption and geocoronal emission, limiting the information about the atmosphere th… ▽ More Observational evidence for escaping exoplanet atmospheres has been obtained for a few exoplanets to date. It comes from strong transit signals detected in the ultraviolet, most notably in the wings of the hydrogen Lyman-$α$ (Ly$α$) line. However, the core of the Ly$α$ line is often heavily affected by interstellar absorption and geocoronal emission, limiting the information about the atmosphere that can be extracted from that part of the spectrum. Transit observations in atomic lines that are (a) sensitive enough to trace the rarefied gas in the planetary wind and (b) do not suffer from significant extinction by the interstellar medium could enable more detailed observations, and thus provide better constraints on theoretical models of escaping atmospheres. The absorption line of a metastable state of helium at 10830 Å could satisfy both of these conditions for some exoplanets. We develop a simple 1D model of escaping planetary atmospheres containing hydrogen and helium. We use it to calculate the density profile of helium in the 2$^3$S metastable excited state and the expected in-transit absorption at 10830 Å for two exoplanets known to have escaping atmospheres. Our results indicate that exoplanets similar to GJ 436b and HD 209458b should exhibit enhanced transit depths at 10830 Å, with $\sim 8\%$ and $\sim 2\%$ excess absorption in the line core, respectively. △ Less

Submitted 7 March, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: New version matches the published version

Journal ref: ApJL 855, L11 (2018)

Dynamics of Quadruple Systems Composed of Two Binaries: Stars, White Dwarfs, and Implications for Ia Supernovae

Authors: Xiao Fang, Todd A. Thompson, Christopher M. Hirata

Abstract: [Abridged]We investigate the long-term secular dynamics and Lidov-Kozai(LK) eccentricity oscillations of quadruple systems composed of two binaries at quadrupole and octupole order in the perturbing Hamiltonian. We show that the fraction of systems reaching high eccentricities is enhanced relative to triple systems, over a broader range of parameter space. We show that this fraction grows with tim… ▽ More [Abridged]We investigate the long-term secular dynamics and Lidov-Kozai(LK) eccentricity oscillations of quadruple systems composed of two binaries at quadrupole and octupole order in the perturbing Hamiltonian. We show that the fraction of systems reaching high eccentricities is enhanced relative to triple systems, over a broader range of parameter space. We show that this fraction grows with time, unlike triple systems evolved at quadrupole order. This is fundamentally because with their additional degrees of freedom, quadruple systems do not have a maximal set of commuting constants of the motion, even in secular theory at quadrupole order. We discuss these results in the context of star-star and white dwarf-white dwarf(WD) binaries, with emphasis on WD-WD mergers and collisions relevant to the Type Ia supernova problem. For star-star systems, we find that more than 30% of systems reach high eccentricity within a Hubble time, potentially forming triple systems via stellar mergers or close binaries. For WD-WD systems, taking into account general relativistic and tidal precession and dissipation, we show that the merger rate is enhanced in quadruple systems relative to triple systems by a factor of 3.5-10, and that the long-term evolution of quadruple systems leads to a delay-time distribution ~1/t for mergers and collisions. In GW-driven mergers of compact objects, we classify the mergers by their evolutionary patterns in phase space and identify a regime in ~8% of orbital shrinking mergers, where eccentricity oscillations occur on the general relativistic precession timescale, rather than the much longer LK timescale. Finally, we generalize previous treatments of oscillations in the inner binary eccentricity to eccentric mutual orbits. We assess the merger rate in quadruple and triple systems and the implications for their viability as progenitors of stellar mergers and SNe Ia. △ Less

Submitted 20 February, 2018; v1 submitted 25 September, 2017; originally announced September 2017.

Comments: 32 pages, 29 figures, matches MNRAS accepted version

Solar system science with the Wide-Field InfraRed Survey Telescope (WFIRST)

Authors: B. J. Holler, S. N. Milam, J. M. Bauer, C. Alcock, M. T. Bannister, G. L. Bjoraker, D. Bodewits, A. S. Bosh, M. W. Buie, T. L. Farnham, N. Haghighipour, P. S. Hardersen, A. W. Harris, C. M. Hirata, H. H. Hsieh, M. S. P. Kelley, M. M. Knight, E. A. Kramer, A. Longobardo, C. A. Nixon, E. Palomba, S. Protopapa, L. C. Quick, D. Ragozzine, V. Reddy , et al. (8 additional authors not shown)

Abstract: We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $μ$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are wher… ▽ More We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $μ$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations. △ Less

Submitted 27 July, 2018; v1 submitted 8 September, 2017; originally announced September 2017.

Comments: 58 pages, 14 figures, 7 tables

How does the shape of the stellar spectrum affect the Raman scattering features in the albedo of exoplanets?

Authors: Antonija Oklopčić, Christopher M. Hirata, Kevin Heng

Abstract: Diagnostic potential of the spectral signatures of Raman scattering imprinted in planetary albedo spectra at short optical wavelengths has been demonstrated in research on Solar System planets and has recently been proposed as a probe of exoplanet atmospheres, complementary to albedo studies at longer wavelengths. Spectral features caused by Raman scattering offer insight into the properties of pl… ▽ More Diagnostic potential of the spectral signatures of Raman scattering imprinted in planetary albedo spectra at short optical wavelengths has been demonstrated in research on Solar System planets and has recently been proposed as a probe of exoplanet atmospheres, complementary to albedo studies at longer wavelengths. Spectral features caused by Raman scattering offer insight into the properties of planetary atmospheres, such as the atmospheric depth, composition, and temperature, as well as the possibility of detecting and spectroscopically identifying spectrally inactive species, such as H$_2$ and N$_2$, in the visible wavelength range. Raman albedo features, however, depend on both the properties of the atmosphere and the shape of the incident stellar spectrum. Identical planetary atmospheres can produce very different albedo spectra depending on the spectral properties of the host star. Here we present a set of geometric albedo spectra calculated for atmospheres with H$_2$/He, N$_2$, and CO$_2$ composition, irradiated by different stellar types ranging from late A to late K stars. Prominent albedo features caused by Raman scattering appear at different wavelengths for different types of host stars. We investigate how absorption due to alkali elements, sodium and potassium, may affect the intensity of Raman features, and discuss the preferred observing strategies for detecting Raman features in future observations. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: 12 pages; accepted for publication in ApJ

Detecting primordial gravitational waves with circular polarization of the redshifted 21 cm line: II. Forecasts

Authors: Abhilash Mishra, Christopher M. Hirata

Abstract: In the first paper of this series, we showed that the CMB quadrupole at high redshifts results in a small circular polarization of the emitted 21 cm radiation. In this paper we forecast the sensitivity of future radio experiments to measure the CMB quadrupole during the era of first cosmic light ($z\sim 20$). The tomographic measurement of 21 cm circular polarization allows us to construct a 3D re… ▽ More In the first paper of this series, we showed that the CMB quadrupole at high redshifts results in a small circular polarization of the emitted 21 cm radiation. In this paper we forecast the sensitivity of future radio experiments to measure the CMB quadrupole during the era of first cosmic light ($z\sim 20$). The tomographic measurement of 21 cm circular polarization allows us to construct a 3D remote quadrupole field. Measuring the $B$-mode component of this remote quadrupole field can be used to put bounds on the tensor-to-scalar ratio $r$. We make Fisher forecasts for a future Fast Fourier Transform Telescope (FFTT), consisting of an array of dipole antennas in a compact grid configuration, as a function of array size and observation time. We find that a FFTT with a side length of 100 km can achieve $σ(r)\sim 4\times 10^{-3}$ after ten years of observation and with a sky coverage $f_{\mathrm{sky}}\sim 0.7$. The forecasts are dependent on the evolution of the Lyman-$α$ flux in the pre-reionization era, that remains observationally unconstrained. Finally, we calculate the typical order of magnitudes for circular polarization foregrounds and comment on their mitigation strategies. We conclude that detection of primordial gravitational waves with 21 cm observations is in principle possible, so long as the primordial magnetic field amplitude is small, but would require a very futuristic experiment with corresponding advances in calibration and foreground suppression techniques. △ Less

Submitted 27 April, 2018; v1 submitted 11 July, 2017; originally announced July 2017.

Comments: 19 pages, matches PRD accepted version

Journal ref: Phys. Rev. D 97, 103522 (2018)