Tau neutrinos are expected to comprise roughly one third of both the astrophysical and cosmogenic neutrino flux, but currently the flavor ratio is poorly constrained and the expected flux at energies above $10^{17}$ eV is low. We present a detector concept aimed at measuring the diffuse flux of tau neutrinos in this energy range via a high-elevation mountaintop detector using the radio technique. The detector searches for radio signals from upgoing air showers generated by Earth-skimming tau neutrinos. Signals from several antennas in a compact array are coherently summed at the trigger level, permitting not only directional masking of anthropogenic backgrounds, but also a low trigger threshold. This design takes advantage of both the large viewing area available at high-elevation sites and the nearly full duty cycle available to radio instruments. We present trade studies that consider the station elevation, frequency band, number of antennas in the array, and the trigger threshold to develop a highly efficient station design. Such a mountaintop detector can achieve a factor of ten improvement in acceptance over existing instruments with 100 independent stations. With 1000 stations and three years of observation, it can achieve a sensitivity to an integrated $\mathcal{E}^{-2}$ flux of $<10^{-9}$ GeV cm$^{-2}$ sr$^{-1}$ s$^{-1}$, in the range of the expected flux of all-flavor cosmogenic neutrinos assuming a pure iron cosmic-ray composition.

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale Cosmic Microwave Background (CMB) polarization from the South Pole in search of a primordial $B$-mode signature. This $B$-mode signal arises from primordial gravitational waves interacting with the CMB, and has amplitude parametrized by the tensor-to-scalar ratio $r$. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier Transform Spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on $r$ induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.

Branches of cosmic inflationary models, such as slow-roll inflation, predict a background of primordial gravitational waves that imprints a unique odd-parity B-mode pattern in the Cosmic Microwave Background (CMB) at amplitudes that are within experimental reach. The BICEP/Keck (BK) experiment targets this primordial signature, the amplitude of which is parameterized by the tensor-to-scalar ratio r, by observing the polarized microwave sky through the exceptionally clean and stable atmosphere at the South Pole. B-mode measurements require an instrument with exquisite sensitivity, tight control of systematics, and wide frequency coverage to disentangle the primordial signal from the Galactic foregrounds. BICEP Array represents the most recent stage of the BK program, and comprises four BICEP3-class receivers observing at 30/40, 95, 150 and 220/270 GHz. The 30/40 GHz receiver will be deployed at the South Pole during the 2019/2020 austral summer. After 3 full years of observations with 30,000+ detectors, BICEP Array will measure primordial gravitational waves to a precision $\sigma(r)$ between 0.002 and 0.004, depending on foreground complexity and the degree of lensing removal. In this paper we give an overview of the instrument, highlighting the design features in terms of cryogenics, magnetic shielding, detectors and readout architecture as well as reporting on the integration and tests that are ongoing with the first receiver at 30/40 GHz.

The BICEP/Keck (BK) experiment aims to detect the imprint of primordial gravitational waves in the Cosmic Microwave Background polarization, which would be direct evidence of the inflation theory. While the tensor-to-scalar ratio has been constrained to be r_0.05 < 0.06 at 95% c.l., further improvements on this upper limit are hindered by polarized Galactic foreground emissions and removal of gravitational lensing polarization. The 30/40 GHz receiver of the BICEP Array (BA) will deploy at the end of 2019 and will constrain the synchrotron foreground with unprecedented accuracy within the BK sky patch. We will show the design of the 30/40 GHz detectors and test results summarizing its performance. The low optical and atmospheric loading at these frequencies requires our TES detectors to have low saturation power in order to be photon-noise dominated. To realize the low thermal conductivity required from a 250 mK base temperature, we developed new bolometer leg designs. We will present the relevant measured detector parameters: G, Tc, Rn, Psat , and spectral bands, and noise spectra. We achieved a per bolometer NEP including all noise components of 2.07E-17 W/sqrt(Hz), including an anticipated photon noise level 1.54E-17 W/sqrt(Hz).

NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions.

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

The past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

RNO is the mid-scale discovery instrument designed to make the first observation of neutrinos from the cosmos at extreme energies, with sensitivity well beyond current instrument capabilities. This new observatory will be the largest ground-based neutrino telescope to date, enabling the measurement of neutrinos above $10^{16}$ eV, determining the nature of the astrophysical neutrino flux that has been measured by IceCube at higher energies, similarly extending the reach of multi-messenger astrophysics to the highest energies, and enabling investigations of fundamental physics at energies unreachable by particle accelerators on Earth.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

A Monte Carlo simulation program for the radio detection of Ultra High Energy (UHE) neutrino interactions in the Antarctic ice as viewed by the Antarctic Impulsive Transient Antenna (ANITA) is described in this article. The program, icemc, provides an input spectrum of UHE neutrinos, the parametrization of the Askaryan radiation generated by their interaction in the ice, and the propagation of the radiation through ice and air to a simulated model of the third and fourth ANITA flights. This paper provides an overview of the icemc simulation, descriptions of the physics models used and of the ANITA electronics processing chain, data/simulation comparisons to validate the predicted performance, and a summary of the impact of published results.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

