From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, microphysics, and chemistry. In order to understand these key microphysical parameters, we present a one-dimensional non-ideal magnetohydrodynamics code and apply it to a model of a time-dependent, oblique, magnetic shock wave. By varying the microphysical ingredients of the model, we find that cosmic rays and dust play a major role, and that, despite the uncertainties, the inclusion of microphysics is essential to obtain a realistic outcome in magnetic astrophysical simulations.

We investigate the formation of protoplanetary disks around nine solar mass stars formed in the context of a (40 pc)$^3$ Giant Molecular Cloud model, using RAMSES adaptive-mesh refinement simulations extending over a scale range of about 4 million, from an outer scale of 40 pc down to cell sizes of 2 AU. Our most important result is that the accretion process is heterogeneous in multiple ways; in time, in space, and among protostars of otherwise similar mass. Accretion is heterogeneous in time, in the sense that accretion rates vary during the evolution, with generally decreasing profiles, whose slopes vary over a wide range, and where accretion can increase again if a protostar enters a region with increased density and low speed. Accretion is heterogeneous in space, because of the mass distribution, with mass approaching the accreting star-disk system in filaments and sheets. Finally, accretion is heterogeneous among stars, since the detailed conditions and dynamics in the neighborhood of each star can vary widely. We also investigate the sensitivity of disk formation to physical conditions, and test their robustness by varying numerical parameters. We find that disk formation is robust even when choosing the least favorable sink particle parameters, and that turbulence cascading from larger scales is a decisive factor in disk formation. We also investigate the transport of angular momentum, finding that the net inward mechanical transport is compensated for mainly by an outward directed magnetic transport, with a contribution from gravitational torques usually subordinate to the magnetic transport.

Dust plays a key role in the evolution of the ISM and its correct modelling in numerical simulations is therefore fundamental. We present a new and self-consistent model that treats grain thermal coupling with the gas, radiation balance, and surface chemistry for molecular hydrogen. This method can be applied to any dust distribution with an arbitrary number of grain types without affecting the overall computational cost. In this paper we describe in detail the physics and the algorithm behind our approach, and in order to test the methodology, we present some examples of astrophysical interest, namely (i) a one-zone collapse with complete gas chemistry and thermochemical processes, (ii) a 3D model of a low-metallicity collapse of a minihalo starting from cosmological initial conditions, and (iii) a turbulent molecular cloud with H-C-O chemistry (277 reactions), together with self-consistent cooling and heating solved on the fly. Although these examples employ the publicly available code KROME, our approach can be easily integrated into any computational framework.

The short-lived $^{26}$Al and $^{60}$Fe radionuclides are synthesized and expelled in the interstellar medium by core-collapse supernova events. The solar system's first solids, calcium-aluminium refractory inclusions (CAIs), contain evidence for the former presence of the $^{26}$ Al nuclide defining the canonical $^{26}$Al/$^{27}$ Al ratio of $\sim5 \times10^{-5}$. A different class of objects temporally related to canonical CAIs are CAIs with fractionation and unidentified nuclear effects (FUN CAIs), which record a low initial $^{26}$Al/$^{27}$Al of $10^{-6}$. The contrasting level of $^{26}$Al between these objects is often interpreted as reflecting the admixing of the $^{26}$Al nuclide during the early formative phase of the Sun. We use giant molecular cloud (GMC) scale adaptive mesh-refinement numerical simulations to trace the abundance of $^{26}$Al and $^{60}$Fe in star-forming gas during the early stages of accretion of individual low mass protostars. We find that the $^{26}$Al/$^{27}$Al and $^{60}$Fe/$^{56}$Fe ratios of accreting gas within a vicinity of 1000 AU of the stars follow the predicted decay curves of the initial abundances at time of star formation without evidence of spatial or temporal heterogeneities for the first 100 kyr of star formation. Therefore, the observed differences in $^{26}$Al/$^{27}$Al ratios between FUN and canonical CAIs are likely not caused by admixing of supernova material during the early evolution of the proto-Sun. Selective thermal processing of dust grains is a more viable scenario to account for the heterogeneity in $^{26}$Al/$^{27}$Al ratios at the time of solar system formation.

