We find that the cosmic microwave background temperature and polarization power spectra measurements from Planck constrain the parameter $\Omega_bh^2$ mostly through: A) the amplitude of Thomson scattering and B) a factor that ensures Thomson scattering does not violate momentum conservation of the baryon-photon fluid. This allows us to obtain two distinct but comparably strong constraints on $\Omega_b h^2$ from the Planck data alone. They are consistent, showing robustness of the Planck $\Omega_b h^2$ constraint. We can alternatively rephrase these constraints as A) the change of the Thomson scattering cross section since recombination is less than $\sim 2\%$ and B) momentum during recombination is conserved to better than $\sim 2\%$ by Thomson scattering. Decoupling the eight various ways in which $\Omega_b h^2$ affects the Planck data leads to $H_0$ only slightly higher than in the standard analysis, $(68.3 \pm 1.6)\,\mathrm{km/s/Mpc}$. The overall consistency of all $\Omega_b h^2$ constraints does not suggest any problem with the standard cosmological model.

We apply a recently proposed cross-correlation power spectrum technique to study relationship between the ultra-high energy cosmic ray flux from the Pierre Auger Observatory and galaxies from the 2MASS Redshift Survey. Using a simple linear bias model relative to the galaxy auto power spectrum we are able to constrain the value of bias to be less than 5%, with the best fit value of about 1% detected with statistical significance of $\sim 1.8\sigma$. We study energy dependence of the bias, but the small sample size does not allow us to reach any statistically significant conclusions. For the cosmic ray events above 52 EeV we discover a curious excess cross-correlation at $\sim 1^\circ$ degree scales. Given similar cross-correlation is not visible at larger angular scales, statistical fluctuation seems like the most plausible explanation.

We apply a tension metric $Q_\textrm{UDM}$, the update difference in mean parameters, to understand the source of the difference in the measured Hubble constant $H_0$ inferred with cosmic microwave background lensing measurements from the Planck satellite ($H_0=67.9^{+1.1}_{-1.3}\, \mathrm{km/s/Mpc}$) and from the South Pole Telescope ($H_0=72.0^{+2.1}_{-2.5}\, \mathrm{km/s/Mpc}$) when both are combined with baryon acoustic oscillation (BAO) measurements with priors on the baryon density (BBN). $Q_\textrm{UDM}$ isolates the relevant parameter directions for tension or concordance where the two data sets are both informative, and aids in the identification of subsets of data that source the observed tension. With $Q_\textrm{UDM}$, we uncover that the difference in $H_0$ is driven by the tension between Planck lensing and BAO+BBN, at probability-to-exceed of 6.6%. Most of this mild tension comes from the galaxy BAO measurements parallel to the line of sight. The redshift dependence of the parallel BAOs pulls both the matter density $\Omega_m$ and $H_0$ high in $\Lambda$CDM, but these parameter anomalies are usually hidden when the BAO measurements are combined with other cosmological data sets with much stronger $\Omega_m$ constraints.

We analyze the final release of the Planck satellite data to constrain the gravitational lensing potential in a model-independent manner. The amount of lensing determined from the smoothing of the acoustic peaks in the temperature and polarization power spectra is 2$\sigma$ too high when compared with the measurements using the lensing reconstruction and 2.8$\sigma$ too high when compared with $\Lambda$CDM expectation based on the "unlensed" portion of the temperature and polarization power spectra. The largest change from the previous data release is the $\Lambda$CDM expectation, driven by improved constraints to the optical depth to reionization. The anomaly still is inconsistent with actual gravitational lensing, given that the lensing reconstruction constraints are discrepant independent of the model. Within the context of $\Lambda$CDM, improvements in its parameter constraints from lensing reconstruction bring this tension to 2.1$\sigma$ and from further adding baryon acoustic oscillation and Pantheon supernova data to a marginally higher 2.2$\sigma$. Once these other measurements are included, marginalizing this lensing-like anomaly cannot substantially resolve tensions with low-redshift measurements of $H_0$ and $S_8$ in $\Lambda$CDM, $\Lambda$CDM+$N_\mathrm{eff}$ or $\Lambda$CDM+$\sum m_\nu$; furthermore the artificial strengthening of constraints on $\sum m_\nu$ is less than 20%.