High-energy cosmic neutrinos carry unique information about the most energetic non-thermal sources in the Universe. This white paper describes the outstanding astrophysics questions that neutrino astronomy can address in the coming decade. A companion white paper discusses how the observation of cosmic neutrinos can address open questions in fundamental physics. Detailed measurements of the diffuse neutrino flux, measurements of neutrinos from point sources, and multi-messenger observations with neutrinos will enable the discovery and characterization of the most energetic sources in the Universe.

High-energy cosmic neutrinos can reveal new fundamental particles and interactions, probing energy and distance scales far exceeding those accessible in the laboratory. This white paper describes the outstanding particle physics questions that high-energy cosmic neutrinos can address in the coming decade. A companion white paper discusses how the observation of cosmic neutrinos can address open questions in astrophysics. Tests of fundamental physics using high-energy cosmic neutrinos will be enabled by detailed measurements of their energy spectrum, arrival directions, flavor composition, and timing.

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

The ANITA experiment, which is designed to detect ultra-high energy neutrinos, has reported the observation of two anomalous events, directed at angles of $27^{\circ}$ and $35^{\circ}$ with respect to the horizontal. At these angles, the Earth is expected to efficiently absorb ultra-high energy neutrinos, making the origin of these events unclear and motivating explanations involving physics beyond the Standard Model. In this study, we consider the possibility that ANITA's anomalous events are the result of Askaryan emission produced by exotic weakly interacting particles scattering elastically with nuclei in the Antarctic ice sheet. Such particles could be produced by superheavy ($\sim 10^{10}-10^{13}$ GeV) dark matter particles decaying in the halo of the Milky Way. Such scenarios can be constrained by existing measurements of the high-latitude gamma-ray background and the ultra-high energy cosmic ray spectrum, along with searches for ultra-high energy neutrinos by IceCube and other neutrino telescopes.

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.

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.

The ANtarctic Impulsive Transient Antenna (ANITA) NASA long-duration balloon payload completed its fourth flight in December 2016, after 28 days of flight time. ANITA is sensitive to impulsive broadband radio emission from interactions of ultra-high-energy neutrinos in polar ice (Askaryan emission). We present the results of two separate blind analyses searching for signals from Askaryan emission in the data from the fourth flight of ANITA. The more sensitive analysis, with a better expected limit, has a background estimate of $0.64^{+0.69}_{-0.45}$ and an analysis efficiency of $82\pm2\%$. The second analysis has a background estimate of $0.34^{+0.66}_{-0.16}$ and an analysis efficiency of $71\pm6\%$. Each analysis found one event in the signal region, consistent with the background estimate for each analysis. The resulting limit further tightens the constraints on the diffuse flux of ultra-high-energy neutrinos at energies above $10^{19.5}$ eV.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

The excess of neutrino candidate events detected by IceCube from the direction of TXS 0506+056 has generated a great deal of interest in blazars as sources of high-energy neutrinos. In this study, we analyze the publicly available portion of the IceCube dataset, performing searches for neutrino point sources in spatial coincidence with the blazars and other active galactic nuclei contained in the Fermi 3LAC and the Roma BZCAT catalogs, as well as in spatial and temporal coincidence with flaring sources identified in the Fermi Collaboration's All-Sky Variability Analysis (FAVA). We find no evidence that blazars generate a significant flux of high-energy neutrinos, and conclude that no more than 5-15% of the diffuse flux measured by IceCube can originate from this class of objects. While we cannot rule out the possibility that TXS 0506+056 has at times generated significant neutrino emission, we find that such behavior cannot be common among blazars, requiring TXS 0506+056 to be a rather extreme outlier and not representative of the overall blazar population. The bulk of the diffuse high-energy neutrino flux must instead be generated by a significantly larger population of less-luminous sources, such as non-blazar active galactic nuclei.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

The compelling science case for the observation of B-mode polarization in the cosmic microwave background (CMB) is driving the CMB community to expand the observed sky fraction, either by extending survey sizes or by deploying receivers to potential new northern sites. For ground-based CMB instruments, poorly-mixed atmospheric water vapor constitutes the primary source of short-term sky noise. This results in short-timescale brightness fluctuations, which must be rejected by some form of modulation. To maximize the sensitivity of ground-based CMB observations, it is useful to understand the effects of atmospheric water vapor over timescales and angular scales relevant for CMB polarization measurements. To this end, we have undertaken a campaign to perform a coordinated characterization of current and potential future observing sites using scanning 183 GHz water vapor radiometers (WVRs). So far, we have deployed two identical WVR units; one at the South Pole, Antarctica, and the other at Summit Station, Greenland. The former site has a long heritage of ground-based CMB observations and is the current location of the Bicep/Keck Array telescopes as well as the South Pole Telescope. The latter site, though less well characterized, is under consideration as a northern-hemisphere location for future CMB receivers. Data collection from this campaign began in January 2016 at South Pole and July 2016 at Summit Station. Data analysis is ongoing to reduce the data to a single spatial and temporal statistic that can be used for one-to-one site comparison.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

