We establish strong gravitational lens systems as robust probes of axion-like particles (ALPs) -- a candidate for dark matter. A tiny interaction of photons with ALPs induces birefringence. Multiple images of gravitationally lensed polarised objects allow differential birefringence measurement, alleviating systematics and astrophysical dependencies. We apply this novel method to the lens system CLASS B1152+199 and constrain ALP-photon coupling $\le 9.2\times 10^{-11}\, {\rm GeV}^{-1} \textrm{ to } 7.7\times 10^{-8}\, {\rm GeV}^{-1}$ ($95\%$ C.L.) for ALP mass between $3.6\times 10^{-21} \,{\rm eV}$ and $4.6 \times 10^{-18} \,{\rm eV}$. A larger sample will improve the constraints.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

We present a concise derivation of geometric optics in the presence of axionic fields in a curved space-time. Whenever light can be described via geometric optics (the eikonal approximation), the only difference to the situation without axionic field is the phenomenon of achromatic birefringence. Consequently, redshift of light and distance estimates based on propagating light rays, as well as shear and magnification due to gravitational lensing are not affected by the interaction of light with an axionic field.

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known.

The new generation of high-resolution broad-band radio telescopes, like the Low Frequency Array (LOFAR), produces, depending on the level of compression, between 1 to 10 TB of data per hour after correlation. Such a large amount of scientific data demand powerful computing resources and efficient data handling strategies to be mastered. The LOFAR Two-metre Sky Survey (LoTSS) is a Key Science Project (KSP) of the LOFAR telescope. It aims to map the entire northern hemisphere at unprecedented sensitivity and resolution. The survey consist of 3 168 pointings, requiring about 30 PBytes of storage space. As a member of the German Long Wavelength Consortioum (GLOW) the Forschungszentrum J\"ulich (FSZ) stores in the Long Term Archive (LTA) about 50% of all LoTSS observations conducted to date. In collaboration with SURFsara in Amsterdam we developed service tools that enable the KSP to process LOFAR data stored in the J\"ulich LTA at the supercomputer JUWELS in an automated and robust fashion. Through our system more than 500 out of 800 existing LoTSS observations have already been processed with the prefactor pipeline. This pipeline calibrates the direction-independent instrumental and ionospheric effects and furthermore reduces the data size significantly. For continuum imaging, this processing pipeline is the standard pipeline that is executed before more advanced processing and image reconstruction methods are applied.

The LOFAR Two-metre Sky Survey (LoTSS) will map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. We study the completeness of the LoTSS first data release (DR1) and find a point-source completeness of 99 % above flux densities of 0.8 mJy and define a suite of quality cuts. We determine the count-in-cell statistics and differential source counts statistic and measure the angular two-point correlation function of the LoTSS radio sources. The counts-in-cell statistic reveals that the distribution of radio sources cannot be described by a spatial Poisson process. Instead, a good fit is provided by a compound Poisson distribution. The differential source counts are in good agreement with previous findings in deep fields at low radio frequencies. Simulated catalogues from the SKA design study sky match with our findings, in contrast to the more recent T-RECS models. For angular separations between 0.1 deg and 1 deg we find self-consistency amongst different estimates of the angular two-point correlation. The angular two-point correlation is $<10^{-2}$ at angular scales $> 1$ deg and up to the largest scales probed. We see indications of unidentified systematic issues at angular scales $>1$ deg, likely due to data calibration issues. Based on the distribution of photometric redshifts of LoTSS sources and the Planck 2018 best-fit cosmological model, the theoretically predicted angular two-point correlation between 0.1 deg and 1 deg agrees with the measured clustering for a subsample of radio sources with redshift information. We find agreement with the expectation of large-scale statistical isotropy of the radio sky at the per cent level and below angular separations of 1 deg the angular two-point correlation agrees with the expectation of the cosmological standard model. (abbreviated)

We propose that gravitational interactions of cosmic neutrinos with the statistically homogeneous and isotropic fluctuations of space-time leads to decoherence. This working hypothesis, which we describe by means of a Lindblad operator, is applied to the system of two- and three-flavour neutrinos undergoing vacuum oscillations and the consequences are investigated. As result of this decoherence we find that the neutrino entropy would increase as a function of initial spectral distortions, mixing angles and CP-violation phase. Subsequently we discuss the chances to discover such an increase observationally (in principle). We also present the expected flavour composition of the cosmic neutrino background after decoherence is completed. The physics of two- or three-flavour oscillation of cosmological neutrinos resembles in many aspects two- or three-level systems in atomic clocks, which were recently proposed by Weinberg for the study of decoherence phenomena.

We describe Cal X-1, a SmallSat mission concept to establish X-ray standard candles in the sky, which will enable absolute calibration of the current and future X-ray observatories. It consists of two CubeSats flying in formation, one containing an absolutely calibrated X-ray source and another, 1-2 km away, a simple X-ray telescope.

Precise measurement of the foreground synchrotron emission, which contaminates the faint polarized cosmic microwave background radiation (CMB), is a major challenge for the next-generation of CMB experiments. To address this, dedicated foreground measurement experiments are being undertaken at radio frequencies between 2 and 40 GHz. Foreground polarized synchrotron emission measurements are particularly challenging, primarily due to the complicated frequency dependence in the presence of Faraday rotation, and are best recovered through broad fractional-bandwidth polarization measurements at frequencies $\lesssim 5$ GHz. A unique opportunity for measuring the foreground polarized synchrotron emission will be provided by the 15-m SKA-MPG telescope operating in the frequency range 1.7 to 3.5~GHz (S-Band). Here, we present the scope of a Southern sky survey in S-Band at 1 degree angular resolution and explore its added advantage for application of powerful techniques, such as, Stokes $Q$, $U$ fitting and RM-synthesis. A full Southern-sky polarization survey with this telescope, when combined with other on-going efforts at slightly higher frequencies, will provide an excellent frequency coverage for modeling and extrapolating the foreground polarized synchrotron emission to CMB frequencies ($\gtrsim80$~GHz) with rms brightness temperature better than 10~nK per 1 degree$^2$. We find that this survey will be crucial for understanding the effects of Faraday depolarization, especially in low Galactic latitude regions. This will allow better foreground cleaning and thus will contribute significantly in further improving component separation analyses and increase usable sky area for cosmological analysis of the \textit{Planck} data, and the \textit{LiteBIRD} mission in the future.

We study the impact of non-standard momentum distributions of cosmic neutrinos on the anisotropy spectrum of the cosmic microwave background and the matter power spectrum of the large scale structure. We show that the neutrino distribution has almost no unique observable imprint, as it is almost entirely degenerate with the effective number of neutrino flavours, $N_{\mathrm{eff}}$, and the neutrino mass, $m_{\nu}$. Performing a Markov chain Monte Carlo analysis with current cosmological data, we demonstrate that the neutrino mass bound heavily depends on the assumed momentum distribution of relic neutrinos. The message of this work is simple and has to our knowledge not been pointed out clearly before: Cosmology allows that neutrinos have larger masses if their average momentum is larger than that of a perfectly thermal distribution. Here we provide an example in which the mass limits are relaxed by a factor of two.

Multipole vectors and pseudoentropies provide powerful tools for a numerically fast and vivid investigation of possible statistically anisotropic, respectively non-Gaussian signs in CMB temperature fluctuations. After reviewing and linking these two conceptions we compare their application to data analysis using the Planck 2015 NILC full sky map.

The dipole anisotropy seen in the {cosmic microwave background radiation} is interpreted as due to our peculiar motion. The Cosmological Principle implies that this cosmic dipole signal should also be present, with the same direction, in the large-scale distribution of matter. Measurement of the cosmic matter dipole constitutes a key test of the standard cosmological model. Current measurements of this dipole are barely above the expected noise and unable to provide a robust test. Upcoming radio continuum surveys with the SKA should be able to detect the dipole at high signal to noise. We simulate number count maps for SKA survey specifications in Phases 1 and 2, including all relevant effects. Nonlinear effects from local large-scale structure contaminate the {cosmic (kinematic)} dipole signal, and we find that removal of radio sources at low redshift ($z\lesssim 0.5$) leads to significantly improved constraints. We forecast that the SKA could determine the kinematic dipole direction in Galactic coordinates with an error of $(\Delta l,\Delta b)\sim(9^\circ,5^\circ)$ to $(8^\circ, 4^\circ)$, depending on the sensitivity. The predicted errors on the relative speed are $\sim 10\%$. These measurements would significantly reduce the present uncertainty on the direction of the radio dipole, and thus enable the first critical test of consistency between the matter and CMB dipoles.

Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by assuming a thin layer model. For this aim, ionospheric maps of the free electron density (based on Global Positioning System data) and semi-empirical geomagnetic models are needed. Through the data of five highly polarized pulsars observed with the individual German LOw-Frequency ARray stations, we investigate the performances of the ionospheric modelling. In addition, we estimate the parameters of the systematics and the correlated noise generated by the residual unmodelled ionospheric effects, and show the comparison of the different free-electron density maps. For the best ionospheric maps, we have found that the rotation measure corrections on one-year timescales after subtraction of diurnal periodicity are accurate to $\sim$ 0.06--0.07 rad m$^{-2}$.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release.