Turbulence is ubiquitous in molecular clouds (MCs), but its origin is still unclear because MCs are usually assumed to live longer than the turbulence dissipation time. Interstellar medium (ISM) turbulence is likely driven by SN explosions, but it has never been demonstrated that SN explosions can establish and maintain a turbulent cascade inside MCs consistent with the observations. In this work, we carry out a simulation of SN-driven turbulence in a volume of (250 pc)$^3$, specifically designed to test if SN driving alone can be responsible for the observed turbulence inside MCs. We find that SN driving establishes a velocity scaling consistent with the usual scaling laws of supersonic turbulence, suggesting that previous idealized simulations of MC turbulence, driven with a random, large-scale volume force, were correctly adopted as appropriate models for MC turbulence, despite the artificial driving. We also find that the same scaling laws extend to the interior of MCs, and that the velocity-size relation of the MCs selected from our simulation is consistent with that of MCs from the Outer-Galaxy Survey, the largest MC sample available. The mass-size relation and the mass and size probability distributions also compare successfully with those of the Outer Galaxy Survey. Finally, we show that MC turbulence is super-Alfv\'{e}nic with respect to both the mean and rms magnetic-field strength. We conclude that MC structure and dynamics are the natural result of SN-driven turbulence.

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

We use the adaptive mesh refinement code RAMSES to model the formation of protoplanetary disks in realistic star formation environments. The resolution scales over up to 29 powers of two ($\sim$ 9 orders of magnitude) covering a range from outer scales of 40 pc to inner scales of 0.015 AU. The accretion rate from a 1.5 solar mass envelope peaks near $10^{-4}$ $\mspy$ about 6 kyr after sink particle formation and then decays approximately exponentially, reaching $10^{-6}$ $\mspy$ in 100 kyr. The models suggest universal scalings of physical properties with radius during the main accretion phase, with kinetic and / or magnetic energy in approximate balance with gravitational energy. Efficient accretion is made possible by the braking action of the magnetic field, which nevertheless allows a near-Keplerian disk to grow to a 100 AU size. The magnetic field strength ranges from more than 10 G at 0.1 AU to less than 1 mG at 100 AU, and drives a time dependent bipolar outflow, with a collimated jet and a broader disk wind.

We present the Photon-Plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.

We show that supersonic MHD turbulence yields a star formation rate (SFR) as low as observed in molecular clouds (MCs), for characteristic values of the free-fall time divided by the dynamical time, $t_{\rm ff}/t_{\rm dyn}$, the alfv\'{e}nic Mach number, ${\cal M}_{\rm a}$, and the sonic Mach number, ${\cal M}_{\rm s}$. Using a very large set of deep adaptive-mesh-refinement simulations, we quantify the dependence of the SFR per free-fall time, $\epsilon_{\rm ff}$, on the above parameters. Our main results are: i) $\epsilon_{\rm ff}$ decreases exponentially with increasing $t_{\rm ff}/t_{\rm dyn}$, but is insensitive to changes in ${\cal M}_{\rm s}$, for constant values of $t_{\rm ff}/t_{\rm dyn}$ and ${\cal M}_{\rm a}$. ii) Decreasing values of ${\cal M}_{\rm a}$ (stronger magnetic fields) reduce $\epsilon_{\rm ff}$, but only to a point, beyond which $\epsilon_{\rm ff}$ increases with a further decrease of ${\cal M}_{\rm a}$. iii) For values of ${\cal M}_{\rm a}$ characteristic of star-forming regions, $\epsilon_{\rm ff}$ varies with ${\cal M}_{\rm a}$ by less than a factor of two. We propose a simple star-formation law, based on the empirical fit to the minimum $\epsilon_{\rm ff}$, and depending only on $t_{\rm ff}/t_{\rm dyn}$: $\epsilon_{\rm ff} \approx \epsilon_{\rm wind} \exp(-1.6 \,t_{\rm ff}/t_{\rm dyn})$. Because it only depends on the mean gas density and rms velocity, this law is straightforward to implement in simulations and analytical models of galaxy formation and evolution.

We study the non-linear gravitational collapse of dark matter into halos through numerical N-body simulations of Lemaitre-Tolman-Bondi void models. We extend the halo mass function formalism to these models in a consistent way. This extension not only compares well with the simulated data at all times and radii, but it also gives interesting clues about the impact of the background shear on the growth of perturbations. Our results give hints about the possibility of constraining the background shear via cluster number counts, which could then give rise to strong constraints on general inhomogeneous models, of any scale.

