We show that the excess in electron recoil events seen by the XENON1T experiment can be explained by relatively low-mass Luminous Dark Matter candidate. The dark matter scatters inelastically in the detector (or the surrounding rock), to produce a heavier dark state with a ~2.75 keV mass splitting. This heavier state then decays within the detector, producing a peak in the electron recoil spectrum which is a good fit to the observed excess. We comment on the ability of future direct detection datasets to differentiate this model from other Beyond the Standard Model scenarios, and from possible tritium backgrounds, including the use of diurnal modulation, multi-channel signals etc.,~as possible distinguishing features of this scenario.

Indirect detection of dark matter via its annihilation products is a key technique in the search for dark matter in the form of weakly interacting massive particles (WIMPs). Strong constraints exist on the annihilation of WIMPs to highly visible Standard Model final states such as photons or charged particles. In the case of s-wave annihilation, this typically eliminates thermal relic cross sections for dark matter of mass below $\cal{O}$(10) GeV. However, such limits typically neglect the possibility that dark matter may annihilate to assumed invisible or hard-to-detect final states, such as neutrinos. This is a difficult paradigm to probe due to the weak neutrino interaction cross section. Considering dark matter annihilation in the Galactic halo, we study the prospects for indirect detection using the Hyper-Kamiokande (HyperK) neutrino experiment, for dark matter of mass below 1 GeV. We undertake a dedicated simulation of the HyperK detector, which we benchmark against results from the similar Super-Kamiokande experiment and HyperK physics projections. We provide projections for the annihilation cross-sections that can be probed by HyperK for annihilation to muon or neutrino final states, and discuss uncertainties associated with the dark matter halo profile. For neutrino final states, we find that HyperK is sensitive to thermal annihilation cross-sections for dark matter with mass around 20 MeV, assuming an NFW halo profile. We also discuss the effects of Gadolinium enhancement, and prospects for improving the reach at low mass.

Neutron stars provide a cosmic laboratory to study the nature of dark matter particles and their interactions. Dark matter can be captured by neutron stars via scattering interactions, in which kinetic energy is transferred to the star. This can have a number of observational consequences, such as the heating of old neutron stars to infra-red temperatures. Previous treatments of the capture process have employed various approximation or simplifications. We present here an improved treatment of dark matter capture, valid for a wide dark matter mass range, that correctly incorporates all relevant physical effects. These include gravitational focusing, a fully relativistic scattering treatment, Pauli blocking, neutron star (NS) opacity and multi-scattering effects. We provide general expressions that enable the exact capture rate to be calculated numerically, and derive simplified expressions that are valid for particular interaction types or mass regimes and that greatly increase the computational efficiency. Our formalism is applicable to the scattering of dark matter from any NS constituents, or to the capture of dark matter in other compact objects.

Dark matter particles will be captured in neutron stars if they undergo scattering interactions with nucleons or leptons. These collisions transfer the dark matter kinetic energy to the star, resulting in appreciable heating that is potentially observable by forthcoming infrared telescopes. While previous work considered scattering only on nucleons, neutron stars contain small abundances of other particle species, including electrons and muons. We perform a detailed analysis of the neutron star kinetic heating constraints on leptophilic dark matter. We also estimate the size of loop induced couplings to quarks, arising from the exchange of photons and Z bosons. Despite having relatively small lepton abundances, we find that an observation of an old, cold, neutron star would provide very strong limits on dark matter interactions with leptons, with the greatest reach arising from scattering off muons. The projected sensitivity is orders of magnitude more powerful than current dark matter-electron scattering bounds from terrestrial direct detection experiments.

Many dark matter interaction types lead to annihilation processes which suffer from p-wave suppression or helicity suppression, rendering them subdominant to unsuppressed s-wave processes. We demonstrate that the natural inclusion of dark initial state radiation can open an unsuppressed s-wave annihilation channel, and thus provide the dominant dark matter annihilation process for particular interaction types. We illustrate this effect with the bremsstrahlung of a dark spin-0 or dark spin-1 particle from fermionic dark matter, $\overline{\chi}\chi\rightarrow \overline{f}f\phi$ or $\overline{f}fZ'$. The dark initial state radiation process, despite having a 3-body final state, proceeds at the same order in the new physics scale $\Lambda$ as the annihilation to the 2-body final state $\overline{\chi}\chi\rightarrow \overline{f}f$. This is lower order in $\Lambda$ than the well-studied lifting of helicity suppression via Standard Model final state radiation, or virtual internal bremsstrahlung. This dark bremsstrahlung process should influence LHC and indirect detection searches for dark matter.

