Search

Recent reports of cosmological parity violation in the 4PCF raises the question of how such violations could be systematically generated. Here we present a constructive procedure to generate arbitrary violations of vectorial and tensorial types on any scale, which is computationally efficient in the squeezed limit. We directly compute their numerical transfer function, and find strong conservation in the linear regime. This procedure spans all squeezed parity violating observables at the 4PCF, following the quadratic estimator classification., Comment: 6 pages, 3 figures

The gravitational coupling between large- and small-scale density perturbations leads to anisotropic distortions to the local small-scale matter fluctuations. Such local anisotropic distortions can be used to reconstruct the large-scale matter distribution, known as tidal reconstruction. In this paper, we apply the tidal reconstruction methods to simulated galaxies in redshift space. We find that redshift-space distortions lead to anisotropic reconstruction results. While the reconstructed radial modes are more noisy mainly due to the small-scale velocity dispersion, the transverse modes are still reconstructed with high fidelity, well correlated with the original large-scale density modes. The bias of the reconstructed field at large scales shows a simple angular dependence which can be described by a form similar to the linear redshift-space distortion. The noise power spectrum is nearly isotropic and scale-independent on large scales. This makes the reconstructed tides field an ideal tracer for cosmic variance cancellation and multi-tracer analysis and has profound implications for future 21~cm intensity mapping surveys., Comment: 18 pages, 11 figures

The angular momentum of galaxies (galaxy spin) contains rich information about the initial condition of the Universe, yet it is challenging to efficiently measure the spin direction for the tremendous amount of galaxies that are being mapped by the ongoing and forthcoming cosmological surveys. We present a machine learning based classifier for the Z-wise vs S-wise spirals, which can help to break the degeneracy in the galaxy spin direction measurement. The proposed Chirality Equivariant Residual Network (CE-ResNet) is manifestly equivariant under a reflection of the input image, which guarantees that there is no inherent asymmetry between the Z-wise and S-wise probability estimators. We train the model with Sloan Digital Sky Survey (SDSS) images, with the training labels given by the Galaxy Zoo 1 (GZ1) project. A combination of data augmentation tricks are used during the training, making the model more robust to be applied to other surveys. We find a $\sim\!30\%$ increase of both types of spirals when Dark Energy Spectroscopic Instrument (DESI) images are used for classification, due to the better imaging quality of DESI. We verify that the $\sim\!7\sigma$ difference between the numbers of Z-wise and S-wise spirals is due to human bias, since the discrepancy drops to $<\!1.8\sigma$ with our CE-ResNet classification results. We discuss the potential systematics that are relevant to the future cosmological applications., Comment: 13+4 pages, 11 figures, 2 tables, accepted by ApJ

Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to discover and localize a large sample of rare, high-fluence, nearby FRBs. This is the population most amenable to multi-messenger, multi-wavelength follow-up, allowing deeper understanding of source mechanisms. BURSTT will provide horizon-to-horizon sky coverage with a half power field-of-view (FoV) of $\sim$10$^{4}$ deg$^{2}$, a 400 MHz effective bandwidth between 300-800 MHz, and sub-arcsecond localization, made possible using outrigger stations hundreds to thousands of km from the main array. Initially, BURSTT will employ 256 antennas. After tests of various antenna designs and optimization of system performance we plan to expand to 2048 antennas. We estimate that BURSTT-256 will detect and localize $\sim$100 bright ($\geq$100 Jy ms) FRBs per year. Another advantage of BURSTT's large FoV and continuous operation will be greatly enhanced monitoring of FRBs for repetition. The current lack of sensitive all-sky observations likely means that many repeating FRBs are currently cataloged as single-event FRBs., Comment: Published version

The reconstruction method was proposed more than a decade ago to boost the signal of baryonic acoustic oscillations measured in galaxy redshift surveys, which is one of key probes for dark energy. After moving the observed overdensities in galaxy surveys back to their initial position, the reconstructed density field is closer to a linear Gaussian field, with higher-order information moved back into the power spectrum. We find that by jointly analysing power spectra measured from the pre- and post-reconstructed galaxy samples, higher-order information beyond the $2$-point power spectrum can be efficiently extracted, which generally yields an information gain upon the analysis using the pre- or post-reconstructed galaxy sample alone. This opens a window to easily use higher-order information when constraining cosmological models., Comment: 38 pages including 4 figures and 1 table, with a Supplement Information document attached; Version accepted to Communications Physics

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