We present measurements of radio transmission in the $\sim$100 MHz range through a $\sim100$ m deep region below the surface of the ice at Summit Station, Greenland, called the firn. In the firn, the index of refraction changes due to the transition from snow at the surface to glacial ice below, affecting the propagation of radio signals in that region. We compare our observations to a finite-difference time-domain (FDTD) electromagnetic wave simulation, which supports the existence of three classes of propagation: a bulk propagation ray-bending mode that leads to so-called "shadowed" regions for certain geometries of transmission, a surface-wave mode induced by the ice/air interface, and an arbitrary-depth horizontal propagation mode that requires perturbations from a smooth density gradient. In the non-shadowed region, our measurements are consistent with the bulk propagation ray-bending mode both in timing and in amplitude. We also observe signals in the shadowed region, in conflict with a bulk-propagation-only ray-bending model, but consistent with FDTD simulations using a variety of firn models for Summit Station. The amplitude and timing of our measurements in all geometries are consistent with the predictions from FDTD simulations. In the shadowed region, the amplitude of the observed signals is consistent with a best-fit coupling fraction value of $2.4$% (0.06% in power) or less to a surface or horizontal propagation mode from the bulk propagation mode. The relative amplitude of observable signals in the two regions is important for experiments that aim to detect radio emission from astrophysical high-energy neutrinos interacting in glacial ice, which rely on a radio propagation model to inform simulations and perform event reconstruction.

The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky.

The Antarctic Impulsive Transient Antenna (ANITA), a NASA long-duration balloon payload, searches for radio emission from interactions of ultra-high-energy neutrinos in polar ice. The third flight of ANITA (ANITA-III) was launched in December 2014 and completed a 22-day flight. We present the results of three analyses searching for Askaryan radio emission of neutrino origin. In the most sensitive of the analyses, we find one event in the signal region on an expected a priori background of $0.7^{+0.5}_{-0.3}$. Though consistent with the background estimate, the candidate event remains compatible with a neutrino hypothesis even after additional post-unblinding scrutiny.

We report on an upward traveling, radio-detected cosmic-ray-like impulsive event with characteristics closely matching an extensive air shower. This event, observed in the third flight of the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload, is consistent with a similar event reported in a previous flight. These events may be produced by the atmospheric decay of an upward-propagating $\tau$-lepton produced by a $\nu_{\tau}$ interaction, although their relatively steep arrival angles create tension with the standard model (SM) neutrino cross section. Each of the two events have $a~posteriori$ background estimates of $\lesssim 10^{-2}$ events. If these are generated by $\tau$-lepton decay, then either the charged-current $\nu_{\tau}$ cross section is suppressed at EeV energies, or the events arise at moments when the peak flux of a transient neutrino source was much larger than the typical expected cosmogenic background neutrinos.

The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.

The Antarctic Impulsive Transient Antenna (ANITA) is a NASA long-duration balloon experiment with the primary goal of detecting ultra-high-energy ($>10^{18}\,\mbox{eV}$) neutrinos via the Askaryan Effect. The fourth ANITA mission, ANITA-IV, recently flew from Dec 2 to Dec 29, 2016. For the first time, the Tunable Universal Filter Frontend (TUFF) boards were deployed for mitigation of narrow-band, anthropogenic noise with tunable, switchable notch filters. The TUFF boards also performed second-stage amplification by approximately 45 dB to boost the $\sim\,\mu\mbox{V-level}$ radio frequency (RF) signals to $\sim$ mV-level for digitization, and supplied power via bias tees to the first-stage, antenna-mounted amplifiers. The other major change in signal processing in ANITA-IV is the resurrection of the $90^{\circ}$ hybrids deployed previously in ANITA-I, in the trigger system, although in this paper we focus on the TUFF boards. During the ANITA-IV mission, the TUFF boards were successfully operated throughout the flight. They contributed to a factor of 2.8 higher total instrument livetime on average in ANITA-IV compared to ANITA-III due to reduction of narrow-band, anthropogenic noise before a trigger decision is made.

The in-ice radio interferometric phased array technique for detection of high energy neutrinos looks for Askaryan emission from neutrinos interacting in large volumes of glacial ice, and is being developed as a way to achieve a low energy threshold and a large effective volume at high energies. The technique is based on coherently summing the impulsive Askaryan signal from multiple antennas, which increases the signal-to-noise ratio for weak signals. We report here on measurements and a simulation of thermal noise correlations between nearby antennas, beamforming of impulsive signals, and a measurement of the expected improvement in trigger efficiency through the phased array technique. We also discuss the noise environment observed with an analog phased array at Summit Station, Greenland, a possible site for an interferometric phased array for radio detection of high energy neutrinos.