In the simplified dark matter models commonly studied, the mass generation mechanism for the dark fields is not typically specified. We demonstrate that the dark matter interaction types, and hence the annihilation processes relevant for relic density and indirect detection, are strongly dictated by the mass generation mechanism chosen for the dark sector particles, and the requirement of gauge invariance. We focus on the class of models in which fermionic dark matter couples to a spin-1 vector or axial-vector mediator. However, in order to generate dark sector mass terms, it is necessary in most cases to introduce a dark Higgs field and thus a spin-0 scalar mediator will also be present. In the case that all the dark sector fields gain masses via coupling to a single dark sector Higgs field, it is mandatory that the axial-vector coupling of the spin-1 mediator to the dark matter is non-zero; the vector coupling may also be present depending on the charge assignments. For all other mass generation options, only pure vector couplings between the spin-1 mediator and the dark matter are allowed. If these coupling restrictions are not obeyed, unphysical results may be obtained such as a violation of unitarity at high energies. These two-mediator scenarios lead to important phenomenology that does not arise in single mediator models. We survey two-mediator dark matter models which contain both vector and scalar mediators, and explore their relic density and indirect detection phenomenology.

We consider the indirect detection signals for a self-consistent hidden $U(1)$ model containing a Majorana dark matter candidate, $\chi$, a dark gauge boson, $Z'$, and a dark Higgs, $s$. Compared with a model containing only a dark matter candidate and $Z'$ mediator, the addition of the scalar provides a mass generation mechanism for the dark sector particles and is required in order to avoid unitarity violation at high energies. We find that the inclusion of the two mediators opens up a new two-body s-wave annihilation channel, $\chi\overline\chi\rightarrow sZ'$. This new process, which is missed in the usual single-mediator simplified model approach, can be the dominant annihilation channel. This provides rich phenomenology for indirect detection searches, allows indirect searches to explore regions of parameter space not accessible with other commonly considered s-wave annihilation processes, and enables both the $Z'$ and scalar couplings to be probed. We examine the phenomenology of the sector with a focus on this new process, and determine the limits on the model parameter space from Fermi data on dwarf spheriodal galaxies and other relevant experiments.

The relic abundance of particle and antiparticle dark matter (DM) need not be vastly different in thermal asymmetric dark matter (ADM) models. By considering the effect of a primordial asymmetry on the thermal Boltzmann evolution of coupled DM and anti-DM, we derive the requisite annihilation cross section. This is used in conjunction with CMB and Fermi-LAT gamma-ray data to impose a limit on the number density of anti-DM particles surviving thermal freeze-out. When the extended gamma-ray emission from the Galactic Center is reanalyzed in a thermal ADM framework, we find that annihilation into $\tau$ leptons prefer anti-DM number densities 1-4$\%$ that of DM while the $b$-quark channel prefers 50-100$\%$.

We consider a scenario where dark matter (DM) interacts exclusively with Standard Model (SM) leptons at tree level. Due to the absence of tree-level couplings to quarks, the constraints on leptophilic dark matter arising from direct detection and hadron collider experiments are weaker than those for a generic WIMP. We study a simple model in which interactions of DM with SM leptons are mediated by a leptophilic $Z'$ boson, and determine constraints on this scenario arising from relic density, direct detection, and other experiments. We then determine current LHC limits and project the future discovery reach. We show that, despite the absence of direct interactions with quarks, this scenario can be strongly constrained.

