We study the mapping from Lagrangian to Eulerian space in the context of the Effective Field Theory (EFT) of Large Scale Structure. We compute Lagrangian displacements with Lagrangian Perturbation Theory (LPT) and perform the full non-perturbative transformation from displacement to density. When expanded up to a given order, this transformation reproduces the standard Eulerian Perturbation Theory (SPT) at the same order. However, the full transformation from displacement to density also includes higher order terms. These terms explicitly resum long wavelength motions, thus making the resulting density field better correlated with the true non-linear density field. As a result, the regime of validity of this approach is expected to extend that of the Eulerian EFT, and match that of the IR-resummed Eulerian EFT. This approach thus effectively enables a test of the IR-resummed EFT at the field level. We estimate the size of stochastic, non-perturbative contributions to the matter density power spectrum. We find that in our highest order calculation, at redshift z = 0 the power spectrum of the density field is reproduced with an accuracy of 1% (10%) up to k = 0.25 hMpc-1 (k = 0.46 hMpc-1). We believe that the dominant source of the remaining error is the stochastic contribution. Unfortunately, on these scales the stochastic term does not yet scale as k4 as it does in the very low k regime. Thus, modeling this contribution might be challenging.

We study Lagrangian Perturbation Theory (LPT) and its regularization in the Effective Field Theory (EFT) approach. We evaluate the LPT displacement with the same phases as a corresponding N-body simulation, which allows us to compare perturbation theory to the non-linear simulation with significantly reduced cosmic variance, and provides a more stringent test than simply comparing power spectra. We reliably detect a non-vanishing leading order EFT coefficient and a stochastic displacement term, uncorrelated with the LPT terms. This stochastic term is expected in the EFT framework, and, to the best of our understanding, is not an artifact of numerical errors or transients in our simulations. This term constitutes a limit to the accuracy of perturbative descriptions of the displacement field and its phases, corresponding to a 1% error on the non-linear power spectrum at k = 0.2 h-1Mpc at z = 0. Predicting the displacement power spectrum to higher accuracy or larger wavenumbers thus requires a model for the stochastic displacement.

Searches for gravitational waves crucially depend on exact signal processing of noisy strain data from gravitational wave detectors, which are known to exhibit significant nonstationary and non-Gaussian behavior. In this paper, we study two distinct effects in the LIGO/Virgo data that reduce the sensitivity of searches: first, variations in the noise power spectral density (PSD) on timescales of more than a few seconds; and second, loud and abrupt transient "glitches"of terrestrial or instrumental origin. We derive a simple procedure to correct, at first order, the effect of the variation in the PSD on the search background. Given the knowledge of the existence of localized glitches, in particular segments of data, we also develop a method to insulate statistical inference from these glitches, so as to cleanly excise them without affecting the search background in neighboring seconds. We show the importance of applying these methods on the publicly available LIGO data and estimate an increase in the detection volume of at least 15% from the PSD-drift correction alone, due to the improved background distribution.

We introduce a new technique to search for gravitational wave events from compact binary mergers that produce a clear signal only in a single gravitational wave detector, and marginal signals in other detectors. Such a situation can arise when the detectors in a network have different sensitivities, or when sources have unfavorable sky locations or orientations. We start with a short list of loud single-detector triggers from regions of parameter space that are empirically unaffected by glitches (after applying signal-quality vetoes). For each of these triggers, we compute evidence for astrophysical origin from the rest of the detector network by coherently combining the likelihoods from all detectors and marginalizing over extrinsic geometric parameters. We report the discovery of two new binary black hole (BBH) mergers in the second observing run of Advanced LIGO and Virgo (O2), in addition to the ones that were reported in [B. P. Abbott (LIGO Scientific and Virgo Collaborations), Phys. Rev. X 9, 031040 (2019)PRXHAE2160-330810.1103/PhysRevX.9.031040 and [T. Venumadhav, Phys. Rev. D 101, 083030 (2020)PRVDAQ2470-001010.1103/PhysRevD.101.083030]. We estimate that the two events have false alarm rates of one in 19 years (60 O2) and one in 11 years (36 O2). One of the events, GW170817A, has primary and secondary masses m1src=56-10+16 M⊙ and m2src=40-11+10 M⊙ in the source frame. The existence of GW170817A should be very informative about the theoretically predicted upper mass gap for stellar mass black holes. Its effective spin parameter is measured to be χeff=0.5±0.2, which is consistent with the tendency of the heavier detected BBH systems to have large and positive effective spin parameters. The other event, GWC170402, will be discussed thoroughly in future work.

The intergalactic medium is expected to be at its coldest point before the formation of the first stars in the Universe. Motivated by recent results from the EDGES experiment, we revisit the standard calculation of the kinetic temperature of the neutral gas through this period. When the first ultraviolet (UV) sources turn on, photons redshift into the Lyman lines of neutral hydrogen and repeatedly scatter within the Lyman-α line. They heat the gas via atomic recoils, and, through the Wouthuysen-Field effect, set the spin temperature of the 21-cm hyperfine (spin-flip) line of atomic hydrogen in competition with the resonant cosmic microwave background (CMB) photons. We show that the Lyman-α photons also mediate energy transfer between the CMB photons and the thermal motions of the hydrogen atoms. In the absence of x-ray heating, this new mechanism is the major correction to the temperature of the adiabatically cooling gas (∼10% at z=17), and is several times the size of the heating rate found in previous calculations. We also find that the effect is more dramatic in nonstandard scenarios that either enhance the radio background above the CMB or invoke new physics to cool the gas in order to explain the EDGES results. The coupling with the radio background can reduce the depth of the 21-cm absorption feature by almost a factor of 2 relative to the case with no sources of heating, and prevent the feature from developing a flattened bottom. As an inevitable consequence of the UV background that generates the absorption feature, this heating should be accounted for in any theoretical prediction.

We present a reanalysis of GW151226, the second binary black hole merger discovered by the LIGO-Virgo Collaboration. Previous analysis showed that the best-fit waveform for this event corresponded to the merger of a ∼14 M⊙ black hole with a ∼7.5 M companion, and the posterior distribution in mass ratio (q≤1) is rather flat. In this work, we perform parameter estimation using a waveform model that includes the effects of orbital precession and higher-order radiative multipole modes, and we find that the source parameters of GW151226 shift toward the low q and high effective spin (χeff) region and that q is better measured. The new solution has a log likelihood roughly two points higher than when either higher multipoles or orbital precession is neglected and can alter the astrophysical interpretation of GW151226. Additionally, we find it useful to use a flat-in-χeff prior, which does not penalize the large |χeff| region, in order to uncover the higher likelihood region for GW151226. Our solution has several interesting properties: (a) the secondary black hole mass is close to the upper limit of the hypothesized lower mass gap of astrophysical black hole population; and (b) orbital precession is driven by the primary black hole spin, which has a dimensionless magnitude as large as ∼0.85 and is tilted away from the orbital angular momentum at an angle of ∼57°. Since GW151226 is a relatively weak signal, an unambiguous claim of the detection of these effects in the signal cannot be made.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.