CMB-S4 is a proposed experiment to map the polarization of the Cosmic Microwave Background (CMB) to nearly the cosmic variance limit for angular scales that are accessible from the ground. The science goals and capabilities of CMB-S4 in illuminating cosmic inflation, measuring the sum of neutrino masses, searching for relativistic relics in the early universe, characterizing dark energy and dark matter, and mapping the matter distribution in the universe have been described in the CMB-S4 Science Book. This Technology Book is a companion volume to the Science Book. The ambitious science goals of CMB-S4, a "Stage-4" experiment, require a step forward in experimental capability from the current Stage=II experiments. To guide this process, we summarize the current state of CMB instrumentation technology, and identify R&D efforts necessary to advance it for use in CMB-S4. The book focuses on technical challenges in four broad areas: Telescope Design; Receiver Optics; Focal-Plane Optical Coupling; and Focal-Plane Sensor and Readout.

The primary science goal of the NASA-sponsored ANITA project is measurement of ultra-high energy neutrinos and cosmic rays, observed via radio-frequency signals resulting from a neutrino- or cosmic ray- interaction with terrestrial matter (atmospheric or ice molecules, e.g.). Accurate inference of the energies of these cosmic rays requires understanding the transmission/reflection of radio wave signals across the ice-air boundary. Satellite-based measurements of Antarctic surface reflectivity, using a co-located transmitter and receiver, have been performed more-or-less continuously for the last few decades. Satellite-based reflectivity surveys, at frequencies ranging from 2--45 GHz and at near-normal incidence, yield generally consistent reflectivity maps across Antarctica. Using the Sun as an RF source, and the ANITA-3 balloon borne radio-frequency antenna array as the RF receiver, we have also measured the surface reflectivity over the interval 200-1000 MHz, at elevation angles of 12-30 degrees, finding agreement with the Fresnel equations within systematic errors. To probe low incidence angles, inaccessible to the Antarctic Solar technique and not probed by previous satellite surveys, a novel experimental approach ("HiCal-1") was devised. Unlike previous measurements, HiCal-ANITA constitute a bi-static transmitter-receiver pair separated by hundreds of kilometers. Data taken with HiCal, between 200--600 MHz shows a significant departure from the Fresnel equations, constant with frequency over that band, with the deficit increasing with obliquity of incidence, which we attribute to the combined effects of possible surface roughness, surface grain effects, radar clutter and/or shadowing of the reflection zone due to Earth curvature effects.

This is a review on radio detection of high energy neutrinos.

BICEP3, the latest telescope in the BICEP/Keck program, started science observations in March 2016. It is a 550mm aperture refractive telescope observing the polarization of the cosmic microwave background at 95 GHz. We show the focal plane design and detector performance, including spectral response, optical efficiency and preliminary sensitivity of the upgraded BICEP3. We demonstrate 9.72$\mu$K$\sqrt{\textrm{s}}$ noise performance of the BICEP3 receiver.

BICEP3 is a 520 mm aperture, compact two-lens refractor designed to observe the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal plane consists of modularized tiles of antenna-coupled transition edge sensors (TESs), similar to those used in BICEP2 and the Keck Array. The increased per-receiver optical throughput compared to BICEP2/Keck Array, due to both its faster f/1.7 optics and the larger aperture, more than doubles the combined mapping speed of the BICEP/Keck program. The BICEP3 receiver was recently upgraded to a full complement of 20 tiles of detectors (2560 TESs) and is now beginning its second year of observation (and first science season) at the South Pole. We report on its current performance and observing plans. Given its high per-receiver throughput while maintaining the advantages of a compact design, BICEP3-class receivers are ideally suited as building blocks for a 3rd-generation CMB experiment, consisting of multiple receivers spanning 35 GHz to 270 GHz with total detector count in the tens of thousands. We present plans for such an array, the new "BICEP Array" that will replace the Keck Array at the South Pole, including design optimization, frequency coverage, and deployment/observing strategies.

BICEP3 is a small-aperture refracting cosmic microwave background (CMB) telescope designed to make sensitive polarization maps in pursuit of a potential B-mode signal from inflationary gravitational waves. It is the latest in the BICEP/Keck Array series of CMB experiments at the South Pole, which has provided the most stringent constraints on inflation to date. For the 2016 observing season, BICEP3 was outfitted with a full suite of 2400 optically coupled detectors operating at 95 GHz. In these proceedings we report on the far field beam performance using calibration data taken during the 2015-2016 summer deployment season in situ with a thermal chopped source. We generate high-fidelity per-detector beam maps, show the array-averaged beam profile, and characterize the differential beam response between co-located, orthogonally polarized detectors which contributes to the leading instrumental systematic in pair differencing experiments. We find that the levels of differential pointing, beamwidth, and ellipticity are similar to or lower than those measured for BICEP2 and Keck Array. The magnitude and distribution of BICEP3's differential beam mismatch - and the level to which temperature-to-polarization leakage may be marginalized over or subtracted in analysis - will inform the design of next-generation CMB experiments with many thousands of detectors.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

We report on four radio-detected cosmic-ray (CR) or CR-like events observed with the Antarctic Impulsive Transient Antenna (ANITA), a NASA-sponsored long-duration balloon payload. Two of the four were previously identified as stratospheric CR air showers during the ANITA-I flight. A third stratospheric CR was detected during the ANITA-II flight. Here we report on characteristics these three unusual CR events, which develop nearly horizontally, 20-30~km above the surface of the Earth. In addition, we report on a fourth steeply upward-pointing ANITA-I CR-like radio event which has characteristics consistent with a primary that emerged from the surface of the ice. This suggests a possible $\tau$-lepton decay as the origin of this event, but such an interpretation would require significant suppression of the Standard Model ${\tau}$-neutrino cross section.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