Galaxy spins can be predicted from the initial conditions in the early Universe through the tidal tensor twist. In simulations, their directions are well preserved through cosmic time, consistent with expectations of angular momentum conservation. We report a $\sim 3 \sigma$ detection of correlation between observed oriented directions of galaxy angular momenta and their predictions based on the initial density field reconstructed from the positions of SDSS galaxies. The detection is driven by a group of spiral galaxies classified by the Galaxy Zoo as (anti-)clockwise, with a modest improvement from adding galaxies from MaNGA and SAMI surveys. This is the first such detection of the oriented galaxy spin direction, which opens a way to use measurements of galaxy spins to probe fundamental physics in the early Universe.

Chiral symmetry is maximally violated in weak interactions, and such microscopic asymmetries in the early Universe might leave observable imprints on astrophysical scales without violating the cosmological principle. In this Letter, we propose a helicity measurement to detect primordial chiral violation. We point out that observations of halo-galaxy angular momentum directions (spins), which are frozen in during the galaxy formation process, provide a fossil chiral observable. From the clustering mode of large scale structure of the Universe, we construct a spin mode in Lagrangian space and show in simulations that it is a good probe of halo-galaxy spins. In standard model, a strong symmetric correlation between the left and right helical components of this spin mode and galaxy spins is expected. Measurements of these correlations will be sensitive to chiral breaking, providing a direct test of chiral symmetry breaking in the early Universe.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

We report on the measurement of coherent radio emission from the electron beam sudden appearance at the Telescope Array Electron Light Source facility. This emission was detected by four independent radio detector setups sensitive to frequencies ranging from 50 MHz up to 12.5 GHz. We show that this phenomenon can be understood as a special case of coherent transition radiation by comparing the observed results with simulations. The in-nature application of this signal is given by the emission of cosmic ray or neutrino induced particle cascades traversing different media such as air, rock and ice.

We investigate correlations induced by gravitational lensing on simulated cosmic microwave background data of experiments with an incomplete sky coverage and their effect on inferences from the South Pole Telescope data. These correlations agree well with the theoretical expectations, given by the sum of super-sample and intra-sample lensing terms, with only a typically negligible $\sim$ 5% discrepancy in the amplitude of the super-sample lensing effect. Including these effects we find that lensing constraints are in $3.0\sigma$ or $2.1\sigma$ tension between the SPT polarization measurements and Planck temperature or lensing reconstruction constraints respectively. If the lensing-induced covariance effects are neglected, the significance of these tensions increases to $3.5\sigma$ or $2.5\sigma$. Using the standard scaling parameter $A_L$ substantially underestimates the significance of the tension once other parameters are marginalized over. By parameterizing the super-sample lensing through the mean convergence in the SPT footprint, we find a hint of underdensity in the SPT region. We also constrain extra sharpening of the CMB acoustic peaks due to missing smoothing of the peaks by super-sample lenses at a level that is much smaller than the lens sample variance. Finally, we extend the usual "shift in the means" statistic for evaluating tensions to non-Gaussian posteriors, generalize an approach to extract correlation modes from noisy simulated covariance matrices, and present a treatment of correlation modes not as data covariances but as auxiliary model parameters.