We study dark matter (DM) models in which there are two dark sector particles, $\chi_1$ and $\chi_2$, of near mass. In such models, co-annihilation of $\chi_1$ and $\chi_2$ may be the dominant process controlling the DM relic density during freezeout in the early universe. In this scenario, there is no significant contribution to direct and indirect detection signals, unless there exists an extreme degeneracy in the masses of the lightest dark sector particles. Therefore, relic density constraints and collider searches provide the most relevant information about these models. We consider Dirac fermion dark matter which couples to standard model (SM) particles via an effective operator. For the collider phenomenology, where an effective field theory may not be valid, we adopt a simple Z' model to provide an appropriate UV completion. We explore the interesting LHC signals that arise from the dark matter production process $pp \rightarrow \overline{\chi_1} + \chi_2 + \textrm{ jet}$, followed by the decay $\chi_2 \rightarrow \chi_1 + SM$.

It has been argued that the existence of old neutron stars excludes the possibility of non-annihilating light bosonic dark matter, such as that arising in asymmetric dark matter scenarios. If non-annihilating dark matter is captured by neutron stars, the density will eventually become sufficient for black hole formation. However, the dynamics of collapse is highly sensitive to dark-matter self-interactions. Repulsive self-interactions, even if extremely weak, can prevent black hole formation. We argue that self-interactions will necessarily be present, and estimate their strength in representative models. We also consider co-annihilation of dark matter with nucleons, which arises naturally in many asymmetric dark matter models, and which again acts to prevent black hole formation. We demonstrate how the excluded region of the dark-matter parameter space shrinks as the strength of such interactions is increased, and conclude that neutron star observations do not exclude most realistic bosonic asymmetric dark matter models.

We revisit the calculation of electroweak bremsstrahlung contributions to dark matter annihilation. Dark matter annihilation to leptons is necessarily accompanied by electroweak radiative corrections, in which a $W$ or $Z$ boson is also radiated. Significantly, while many dark matter models feature a helicity suppressed annihilation rate to fermions, bremsstrahlung process can remove this helicity suppression such that the branching ratios Br($\ell \nu W $), Br($\ell^+\ell^-Z$), and Br($\bar\nu \nu Z$) dominate over Br($\ell^+\ell^-$) and Br($\bar\nu \nu$). We find this is most significant in the limit where the dark matter mass is nearly degenerate with the mass of the boson which mediates the annihilation process. Electroweak bremsstrahlung has important phenomenological consequences both for the magnitude of the total dark matter annihilation cross section and for the character of the astrophysical signals for indirect detection. Given that the $W$ and $Z$ gauge bosons decay dominantly via hadronic channels, it is impossible to produce final state leptons without accompanying protons, antiprotons, and gamma rays.

We examine observational signatures of dark matter annihilation in the Milky Way arising from electroweak bremsstrahlung contributions to the annihilation cross section. It has been known for some time that photon bremsstrahlung may significantly boost DM annihilation yields. Recently, we have shown that electroweak bremsstrahlung of W and Z gauge bosons can be the dominant annihilation channel in some popular models with helicity-suppressed 2 --> 2 annihilation. W/Z-bremsstrahlung is particularly interesting because the gauge bosons produced via annihilation subsequently decay to produce large correlated fluxes of electrons, positrons, neutrinos, hadrons (including antiprotons) and gamma rays, which are all of importance in indirect dark matter searches. Here we calculate the spectra of stable annihilation products produced via gamma/W/Z-bremsstrahlung. After modifying the fluxes to account for the propagation through the Galaxy, we set upper bounds on the annihilation cross section via a comparison with observational data. We show that stringent cosmic ray antiproton limits preclude a sizable dark matter contribution to observed cosmic ray positron fluxes in the class of models for which the bremsstrahlung processes dominate.

The similarity of the visible and dark matter abundances indicates that they may originate via the same mechanism. If both the dark and the visible matter are charged under a generalized baryon number which remains always conserved, then the asymmetry of the visible sector may be compensated by an asymmetry in the dark sector. We show how the separation of baryonic and antibaryonic charge can originate in the vacuum, via the Affleck-Dine mechanism, due to the breaking of a symmetry orthogonal to the baryon number. Symmetry restoration in the current epoch guarantees the individual stability of the two sectors.