Potential features in the primordial power spectrum have been searched for in galaxy surveys in recent years since these features can assist in understanding the nature of inflation. The null detection to date suggests that any such features should be fairly weak, and next-generation galaxy surveys, with their unprecedented sizes and precisions, are in a position to place stronger constraints than before. However, even if such primordial features once existed in the early Universe, they would have been significantly damped in the nonlinear regime at low redshift due to structure formation, which makes them difficult to be directly detected in real observations. A potential way to tackle this challenge for probing the features is to undo the cosmological evolution, i.e., using reconstruction to obtain an approximate linear density field. By employing a set of N-body simulations, here we show that a recently-proposed nonlinear reconstruction algorithm can effectively retrieve damped oscillatory features from halo catalogues and improve the accuracy of the measurement of feature parameters (assuming that such primordial features do exist). We do a Fisher analysis to forecast how nonlinear reconstruction affects the constraining power, and find that it can lead to significantly more robust constraints on the feature amplitude for a DESI-like survey. Comparing nonlinear reconstruction with other ways of improving constraints, such as increasing the survey volume and range of scales, this shows that it is possible to achieve what the latter do, but at a lower cost.

We present the optimal convergence and shear estimators for lensing reconstruction from the cosmic microwave background temperature field. This generalizes the deflection estimator, is sensitive to non-lensing modes, provides internal consistency checks, and is always at least as optimal. Previously, these estimators were only known in the squeezed limit. This paper decomposes convergence and shear fields into cosine and sine waves and the lensed correlation function is then Taylor expanded in the wave amplitudes. Maximizing the likelihood function gives the optimal estimators for the convergence and shear fields. This method has the potential to improve the lensing reconstruction of the cosmic microwave background polarization field: the shear and convergence can be optimally combined to form a deflection estimator, or used separately to separate non-lensing modes, or utilize lensing of non-Gaussian secondary foregrounds., Comment: 9 pages, 3 figures

Galaxy-galaxy strong gravitational lenses have become a popular probe of dark matter (DM) by providing a window into structure formation on the smallest scales. In particular, the convergence power spectrum of subhalos within lensing galaxies has been suggested as a promising observable to study DM. However, the distances involved in strong-lensing systems are vast, and we expect the relevant volume to contain line-of-sight (LOS) halos that are not associated with the main lens. We develop a formalism to calculate the effect of LOS halos as an effective convergence power spectrum. The multi-lens plane equation couples the angular deflections of consecutive lens planes, but by assuming that the perturbations due to the LOS halos are small, we show that they can be projected onto the main-lens plane as effective subhalos. We test our formalism by simulating lensing systems using the full multi-plane lens equation and find excellent agreement. We show how the relative contribution of LOS halos and subhalos depends on the source and lens redshift, as well as the assumed halo and subhalo mass functions. For a fiducial system with fraction of DM halo mass in substructure $f_{\rm sub}=0.4\%$ for subhalo masses $[10^5-10^8]\rm{M}_{\odot}$, the interloper contribution to the power spectrum is at least several times greater than that of subhalos for source redshifts $z_s\gtrsim0.5$. Furthermore, it is likely that for the SLACS and BELLS lenses the interloper contribution dominates: $f_{\rm sub}\gtrsim2\%$ ($4\%$) is needed for subhalos to dominate in SLACS (BELLS), which is higher than current upper bounds on $f_{\rm sub}$ for our mass range. Since the halo mass function is better understood from first principles, the dominance of interlopers in galaxy-galaxy lenses with high-quality imaging can be seen as a significant advantage when translating this observable into a constraint on DM., Comment: 24 pages, 8 figures, Matches accepted version, v3: Fixed a typo in Eq. C4 in the appendix. No changes in the derivations

Reconstruction is becoming a crucial procedure of galaxy clustering analysis for future spectroscopic redshift surveys to obtain subpercent level measurement of the baryon acoustic oscillation scale. Most reconstruction algorithms rely on an estimation of the displacement field from the observed galaxy distribution. However, the displacement reconstruction degrades near the survey boundary due to incomplete data and the boundary effects extend to $\sim100\,\mathrm{Mpc}/h$ within the interior of the survey volume. We study the possibility of using radial velocities measured from the cosmic microwave background observation through the kinematic Sunyaev-Zel'dovich effect to improve performance near the boundary. We find that the boundary effect can be reduced to $\sim30-40\,\mathrm{Mpc}/h$ with the velocity information from Simons Observatory. This is especially helpful for dense low redshift surveys where the volume is relatively small and a large fraction of total volume is affected by the boundary., Comment: 12 pages, 14 figures, in press