We apply a recently developed method to directly measure the gravitational lensing power spectrum from CMB power spectra to the Planck satellite data. This method allows us to analyze the tension between the temperature power spectrum and lens reconstruction in a model independent way. Even when allowing for arbitrary variations in the lensing power spectrum, the tension remains at the 2.4$\sigma$ level. By separating the lensing and unlensed high redshift information in the CMB power spectra, we also show that under $\Lambda$CDM the two are in tension at a similar level whereas the unlensed information is consistent with lensing reconstruction. These anomalies are driven by the smoother acoustic peaks relative to $\Lambda$CDM at $\ell \sim 1250 - 1500$. Both tensions relax slightly when polarization data are considered. This technique also isolates the one aspect of the lensing power spectrum that the Planck CMB power spectra currently constrain and can be straightforwardly generalized to future data when CMB power spectra constrain multiple aspects of lensing which are themselves correlated with lensing reconstruction.

Non-Gaussian correlations induced in CMB power spectra by gravitational lensing must be included in likelihood analyses for future CMB experiments. We present a simple but accurate likelihood model which includes these correlations and use it for Markov Chain Monte Carlo parameter estimation from simulated lensed CMB maps in the context of $\Lambda$CDM and extensions which include the sum of neutrino masses or the dark energy equation of state $w$. If lensing-induced covariance is not taken into account for a CMB-S4 type experiment, the errors for one combination of parameters in each case would be underestimated by more then a factor of two and lower limits on $w$ could be misestimated substantially. The frequency of falsely ruling out the true model or finding tension with other data sets would also substantially increase. Our analysis also enables a separation of lens and unlensed information from CMB power spectra, which provides for consistency tests of the model and, if combined with other such measurements, a nearly lens-sample-variance free test for systematics and new physics in the unlensed spectrum. This parameterization also leads to a simple effective likelihood that can be used to assist model building in case consistency tests of $\Lambda$CDM fail.

Gravitational lensing information from the two and higher point statistics of the CMB temperature and polarization fields are intrinsically correlated because they are lensed by the same realization of structure between last scattering and observation. Using an analytic model for lens sample covariance, we show that there is one mode, separately measurable in the lensed CMB power spectra and lensing reconstruction, that carries most of this correlation. Once these measurements become lens sample variance dominated, this mode should provide a useful consistency check between the observables that is largely free of sampling and cosmological parameter errors. Violations of consistency could indicate systematic errors in the data and lens reconstruction or new physics at last scattering, any of which could bias cosmological inferences and delensing for gravitational waves. A second mode provides a weaker consistency check for a spatially flat universe. Our analysis isolates the additional information supplied by lensing in a model independent manner but is also useful for understanding and forecasting CMB cosmological parameter errors in the extended $\Lambda$CDM parameter space of dark energy, curvature and massive neutrinos. We introduce and test a simple but accurate forecasting technique for this purpose that neither double counts lensing information nor neglects lensing in the observables.

We investigate whether transition radiation from a particle shower crossing the interface between Earth and air and induced by an Earth-skimming neutrino can explain the upward event announced recently by the ANITA Collaboration. While the properties of the observed signal can in principle be explained with transition radiation, the flux necessary for a successful explanation is in tension with the current best limits from the Pierre Auger Observatory, the IceCube neutrino detector and the ANITA balloon. We also show that in this scenario, the direction of the incoming neutrino is determined precisely to within a few degrees, from the polarization properties of the observed events and thanks to the Earth opacity to ultra high energy neutrinos.

Coherent radiation at radio frequencies from high-energy showers fully contained in a dense radio-transparent medium - like ice, salt or regolith - has been extensively investigated as a promising technique to search for ultra-high energy (UHE) neutrinos. Additional emission in the form of transition radiation may occur when a neutrino-induced shower produced close to the Earth surface emerges from the ground into atmospheric air. We present the first detailed evaluation of transition radiation from high-energy showers crossing the boundary between two different media. We found that transition radiation is sizable over a wide solid angle and coherent up to $\sim$ 1 GHz. These properties encourage further work to evaluate the potential of a large-aperture UHE neutrino experiment based on detection of transition radiation.