Dark matter annihilation to leptons, $\chi \chi \rightarrow l \bar{l}$, is necessarily accompanied by electroweak radiative corrections, in which a W or Z boson is radiated from a final state particle. Given that the W and Z gauge bosons decay dominantly via hadronic channels, it is thus impossible to produce final state leptons without accompanying protons, antiprotons, and gamma rays. Significantly, while many dark matter models feature a helicity suppressed annihilation rate to fermions, radiating a massive gauge boson from a final state fermion removes this helicity suppression, such that the branching ratios Br(l \nu W), Br(l^+l^- Z), and Br(\nu\nubar Z) dominate over Br(l\bar{l}). W/Z-bremsstrahlung thus allows indirect detection of many WIMP models that would otherwise be helicity-suppressed, or v^2 suppressed. Antiprotons and even anti-deuterons become consequential final state particles. This is an important result for future DM searches. We discuss the implications of W/Z-bremsstrahlung for "leptonic" DM models which aim to fit recent cosmic ray positron and antiproton data.

We calculate the neutrino signal resulting from annihilation of secluded dark matter in the Sun. In this class of models, dark matter annihilates first into metastable mediators, which subsequently decay into Standard Model particles. If the mediators are long lived, they will propagate out from the dense solar core before decaying. High energy neutrinos undergo absorption in the Sun. In the standard scenario in which neutrinos are produced directly in the centre of the Sun, absorption is relevant for E > 100 GeV, resulting in a significant suppression of the neutrino spectrum beyond E ~ 1 TeV. In the secluded dark matter scenario, the neutrino signal is greatly enhanced because neutrinos are injected away from the core, at lower density. Since the solar density falls exponentially with radius, metastable mediators have a significant effect on the neutrino flux, even for decay lengths which are small compared to the solar radius. Moreover, since neutrino detection cross sections grow with energy, this enhancement of the high energy region of the neutrino spectrum would have a large effect on overall event rates.

The standard cold dark matter cosmological model, while successful in explaining the observed large scale structure of the Universe, tends to overpredict structure on small scales. It has been proposed this problem may be alleviated in a class of late-decaying dark matter models, in which the parent dark matter particle decays to an almost degenerate daughter, plus a relativistic final state. We construct explicit particle physics models that realize this goal while obeying observational constraints. To achieve this, we introduce a pair of fermionic dark matter candidates and a new scalar field, which obey either a Z4 or a U(1) symmetry. Through the spontaneous breaking of these symmetries, and coupling of the new fields to standard model particles, we demonstrate that the desired decay process may be obtained. We also discuss the dark matter production processes in these models.

We consider a class of late-decaying dark-matter models, in which a dark matter particle decays to a heavy stable daughter of approximately the same mass, together with one or more relativistic particles which carry away only a small fraction of the parent rest mass. Such decays can affect galactic halo structure and evolution, and have been invoked as a remedy to some of the small scale structure-formation problems of cold dark matter. There are existing stringent limits on the dark matter lifetime if the decays produce photons. By considering examples in which the relativistic decay products instead consist of neutrinos or electron-position pairs, we derive stringent limits on these scenarios for a wide range of dark matter masses. We thus eliminate a sizable portion of the parameter space for these late decay models if the dominant decay channel involves Standard Model final states.

We determine upper limits on the dark matter (DM) self-annihilation cross section for scenarios in which annihilation leads to the production of electron--positron pairs. In the Galactic centre (GC), relativistic electrons and positrons produce a radio flux via synchroton emission, and a gamma ray flux via bremsstrahlung and inverse Compton scattering. On the basis of archival, interferometric and single-dish radio data, we have determined the radio spectrum of an elliptical region around the Galactic centre of extent 3 degrees semi-major axis (along the Galactic plane) and 1 degree semi-minor axis and a second, rectangular region, also centered on the GC, of extent 1.6 degrees x 0.6 degrees. The radio spectra of both regions are non-thermal over the range of frequencies for which we have data: 74 MHz -- 10 GHz. We also consider gamma-ray data covering the same region from the EGRET instrument (about GeV) and from HESS (around TeV). We show how the combination of these data can be used to place robust constraints on DM annihilation scenarios, in a way which is relatively insensitive to assumptions about the magnetic field amplitude in this region. Our results are approximately an order of magnitude more constraining than existing Galactic centre radio and gamma ray limits. For a DM mass of m_\chi =10 GeV, and an NFW profile, we find that the velocity-averaged cross-section must be less than a few times 10^-25 cm^3 s^-1.

