We propose a new constraint on light (sub-GeV) particles beyond the Standard Model that can be produced inside the proto-neutron star core resulting from the core-collapse supernova explosion. It is derived by demanding that the energy carried by exotic particles being transferred to the progenitor stellar envelopes must not exceed the explosion energy of $\lesssim 2\cdot 10^{51}$ erg of observed supernovae. We show specifically that for the case of a dark photon which kinetically mixes with the SM photon and decays predominantly to an $e^\pm$ pair, a smaller mixing parameter of one order of magnitude below the well-established supernova cooling bound can be excluded. Furthermore, our bound fills the gap between the cooling bound and the region constrained by (non)observation of gamma rays produced from supernovae for dark photons lighter than ~ 20 MeV. Our result also rules out the possibility of aiding successful supernova explosions by transferring energy from the supernova core to the shock with exotic particles.

When a neutron star cools to below the critical temperature for the onset of superfluidity, nucleon pair breaking and formation (PBF) processes become the dominant mechanism for neutrino emission, while the modified URCA and the nuclear bremsstrahlung processes are suppressed. The PBF processes in neutron stars have also been used to set upper limits on the properties of axions, which are comparable to those set by supernova SN 1987A. We apply this constraint on Weinberg's Higgs portal model, in which the dark radiation particles (the Goldstone bosons) and the dark matter candidate (a Majorana fermion) interact with the Standard Model (SM) fields solely through the mixing of the SM Higgs boson and a light Higgs boson. We compare the Goldstone boson emissivity with that of the neutrinos by considering several superfluid gap models for the neutron singlet-state pairing in the neutron star inner crust, as well as in the core region. We find that the PBF processes in the superfluid neutron star interior can indeed probe Weinberg's Higgs portal model in a new parameter space region. Together with our previous works on the constraints from supernovae and gamma-ray bursts, this study demonstrates further the competitiveness and complementarity of astrophysics to laboratory particle physics experiments.

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

The observed burst duration and energies of the neutrinos from Supernova 1987A strongly limit the possibility of any weakly-interacting light particle species being produced in the proto-neutron star (PNS) core and leading to efficient energy loss. We reexamine this constraint on Weinberg's Higgs portal model, in which the dark radiation particles (the Goldstone bosons) and the dark matter candidate (a Majorana fermion) interact with Standard Model (SM) fields solely through the mixing of the SM Higgs boson and a light Higgs boson. In order for the Goldstone bosons to freely stream out of the PNS core region, the Higgs portal coupling has to be about a factor of $4$--$9$ smaller than the current collider bound inferred from the SM Higgs invisible decay width. We find that in the energy loss rate calculations, results obtained by using the one-pion exchange (OPE) approximation and the SP07 global fits for the nucleon-nucleon total elastic cross section differ only by a factor $\lesssim 3$. The SN 1987A constraints surpass those set by laboratory experiments or by the energy loss arguments in other astrophysical objects such as the gamma-ray bursts, even with other nuclear uncertainties taken into account. Furthermore, the SN 1987A constraints are comparable to bounds from the latest dark matter direct search for low-mass WIMPs ($\lesssim 10~\Gev$.)

We report on a search for ultra-high-energy (UHE) neutrinos from gamma-ray bursts (GRBs) in the data set collected by the Testbed station of the Askaryan Radio Array (ARA) in 2011 and 2012. From 57 selected GRBs, we observed no events that survive our cuts, which is consistent with 0.12 expected background events. Using NeuCosmA as a numerical GRB reference emission model, we estimate upper limits on the prompt UHE GRB neutrino fluence and quasi-diffuse flux from $10^{7}$ to $10^{10}$ GeV. This is the first limit on the prompt UHE GRB neutrino quasi-diffuse flux above $10^{7}$ GeV.