Spontaneous rapid growth of strong magnetic fields is rather ubiquitous in high-energy density environments ranging from astrophysical sources (e.g., gamma-ray bursts and relativistic shocks), to reconnection, to laser-plasma interaction laboratory experiments, where they are produced by kinetic streaming instabilities of the Weibel type. Relativistic electrons propagating through these sub-Larmor-scale magnetic fields radiate in the jitter regime, in which the anisotropy of the magnetic fields and the particle distribution have a strong effect on the produced radiation. Here we develop the general theory of jitter radiation, which includes (i) anisotropic magnetic fields and electron velocity distributions, (ii) the effects of trapped electrons and (iii) extends the description to large deflection angles of radiating particles thus establishing a cross-over between the classical jitter and synchrotron regimes. Our results are in remarkable agreement with the radiation spectra obtained from particle-in-cell simulations of the classical Weibel instability. Particularly interesting is the onset of the field growth, when the transient hard synchrotron-violating spectra are common as a result of the dominant role of the trapped population. This effect can serve as a distinct observational signature of the violent field growth in astrophysical sources and lab experiments. It is also interesting that a system with small-scale fields tends to evolve toward the small-angle jitter regime, which can, under certain conditions, dominate the overall emission of a source.

We perform numerical simulations of large scale structure evolution in an inhomogeneous Lemaitre-Tolman-Bondi (LTB) model of the Universe. We follow the gravitational collapse of a large underdense region (a void) in an otherwise flat matter-dominated Einstein-deSitter model. We observe how the (background) density contrast at the centre of the void grows to be of order one, and show that the density and velocity profiles follow the exact non-linear LTB solution to the full Einstein equations for all but the most extreme voids. This result seems to contradict previous claims that fully relativistic codes are needed to properly handle the non-linear evolution of large scale structures, and that local Newtonian dynamics with an explicit expansion term is not adequate. We also find that the (local) matter density contrast grows with the scale factor in a way analogous to that of an open universe with a value of the matter density OmegaM(r) corresponding to the appropriate location within the void.

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to arise from relativistically shocked collisionless plasmas. Such sources often exhibit highly transient spectra evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such circumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we present synthetic radiation spectra representing the early stage development of the filamentation (streaming) instability of an initially unmagnetized plasma, which is relevant for both collisionless shock formation and reconnection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. Results were obtained using a highly efficient "in situ" diagnostics method, based on detailed particle-in-cell modeling of collisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra exhibit dependence on the presence - or absence - of an inert plasma constituent, when comparing baryonic plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be taken when using lower-dimensional models to obtain information about the astrophysical phenomena generating observed spectra.

Two dimensional modeling of collisionless shocks has been of tremendous importance in understanding the physics of the non-linear evolution, momentum transfer and particle acceleration, but current computer capacities have now reached a point where three dimensional modeling is becoming feasible. We present the first three dimensional model of a fully developed and relaxed relativistic ion-electron shock, and analyze and compare it to similar 2D models. Quantitative and qualitative differences are found with respect to the two-dimensional models. The shock jump conditions are of course different, because of the extra degree of freedom, but in addition it is found that strong parallel electric fields develop at the shock interface, the level of magnetic field energy is lower, and the non-thermal particle distribution is shallower with a powerlaw index of ~2.2.

We use N-body simulations to find the effect of neutrino masses on halo properties, and investigate how the density profiles of both the neutrino and the dark matter components change as a function of the neutrino mass. We compare our neutrino density profiles with results from the N-one-body method and find good agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos is not saturated. Finally we study how the halo mass function changes as a function of the neutrino mass and compare our results with the Sheth-Tormen semi-analytic formulae. Our results are important for surveys which aim at probing cosmological parameters using clusters, as well as future experiments aiming at measuring the cosmic neutrino background directly.

We study the primordial non-Gaussinity predicted from simple models of inflation with a linear potential and superimposed oscillations. This generic form of the potential is predicted by the axion monodromy inflation model, that has recently been proposed as a possible realization of chaotic inflation in string theory, where the monodromy from wrapped branes extends the range of the closed string axions to beyond the Planck scale. The superimposed oscillations in the potential can lead to new signatures in the CMB spectrum and bispectrum. In particular the bispectrum will have a new distinct shape. We calculate the power spectrum and bispectrum of curvature perturbations in the model, as well as make analytic estimates in various limiting cases. From the numerical analysis we find that for a wide range of allowed parameters the model produces a feature in the bispectrum with fnl ~ 50 or larger while the power spectrum is almost featureless. This model is therefore an example of a string-inspired inflationary model which is testable mainly through its non-Gaussian features. Finally we provide a simple analytic fitting formula for the bispectrum which is accurate to approximately 5% in all cases, and easily implementable in codes designed to provide non-Gaussian templates for CMB analyses.