Dark matter annihilation into charged particles is necessarily accompanied by gamma rays, produced via radiative corrections. Internal bremsstrahlung from the final state particles can produce hard gamma rays up to the dark matter mass, with an approximately model-independent spectrum. Focusing on annihilation into electrons, we compute robust upper bounds on the dark matter self annihilation cross section $<\sigma_A v >_{e^+e^-}$ using gamma ray data from the Milky Way spanning a wide range of energies, $\sim10^{-3} - 10^4$ GeV. We also compute corresponding bounds for the other charged leptons. We make conservative assumptions about the astrophysical inputs, and demonstrate how our derived bounds would be strengthened if stronger assumptions about these inputs are adopted. The fraction of hard gamma rays near the endpoint accompanying annihilation to $e^+e^-$ is only a factor of $\alt 10^2$ lower than for annihilation directly to monoenergetic gamma rays. The bound on $<\sigma_A v >_{e^+e^-}$ is thus weaker than that for $<\sigma_A v >_{\gamma\gamma}$ by this same factor. The upper bounds on the annihilation cross sections to charged leptons are compared with an upper bound on the {\it total} annihilation cross section defined by neutrinos.

A conservative upper bound on the total dark matter (DM) annihilation rate can be obtained by constraining the appearance rate of the annihilation products which are hardest to detect. The production of neutrinos, via the process $\chi \chi \to \bar\nu \nu $, has thus been used to set a strong general bound on the dark matter annihilation rate. However, Standard Model radiative corrections to this process will inevitably produce photons which may be easier to detect. We present an explicit calculation of the branching ratios for the electroweak bremsstrahlung processes $\chi \chi \to \bar\nu \nu Z$ and $\chi \chi \to \bar\nu e W$. These modes inevitably lead to electromagnetic showers and further constraints on the DM annihilation cross-section. In addition to annihilation, our calculations are also applicable to the case of dark matter decay.

Using gamma-ray data from observations of the Milky Way, Andromeda (M31), and the cosmic background, we calculate conservative upper limits on the dark matter self-annihilation cross section to monoenergetic gamma rays, <sigma_A v>_{gamma gamma}, over a wide range of dark matter masses. (In fact, over most of this range, our results are unchanged if one considers just the branching ratio to gamma rays with energies within a factor of a few of the endpoint at the dark matter mass.) If the final-state branching ratio to gamma rays, Br(gamma gamma), were known, then <sigma_A v>_{gamma gamma} / Br(gamma gamma) would define an upper limit on the total cross section; we conservatively assume Br(gamma gamma) > 10^{-4}. An upper limit on the total cross section can also be derived by considering the appearance rates of any Standard Model particles; in practice, this limit is defined by neutrinos, which are the least detectable. For intermediate dark matter masses, gamma-ray-based and neutrino-based upper limits on the total cross section are comparable, while the gamma-ray limit is stronger for small masses and the neutrino limit is stronger for large masses. We comment on how these results depend on the assumptions about astrophysical inputs and annihilation final states, and how GLAST and other gamma-ray experiments can improve upon them.

We consider dark matter annihilation into Standard Model particles and show that the least detectable final states, namely neutrinos, define an upper bound on the total cross section. Calculating the cosmic diffuse neutrino signal, and comparing it to the measured terrestrial atmospheric neutrino background, we derive a strong and general bound. This can be evaded if the annihilation products are dominantly new and truly invisible particles. Our bound is much stronger than the unitarity bound at the most interesting masses, shows that dark matter halos cannot be significantly modified by annihilations, and can be improved by a factor of 10--100 with existing neutrino experiments.