Self-accelerating backgrounds in massive gravity provide an arena to explore the Cauchy problem for derivatively coupled fields that obey complex constraints which reduce the phase space degrees of freedom. We present here an algorithm based on the Kronecker form of a matrix pencil that finds all hidden constraints, for example those associated with derivatives of the equations of motion, and characteristic curves for any 1+1 dimensional system of linear partial differential equations. With the Regge-Wheeler-Zerilli decomposition of metric perturbations into angular momentum and parity states, this technique applies to fully 3+1 dimensional perturbations of massive gravity around any isotropic self-accelerating background. Five spin modes of the massive graviton propagate once the constraints are imposed: two spin-2 modes with luminal characteristics present in the massless theory as well as two spin-1 modes and one spin-0 mode. Although the new modes all possess the same - typically spacelike - characteristic curves, the spin-1 modes are parabolic while the spin-0 modes are hyperbolic. The joint system, which remains coupled by non-derivative terms, cannot be solved as a simple Cauchy problem from a single non-characteristic surface. We also illustrate the generality of the algorithm with other cases where derivative constraints reduce the number of propagating degrees of freedom or order of the equations.

We present a concept for large-area, low-cost detection of ultra-high energy cosmic rays (UHECRs) with a Fluorescence detector Array of Single-pixel Telescopes (FAST), addressing the requirements for the next generation of UHECR experiments. In the FAST design, a large field of view is covered by a few pixels at the focal plane of a mirror or Fresnel lens. We report first results of a FAST prototype installed at the Telescope Array site, consisting of a single 200 mm photomultiplier tube at the focal plane of a 1 m$^2$ Fresnel lens system taken from the prototype of the JEM-EUSO experiment. The FAST prototype took data for 19 nights, demonstrating remarkable operational stability. We detected laser shots at distances of several kilometres as well as 16 highly significant UHECR shower candidates.

Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the general interplay between superluminality, the well-posedness of the Cauchy problem, and strong coupling. For three particular classes of vacuum solutions, one of which is new to this work, we construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy perturbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a pathological frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike characteristics do and do not originate from a point where perturbation theory breaks down and where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given observer. The global structure of spacetime also reveals issues that are unique to theories with two metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singularities. Characteristics run tangent to these surfaces requiring {\it ad hoc} rules to establish continuity across singularities.

We undertake a complete and covariant treatment for the quadratic Lagrangian of all of the degrees of freedom of massive gravity with a fixed flat fiducial metric for arbitrary massive gravity parameters around any isotropic self-accelerating background solution. Generically, 3 out of 4 Stuckelberg degrees of freedom propagate in addition to the usual 2 tensor degrees of freedom of general relativity. The complete kinetic structure typically is only revealed at an order in the graviton mass that is equivalently to retaining curvature terms in a locally flat expansion. These results resolve several apparent discrepancies in the literature where zero degrees of freedom propagate in either special cases or approximate treatments as well as decoupling limit analyses which attempt to count longitudinal degrees of freedom.

In a recent Letter (Rimmer et al. 2014), Jupiter is presented as an efficient detector for Ultra-High Energy Cosmic Rays (UHECRs), through measurement by an Earth-orbiting satellite of gamma rays from UHECRs showers produced in Jupiter's atmosphere. We show that this result is incorrect, due to erroneous assumptions on the angular distribution of shower particles. We evaluated other Solar System objects as potential targets for UHECRs detection, and found that the proposed technique is either not viable or not competitive with traditional ground-based UHECRs detectors.

Radio emission from Ultra-High Energy Cosmic Rays (UHECR) showers detected after specular reflection off the Antarctic ice surface has been recently demonstrated by the ANITA balloon-borne experiment. An antenna observing a large area of ice or water from a mountaintop, a balloon or a satellite may be competitive with more conventional techniques. We present an estimate of the exposure of a high altitude antenna, which provides insight on the prospects of this technique for UHECR detection. We find that a satellite antenna may reach a significantly larger exposure than existing UHECR observatories, but an experimental characterization of the radio reflected signal is required to establish the potential of this approach. A balloon-borne or a mountaintop antenna are found not to be competitive under any circumstances.