We present results from an analysis of all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season. This includes the first Keck Array observations at 95 GHz. The maps reach a depth of 50 nK deg in Stokes $Q$ and $U$ in the 150 GHz band and 127 nK deg in the 95 GHz band. We take auto- and cross-spectra between these maps and publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353 GHz. An excess over lensed-LCDM is detected at modest significance in the 95x150 $BB$ spectrum, and is consistent with the dust contribution expected from our previous work. No significant evidence for synchrotron emission is found in spectra such as 23x95, or for correlation between the dust and synchrotron sky patterns in spectra such as 23x353. We take the likelihood of all the spectra for a multi-component model including lensed-LCDM, dust, synchrotron and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio $r$), using priors on the frequency spectral behaviors of dust and synchrotron emission from previous analyses of WMAP and Planck data in other regions of the sky. This analysis yields an upper limit $r_{0.05}<0.09$ at 95% confidence, which is robust to variations explored in analysis and priors. Combining these $B$-mode results with the (more model-dependent) constraints from Planck analysis of CMB temperature plus BAO and other data, yields a combined limit $r_{0.05}<0.07$ at 95% confidence. These are the strongest constraints to date on inflationary gravitational waves.

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. For 14 of these events, this radiation was reflected from the ice. The dominant contribution to the radiation from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. This radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of 36km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 EeV. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations. In addition, we find that the Monte Carlo simulation of the ANITA data set is in agreement with the total number of observed events and with the properties of those events.

The detection of high energy neutrinos ($10^{15}-10^{20}$ eV or $1-10^{5}$ PeV) is an important step toward understanding the most energetic cosmic accelerators and would enable tests of fundamental physics at energy scales that cannot easily be achieved on Earth. In this energy range, there are two expected populations of neutrinos: the astrophysical flux observed with IceCube at lower energies ($\sim1$ PeV) and the predicted cosmogenic flux at higher energies ($\sim10^{18}$ eV). Radio detector arrays such as RICE, ANITA, ARA, and ARIANNA exploit the Askaryan effect and the radio transparency of glacial ice, which together enable enormous volumes of ice to be monitored with sparse instrumentation. We describe here the design for a phased radio array that would lower the energy threshold of radio techniques to the PeV scale, allowing measurement of the astrophysical flux observed with IceCube over an extended energy range. Meaningful energy overlap with optical Cherenkov telescopes could be used for energy calibration. The phased radio array design would also provide more efficient coverage of the large effective volume required to discover cosmogenic neutrinos.

BICEP3 is a $550~mm$ aperture telescope with cold, on-axis, refractive optics designed to observe at the $95~GHz$ band from the South Pole. It is the newest member of the BICEP/Keck family of inflationary probes specifically designed to measure the polarization of the cosmic microwave background (CMB) at degree-angular scales. BICEP3 is designed to house 1280 dual-polarization pixels, which, when fully-populated, totals to $\sim$9$\times$ the number of pixels in a single Keck $95~GHz$ receiver, thus further advancing the BICEP/Keck program's $95~GHz$ mapping speed. BICEP3 was deployed during the austral summer of 2014-2015 with 9 detector tiles, to be increased to its full capacity of 20 in the second season. After instrument characterization measurements were taken, CMB observation commenced in April 2015. Together with multi-frequency observation data from Planck, BICEP2, and the Keck Array, BICEP3 is projected to set upper limits on the tensor-to-scalar ratio to $r$ $\lesssim 0.03$ at $95\%$ C.L..

In a companion paper we have reported a $>5\sigma$ detection of degree scale $B $-mode polarization at 150 GHz by the BICEP2 experiment. Here we provide a detailed study of potential instrumental systematic contamination to that measurement. We focus extensively on spurious polarization that can potentially arise from beam imperfections. We present a heuristic classification of beam imperfections according to their symmetries and uniformities, and discuss how resulting contamination adds or cancels in maps that combine observations made at multiple orientations of the telescope about its boresight axis. We introduce a technique, which we call "deprojection", for filtering the leading order beam-induced contamination from time ordered data, and show that it removes power from BICEP2's $BB$ spectrum consistent with predictions using high signal-to-noise beam shape measurements. We detail the simulation pipeline that we use to directly simulate instrumental systematics and the calibration data used as input to that pipeline. Finally, we present the constraints on $BB$ contamination from individual sources of potential systematics. We find that systematics contribute $BB$ power that is a factor $\sim10\times$ below BICEP2's 3-year statistical uncertainty, and negligible compared to the observed $BB$ signal. The contribution to the best-fit tensor/scalar ratio is at a level equivalent to $r=(3-6)\times10^{-3}$.