We perform an autocorrelation study of the Auger data with the aim to constrain the number density n_s of ultrahigh energy cosmic ray (UHECR) sources, estimating at the same time the effect on $n_s$ of the systematic energy scale uncertainty and of the distribution of UHECR. The use of global analysis has the advantage that no biases are introduced, either in n_s or in the related error bar, by the a priori choice of a single angular scale. The case of continuous, uniformly distributed sources is nominally disfavored at 99% C.L. and the fit improves if the sources follow the large-scale structure of matter in the universe. The best fit values obtained for the number density of proton sources are within a factor ~2 around n_s ~ 10^{-4}/Mpc^3 and depend mainly on the Auger energy calibration scale, with lower densities being preferred if the current scale is correct. The data show no significant small-scale clustering on scales smaller than a few degrees. This might be interpreted as a signature of magnetic smearing of comparable size, comparable with the indication of a ~3 degrees magnetic deflection coming from cross-correlation results. The effects on the above results of some approximations done is also discussed.

As an alternative explanation of the dimming of distant supernovae it has recently been advocated that we live in a special place in the Universe near the centre of a large spherical void described by a Lemaitre-Tolman-Bondi (LTB) metric. In this scenario, the Universe is no longer homogeneous and isotropic, and the apparent late time acceleration is actually a consequence of spatial gradients. We propose in this paper a new observable, the normalized cosmic shear, written in terms of directly observable quantities, and calculable in arbitrary inhomogeneous cosmologies. This will allow future surveys to determine whether we live in a homogeneous universe or not. In this paper we also update our previous observational constraints from geometrical measures of the background cosmology. We include the Union Supernovae data set of 307 Type Ia supernovae, the CMB acoustic scale and the first measurement of the radial baryon acoustic oscillation scale. Even though the new data sets are significantly more constraining, LTB models -- albeit with slightly larger voids -- are still in excellent agreement with observations, at chi^2/d.o.f. = 307.7/(310-4)=1.005. Together with the paper we also publish the updated easyLTB code used for calculating the models and for comparing them to the observations.

The possibility that we live in a special place in the universe, close to the centre of a large void, seems an appealing alternative to the prevailing interpretation of the acceleration of the universe in terms of a LCDM model with a dominant dark energy component. In this paper we confront the asymptotically flat Lemaitre-Tolman-Bondi (LTB) models with a series of observations, from Type Ia Supernovae to Cosmic Microwave Background and Baryon Acoustic Oscillations data. We propose two concrete LTB models describing a local void in which the only arbitrary functions are the radial dependence of the matter density Omega_M and the Hubble expansion rate H. We find that all observations can be accommodated within 1 sigma, for our models with 4 or 5 independent parameters. The best fit models have a chi^2 very close to that of the LCDM model. We perform a simple Bayesian analysis and show that one cannot exclude the hypothesis that we live within a large local void of an otherwise Einstein-de Sitter model.

Modeling plasma physical processes in astrophysical context demands for both detailed kinetics and large scale development of the electromagnetic field densities. We present a new framework for modeling plasma physics of hot tenuous plasmas by a two-split scheme, in which the large scale fields are modeled by means of a particle-in-cell (PIC) code, and in which binary collision processes and single-particle processes are modeled through a Monte-Carlo approach. Our novel simulation tool -- the PhotonPlasma code -- is a unique hybrid model; it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting of particles and a sub-Debye Monte-Carlo binary particle interaction framework. As an illustration of the capabilities we present results from a numerical study of Gamma-Ray Burst - Circumburst Medium interaction and plasma preconditioning via Compton scattering. We argue that important microphysical processes can only viably be investigated by means of hybrid codes such as the PhotonPlasma code. Our first results from 3D simulations with this new simulation tool suggest that magnetic fields and plasma filaments are created in the wakefield of prompt gamma-ray bursts. Furthermore, the photon flux density gradient impacts on particle acceleration in the burst head and wakefield. We discuss some possible implications of the circumburst medium being preconditioned for a trailing afterglow shock front. We also discuss important improvements for future studies of GRB wakefields processes, using the PhotonPlasma code.

