Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

We explore data analysis techniques for signatures from heavy particle production during inflation. Heavy particules can be produced by time dependent masses and couplings, which are ubiquitous in string theory. These localized excitations induce curvature perturbations with non-zero correlation functions at all orders. In particular, Flauger et. al. 2016 has shown that the signal-to-noise as a function of the order $N$ of the correlation function can peak for $N$ of order $\mathcal{O}(1)$ to $\mathcal{O}(100)$ for an interesting space of models. As previous non-Gaussianity analyses have focused on $N=\{3,4\}$, in principle this provides an unexplored data analysis window with new discovery potential. We derive estimators for arbitrary $N$-point functions in this model and discuss their properties and covariances. To lowest order, the heavy particle production phenomenology reduces to a classical Poisson process, which can be implemented as a search for spherically symmetric profiles in the curvature perturbations. We explicitly show how to recover this result from the $N$-point functions and their estimators. Our focus in this paper is on method development, but we provide an initial data analysis using WMAP data, which illustrates the particularities of higher $N$-point function searches.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

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.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

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

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged).

We show how a neural network can be trained to Wiener filter masked CMB maps to high accuracy. We propose an innovative neural network architecture, the WienerNet, which guarantees linearity in the data map. Our method does not require Wiener filtered training data, but rather learns Wiener filtering from tailored loss functions which are mathematically guaranteed to be minimized by the exact solution. Once trained, the neural network Wiener filter is extremely fast, about a factor of 1000 faster than the standard conjugate gradient method. Wiener filtering is the computational bottleneck in many optimal CMB analyses, including power spectrum estimation, lensing and non-Gaussianities, and our method could potentially be used to speed them up by orders of magnitude with minimal loss of optimality. The method should also be useful to analyze other statistical fields in cosmology.

Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity.

Kinetic Sunyaev Zel'dovich (kSZ) tomography provides a powerful probe of the radial velocity field of matter in the Universe. By cross-correlating a high resolution CMB experiment like CMB S4 and a galaxy survey like DESI or LSST, one can measure the radial velocity field with very high signal to noise over a large volume of the universe. In this paper we show how this measurement can be used to improve constraints on primordial non-Gaussianities of the local type. The velocity field provides a measurement of the unbiased matter perturbations on large scales, which can be cross-correlated with the biased large-scale galaxy density field. This results in sample variance cancellation for a measurement of scale-dependent bias due to a non-zero $f_{NL}$. Using this method we forecast that CMB S4 and LSST combined reach a sensitivity $\sigma_{f_{NL}} \sim 0.5$, which is a factor of three improvement over the sensitivity using LSST alone (without internal sample variance cancellation). We take into account critical systematics like photometric redshifts, the kSZ optical depth degeneracy, and systematics affecting the shape of the galaxy auto-power spectrum and find that these have negligible impact, thus making kSZ tomography a robust probe for primordial non-Gaussianities. We also forecast the impact of mass binning on our constraints. The techniques proposed in this paper could be an important component of achieving the theoretically important threshold of $\sigma_{f_{NL}} \lesssim 1$ with future surveys.

Several statistics have been proposed for measuring the kSZ effect by combining the small-scale CMB with galaxy surveys. We review five such statistics, and show that they are all mathematically equivalent to the optimal bispectrum estimator of type $\langle ggT \rangle$. Reinterpreting these kSZ statistics as special cases of bispectrum estimation makes many aspects transparent, for example optimally weighting the estimator, or incorporating photometric redshift errors. We analyze the information content of the bispectrum and show that there are two observables: the small-scale galaxy-electron power spectrum $P_{ge}(k_S)$, and the large-scale galaxy-velocity power spectrum $P_{gv}(k)$. The cosmological constraining power of the kSZ arises from its sensitivity to fluctuations on large length scales, where its effective noise level can be much better than galaxy surveys.

The detection and characterization of primordial gravitational waves through their impact on the polarization anisotropies of the cosmic microwave background (CMB) is a primary science goal of current and future observations of the CMB. An ancillary dataset that will become accessible with the great leaps in sensitivity of CMB experiments is the polarized Sunyaev Zel'dovich (pSZ) effect, small-scale CMB polarization anisotropies induced by scattering from free electrons in the post-reionization Universe. The cross correlation of the pSZ effect with galaxy surveys, a technique known as pSZ tomography, can be used to reconstruct the remote quadrupole field: the CMB quadrupole observed from different locations in the Universe. Primordial gravitational waves leave a distinct imprint on the remote quadrupole field, making pSZ tomography a potential new method to characterize their properties. Building on previous work, we explore the utility of the full set of correlations between the primary CMB and the reconstructed remote quadrupole field to both provide exclusion limits on the amplitude of primordial gravitational waves, as well as to provide constraints on several phenomenological models of the tensor sector: axion gauge field inflation, general models with chiral tensors, and models with modified late-time decay of tensors. We find that relatively futuristic experimental requirements are necessary to provide competitive exclusion limits compared with the primary CMB. However, pSZ tomography can be a powerful probe of the late-time evolution of tensors and, through cross-correlations with the primary CMB, can provide mild improvements on parameter constraints in various models with chiral primordial gravitational waves.