We report an in situ measurement of the electric field attenuation length at radio frequencies for the bulk ice at Summit Station, Greenland, made by broadcasting radio-frequency signals vertically through the ice and measuring the relative power in the return ground bounce signal. We find the depth-averaged field attenuation length to be 947 +92/-85 meters at 75 MHz. While this measurement has clear radioglaciological applications, the radio clarity of the ice also has implications for the detection of ultra-high energy (UHE) astrophysical particles via their radio emission in dielectric media such as ice. Assuming a reliable extrapolation to higher frequencies, the measured attenuation length at Summit Station is comparable to previously measured radio-frequency attenuation lengths at candidate particle detector sites around the world, and strengthens the case for Summit Station as a promising northern site for UHE neutrino detection.

The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.

BICEP2 and the Keck Array are polarization-sensitive microwave telescopes that observe the cosmic microwave background (CMB) from the South Pole at degree angular scales in search of a signature of inflation imprinted as B-mode polarization in the CMB. BICEP2 was deployed in late 2009, observed for three years until the end of 2012 at 150 GHz with 512 antenna-coupled transition edge sensor bolometers, and has reported a detection of B-mode polarization on degree angular scales. The Keck Array was first deployed in late 2010 and will observe through 2016 with five receivers at several frequencies (95, 150, and 220 GHz). BICEP2 and the Keck Array share a common optical design and employ the field-proven BICEP1 strategy of using small-aperture, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. Here we describe the optical design of both instruments and report a full characterization of the optical performance and beams of BICEP2 and the Keck Array at 150 GHz.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.

We have developed antenna-coupled transition-edge sensor (TES) bolometers for a wide range of cosmic microwave background (CMB) polarimetry experiments, including BICEP2, Keck Array, and the balloon borne SPIDER. These detectors have reached maturity and this paper reports on their design principles, overall performance, and key challenges associated with design and production. Our detector arrays repeatedly produce spectral bands with 20%-30% bandwidth at 95, 150, or 220~GHz. The integrated antenna arrays synthesize symmetric co-aligned beams with controlled side-lobe levels. Cross-polarized response on boresight is typically ~0.5%, consistent with cross-talk in our multiplexed readout system. End-to-end optical efficiencies in our cameras are routinely 35% or higher, with per detector sensitivities of NET~300 uKrts. Thanks to the scalability of this design, we have deployed 2560 detectors as 1280 matched pairs in Keck Array with a combined instantaneous sensitivity of ~9 uKrts, as measured directly from CMB maps in the 2013 season. Similar arrays have recently flown in the SPIDER instrument, and development of this technology is ongoing.

From an elevation of ~38 km, the balloon-borne ANtarctic Impulsive Transient Antenna (ANITA) is designed to detect the up-coming radio frequency (RF) signal resulting from a sub-surface neutrino-nucleon collision. Although no neutrinos have been discovered thus far, ANITA is nevertheless the only experiment to self-trigger on radio frequency emissions from cosmic-ray induced atmospheric air showers. In the majority of those cases, down-coming RF signals are observed via their reflection from the Antarctic ice sheet and back up to the ANITA interferometer. Estimating the energy scale of the incident cosmic rays therefore requires an estimate of the fractional power reflected at the air-ice interface. Similarly, inferring the energy of neutrinos interacting in-ice from observations of the upwards-directed signal refracting out to ANITA also requires consideration of signal coherence across the interface. By comparing the direct Solar RF signal intensity measured with ANITA to the surface-reflected Solar signal intensity, as a function of incident elevation angle relative to the surface {\Theta}, we estimate the power reflection coefficients R({\Theta}). We find general consistency between our average measurements and the values of R({\Theta}) expected from the Fresnel equations, separately for horizontal- vs. vertical-polarizations.