We have performed detailed studies of non-linear structure formation in cosmological models with light neutrinos. For the first time the effect of neutrino thermal velocities has been included in a consistent way, and the effect on the matter power spectrum is found to be significant. The effect is large enough to be measured in future, high precision surveys. Additionally, we provide a simple but accurate analytic expression for the suppression of fluctuation power due to massive neutrinos. Finally, we describe a simple and fast method for including the effect of massive neutrinos in large-scale N-body simulations which is accurate at the 1% level for \sum m_nu < 0.15 eV.

As an alternative explanation of the dimming of distant supernovae it has recently been advocated that we live in a special place in the Universe near the centre of a large void described by a Lemaitre-Tolman-Bondi (LTB) metric. The Universe is no longer homogeneous and isotropic and the apparent late time acceleration is actually a consequence of spatial gradients in the metric. If we did not live close to the centre of the void, we would have observed a Cosmic Microwave Background (CMB) dipole much larger than that allowed by observations. Hence, until now it has been argued, for the model to be consistent with observations, that by coincidence we happen to live very close to the centre of the void or we are moving towards it. However, even if we are at the centre of the void, we can observe distant galaxy clusters, which are off-centre. In their frame of reference there should be a large CMB dipole, which manifests itself observationally for us as a kinematic Sunyaev-Zeldovich (kSZ) effect. kSZ observations give far stronger constraints on the LTB model compared to other observational probes such as Type Ia Supernovae, the CMB, and baryon acoustic oscillations. We show that current observations of only 9 clusters with large error bars already rule out LTB models with void sizes greater than approximately 1.5 Gpc and a significant underdensity, and that near future kSZ surveys like the Atacama Cosmology Telescope, South Pole Telescope, APEX telescope, or the Planck satellite will be able to strongly rule out or confirm LTB models with giga parsec sized voids. On the other hand, if the LTB model is confirmed by observations, a kSZ survey gives a unique possibility of directly reconstructing the expansion rate and underdensity profile of the void.

The extragalactic cosmic gamma-ray background (CGB) is an interesting channel to look for signatures of dark matter annihilation. In particular, besides the imprint in the energy spectrum, peculiar anisotropy patterns are expected compared to the case of a pure astrophysical origin of the CGB. We take into account the uncertainties in the dark matter clustering properties on sub-galactic scales, deriving two possible anisotropy scenarios. A clear dark matter angular signature is achieved when the annihilation signal receives only a moderate contribution from sub-galactic clumps and/or cuspy haloes. Experimentally, if galactic foregrounds systematics are efficiently kept under control, the angular differences are detectable with the forthcoming GLAST observatory, provided that the annihilation signal contributes to the CGB for a fraction >10-20%. If, instead, sub-galactic structures have a more prominent role, the astrophysical and dark matter anisotropies become degenerate, correspondingly diluting the DM signature. As complementary observables we also introduce the cross-correlation between surveys of galaxies and the CGB and the cross-correlation between different energy bands of the CGB and we find that they provide a further sensitive tool to detect the dark matter angular signatures.

The arrival directions of ultrahigh energy cosmic rays (UHECRs) may show anisotropies on all scales, from just above the experimental angular resolution up to medium scales and dipole anisotropies. We find that a global comparison of the two-point auto-correlation function of the data with the one of catalogues of potential sources is a powerful diagnostic tool. In particular, this method is far less sensitive to unknown deflections in magnetic fields than cross-correlation studies while keeping a strong discrimination power among source candidates. We illustrate these advantages by considering ordinary galaxies, gamma ray bursts and active galactic nuclei as possible sources. Already the sparse publicly available data suggest that the sources of UHECRs may be a strongly clustered sub-sample of galaxies or of active galactic nuclei. We present forecasts for various cases of source distributions which can be checked soon by the Pierre Auger Observatory.