We give the three-dimensional dynamical autonomous systems for most of the popular scalar field dark energy models including (phantom) quintessence, (phantom) tachyon, k-essence and general non-canonical scalar field models, change the dynamical variables from variables $(x, y, \lambda)$ to observable related variables $(w_{\phi}, \Omega_{\phi}, \lambda)$, and show the intimate relationships between those scalar fields that the three-dimensional system of k-essence can reduce to (phantom) tachyon, general non-canonical scalar field can reduce to (phantom) quintessence and k-essence can also reduce to (phantom) quintessence for some special cases. For the applications of the three-dimensional dynamical systems, we investigate several special cases and give the exactly dynamical solutions in detail. In the end of this paper, we argue that, it is more convenient and also has more physical meaning to express the differential equations of dynamical systems in $(w_{\phi}, \Omega_{\phi}, \lambda)$ instead of variables $(x, y, \lambda)$ and to investigate the dynamical system in 3-Dimension instead of 2-Dimension. We also raise a question about the possibility of the chaotic behavior in the spatially flat single scalar field FRW cosmological models in the presence of ordinary matter.

Gamma-ray bursts (GRBs) are the most energetic explosion events in the universe. An amount of gravitational energy of the order of the rest-mass energy of the Sun is released from a small region, within seconds or longer. This should lead to the formation of a fireball of temperature in the MeV range, consisting of electrons/positrons, photons, and a small fraction of baryons. We exploit the potential of GRB fireballs for being a laboratory for testing particle physics beyond the Standard Model, where we find that Weinberg's Higgs portal model serves as a good candidate for this purpose. Due to the resonance effects, the Goldstone bosons can be rapidly produced by electron-positron annihilation process in the initial fireballs of the gamma-ray bursts. On the other hand, the mean free path of the Goldstone bosons is larger than the size of the GRB initial fireballs, so they are not coupled to the GRB's relativistic flow and can lead to significant energy loss. Using generic values for the GRB initial fireball energy, temperature, radius, expansion rate, and baryon number density, we find that the GRB bounds on the parameters of Weinberg's Higgs portal model are indeed competitive to current laboratory constraints.

TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade.

$f(R)$ gravity is one of the simplest theories of modified gravity to explain the accelerated cosmic expansion. Although it is usually assumed that the quasi-Newtonian approach (a combination of the quasi-static approximation and sub-Hubble limit) for cosmic perturbations is good enough to describe the evolution of large scale structure in $f(R)$ models, some studies have suggested that this method is not valid for all $f(R)$ models. Here, we show that in the matter-dominated era, the pressure and shear equations alone, which can be recast into four first-order equations to solve for cosmological perturbations exactly, are sufficient to solve for the Newtonian potential, $\Psi$, and the curvature potential, $\Phi$. Based on these two equations, we are able to clarify how the exact linear perturbations fit into different limits. We find that the Compton length controls the quasi-static behaviours in $f(R)$ gravity. In addition, regardless the validity of quasi-static approximation, a strong version of the sub-Hubble limit alone is sufficient to reduce the exact linear perturbations in any viable $f(R)$ gravity to second order. Our findings disagree with some previous studies where we find little difference between our exact and quasi-Newtonian solutions even up to $k=10 c^{-1} \mathcal{H}_0$.

Unlike the random radial orientation distribution of field elliptical galaxies, galaxies in a cluster are expected to point preferentially towards the center of the cluster, as a result of the cluster's tidal force on its member galaxies. In this work an analytic model is formulated to simulate this effect. The deformation time scale of a galaxy in a cluster is usually much shorter than the time scale of change of the tidal force; the dynamical process of the tidal interaction within the galaxy can thus be ignored. An equilibrium shape of a galaxy is then assumed to be the surface of equipotential, which is the sum of the self-gravitational potential of the galaxy and the tidal potential of the cluster at this location. We use a Monte-Carlo method to calculate the radial orientation distribution of these galaxies, by assuming the NFW mass profile of the cluster and the initial ellipticity of field galaxies. The radial angles show a single peak distribution centered at zero. The Monte-Carlo simulations also show that a shift of the reference center from the real cluster center weakens the anisotropy of the radial angle distribution. Therefore, the expected radial alignment cannot be revealed if the distribution of spatial position angle is used instead of that of radial angle. The observed radial orientations of elliptical galaxies in cluster Abell~2744 are consistent with the simulated distribution.

The Askaryan Radio Array (ARA) is an ultra-high energy ($>10^{17}$ eV) cosmic neutrino detector in phased construction near the South Pole. ARA searches for radio Cherenkov emission from particle cascades induced by neutrino interactions in the ice using radio frequency antennas ($\sim150-800$ MHz) deployed at a design depth of 200 m in the Antarctic ice. A prototype ARA Testbed station was deployed at $\sim30$ m depth in the 2010-2011 season and the first three full ARA stations were deployed in the 2011-2012 and 2012-2013 seasons. We present the first neutrino search with ARA using data taken in 2011 and 2012 with the ARA Testbed and the resulting constraints on the neutrino flux from $10^{17}-10^{21}$ eV.