We investigate signatures of neutrino scattering in the Cosmic Microwave Background (CMB) and matter power spectra, and the extent to which present cosmological data can distinguish between a free streaming or tightly coupled fluid of neutrinos. If neutrinos have strong non-standard interactions, for example, through the coupling of neutrinos to a light boson, they may be kept in equilibrium until late times. We show how the power spectra for these models differ from more conventional neutrino scenarios, and use CMB and large scale structure data to constrain these models. CMB polarization data improves the constraints on the number of massless neutrinos, while the Lyman--$\alpha$ power spectrum improves the limits on the neutrino mass. Neutrino mass limits depend strongly on whether some or all of the neutrino species interact and annihilate. The present data can accommodate a number of tightly-coupled relativistic degrees of freedom, and none of the interacting-neutrino scenarios considered are ruled out by current data -- although considerations regarding the age of the Universe disfavor a model with three annihilating neutrinos with very large neutrino masses.

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.

We discuss the prospects for next generation neutrino telescopes, such as IceCube, to measure the flavor ratios of high-energy astrophysical neutrinos. The expected flavor ratios at the sources are $\phi_{\nu_e}:\phi_{\nu_{\mu}}:\phi_{\nu_{\tau}} = 1:2:0$, and neutrino oscillations quickly transform these to $1:1:1$. The flavor ratios can be deduced from the relative rates of showers ($\nu_e$ charged-current, most $\nu_\tau$ charged-current, and all flavors neutral-current), muon tracks ($\nu_\mu$ charged-current only), and tau lepton lollipops and double-bangs ($\nu_\tau$ charged-current only). The peak sensitivities for these interactions are at different neutrino energies, but the flavor ratios can be reliably connected by a reasonable measurement of the spectrum shape. Measurement of the astrophysical neutrino flavor ratios tests the assumed production mechanism and also provides a very long baseline test of a number of exotic scenarios, including neutrino decay, CPT violation, and small-$\delta m^2$ oscillations to sterile neutrinos.