The Sunyaev-Zeldovich (SZ) effect gives a measure of the thermal energy and electron pressure in groups and clusters of galaxies. In the near future SZ surveys will map hundreds of systems, shedding light on the pressure distribution in the systems. The thermal energy is related to the total mass of a system of galaxies, but it is only a projection that is observed through the SZ effect. A model for the 3D distribution of pressure is needed to link the SZ signal to the total mass of the system. In this work we construct an empirical model for the 2D and 3D SZ profile, and compare it to a set of realistic high resolution SPH simulations of galaxy clusters and groups, and to a stacked SZ profile for massive clusters derived from WMAP data. Furthermore, we combine observed temperature profiles with dark matter potentials to yield an additional constraint, under the assumption of hydrostatic equilibrium. We find a very tight correlation between the characteristic scale in the model, the integrated SZ signal, and the total mass in the systems with a scatter of only 4%. The model only contains two free parameters, making it readily applicable even to low resolution SZ observations of galaxy clusters. A fitting routine for the model that can be applied to observed or simulated data can be found at http://www.phys.au.dk/~haugboel/software.shtml

Type Ia supernovae are currently the best known standard candles at cosmological distances. In addition to providing a powerful probe of dark energy they are an ideal source of information about the peculiar velocity field of the local universe. Even with the very small number of supernovae presently available it has been possible to measure the dipole and quadrupole of the local velocity field out to z~0.025. With future continuous all-sky surveys like the LSST project the luminosity distances of tens of thousands of nearby supernovae will be measured accurately. This will allow for a determination of the local velocity structure of the universe as a function of redshift with unprecedented accuracy, provided the redshifts of the host galaxies are known. Using catalogues of mock surveys we estimate that future low redshift supernova surveys will be able to probe sigma-8 to a precision of roughly 5% at 95% C.L. This is comparable to the precision in future galaxy and weak lensing surveys and with a relatively modest observational effort it will provide a crucial cross-check on future measurements of the matter power spectrum.

We present a measurement of the velocity flow of the local universe relative to the CMB rest frame, based on the Jha, Riess & Kirshner (2007) sample of 133 low redshift type Ia supernovae. At a depth of 4500 km/s we find a dipole amplitude of 279+-68 km/s in the direction (l,b) = (285+-18,-10+-15), consistent with earlier measurements and with the assumption that the local velocity field is dominated by the Great Attractor region. At a larger depth of 5900 km/s we find a shift in the dipole direction towards the Shapley concentration. We also present the first measurement of the quadrupole term in the local velocity flow at these depths. Finally, we have performed detailed studies based on N-body simulations of the expected precision with which the lowest multipoles in the velocity field can be measured out to redshifts of order 0.1. Our mock catalogues are in good agreement with current observations, and demonstrate that our results are robust with respect to assumptions about the influence of local environment on the type Ia supernova rate.

In this thesis different numerical methods, as well as applications of the methods to a number of current problems in relativistic astrophysics, are presented. In the first part the theoretical foundation and numerical implementation of a new general relativistic magnetohydrodynamics code is discussed. A new form of the equations of motion using global coordinates, but evolving the dynamical variables from the point of view of a local observer is presented. No assumptions are made about the background metric and the design is ready to be coupled with methods solving the full Einstein equations. In the second part of the thesis important results concerning the understanding of collisionless shocks, obtained from experiments with a relativistic charged particle code, are presented. Relativistic collisionless shocks are important in a range of astrophysical objects; in particular in gamma ray burst afterglows and other relativistic jets. It is shown that a strong small scale, fluctuating, and predominantly transversal magnetic field is unavoidably generated by a two-stream instability. The magnetic energy density reaches a few percent of equipartition. A new acceleration mechanism for electrons in ion-electron collisionless shocks is proposed. The mechanism is capable of creating a powerlaw electron distribution in a collisionless shocked region. The non-thermal acceleration of the electrons is directly related to the ion current channels generated by the two-stream instability and is local in nature. Thus the observed radiation field may be tied directly to the local conditions of the plasma and could be a strong handle on the physical processes. (abridged)

Gamma ray bursts are among the most energetic events in the known universe. A highly relativistic fireball is ejected. In most cases the burst itself is followed by an afterglow, emitted under deceleration as the fireball plunges through the circum-stellar media. To interpret the observations of the afterglow emission, two physical aspects need to be understood: 1) The origin and nature of the magnetic field in the fireball and 2) the particle velocity distribution function behind the shock. Both are necessary in existing afterglow models to account for what is believed to be synchrotron radiation. To answer these questions, we need to understand the microphysics at play in collisionless shocks. Using 3D particle-in-cell simulations we can gain insight in the microphysical processes that take place in such shocks. We discuss the results of such computer experiments. It is shown how a Weibel-like two-stream plasma instability is able to create a strong transverse intermittent magnetic field and points to a connected mechanism for in situ particle acceleration in the shock region.

