Search

Active galactic nuclei (AGN) are the signposts of black hole growth, and likely play an important role in galaxy evolution. An outstanding question is whether AGN of different spectral types indicate different evolutionary stages in the coevolution of black holes and galaxies. We present the angular correlation function between an AGN sample selected from the Hyper Suprime Camera Subaru Strategic Program (HSC-SSP) optical + Wide-field Infrared Survey Explorer (WISE) mid-IR photometry, and a luminous red galaxy (LRG) sample from HSC-SSP. We investigate AGN clustering strength as a function of their luminosity and spectral features across three independent HSC fields totaling $\sim600\,{\rm deg^{2}}$, for $z\in0.6-1.2$ and AGN with $L_{6\mu m}>3\times10^{44}{\rm\,erg\,s^{-1}}$. There are $\sim28,500$ AGN and $\sim1.5$ million LRGs in our primary analysis. We determine the inferred average halo mass for the full AGN sample ($M_h \approx 10^{12.9}h^{-1}M_\odot$), and note that it does not evolve significantly as a function of redshift (over this narrow range) or luminosity. We find that, on average, unobscured AGN ($M_h \approx10^{13.3}h^{-1}M_\odot$) occupy $\sim4.5\times$ more massive halos than obscured AGN ($M_h \approx10^{12.6}h^{-1}M_\odot$), at $5\sigma$ statistical significance using 1-D uncertainties, and at $3\sigma$ using the full covariance matrix, suggesting a physical difference between unobscured and obscured AGN, beyond the line-of-sight viewing angle. Furthermore, we find evidence for a halo mass dependence on reddening level within the Type I AGN population, which could support the existence of a previously claimed dust-obscured phase in AGN-host galaxy coevolution. However, we also find that even quite small systematic shifts in the redshift distributions of the AGN sample could plausibly explain current and previously observed differences in $M_{h}$., Comment: 31 pages, 14 figures, submitted to ApJ

The sum of cosmic neutrino masses can be measured cosmologically, as the sub-eV particles behave as `hot' dark matter whose main effect is to suppress the clustering of matter compared to a universe with the same amount of purely cold dark matter. Current astronomical data provide an upper limit on $\Sigma m_{\nu}$ between 0.07 - 0.12 eV at 95% confidence, depending on the choice of data. This bound assumes that the cosmological model is $\Lambda$CDM, where dark energy is a cosmological constant, the spatial geometry is flat, and the primordial fluctuations follow a pure power-law. Here, we update studies on how the mass limit degrades if we relax these assumptions. To existing data from the Planck satellite we add new gravitational lensing data from the Atacama Cosmology Telescope, the new Type Ia Supernova sample from the Pantheon+ survey, and baryonic acoustic oscillation (BAO) measurements from the Sloan Digital Sky Survey and the Dark Energy Spectrosopic Instrument. We find the neutrino mass limit is stable to most model extensions, with such extensions degrading the limit by less than 10%. We find a broadest bound of $\Sigma m_{\nu} < 0.19 ~\rm{eV}$ at 95% confidence for a model with dynamical dark energy, although this scenario is not statistically preferred over the simpler $\Lambda$CDM model., Comment: 8 pages

The large-angular-scale falloff in the autocorrelation function for the cosmic microwave background (CMB) temperature has long intrigued cosmologists and fueled speculation about suppressed superhorizon power. Here we highlight an inconsistency between the temperature quadrupole and the more recently obtained E-mode polarization quadrupole from Planck PR3. The temperature quadrupole arises primarily at the CMB surface of last scatter, while the polarization primarily from the epoch of reionization, but the two still probe comparable distance scales. Although the temperature quadrupole is intriguingly low (much greater than a $1\sigma$ fluctuation) compared with that expected in the standard $\Lambda$CDM cosmological model, the polarization quadrupole turns out to be somewhat high, at the $1\sigma$ level. We calculate the joint probability distribution function for both and find a slight tension: the observed pair of quadrupoles is inconsistent at a $2.3\sigma$ confidence level. The problem is robust to simple changes to the cosmological model. If the high polarization quadrupole survives further scrutiny, then this result disfavors, at comparable significance, new superhorizon physics. The full-sky coverage and pristine foreground subtraction of the LiteBIRD satellite will be ideal to help resolve this question., Comment: 7 pages with references, 4 figures

Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB., Comment: The Compton-y map and associated products will be made publicly available upon publication of the paper. The CMB T and E mode maps will be made available when the DR6 maps are made public

We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $\sigma_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equiv\sigma_8({\Omega_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $\sigma_8 = 0.812 \pm 0.013$, $S_8\equiv\sigma_8({\Omega_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $\Lambda$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$\sigma$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $\Lambda$CDM, limiting the sum of the neutrino masses to $\sum m_{\nu} < 0.13$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $\Lambda$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys., Comment: 32 pages, 17 figures, replaced with version accepted in ApJ (Feb 2024). Cosmological likelihood data and mass maps are public here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Qu et al and MacCrann et al

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., Comment: 25 pages, 20 figures. Submitted to MNRAS

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\'e 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., Comment: 28 pages, 19 figures, matches version accepted by MNRAS

Three-dimensional correlations of the Lyman-$\alpha$ (Ly$\alpha$) forest and cross correlations between the Ly$\alpha$ forest and quasars have been measured on large scales, allowing a precise measurement of the baryon acoustic oscillation (BAO) feature at redshifts $z>2$. These 3D correlations are often modelled using linear perturbation theory, but full-shape analyses to extract cosmological information beyond BAO will require more realistic models capable of describing non-linearities present at smaller scales. We present a measurement of the Ly$\alpha$ forest flux power spectrum from large hydrodynamic simulations -- the Sherwood simulations -- and compare it to different models describing the small-scale deviations from linear theory. We confirm that the model presented in Arinyo-i-Prats et al. (2015) fits the measured 3D power up to $k=10\, h\rm{Mpc^{-1}}$ with an accuracy better than 5%, and show that the same model can also describe the 1D correlations with similar precision. We also present, for the first time, an equivalent study for the cross-power spectrum of halos with the Ly$\alpha$ forest, and we discuss different challenges we face when modelling the cross-power spectrum beyond linear scales. We make all our measured power spectra public in \url{https://github.com/andreufont/sherwood_p3d}. This study is a step towards joint analyses of 1D and 3D flux correlations, and towards using the quasar-Ly$\alpha$ cross-correlation beyond BAO analyses., Comment: 31 pages, 19 figures

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 $\mu$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 $\xi_+ \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., Comment: 26 pages, 9 figures

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\sigma$ per $\Delta z = 0.5$ bin)., Comment: 12 pages, 5 figures, accepted for publication in MNRAS

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.

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-$\alpha$ 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-$\alpha$ 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-$\alpha$ 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., Comment: 25 pages, 5 figures. This version is the closest match to the PRD final print