The LISA Pathfinder mission will demonstrate the technology of drag-free test masses for use as inertial references in future space-based gravitational wave detectors. To accomplish this, the Pathfinder spacecraft will perform drag-free flight about a test mass while measuring the acceleration of this primary test mass relative to a second reference test mass. Because the reference test mass is contained within the same spacecraft, it is necessary to apply forces on it to maintain its position and attitude relative to the spacecraft. These forces are a potential source of acceleration noise in the LISA Pathfinder system that are not present in the full LISA configuration. While LISA Pathfinder has been designed to meet it's primary mission requirements in the presence of this noise, recent estimates suggest that the on-orbit performance may be limited by this `suspension noise'. The drift-mode or free-flight experiments provide an opportunity to mitigate this noise source and further characterize the underlying disturbances that are of interest to the designers of LISA-like instruments. This article provides a high-level overview of these experiments and the methods under development to analyze the resulting data.

Recently, Weinberg proposed a Higgs portal model with a spontaneously broken global $U(1)$ symmetry in which Goldstone bosons may be masquerading as fractional cosmic neutrinos. We extend the model by gauging the $U(1)$ symmetry. This gives rise to the so-called dark photon and dark Higgs. The dark photons can constitute about 0.912 (0.167) to the effective number of light neutrino species if they decouple from the thermal bath before the pions become non-relativistic and after (before) the QCD transition. Restriction on the parameter space of the portal coupling and the dark Higgs mass is obtained from the freeze-out condition of the dark photons. Combining with the collider data constraints on the invisible width of the standard model Higgs requires the dark Higgs mass to be less than a few GeV.

We give a full investigation on the dynamics of power-law kinetic quintessence $L(X, \phi)=V(\phi)(-X+X^2)$ by considering the potential related parameter $\Gamma$($=\frac{V V''}{V'^2}$) as a function of another potential parameter $\lambda$($=\frac{V'}{\kappa V^{3/2}}$), which correspondingly extends the analysis of the dynamical system of our universe from two-dimension to three-dimension. Beside the critical points found in previous papers, we find a new de-Sitter-like dominant attractor(\textit{cp$6$}) and give its stable condition using the center manifold theorem. For the dark energy dominant solution(\textit{cp$6$} and \textit{cp$7$}), it could be distinguished from canonical quintessence and tachyon models since the sound speed $c_s^2=0$ or $c_s^2\ll 1$. For the scaling solution (\textit{cp$8$}), it is very interesting that the sound speed $c_s^2=1/5$ while it behaves as ordinary matter. We therefore point out that the power-law kinetic quintessence should have different signatures on cold dark matter power spectrum and cosmic microwave background both at early time when this scalar field is an early dark energy with $\Omega_\phi$ being non-negligible at high redshift and at late time when it drives the accelerating expansion. We even do not know whether there are any degeneracies of the impacts between these two epoches. They are expected to be investigated in future.

Recently, Weinberg proposed a scenario where Goldstone bosons may be masquerading as fractional cosmic neutrinos. We calculate the energy loss rates through the emission of these Goldstone bosons in a post-collapse supernova core. Invoking the well established emissivity bound from the Supernova 1987A observations and simulations, we find that nuclear bremsstrahlung processes can notably impose a bound on the Goldstone boson coupling to the Standard Model Higgs, $g$, dependent on the mass of the associated radial field, $m_r$. For $m_r$ large enough compared with the temperature in the post-collapse supernova core, our bound is $|g| \lesssim 0.011\, (m_r / 500~{\rm MeV})^2$, very competitive to that derived from collider experiments.

LISA Pathfinder is a mission planned by the European Space Agency to test the key technologies that will allow the detection of gravitational waves in space. The instrument on-board, the LISA Technology package, will undergo an exhaustive campaign of calibrations and noise characterisation campaigns in order to fully describe the noise model. Data analysis plays an important role in the mission and for that reason the data analysis team has been developing a toolbox which contains all the functionalities required during operations. In this contribution we give an overview of recent activities, focusing on the improvements in the modelling of the instrument and in the data analysis campaigns performed both with real and simulated data.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

