Numerical simulations of massive neutrino cosmologies consistently find a spoon-like feature in the non-linear matter power spectrum ratios of cosmological models that differ only in the neutrino mass fraction f_N. Typically, the ratio approaches unity at low wave numbers k, decreases by ~ 10 f_N at k ~ 1 h/Mpc, and turns up again at large k. Using the halo model of large-scale structure, we show that this spoon feature originates in the transition from the two-halo power spectrum to the one-halo power spectrum. The former's sensitivity to f_N rises with k, while that of the latter decreases with k. The presence of this spoon feature is robust with respect to different choices of the halo mass function and the halo density profile, and does not require any parameter tuning within the halo model. We demonstrate that a standard halo model calculation is already able to predict the depth, width, and position of this spoon as well as its evolution with redshift z with remarkable accuracy. Predictions at z >= 1 can be further improved using non-linear perturbative inputs.

Investigating the Gunn-Peterson trough of high redshift quasars (QSOs) is a powerful way to reveal the cosmic reionization. As one of such attempts, we perform a series of analyses to examine the absorption lines observed with one of the highest redshift QSOs, PSO J006.1240+39.2219, which we previously discovered at z = 6.62. Using the Subaru telescope, we obtained medium-resolution spectrum with a total exposure time of 7.5 hours. We calculate the Ly$\alpha$ transmission in different redshift bins to determine the near zone radius and the optical depth at 5.6$<$z$<$6.5. We find a sudden change in the Ly$\alpha$ transmission at 5.75$<$z$<$5.86, which is consistent with the result from the literature. The near zone radius of the QSO is 5.79$\pm$0.09 $p$Mpc, within the scatter of the near zone radii of other QSOs measured in previous studies. We also analyze the dark gap distribution to probe the neutral hydrogen fractions beyond the saturation limit of the Gunn-Peterson trough. We extend the measurement of the dark gaps to 5.7$<$z$<$6.3. We find that the gap widths increase with increasing redshifts, suggesting more neutral Universe at higher redshifts. However, these measurements strongly depend on the continuum modeling. As a continuum model-free attempt, we also perform the dark-pixel counting analysis, to find the upper limit of $\langle x_{\rm H I} \rangle \sim$0.6 (0.8) at $z<$5.8 ($z>$5.8). All three analyses based on this QSO show increasingly neutral hydrogen towards higher redshifts, adding precious measurements up to z$\sim$6.5.

We propose an alternative approach to the construction of fitting functions to the nonlinear matter power spectrum extracted from $N$-body simulations based on the relative matter power spectrum $\delta(k,a)$, defined as the fractional deviation in the absolute matter power spectrum produced by a target cosmology away from a reference $\Lambda$CDM prediction. From the computational perspective, $\delta(k,a)$ is fairly insensitive to the specifics of the simulation settings, and numerical convergence at the 1%-level can be readily achieved without the need for huge computing capacity. Furthermore, $\delta(k,a)$ exhibits several interesting properties that enable a piece-wise construction of the full fitting function, whereby component fitting functions are sought for single-parameter variations and then multiplied together to form the final product. Then, to obtain 1%-accurate absolute power spectrum predictions for any target cosmology only requires that the community as a whole invests in producing one single ultra-precise reference $\Lambda$CDM absolute power spectrum, to be combined with the fitting function to produce the desired result. To illustrate the power of this approach, we have constructed the fitting function RelFit using only five relatively inexpensive $w$CDM simulations (box length $L=256 h^{-1}$Mpc, $N=1024^3$ particles, initialised at $z_i=49$). In a 6-parameter space spanning $\{\omega_m,A_s,n_s,w,\omega_b,h\}$, the output relative power spectra of RelFit are consistent with the predictions of the CosmicEmu emulator to 1% or better for a wide range of cosmologies up to $k\simeq 10$/Mpc. Thus, our approach could provide an inexpensive and democratically accessible route to fulfilling the 1%-level accuracy demands of the upcoming generation of large-scale structure probes, especially in the exploration of "non-standard" or "exotic" cosmologies on nonlinear scales.

We revisit several aspects of Standard Model physics at finite temperature that drive the theoretical value of the cosmological parameter $N_{\rm eff}$, the effective number of neutrinos in the early universe, away from 3. Our chief focus is finite-temperature corrections to the equation of state of the QED plasma in the vicinity of neutrino decoupling at $T \sim 1$ MeV, where $T$ is the photon temperature. Working in the instantaneous decoupling approximation, we recover at ${\cal O}(e^2)$, where $e$ is the elementary electric charge, the well-established correction of $\delta N_{\rm eff}^{(2)} \simeq 0.010$ across a range of plausible neutrino decoupling temperatures, in contrast to an erroneous claim in the recent literature which found twice as large an effect. At ${\cal O}(e^3)$ we find a new and significant correction of $\delta N_{\rm eff}^{(3)} \simeq -0.001$ that has so far not been accounted for in any precision calculation of $N_{\rm eff}$, significant because this correction is potentially larger than the change in $N_{\rm eff}$ induced between including and excluding neutrino oscillations in the transport modelling. In addition to the QED equation of state, we make a first pass at quantifying finite-temperature QED corrections to the weak interaction rates that directly affect the neutrino decoupling process, and find that the ${\cal O}(e^2)$ thermal electron mass correction induces a change of $\delta N_{\rm eff}^{m_{\rm th}} \lesssim 10^{-4}$. A complete assessment of the various effects considered in this work on the final value of $N_{\rm eff}$ will necessitate an account of neutrino energy transport beyond the instantaneous decoupling approximation. However, relative to $N_{\rm eff} = 3.044$ obtained in the most recent such calculation, we expect the new effects found in this work to lower the number to $N_{\rm eff} = 3.043$.