The kinetic Sunyaev Zel'dovich (kSZ) and polarized Sunyaev Zel'dovich (pSZ) effects are temperature and polarization anisotropies induced by the scattering of CMB photons from structure in the post-reionization Universe. In the case of the kSZ effect, small angular scale anisotropies in the optical depth are modulated by the cosmic microwave background (CMB) dipole field, i.e. the CMB dipole observed at each spacetime point, which is sourced by the primordial dipole and especially the local peculiar velocity. In the case of the pSZ effect, similar small-scale anisotropies are modulated by the CMB quadrupole field, which receives contributions from both scalar and tensor modes. Statistical anisotropies in the cross correlations of CMB temperature and polarization with tracers of the inhomogeneous distribution of electrons provide a means of isolating and reconstructing the dipole and quadrupole fields. In this paper, we present a set of unbiased minimum variance quadratic estimators for the reconstruction of the dipole and quadrupole fields, and forecast the ability of future CMB experiments and large scale structure surveys to perform this reconstruction. Consistent with previous work, we find that a high fidelity reconstruction of the dipole and quadrupole fields over a variety of scales is indeed possible, and demonstrate the sensitivity of the pSZ effect to primordial tensor modes. Using a principle component analysis, we estimate how many independent modes could be accessed in such a reconstruction. We also comment on a few first applications of a detection of the dipole and quadrupole fields, including a reconstruction of the primordial contribution to our locally observed CMB dipole, a test of statistical homogeneity on large scales from the first modes of the quadrupole field, and a reconstruction technique for the primordial potential on the largest scales.

Secondary CMB polarization is induced by the late-time scattering of CMB photons by free electrons on our past light cone. This polarized Sunyaev Zel'dovich (pSZ) effect is sensitive to the electrons' locally observed CMB quadrupole, which is sourced primarily by long wavelength inhomogeneities. By combining the remote quadrupoles measured by free electrons throughout the Universe after reionization, the pSZ effect allows us to obtain additional information about large scale modes beyond what can be learned from our own last scattering surface. Here we determine the power of pSZ tomography, in which the pSZ effect is cross-correlated with the density field binned at several redshifts, to provide information about the long wavelength Universe. The signal we explore here is a power asymmetry in the cross-correlation between $E$ or $B$ mode CMB polarization and the density field. We compare this to the cosmic variance limited noise: the random chance to get a power asymmetry in the absence of a large scale quadrupole field. By computing the necessary transfer functions and cross-correlations, we compute the signal-to-noise ratio attainable by idealized next generation CMB experiments and galaxy surveys. We find that a signal-to-noise ratio of $\sim 1-10$ is in principle attainable over a significant range of power multipoles, with the strongest signal coming from the first multipoles in the lowest redshift bins. These results prompt further assessment of realistically measuring the pSZ signal and the potential impact for constraining cosmology on large scales.

It is generally expected that heavy fields are present during inflation, which can leave their imprint in late-time cosmological observables. The main signature of these fields is a small amount of distinctly shaped non-Gaussianity, which if detected, would provide a wealth of information about the particle spectrum of the inflationary Universe. Here we investigate to what extent these signatures can be detected or constrained using futuristic 21-cm surveys. We construct model-independent templates that extract the squeezed-limit behavior of the bispectrum, and examine their overlap with standard inflationary shapes and secondary non-Gaussianities. We then use these templates to forecast detection thresholds for different masses and couplings using a 3D reconstruction of modes during the dark ages ($z\sim 30-100$). We consider interactions of several broad classes of models and quantify their detectability as a function of the baseline of a dark ages interferometer. Our analysis shows that there exists the tantalizing possibility of discovering new particles with different masses and interactions with future 21-cm surveys.