We propose a new velocity reconstruction method based on the displacement estimation by recently developed methods. The velocity is first reconstructed by transfer functions in Lagrangian space and then mapped into Eulerian space. High resolution simulations are used to test the performance. We find that the new reconstruction method outperforms the standard velocity reconstruction in the sense of better cross-correlation coefficient, less velocity misalignment and smaller amplitude difference. We conclude that this new method has the potential to improve the large-scale structure sciences involving a velocity reconstruction, such as kinetic Sunyaev-Zel'dovich measurement and supernova cosmology., Comment: 14 pages, 12 figures, matched to the published version

We present a new method to measure neutrino masses using the dark matter-neutrino relative velocity. The relative motion between dark matter and neutrinos results in a parity-odd bispectrum which can be measured from cross-correlation of different cosmic fields. This new method is not affected by most systematics which are parity even and not limited by the knowledge of optical depth to the cosmic microwave background. We estimate the detectability of the relative velocity effect and find that the minimal sum of neutrino masses could be detected at high significance with upcoming surveys., Comment: 6 pages, 2 figures, 1 table, published version

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space., Comment: Science white paper submitted to the Astro2020 Decadal Survey

We discuss the possibility of studying diffuse baryon distributions with kinematic Sunyaev-Zel'dovich (kSZ) effect by correlating cosmic microwave background (CMB) temperature fluctuations with density fluctuations from 21\,cm intensity mapping (IM). The biggest challenge for the cross-correlation is the loss of large-scale information in IM, due to foregrounds and the zero spacing problem of interferometers. We apply the tidal reconstruction algorithm to restore the lost large-scale modes, which increases the correlation by more than a factor of three. With the predicted foreground level, we expect a $\sim20\,\sigma$ detection of kSZ signal for $0.8\lesssim z\lesssim2.5$ with CHIME and Planck, and a $\sim40\, \sigma$ detection with HIRAX and Planck. The significance can be greatly increased with next-generation facilities of higher spatial resolutions., Comment: 8 pages, 4 figures, submit to PRD

We apply nonlinear reconstruction to the dark matter density field in redshift space and solve for the nonlinear mapping from the initial Lagrangian position to the final redshift space position. The reconstructed anisotropic field inferred from the nonlinear displacement correlates with the linear initial conditions to much smaller scales than the redshift space density field. The number of linear modes in the density field is improved by a factor of 30-40 after reconstruction. We thus expect this reconstruction approach to substantially expand the cosmological information including baryon acoustic oscillations and redshift space distortions for dense low-redshift large scale structure surveys including for example SDSS main sample, DESI BGS, and 21 cm intensity mapping surveys., Comment: 18 pages, 21 figures, published version. The nonlinear reconstruction code is available at https://github.com/ColdThunder/NR-code

In this paper, we report a significant recovery of the linear baryonic acoustic oscillation (BAO) signature by applying the isobaric reconstruction algorithm to the non-linear matter density field. Assuming only the longitudinal component of the displacement being cosmologically relevant, this algorithm iteratively solves the coordinate transform between the Lagrangian and Eulerian frames without requiring any specific knowledge of the dynamics. For dark matter field, it produces the non-linear displacement potential with very high fidelity. The reconstruction error at the pixel level is within a few percent, and is caused only by the emergence of the transverse component after the shell-crossing. As it circumvents the strongest non-linearity of the density evolution, the reconstructed field is well-described by linear theory and immune from the bulk-flow smearing of the BAO signature. Therefore this algorithm could significantly improve the measurement accuracy of the sound horizon scale. For a perfect large-scale structure survey at redshift zero without Poisson or instrumental noise, the fractional error is reduced by a factor of 2.7, very close to the ideal limit with linear power spectrum and Gaussian covariance matrix., Comment: 5 pages, 3 figures, accepted version

We apply the nonlinear reconstruction method to simulated halo fields. For halo number density $2.77\times 10^{-2}$ $(h^{-1} {\rm Mpc})^{-3}$ at $z=0$, corresponding to the SDSS main sample density, we find the scale where the noise saturates the linear signal is improved to $k\gtrsim0.36\ h {\rm Mpc}^{-1}$, a factor of $2.29$ improvement in scale, or $12$ in number of linear modes. The improvement is less for higher redshift or lower halo density. We expect this to substantially improve the BAO accuracy of dense, low redshift surveys, including the SDSS main sample, 6dFGS and 21cm intensity mapping initiatives., Comment: 10 pages, 7 figures, published on ApJ. Match to published version. Two subsections and two figures are added