We explore the issue of prior dependence in the context of one-sided constraints on the dark matter-photon and dark matter-neutrino elastic scattering cross-sections derived from cosmic microwave background (CMB) anisotropies measurements. Testing in particular the linear flat, Jeffreys, and logarithmic flat priors, we find that the former two yield upper limits on the cross-sections that are mutually consistent to within 20%. In contrast, bounds derived under the assumption of the logarithmic flat prior are strongly sensitive to the choice of the lower prior boundary. Indeed, surveying the recent literature, we find that this pathology of the logarithmic prior has resulted in published constraints that are up to an order of magnitude artificially tighter than they should objectively be. Our revised `objective' constraints from the 2015 data of the Planck CMB mission on the present-day scattering cross-sections are $\sigma_{\rm DM-\gamma} < 1.72 \times 10^{-6} \, \sigma_{\rm T} \, (m_{\rm DM}/{\rm GeV})$ and $\sigma_{\rm DM-\gamma} < 2.74 \times 10^{-15} \, \sigma_{\rm T} \, (m_{\rm DM}/{\rm GeV})$ for dark matter-photon interactions scaling as $a^0$ and $a^{-2}$ respectively, where $a$ is the scale factor, and $\sigma_{\rm T}$ the total Thomson scattering cross-section. Their dark matter-neutrino counterparts read $\sigma_{\rm DM-\nu} < 2.14 \times 10^{-6} \, \sigma_{\rm T} \, (m_{\rm DM}/{\rm GeV})$ and $\sigma_{\rm DM-\nu}< 2.46 \times 10^{-15} \, \sigma_{\rm T} \, (m_{\rm DM}/{\rm GeV})$. All have been computed assuming the Jeffreys prior.

Elastic scattering between dark matter particles and a relativistic species such as photons or neutrinos leads to a transfer of energy from the latter due to their intrinsically different temperature scaling relations. In this work, we point out that this siphoning of energy from the radiation bath manifests as a change in the effective number of neutrinos $N_{\rm eff}$, and compute the expected shift $\Delta N_{\rm eff}$ for dark matter-photon and dark matter-neutrino elastic scattering as a function of the dark matter mass $m_\psi$ and scattering cross section $\sigma_{\psi-X}$. For $(m_\psi,\sigma_{\psi-X})$-parameter regions already explored by nonlinear probes such as the Lyman-$\alpha$ forest through collisional and/or free-streaming damping, we find shifts of $|\Delta N_{\rm eff}| \simeq O(10^{-2})$, which may be within the reach of the proposed CMB-S4 experiment. For most of the as-yet-unexplored parameter space, however, we expect $|\Delta N_{\rm eff}| \lesssim O(10^{-3})$. An ideal 21 cm tomography survey of the dark ages limited only by cosmic variance is potentially sensitive to $|\Delta N_{\rm eff}| \simeq O(10^{-6})$, in which case dark matter masses up to $m_{\psi} \sim 10 \, \textrm{MeV}$ may be probed via their effect on $N_{\rm eff}$.

We consider the impact of neutrino self-interactions described by an effective four-fermion coupling on cosmological observations. Implementing the exact Boltzmann hierarchy for interacting neutrinos first derived in [arxiv:1409.1577] into the Boltzmann solver CLASS, we perform a detailed numerical analysis of the effects of the interaction on the cosmic microwave background (CMB) anisotropies, and compare our results with known approximations in the literature. While we find good agreement between our exact approach and the relaxation time approximation used in some recent studies, the popular $\left( c_{\text{eff}}^2,c_{\text{vis}}^2 \right)$-parameterisation fails to reproduce the correct scale dependence of the CMB temperature power spectrum. We then proceed to derive constraints on the effective coupling constant $G_{\text{eff}}$ using currently available cosmological data via an MCMC analysis. Interestingly, our results reveal a bimodal posterior distribution, where one mode represents the standard $\Lambda$CDM limit with $G_{\rm eff} \lesssim 10^8 \, G_{\rm F}$, and the other a scenario in which neutrinos self-interact with an effective coupling constant $G_{\rm eff} \simeq 3 \times 10^9 \, G_{\rm F}$.

In this work we consider CMB spectral distortions as a probe of dark matter microphysics in the early universe. We demonstrate that future experiments such as PRISM have the potential to distinguish between scenarios which offer solutions to the small-scale problems of CDM cosmology.

The dissipation of small-scale perturbations in the early universe produces a distortion in the blackbody spectrum of cosmic microwave background photons. In this work, we propose to use these distortions as a probe of the microphysics of dark matter on scales $1\,\textrm{Mpc}^{-1}\lesssim k \lesssim 10^{4}\,\textrm{Mpc}^{-1}$. We consider in particular models in which the dark matter is kinetically coupled to either neutrinos or photons until shortly before recombination, and compute the photon heating rate and the resultant $\mu$-distortion in both cases. We show that the $\mu$-parameter is generally enhanced relative to $\Lambda$CDM for interactions with neutrinos, and may be either enhanced or suppressed in the case of interactions with photons. The deviations from the $\Lambda$CDM signal are potentially within the sensitivity reach of a PRISM-like experiment if $\sigma_{\textrm{DM}-\gamma} \gtrsim 1.1\times10^{-30} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $\sigma_{\textrm{DM}-\nu} \gtrsim 4.8\times 10^{-32} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for time-independent cross sections, and $\sigma^{0}_{\textrm{DM}-\gamma} \gtrsim 1.8 \times 10^{-40} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ and $\sigma^{0}_{\textrm{DM}-\nu} \gtrsim 2.5 \times 10^{-47} \left(m_{\textrm{DM}}/\textrm{GeV}\right) \textrm{cm}^{2}$ for cross sections scaling as temperature squared, coinciding with the parameter regions in which late kinetic decoupling may serve as a solution to the small-scale crisis. Furthermore, these $\mu$-distortion signals differ from those of warm dark matter (no deviation from $\Lambda$CDM) and a suppressed primordial power spectrum (strongly suppressed or a negative $\mu$-parameter), demonstrating that CMB spectral distortion can potentially be used to distinguish between solutions to the small-scale crisis.

Grain growth during star formation affects the physical and chemical processes in the evolution of star-forming clouds. We investigate the origin of the millimeter (mm)-sized grains recently observed in Class I protostellar envelopes. We use the coagulation model developed in our previous paper and find that a hydrogen number density of as high as $10^{10}~{\rm cm^{-3}}$, instead of the typical density $10^5~{\rm cm^{-3}}$, is necessary for the formation of mm-sized grains. Thus, we test a hypothesis that such large grains are transported to the envelope from the inner, denser parts, finding that gas drag by outflow efficiently "launches" the large grains as long as the central object has not grown to $\gtrsim 0.1$ M$_{\odot}$. By investigating the shattering effect on the mm-sized grains, we ensure that the large grains are not significantly fragmented after being injected in the envelope. We conclude that the mm-sized grains observed in the protostellar envelopes are not formed in the envelopes but formed in the inner parts of the star-forming regions and transported to the envelopes before a significant mass growth of the central object, and that they survive in the envelopes.

Sterile neutrinos are thermalised in the early Universe via oscillations with the active neutrinos for certain mixing parameters. The most detailed calculation of this thermalisation process involves the solution of the momentum-dependent quantum kinetic equations, which track the evolution of the neutrino phase space distributions. Until now the collision terms in the quantum kinetic equations have always been approximated using equilibrium distributions, but this approximation has never been checked numerically. In this work we revisit the sterile neutrino thermalisation calculation using the full collision term, and compare the results with various existing approximations in the literature. We find a better agreement than would naively be expected, but also identify some issues with these approximations that have not been appreciated previously. These include an unphysical production of neutrinos via scattering and the importance of redistributing momentum through scattering, as well as details of Pauli blocking. Finally, we devise a new approximation scheme, which improves upon some of the shortcomings of previous schemes.

We explore a chameleon type of interacting dark matter-dark energy scenario in which a scalar field adiabatically traces the minimum of an effective potential sourced by the dark matter density. We discuss extensively the effect of this coupling on cosmological observables, especially the parameter degeneracies expected to arise between the model parameters and other cosmological parameters, and then test the model against observations of the cosmic microwave background (CMB) anisotropies and other cosmological probes. We find that the chameleon parameters $\alpha$ and $\beta$, which determine respectively the slope of the scalar field potential and the dark matter-dark energy coupling strength, can be constrained to $\alpha < 0.17$ and $\beta < 0.19$ using CMB data alone. The latter parameter in particular is constrained only by the late Integrated Sachs-Wolfe effect. Adding measurements of the local Hubble expansion rate $H_0$ tightens the bound on $\alpha$ by a factor of two, although this apparent improvement is arguably an artefact of the tension between the local measurement and the $H_0$ value inferred from Planck data in the minimal $\Lambda$CDM model. The same argument also precludes chameleon models from mimicking a dark radiation component, despite a passing similarity between the two scenarios in that they both delay the epoch of matter-radiation equality. Based on the derived parameter constraints, we discuss possible signatures of the model for ongoing and future large-scale structure surveys.

Starting from the collisional Boltzmann equation, we derive for the first time and from first principles the Boltzmann hierarchy for neutrinos including interactions with a scalar particle. Such interactions appear, for example, in majoron-like models of neutrino mass generation. We study two limits of the scalar mass: (i) An extremely massive scalar whose only role is to mediate an effective 4-fermion neutrino-neutrino interaction, and (ii) a massless scalar that can be produced in abundance and thus demands its own Boltzmann hierarchy. In contrast to, e.g., the first-order Boltzmann hierarchy for Thomson-scattering photons, our interacting neutrino/scalar Boltzmann hierarchies contain additional momentum-dependent collision terms arising from a non-negligible energy transfer in the neutrino-neutrino and neutrino-scalar interactions. This necessitates that we track each momentum mode of the phase space distributions individually, even if the particles were massless. Comparing our hierarchy with the commonly used $(c_{\rm eff}^2,c_{\rm vis}^2)$-parameterisation, we find no formal correspondence between the two approaches, which raises the question of whether the latter parameterisation even has an interpretation in terms of particle scattering. Lastly, although we have invoked majoron-like models as a motivation for our study, our treatment is in fact generally applicable to all scenarios in which the neutrino and/or other ultrarelativistic fermions interact with scalar particles.

We develop a higher-order perturbation theory for large-scale structure formation involving a free-streaming hot or warm dark matter species. We focus on the case of mixed cold dark matter and massive neutrinos, although our approach is applicable also to a single warm dark matter species. In order to capture the suppressed growth of neutrino density perturbations on small scales, we account for the full momentum dependence of the phase space distribution using the Vlaslov equation, and derive from it a formal closed-form nonlinear equation for the neutrino density. Using a systematic perturbative expansion of this equation we compute high-order corrections to the neutrino density contrast without the explicit need to track the perturbed neutrino momentum distribution. We calculate the leading-order total matter bispectrum for several neutrino masses. Using our result as a benchmark, we test the accuracy of the fluid approximation and a linear approximation used in perturbative and N-body analyses, as well as a new hybrid approach that combines the exact linear evolution with the nonlinear structure of the fluid equations. Aiming at $\lesssim1%$ accuracy, we find that the total matter bispectrum with a low neutrino mass m=0.046 eV can be reproduced by all but the fluid approximation, while for larger neutrino masses m=0.46-0.93 eV only the hybrid approach has the desired accuracy on a large range of scales. This result serves as a cautionary note that approximate nonlinear models of neutrino clustering that reproduce the gross features of some observables may not suffice for precision calculations, nor are they guaranteed to apply to other observables. All of the approximation schemes fail to reproduce the bispectrum of the neutrino density perturbations at better than 20% accuracy across all scales, indicating that an exact treatment of nonlinear neutrino perturbations is necessary.

We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch, assumed here to occur at a temperature $T_{\rm QCD} \sim 170$ MeV; this leaves an axion population of such low density that its late-time cosmological impact is negligible. For larger axion masses, $m_a \gtrsim 0.15$ eV, where axions remain in equilibrium until after the QCD phase transition, we find that a Euclid-like survey combined with Planck CMB data can detect $m_a$ at very high significance. Our conclusions are robust against assumptions about prior knowledge of the neutrino mass. Given that the proposed IAXO solar axion search is sensitive to $m_a\lesssim 0.2$ eV, the axion mass range probed by cosmology is nicely complementary.

We assess the potential of a future large-volume photometric redshift survey to constrain observational inflationary parameters using three large-scale structure observables: the angular shear and galaxy power spectra, and the cluster mass function measured through weak lensing. When used in combination with Planck-like CMB measurements, we find that the spectral index n_s can be constrained to a 1 sigma precision of up to 0.0025. The sensitivity to the running of the spectral index can potentially improve to 0.0017, roughly a factor of five better than the present 1 sigma~constraint from Planck and auxiliary CMB data, allowing us to test the assumptions of the slow-roll scenario with unprecedented accuracy. Interestingly, neither CMB+shear nor CMB+galaxy nor CMB+clusters alone can achieve this level of sensitivity; it is the combined power of all three probes that conspires to break the different parameter degeneracies inherent in each type of observations. We make our forecast software publicly available via download or upon request from the authors.

We perform a detailed forecast on how well a {\sc Euclid}-like survey will be able to constrain dark energy and neutrino parameters from a combination of its cosmic shear power spectrum, galaxy power spectrum, and cluster mass function measurements. We find that the combination of these three probes vastly improves the survey's potential to measure the time evolution of dark energy. In terms of a dark energy figure-of-merit defined as $(\sigma(w_{\mathrm p}) \sigma(w_a))^{-1}$, we find a value of 690 for {\sc Euclid}-like data combined with {\sc Planck}-like measurements of the cosmic microwave background (CMB) anisotropies in a 10-dimensional cosmological parameter space, assuming a $\Lambda$CDM fiducial cosmology. For the more commonly used 7-parameter model, we find a figure-of-merit of 1900 for the same data combination. We consider also the survey's potential to measure dark energy perturbations in models wherein the dark energy is parameterised as a fluid with a nonstandard non-adiabatic sound speed, and find that in an \emph{optimistic} scenario in which $w_0$ deviates by as much as is currently observationally allowed from $-1$, models with $\hat{c}_\mathrm{s}^2 = 10^{-6}$ and $\hat{c}_\mathrm{s}^2 = 1$ can be distinguished at more than $2\sigma$ significance. We emphasise that constraints on the dark energy sound speed from cluster measurements are strongly dependent on the modelling of the cluster mass function; significantly weaker sensitivities ensue if we modify our model to include fewer features of nonlinear dark energy clustering. Finally, we find that the sum of neutrino masses can be measured with a $1 \sigma$ precision of 0.015~eV, (abridged)

We present a new analysis on how to distinguish between isotropic and anisotropic cosmological models based on tracking the angular displacements of a large number of distant quasars over an extended period of time, and then performing a multipole-vector decomposition of the resulting displacement maps. We find that while the GAIA mission operating at its nominal specifications does not have sufficient angular resolution to resolve anisotropic universes from isotropic ones using this method within a reasonable timespan of ten years, a next-generation GAIA-like survey with a resolution ten times better should be equal to the task. Distinguishing between different anisotropic models is however more demanding. Keeping the observational timespan to ten years, we find that the angular resolution of the survey will need to be of order 0.1 micro-arcsecond in order for certain rotating anisotropic models to produce a detectable signature that is also unique to models of this class. However, should such a detection become possible, it would immediately allow us to rule out large local void models.

We use cosmological observations in the post-Planck era to derive limits on thermally produced cosmological axions. In the early universe such axions contribute to the radiation density and later to the hot dark matter fraction. We find an upper limit m_a < 0.67 eV at 95% C.L. after marginalising over the unknown neutrino masses, using CMB temperature and polarisation data from Planck and WMAP respectively, the halo matter power spectrum extracted from SDSS-DR7, and the local Hubble expansion rate H_0 released by the Carnegie Hubble Program based on a recalibration of the Hubble Space Telescope Key Project sample. Leaving out the local H_0 measurement relaxes the limit somewhat to 0.86 eV, while Planck+WMAP alone constrain the axion mass to 1.01 eV, the first time an upper limit on m_a has been obtained from CMB data alone. Our axion limit is therefore not very sensitive to the tension between the Planck-inferred H_0 and the locally measured value. This is in contrast with the upper limit on the neutrino mass sum, which we find here to range from 0.27 eV at 95% C.L. combining all of the aforementioned observations, to 0.84 eV from CMB data alone.

We perform a detailed forecast on how well a Euclid-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of sum m_nu ~ 0.06 eV in the normal hierarchy can be detected at 1.5 sigma to 2.5 sigma significance, depending on the model complexity, using a combination of galaxy and cosmic shear power spectrum measurements in conjunction with CMB temperature and polarisation observations from Planck. With better knowledge of the galaxy bias, the significance of the detection could potentially reach 5.4 sigma. Interestingly, neither Planck+shear nor Planck+galaxy alone can achieve this level of sensitivity; it is the combined effect of galaxy and cosmic shear power spectrum measurements that breaks the persistent degeneracies between the neutrino mass, the physical matter density, and the Hubble parameter. Notwithstanding this remarkable sensitivity to sum m_nu, Euclid-like shear and galaxy data will not be sensitive to the exact mass spectrum of the neutrino sector; no significant bias (< 1 sigma) in the parameter estimation is induced by fitting inaccurate models of the neutrino mass splittings to the mock data, nor does the goodness-of-fit of these models suffer any significant degradation relative to the true one (Delta chi_eff ^2< 1).

This is part two in a series of papers in which we investigate an approach based on Lagrangian perturbation theory (LPT) to study the non-linear evolution of the large-scale structure distribution in the universe. Firstly, we compute the matter bispectrum in real space using LPT up one-loop order, for both Gaussian and non-Gaussian initial conditions. In the initial position limit, we find that the one-loop bispectrum computed in this manner is identical to its counterpart obtained from standard Eulerian perturbation theory (SPT). Furthermore, the LPT formalism allows for a simple reorganisation of the perturbative series corresponding to the resummation of an infinite series of perturbations in SPT. Applying this method, we find a resummed one-loop bispectrum that compares favourably with results from N-body simulations. We generalise the resummation method also to the computation of the redshift-space bispectrum up to one loop.

Future cluster surveys will observe galaxy clusters numbering in the hundred thousands. We consider this work how these surveys can be used to constrain dark energy parameters: in particular, the equation of state parameter w and the non-adiabatic sound speed c_s^2. We demonstrate that, in combination with Cosmic Microwave Background (CMB) observations from Planck, cluster surveys such as that in the ESA Euclid project will be able to determine a time-independent w with subpercent precision. Likewise, if the dark energy sound horizon falls within the length scales probed by the cluster survey, then c_s^2 can be pinned down to within an order of magnitude. In the course of this work, we also investigate the process of dark energy virialisation in the presence of an arbitrary sound speed. We find that dark energy clustering and virialisation can lead to dark energy contributing to the total cluster mass at approximately the 0.1% level at maximum.

We consider a generic type of dark energy fluid, characterised by a constant equation of state parameter w and sound speed c_s, and investigate the impact of dark energy clustering on cosmic structure formation using the spherical collapse model. Along the way, we also discuss in detail the evolution of dark energy perturbations in the linear regime. We find that the introduction of a finite sound speed into the picture necessarily induces a scale-dependence in the dark energy clustering, which in turn affects the dynamics of the spherical collapse in a scale-dependent way. As with other, more conventional fluids, we can define a Jeans scale for the dark energy clustering, and hence a Jeans mass M_J for the dark matter which feels the effect of dark energy clustering via gravitational interactions. For bound objects (halos) with masses M >> M_J, the effect of dark energy clustering is maximal. For those with M << M_J, the dark energy component is effectively homogeneous, and its role in the formation of these structures is reduced to its effects on the Hubble expansion rate. To compute quantitatively the virial density and the linearly extrapolated threshold density, we use a quasi-linear approach which is expected to be valid up to around the Jeans mass. We find an interesting dependence of these quantities on the halo mass M, given some w and c_s. The dependence is the strongest for masses lying in the vicinity of M ~ M_J. Observing this M-dependence will be a tell-tale sign that dark energy is dynamic, and a great leap towards pinning down its clustering properties.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We investigate the question of whether the recent modulation signal claimed by CoGeNT is best explained by the dark matter (DM) hypothesis from a Bayesian model comparison perspective. We consider five phenomenological explanations for the data: no modulation signal, modulation due to DM, modulation due to DM compatible with the total CoGeNT rate, and a signal coming from other physics with a free phase but annual period, or with a free phase and a free period. In each scenario, we assign to the free parameters physically motivated priors. We find that while the DM models are weakly preferred to the no modulation model, but when compared to models where the modulation is due to other physics, the DM hypothesis is favoured with odds ranging from 185:1 to 560:1. This result is robust even when astrophysical uncertainties are taken into account and the impact of priors assessed. Interestingly, the odds for the DM model in which the modulation signal is compatible with the total rate against a DM model in which this prior is not implemented is only 5:8, in spite of the former's prediction of a modulation amplitude in the energy range 0.9 to 3.0 keVee that is significantly smaller than the value observed by CoGeNT. Classical hypothesis testing also rules out the null hypothesis of no modulation at the 1.6 sigma to 2.3 sigma level, depending on the details of the alternative. Lastly, we investigate whether anisotropic velocity distributions can help to mitigate the tension between the CoGeNT total and modulated rates, and find encouraging results.

I give an overview of the effects of neutrino masses in cosmology, focussing on the role they play in the evolution of cosmological perturbations. I discuss how recent observations of the cosmic microwave background anisotropies and the large-scale matter distribution can probe neutrino masses with greater precision than current laboratory experiments. I describe several new techniques that will be used to probe cosmology in the future, as well as recent advances in the computation of the nonlinear matter power spectrum and related observables.

We study cosmological models that contain sterile neutrinos with eV-range masses as suggested by reactor and short-baseline oscillation data. We confront these models with both precision cosmological data (probing the CMB decoupling epoch) and light-element abundances (probing the BBN epoch). In the minimal LambdaCDM model, such sterile neutrinos are strongly disfavoured by current data because they contribute too much hot dark matter. However, if the cosmological framework is extended to include also additional relativistic degrees of freedom -- beyond the three standard neutrinos and the putative sterile neutrinos, then the hot dark matter constraint on the sterile states is considerably relaxed. A further improvement is achieved by allowing a dark energy equation of state parameter w<-1. While BBN strongly disfavours extra radiation beyond the assumed eV-mass sterile neutrino, this constraint can be circumvented by a small nu_e degeneracy. Any model containing eV-mass sterile neutrinos implies also strong modifications of other cosmological parameters. Notably, the inferred cold dark matter density can shift up by 20 to 75% relative to the standard LambdaCDM value.

Bayesian statistical methods offer a simple and consistent framework for incorporating uncertainties into a multi-parameter inference problem. In this work we apply these methods to a selection of current direct dark matter searches. We consider the simplest scenario of spin-independent elastic WIMP scattering, and infer the WIMP mass and cross-section from the experimental data with the essential systematic uncertainties folded into the analysis. We find that when uncertainties in the scintillation efficiency of Xenon100 have been accounted for, the resulting exclusion limit is not sufficiently constraining to rule out the CoGeNT preferred parameter region, contrary to previous claims. In the same vein, we also investigate the impact of astrophysical uncertainties on the preferred WIMP parameters. We find that within the class of smooth and isotropic WIMP velocity distributions, it is difficult to reconcile the DAMA and the CoGeNT preferred regions by tweaking the astrophysics parameters alone. If we demand compatibility between these experiments, then the inference process naturally concludes that a high value for the sodium quenching factor for DAMA is preferred.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

Precision cosmology and big-bang nucleosynthesis mildly favor extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to derive credible regions for the assumed common mass scale m_s and effective number N_s of thermally excited sterile neutrino states. Our results are compatible with the existence of one or perhaps two sterile neutrinos, as suggested by LSND and MiniBooNE, if m_s is in the sub-eV range.

The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape information currently contained in the SDSS DR7 halo power spectrum (HPS) is by far superseded by geometric information derived from the baryonic features. An immediate corollary is that contrary to popular beliefs, the upper limit on the neutrino mass m_\nu presently derived from LSS combined with cosmic microwave background (CMB) data does not in fact arise from the possible small-scale power suppression due to neutrino free-streaming, if we limit the model framework to minimal LambdaCDM+m_\nu. However, in more complicated models, such as those extended with extra light degrees of freedom and a dark energy equation of state parameter w differing from -1, shape information becomes crucial for the resolution of parameter degeneracies. This conclusion will remain true even when data from the Planck surveyor become available. In the course of our analysis, we introduce a new dewiggling procedure that allows us to extend consistently the use of the SDSS HPS to models with an arbitrary sound horizon at decoupling. All the cases considered here are compatible with the conservative 95%-bounds \sum m_\nu < 1.16 eV, N_eff = 4.8 \pm 2.0.

We update cosmological hot dark matter constraints on neutrinos and hadronic axions. Our most restrictive limits use 7-year data from the Wilkinson Microwave Anisotropy Probe for the cosmic microwave background anisotropies, the halo power spectrum (HPS) from the 7th data release of the Sloan Digital Sky Survey, and the Hubble constant from Hubble Space Telescope observations. We find 95% C.L. upper limits of \sum m_\nu<0.44 eV (no axions), m_a<0.91 eV (assuming \sum m_\nu=0), and \sum m_\nu<0.41 eV and m_a<0.72 eV for two hot dark matter components after marginalising over the respective other mass. CMB data alone yield \sum m_\nu<1.19 eV (no axions), while for axions the HPS is crucial for deriving m_a constraints. This difference can be traced to the fact that for a given hot dark matter fraction axions are much more massive than neutrinos.

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 discuss current cosmological constraints on axions, as well as future sensitivities. Bounds on axion hot dark matter are discussed first, and subsequently we discuss both current and future sensitivity to models in which axions play the role as cold dark matter, but where the Peccei-Quinn symmetry is not restored during reheating.

We explore the cosmological sensitivity to the amplitude of isocurvature fluctuations that would be caused by axions in the "anthropic window" where the axion decay constant f_a >> 10^12 GeV and the initial misalignment angle Theta_i << 1. In a minimal Lambda-CDM cosmology extended with subdominant scale-invariant isocurvature fluctuations, existing data constrain the isocurvature fraction to alpha < 0.09 at 95% C.L. If no signal shows up, Planck can improve this constraint to 0.042 while an ultimate CMB probe limited only by cosmic variance in both temperature and E-polarisation can reach 0.017, about a factor of five better than the current limit. In the parameter space of f_a and H_I (Hubble parameter during inflation) we identify a small region where axion detection remains within the reach of realistic cosmological probes.

We present the first systematic derivation of the one-loop correction to the large scale matter power spectrum in a mixed cold+hot dark matter cosmology with subdominant massive neutrino hot dark matter. Starting with the equations of motion for the density and velocity fields, we derive perturbative solutions to these quantities and construct recursion relations for the interaction kernels, noting and justifying all approximations along the way. We find interaction kernels similar to those for a cold dark matter-only universe, but with additional dependences on the neutrino energy density fraction f_nu and the linear growth functions of the incoming wavevectors. Compared with the f_nu=0 case, the one-loop corrected matter power spectrum for a mixed dark matter cosmology exhibits a decrease in small scale power exceeding the canonical ~8 f_nu suppression predicted by linear theory, a feature also seen in multi-component N-body simulations.

We revisit the issue of ripples in the primordial power spectra caused by trans-Planckian physics, and the potential for their detection by future cosmological probes. We find that for reasonably large values of the first slow-roll parameter epsilon (> 0.001), a positive detection of trans-Planckian ripples can be made even if the amplitude is as low as 10^-4. Data from the Large Synoptic Survey Telescope (LSST) and the proposed future 21 cm survey with the Fast Fourier Transform Telescope (FFTT) will be particularly useful in this regard. If the scale of inflation is close to its present upper bound, a scale of new physics as high as 0.2 M_Planck could lead to observable signatures.

We explore and compare the performances of two nonlinear correction and scale-dependent biasing models for the extraction of cosmological information from galaxy power spectrum data, especially in the context of beyond-LCDM cosmologies. The first model is the well known Q model, first applied in the analysis of 2dFGRS data. The second, the P model, is inspired by the halo model, in which nonlinear evolution and scale-dependent biasing are encapsulated in a single non-Poisson shot noise term. We find that while both models perform equally well in providing adequate correction for a range of galaxy clustering data in standard LCDM cosmology and in extensions with massive neutrinos, the Q model can give unphysical results in cosmologies containing a subdominant free-streaming dark matter whose temperature depends on the particle mass, e.g., relic thermal axions, unless a suitable prior is imposed on the correction parameter. This last case also exposes the danger of analytic marginalisation, a technique sometimes used in the marginalisation of nuisance parameters. In contrast, the P model suffers no undesirable effects, and is the recommended nonlinear correction model also because of its physical transparency.

We update our previous constraints on two-component hot dark matter (axions and neutrinos), including the recent WMAP 5-year data release. Marginalising over sum m_nu provides m_a < 1.02 eV (95% C.L.) for the axion mass. In the absence of axions we find sum m_nu < 0.63 eV (95% C.L.).

We estimate the effect of the experimental uncertainty in the measurement of the temperature of the cosmic microwave background (CMB) on the extraction of cosmological parameters from future CMB surveys. We find that even for an ideal experiment limited only by cosmic variance up to l = 2500 for both the temperature and polarisation measurements, the projected cosmological parameter errors are remarkably robust against the uncertainty of 1 mK in the FIRAS instrument's CMB temperature monopole measurement. The maximum degradation in sensitivity is 20%, for the baryon density estimate, relative to the case in which the monopole is known infinitely well. While this degradation is acceptable, we note that reducing the uncertainty in the current temperature measurement by a factor of five will bring it down to the per cent level. We also estimate the effect of the uncertainty in the dipole temperature measurement. Assuming the overall calibration of the data to be dominated by the dipole error of 0.2% from FIRAS, the sensitivity degradation is insignificant and does not exceed 10% in any parameter direction.

We consider the inference of the cosmic radiation density, traditionally parameterised as the effective number of neutrino species N_eff, from precision cosmological data. Paying particular attention to systematic effects, notably scale-dependent biasing in the galaxy power spectrum, we find no evidence for a significant deviation of N_eff from the standard value of N_eff^0=3.046 in any combination of cosmological data sets, in contrast to some recent conclusions of other authors. The combination of all available data in the linear regime prefers, in the context of a ``vanilla+N_eff'' cosmological model, 1.1<N_eff<4.8 (95% C.L.) with a best-fit value of 2.6. Adding data at smaller scales, notably the Lyman-alpha forest, we find 2.2<N_eff<5.8 (95% C.L.) with 3.8 as the best fit. Inclusion of the Lyman-alpha data shifts the preferred N_eff upwards because the sigma_8 value derived from the SDSS Lyman-alpha data is inconsistent with that inferred from CMB. In an extended cosmological model that includes a nonzero mass for N_eff neutrino flavours, a running scalar spectral index and a w parameter for the dark energy, we find 0.8<N_eff<6.1 (95% C.L.) with 3.0 as the best fit.

The existence of an unparticle sector, weakly coupled to the standard model, would have a profound impact on supernova (SN) physics. Emission of energy into the unparticle sector from the core of SN1987A would have significantly shortened the observed neutrino burst. The unparticle interaction with nucleons, neutrinos, electrons and muons is constrained to be so weak that it is unlikely to provide any missing-energy signature at colliders. One important exception are models where scale invariance in the hidden sector is broken by the Higgs vacuum expectation value. In this case the SN emission is suppressed by threshold effects.

We calculate the sensitivity of future cosmic microwave background probes and large scale structure measurements from galaxy redshift surveys to the neutrino mass. We find that, for minimal models with few parameters, a measurement of the matter power spectrum over a large range of redshifts has more constraining power than a single measurement at low redshifts. However, this improvement in sensitivity does not extend to larger models. We also quantify how the non-Gaussian nature of the posterior distribution function with respect to the individual cosmological parameter influences such quantities as the sensitivity and the detection threshold. For realistic assumptions about future large scale structure data, the minimum detectable neutrino mass at 95 % C.L. is about 0.05 eV in the context of a minimal 8-parameter cosmological model. In a more general model framework, however, the detection threshold can increase by as much as a factor of three.

We use observations of the cosmological large-scale structure to derive limits on two-component hot dark matter consisting of mass-degenerate neutrinos and hadronic axions, both components having velocity dispersions corresponding to their respective decoupling temperatures. We restrict the data samples to the safely linear regime, in particular excluding the Lyman-alpha forest. Using standard Bayesian inference techniques we derive credible regions in the two-parameter space of m_a and sum(m_nu). Marginalising over sum(m_nu) provides m_a < 1.2 eV (95% C.L.). In the absence of axions the same data and methods give sum(m_nu) < 0.65 eV (95% C.L.). We also derive limits on m_a for a range of axion-pion couplings up to one order of magnitude larger or smaller than the hadronic value.

Neutrino-neutrino interactions can lead to collective flavour conversion effects in supernovae and in the early universe. We demonstrate that the case of "bipolar" oscillations, where a dense gas of neutrinos and antineutrinos in equal numbers completely converts from one flavour to another even if the mixing angle is small, is equivalent to a pendulum in flavour space. Bipolar flavour conversion corresponds to the swinging of the pendulum, which begins in an unstable upright position (the initial flavour), and passes through momentarily the vertically downward position (the other flavour) in the course of its motion. The time scale to complete one cycle of oscillation depends logarithmically on the vacuum mixing angle. Likewise, the presence of an ordinary medium can be shown analytically to contribute to a logarithmic increase in the bipolar conversion period. We further find that a more complex (and realistic) system of unequal numbers of neutrinos and antineutrinos is analogous to a spinning top subject to a torque. This analogy easily explains how such a system can oscillate in both the bipolar and the synchronised mode, depending on the neutrino density and the size of the neutrino-antineutrino asymmetry. Our simple model applies strictly only to isotropic neutrino gasses. In more general cases, and especially for neutrinos streaming from a supernova core, different modes couple to each other with unequal strength, an effect that can lead to kinematical decoherence in flavour space rather than collective oscillations. The exact circumstances under which collective oscillations occur in non-isotropic media remain to be understood.

High-precision data from observation of the cosmic microwave background and the large scale structure of the universe provide very tight constraints on the effective parameters that describe cosmological inflation. Indeed, within a constrained class of LambdaCDM models, the simple lambda phi^4 chaotic inflation model already appears to be ruled out by cosmological data. In this paper, we compute constraints on inflationary parameters within a more general framework that includes other physically motivated parameters such as a nonzero neutrino mass. We find that a strong degeneracy between the tensor-to-scalar ratio r and the neutrino mass prevents lambda phi^4 from being excluded by present data. Reversing the argument, if lambda phi^4 is the correct model of inflation, it predicts a sum of neutrino masses at 0.3-0.5 eV, a range compatible with present experimental limits and within the reach of the next generation of neutrino mass measurements. We also discuss the associated constraints on the dark matter density, the dark energy equation of state, and spatial curvature, and show that the allowed regions are significantly altered. Importantly, we find an allowed range of 0.094 < Omega_c h^2 < 0.136 for the dark matter density, a factor of two larger than that reported in previous studies. This expanded parameter space may have implications for constraints on SUSY dark matter models.

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 study post weak decoupling coherent active-sterile and active-active matter-enhanced neutrino flavor transformation in the early universe. We show that flavor conversion efficiency at Mikheyev-Smirnov-Wolfenstein resonances is likely to be high (adiabatic evolution) for relevant neutrino parameters and energies. However, we point out that these resonances cannot sweep smoothly and continuously with the expansion of the universe. We show how neutrino flavor conversion in this way can leave both the active and sterile neutrinos with non-thermal energy spectra, and how, in turn, these distorted energy spectra can affect the neutron-to-proton ratio, primordial nucleosynthesis, and cosmological mass/closure constraints on sterile neutrinos. We demonstrate that the existence of a light sterile neutrino which mixes with active neutrinos can change fundamentally the relationship between the cosmological lepton numbers and the primordial nucleosynthesis He-4 yield.

Resonant annihilation of ultra-high energy cosmic neutrinos (UHECnu) on the cosmic neutrino background (CnuB) into Z bosons--the Z-burst mechanism--and its associated absorption and emission phenomenology provide a unique, albeit indirect, probe of the CnuB in its present state. In this paper, we examine the implications of gravitational clustering of the CnuB in nearby galaxy clusters for the Z-burst phenomenology. In particular, we study the emission features of the Z-decay products originating from the Virgo cluster, and the potential of future cosmic ray experiments to observe clustering-enhanced Z-burst rates. We find that a detector with an exposure equivalent to three years of observations at the Extreme Universe Space Observatory (EUSO) will very likely measure these enhanced rates together with the associated UHECnu flux, provided that the latter saturates current observational limits and the neutrino masses are quasi-degenerate, m_nu_i < 0.1 eV. In the case of UHECnu fluxes below the electromagnetic cascade limit, or a hierarchical neutrino mass spectrum, an experimental sensitivity exceeding that of EUSO by at least two orders of magnitude is required to detect the clustering enhancements with any certainty.

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.

We study the gravitational clustering of big bang relic neutrinos onto existing cold dark matter (CDM) and baryonic structures within the flat $\Lambda$CDM model, using both numerical simulations and a semi-analytical linear technique, with the aim of understanding the neutrinos' clustering properties for direct detection purposes. In a comparative analysis, we find that the linear technique systematically underestimates the amount of clustering for a wide range of CDM halo and neutrino masses. This invalidates earlier claims of the technique's applicability. We then compute the exact phase space distribution of relic neutrinos in our neighbourhood at Earth, and estimate the large scale neutrino density contrasts within the local Greisen--Zatsepin--Kuzmin zone. With these findings, we discuss the implications of gravitational neutrino clustering for scattering-based detection methods, ranging from flux detection via Cavendish-type torsion balances, to target detection using accelerator beams and cosmic rays. For emission spectroscopy via resonant annihilation of extremely energetic cosmic neutrinos on the relic neutrino background, we give new estimates for the expected enhancement in the event rates in the direction of the Virgo cluster.

A recent revival of interest in synchronised oscillations due to neutrino--neutrino forward scattering in dense gases has led to two interesting applications with notable outcomes: (i) cosmological bounds on neutrino--antineutrino asymmetries are improved owing to flavour equilibration prior to the onset of big bang nucleosynthesis, and (ii) a neutron-rich environment required for $r$-process nucleosynthesis is shown to be always maintained in a supernova hot bubble irrespective of flavour oscillations, contrary to results from previous studies. I present in this talk a pedagogical review of these works.

The exact Quantum Kinetic Equations for a two-flavour active-sterile neutrino system are used to provide a systematic derivation of approximate evolution equations for the relic neutrino asymmetry. An extension of the adiabatic approximation for matter-affected neutrino oscillations is developed which incorporates decoherence due to collisions. Exact and approximate expressions for the decoherence and repopulation functions are discussed. A first pass is made over the exact treatment of multi-flavour partially incoherent oscillations.

It has been proposed that the annihilation process $\nu\bar{\nu} \to e^{-}e^{+} \to \gamma\gamma$ may be responsible for the generation of gamma-ray bursts (GRBs). The incipient neutrino--antineutrino pairs carry virtually all of the gravitational binding energy available from the central engine. However, gamma-ray bursters proposed to date are inevitably surrounded by an excess of baryons, leading to the "baryon-loading problem". In the light of growing evidence for neutrino oscillations, we discuss the implications of matter-affected oscillations for GRB energetics, and on the viability of "mirror" stars as GRB progenitors.