[ABRIDGED] The new generation of wide field optical imaging like the Canada France Hawaii Telescope Legacy Survey (CFHTLS) enables discoveries of all types of gravitational lenses present in the sky. The Strong Lensing Legacy Survey (SL2S) project has started an inventory, respectively for clusters or groups of galaxies lenses, and for Einstein rings around distant massive ellipticals. Here we attempt to extend this inventory by finding lensing events produced by massive edge-on disk galaxies which remains a poorly documented class of lenses. We implement and test an automated search procedure of edge-on galaxy lenses in the CFHTLS Wide fields with magnitude 18

Can dark matter be stabilized by charge conservation, just as the electron is in the standard model? We examine the possibility that dark matter is hidden, that is, neutral under all standard model gauge interactions, but charged under an exact U(1) gauge symmetry of the hidden sector. Such candidates are predicted in WIMPless models, supersymmetric models in which hidden dark matter has the desired thermal relic density for a wide range of masses. Hidden charged dark matter has many novel properties not shared by neutral dark matter: (1) bound state formation and Sommerfeld-enhanced annihilation after chemical freeze out may reduce its relic density, (2) similar effects greatly enhance dark matter annihilation in protohalos at redshifts of z ~ 30, (3) Compton scattering off hidden photons delays kinetic decoupling, suppressing small scale structure, and (4) Rutherford scattering makes such dark matter self-interacting and collisional, potentially impacting properties of the Bullet Cluster and the observed morphology of galactic halos. We analyze all of these effects in a WIMPless model in which the hidden sector is a simplified version of the minimal supersymmetric standard model and the dark matter is a hidden sector stau. We find that charged hidden dark matter is viable and consistent with the correct relic density for reasonable model parameters and dark matter masses in the range 1 GeV < m_X < 10 TeV. At the same time, in the preferred range of parameters, this model predicts cores in the dark matter halos of small galaxies and other halo properties that may be within the reach of future observations. These models therefore provide a viable and well-motivated framework for collisional dark matter with Sommerfeld enhancement, with novel implications for astrophysics and dark matter searches.