We report on the design and performance of the BICEP2 instrument and on its three-year data set. BICEP2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1 to 5 degrees ($\ell$=40-200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The BICEP2 telescope observed from the South Pole with a 26~cm aperture and cold, on-axis, refractive optics. BICEP2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed BICEP2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, BICEP2 achieved an instrument noise-equivalent temperature of 15.8 $\mu$K sqrt(s). The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (5.2 $\mu$K arcmin) over an effective area of 384 square degrees within a 1000 square degree field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales.

BICEP3 is a 550 mm-aperture refracting telescope for polarimetry of radiation in the cosmic microwave background at 95 GHz. It adopts the methodology of BICEP1, BICEP2 and the Keck Array experiments - it possesses sufficient resolution to search for signatures of the inflation-induced cosmic gravitational-wave background while utilizing a compact design for ease of construction and to facilitate the characterization and mitigation of systematics. However, BICEP3 represents a significant breakthrough in per-receiver sensitivity, with a focal plane area 5$\times$ larger than a BICEP2/Keck Array receiver and faster optics ($f/1.6$ vs. $f/2.4$). Large-aperture infrared-reflective metal-mesh filters and infrared-absorptive cold alumina filters and lenses were developed and implemented for its optics. The camera consists of 1280 dual-polarization pixels; each is a pair of orthogonal antenna arrays coupled to transition-edge sensor bolometers and read out by multiplexed SQUIDs. Upon deployment at the South Pole during the 2014-15 season, BICEP3 will have survey speed comparable to Keck Array 150 GHz (2013), and will significantly enhance spectral separation of primordial B-mode power from that of possible galactic dust contamination in the BICEP2 observation patch.

(abridged for arXiv) We report results from the BICEP2 experiment, a cosmic microwave background (CMB) polarimeter specifically designed to search for the signal of inflationary gravitational waves in the B-mode power spectrum around $\ell\sim80$. The telescope comprised a 26 cm aperture all-cold refracting optical system equipped with a focal plane of 512 antenna coupled transition edge sensor 150 GHz bolometers each with temperature sensitivity of $\approx300\mu\mathrm{K}_\mathrm{CMB}\sqrt{s}$. BICEP2 observed from the South Pole for three seasons from 2010 to 2012. A low-foreground region of sky with an effective area of 380 square deg was observed to a depth of 87 nK deg in Stokes $Q$ and $U$. We find an excess of $B$-mode power over the base lensed-LCDM expectation in the range $30< \ell< 150$, inconsistent with the null hypothesis at a significance of $> 5\sigma$. Through jackknife tests and simulations we show that systematic contamination is much smaller than the observed excess. We also examine a number of available models of polarized dust emission and find that at their default parameter values they predict power $\sim(5-10)\times$ smaller than the observed excess signal. However, these models are not sufficiently constrained to exclude the possibility of dust emission bright enough to explain the entire excess signal. Cross correlating BICEP2 against 100 GHz maps from the BICEP1 experiment, the excess signal is confirmed and its spectral index is found to be consistent with that of the CMB, disfavoring dust at $1.7\sigma$. The observed $B$-mode power spectrum is well fit by a lensed-LCDM + tensor theoretical model with tensor-to-scalar ratio $r=0.20^{+0.07}_{-0.05}$, with $r=0$ disfavored at $7.0\sigma$. Accounting for the contribution of foreground dust will shift this value downward by an amount which will be better constrained with upcoming data sets.

This is a report on the status and prospects of the quantification of neutrino properties through the cosmological neutrino background for the Cosmic Frontier of the Division of Particles and Fields Community Summer Study long-term planning exercise. Experiments planned and underway are prepared to study the cosmological neutrino background in detail via its influence on distance-redshift relations and the growth of structure. The program for the next decade described in this document, including upcoming spectroscopic galaxy surveys eBOSS and DESI and a new Stage-IV CMB polarization experiment CMB-S4, will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass measurement will produce a high significance detection of non-zero sum m_nu, whose lower bound derived from atmospheric and solar neutrino oscillation data is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this measurement will definitively rule out the inverted neutrino mass hierarchy, shedding light on one of the most puzzling aspects of the Standard Model of particle physics --- the origin of mass. This precise a measurement of N_eff will allow for high sensitivity to any light and dark degrees of freedom produced in the big bang and a precision test of the standard cosmological model prediction that N_eff = 3.046.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

We present an interferometric technique for the reconstruction of ultra-wide band impulsive signals from point sources. This highly sensitive method was developed for the search for ultra-high energy neutrinos with the ANITA experiment but is fully generalizable to any antenna array detecting radio impulsive events. Applications of the interferometric method include event reconstruction, thermal noise and anthropogenic background rejection, and solar imaging for calibrations. We illustrate this technique with applications from the analysis of the ANITA-I and ANITA-II data in the 200-1200 MHz band.

We consider the implications of Lorentz-invariance violation (LIV) on cosmogenic neutrino observations, with particular focus on the constraints imposed on several well-developed models for ultra-high energy cosmogenic neutrino production by recent results from the Antarctic Impulsive Transient Antenna (ANITA) long-duration balloon payload, and Radio Ice Cherenkov Experiment (RICE) at the South Pole. Under a scenario proposed originally by Coleman and Glashow, each lepton family may attain maximum velocities that can exceed the speed of light, leading to energy-loss through several interaction channels during propagation. We show that future observations of cosmogenic neutrinos will provide by far the most stringent limit on LIV in the neutrino sector. We derive the implied level of LIV required to suppress observation of predicted fluxes from several mainstream cosmogenic neutrino models, and specifically those recently constrained by the ANITA and RICE experiments. We simulate via detailed Monte Carlo code the propagation of cosmogenic neutrino fluxes in the presence of LIV-induced energy losses. We show that this process produces several detectable effects in the resulting attenuated neutrino spectra, even at LIV-induced neutrino superluminality of (u_{\nu}-c)/c ~ 10^{-26}, about 13 orders of magnitude below current bounds.

Between the BICEP2 and Keck Array experiments, we have deployed over 1500 dual polarized antenna coupled bolometers to map the Cosmic Microwave Background's polarization. We have been able to rapidly deploy these detectors because they are completely planar with an integrated phased-array antenna. Through our experience in these experiments, we have learned of several challenges with this technology- specifically the beam synthesis in the antenna- and in this paper we report on how we have modified our designs to mitigate these challenges. In particular, we discus differential steering errors between the polarization pairs' beam centroids due to microstrip cross talk and gradients of penetration depth in the niobium thin films of our millimeter wave circuits. We also discuss how we have suppressed side lobe response with a Gaussian taper of our antenna illumination pattern. These improvements will be used in Spider, Polar-1, and this season's retrofit of Keck Array.

The Keck Array (SPUD) is a set of microwave polarimeters that observes from the South Pole at degree angular scales in search of a signature of Inflation imprinted as B-mode polarization in the Cosmic Microwave Background (CMB). The first three Keck Array receivers were deployed during the 2010-2011 Austral summer, followed by two new receivers in the 2011-2012 summer season, completing the full five-receiver array. All five receivers are currently observing at 150 GHz. The Keck Array employs the field-proven BICEP/BICEP2 strategy of using small, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. We describe our efforts to characterize the main beam shape and beam shape mismatch between co-located orthogonally-polarized detector pairs, and discuss the implications of measured differential beam parameters on temperature to polarization leakage in CMB analysis.

The Keck Array (SPUD) began observing the cosmic microwave background's polarization in the winter of 2011 at the South Pole. The Keck Array follows the success of the predecessor experiments Bicep and Bicep2, using five on-axis refracting telescopes. These have a combined imaging array of 2500 antenna-coupled TES bolometers read with a SQUID-based time domain multiplexing system. We will discuss the detector noise and the optimization of the readout. The achieved sensitivity of the Keck Array is 11.5 {\mu}K_(CMB)*sqrt{s} in the 2012 configuration.

The BICEP2 and Keck Array experiments are designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 2-4 degrees (l=50-100). This is the region in which the B-mode signal, a signature prediction of cosmic inflation, is expected to peak. BICEP2 was deployed to the South Pole at the end of 2009 and is in the middle of its third year of observing with 500 polarization-sensitive detectors at 150 GHz. The Keck Array was deployed to the South Pole at the end of 2010, initially with three receivers--each similar to BICEP2. An additional two receivers have been added during the 2011-12 summer. We give an overview of the two experiments, report on substantial gains in the sensitivity of the two experiments after post-deployment optimization, and show preliminary maps of CMB polarization from BICEP2.

We set the first limits on the ultra-high energy (UHE) neutrino fluence at energies greater than 10^9 GeV from gamma-ray bursts (GRBs) based on data from the second flight of the ANtarctic Impulsive Transient Antenna (ANITA). During the 31 day flight of ANITA-II, 26 GRBs were recorded by Swift or Fermi. Of these, we analyzed the 12 GRBs which occurred during quiet periods when the payload was away from anthropogenic activity. In a blind analysis, we observe 0 events on a total background of 0.0044 events in the combined prompt window for all 12 low-background bursts. We also observe 0 events from the remaining 14 bursts. We place a 90% confidence level limit on the E^-4 prompt neutrino fluence of 2.5x10^17 GeV^3/cm^2 between 10^8 and 10^12 GeV from GRB090107A. This is the first reported limit on the UHE neutrino fluence from GRBs above 10^9 GeV, and the strongest limit above 10^8 GeV.

We have conducted a search for extended energy deposition trails left by ultra-relativistic magnetic monopoles interacting in Antarctic ice. The non-observation of any satisfactory candidates in the 31 days of accumulated ANITA-II flight data results in an upper limit on the diffuse flux of relativistic monopoles. We obtain a 90% C.L. limit of order 10^{-19}/(cm^2-s-sr) for values of Lorentz boost factor 10^{10}<gamma at the anticipated energy E=10^{16} GeV. This bound is stronger than all previously published experimental limits for this kinematic range.

This is an erratum to our paper in Physical Review D82:022004,2010, corresponding to preprint: arXiv:1003.2961 .

The Antarctic Impulsive Transient Antenna (ANITA) completed its second long-duration balloon flight in January 2009, with 31 days aloft (28.5 live days) over Antarctica. ANITA searches for impulsive coherent radio Cherenkov emission from 200 to 1200 MHz, arising from the Askaryan charge excess in ultra-high energy neutrino-induced cascades within Antarctic ice. This flight included significant improvements over the first flight in the payload sensitivity, efficiency, and a flight trajectory over deeper ice. Analysis of in-flight calibration pulses from surface and sub-surface locations verifies the expected sensitivity. In a blind analysis, we find 2 surviving events on a background, mostly anthropogenic, of 0.97+-0.42 events. We set the strongest limit to date for 1-1000 EeV cosmic neutrinos, excluding several current cosmogenic neutrino models.

We report the observation of sixteen cosmic ray events of mean energy of 1.5 x 10^{19} eV, via radio pulses originating from the interaction of the cosmic ray air shower with the Antarctic geomagnetic field, a process known as geosynchrotron emission. We present the first ultra-wideband, far-field measurements of the radio spectral density of geosynchrotron emission in the range from 300-1000 MHz. The emission is 100% linearly polarized in the plane perpendicular to the projected geomagnetic field. Fourteen of our observed events are seen to have a phase-inversion due to reflection of the radio beam off the ice surface, and two additional events are seen directly from above the horizon.