We present a direct approach to nonparametrically reconstruct the linear density field from an observed nonlinear map. We solve for the unique displacement potential consistent with the nonlinear density and positive definite coordinate transformation using a multigrid algorithm. We show that we recover the linear initial conditions up to the nonlinear scale ($r_{\delta_r\delta_L}>0.5$ for $k\lesssim1\ h/\mathrm{Mpc}$) with minimal computational cost. This reconstruction approach generalizes the linear displacement theory to fully nonlinear fields, potentially substantially expanding the baryon acoustic oscillations and redshift space distortions information content of dense large scale structure surveys, including for example SDSS main sample and 21cm intensity mapping initiatives., Comment: 7 pages, 7 figures, published version

Constraining neutrino mass remains an elusive challenge in modern physics. Precision measurements are expected from several upcoming cosmological probes of large-scale structure. Achieving this goal relies on an equal level of precision from theoretical predictions of neutrino clustering. Numerical simulations of the non-linear evolution of cold dark matter and neutrinos play a pivotal role in this process. We incorporate neutrinos into the cosmological N-body code CUBEP3M and discuss the challenges associated with pushing to the extreme scales demanded by the neutrino problem. We highlight code optimizations made to exploit modern high performance computing architectures and present a novel method of data compression that reduces the phase-space particle footprint from 24 bytes in single precision to roughly 9 bytes. We scale the neutrino problem to the Tianhe-2 supercomputer and provide details of our production run, named TianNu, which uses 86% of the machine (13,824 compute nodes). With a total of 2.97 trillion particles, TianNu is currently the world's largest cosmological N-body simulation and improves upon previous neutrino simulations by two orders of magnitude in scale. We finish with a discussion of the unanticipated computational challenges that were encountered during the TianNu runtime., Comment: Accepted for publication in Research in Astronomy and Astrophysics

Measurements of neutrino mass in cosmological observations rely on two point statistics that are hindered by significant degeneracies with the optical depth and galaxy bias. The relative velocity effect between cold dark matter and neutrinos induces a large scale dipole into the matter density field and may be able to provide orthogonal constraints to standard techniques. We numerically investigate this dipole in the TianNu Simulation, which contains cold dark matter and 50 meV neutrinos. We first compute the dipole using a new linear response technique where we treat the displacement caused by the relative velocity as a phase in Fourier space and then integrate the matter power spectrum over redshift. Then, we compute the dipole numerically in real space using the simulation density and velocity fields. We find excellent agreement between the linear response and N-body methods. Utilizing the dipole as an observational tool will require two tracers of the matter distribution that are differently biased with respect to the neutrino density., Comment: 6 pages, 5 figures

We study the nonlinear $E$-mode clustering in Lagrangian space by using large scale structure $N$-body simulations and use the displacement field information in Lagrangian space to recover the primordial linear density field. We find that, compared to Eulerian nonlinear density fields, the $E$-mode displacement fields in Lagrangian space improves the cross-correlation scale $k$ with initial density field by a factor of 6-7, containing two orders of magnitude more primordial information. This illustrates ability of potential density reconstruction algorithms, to improve the baryonic acoustic oscillation measurements from current and future large scale structure surveys., Comment: 5 pages, 4 figures

21 cm intensity mapping has emerged as a promising technique to map the large-scale structure of the Universe, at redshifts $z$ from 1 to 10. Unfortunately, many of the key cross-correlations with the cosmic microwave background and photo-$z$ galaxies have been thought to be impossible due to the foreground contamination for radial modes with small wave numbers. In this paper, we apply tidal reconstruction to the simulated 21 cm fields and recover the lost large-scale radial modes successfully. We estimate the detectability of the cross-correlation signals and find they can be detected at high significance with current 21 cm experiments. The tidal field reconstruction method opens up a new set of possibilities to probe the Universe and is extremely valuable not only for 21 cm surveys but also for cosmic microwave background and photometric-redshift observations., Comment: 9 pages, 6 figures, published version

Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys., Comment: 15 pages, 6 figures

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

The gravitational coupling of a long-wavelength tidal field with small-scale density fluctuations leads to anisotropic distortions of the locally measured small-scale matter correlation function. Since the local correlation function is known to be statistically isotropic in the absence of such tidal interactions, the tidal distortions can be used to reconstruct the long-wavelength tidal field and large-scale density field in analogy with the cosmic microwave background lensing reconstruction. In this paper we present the theoretical framework of cosmic tidal reconstruction and test the reconstruction in numerical simulations. We find that the density field on large scales can be reconstructed with good accuracy and the cross-correlation coefficient between the reconstructed density field and the original density field is greater than 0.9 on large scales $(k\lesssim0.1\ h/\mathrm{Mpc})$, with the filter scale $\sim1.25\ \mathrm{Mpc}/h$. This is useful in the 21cm intensity mapping survey, where the long-wavelength radial modes are lost due to a foreground subtraction process., Comment: 12 pages, 8 figures, published version

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

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