The Galactic positrons, as observed by their annihilation gamma-ray line at 0.511 MeV, are difficult to account for with astrophysical sources. It has been proposed that they are produced instead by dark matter annihilation or decay in the inner Galactic halo. To avoid other constraints, these processes are required to occur "invisibly," such that the eventual positron annihilation is the only detectable signal. However, electromagnetic radiative corrections to these processes inevitably produce real gamma rays (``internal bremsstrahlung"); this emission violates COMPTEL and EGRET constraints unless the dark matter mass is less than about 20 MeV.

We consider the consequences for the relic neutrino abundance if extra neutrino interactions are allowed, e.g., the coupling of neutrinos to a light (compared to $m_\nu$) boson. For a wide range of couplings not excluded by other considerations, the relic neutrinos would annihilate to bosons at late times, and thus make a negligible contribution to the matter density today. This mechanism evades the neutrino mass limits arising from large scale structure.

Neutrinos may be pseudo-Dirac states, such that each generation is actually composed of two maximally-mixed Majorana neutrinos separated by a tiny mass difference. The usual active neutrino oscillation phenomenology would be unaltered if the pseudo-Dirac splittings are $\delta m^2 \alt 10^{-12}$ eV$^2$; in addition, neutrinoless double beta decay would be highly suppressed. However, it may be possible to distinguish pseudo-Dirac from Dirac neutrinos using high-energy astrophysical neutrinos. By measuring flavor ratios as a function of $L/E$, mass-squared differences down to $\delta m^2 \sim 10^{-18}$ eV$^2$ can be reached. We comment on the possibility of probing cosmological parameters with neutrinos.

We have previously shown that the decay of high-energy neutrinos from distant astrophysical sources would be revealed by flavor ratios that deviate strongly from the $\phi_{\nu_e}:\phi_{\nu_\mu}:\phi_{\nu_\tau} = 1:1:1$ expected from oscillations alone. Here we show that the deviations are significantly larger when the mixing angle $\theta_{13}$ and the CP phase $\delta$ are allowed to be nonzero. If neutrinos decay, this could allow measurement of $\theta_{13}$ and $\delta$ in IceCube and other near-term neutrino telescopes.

We study a system containing many particles of identical kinematics with a zero range interaction that scatters one from the other, and with the possible exchange of an attribute. Taking an initial condition in which the attribute is asymmetrically distributed in the regions of momentum space occupied by the particles, we study the rate at which it becomes uniformly distributed, through collisions. We find, in some circumstances, a rate that is much faster than that which would be estimated from cross-sections. This behavior is attributable in some general sense to N-particle entanglement. We suggest applications to neutrino physics, where the attribute is neutrino flavor.

Existing limits on the non-radiative decay of one neutrino to another plus a massless particle (e.g., a singlet Majoron) are very weak. The best limits on the lifetime to mass ratio come from solar neutrino observations, and are $\tau/m \agt 10^{-4}$ s/eV for the relevant mass eigenstate(s). For lifetimes even several orders of magnitude longer, high-energy neutrinos from distant astrophysical sources would decay. This would strongly alter the flavor ratios from the $\phi_{\nu_e}:\phi_{\nu_{\mu}}:\phi_{\nu_{\tau}} = 1:1:1$ expected from oscillations alone, and should be readily visible in the near future in detectors such as IceCube.

We assess a mechanism which can transform neutrino-antineutrino asymmetries between flavors in the early universe, and confirm that such transformation is unavoidable in the near bi-maximal framework emerging for the neutrino mixing matrix. We show that the process is a standard Mikheyev-Smirnov-Wolfenstein flavor transformation dictated by a synchronization of momentum states. We also show that flavor ``equilibration'' is a special feature of maximal mixing, and carefully examine new constraints placed on neutrino asymmetries. In particular, the big bang nucleosynthesis limit on electron neutrino degeneracy xi_e < 0.04 does not apply directly to all flavors, yet confirmation of the large-mixing-angle solution to the solar neutrino problem will eliminate the possibility of degenerate big bang nucleosynthesis.

Despite the fact that the solar neutrino flux is now well-understood in the context of matter-affected neutrino mixing, we find that it is not yet possible to set a strong and model-independent bound on solar neutrino decays. If neutrinos decay into truly invisible particles, the Earth-Sun baseline defines a lifetime limit of $\tau/m \agt 10^{-4}$ s/eV. However, there are many possibilities which must be excluded before such a bound can be established. There is an obvious degeneracy between the neutrino lifetime and the mixing parameters. More generally, one must also allow the possibility of active daughter neutrinos and/or antineutrinos, which may partially conceal the characteristic features of decay. Many of the most exotic possibilities that presently complicate the extraction of a decay bound will be removed if the KamLAND reactor antineutrino experiment confirms the large-mixing angle solution to the solar neutrino problem and measures the mixing parameters precisely. Better experimental and theoretical constraints on the $^8$B neutrino flux will also play a key role, as will tighter bounds on absolute neutrino masses. Though the lifetime limit set by the solar flux is weak, it is still the strongest direct limit on non-radiative neutrino decay. Even so, there is no guarantee (by about eight orders of magnitude) that neutrinos from astrophysical sources such as a Galactic supernova or distant Active Galactic Nuclei will not decay.

We investigate the mixing of heavy gauge singlet neutrinos in a mirror matter model employing the seesaw mechanism. The parameter constraints that must be satisfied to prevent the overproduction of mirror matter in the early universe are deduced. We find that no fine tuning in the heavy neutral fermion sector is required for this mirror matter model to satisfy cosmological constraints. Baryogenesis scenarios are briefly discussed in the context of the mirror model.

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.

A consequence of the evaporation of primordial black holes in the early universe may be the generation of mirror matter. This would have implications with regard to dark matter, and the number of light particle species in equilibrium at the time of big bang nucleosynthesis. The possibilities for the production of mirror matter by this mechanism are explored.

Oscillations between ordinary and sterile neutrinos can generate large neutrino asymmetries in the early universe. These asymmetries can significantly affect big bang nucleosynthesis (BBN) through modification of nuclear reaction rates. We study this phenomenon within a model consisting of the three ordinary neutrinos plus one sterile neutrino that can be motivated by the neutrino anomalies and the dark matter problem. We calculate how the lepton asymmetries produced evolve at temperatures where they impact on BBN. The effect of the asymmetries on primordial helium production is determined, leading to an effective number of neutrino flavours during BBN of either about 2.7 or 3.1 depending on the sign of the lepton asymmetry.