Determining the large-scale distribution of baryons in the late universe is a long-standing challenge in cosmology. To gain insight into this problem, we present a new approach for extracting the kinematic Sunyaev-Zel'dovich (kSZ) effect from observations, ARF-kSZ tomography. This technique involves the cross-correlation of maps of Angular Redshift Fluctuations (ARF), which contain precise information about the cosmic density and velocity fields, and cosmic microwave background (CMB) temperature maps high-pass filtered using aperture photometry. To produce the first and second, in this work we resort to galaxies and quasars from 6dF and SDSS and foreground-cleaned CMB maps from Planck, respectively. We detect statistically significant cross-correlation between ARF and filtered CMB maps for a wide range of redshifts and filter apertures, yielding a joint detection of the kSZ effect at the $>10\sigma$ level. Using measurements of the kSZ optical depth extracted from these cross-correlations, we then set constraints on the properties of the gas responsible for the kSZ effect, finding that the kSZ gas resides mostly outside haloes and presents densities from 10 to 250 times the cosmic average, which is the density of baryons in filaments and sheets according to cosmological hydrodynamical simulations. Finally, we conduct a tomographic census of baryons from redshift z=0 to 5, finding that ARF-kSZ tomography is sensitive to approximately half of the baryons in the Universe.

We analytically model the non-linear effects induced by massive neutrinos on the total matter power spectrum using the halo model reaction framework of Cataneo et al. 2019. In this approach the halo model is used to determine the relative change to the matter power spectrum caused by new physics beyond the concordance cosmology. Using standard fitting functions for the halo abundance and the halo mass-concentration relation, the total matter power spectrum in the presence of massive neutrinos is predicted to percent-level accuracy, out to $k=10 \, h \, {\rm Mpc}^{-1}$. We find that refining the prescriptions for the halo properties using $N$-body simulations improves the recovered accuracy to better than 1%. This paper serves as another demonstration for how the halo model reaction framework, in combination with a single suite of standard $\Lambda$CDM simulations, can recover percent-level accurate predictions for beyond-$\Lambda$CDM matter power spectra, well into the non-linear regime.

A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We present results from the first cosmological simulation using Conservative Reproducing Kernel Smoothed Particle Hydrodynamics (CRK-SPH). We employ two simulations: one evolved purely under gravity and the other with non-radiative hydrodynamics. Each contains 2x2304^3 cold dark matter plus baryon particles in an 800 Mpc/h box. We compare statistics to previous non-radiative simulations including power spectra, mass functions, baryon fractions, and concentration. We find self-similar radial profiles of gas temperature, entropy, and pressure and show that a simple analytic model recovers these results to better than 40% over two orders of magnitude in mass. We quantify the level of non-thermal pressure support in halos and demonstrate that hydrostatic mass estimates are biased low by 24% (10%) for halos of mass 10^15 (10^13) Msun/h. We compute angular power spectra for the thermal and kinematic Sunyaev-Zel'dovich effects and find good agreement with the low-l Planck measurements. Finally, artificial scattering between particles of unequal mass is shown to have a large impact on the gravity-only run and we highlight the importance of better understanding this issue in hydrodynamic applications. This is the first in a simulation campaign using CRK-SPH with future work including subresolution gas treatments.

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.

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.

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.

We carry out a detailed study of the orbital dynamics and structural evolution of over 6000 subhalos in the Via Lactea II simulation, from infall to present. By analyzing subhalos with masses down to m = 4e5 Msun, we find that lower mass subhalos, which are not strongly affected by dynamical friction, exhibit behaviors qualitatively different from those found previously for more massive ones. Furthermore, there is a clear trend of subhalos that fell into the host earlier being less concentrated. We show that the concentration at infall characterizes various aspects of subhalo evolution. In particular, tidal effects truncate the growth of less concentrated subhalos at larger distances from the host; subhalos with smaller concentrations have larger infall radii. The concentration at infall is further shown to be a determining factor for the subsequent mass loss of subhalos within the host, and also for the evolution of their internal structure in the v_max-r_max plane. Our findings raise the prospects of using the concentration to predict the tidal evolution of subhalos, which will be useful for obtaining analytic models of galaxy formation, as well as for near field cosmology.

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

Early in the reionization process, the intergalactic medium (IGM) would have been quite inhomogeneous on small scales, due to the low Jeans mass in the neutral IGM and the hierarchical growth of structure in a cold dark matter Universe. This small-scale structure acted as an important sink during the epoch of reionization, impeding the progress of the ionization fronts that swept out from the first sources of ionizing radiation. Here we present results of high-resolution cosmological hydrodynamics simulations that resolve the cosmological Jeans mass of the neutral IGM in representative volumes several Mpc across. The adiabatic hydrodynamics we follow are appropriate in an unheated IGM, before the gas has had a chance to respond to the photoionization heating. Our focus is determination of the resolution required in cosmological simulations in order to sufficiently sample and resolve small-scale structure regulating the opacity of an unheated IGM. We find that a dark matter particle mass of m_dm < 50 Msun and box size of L > 1 Mpc are required. With our converged results we show how the mean free path of ionizing radiation and clumping factor of ionized hydrogen depends upon the ultraviolet background (UVB) flux and redshift. We find, for example at z = 10, clumping factors typically of 10 to 20 for an ionization rate of Gamma ~ 0.3 - 3 x 1e-12 s^-1, with corresponding mean free paths of ~ 3 - 15 Mpc, extending previous work on the evolving mean free path to considerably smaller scales and earlier times.

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