The outflows from gamma ray bursts, active galactic nuclei and relativistic jets in general interact with the surrounding media through collisionless shocks. With three dimensional relativistic particle-in-cell simulations we investigate such shocks. The results from these experiments show that small--scale magnetic filaments with strengths of up to percents of equipartition are generated and that electrons are accelerated to power law distributions N(E)~E^{-p} in the vicinity of the filaments through a new acceleration mechanism. The acceleration is locally confined, instantaneous and differs from recursive acceleration processes such as Fermi acceleration. We find that the proposed acceleration mechanism competes with thermalization and becomes important at high Lorentz factors.

Collisionless plasma shock theory, which applies for example to the afterglow of gamma ray bursts, still contains key issues that are poorly understood. In this paper we study charged particle dynamics in a highly relativistic collisionless shock numerically using ~10^9 particles. We find a power law distribution of accelerated electrons, which upon detailed investigation turns out to originate from an acceleration mechanism that is decidedly different from Fermi acceleration. Electrons are accelerated by strong filamentation instabilities in the shocked interpenetrating plasmas and coincide spatially with the power law distributed current filamentary structures. These structures are an inevitable consequence of the now well established Weibel-like two-stream instability that operates in relativistic collisionless shocks. The electrons are accelerated and decelerated instantaneously and locally; a scenery that differs qualitatively from recursive acceleration mechanisms such as Fermi acceleration. The slopes of the electron distribution power laws are in concordance with the particle power law spectra inferred from observed afterglow synchrotron radiation in gamma ray bursts, and the mechanism can possibly explain more generally the origin of non-thermal radiation from shocked inter- and circum-stellar regions and from relativistic jets.

Using a three dimensional relativistic particle-in-cell code we have performed numerical experiments of plasma shells colliding at relativistic velocities. Such scenarios are found in many astrophysical objects e.g. the relativistic outflow from gamma ray bursts, active galactic nuclei jets and supernova remnants. We show how a Weibel-like two-stream instability is capable of generating small-scale magnetic filaments with strength up to percents of equipartition. Such field topology is ideal for the generation of jitter radiation as opposed to synchrotron radiation. We also explain how the field generating mechanism involves acceleration of electrons to power law distributions (N(E)~E^(-p))through a non-Fermi acceleration mechanism. The results add to our understanding of collisionless shocks.

We present results from three-dimensional particle simulations of collisionless shocks with relativistic counter-streaming ion-electron plasmas. Particles are followed over many skin depths downstream of the shock. Open boundaries allow the experiments to be continued for several particle crossing times. The experiments confirm the generation of strong magnetic and electric fields by a Weibel-like kinetic streaming instability, and demonstrate that the electromagnetic fields propagate far downstream of the shock. The magnetic fields are predominantly transversal, and are associated with merging ion current channels. The total magnetic energy grows as the ion channels merge, and as the magnetic field patterns propagate down stream. The electron populations are quickly thermalized, while the ion populations retain distinct bulk speeds in shielded ion channels and thermalize much more slowly. These results may help explain the origin of the magnetic fields responsible for afterglow synchrotron/jitter radiation from Gamma-Ray Bursts.

Resonant particle production, along with many other physical processes which change the effective equation of state (EOS) during inflation, introduces a step-like feature in the primordial power spectrum. We calculate observational constraints on resonant particle production, parameterised in form of an effective step height, N_eff and location in k-space, k_break. Combining data from the cosmic microwave background and the 2dF Galaxy Redshift Survey yields strong constraints in some regions of parameter space, although the range in k-space which can be probed is restricted to k ~ 0.001 - 0.1 h Mpc^-1. We also discuss the implications of our findings for general models which change the effective EOS during inflation.

We present numerical results from plasma particle simulations of collisionless shocks and ultra-relativistic counter-streaming plasmas. We demonstrate how the field-particle interactions lead to particle acceleration behind the shock-front. Further, we demonstrate how ultra relativistic counter-streaming plasmas create large scale patchy magnetic field structures and that these field structures propagate down-stream of the shock front. These results may help explain the origin of the magnetic fields and accelerated electrons responsible for afterglow synchrotron radiation from gamma ray bursts.