Detecting a deviation from a featureless primordial power spectrum of fluctuations would give profound insight into the physics of the primordial Universe. Depending on their nature, primordial features can either provide direct evidence for the inflation scenario or pin down details of the inflation model. Thus far, using the cosmic microwave background (CMB) we have only been able to put stringent constraints on the amplitude of features, but no significant evidence has been found for such signals. Here we explore the limit of the experimental reach in constraining such features using 21 cm tomography at high redshift. A measurement of the 21 cm power spectrum from the Dark Ages is generally considered as the ideal experiment for early Universe physics, with potentially access to a large number of modes. We consider three different categories of theoretically motivated models: the sharp feature models, resonance models, and standard clock models. We study the improvements on bounds on features as a function of the total number of observed modes and identify parameter degeneracies. The detectability depends critically on the amplitude, frequency and scale-location of the features, as well as the angular and redshift resolution of the experiment. We quantify these effects by considering different fiducial models. Our forecast shows that a cosmic variance limited 21 cm experiment measuring fluctuations in the redshift range $30\leq z \leq 100$ with a 0.01-MHz bandwidth and sub-arcminute angular resolution could potentially improve bounds by several orders of magnitude for most features compared to current Planck bounds. At the same time, 21 cm tomography also opens up a unique window into features that are located on very small scales.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

We develop the tools necessary to assess the statistical significance of resonant features in the CMB correlation functions, combining power spectrum and bispectrum measurements. This significance is typically addressed by running a large number of simulations to derive the probability density function (PDF) of the feature-amplitude in the Gaussian case. Although these simulations are tractable for the power spectrum, for the bispectrum they require significant computational resources. We show that, by assuming that the PDF is given by a multi-variate Gaussian where the covariance is determined by the Fisher matrix of the sine and cosine terms, we can efficiently produce spectra that are statistically close to those derived from full simulations. By drawing a large number of spectra from this PDF, both for the power spectrum and the bispectrum, we can quickly determine the statistical significance of candidate signatures in the CMB, considering both single frequency and multi-frequency estimators. We show that for resonance models, cosmology and foreground parameters have little influence on the estimated amplitude, which allows to simplify the analysis considerably. A more precise likelihood treatment can then be applied to candidate signatures only. We also discuss a modal expansion approach for the power spectrum, aimed at quickly scanning through large families of oscillating models.

We propose optimal estimators for bispectra from excited states. Two common properties of such bispectra are the enhancement in the collinear limit, and the prediction of oscillating features. We review the physics behind excited states and some of the choices made in the literature. We show that the enfolded template is a good template in the collinear limit, but does poorly elsewhere, establishing a strong case for an improved estimator. Although the detailed scale dependence of the bispectra differs depending on various assumptions, generally the predicted bispectra are either effectively 1 or 2-dimensional and a simple Fourier basis suffices for accurate reconstruction. For an optimal CMB data analysis, combining all n-point functions, the choice for the excited state needs to be the same when computing power spectrum, bispectrum and higher order correlation functions. This has not always been the case, which could lead to wrong conclusions. We calculate the bispectrum for different choices previously discussed for the power spectrum, setting up a consistent framework to search for evidence of excited states in the CMB data.

We propose an (optimal) estimator for a CMB bispectrum containing logarithmically spaced oscillations. There is tremendous theoretical interest in such bispectra, and they are predicted by a plethora of models, including axion monodromy models of inflation and initial state modifications. The number of resolved logarithmical oscillations in the bispectrum is limited due to the discrete resolution of the multipole bispectrum. We derive a simple relation between the maximum number of resolved oscillations and the frequency. We investigate several ways to factorize the primordial bispectrum, and conclude that a one dimensional expansion in the sum of the momenta $\sum k_i = k_t$ is the most efficient and flexible approach. We compare the expansion to the exact result in multipole space and show for $\omega_{\rm eff}=100$ that $\mathcal{O}(10^3)$ modes are sufficient for an accurate reconstruction. We compute the expected $\sigma_{f_{\rm NL}}$ and find that within an effective field theory (EFT) the overall signal to noise scales as $S/N\propto \omega^{3/2}$. Using only the temperature data we find $S/N\sim\mathcal{O}(1-10^2)$ for the frequency domain set by the EFT.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

Oscillating shapes of the primordial bispectrum are present in many inflationary models. The Planck experiment has recently published measurements of oscillating shapes, which were however limited to the efficient frequency range of the used analysis method. Here we study the KSW estimator for oscillations in the CMB bispectrum, that allows to examine arbitrary oscillation frequencies for separable oscillating bispectrum shapes. We study the precision with which amplitude, phase and frequency can be determined with our estimator. An examination of the 3-point function in real space gives further insight into the estimator.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.