Marcel Schmittfull, Tobias Baldauf, Matias Zaldarriaga

Motivated by recent developments in perturbative calculations of the nonlinear evolution of large-scale structure, we present an iterative algorithm to reconstruct the initial conditions in a given volume starting from the dark matter distribution in real space. In our algorithm, objects are first iteratively moved back along estimated potential gradients until a nearly uniform catalog is obtained. The linear initial density is then estimated as the divergence of the cumulative displacement, with an optional second-order correction. This algorithm should undo non-linear effects up to one-loop order, including the higher-order infrared resummation piece. We test the method using dark matter simulations in real space. At redshift $z=0$, we find that after eight iterations the reconstructed density is more than $95\%$ correlated with the initial density at $k\le 0.35\; h\mathrm{Mpc}^{-1}$. The reconstruction also reduces the power in the difference between reconstructed and initial fields by more than two orders of magnitude at $k\le 0.2\; h\mathrm{Mpc}^{-1}$, and it extends the range of scales where the full broad-band shape of the power spectrum matches linear theory by a factor 2-3. As a specific application, we consider measurements of the Baryonic Acoustic Oscillation (BAO) scale that can be improved by reducing the degradation effects of large-scale flows. We find that the method improves the BAO signal-to-noise by a factor 2.7 at redshift $z=0$ and by a factor 2.5 at $z=0.6$ in our idealistic simulations. This improves standard BAO reconstruction by $70\%$ at $z=0$ and $30\%$ at $z=0.6$, and matches the optimal BAO signal and signal-to-noise of the linear density in the same volume.

Lachlan Lancaster, Francis-Yan Cyr-Racine, Lloyd Knox, Zhen Pan

We present updated constraints on the free-streaming nature of cosmological neutrinos from cosmic microwave background (CMB) power spectra, baryonic acoustic oscillation data, and local measurements of the Hubble constant. Specifically, we consider a Fermi-like four-fermion interaction between massless neutrinos, characterized by an effective coupling constant $ G_{\rm eff}$, and resulting in a neutrino opacity $\dot{\tau}_\nu\propto G_{\rm eff}^2 T_\nu^5$. Using a conservative prior on the parameter $\log_{10}\left(G_{\rm eff} {\rm MeV}^2\right)$, we find a bimodal posterior distribution. The first of these modes is consistent with the standard $\Lambda$CDM cosmology and corresponds to neutrinos decoupling at redshift $z_{\nu,{\rm dec}} > 1.3\times10^5$. The other mode of the posterior, dubbed the "interacting neutrino mode", corresponds to neutrino decoupling occurring within a narrow redshift window centered around $z_{\nu,{\rm dec}}\sim8300$. This mode is characterized by a high value of the effective neutrino coupling constant, together with a lower value of the scalar spectral index and amplitude of fluctuations, and a higher value of the Hubble parameter. Using both a maximum likelihood analysis and the ratio of the two mode's Bayesian evidence, we find the interacting neutrino mode to be statistically disfavored compared to the standard $\Lambda$CDM cosmology. Interestingly, the addition of CMB polarization and direct Hubble constant measurements significantly raises the statistical significance of this secondary mode, indicating that new physics in the neutrino sector could help explain the difference between local measurements of $H_0$, and those inferred from CMB data. A robust consequence of our results is that neutrinos must be free streaming long before the epoch of matter-radiation equality.

Benjamin M. Roberts, Geoffrey Blewitt, Conner Dailey, Mac Murphy, Maxim Pospelov, Alex Rollings, Jeff Sherman, Wyatt Williams, Andrei Derevianko

Multiple cosmological observations indicate that dark matter (DM) constitutes 85% of all matter in the Universe [1]. All the current evidence for DM comes from galactic or larger scale observations through the gravitational pull of DM on ordinary matter [1], leaving the microscopic nature of DM a mystery. Ambitious programs in particle physics have mostly focused on searches for WIMPs (Weakly Interacting Massive Particles) as DM candidates [2]. As WIMPs remain elusive, there is a growing interest in alternatives. Some models [3-7] predict DM in the form of spatially large objects ("clumps") that may cause glitches in atomic clocks [6]. Here we use the network of atomic clocks on board the GPS satellites as a 50,000 km aperture detector to search for DM clumps. As DM clumps sweep through the GPS satellite constellation at galactic velocities ~300 km/s, their predicted signature is a correlated propagation of clock glitches through the constellation over a period of a few minutes [6]. By mining 16 years of archival GPS data, we find no evidence for DM clumps in the form of domain walls. This enables us to improve limits on DM couplings to atomic clocks by several orders of magnitude. Our work demonstrates the use of a global network of precision measurement devices in the search for DM. Several global networks including magnetometers, laboratory atomic clocks, and other precision devices are being developed [5,8,9]. We anticipate that our methods will be valuable for probing new physics in this emerging area.

Abhirup Ghosh, Nathan K. Johnson-McDaniel, Archisman Ghosh, Chandra Kant Mishra, Parameswaran Ajith, Walter Del Pozzo, Christopher P. L. Berry, Alex B. Nielsen, Lionel London

Advanced LIGO's recent observations of gravitational waves (GWs) from merging binary black holes have opened up a unique laboratory to test general relativity (GR) in the highly relativistic regime. One of the tests used to establish the consistency of the first LIGO event with a binary black hole merger predicted by GR was the inspiral-merger-ringdown consistency test. This involves inferring the mass and spin of the remnant black hole from the inspiral (low-frequency) part of the observed signal and checking for the consistency of the inferred parameters with the same estimated from the post-inspiral (high-frequency) part of the signal. Based on the observed rate of binary black hole mergers, we expect the advanced GW observatories to observe hundreds of binary black hole mergers every year when operating at their design sensitivities, most of them with modest signal to noise ratios (SNRs). Anticipating such observations, this paper shows how constraints from a large number of events with modest SNRs can be combined to produce strong constraints on deviations from GR. Using kludge modified GR waveforms, we demonstrate how this test could identify certain types of deviations from GR if such deviations are present in the signal waveforms. We also study the robustness of this test against reasonable variations of a variety of different analysis parameters.