The existence of strong lensing systems with Einstein radii (Re) covering the full mass spectrum, from ~1-2" (produced by galaxy scale dark matter haloes) to >10" (produced by galaxy cluster scale haloes) have long been predicted. Many lenses with Re around 1-2" and above 10" have been reported but very few in between. In this article, we present a sample of 13 strong lensing systems with Re in the range 3"- 8", i.e. systems produced by galaxy group scale dark matter haloes, spanning a redshift range from 0.3 to 0.8. This opens a new window of exploration in the mass spectrum, around 10^{13}- 10^{14} M_{sun}, which is a crucial range for understanding the transition between galaxies and galaxy clusters. Our analysis is based on multi-colour CFHTLS images complemented with HST imaging and ground based spectroscopy. Large scale properties are derived from both the light distribution of the elliptical galaxies group members and weak lensing of the faint background galaxy population. On small scales, the strong lensing analysis yields Einstein radii between 2.5" and 8". On larger scales, the strong lenses coincide with the peak of the light distribution, suggesting that mass is traced by light. Most of the luminosity maps have complicated shapes, indicating that these intermediate mass structures are dynamically young. Fitting the reduced shear with a Singular Isothermal Sphere, we find sigma ~ 500 km/s and an upper limit of ~900 km/s for the whole sample. The mass to light ratio for the sample is found to be M/L_i ~ 250 (solar units, corrected for evolution), with an upper limit of 500. This can be compared to mass to light ratios of small groups (with sigma ~ 300 km/s and galaxy clusters with sigma > 1000 km/s, thus bridging the gap between these mass scales.

SL2SJ02176-0513 is a remarkable lens for the presence of two multiply-imaged systems at different redshifts lensed by a foreground massive galaxy at $z_{\rm lens}=0.656$: a bright cusp arc at $z_{\rm arc}=1.847$ and an additional double-image system at an estimated redshift of $z_{\rm dbl}\sim2.9$ based on photometry and lensing geometry. The system is located about 400 kpc away from the center of a massive group of galaxies. Mass estimates for the group are available from X-ray observations and satellite kinematics. Multicolor photometry provides an estimate of the stellar mass of the main lens galaxy. The lensing galaxy is modeled with two components (stars and dark matter), and we include the perturbing effect of the group environment, and all available constraints. We find that classic lensing degeneracies, e.g. between external convergence and mass density slope, are significantly reduced with respect to standard systems and infer tight constraints on the mass density profile: (i) the dark matter content of the main lens galaxy is in line with that of typical galaxies $f_{\rm dm}(<R_{\rm e})=0.41^{+0.09}_{-0.06}$; (ii) the required mass associated with the dark matter halo of the nearby group is consistent with X-ray and weak-lensing estimates ($\sigma_{\rm grp}=550^{+130}_{-240}$); (iii) accounting for the group contribution in the form of an external convergence, the slope of the mass density profile of the main lens galaxy alone is found to be $\alpha=-1.03^{+0.22}_{-0.16}$, consistent with the isothermal ($\alpha=-1$) slope. We demonstrate that multiple source plane systems together with good ancillary dataset can be used to disentangle local and environmental effects.

Dark matter may be hidden, with no standard model gauge interactions. At the same time, in WIMPless models with hidden matter masses proportional to hidden gauge couplings squared, the hidden dark matter's thermal relic density may naturally be in the right range, preserving the key quantitative virtue of WIMPs. We consider this possibility in detail. We first determine model-independent constraints on hidden sectors from Big Bang nucleosynthesis and the cosmic microwave background. Contrary to conventional wisdom, large hidden sectors are easily accommodated. A flavour-free version of the standard model is allowed if the hidden sector is just 30% colder than the observable sector after reheating. Alternatively, if the hidden sector contains a 1-generation version of the standard model with characteristic mass scale below 1 MeV, even identical reheating temperatures are allowed. We then analyze hidden sector freezeout in detail for a concrete model, solving the Boltzmann equation numerically and understanding the results from both observable and hidden sector points of view. We find that WIMPless dark matter indeed obtains the correct relic density for masses in the range keV < m_X < TeV. The upper bound results from the requirement of perturbativity, and the lower bound assumes that the observable and hidden sectors reheat to the same temperature and is raised to the MeV scale if the hidden sector is 10 times colder. WIMPless dark matter therefore generalizes the WIMP paradigm to the largest mass range possible for viable thermal relics and provides a unified framework for exploring dark matter signals across nine orders of magnitude in dark matter mass.

In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' formed by 4 bright images at z_{spec}=0.888 located very close to the cD galaxy, providing observational constraints that are potentially very interesting to probe the central mass distribution. The stellar contribution from the cD galaxy (~1.25 10^{12} M_{sun} within 7") is accounted for in our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized NFW profile parametrized with a central logarithmic slope \alpha. We find that within the range where observational constraints are present (from ~5" to ~50"), the slope of the dark matter distribution in Abell1703 is equal to 1.09^{+0.05}_{-0.11} (3\sigma confidence level). The concentration parameter is equal to c_{200} ~ 3.5, and the scale radius is constrained to be larger than the region where observational constraints are available. Within this radius, the 2D mass is equal to M(50")=2.4 10^{14} M_{\sun}. We cannot draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing baryons to make a proper comparison. We advocate the need for a sample of observed and simulated unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models.

The strong lensing modelling of gravitational ``rings'' formed around massive galaxies is sensitive to the amplitude of the external shear and convergence produced by nearby mass condensations. In current wide field surveys, it is now possible to find out a large number of rings, typically 10 gravitational rings per square degree. We propose here, to systematically study gravitational rings around galaxy clusters to probe the cluster mass profile beyond the cluster strong lensing regions. For cluster of galaxies with multiple arc systems, we show that rings found at various distances from the cluster centre can improve the modelling by constraining the slope of the cluster mass profile. We outline the principle of the method with simple numerical simulations and we apply it to 3 rings discovered recently in Abell~1689. In particular, the lens modelling of the 3 rings confirms that the cluster is bimodal, and favours a slope of the mass profile steeper than isothermal at a cluster radius $\sim 300 \kpc$. These results are compared with previous lens modelling of Abell~1689 including weak lensing analysis. Because of the difficulty arising from the complex mass distribution in Abell~1689, we argue that the ring method will be better implemented on simpler and relaxed clusters.

If the diffuse extragalactic gamma ray emission traces the large scale structures of the universe, peculiar anisotropy patterns are expected in the gamma ray sky. In particular, because of the cutoff distance introduced by the absorption of 0.1-10 TeV photons on the infrared/optical background, prominent correlations with the local structures within a range of few hundreds Mpc should be present. We provide detailed predictions of the signal based on the PSCz map of the local universe. We also use mock N-body catalogues complemented with the halo model of structures to study some statistical features of the expected signatures. The results are largely independent from cosmological details, and depend mostly on the index of correlation (or bias) of the sources with respect to the large scale distribution of galaxies. For instance, the predicted signal in the case of a quadratic correlation (as it may happen for a dark matter annihilation contribution to the diffuse gamma flux) differs substantially from a linear correlation case, providing a complementary tool to unveil the nature of the sources of the diffuse gamma ray emission. The chances of the present and future space and ground based observatories to measure these features are discussed.

Surveys of weak gravitational lensing of distant galaxies will be one of the key cosmological probes in the future. We study the ability of such surveys to constrain neutrino masses and the equation of state parameter of the dark energy, focussing on how tomographic information can improve the sensitivity to these parameters. We also provide a detailed discussion of systematic effects pertinent to weak lensing surveys, and the possible degradation of sensitivity to cosmological parameters due to these effects. For future probes such as the Large Synoptic Survey Telescope survey, we find that, when combined with cosmic microwave background data from the Planck satellite, a sensitivity to neutrino masses of sigma(sum m_nu) < 0.05 eV can be reached. This results is robust against variations in the running of the scalar spectral index, the time-dependence of dark energy equation of state, and/or the number of relativistic degrees of freedom.

The Fisher matrix formalism has in recent times become the standard method for predicting the precision with which various cosmological parameters can be extracted from future data. This approach is fast, and generally returns accurate estimates for the parameter errors when the individual parameter likelihoods approximate a Gaussian distribution. However, where Gaussianity is not respected (due, for instance, to strong parameter degeneracies), the Fisher matrix formalism loses its reliability. In this paper, we compare the results of the Fisher matrix approach with those from Monte Carlo simulations. The latter method is based on the publicly available CosmoMC code, but uses synthetic realisations of data sets anticipated for future experiments. We focus on prospective cosmic microwave background (CMB) data from the Planck satellite, with or without CMB lensing information, and its implications for a minimal cosmological scenario with eight parameters and an extended model with eleven parameters. We show that in many cases, the projected sensitivities from the Fisher matrix and the Monte Carlo methods differ significantly, particularly in models with many parameters. Sensitivities to the neutrino mass and the dark matter fraction are especially susceptible to change.

We have studied bounds on the neutrino mass using new data from the WMAP 3 year data, the Sloan Digital Sky Survey measurement of the baryon acoustic peak, the Type Ia supernovae from SNLS, and the Lyman-alpha forest. We find that even in the most general models with a running spectral index where the number of neutrinos and the dark energy equation of state are allowed to vary, the 95% C.L. bound on the sum of neutrino masses is sum m_nu < 0.62 eV (95% C.L.), a bound which we believe to be robust. In the more often used constrained analysis with N_nu =3, w = -1, and alpha_s = 0, we find a bound of 0.48 eV without using the Lyman-alpha data. If Lyman-alpha data is used, the bound shrinks to \sum m_nu < 0.2-0.4 eV (95% C.L.), depending strongly on the Lyman-alpha analysis used.

Scattering processes in the cosmic microwave background limit the propagation of ultra high energy charged particles in our Universe. For extragalactic proton sources resonant photopion production results in the famous Greisen-Zatsepin-Kuzmin (GZK) cutoff at about 4x10^10 GeV expected in the spectrum observed on Earth. The faint flux of ultra high energy cosmic rays of less than one event per year and cubic kilometer and the large systematic uncertainties in the energy calibration of cosmic ray showers is a challenge for cosmic ray observatories and so far the GZK cutoff has not been unambiguously confirmed. We have investigated the possibility that the primaries of super-GZK events are strongly interacting neutrinos which are not subject to the GZK cutoff. For the flux of protons and neutrinos from extragalactic optically thin sources and a flexible parameterization of the neutrino-nucleon cross section we have analyzed the cosmic ray spectra observed at AGASA and HiRes taking also into account results from horizontal events at AGASA and contained events at RICE. We find that scenarios of strongly interacting neutrinos are still compatible with the data requiring a steep increase of the inelastic neutrino-nucleon cross section by four order of magnitude within one energy decade compared to the Standard Model predictions. We also discuss the impact of the preliminary cosmic ray spectrum observed by the Pierre Auger Observatory.

We study the difference between thermally produced fermionic and bosonic hot dark matter in detail. In the linear regime of structure formation, their distinct free-streaming behaviours can lead to pronounced differences in the matter power spectrum. While not detectable with current cosmological data, such differences will be clearly observable with upcoming large scale weak lensing surveys for particles as light as m_HDM ~ 0.2 eV. In the nonlinear regime, bosonic hot dark matter is not subject to the same phase space constraints that severely limit the amount of fermionic hot dark matter infall into cold dark matter halos. Consequently, the overdensities in fermionic and bosonic hot dark matter of equal particle mass can differ by more than a factor of five in the central part of a halo. However, this unique manifestation of quantum statistics may prove very difficult to detect unless the mass of the hot dark matter particle and its decoupling temperature fall within a very narrow window, 1 < m_HDM/eV < 4 and g* < 30. In this case, hot dark matter infall may have some observable consequences for the nonlinear power spectrum and hence the weak lensing convergence power spectrum at l ~ 10^3 - 10^4 at the percent level.

Resonant photopion production with the cosmic microwave background predicts a suppression of extragalactic protons above the famous Greisen-Zatsepin-Kuzmin cutoff at about E_{GZK} ~ 5 x 10^10 GeV. Current cosmic ray data measured by the AGASA and HiRes Collaborations do not unambiguously confirm the GZK cutoff and leave a window for speculations about the origin and chemical composition of the highest energy cosmic rays. In this work we analyze the possibility of strongly interacting neutrino primaries and derive model-independent quantitative requirements on the neutrino-nucleon inelastic cross section for a viable explanation of the cosmic ray data. Search results on weakly interacting cosmic particles from the AGASA and RICE experiments are taken into account simultaneously. Using a flexible parameterization of the inelastic neutrino-nucleon cross section we find that a combined fit of the data does not favor the Standard Model neutrino-nucleon inelastic cross section, but requires, at 90 % confidence level, a steep increase within one energy decade around E_{GZK} by four orders of magnitude. We illustrate such an enhancement within some extensions of the Standard Model. The impact of new cosmic ray data or cosmic neutrino search results on this scenario, notably from the Pierre Auger Observatory soon, can be immediately evaluated within our approach.

Extraterrestrial neutrinos can initiate deeply developing air showers, and those that traverse the atmosphere unscathed may produce cascades in the ice or water. Up to now, no such events have been observed. This can be translated into upper limits on the diffuse neutrino flux. On the other hand, the observation of cosmic rays with primary energies > 10^{10} GeV suggests that there is a guaranteed flux of cosmogenic neutrinos, arising from the decay of charged pions (and their muon daughters) produced in proton interactions with the cosmic microwave background. In this work, armed with these cosmogenic neutrinos and the increased exposure of neutrino telescopes we bring up-to-date model-independent upper bounds on the neutrino-nucleon inelastic cross section. Uncertainties in the cosmogenic neutrino flux are discussed and taken into account in our analysis. The prospects for improving these bounds with the Pierre Auger Observatory are also estimated. The unprecedented statistics to be collected by this experiment in 6 yr of operation will probe the neutrino-nucleon inelastic cross section at the level of Standard Model predictions.

After a short review of the ultrahigh energy cosmic ray puzzle - the apparent observation of cosmic rays originating from cosmological distances with energies above the expected Greisen-Zatsepin-Kuzmin cutoff 4x10^{19} eV - we consider strongly interacting neutrino scenarios as an especially interesting solution. We show that all features of the ultrahigh energy cosmic ray spectrum from 10^{17} eV to 10^{21} eV can be described to originate from a simple power-like injection spectrum of protons, under the assumption that the neutrino-nucleon cross-section is significantly enhanced at center of mass energies above ~ 100 TeV. In such a scenario, the cosmogenic neutrinos produced during the propagation of protons through the cosmic microwave background initiate air showers in the atmosphere, just as the protons. The total air shower spectrum induced by protons and neutrinos shows excellent agreement with the observations. We shortly discuss TeV-scale extensions of the Standard Model which may lead to a realization of a strongly interacting neutrino scenario. We emphasize, however, that such a scenario may even be realized within the standard electroweak model: electroweak instanton/sphaleron induced processes may get strong at ultrahigh energies. Possible tests of strongly interacting neutrino scenarios range from observations at cosmic ray facilities and neutrino telescopes to searches at lepton nucleon scattering experiments.

We show that all features of the ultrahigh energy cosmic ray spectrum from 10^{17} eV to 10^{21} eV can be described with a simple power-like injection spectrum of protons under the assumption that the neutrino-nucleon cross-section is significantly enhanced at center of mass energies above \approx 100 TeV. In our scenario, the cosmogenic neutrinos produced during the propagation of protons through the cosmic microwave background initiate air showers in the atmosphere, just as the protons. The total air shower spectrum induced by protons and neutrinos shows excellent agreement with the observations. A particular possibility for a large neutrino-nucleon cross-section exists within the Standard Model through electroweak instanton-induced processes.

Under the assumption that some part of the observed highest energy cosmic rays consists of protons originating from cosmological distances, we derive bounds on the associated flux of neutrinos generated by inelastic processes with the cosmic microwave background photons. We exploit two methods. First, a power-like injection spectrum is assumed. Then, a model-independent technique, based on the inversion of the observed proton flux, is presented. The inferred lower bound is quite robust. As expected, the upper bound depends on the unknown composition of the highest energy cosmic rays. Our results represent benchmarks for all ultrahigh energy neutrino telescopes.

We propose a scenario in which a simple power-like primary spectrum for protons with sources at cosmological distances leads to a quantitative description of all the details of the observed cosmic ray spectrum for energies from 10^{17} eV to 10^{21} eV. As usual, the ultrahigh energy protons with energies above E_{GZK} ~ 4 x 10^{19} eV loose a large fraction of their energies by the photoproduction of pions on the cosmic microwave background, which finally decay mainly into neutrinos. In our scenario, these so-called cosmogenic neutrinos interact with nucleons in the atmosphere through Standard Model electroweak instanton-induced processes and produce air showers which are hardly distinguishable from ordinary hadron-initiated air showers. In this way, they give rise to a second contribution to the observed cosmic ray spectrum -- in addition to the one from above mentioned protons -- which reaches beyond E_{GZK}. Since the whole observed spectrum is uniquely determined by a single primary injection spectrum, no fine tuning is needed to fix the ratio of the spectra below and above E_{GZK}. The statistical analysis shows an excellent goodness of this scenario. Possible tests of it range from observations at cosmic ray facilities and neutrino telescopes to searches for QCD instanton-induced processes at HERA.

In scenarios with extra dimensions and TeV-scale quantum gravity, black holes are expected to be produced copiously at center-of-mass energies above the fundamental Planck scale. The Large Hadron Collider (LHC) may thus turn into a factory of black holes, at which their production and evaporation may be studied in detail. But even before the LHC starts operating, the Pierre Auger Observatory for cosmic rays, presently under construction, has an opportunity to search for black hole signatures. Black hole production in the scattering of ultrahigh energy cosmic neutrinos on nucleons in the atmosphere may initiate quasi-horizontal air showers far above the Standard Model rate. In this letter, we compare the sensitivity of LHC and Auger to black hole production by studying their respective reach in black hole production parameter space. Moreover, we present constraints in this parameter space from the non-observation of horizontal showers by the Fly's Eye collaboration. We find that if the ultrahigh energy neutrino flux is at the level expected from cosmic ray interactions with the cosmic microwave background radiation, Auger has only a small window of opportunity to detect black holes before the start of the LHC. If, on the other hand, larger ultrahigh energy neutrino fluxes on the level of the upper limit from ``hidden'' hadronic astrophysical sources are realized in nature, then the first signs of black hole production may be observed at Auger. Moreover, in this case, the Fly's Eye constraints, although more model dependent, turn out to be competitive with other currently available constraints on TeV-scale gravity which are mainly based on interactions associated with Kaluza-Klein gravitons.