The LOFAR Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey of the Northern sky. The LoTSS First Data Release (DR1) presents 424 square degrees of radio continuum observations over the HETDEX Spring Field (10h45m00s $<$ right ascension $<$ 15h30m00s and 45$^\circ$00$'$00$'$ $<$ declination $<$ 57$^\circ$00$'$00$''$) with a median sensitivity of 71$\mu$Jy/beam and a resolution of 6$''$. In this paper we present photometric redshifts (photo-$z$) for 94.4% of optical sources over this region that are detected in the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 3$\pi$ steradian survey. Combining the Pan-STARRS optical data with mid-infrared photometry from the Wide-field Infrared Survey Explorer, we estimate photo-$z$s using a novel hybrid photometric redshift methodology optimised to produce the best possible performance for the diverse sample of radio continuum selected sources. For the radio-continuum detected population, we find an overall scatter in the photo-$z$ of 3.9% and an outlier fraction ($\left | z_{\rm{phot}} - z_{\rm{spec}} \right | / (1+z_{\rm{spec}}) > 0.15$) of 7.9%. We also find that, at a given redshift, there is no strong trend in photo-$z$ quality as a function of radio luminosity. However there are strong trends as a function of redshift for a given radio luminosity, a result of selection effects in the spectroscopic sample and/or intrinsic evolution within the radio source population. Additionally, for the sample of sources in the LoTSS First Data Release with optical counterparts, we present rest-frame optical and mid-infrared magnitudes based on template fits to the consensus photometric (or spectroscopic when available) redshift.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

We investigate how a lepton asymmetry impacts the cosmic trajectory in the QCD phase diagram. We study the evolution of chemical potentials during the QCD epoch of the early Universe using susceptibilities from lattice QCD to interpolate between an ideal quark gas and an ideal hadron resonance gas. The lepton asymmetry affects the evolution of all chemical potentials. The standard cosmic trajectory is obtained assuming tiny lepton and baryon asymmetries. For larger lepton asymmetry, the charge chemical potential exceeds the baryon chemical potential before pion annihilation.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

The statistical cosmological principle states that observables on the celestial sphere are sampled from a rotationally invariant distribution. Previously certain large scale anomalies which violate this principle have been found, for example an alignment of the lowest multipoles with the cosmic dipole direction. In this work we continue the search for possible anomalies using multipole vectors which represent a convenient tool for this purpose. In order to study the statistical behavior of multipole vectors, we revisit several construction methods. We investigate all four full-sky foreground-cleaned maps from the Planck 2015 release with respect to four meaningful physical directions using computationally cheap statistics that have a simple geometric interpretation. We find that the full-sky SEVEM map deviates from all the other cleaned maps, as it shows a strong correlation with the Galactic Pole and Galactic Center. The other three maps COMMANDER, NILC and SMICA show a consistent behavior. On the largest angular scales, $l \leq 5$, as well as on intermediate scales, $l = 20, 21, 22, 23, 24$, all of them are unusually correlated with the cosmic dipole direction. These scales coincide with the scales on which the angular power spectrum deviates from the Planck 2015 best-fit {\Lambda}CDM model. In the range $2 \leq l \leq 50$ as a whole there is no unusual behavior visible globally. We do not find abnormal intramultipole correlation, i.e. correlation of multipole vectors inside a given multipole without reference to any outer direction.

We study the prospects to measure the cosmic radio dipole by means of continuum surveys with the Square Kilometre Array. Such a measurement will allow a critical test of the cosmological principle. It will test whether the cosmic rest frame defined by the cosmic microwave background at photon decoupling agrees with the cosmic rest frame of matter at late times.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

PSR J0922+0638 (B0919+06) shows unexplained anomalous variations in the on-pulse phase, where the pulse appears to episodically move to an earlier longitude for a few tens of rotations before reverting to the usual phase for approximately several hundred to more than a thousand rotations. These events, where the pulse moves in phase by up to 5$^{\circ}$, have been previously detected in observations from $\sim$300 to 2000 MHz. We present simultaneous observations from the Effelsberg 100-m radio telescope at 1350 MHz and the Bornim (Potsdam) station of the LOw Frequency ARray at 150 MHz. Our observations present the first evidence for an absence of the anomalous phase-shifting behaviour at 150 MHz. Instead, the observed intensity at the usual pulse-phase typically decreases, often showing a pseudo-nulling feature corresponding to the times when phase shifts are observed at 1350 MHz. The presence of weak emission at the usual pulse-phase supports the theory that these shifts may result from processes similar to the 'profile-absorption' expected to operate for PSR J0814+7429 (B0809+74). A possible mechanism for this could be intrinsic variations of the emission within the pulsar's beam combined with absorption by expanding shells of electrons in the line of sight.

Radio astronomy has experienced phenomenal progress in recent years due to advances in digital technologies and processing speed, the development of new technologies, and the prospect for new powerful facilities A new generation of radio interferometers is opening new windows to the Universe. LOFAR, the world largest telescope, extended the classical radio window at low frequencies. ALMA started a new era at millimetre wavelengths and the planned SKA will revolutionize the sciences of the Universe, well beyond the traditional limits of astronomy. Radio astronomy research is making leaps with enhancements in resolution, sensitivity, and image fidelity. In Germany, the community has access to excellent facilities and training opportunities. The Effelsberg telescope remains the flagship of radio astronomical research at centimeter wavelengths and serves as a test bed for new technologies, LOFAR has notably expanded the German community and now includes six German stations and a LOFAR long term archive coordinated by the GLOW consortium. The SKA will be a transformational astronomical facility in the coming decade(s), and the German community is looking forward to broadly participate in SKA-enabled research.

We study the propagation of light bundles in non-empty spacetime, as most of the Universe is filled by baryonic matter in the form of a (dilute) plasma. Here we restrict to the case of a cold (i.e., pressureless) and non-magnetised plasma. Then the influence of the medium on the light rays is encoded in the spacetime dependent plasma frequency. Our result for a general spacetime generalises the Sachs equations to the case of a cold plasma Universe. We find that the reciprocity law (Etherington theorem), the relation that connects area distance with luminosity distance, is modified. Einstein's field equation is not used, i.e., our results apply independently of whether or not the plasma is self-gravitating. As an example, our findings are applied to a homogeneous plasma in a Robertson-Walker spacetime. We find small modifications of the cosmological redshift of frequencies and of the Hubble law.

A cosmic lepton asymmetry $\eta_{\text{l}}=(n_{\text{l}}-n_{\bar{\text{l}}})/n_{\gamma}$ affects the primordial helium abundance and the expansion rate of the early Universe. Both of these effects have an impact on the anisotropies of the cosmic microwave background (CMB). We derive constraints on the neutrino chemical potentials from the Planck 2015 data, assuming equal lepton flavour asymmetries and negligible neutrino masses. We find $\xi=-0.002 ^{+0.114}_{-0.111}$ (95\% CL) for the chemical potentials, which corresponds to $ -0.085 \leq \eta_{\text{l}} \leq 0.084$. Our constraints on the lepton asymmetry are significantly stronger than previous constraints from CMB data analysis and we argue that they are more robust than those from primordial light element abundances. The resulting constraints on the primordial helium and deuterium abundances are consistent with those from direct measurements.

Context. The Sun is an active star that produces large-scale energetic events such as solar flares and coronal mass ejections and numerous smaller-scale events such as solar jets. These events are often associated with accelerated particles that can cause emission at radio wavelengths. The reconfiguration of the solar magnetic field in the corona is believed to be the cause of the majority of solar energetic events and accelerated particles. Aims. Here, we investigate a bright J-burst that was associated with a solar jet and the possible emission mechanism causing these two phenomena. Methods. We used data from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph (NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July 2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic components and it was associated with a solar jet observed at extreme ultraviolet wavelengths with SDO. The solar jet occurred at a time and location coincident with the radio burst, in the northern hemisphere, and not inside a group of complex active regions in the southern hemisphere. The jet occurred in the negative polarity region of an area of bipolar plage. Newly emerged positive flux in this region appeared to be the trigger of the jet. Conclusions. Magnetic reconnection between the overlying coronal field lines and the newly emerged positive field lines is most likely the cause of the solar jet. Radio imaging provides a clear association between the jet and the J-burst which shows the path of the accelerated electrons.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

Cosmic reionization and dark matter decay can impact observations of the cosmic microwave sky in a similar way. A simultaneous study of both effects is required to constrain unstable dark matter from cosmic microwave background observations. We compare two reionization models with and without dark matter decay. We find that a reionization model that fits also data from quasars and star forming galaxies results in tighter constraints on the reionization optical depth $\tau_{\text{reio}}$, but weaker constraints on the spectral index $n_{\text{s}}$ than the conventional parametrization. We use the Planck 2015 data to constrain the effective decay rate of dark matter to $\Gamma_{\rm eff} < 2.9 \times 10^{-25}/$s at $95$\% C.L. This limit is robust and model independent. It holds for any type of decaying dark matter and it depends only weakly on the chosen parametrization of astrophysical reionization. For light dark matter particles that decay exclusively into electromagnetic components this implies a limit of $\Gamma < 5.3 \times 10^{-26}/$s at $95$\% C.L. Specifying the decay channels, we apply our result to the case of keV-mass sterile neutrinos as dark matter candidates and obtain constraints on their mixing angle and mass, which are comparable to the ones from the diffuse X-ray background.

We measure the angular two-point correlation and angular power spectrum from the NRAO VLA Sky Survey (NVSS) of radio galaxies. They are found to be consistent with the best-fit cosmological model from the Planck analysis, and with the redshift distribution obtained from the Combined EIS-NVSS Survey Of Radio Sources (CENSORS). Our analysis is based on an optimal estimation of the two-point correlation function and makes use of a new mask, that takes into account direction dependent effects of the observations, sidelobe effects of bright sources and galactic foreground. We also set a flux threshold and take the cosmic radio dipole into account. The latter turns out to be an essential step in the analysis. This improved cosmological analysis of the NVSS emphasizes the importance of a flux calibration that is robust and stable on large angular scales for future radio continuum surveys.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

Several unexpected features have been observed in the microwave sky at large angular scales, both by WMAP an by Planck. Among those features is a lack of both variance and correlation on the largest angular scales, alignment of the lowest multipole moments with one another and with the motion and geometry of the Solar System, a hemispherical power asymmetry or dipolar power modulation, a preference for odd parity modes and an unexpectedly large cold spot in the Southern hemisphere. The individual p-values of the significance of these features are in the per mille to per cent level, when compared to the expectations of the best-fit inflationary $\Lambda$CDM model. Some pairs of those features are demonstrably uncorrelated, increasing their combined statistical significance and indicating a significant detection of CMB features at angular scales larger than a few degrees on top of the standard model. Despite numerous detailed investigations, we still lack a clear understanding of these large-scale features, which seem to imply a violation of statistical isotropy and scale invariance of inflationary perturbations. In this contribution we present a critical analysis of our current understanding and discuss several ideas of how to make further progress.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

We revisit the alignments of the largest structures observed in the cosmic microwave background (CMB) using the seven and nine-year WMAP and first-year Planck data releases. The observed alignments -- the quadrupole with the octopole and their joint alignment with the direction of our motion with respect to the CMB (the dipole direction) and the geometry of the Solar System (defined by the Ecliptic plane) -- are generally in good agreement with results from the previous WMAP data releases. However, a closer look at full-sky data on the largest scales reveals discrepancies between the earlier WMAP data releases (three to seven-year) and the final nine-year release. There are also discrepancies between all the WMAP data releases and the first-year Planck release. Nevertheless, both the WMAP and Planck data confirm the alignments of the largest observable CMB modes in the Universe. In particular, the p-values for the mutual alignment between the quadrupole and octopole, and the alignment of the plane defined by the two with the dipole direction, are both at the greater than 3-sigma level for all three Planck maps studied. We also calculate conditional statistics on the various alignments and find that it is currently difficult to unambiguously identify a leading anomaly that causes the others or even to distinguish correlation from causation.

The lack of large-angle correlations in the observed microwave background temperature fluctuations persists in the final-year maps from WMAP and the first cosmological data release from Planck. We find a statistically robust and significant result: p-values for the missing correlations lying below 0.24 per cent (i.e. evidence at more than 3 sigma) for foreground cleaned maps, in complete agreement with previous analyses based upon earlier WMAP data. A conservative cut-sky analysis of the Planck HFI 100 GHz frequency band, the `cleanest CMB channel', returns a p-value as small as 0.03 per cent, based on the conservative mask defined by WMAP. These findings are in stark contrast to expectations from the inflationary Lambda cold dark matter model and still lack a convincing explanation. If this lack of large-angle correlations is a true feature of our Universe, and not just a statistical fluke, then the cosmological dipole must be considerably smaller than that predicted in the best-fitting model.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

We investigate the differential number counts of sources in radio continuum surveys, including all terms at linear order in cosmological perturbations. Our framework does not assume a specific gauge condition. This general approach allows us to recover gauge invariance explicitly. With the complete derivations of the covariant volume integral on the past light cone, we have identified several contributions in the number counts. To clarify their underlying physics, we present each contributions in terms of scalar, vector and tensor modes. This theoretical framework promises to be widely applicable to continuum radio galaxy surveys to model the expected angular power spectrum and two-point correlation.

Predictions of the microwave thermal emission from the interplanetary dust cloud are made using several contemporary meteoroid models to construct the distributions of cross-section area of dust in space, and applying the Mie light-scattering theory to estimate the temperatures and emissivities of dust particles in broad size and heliocentric distance ranges. In particular, the model of the interplanetary dust cloud by Kelsall et al. (1998, ApJ 508: 44-73), the five populations of interplanetary meteoroids of Divine (1993, JGR 98(E9): 17,029-17,048) and the Interplanetary Meteoroid Engineering Model (IMEM) by Dikarev et al. (2004, EMP 95: 109-122) are used in combination with the optical properties of olivine, carbonaceous and iron spherical particles. The Kelsall model has been widely accepted by the Cosmic Microwave Background (CMB) community. We show, however, that it predicts the microwave emission from interplanetary dust remarkably different from the results of application of the meteoroid engineering models. We make maps and spectra of the microwave emission predicted by the three models assuming variant composition of dust particles. Predictions can be used to look for the emission from interplanetary dust in CMB experiments as well as to plan new observations.

The imprint of baryon acoustic oscillations (BAO) in large-scale structure can be used as a standard ruler for mapping out the cosmic expansion history, and hence for testing cosmological models. In this article we briefly describe the scientific background to the BAO technique, and forecast the potential of the Phase 1 and 2 SKA telescopes to perform BAO surveys using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys. We find that a 25,000 sq. deg. intensity mapping survey on a Phase 1 array will preferentially constrain the radial BAO, providing a highly competitive 2% constraint on the expansion rate at z ~ 2. A 30,000 sq. deg. galaxy redshift survey on SKA2 will outperform all other planned experiments for z < 1.4.

Radio continuum surveys have, in the past, been of restricted use in cosmology. Most studies have concentrated on cross-correlations with the cosmic microwave background to detect the integrated Sachs-Wolfe effect, due to the large sky areas that can be surveyed. As we move into the SKA era, radio continuum surveys will have sufficient source density and sky area to play a major role in cosmology on the largest scales. In this chapter we summarise the experiments that can be carried out with the SKA as it is built up through the coming decade. We show that the SKA can play a unique role in constraining the non-Gaussianity parameter to \sigma(f_NL) ~ 1, and provide a unique handle on the systematics that inhibit weak lensing surveys. The SKA will also provide the necessary data to test the isotropy of the Universe at redshifts of order unity and thus evaluate the robustness of the cosmological principle.Thus, SKA continuum surveys will turn radio observations into a central probe of cosmological research in the coming decades.

Continuum and HI surveys with the Square Kilometre Array (SKA) will allow us to probe some of the most fundamental assumptions of modern cosmology, including the Cosmological Principle. SKA all-sky surveys will map an enormous slice of space-time and reveal cosmology at superhorizon scales and redshifts of order unity. We illustrate the potential of these surveys and discuss the prospects to measure the cosmic radio dipole at high fidelity. We outline several potentially transformational tests of cosmology to be carried out by means of SKA all-sky surveys.

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

Cosmic structures determine how light propagates through the Universe and consequently must be taken into account in the interpretation of observations. In the standard cosmological model at the largest scales, such structures are either ignored or treated as small perturbations to an isotropic and homogeneous Universe. This isotropic and homogeneous model is commonly assumed to emerge from some averaging process at the largest scales. We assume that there exists an averaging procedure that preserves the causal structure of space-time. Based on that assumption, we study the effects of averaging the geometry of space-time and derive an averaged version of the null geodesic equation of motion. For the averaged geometry we then assume a flat Friedmann-Lema\^{i}tre (FL) model and find that light propagation in this averaged FL model is not given by null geodesics of that model, but rather by a modified light propagation equation that contains an effective Hubble expansion rate, which differs from the Hubble rate of the averaged space-time.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

We investigate the contribution that a local over- or under-density can have on linear cosmic dipole estimations. We focus here on radio surveys, such as the NRAO VLA Sky Survey (NVSS), and forthcoming surveys such as those with the LOw Frequency ARray (LOFAR), the Australian Square Kilometre Array Pathfinder (ASKAP) and the Square Kilometre Array (SKA). The NVSS has already been used to estimate the cosmic radio dipole; it was shown recently that this radio dipole amplitude is larger than expected from a purely kinematic effect, assuming the velocity inferred from the dipole of the cosmic microwave background. We show here that a significant contribution to this excess could come from a local void or similar structure. In contrast to the kinetic contribution to the radio dipole, the structure dipole depends on the flux threshold of the survey and the wave band, which opens the chance to distinguish the two contributions.

Local measurements of the Hubble expansion rate are affected by structures like galaxy clusters or voids. Here we present a fully relativistic treatment of this effect, studying how clustering modifies the mean distance (modulus)-redshift relation and its dispersion in a standard $\Lambda$CDM universe. The best estimates of the local expansion rate stem from supernova observations at small redshifts (0.01<z<0.1). It is interesting to compare these local measurements with global fits to data from cosmic microwave background anisotropies. In particular, we argue that cosmic variance (i.e. the effects of the local structure) is of the same order of magnitude as the current observational errors and must be taken into account in local measurements of the Hubble expansion rate.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

We assume cold dark matter to possess a small bulk-viscous pressure which typically attenuates the growth of inhomogeneities. Explicit calculations, based on Eckart's theory of dissipative processes, reveal that for viscous cold dark matter the usual Newtonian approximation for perturbation scales smaller than the Hubble scale is no longer valid. We advocate the use of a neo-Newtonian approach which consistently incorporates pressure effects into the fluid dynamics and correctly reproduces the general relativistic dynamics. This result is of interest for numerical simulations of nonlinear structure formation involving nonstandard dark-matter fluids. We obtain upper limits on the magnitude of the viscous pressure by requiring that relevant perturbation amplitudes should grow sufficiently to enter the nonlinear stage.

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

We investigate the angular two-point correlation function of temperature in the WMAP maps. Updating and extending earlier results, we confirm the lack of correlations outside the Galaxy on angular scales greater than about 60 degrees at a level that would occur in 0.025 per cent of realizations of the concordance model. This represents a dramatic increase in significance from the original observations by the COBE-DMR and a marked increase in significance from the first-year WMAP maps. Given the rest of the reported angular power spectrum C_\ell, the lack of large-angle correlations that one infers outside the plane of the Galaxy requires covariance among the C_\ell up to \ell=5. Alternately, it requires both the unusually small (5 per cent of realizations) full-sky large-angle correlations, and an unusual coincidence of alignment of the Galaxy with the pattern of cosmological fluctuations (less than 2 per cent of those 5 per cent). We argue that unless there is some undiscovered systematic error in their collection or reduction, the data point towards a violation of statistical isotropy. The near-vanishing of the large-angle correlations in the cut-sky maps, together with their disagreement with results inferred from full-sky maps, remain open problems, and are very difficult to understand within the concordance model.

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

We use linear estimators to determine the magnitude and direction of the cosmic radio dipole from the NRAO VLA Sky Survey (NVSS) and the Westerbork Northern Sky Survey (WENSS). We show that special attention has to be given to the issues of bias due to shot noise, incomplete sky coverage and masking of the Milky Way. We compare several different estimators and show that conflicting claims in the literature can be attributed to the use of different estimators. We find that the NVSS and WENSS estimates of the cosmic radio dipole are consistent with each other and with the direction of the cosmic microwave background (CMB) dipole. We find from the NVSS a dipole amplitude of $(1.6 \pm 0.6) \times 10^{-2}$ in direction $(\mathrm{RA}, \mathrm{dec})= (154^\circ \pm 21^\circ, -2^\circ \pm 21^\circ)$. This amplitude exceeds the one expected from the CMB by a factor of about 3 and is inconsistent with the assumption of a pure kinetic origin of the radio dipole at 99.5% CL.

The anomalous lack of large angle temperature correlations has been a surprising feature of the cosmic microwave background (CMB) since first observed by COBE-DMR and subsequently confirmed and strengthened by the Wilkinson Microwave Anisotropy Probe. This anomaly may point to the need for modifications of the standard model of cosmology or may indicate that our Universe is a rare statistical fluctuation within that model. Further observations of the temperature auto-correlation function will not elucidate the issue; sufficiently high precision statistical observations already exist. Instead, alternative probes are required. In this work we explore the expectations for forthcoming polarization observations. We define a prescription to test the hypothesis that the large-angle CMB temperature perturbations in our Universe represent a rare statistical fluctuation within the standard cosmological model. These tests are based on the temperature-$Q$ Stokes parameter correlation. Unfortunately these tests cannot be expected to be definitive. However, we do show that if this $TQ$-correlation is observed to be sufficiently large over an appropriately chosen angular range, then the hypothesis can be rejected at a high confidence level. We quantify these statements and optimize the statistics we have constructed to apply to the anticipated polarization data. We find that we can construct a statistic that has a 25 per cent chance of excluding the hypothesis that we live in a rare realization of LCDM at the 99.9 per cent confidence level.

We test the isotropy of the expansion of the Universe by estimating the hemispherical anisotropy of supernova type Ia (SN Ia) Hubble diagrams at low redshifts (z<0.2). We compare the best fit Hubble diagrams in pairs of hemispheres and search for the maximal asymmetric orientation. For an isotropic Universe, we expect only a small asymmetry due to noise and the presence of nearby structures. This test does not depend on the assumed content of the Universe, the assumed model of gravity, or the spatial curvature of the Universe. The expectation for possible fluctuations due to large scale structure is evaluated for the \Lambda cold dark matter (\Lambda CDM) model and is compared to the supernova data from the Constitution set for four different light curve fitters, thus allowing a study of the systematic effects. The expected order of magnitude of the hemispherical asymmetry of the Hubble expansion agrees with the observed one. The direction of the Hubble asymmetry is established at 95% confidence level (C.L.) using both, the MLCS2k2 and the SALT II light curve fitter. The highest expansion rate is found towards (l, b) ~ (-35{\deg},-19{\deg}), which agrees with directions reported by other studies. Its amplitude is not in contradiction to expectations from the \Lambda CDM model. The measured Hubble anisotropy is \Delta H/H ~ 0.026. With 95% C.L. the expansion asymmetry is \Delta H/H<0.038.

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

Recent observations of the cosmic microwave background (CMB) at smallest angular scales and updated abundances of primordial elements, indicate an increase of the energy density and the helium-4 abundance with respect to standard big bang nucleosynthesis with three neutrino flavour. This calls for a reanalysis of the observational bounds on neutrino chemical potentials, which encode the number asymmetry between cosmic neutrinos and anti-neutrinos and thus measures the lepton asymmetry of the Universe. We compare recent data with a big bang nucleosynthesis code, assuming neutrino flavour equilibration via neutrino oscillations before the onset of big bang nucleosynthesis. We find a slight preference for negative neutrino chemical potentials, which would imply an excess of anti-neutrinos and thus a negative lepton number of the Universe. This lepton asymmetry could exceed the baryon asymmetry by orders of magnitude.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

Fluids often display dissipative properties. We explore dissipation in the form of bulk viscosity in the cold dark matter fluid. We constrain this model using current data from supernovae, baryon acoustic oscillations and the cosmic microwave background. Considering the isotropic and homogeneous background only, viscous dark matter is allowed to have a bulk viscosity $\lesssim 10^7$ Pa$\cdot$s, also consistent with the expected integrated Sachs-Wolfe effect (which plagues some models with bulk viscosity). We further investigate the small-scale formation of viscous dark matter halos, which turns out to place significantly stronger constraints on the dark matter viscosity. The existence of dwarf galaxies is guaranteed only for much smaller values of the dark matter viscosity, $\lesssim 10^{-3}$ Pa$\cdot$s.

We show how the relativistic matter and velocity power spectra behave in different gauges. We construct a new gauge where both spectra coincide with Newtonian theory on all scales. However, in this gauge there are geometric quantities present which do not exist in Newtonian theory, for example the local variation of the Hubble parameter. Comparing this quantity to second order Newtonian quantities, we find that Newtonian theory is inaccurate on scales larger than 10 Mpc. This stresses the importance of relativistic corrections to Newtonian cosmological N-body simulations on these scales.

We find the best-fit cosmological parameters for a scenario of inflation with only the sufficient amount of accelerated expansion for the $\lambda\phi^4$ potential. While for the simplest scenario of chaotic inflation all observable primordial fluctuations cross the Hubble horizon during the slow-roll epoch, for the scenario of just-enough inflation the slow-roll conditions are violated at the largest length scales. Performing a numerical mode-by-mode integration for the perturbations on the largest scales and comparing the predicted anisotropies of the cosmic microwave background to results from the WMAP 7-yr data analysis, we find the initial conditions in agreement with current cosmological data. In contrast to the simplest chaotic model for the quartic potential, the just-enough inflation scenario is not ruled out. Although this scenario naturally gives rise to a modification of the first multipoles, for a quartic potential it cannot explain the lack of power at the largest angular scales.

We investigate how large lepton asymmetries affect the evolution of the early universe at times before big bang nucleosynthesis and in particular how they influence the relic density of WIMP dark matter. In comparison to the standard calculation of the relic WIMP abundance we find a decrease, depending on the lepton flavour asymmetry. We find an effect of up to 20 per cent for lepton flavour asymmetries $l_f= {\cal O}(0.1)$.

A poorly constrained parameter in the Standard Model of Cosmology is the lepton asymmetry l = \sum_f l_f=\sum_f(n_f+n_{\nu_f})/s. Each flavour asymmetry l_f with f=e, \mu, {\tau} is the sum of the net particle density of the charged leptons n_f and their corresponding neutrinos, normalized with the entropy density s. Constraints on l_f \leq O(0.1) from BBN and CMB allow for lepton flavour asymmetries orders of magnitudes larger then the baryon asymmetry b ~ 10^{-10}. In this article we show how such large lepton (flavour) asymmetries influence the early universe, in particular the freeze out of WIMPs and the cosmic QCD transition.

Modern cosmology relies on the assumption of large-scale isotropy and homogeneity of the Universe. However, locally the Universe is inhomogeneous and anisotropic. So, how can local measurements (at the 100 Mpc scale) be used to determine global cosmological parameters (defined at the 10 Gpc scale)? We use Buchert's averaging formalism and determine a set of locally averaged cosmological parameters in the context of the flat Lambda cold dark matter model. We calculate their ensemble means (i.e. their global values) and variances (i.e. their cosmic variances). We apply our results to typical survey geometries and focus on the study of the effects of local fluctuations of the curvature parameter. By this means we show, that in the linear regime cosmological backreaction and averaging can be reformulated as the issue of cosmic variance. The cosmic variance is found largest for the curvature parameter and discuss some of its consequences. We further propose to use the observed variance of cosmological parameters to measure the growth factor. [abbreviated]

We discuss selected large-scale anomalies in the maps of temperature anisotropies in the cosmic microwave background. Specfically, these include alignments of the largest modes of CMB anisotropy with one another and with the geometry and direction of motion of the Solar System, and the unexpected absence of two-point angular corellations especially outside the region of the sky most contaminated by the Galaxy. We discuss these findings in relation to expectations from standard inflationary cosmology. This paper is adapted from a talk given by one of us (GDS) at the SEENET-2011 meeting in August 2011 on the Serbian bank of the Danube River.

Modern cosmology suggests that the Universe contains two dark components -- dark matter and dark energy -- both unkown in laboratory physics and both lacking direct evidence. Alternatively, a unified dark sector, described by a single fluid, has been proposed. Dissipation is a common phenomenon in nature and it thus seems natural to consider models dominated by a viscous dark fluid. We focus on the study of bulk viscosity, as isotropy and homogeneity at large scales implies the suppression of shear viscosity, heat flow and diffusion. The generic ansatz $\xi \propto \rho^{\nu}$ for the coefficient of bulk viscosity ($\rho$ denotes the mass/energy density), which for $\nu = - 1/2$ mimics the $\Lambda$CDM background evolution, offers excellent fits to supernova and H(z) data. We show that viscous dark fluids suffer from large contributions to the integrated Sachs-Wolfe effect (generalising a previous study by Li & Barrow) and a suppression of structure growth at small-scales (as seen from a generalized Meszaros equation). Based on recent observations, we conclude that viscous dark fluid models (with $\xi \propto \rho^{\nu}$ and neglecting baryons) are strongly challenged.

Measuring the deviation from scale invariance of the primordial power spectrum is a critical test of inflation. In this paper we reconstruct the shape of the primordial power spectrum of curvature perturbations from the cosmic microwave background data, including the 7-year Wilkinson Microwave Anisotropy Probe data and the Atacama Cosmology Telescope 148 GHz data, by using a binning method of a cubic spline interpolation in log-log space. We find that the power-law spectrum is preferred by the data and that the Harrison-Zel'dovich spectrum is disfavored at 95% confidence level. These conclusions hold with and without allowing for tensor modes, however the simpler model without tensors is preferred by the data. We do not find evidence for a feature in the primordial power spectrum - in full agreement with generic predictions from cosmological inflation.

We present forecasts for constraints on cosmological models which can be obtained by forthcoming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the planned surveys to predict measurements obtained with the source auto-correlation, the cross-correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the cross-correlation of radio sources with foreground objects due to cosmic magnification, and a joint analysis together with the CMB power spectrum and supernovae. We show that near future radio surveys will bring complementary measurements to other experiments, probing different cosmological volumes, and having different systematics. Our results show that the unprecedented sky coverage of these surveys combined should provide the most significant measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the ISW effect will significantly tighten constraints on modified gravity parameters, while the best measurements of dark energy models will come from galaxy auto-correlation function analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will be possible to measure the dark energy parameters with an uncertainty of \{$\sigma (w_0) = 0.05$, $\sigma (w_a) = 0.12$\} and the modified gravity parameters \{$\sigma (\eta_0) = 0.10$, $\sigma (\mu_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors. Finally, we show that radio surveys would detect a primordial non-Gaussianity of $f_{\rm NL}$ = 8 at 1-$\sigma$ and we briefly discuss other promising probes.

The large-angle, low multipole cosmic microwave background (CMB) provides a unique view of the largest angular scales in the Universe. Study of these scales is hampered by the facts that we have only one Universe to observe, only a few independent samples of the underlying statistical distribution of these modes, and an incomplete sky to observe due to the interposing Galaxy. Techniques for reconstructing a full sky from partial sky data are well known and have been applied to the large angular scales. In this work we critically study the reconstruction process and show that, in practise, the reconstruction is biased due to leakage of information from the region obscured by foregrounds to the region used for the reconstruction. We conclude that, despite being suboptimal in a technical sense, using the unobscured region without reconstructing is the most robust measure of the true CMB sky. We also show that for noise free data reconstructing using the usual optimal, unbiased estimator may be employed without smoothing thus avoiding the leakage problem. Unfortunately, directly applying this to real data with noise and residual, unmasked foregrounds yields highly biased reconstructions requiring further care to apply this method successfully to real-world CMB.

Determining the energy scale of inflation is crucial to understand the nature of inflation in the early Universe. We place observational constraints on the energy scale of the observable part of the inflaton potential by combining the 7-year Wilkinson Microwave Anisotropy Probe data with distance measurements from the baryon acoustic oscillations in the distribution of galaxies and the Hubble constant measurement. Our analysis provides an upper limit on this energy scale, 2.3 \times 10^{16} GeV at 95% confidence level. Moreover, we forecast the sensitivity and constraints achievable by the Planck experiment by performing Monte Carlo studies on simulated data. Planck could significantly improve the constraints on the energy scale of inflation and on the shape of the inflaton potential.

We show non-perturbatively that the power spectrum of a self-interacting scalar field in de Sitter space-time is strongly suppressed on large scales. The cut-off scale depends on the strength of the self-coupling, the number of e-folds of quasi-de Sitter evolution, and its expansion rate. As a consequence, the two-point correlation function of field fluctuations is free from infra-red divergencies.

We present a short (and necessarily incomplete) review of the evidence for the accelerated expansion of the Universe. The most direct probe of acceleration relies on the detailed study of supernovae (SN) of type Ia. Assuming that these are standardizable candles and that they fairly sample a homogeneous and isotropic Universe, the evidence for acceleration can be tested in a model- and calibration-independent way. Various light-curve fitting procedures have been proposed and tested. While several fitters give consistent results for the so-called Constitution set, they lead to inconsistent results for the recently released SDSS SN. Adopting the SALT fitter and relying on the Union set, cosmic acceleration is detected by a purely kinematic test at 7 sigma when spatial flatness is assumed and at 4 sigma without assumption on the spatial geometry. A weak point of the described method is the local set of SN (at z < 0.2), as these SN are essential to anchor the Hubble diagram. These SN are drawn from a volume much smaller than the Hubble volume and could be affected by local structure. Without the assumption of homogeneity, there is no evidence for acceleration, as the effects of acceleration are degenerate with the effects of inhomogeneities. Unless we sit in the centre of the Universe, such inhomogeneities can be constrained by SN observations by means of tests of the isotropy of the Hubble flow.

We review the recently found large-scale anomalies in the maps of temperature anisotropies in the cosmic microwave background. These include alignments of the largest modes of CMB anisotropy with each other and with geometry and direction of motion of the Solar System, and the unusually low power at these largest scales. We discuss these findings in relation to expectation from standard inflationary cosmology, their statistical significance, the tools to study them, and the various attempts to explain them.

We study the influence of lepton asymmetry on the evolution of the early Universe. The lepton asymmetry $l$ is poorly constrained by observations and might be orders of magnitude larger than the baryon asymmetry $b$, $|l|/b \leq 2\times 10^8$. We find that lepton asymmetries that are large compared to the tiny baryon asymmetry, can influence the dynamics of the QCD phase transition significantly. The cosmic trajectory in the $\mu_B-T$ phase diagram of strongly interacting matter becomes a function of lepton (flavour) asymmetry. Large lepton asymmetry could lead to a cosmic QCD phase transition of first order.

The Wilkinson Microwave Anisotropy Probe (WMAP) experiment has detected reionization at the $5.5 \sigma$ level and has reported a mean optical depth of $0.088 \pm 0.015$. A powerful probe of reionization is the large-angle $EE$ polarization power spectrum, which is now (since the first five years of data from WMAP) cosmic variance limited for $2\le l \le6$. Here we consider partial reionization caused by WIMP dark matter annihilation, and calculate the expected polarization power spectrum. We compare the dark matter models with a standard 2-step reionization theory, and examine whether the models may be distinguished using current, and future CMB observations. We consider dark matter annihilation at intermediate redshifts ($z<60$) due to halos, as well as annihilation at higher redshifts due to free particles. In order to study the effect of high redshift dark matter annihilation on CMB power spectra, it is essential to include the contribution of residual electrons (left over from recombination) to the ionization history. Dark matter halos at redshifts $z<60$ influence the low multipoles $l<20$ in the $EE$ power spectrum, while the annihilation of free particle dark matter at high redshifts $z>100$ mainly affects multipoles $l>10$.

We consider slow-roll inflation for a single scalar field with an arbitrary potential and an arbitrary nonminimal coupling to the Gauss-Bonnet term. By introducing a combined hierarchy of Hubble and Gauss-Bonnet flow functions, we analytically derive the power spectra of scalar and tensor perturbations. The standard consistency relation between the tensor-to-scalar ratio and the spectral index of tensor perturbations is broken. We apply this formalism to a specific model with a monomial potential and an inverse monomial Gauss-Bonnet coupling and constrain it by the 7-year Wilkinson Microwave Anisotropy Probe data. The Gauss-Bonnet term with a positive (or negative) coupling may lead to a reduction (or enhancement) of the tensor-to-scalar ratio and hence may revive the quartic potential ruled out by recent cosmological data.

On large angular scales (greater than about 60 degrees), the two-point angular correlation function of the temperature of the cosmic microwave background (CMB), as measured (outside of the plane of the Galaxy) by the Wilkinson Microwave Anisotropy Probe, shows significantly lower large-angle correlations than expected from the standard inflationary cosmological model. Furthermore, when derived from the full CMB sky, the two lowest cosmologically interesting multipoles, the quadrupole (l=2) and the octopole (l=3), are unexpectedly aligned with each other. Using randomly generated full-sky and cut-sky maps, we investigate whether these anomalies are correlated at a statistically significant level. We conclusively demonstrate that, assuming Gaussian random and statistically isotropic CMB anisotropies, there is no statistically significant correlation between the missing power on large angular scales in the CMB and the alignment of the l=2 and l=3 multipoles. The chance to measure the sky with both such a lack of large-angle correlation and such an alignment of the low multipoles is thus quantified to be below 10^{-6}.

This work summarises some of the attempts to explain the phenomenon of dark energy as an effective description of complex gravitational physics and the proper interpretation of observations. Cosmological backreaction has been shown to be relevant for observational (precision) cosmology, nevertheless no convincing explanation of dark energy by means of backreaction has been given so far.

As the Einstein equations are non-linear, spatial averaging and temporal evolution do not commute. Therefore, the evolution of the averaged universe is affected by inhomogeneities. It is, however, highly controversial how large these cosmological backreaction effects are. We use the supernova data of the Constitution set up to a redshift of 0.1 in order to analyse to what extent the measurement of the Hubble constant is affected. The size of the effect depends on the size of the volume that is averaged over. The observational results are then compared to the theory of the backreaction mechanism.

We consider power-law inflation with a Gauss-Bonnet correction inspired by string theory. We analyze the stability of cosmological perturbations and obtain the allowed parameter space. We find that for GB-dominated inflation ultra-violet instabilities of either scalar or tensor perturbations show up on small scales. The Gauss-Bonnet correction with a positive (or negative) coupling may lead to a reduction (or enhancement) of the tensor-to-scalar ratio in the potential-dominated case. We place tight constraints on the model parameters by making use of the WMAP 5-year data.

We present counter examples to the claim that the $\lambda \phi^4$ inflaton potential is excluded by recent cosmological data. Finding counter examples requires that the actually observed primordial fluctuations are generated at the onset of the slow-roll regime of inflation. This set up for the initial conditions is therefore different from the usual scenario of chaotic inflation where inflation starts long before the observed fluctuations are created. The primordial power spectrum of "just enough" chaotic inflation violates scale-invariance in a way consistent with observations.

If dark matter is made up of Weakly Interacting Massive Particles, the annihilation of these particles in halos results in energy being released, some of which is absorbed by gas, causing partial ionization and heating. Dark matter annihilation may result in partial ionization and gas heating at high redshifts, even before the formation of the first stars. It is shown that early ionization results in a transfer of power to higher multipoles in the large angle CMB polarization power spectra. Future CMB experiments may be able to place constraints on certain light dark matter models. We also investigate the effect of gas heating on the expected H21 cm power spectrum. Heating by particle annihilation results in a decrease in the amplitude of the H21 cm power spectrum as the gas temperature $T$ becomes comparable to the CMB temperature $T_\gamma$, and then an increase as $T > T_\gamma$. The result is a minimum in the power spectrum at the redshift for which $T \approx T_\gamma$. Only certain models (low particle masses $\sim$ 10 GeV, or favorable halo parameters) show this effect. Within these models, observations of the H21 cm power spectrum at multiple redshifts can help us obtain constraints on dark matter particle and halo properties.

The cosmological principle says that the Universe is spatially homogeneous and isotropic. It predicts, among other phenomena, the cosmic redshift of light and the Hubble law. Nevertheless, the existence of structure in the Universe violates the (exact) cosmological principle. A more precise formulation of the cosmological principle must allow for the formation of structure and must therefore incorporate probability distributions. In this contribution to the Memorial Volume for Wolfgang Kummer, a great teacher and mentor to me, I discuss how we could formulate a new version of the cosmological principle, how to test it, and how to possibly justify it by fundamental physics. My contribution starts with some of my memories of Wolfgang.

We present a calibration-independent test of the accelerated expansion of the universe using supernova type Ia data. The test is also model-independent in the sense that no assumptions about the content of the universe or about the parameterization of the deceleration parameter are made and that it does not assume any dynamical equations of motion. Yet, the test assumes the universe and the distribution of supernovae to be statistically homogeneous and isotropic. A significant reduction of systematic effects, as compared to our previous, calibration-dependent test, is achieved. Accelerated expansion is detected at significant level (4.3 sigma in the 2007 Gold sample, 7.2 sigma in the 2008 Union sample) if the universe is spatially flat. This result depends, however, crucially on supernovae with a redshift smaller than 0.1, for which the assumption of statistical isotropy and homogeneity is less well established.

The annihilation of dark matter particles releases energy, ionizing some of the gas in the Universe. We investigate the effect of dark matter halos on reionization. We show that the effect depends on the assumed density profile, the particle mass, and the assumed minimum halo mass. For NFW halos and typical WIMPs, we find the effect to be quite small. However, light dark matter candidates in the MeV range can contribute significantly to reionization and can make an important contribution to the measured optical depth. This effect may be used to constrain light dark matter models. We also study the effect of varying the halo density profile on reionization.

Many cosmological observations call for the existence of dark matter. The most direct evidence for dark matter is inferred from the measured flatness of galactic rotation curves. The latter is based on Newtonian gravity. Alternative approaches to the rotation curve problem by means of general relativity have recently been put forward. The class of models of interest is a subset of the axially symmetric and stationary solutions of Einstein's equations with rotating dust. As a step toward the understanding of general relativistic galaxy models, we analyse rigidly as well as non-rigidly rotating (Post-)Newtonian spacetimes. We find that the Newtonian limit of the considered general relativistic galaxy model leads to Post-Newtonian terms in the metric.

The missing of a Keplerian fall-off in the observed galaxy rotation curves represents classical evidence for the existence of dark matter on galactic scales. There has been some recent activity concerning the potential of modelling galactic systems with the help of general relativity. This was motivated by claims that by the use of full general relativity dark matter could be made superfluous. Here we focus on possible axisymmetric and stationary solutions of Einstein's equations with rotating dust. After a short review of the current debate we pursue the idea of approaching such relativistic models in a Newtonian language. We analyse rigidly as well as differentially rotating Newtonian and Post-Newtonian spacetimes and find that it is necessary to incorporate a Post-Newtonian term in order to make physical sense.

We present a replica field-theoretic approach to stochastic inflation in which a manifestation of dimensional reduction is found. The scale above which the latter dominates grows exponentially fast with time and thus affects largest super-horizon scales. We find inevitable modifications of the spectral index on those scales. An explicit relation between the noise correlator and the non-gaussianity parameter f_NL is found.

We adopt methods from statistical field theory to stochastic inflation. For the example of a free test field in de Sitter and power-law inflation, the power spectrum of long-wavelength fluctuations is computed. We study its dependence on the shape of the filter that separates long from short wavelength modes. While for filters with infinite support the phenomenon of dimensional reductions is found on large super-horizon scales, filters with compact support return a scale-invariant power spectrum in the infra-red. Features of the power spectrum, induced by the filter, decay within a few e-foldings. Thus the late-time power spectrum is independent of the filter details.

Due to the non-commutation of spatial averaging and temporal evolution, inhomogeneities and anisotropies (cosmic structures) influence the evolution of the averaged Universe via the cosmological backreaction mechanism. We study the backreaction effect as a function of averaging scale in a perturbative approach up to higher orders. We calculate the hierarchy of the critical scales, at which 10% effects show up from averaging at different orders. The dominant contribution comes from the averaged spatial curvature, observable up to scales of 200 Mpc. The cosmic variance of the local Hubble rate is 10% (5%) for spherical regions of radius 40 (60) Mpc. We compare our result to the one from Newtonian cosmology and Hubble Space Telescope Key Project data.

The present standard model of cosmology states that the known particles carry only a tiny fraction of total mass and energy of the Universe. Rather, unknown dark matter and dark energy are the dominant contributions to the cosmic energy budget. We review the logic that leads to the postulated dark energy and present an alternative point of view, in which the puzzle may be solved by properly taking into account the influence of cosmic structures on global observables. We illustrate the effect of averaging on the measurement of the Hubble constant.

We test the present expansion of the universe using supernova type Ia data without making any assumptions about the matter and energy content of the universe or about the parameterization of the deceleration parameter. We assume the cosmological principle to apply in a strict sense. The result strongly depends on the data set, the light-curve fitting method and the calibration of the absolute magnitude used for the test, indicating strong systematic errors. Nevertheless, in a spatially flat universe there is at least a 5 sigma evidence for acceleration which drops to 1.8 sigma in an open universe.

Cosmological backreaction has been suggested as an explanation of dark energy and is heavily disputed since. We combine cosmological perturbation theory with Buchert's non-perturbative framework, calculate the relevant averaged observables up to second order in the comoving synchronous gauge, and discuss their gauge dependence. With the help of an integrability condition the leading second order contributions follow from the first order calculation. We focus on the onset of cosmological backreaction, as a perturbative treatment is necessarily restricted to the era when the effect is still small. We demonstrate that the leading contributions to all averaged physical observables are completely specified on the boundary of the averaged domain. For any finite domain these surface terms are nonzero in general and thus backreaction is for real. We map the backreaction effect on an effectively homogeneous and isotropic fluid. The generic effective equation of state is not only time dependent, but also depends on scale.

We test the isotropy of the Hubble diagram. At small redshifts, this is possible without assumptions on the cosmic inventory and provides a fundamental test of the cosmological principle. At higher redshift we check for the self-consistency of the LambdaCDM model. At small redshifts, we use public supernovae (SNe) Ia data to determine the deceleration parameter q_0 and the SN calibration on opposite hemispheres. For the complete data sets we fit Omega_M and the SN calibration on opposite hemispheres. A statistically significant anisotropy of the Hubble diagram at redshifts z < 0.2 is discovered (> 95% C.L.). While data from the North Galactic hemisphere favour the accelerated expansion of the Universe, data from the South Galactic hemisphere are not conclusive. The hemispheric asymmetry is maximal toward a direction close to the equatorial poles. The discrepancy between the equatorial North and South hemispheres shows up in the SN calibration. For the LambdaCDM model fitted to all available SNe, we find the same asymmetry. The alignment of discrepancies between hemispheric Hubble diagrams with the equatorial frame seems to point toward a systematic error in the SN search, observation, analysis or data reduction. We also find that our model independent test cannot exclude the case of the deceleration of the expansion at a statistically significant level.

At the largest angular scales the presence of a number of unexpected features has been confirmed by the latest measurements of the cosmic microwave background (CMB). Among them are the anomalous alignment of the quadrupole and octopole with each other as well as the stubborn lack of angular correlation on scales >60deg. We search for correlations between these two phenomena and demonstrate their absence. A Monte Carlo likelihood analysis confirms previous studies and shows that the joint likelihood of both anomalies is incompatible with the best-fit Lambda Cold Dark Matter model at >99.95% C.L. At the same time, a presumed special axis (the `Axis of Evil') identified on the microwave sky demands additional contributions to multipole power on top of the primordial standard inflationary ones. We find that the notion of a preferred axis in the CMB is misleading and inconsistent with three-year data from the Wilkinson Microwave Anisotropy Probe (WMAP). Rather the data require a preferred plane, whereupon the axis is just the normal direction to that plane. Rotational symmetry within that plane is inconsistent with the observations and is ruled out at high confidence.

We address the effect of an extended local foreground on the low-l anomalies found in the CMB. Recent X-ray catalogues point us to the existence of very massive superstructures at the 100 h^(-1) Mpc scale that contribute significantly to the dipole velocity profile. Being highly non-linear, these structures provide us a natural candidate to leave an imprint on the CMB sky via a local Rees-Sciama effect. We show that the Rees-Sciama effect of local foregrounds can induce CMB anisotropy of DeltaT/T ~ 10^(-5) and we analyse its impact on multipole power as well as the induced phase pattern on largest angular scales.

The size and time of formation of the first gravitationally bound objects in the Universe is set by the microphysical properties of the dark matter. It is argued that observations seem to favour cold and thermal candidates for the main contribution to the dark matter. For that type of dark matter, the size and time of formation of the first halos is determined by the elastic cross sections and mass of the CDM particles. Consequently, the astrophysics of CDM might allow us to measure some of the fundamental parameters of CDM particles. Essential for observations is the survival rate and spatial distribution of the very first objetcs, which are currently under debate.

The evolution of the Universe is the ultimate laboratory to study fundamental physics across energy scales that span about 25 orders of magnitude: from the grand unification scale through particle and nuclear physics scales down to the scale of atomic physics. The standard models of cosmology and particle physics provide the basic understanding of the early and present Universe and predict a series of phase transitions that occurred in succession during the expansion and cooling history of the Universe. We survey these phase transitions, highlighting the equilibrium and non-equilibrium effects as well as their observational and cosmological consequences. We discuss the current theoretical and experimental programs to study phase transitions in QCD and nuclear matter in accelerators along with the new results on novel states of matter as well as on multi- fragmentation in nuclear matter. A critical assessment of similarities and differences between the conditions in the early universe and those in ultra- relativistic heavy ion collisions is presented. Cosmological observations and accelerator experiments are converging towards an unprecedented understanding of the early and present Universe.

The large-angle (low-ell) correlations of the Cosmic Microwave Background (CMB) as reported by the Wilkinson Microwave Anisotropy Probe (WMAP) after their first year of observations exhibited statistically significant anomalies compared to the predictions of the standard inflationary big-bang model. We suggested then that these implied the presence of a solar system foreground, a systematic correlated with solar system geometry, or both. We re-examine these anomalies for the data from the first three years of WMAP's operation. We show that, despite the identification by the WMAP team of a systematic correlated with the equinoxes and the ecliptic, the anomalies in the first-year Internal Linear Combination (ILC) map persist in the three-year ILC map, in all-but-one case at similar statistical significance. The three-year ILC quadrupole and octopole therefore remain inconsistent with statistical isotropy -- they are correlated with each other (99.6% C.L.), and there are statistically significant correlations with local geometry, especially that of the solar system. The angular two-point correlation function at scales >60 degrees in the regions outside the (kp0) galactic cut, where it is most reliably determined, is approximately zero in all wavebands and is even more discrepant with the best fit LambdaCDM inflationary model than in the first-year data - 99.97% C.L. for the new ILC map. The full-sky ILC map, on the other hand, has a non-vanishing angular two-point correlation function, apparently driven by the region inside the cut, but which does not agree better with LambdaCDM. The role of the newly identified low-ell systematics is more puzzling than reassuring.

The microwave sky shows unexpected features at the largest angular scales, among them the alignments of the dipole, quadrupole and octopole. Motivated by recent X-ray cluster studies, we investigate the possibility that local structures at the 100 h^(-1) Mpc scale could be responsible for such correlations. These structures give rise to a local Rees-Sciama contribution to the microwave sky that may amount to Delta T/T ~ 10^(-5) at the largest angular scales. We model local structures by a spherical overdensity (Lemaitre-Tolman-Bondi model) and assume that the Local Group is falling toward the centre. We superimpose the local Rees-Sciama effect on a statistically isotropic, gaussian sky. As expected we find alignments among low multipoles, but a closer look reveals that they do not agree with the type of correlations revealed by the data.

The dark matter distribution on small scales may depend on the the properties of the first generation of dark matter halos to form, which is in turn determined by the microphysics of the dark matter particles. We overview the microphysics of WIMPs and our calculations of the collisional damping and free streaming scales. We then plot the resulting density perturbation power spectrum and red-shift at which typical halos form, taking into account the effect of uncertainties in the WIMP properties (mass and interaction channel) and the primordial power spectrum. Finally we review recent developments regarding the properties and fate of the first WIMPy halos.

[Abridged] We apply the multipole vector framework to full-sky maps derived from the first year WMAP data. We significantly extend our earlier work showing that the two lowest cosmologically interesting multipoles, l=2 and 3, are not statistically isotropic. These results are compared to the findings obtained using related methods. In particular, the planes of the quadrupole and the octopole are unexpectedly aligned. Moreover, the combined quadrupole plus octopole is surprisingly aligned with the geometry and direction of motion of the solar system: the plane they define is perpendicular to the ecliptic plane and to the plane defined by the dipole direction, and the ecliptic plane carefully separates stronger from weaker extrema, running within a couple of degrees of the null-contour between a maximum and a minimum over more than 120deg of the sky. Even given the alignment of the quadrupole and octopole with each other, we find that their alignment with the ecliptic is unlikely at >98% C.L., and argue that it is in fact unlikely at >99.9% C.L. We explore the role of foregrounds showing that the known Galactic foregrounds are unlikely to lead to these correlations. Multipole vectors, like individual a_lm, are very sensitive to sky cuts, and we demonstrate that analyses using cut skies induce relatively large errors, thus weakening the observed correlations but preserving their consistency with the full-sky results. Finally we apply our tests to COBE cut-sky maps and briefly extend the analysis to higher multipoles. If the correlations we observe are indeed a signal of non-cosmic origin, then the lack of low-l power will very likely be exacerbated, with important consequences for our understanding of cosmology on large scales.

Dark matter direct and indirect detection signals depend crucially on the dark matter distribution. While the formation of large scale structure is independent of the nature of the cold dark matter (CDM), the fate of inhomogeneities on sub-galactic scales, and hence the present day CDM distribution on these scales, depends on the micro-physics of the CDM particles. We study the density contrast of Weakly Interacting Massive Particles (WIMPs) on sub-galactic scales. We calculate the damping of the primordial power spectrum due to collisional damping and free-streaming of WIMPy CDM and show that free-streaming leads to a CDM power spectrum with a sharp cut-off at about $10^{-6} M_\odot$. We also calculate the transfer function for the growth of the inhomogeneities in the linear regime, taking into account the suppression in the growth of the CDM density contrast after matter-radiation equality due to baryons and show that our analytic results are in good agreement with numerical calculations. Combining the transfer function with the damping of the primordial fluctuations we produce a WMAP normalized primordial CDM power spectrum, which can serve as an input for high resolution CDM simulations. We find that the smallest inhomogeneities typically have co-moving radius of about 1 pc and enter the non-linear regime at a redshift of $60 \pm 20$. We study the effect of scale dependence of the primordial power spectrum on these numbers and also use the spherical collapse model to make simple estimates of the properties of the first generation of WIMP halos to form. We find that the very first WIMPy halos may have a significant impact on indirect dark matter searches.

The large-angle (low-l) correlations of the Cosmic Microwave Background exhibit several statistically significant anomalies compared to the standard inflationary big-bang model, however no connection has hitherto been drawn between them. Here we show that the quadrupole and octopole are far more correlated (99.97% C.L.) than previously thought. The quadrupole plane and the three octopole planes are remarkably aligned. Three of these planes are orthogonal to the ecliptic at a level inconsistent with gaussian random statistically isotropic skies at 99.8% C.L., and the normals to these planes are aligned at 99.9% C.L. with the direction of the cosmological dipole and with the equinoxes. The remaining octopole plane is orthogonal to the supergalactic plane at >99.9% C.L. In a combined quadrupole-octopole map, the ecliptic plane narrowly threads between a hot spot and a cold spot over approximately 1/3 of the sky, and separates the three strongest extrema (in the south ecliptic hemisphere) from the three weakest extrema (in the north ecliptic hemisphere).

The formation of large scale structure is independent of the nature of the cold dark matter (CDM), however the fate of very small scale inhomogeneities depends on the micro-physics of the CDM particles. We investigate the matter power spectrum for scales that enter the Hubble radius well before matter-radiation equality, and follow its evolution until the time when the first inhomogeneities become non-linear. Our focus lies on weakly interacting massive particles (WIMPs), and as a concrete example we analyze the case when the lightest supersymmetric particle is a bino. We show that collisional damping and free-streaming of WIMPs lead to a matter power spectrum with a sharp cut-off at about 10^(-6) M_odot and a maximum close to that cut-off. We also calculate the transfer function for the growth of the inhomogeneities in the linear regime. These three effects (collisional damping, free-streaming and gravitational growth) are combined to provide a WMAP normalized primordial CDM power spectrum, which could serve as an input for high resolution CDM simulations. The smallest inhomogeneities typically enter the non-linear regime at a redshift of about 60.

The observational constraints on the primordial power spectrum have tightened considerably with the release of the first year analysis of the WMAP observations, especially when combined with the results from other CMB experiments and galaxy redshift surveys. These observations allow us to constrain the physics of cosmological inflation: (i) The data show that the Hubble distance is almost constant during inflation. While observable modes cross the Hubble scale, it changes by less than 3% during one e-folding: d(d_H)/dt < 0.032 at 2 sigma. The distance scale of inflation itself remains poorly constrained: 1.2 x 10^{-28} cm < d_H < 1 cm. (ii) We present a new classification of single-field inflationary scenarios (including scenarios beyond slow-roll inflation), based on physical criteria, namely the behaviour of the kinetic and total energy densities of the inflaton field. The current data show no preference for any of the scenarios. (iii) For the first time the slow-roll assumption could be dropped from the data analysis and replaced by the more general assumption that the Hubble scale is (almost) constant during the observable part of inflation. We present simple analytic expressions for the scalar and tensor power spectra for this very general class of inflation models and test their accuracy.

A curvature self-interaction of the cosmic gas is shown to mimic a cosmological constant or other forms of dark energy, such as a rolling tachyon condensate or a Chaplygin gas. Any given Hubble rate and deceleration parameter can be traced back to the action of an effective curvature force on the gas particles. This force self-consistently reacts back on the cosmological dynamics. The links between an imperfect fluid description, a kinetic description with effective antifriction forces, and curvature forces, which represent a non-minimal coupling of gravity to matter, are established.

The history of the Universe after its first second is now tested by high quality observations of light element abundances and temperature anisotropies of the cosmic microwave background. The epoch of the first second itself has not been tested directly yet; however, it is constrained by experiments at particle and heavy ion accelerators. Here I attempt to describe the epoch between the electroweak transition and the primordial nucleosynthesis. The most dramatic event in that era is the quark--hadron transition at 10 $\mu$s. Quarks and gluons condense to form a gas of nucleons and light mesons, the latter decay subsequently. At the end of the first second, neutrinos and neutrons decouple from the radiation fluid. The quark--hadron transition and dissipative processes during the first second prepare the initial conditions for the synthesis of the first nuclei. As for the cold dark matter (CDM), WIMPs (weakly interacting massive particles) -- the most popular candidates for the CDM -- decouple from the presently known forms of matter, chemically (freeze-out) at 10 ns and kinetically at 1 ms. The chemical decoupling fixes their present abundances and dissipative processes during and after thermal decoupling set the scale for the very first WIMP clouds.

We consider approaches to cosmological parameter estimation in the inflationary cosmology, focussing on the required accuracy of the initial power spectra. Parametrizing the spectra, for example by power-laws, is well suited to testing the inflationary paradigm but will only correctly estimate cosmological parameters if the parametrization is sufficiently accurate, and we investigate conditions under which this is achieved both for present data and for upcoming satellite data. If inflation is favoured, reliable estimation of its physical parameters requires an alternative approach adopting its detailed predictions. For slow-roll inflation, we investigate the accuracy of the predicted spectra at first and second order in the slow-roll expansion (presenting the complete second-order corrections for the tensors for the first time). We find that within the presently-allowed parameter space, there are regions where it will be necessary to include second-order corrections to reach the accuracy requirements of MAP and Planck satellite data. We end by proposing a data analysis pipeline appropriate for testing inflation and for cosmological parameter estimation from high-precision data.

A new method for predicting inflationary cosmological perturbations, based on the Wentzel-Kramers-Brillouin (WKB) approximation, is presented. A general expression for the WKB scalar and tensor power spectra is derived. The main advantage of the new scheme of approximation is that it is valid even if the slow-roll conditions are violated. The method is applied to power-law inflation, which allows a comparison with an exact result. It is demonstrated that the WKB approximation predicts the spectral indices exactly and the amplitude with an error lower than 10%, even in regimes far from scale-invariance. The new method of approximation is also applied to a situation where the slow-roll conditions hold. It is shown that the result obtained bears close resemblance with the standard slow-roll calculation. Finally, some possible improvements are briefly mentioned.

The lightest supersymmetric particle, most likely the lightest neutralino, is one of the most prominent particle candidates for cold dark matter (CDM). We show that the primordial spectrum of density fluctuations in neutralino CDM has a sharp cut-off, induced by two different damping mechanisms. During the kinetic decoupling of neutralinos, non-equilibrium processes constitute viscosity effects, which damp or even absorb density perturbations in CDM. After the last scattering of neutralinos, free streaming induces neutralino flows from overdense to underdense regions of space. Both damping mechanisms together define a minimal mass scale for perturbations in neutralino CDM, before the inhomogeneities enter the nonlinear epoch of structure formation. We find that the very first gravitationally bound neutralino clouds ought to have masses above 10^{-6} solar masses, which is six orders of magnitude above the mass of possible axion miniclusters.

Specific internal self-interactions in gaseous cosmic fluids are shown to give rise to effective negative pressures which may violate the strong energy condition. On this basis we discuss the transition from an initial de Sitter phase to a subsequent Friedmann-Lema\^{\i}tre-Robertson-Walker period as a non-equilibrium configuration of an ultrarelativistic gas. An accelerated expansion of the present universe is obtained as the consequence of a negative internal friction force which is self-consistently exerted on the non-relativistic cold dark matter (CDM) type microscopic constituents of the cosmic gas.

The fact that the LambdaCDM model fits the observations does not necessarily imply the physical existence of `dark energy'. Dropping the assumption that cold dark matter (CDM) is a perfect fluid opens the possibility to fit the data without dark energy. For imperfect CDM, negative bulk pressure is favoured by thermodynamical arguments and might drive the cosmic acceleration. The coincidence between the onset of accelerated expansion and the epoch of structure formation at large scales might suggest that the two phenomena are linked. A specific example is considered in which effective (anti-frictional) forces, which may be due to dissipative processes during the formation of inhomogeneities, give rise to accelerated expansion of a CDM universe.

In the paper ``Density perturbations in the ekpyrotic scenario'', it is argued that the expected spectrum of primordial perturbations should be scale invariant in this scenario. Here we show that, contrary to what is claimed in that paper, the expected spectrum depends on an arbitrary choice of matching variable. As no underlying (microphysical) principle exists at the present time that could lift the arbitrariness, we conclude that the ekpyrotic scenario is not yet a predictive model.

We investigate the chances that neutron stars act as the lense in a gravitational microlensing event towards the galactic bulge or a spiral arm. The observation of neutron stars by means of gravitational microlensing would allow the estimation of neutron star masses independently of the property of being a pulsar in a binary system. We estimate the contribution of neutron stars to the optical depth and the lensing rate based on two different models of the pulsar distribution in the galaxy. Since only a small fraction of all neutron stars are pulsars, it is unlikely to find a pulsar that acts as a microlense by chance. A position comparison of known radio pulsars with observed microlensing candidates towards the galactic bulge and spiral arms shows no candidate pair, which is consistent with the theoretical expectation. To improve the probability of microlensing a pulsar, we suggest to search for gravitational microlensing events of known nearby high proper motion pulsars. The pulsar PSR J1932+1059 is a good candidate for an astrometric detection of gravitational lensing.

By considering a simplified but exact model for realizing the ekpyrotic scenario, we clarify various assumptions that have been used in the literature. In particular, we discuss the new ekpyrotic prescription for passing the perturbations through the singularity which we show to provide a spectrum depending on a non physical normalization function. We also show that this prescription does not reproduce the exact result for a sharp transition. Then, more generally, we demonstrate that, in the only case where a bounce can be obtained in Einstein General Relativity without facing singularities and/or violation of the standard energy conditions, the bounce cannot be made arbitrarily short. This contrasts with the standard (inflationary) situation where the transition between two eras with different values of the equation of state can be considered as instantaneous. We then argue that the usually conserved quantities are not constant on a typical bounce time scale. Finally, we also examine the case of a test scalar field (or gravitational waves) where similar results are obtained. We conclude that the full dynamical equations of the underlying theory should be solved in a non singular case before any conclusion can be drawn.

The lightest supersymmetric particle, most likely the neutralino, might account for a large fraction of dark matter in the Universe. We show that the primordial spectrum of density fluctuations in neutralino cold dark matter (CDM) has a sharp cut-off due to two damping mechanisms: collisional damping during the kinetic decoupling of the neutralinos at O(10 MeV) and free streaming after last scattering of neutralinos. The cut-off in the primordial spectrum defines a minimal mass for CDM objects in hierarchical structure formation. For typical neutralino and sfermion masses the first gravitationally bound neutralino clouds have masses above 10^(-6) M_\odot.

Accelerated expansion of the Universe may result from an anti-frictional force that is self-consistently exerted on cold dark matter (CDM). Cosmic anti-friction is shown to give rise to an effective negative pressure of the cosmic medium. While other models introduce a component of dark energy besides ``standard'' CDM, we resort to a phenomenological one-component model of CDM with internal self-interactions. We demonstrate how the dynamics of the LambdaCDM model may be recovered as a special case of cosmic anti-friction.

The lightest supersymmetric particle, most likely the neutralino, might account for a large fraction of dark matter in the Universe. We show that the primordial spectrum of density fluctuations in neutralino cold dark matter (CDM) has a sharp cut-off due to two damping mechanisms: collisional damping during the kinetic decoupling of the neutralinos at about 30 MeV (for typical neutralino and sfermion masses) and free streaming after last scattering of neutralinos. The last scattering temperature is lower than the kinetic decoupling temperature by one order of magnitude. The cut-off in the primordial spectrum defines a minimal mass for CDM objects in hierarchical structure formation. For typical neutralino and sfermion masses the first gravitationally bound neutralino clouds have to have masses above 10^(-7) solar masses.

We calculate power spectra of cosmological perturbations at high accuracy for two classes of inflation models. We classify the models according to the behaviour of the Hubble distance during inflation. Our approximation works if the Hubble distance can be approximated either to be a constant or to grow linearly with cosmic time. Many popular inflationary models can be described in this way, e.g., chaotic inflation with a monomial potential, power-law inflation and inflation at a maximum. Our scheme of approximation does not rely on a slow-roll expansion. Thus we can make accurate predictions for some of the models with large slow-roll parameters.

We explain an accelerated expansion of the present universe, suggested from observations of supernovae of type Ia at high redshift, by introducing an anti-frictional force that is self-consistently exerted on the particles of the cosmic substratum. Cosmic anti-friction, which is intimately related to ``particle production'', is shown to give rise to an effective negative pressure of the cosmic medium. While other explanations for an accelerated expansion (cosmological constant, quintessence) introduce a component of dark energy besides ``standard'' cold dark matter (CDM) we resort to a phenomenological one-component model of CDM with internal self-interactions. We demonstrate how the dynamics of the LambdaCDM model may be recovered as a special case of cosmic anti-friction. We discuss the connection with two-component models and obtain an attractor behavior for the ratio of the energy densities of both components which provides a possible phenomenological solution to the coincidence problem.

We present a new approximation scheme that allows us to increase the accuracy of analytical predictions of the power spectra of inflationary perturbations for two specific classes of inflationary models. Among these models are chaotic inflation with a monomial potential, power-law inflation and natural inflation (inflation at a maximum). After reviewing the established first order results we calculate the amplitudes and spectral indices for these classes of models at higher orders in the slow-roll parameters for scalar and tensorial perturbations.

We investigate a new mechanism for the cosmological QCD phase transition: inhomogeneous nucleation. The primordial temperature fluctuations, measured to be $\delta T/T \sim 10^{-5}$, are larger than the tiny temperature interval, in which bubbles would form in the standard picture of homogeneous nucleation. Thus the bubbles nucleate at cold spots. We find the typical distance between bubble centers to be a few meters. This exceeds the estimates from homogeneous nucleation by two orders of magnitude. The resulting baryon inhomogeneities may affect primordial nucleosynthesis.

We emphasize that the estimation of cosmological parameters from cosmic microwave background (CMB) anisotropy data, such as the recent high resolution maps from BOOMERanG and MAXIMA-1, requires assumptions about the primordial spectra. The latter are predicted from inflation. The physically best-motivated scenario is that of slow-roll inflation. However, very often, the unphysical power-law inflation scenario is (implicitly) assumed in the CMB data analysis. We show that the predicted multipole moments differ significantly in both cases. We identify several misconceptions present in the literature (and in the way inflationary relations are often combined in popular numerical codes). For example, we do not believe that, generically, inflation predicts the relation n_S - 1 = n_T for the spectral indices of scalar and tensor perturbations or that gravitational waves are negligible. We calculate the CMB multipole moments for various values of the slow-roll parameters and demonstrate that an important part of the space of parameters (n_S, n_T) has been overlooked in the CMB data analysis so far.

From a general ansatz for the effective potential of cosmological perturbations we find new, exact solutions in single-scalar-field inflation: a three parameter family of exact inflationary solutions that encompasses all exact solutions that have been known previously (power-law inflation, Easther's model, and a generalised version of Starobinsky's solution). The main feature of this new family is that the spectral indices are scale dependent.

We analyze the effect of primordial density perturbations on the cosmic QCD phase transition. According to our results hadron bubbles nucleate at the cold perturbations. We call this mechanism inhomogeneous nucleation. We find the typical distance between bubble centers to be a few meters. This exceeds the estimates from homogeneous nucleation by two orders of magnitude. The resulting baryon inhomogeneities may affect primordial nucleosynthesis.

We analyze a new mechanism for the cosmological QCD first-order phase transition: inhomogeneous nucleation. The primordial temperature fluctuations are larger than the tiny temperature interval, in which bubbles would form in the standard picture of homogeneous nucleation. Thus the bubbles nucleate at cold spots. We find the typical distance between bubble centers to be a few meters. This exceeds the estimates from homogeneous nucleation by two orders of magnitude. The resulting baryon inhomogeneities may affect primordial nucleosynthesis.

Inflationary predictions of the cosmic microwave background anisotropy are often based on the slow-roll approximation. We study the precision of these predictions and compare them with the recent data from BOOMERanG and MAXIMA-1.

Inflationary predictions for the anisotropy of the cosmic microwave background radiation (CMBR) are often based on the slow-roll approximation. We study the precision with which the multipole moments of the temperature two-point correlation function can be predicted by means of the slow-roll approximation. We ask whether this precision is good enough for the forthcoming high precision observations by means of the MAP and Planck satellites. The error in the multipole moments due to the slow-roll approximation is demonstrated to be bigger than the error in the power spectrum. For power-law inflation with $n_S=0.9$ the error from the leading order slow-roll approximation is $\approx 5%$ for the amplitudes and $\approx 20%$ for the quadrupoles. For the next-to-leading order the errors are within a few percent. The errors increase with $|n_S - 1|$. To obtain a precision of 1% it is necessary, but in general not sufficient, to use the next-to-leading order. In the case of power-law inflation this precision is obtained for the spectral indices if $|n_S - 1| < 0.02$ and for the quadrupoles if $|n_S - 1| < 0.15$ only. The errors in the higher multipoles are even larger than those for the quadrupole, e.g. $\approx 15%$ for l=100, with $n_S = 0.9$ at the next-to-leading order. We find that the accuracy of the slow-roll approximation may be improved by shifting the pivot scale of the primordial spectrum (the scale at which the slow-roll parameters are fixed) into the regime of acoustic oscillations. Nevertheless, the slow-roll approximation cannot be improved beyond the next-to-leading order in the slow-roll parameters.

We investigate the clumping of cold dark matter (CDM) at small scales. If the CDM particle is the neutralino, we find that collisional damping during its kinetic decoupling from the radiation fluid and free streaming introduce a small-scale cut-off in the primordial power spectrum of CDM. This cut-off sets the scale for the very first CDM objects in the Universe, which we expect to have a mass of $\sim 10^{-12} M_\odot$. For non-thermal CDM candidates, such as axions, wimpzillas, or primordial black holes, the cosmological QCD transition might induce features in the primordial spectrum at similar mass scales.

The cosmological QCD transition affects primordial density perturbations. If the QCD transition is first order, the sound speed vanishes during the transition and density perturbations fall freely. For scales below the Hubble radius at the transition the primordial Harrison-Zel'dovich spectrum of density fluctuations develops large peaks and dips. These peaks grow with wave number for both the radiation fluid and for cold dark matter (CDM). The peaks in the radiation fluid are wiped out during neutrino decoupling. For cold dark matter that is kinetically decoupled at the QCD transition (e.g. axions) these peaks lead to the formation of CDM clumps of masses $10^{-20} M_\odot< M < 10^{-10} M_\odot$.

The spectrum of gravitational waves that have been produced in inflation is modified during cosmological transitions. Large drops in the number of relativistic particles, like during the QCD transition or at $e^+e^-$ annihilation, lead to steps in the spectrum of gravitational waves. We calculate the transfer function for the differential energy density of gravitational waves for a first-order and for a crossover QCD transition.

The cosmological QCD transition affects primordial density perturbations. If the QCD transition is first order, the sound speed vanishes during the transition and density perturbations fall freely. For scales below the Hubble radius at the transition the primordial Harrison-Zel'dovich spectrum of density fluctuations develops large peaks and dips. These peaks grow with wave number for both the hadron-photon-lepton fluid and for cold dark matter. At the horizon scale the enhancement is small. This by itself does not lead to the formation of black holes at the QCD transition. The peaks in the hadron-photon-lepton fluid are wiped out during neutrino decoupling. For cold dark matter that is kinetically decoupled at the QCD transition (e.g., axions or primordial black holes) these peaks lead to the formation of CDM clumps of masses $10^{-20} M_\odot< M_{\rm clump} < 10^{-10} M_\odot$.

The different phases of QCD at finite temperature and density lead to interesting effects in cosmology and astrophysics. In this work I review some aspects of the cosmological QCD transition and of astrophysics at high baryon density.

Inflation predicts the generation of cosmological perturbations. Usually, the power spectra for the scalar and tensor modes are calculated with help of the slow roll approximation. In the case of power law inflation an exact result is available. We compare the predictions for the cosmic microwave background anisotropies from the slow roll approximation with the exact results from power law inflation. We find that the so-called consistency check from the slow roll approximation, $C_2^T/C_2^S \approx -6.93 n_T$, may differ considerably from the exact result.

We study the influence of reheating on super-horizon density perturbations and gravitational waves. We correct wrong claims [L. P. Grishchuk Phys. Rev. D 50, 7154 (1994) and gr-qc/9801011] about the joining of perturbations at cosmological transitions and about the quantization of cosmological perturbations.

We reply to L.P. Grishchuk's comment (gr-qc/9801011) on our paper ``The Influence of Cosmological Transitions on the Evolution of Density Perturbations'' (Phys. Rev. D 57, 3302, 1998). We show that his points of criticism are not correct.

We study the influence of the reheating and equality transitions on superhorizon density perturbations and gravitational waves. Recent criticisms of the `standard result' for large-scale perturbations in inflationary cosmology are rectified. The claim that the `conservation law' for the amplitude of superhorizon modes was empty is shown to be wrong. For sharp transitions, i.e. the pressure jumps, we rederive the Deruelle-Mukhanov junction conditions. For a smooth transition we correct a result obtained by Grishchuk recently. We show that the junction conditions are not crucial, because the pressure is continuous during the reheating transition. The problem occurred, because the perturbed metric was not evolved correctly through the smooth reheating transition. Finally, we derive the `standard result' within Grishchuks's smooth (reheating) transition.

We investigate the effect of the cosmological QCD transition on the evolution of primordial density perturbations. If the phase transition is first order, the sound speed vanishes during the transition and density perturbations fall freely. The Harrison-Zel'dovich spectrum of density fluctuations develops large peaks on scales below the Hubble radius at the transition. These peaks above the primordial spectrum grow with wavenumber and produce cold dark matter clumps of masses less than $10^{-10} M_\odot$. At the horizon scale the amplification of overdensities is of order unity, and thus no $1 M_\odot$ black hole production is possible for a COBE normalized scale-invariant spectrum.

The quark-gluon to hadron transition affects the evolution of cosmological perturbations. If the phase transition is first order, the sound speed vanishes during the transition, and density perturbations fall freely. This distorts the primordial Harrison-Zel'dovich spectrum of density fluctuations below the Hubble scale at the transition. Peaks are produced, which grow at most linearly in wavenumber, both for the hadron-photon-lepton fluid and for cold dark matter. For cold dark matter which is kinetically decoupled well before the QCD transition clumps of masses below $10^{-10} M_\odot$ are produced.

We investigate the effect of the quark-gluon to hadron transition on the evolution of cosmological perturbations. If the phase transition is first order, the sound speed vanishes during the transition, and density perturbations fall freely. The primordial Harrison-Zel'dovich density fluctuations for scales below the Hubble radius at the transition develop peaks, which grow linearly with the wavenumber, both for the hadron-photon-lepton fluid and for cold dark matter. The large peaks in the spectrum produce cold dark matter clumps of $10^{-8}$ to $10^{-11} M_\odot$.

Scalar cosmological perturbations of a weakly self-interacting plasma mixed with a perfect radiation fluid are investigated. Effects of this plasma are considered through order $\lambda^{3/2}$ of perturbative thermal-field-theory in the radiation dominated universe. The breakdown of thermal perturbation theory at vastly subhorizon scales is circumvented by a Pad\'e approximant solution. Compared to collisionless plasmas the phase speed and subhorizon damping of the plasma density perturbations are changed. An example for a self-interacting thermal field is provided by the neutrinos with effective 4-fermion interactions.

Using a recent thermal-field-theory approach to cosmological perturbations, the exact solutions that were found for collisionless ultrarelativistic matter are generalized to include the effects from weak self-interactions in a $\lambda\phi^4$ model through order $\lambda^{3/2}$. This includes the effects of a resummation of thermal masses and associated nonlocal gravitational vertices, thus going far beyond classical kinetic theory. Explicit solutions for all the scalar, vector, and tensor modes are obtained for a radiation-dominated Einstein-de Sitter model containing a weakly interacting scalar plasma with or without the admixture of an independent component of perfect radiation fluid.

The exact solutions for linear cosmological perturbations which have been obtained for collisionless relativistic matter within thermal field theory are extended to a self-interacting case. The two-loop contributions of scalar $\lambda\phi^4$ theory to the thermal graviton self-energy are evaluated, which give the $O(\lambda)$ corrections in the perturbation equations. The changes are found to be perturbative on scales comparable to or larger than the Hubble horizon, but the determination of the large-time damping behavior of subhorizon perturbations requires a resummation of thermally induced masses.

All regular and singular cosmological perturbations in a radiation dominated Einstein-de Sitter Universe with collisionless particles can be found by a generalized power series ansatz. Talk given at "Birth of the Universe and Fundamental Physics", May 1994 (Rome).

A closed set of equations for the evolution of linear perturbations of homogeneous, isotropic cosmological models can be obtained in various ways. The simplest approach is to assume a macroscopic equation of state, e.g.\ that of a perfect fluid. For a more refined description of the early universe, a microscopic treatment is required. The purpose of this paper is to compare the approach based on classical kinetic theory to the more recent thermal-field-theory approach. It is shown that in the high-temperature limit the latter describes cosmological perturbations supported by collisionless, massless matter, wherein it is equivalent to the kinetic theory approach. The dependence of the perturbations in a system of a collisionless gas and a perfect fluid on the initial data is discussed in some detail. All singular and regular solutions are found analytically.

We obtain analytic solutions for the density contrast and the anisotropic pressure in a multi-dimensional FRW cosmology with collisionless, massless matter. These are compared with perturbations of a perfect fluid universe. To describe the metric perturbations we use manifest gauge invariant metric potentials. The matter perturbations are calculated by means of (automatically gauge invariant) finite temperature field theory, instead of kinetic theory. (Talk given at the Journ\'ees Relativistes '93, 5 -- 7 April, Brussels, Belgium)