An old neutron star (NS) may capture halo dark matter (DM) and get heated up by the deposited kinetic energy, thus behaving like a thermal DM detector with sensitivity to a wide range of DM masses and a variety of DM-quark interactions. Near future infrared telescopes will measure NS temperatures down to a few thousand Kelvin and probe NS heating by DM capture. We focus on GeV-mass Dirac fermion DM (which is beyond the reach of current DM direct detection experiments) in scenarios in which the DM capture rate can saturate the geometric limit. For concreteness, we study (1) a model that invokes dark decays of the neutron to explain the neutron lifetime anomaly, and (2) a framework of DM coupled to quarks through a vector current portal. In the neutron dark decay model a NS can have a substantial DM population, so that the DM capture rate can reach the geometric limit through DM self-interactions even if the DM-neutron scattering cross section is tiny. We find NS heating to have greater sensitivity than multi-pion signatures in large underground detectors for the neutron dark decay model, and sub-GeV gamma-ray signatures for the quark vector portal model.

We study the stochastic background of gravitational waves which accompany the sudden freeze-out of dark matter triggered by a cosmological first order phase transition that endows dark matter with mass. We consider models that produce the measured dark matter relic abundance via (1) bubble filtering, and (2) inflation and reheating, and show that gravitational waves from these mechanisms are detectable at future interferometers.

We introduce a model in which the genesis of dark matter and neutrino masses is associated with a first order phase transition of a scalar singlet field. During the phase transition a source right-handed neutrino acquires a spacetime-dependent mass dynamically, a small fraction of which is converted via resonant oscillations into a very weakly mixed dark right-handed neutrino with the observed dark matter relic abundance. Neutrino masses are generated via a traditional two right-handed neutrino type-I seesaw between a third right-handed neutrino and the source neutrino. The gravitational waves produced during the phase transition have a peak frequency that increases with the dark matter mass, and are detectable at future gravitational wave interferometers for dark matter masses in the 1 MeV - 20 GeV range. For source right-handed neutrinos heavier than a GeV, successful leptogenesis is also possible.

The Magnificent CE$\nu$NS Workshop (2018) was held November 2 & 3 of 2018 on the University of Chicago campus and brought together theorists, phenomenologists, and experimentalists working in numerous areas but sharing a common interest in the process of coherent elastic neutrino-nucleus scattering (CE$\nu$NS). This is a collection of abstract-like summaries of the talks given at the meeting, including links to the slides presented. This document and the slides from the meeting provide an overview of the field and a snapshot of the robust CE$\nu$NS-related efforts both planned and underway.

The Antarctic Impulsive Transient Antenna (ANITA) experiment has observed two air shower events with energy $\sim 500~{\rm PeV}$ emerging from the Earth with exit angles $\sim 30^\circ$ above the horizon. As was immediately noted by the ANITA Collaboration, these events (in principle) could originate in the atmospheric decay of an upgoing $\tau$-lepton produced through a charged current interaction of a $\nu_\tau$ inside the Earth. However, the relatively steep arrival angles of these perplexing events create tension with the standard model (SM) neutrino-nucleon interaction cross section. Deepening the conundrum, the IceCube neutrino telescope and the Pierre Auger Observatory with substantially larger exposures to cosmic $\nu_\tau$'s in this energy range have not observed any events. This lack of observation implies that the messenger particle (MP) giving rise to ANITA events must produce an air shower event rate at least a factor of 40 larger than that produced by a flux of $\tau$-neutrinos to avoid conflicts with the upper limits reported by the IceCube and the Pierre Auger collaborations. In addition, the sensitivity of ANITA to MP-induced events must be comparable to or larger than those of IceCube and Auger to avoid conflict with the non-observation of any signal at these facilities. Beyond SM interpretations of ANITA events can be classified according to whether the MPs: (i) live inside the Earth, (ii) originate in neutrino-nucleon collisions inside the Earth, (iii) come from cosmological distances. In this communication we investigate the positive and negative facets of these three classes of models.

The 4.8$\sigma$ anomaly in MiniBooNE data cannot be reconciled with MINOS+ and IceCube data within the vanilla framework of neutrino oscillations involving an eV-mass sterile neutrino. We show that an apparently consistent picture can be drawn if charged-current and neutral-current nonstandard neutrino interactions are at work in the 3+1 neutrino scheme. It appears that either the neutrino sector is more elaborate than usually envisioned, or one or more datasets needs revision.

The DAMPE experiment has recently reported an electron spectrum that can be explained by dark matter annihilation into charged lepton pairs in a nearby dark matter clump. The accompanying bremsstrahlung may yield a gamma-ray excess with a known spectral shape that extends over an angular scale of $O(10^\circ)$. We show that such an excess is not present in Fermi-LAT data.

We present a general, model-independent formalism for determining bounds on the production of photons in dwarf spheroidal galaxies via dark matter annihilation, applicable to any set of assumptions about dark matter particle physics or astrophysics. As an illustration, we analyze gamma-ray data from the Fermi Large Area Telescope to constrain a variety of nonstandard dark matter models, several of which have not previously been studied in the context of dwarf galaxy searches.

We explain the two upgoing ultra-high energy shower events observed by ANITA as arising from the decay in the Earth's interior of the quasi-stable dark matter candidate in the CPT symmetric universe. The dark matter particle is a 480 PeV right-handed neutrino that decays into a Higgs boson and a light Majorana neutrino. The latter interacts in the Earth's crust to produce a tau lepton that in turn initiates an atmospheric upgoing shower. The fact that both events emerge at the same angle from the Antarctic ice-cap suggests an atypical dark matter density distribution in the Earth.

Couplings between light scalar dark matter (DM) and neutrinos induce a perturbation to the neutrino mass matrix. If the DM oscillation period is smaller than ten minutes (or equivalently, the DM particle is heavier than $0.69\times10^{-17}$ eV), the fast-averaging over an oscillation cycle leads to a modification of the measured oscillation parameters. We present a specific $\mu-\tau$ symmetric model in which the measured value of $\theta_{13}$ is entirely generated by the DM interaction, and which reproduces the other measured oscillation parameters. For a scalar DM particle lighter than $10^{-15}$ eV, adiabatic solar neutrino propagation is maintained. A suppression of the sensitivity to CP violation at long-baseline neutrino experiments is predicted in this model. We find that DUNE cannot exclude the DM scenario at more than $3\sigma$ C.L. for bimaximal, tribimaximal and hexagonal mixing, while JUNO can rule it out at more than $6\sigma$ C.L. by precisely measuring both $\theta_{12}$ and $\theta_{13}$.

We consider the prospects for multiple dark matter direct detection experiments to determine if the interactions of a dark matter candidate are isospin-violating. We focus on theoretically well-motivated examples of isospin-violating dark matter (IVDM), including models in which dark matter interactions with nuclei are mediated by a dark photon, a Z, or a squark. We determine that the best prospects for distinguishing IVDM from the isospin-invariant scenario arise in the cases of dark photon- or Z-mediated interactions, and that the ideal experimental scenario would consist of large exposure xenon- and neon-based detectors. If such models just evade current direct detection limits, then one could distinguish such models from the standard isospin-invariant case with two detectors with of order 100 ton-year exposure.

The 6-year dataset of high-energy starting events (HESE) at IceCube indicates a spectrum of astrophysical neutrinos much softer than expected from the Fermi shock acceleration mechanism. On the other hand, IceCube's up-going muon neutrino dataset and Fermi-LAT's gamma-ray spectrum point to an $E^{-2}$ neutrino spectrum. If the HESE data above 200 TeV are fit with the latter flux, an excess at lower energies ensues, which then suggests a multicomponent spectrum. We show that the HESE dataset can be explained by a single $E^{-2}$ power-law neutrino flux from a muon-damped $p\gamma$ source if neutrino interactions are modified by CPT violation. The low-energy excess is naturally explained by the pileup of events from superluminal neutrino decay, and there is no cutoff at high energies due to the contribution of subluminal antineutrinos. The best-fit scenario with CPT violation also predicts the observation of Glashow resonance events in the near future.

High energy neutrinos have been detected by IceCube, but their origin remains a mystery. Determining the sources of this flux is a crucial first step towards multi-messenger studies. In this work we systematically compare two classes of sources with the data: Galactic and extragalactic. We assume that the neutrino sources are distributed according to a class of Galactic models. We build a likelihood function on an event by event basis including energy, event topology, absorption, and direction information. We present the probability that each high energy event with deposited energy $E_{\rm dep}>60$ TeV in the HESE sample is Galactic, extragalactic, or background. For Galactic models considered the Galactic fraction of the astrophysical flux has a best fit value of $1.3\%$ and is $<9.5\%$ at 90\% CL. A zero Galactic flux is allowed at $<1\sigma$.

Measurements of the solar neutrino mass-squared difference from KamLAND and solar neutrino data are somewhat discrepant, perhaps due to nonstandard neutrino interactions in matter. We show that the zenith angle distribution of solar neutrinos at Hyper-Kamiokande and the energy spectrum of reactor antineutrinos at JUNO can conclusively confirm the discrepancy and detect new neutrino interactions.

We consider scenarios in which the annihilation of self-conjugate spin-1 dark matter to a Standard Model fermion-antifermion final state is chirality suppressed, but where this suppression can be lifted by the emission of an additional photon via internal bremsstrahlung. We find that this scenario can only arise if the initial dark matter state is polarized, which can occur in the context of self-interacting dark matter. In particular, this is possible if the dark matter pair forms a bound state that decays to its ground state before the constituents annihilate. We show that the shape of the resulting photon spectrum is the same as for self-conjugate spin-0 and spin-1/2 dark matter, but the normalization is less heavily suppressed in the limit of heavy mediators.

Cosmological parameters deduced from the Planck measurements of anisotropies in the cosmic microwave background are at some tension with direct astronomical measurements of various parameters at low redshifts. Very recently, it has been conjectured that this discrepancy can be reconciled if a certain fraction of dark matter is unstable and decays between recombination and the present epoch. Herein we show that if the superheavy relics have a branching into neutrinos B (X \to \nu \bar \nu) \sim 3 \times 10^{-9}, then this scenario can also accommodate the recently discovered extraterrestrial flux of neutrinos, relaxing the tension between IceCube results and Fermi LAT data. The model is fully predictive and can be confronted with future IceCube data. We demonstrate that in 10 years of observation IceCube will be able to distinguish the mono-energetic signal from X decay at the 3\sigma level. In a few years of data taking with the upgraded IceCube-Gen2 enough statistics will be gathered to elucidate the dark matter--neutrino connection at the 5\sigma level.

We analyze the energy and zenith angle distributions of the latest 2-year IceCube dataset of upward going atmospheric neutrinos to constrain sterile neutrinos at the eV scale in the 3+1 scenario. We find that the parameters favored by a combination of LSND and MiniBooNE data are excluded at more than the 99% C.L. We explore the impact of nonstandard matter interactions on this exclusion, and find that the exclusion holds for nonstandard interactions (NSI) that are within the stringent model-dependent bounds set by collider and neutrino scattering experiments. However, for large NSI parameters subject only to model-independent bounds from neutrino oscillation experiments, the LSND and MiniBooNE data are consistent with IceCube.

We study the scenario of fermionic dark matter that annihilates to standard model fermions through an s-channel axial vector mediator. We point out that the well-known chirality suppression of the annihilation cross section can be alleviated by s-channel internal Higgsstrahlung. The shapes of the cosmic ray spectra are identical to that of t-channel internal Higgsstrahlung in the limit of a heavy mediating particle. Unlike the general case of t-channel bremsstrahlung, s-channel Higgsstrahlung can be the dominant annihilation process even for Dirac dark matter. Since the s-channel mediator can be a standard model singlet, collider searches for the mediator are easily circumvented.

We consider monojet searches at the Large Hadron Collider (LHC) for spin-1 dark matter that interacts with quarks through a contact operator. If the dark matter particles are produced with longitudinal polarizations, then the production matrix element is enhanced by factors of the energy. We show that this particularly effective search strategy can test models for which the energy suppression scale of the operator is as large as 10^5 TeV. As such, these searches can probe a large class of models for which the contact-operator approximation is valid. We find that for contact operators that permit velocity-independent dark matter-nucleon scattering, LHC monojet searches for spin-1 dark matter are competitive with or far surpass direct-detection searches depending on whether the scattering is spin independent or spin dependent, respectively.

We provide a compendium of results at the level of matrix elements for a systematic study of dark matter scattering and annihilation. We identify interactions that yield spin-dependent and spin-independent scattering and specify whether the interactions are velocity- and/or momentum-suppressed. We identify the interactions that lead to s-wave or p-wave annihilation, and those that are chirality-suppressed. We also list the interaction structures that can interfere in scattering and annihilation processes. Using these results, we point out situations in which deviations from the standard lore are obtained.

We study the vacuum stability of a minimal Higgs portal model in which the standard model (SM) particle spectrum is extended to include one complex scalar field and one Dirac fermion. These new fields are singlets under the SM gauge group and are charged under a global U(1) symmetry. Breaking of this U(1) symmetry results in a massless Goldstone boson, a massive CP-even scalar, and splits the Dirac fermion into two new mass-eigenstates, corresponding to Majorana fermions. The lightest Majorana fermion (w) is absolutely stable, providing a plausible dark matter (DM) candidate. We show that interactions between the Higgs sector and the lightest Majorana fermion which are strong enough to yield a thermal relic abundance consistent with observation can easily destabilize the electroweak vacuum or drive the theory into a non-perturbative regime at an energy scale well below the Planck mass. However, we also demonstrate that there is a region of the parameter space which develops a stable vacuum (up to the Planck scale), satisfies the relic abundance, and is in agreement with direct DM searches. Such an interesting region of the parameter space corresponds to DM masses 350 GeV \alt m_w \alt 1 TeV. The region of interest is within reach of second generation DM direct detection experiments.

Now that PeV neutrinos have been discovered by IceCube, we optimistically entertain the possibility that neutrinos with energy above 100 PeV exist. We evaluate the dependence of event rates of such neutrinos on the neutrino-nucleon cross section at observatories that detect particles, atmospheric fluorescence, or Cherenkov radiation, initiated by neutrino interactions. We consider how (i) a simple scaling of the total standard model neutrino-nucleon cross section, (ii) a new elastic neutral current interaction, and (iii) a new completely inelastic interaction, individually impact event rates.

The Glashow resonance at E_\nu=6.3 PeV is a measure of the \bar\nu_e content of the astrophysical neutrino flux. The fractional \bar\nu_e content depends on the neutrino production model at the cosmic neutrino source, and the environment at the source. Thus, the strength of the Glashow resonance event rate is a potential window into astrophysical sources. We quantify the "Glashow resonometer" and comment on the significance that no Glashow events are observed in the IceCube three-year data.

There may be a high-energy cutoff of neutrino events in IceCube data. In particular, IceCube does not observe either continuum events above 2 PeV, or the Standard Model Glashow-resonance events expected at 6.3 PeV. There are also no higher energy neutrino signatures in the ANITA and Auger experiments. This absence of high-energy neutrino events motivates a fundamental restriction on neutrino energies above a few PeV. We postulate a simple scenario to terminate the neutrino spectrum that is Lorentz-invariance violating, but with a limiting neutrino velocity that is always smaller than the speed of light. If the limiting velocity of the neutrino applies also to its associated charged lepton, then a significant consequence is that the two-body decay modes of the charged pion are forbidden above two times the maximum neutrino energy, while the radiative decay modes are suppressed at higher energies. Such stabilized pions may serve as cosmic ray primaries.

We show that S-dual inflationary potentials solve the unlikeliness problem manifested in Planck data and explain the excess B-mode power observed by the BICEP2 experiment as arising from primordial tensor fluctuations.

The anomaly detected by AMS-02 and PAMELA in the cosmic-ray positron flux when interpreted as arising from dark matter annihilation suggests that dark matter may interact differently with hadrons and leptons so as to remain compatible with cosmic-ray antiproton data. Such a scenario is readily accommodated in models with extra spatial dimensions. We study indirect detection of Kaluza-Klein (KK) dark matter in Universal Extra Dimensions with brane-localized terms and fermion bulk masses: Next-to-Minimal Universal Extra Dimensions. So that an excess of antiprotons is not produced in explaining the positron anomaly, it is necessary that the KK bulk masses in the lepton and hadron sectors be distinct. Even so, we find that cosmic-ray data disfavor a heavy KK photon dark matter scenario. Also, we find these scenarios with flavor-universal bulk masses to be in conflict with dijet and dilepton searches at the LHC.

The announcement by the IceCube Collaboration of the observation of 28 cosmic neutrino candidates has been greeted with a great deal of justified excitement. The data reported so far depart by 4.3\sigma from the expected atmospheric neutrino background, which raises the obvious question: "Where in the Cosmos are these neutrinos coming from?" We review the many possibilities which have been explored in the literature to address this question, including origins at either Galactic or extragalactic celestial objects. For completeness, we also briefly discuss new physics processes which may either explain or be constrained by IceCube data.

Isospin-violating dark matter (IVDM) generalizes the standard spin-independent scattering parameter space by introducing one additional parameter, the neutron-to-proton coupling ratio f_n/f_p. In IVDM the implications of direct detection experiments can be altered significantly. We review the motivations for considering IVDM and present benchmark models that illustrate some of the qualitatively different possibilities. IVDM strongly motivates the use of a variety of target nuclei in direct detection experiments.

We explore the sensitivity of an experiment at the Daya Bay site, with a point radioactive source and a few meter baseline, to neutrino oscillations involving one or more eV mass sterile neutrinos. We find that within a year, the entire 3+2 and 1+3+1 parameter space preferred by global fits can be excluded at the 3\sigma level, and if an oscillation signal is found, the 3+1 and 3+2 scenarios can be distinguished from each other at more than the 3\sigma level provided one of the sterile neutrinos is lighter than 0.5 eV.

The nonrelativistic annihilation of Majorana dark matter in the Sun to a pair of light fermions is chirality-suppressed. Annihilation to 3-body final states $\ell^+f^-V$, where $V=W,Z,\gamma$, and $\ell$ and $f$ are light fermions (that may be the same), becomes dominant since bremsstrahlung relaxes the chirality suppression. We evaluate the neutrino spectra at the source, including spin and helicity dependent effects, and assess the detectability of each significant bremsstrahlung channel at IceCube/DeepCore. We also show how to combine the sensitivities to the dark matter-nucleon scattering cross section in individual channels, since typically several channels contribute in models.

Monojet and monophoton final states with large missing transverse energy (${\not E}_T$) are important for dark matter (DM) searches at colliders. We present analytic expressions for the differential cross sections for the parton-level processes, $q\bar{q}(qg)\to g(q)\chi\bar{\chi}$ and $q\bar{q}\to \gamma\chi\bar{\chi}$, for a neutral DM particle with a magnetic dipole moment (MDM) or an electric dipole moment (EDM). We collectively call such DM candidates dipole moment dark matter (DMDM). We also provide monojet cross sections for scalar, vector and axial-vector interactions. We then use ATLAS/CMS monojet${+\not E}_T$ data and CMS monophoton$+{\not E}_T$ data to constrain DMDM. We find that 7 TeV LHC bounds on the MDM DM-proton scattering cross section are about six orders of magnitude weaker than on the conventional spin-independent cross section.

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

We show that the energy spectra from dark matter (DM) annihilation into $f^+f^-\gamma$ (where $f$ is a light fermion) via chirality preserving interactions are identical in 3 scenarios: (a) DM is a Majorana fermion and the particle exchanged is a scalar, (b) DM is a Majorana fermion and the particle exchanged is a vector and (c) DM is a scalar and the particle exchanged is a fermion. For cases (a) and (c), we also calculate the differential cross section to $\ell^+f^-V$, where $V=W,Z,\gamma$, and $\ell$ is a light fermion that may or not be the same as $f$. The form of the cross section depends on whether the DM is a Majorana fermion or a scalar but in both cases its form is independent of $V$.

Searches for dark matter scattering off nuclei are typically compared assuming that the dark matter's spin-independent couplings are identical for protons and neutrons. This assumption is neither innocuous nor well motivated. We consider isospin-violating dark matter (IVDM) with one extra parameter, the ratio of neutron to proton couplings, and include the isotope distribution for each detector. For a single choice of the coupling ratio, the DAMA and CoGeNT signals are consistent with each other and with current XENON constraints, and they unambiguously predict near future signals at XENON and CRESST. We provide a quark-level realization of IVDM as WIMPless dark matter that is consistent with all collider and low-energy bounds.

Data from the LSND and MiniBooNE experiments, and revised expectations of the antineutrino flux from nuclear reactors suggest the existence of eV-mass sterile neutrinos. 3+2 and 1+3+1 scenarios accommodate all relevant short-baseline neutrino data except for the low-energy MiniBooNE anomaly. We analyze the angular distribution of upward going atmospheric neutrino events in the IceCube-40 dataset for evidence of sterile neutrinos within these scenarios. Depending on how systematic uncertainties are handled, we find strong evidence for, or weak evidence against sterile neutrinos. We show that future IceCube data will definitively settle the issue.

We consider the prospects for studying spin-independent isospin-violating dark matter-nucleon interactions with neutrinos from dark matter annihilation in the Sun, with a focus on IceCube/DeepCore (IC/DC). If dark matter-nucleon interactions are isospin violating, IC/DC's reach in the spin-independent cross section may be competitive with current direct detection experiments for a wide range of dark matter masses. We also compare IC/DC's sensitivity to that of next generation argon, germanium, neon and xenon based detectors.

With the augmentation of IceCube by DeepCore, the prospect for detecting dark matter annihilation in the Sun is much improved. To complement this experimental development, we provide a thorough template analysis of the particle physics issues that are necessary to precisely interpret the data. Our study is about nitty-gritty and is intended as a framework for detailed work on a variety of dark matter candidates. To accurately predict the source neutrino spectrum, we account for spin correlations of the final state particles and the helicity-dependence of their decays, and absorption effects at production. We fully treat the propagation of neutrinos through the Sun, including neutrino oscillations, energy losses and tau regeneration. We simulate the survival probability of muons produced in the Earth by using the Muon Monte Carlo program, reproduce the published IceCube effective area, and update the parameters in the differential equation that approximates muon energy losses. To evaluate the zenith-angle dependent atmospheric background event rate, we track the Sun and determine the time it spends at each zenith angle. Throughout, we employ neutralino dark matter as our example.

A neutral dark matter particle may possess an electric dipole moment (EDM) or a magnetic dipole moment (MDM), so that its scattering with nuclei is governed by electromagnetic interactions. If the moments are associated with relevant operators of dimension-5, they may be detectable in direct search experiments. We calculate complete expressions of the scattering cross sections and the recoil energy spectra for dark matter with these attributes. We also provide useful formulae pertinent to dark matter that interacts via an electric charge form factor (CFF) which is related to the charge radius defined by an effective dimension-6 operator. We show that a 7 GeV dark matter particle with an EDM, MDM or CFF easily reproduces the CoGeNT excess while remaining consistent with null searches.

We investigate the prospects for indirect detection of fermion WIMPless dark matter at the neutrino telescopes IceCube and DeepCore. The dark matter annihilating in the Sun is a hidden sector Majorana fermion that couples through Yukawa couplings to a connector particle and a visible sector particle, and it exhibits only spin-dependent scattering with nuclei via couplings to first generation quarks. We consider cases where the annihilation products are taus, staus, or sneutrinos of the three generations. To evaluate the muon fluxes incident at the detector, we propagate the neutrino spectra through the solar medium and to the Earth and account for the effects of neutrino oscillations, energy losses due to neutral- and charged-current interactions, and tau regeneration. We find that for the stau and sneutrino channels, a 5 yr 3$\sigma$ detection of dark matter lighter than about 300 GeV is possible at IceCube for large Yukawa couplings or for dark matter and connector particles with similar masses. The tau channel offers far better detection prospects. However, due to its lower energy threshold and better muon background rejection capability, DeepCore is able to detect signals in all annihilation channels and for a wider range of dark matter masses.

Within a supersymmetric (SUSY) type-I seesaw framework with flavor-blind universal boundary conditions, we study the consequences of requiring that the observed baryon asymmetry of the Universe be explained by either thermal or non-thermal leptogenesis. In the former case, we find that the parameter space is very constrained. In the bulk and stop-coannihilation regions of mSUGRA parameter space (that are consistent with the measured dark matter abundance), lepton flavor-violating (LFV) processes are accessible at MEG and future experiments. However, the very high reheat temperature of the Universe needed after inflation (of about 10^{12} GeV) leads to a severe gravitino problem, which disfavors either thermal leptogenesis or neutralino dark matter. Non-thermal leptogenesis in the preheating phase from SUSY flat directions relaxes the gravitino problem by lowering the required reheat temperature. The baryon asymmetry can then be explained while preserving neutralino dark matter, and for the bulk or stop-coannihilation regions LFV processes should be observed in current or future experiments.

We consider extensions of the standard model based on open strings ending on D-branes, in which gauge bosons and their associated gauginos exist as strings attached to stacks of D-branes, and chiral matter exists as strings stretching between intersecting D-branes. Under the assumptions that the fundamental string scale is in the TeV range and the theory is weakly coupled, we study models of supersymmetry for which signals of annihilating neutralino dark matter are observable. In particular, we construct a model with a supersymmetric R-symmetry violating (but R-parity conserving) effective Lagrangian that allows for the s-wave annihilation of neutralinos, once gauginos acquire mass through an unspecified mechanism. The model yields bino-like neutralinos (with the measured relic abundance) that annihilate to a gamma-gamma final state with a substantial branching fraction (~ 10%) that is orders of magnitude larger than in the minimal supersymmetric standard model. A very bright gamma-ray spectral line could be observed by gamma-ray telescopes.

We study lepton flavor-violating (LFV) processes within a supersymmetric type-I seesaw framework with flavor-blind universal boundary conditions, properly accounting for the effect of the neutrino sector on the dark matter relic abundance. We consider several possibilities for the neutrino Yukawa coupling matrix and show that in regions of SUSY parameter space that yield the correct neutralino relic density, LFV rates can differ from naive estimates by up to two orders of magnitude. Contrary to common belief, we find that current LFV limits do not exclude neutrino Yukawa couplings larger than top Yukawa couplings. We introduce the ISAJET-M program that was used for the computations.

We describe a characteristic signature of dark matter (DM) annihilation or decay into gamma-rays. We show that if the total angular momentum of the initial DM particle(s) vanishes, and helicity suppression operates to prevent annihilation/decay into light fermion pairs, then the amplitude for the dominant 3-body final state f^+f^-\gamma has a unique form dictated by gauge invariance. This amplitude and the corresponding energy spectra hold for annihilation of DM Majorana fermions or self-conjugate scalars, and for decay of DM scalars, thus encompassing a variety of possibilities. Within this scenario, we analyze Fermi LAT, PAMELA and HESS data, and predict a hint in future Fermi gamma-ray data that portends a striking signal at atmospheric Cherenkov telescopes (ACTs).

We analyze new diffuse gamma-ray data from the Fermi Gamma-ray Space Telescope, which do not confirm an excess in the EGRET data at galactic mid-latitudes, in combination with measurements of electron and positron fuxes from PAMELA, Fermi and HESS within the context of three possible sources: dark matter (DM) annihilation or decay into charged leptons, and a continuum distribution of pulsars. We allow for variations in the backgrounds, consider several DM halo profiles, and account for systematic uncertainties in data where possible. We find that all three scenarios represent the data well. The pulsar description holds for a wide range of injection energy spectra. We compare with ATIC data and the WMAP haze where appropriate, but do not fit these data since the former are discrepant with Fermi data and the latter are subject to large systematic uncertainties. We show that for cusped halo profiles, Fermi could observe a spectacular gamma-ray signal of DM annihilation from the galactic center while seeing no excess at mid-latitudes.

Three-dimensional neutral hydrogen mapping using the redshifted 21 cm line has recently emerged as a promising cosmological probe. Within the framework of slow-roll reconstruction, we analyze how well the inflationary potential can be reconstructed by combining data from 21 cm experiments and cosmic microwave background data from the Planck satellite. We consider inflationary models classified according to the amplitude of their tensor component, and show that 21 cm measurements can significantly improve constraints on the slow-roll parameters and determine the shape of the inflationary potential.

Assuming that the positron excess in PAMELA satellite data is a consequence of annihilations of cold dark matter, we consider from a model-independent perspective if the data show a preference for the spin of dark matter. We then perform a general analysis of annihilations into two-body final states to determine what weighted combination of channels best describes the data.

The NSF has chosen the site for the Deep Underground Science and Engineering Laboratory (DUSEL) to be in Lead, South Dakota. In fact, the state of South Dakota has already stepped up to the plate and contributed its own funding for the proposed lab, see http://www.sanfordlaboratoryathomestake.org/index.html. The final decision by NSF for funding the Initial Suite of Experiments for DUSEL will be made early in 2009. At that time the NSF Science Board must make a decision. Of order 200 experimentalists have already expressed an interest in performing experiments at DUSEL. In order to assess the interest of the theoretical community, the Center for Cosmology and Astro-Particle Physics (CCAPP) at The Ohio State University (OSU) organized a 3-day DUSEL Theory Workshop in Columbus, Ohio from April 4 - 6, 2008. The workshop focused on the scientific case for six proposed experiments for DUSEL: long baseline neutrino oscillations, proton decay, dark matter, astrophysical neutrinos, neutrinoless double beta decay and N-Nbar oscillations. The outcome of this workshop is the DUSEL Theory White paper addressing the scientific case at a level which may be useful in the decision making process for policy makers at the NSF and in the U.S. Congress. In order to assess the physics interest in the DUSEL project we have posted the DUSEL Theory White paper on the following CCAPP link http://ccapp.osu.edu/whitepaper.html . Please read the white paper and, if you are interested, use the link to show your support by co-signing the white paper.

Motivated by the fact that neutrinos are massive, we study the effect of neutrino Yukawa couplings on neutralino dark matter observables within the framework of a supersymmetric seesaw. We find that neutrino couplings significantly affect the neutralino relic density in regions of parameter space where soft SUSY-breaking slepton masses and/or trilinear couplings are large. Depending on the size of the couplings, the neutralino relic density spans over an order of magnitude in the A-funnel, focus point and stop-coannihilation regions of mSUGRA. We also show that dark matter detection rates can be modified by up to several orders of magnitude.

We show that the rate for upward showers from an isotropic cosmic neutrino flux at neutrino telescopes like IceCube is independent of the neutrino-nucleon cross section. For bins that span a relatively narrow range in energy, neither scaling the cross section, nor changing its power-law energy behavior affects the upward shower rate, which depends only on the flux. The neutrino flux can be completely known since its spectral shape can be determined by comparing the rates in neighboring bins. We also show that the downward shower rate varies linearly with cross section with a proportionality constant determined by the energy-dependence of the cross section, independent of the power-law behavior of the flux. The normalization and energy dependence of the cross section can be known by comparing the downward rates in neighboring bins.

Within the context of mass-varying neutrinos, we construct a cosmological model that has a phase transition in the electromagnetic fine structure constant \alpha at a redshift of 0.5. The model accommodates hints of a time variable \alpha in quasar spectra and the nonobservance of such an effect at very low redshifts. It is consistent with limits from the recombination and primordial nucleosynthesis eras and is free of instabilities.

We test if the latest Gold set of 182 SNIa or the combined "Platinum" set of 192 SNIa from the ESSENCE and Gold sets, in conjunction with the CMB shift parameter show a preference between the LambdaCDM model, three wCDM models, and the DGP model of modified gravity as an explanation for the current accelerating phase of the universe's expansion. We consider flat wCDM models with an equation of state w(a) that is (i) constant with scale factor $a$, (ii) varies as w(a)=w_0+w_a(1-a) for redshifts probed by supernovae but is fixed at -1 at earlier epochs and (iii) varies as w_0+w_a(1-a) since recombination. We find that all five models explain the data with comparable success.

We examine the dependence of event rates at neutrino telescopes on the neutrino-nucleon cross section for neutrinos with energy above 1 PeV, and contrast the results with those for cosmic ray experiments. Simple scaling of the Standard Model cross sections leaves the rates of upward events essentially unchanged. Details, such as detector depth and cross section inelasticity, can influence rates. Numerical estimates of upward shower, muon and tau event rates in the IceCube detector confirm these results.

We study the proposal that mass-varying neutrinos could provide an explanation for the LSND signal for \bar\nu_mu to \bar\nu_e oscillations. We first point out that all positive oscillation signals occur in matter and that three active mass-varying neutrinos are insufficient to describe all existing neutrino data including LSND. We then examine the possibility that a model with four mass-varying neutrinos (three active and one sterile) can explain the LSND effect and remain consistent with all other neutrino data. We find that such models with a 3+1 mass structure in the neutrino sector may explain the LSND data and a null MiniBooNE result for 0.10 < \sin^2 2\theta_x < 0.30. Predictions of the model include a null result at Double-CHOOZ, but positive signals for underground reactor experiments and for \nu_\mu to \nu_e oscillations in long-baseline experiments.

We deduce the cosmogenic neutrino flux by jointly analysing ultra high energy cosmic ray data from HiRes-I and II, AGASA and the Pierre Auger Observatory. We make two determinations of the neutrino flux by using a model-dependent method and a model-independent method. The former is well-known, and involves the use of a power-law injection spectrum. The latter is a regularized unfolding procedure. We then use neutrino flux bounds obtained by the RICE experiment to constrain the neutrino-nucleon inelastic cross section at energies inaccessible at colliders. The cross section bounds obtained using the cosmogenic fluxes derived by unfolding are the most model-independent bounds to date.

We classify dark energy models in a plane of observables that correspond to the common parameterization of a non-constant equation of state, w(a)=w_0 + w_a(1-a), where $a$ is the scale factor of the universe. The models fall into four classes and only two of these classes have a region of overlap in the observable plane. We perform a joint analysis of all Type Ia supernova (SNIa) data compiled by the High-Z SN Search Team (HZT) and the Supernova Legacy Survey (SNLS) and find that no class of models is excluded by current SNIa data. However, an analysis of large scale structure, Ly$\alpha$ forest and bias constraints from SDSS, the Gold SNIa data and WMAP data indicates that non-phantom barotropic models with a a positive sound speed are excluded at the 95% C.L.

We propose that the solar neutrino deficit may be due to oscillations of mass-varying neutrinos (MaVaNs). This scenario elucidates solar neutrino data beautifully while remaining comfortably compatible with atmospheric neutrino and K2K data and with reactor antineutrino data at short and long baselines (from CHOOZ and KamLAND). We find that the survival probability of solar MaVaNs is independent of how the suppression of neutrino mass caused by the acceleron-matter couplings varies with density. Measurements of MeV and lower energy solar neutrinos will provide a rigorous test of the idea.

We demonstrate that the detection of shock modulations of the neutrino spectra from a galactic core-collapse supernova is sufficient to obtain a high significance determination of the neutrino mass hierarchy if the supernova event is observed in both a Mton-class water Cherenkov detector and a 100 kton-class liquid argon detector. Neither detailed supernova neutrino flux modelling nor observation of Earth matter effects is needed for this determination. As a corollary, a nonzero value of \theta_x will be established.

We perform a comparative study of two methods of determining the survival probabilities of low, intermediate, and high energy solar neutrinos that emphasizes the general agreement between the Large Mixing Angle (LMA) solution and extant solar neutrino data. The first analysis is oscillation parameter-independent and the second analysis involves an approximate calculation of the survival probabilities in the three energy ranges that depends only on oscillation parameters. We show that future experiments like BOREXino, CLEAN, Heron, LENS and MOON, that measure $pp$ and $^7$Be neutrinos, will facilitate a stringent test of the LMA solution independently of the Standard Solar Model (SSM), without recourse to earth-matter effects. Throughout, we describe the role of SSM assumptions on our results. If the LMA solution passes the test without needing to be modified, it may be possible to establish that $\theta_x$ is nonzero at more than $2\sigma$ assuming the SSM prediction for the $pp$ flux is correct.

We investigate whether cosmological data suggest the need for massive neutrinos. We employ galaxy power spectrum measurements from the Sloan Digital Sky Survey (SDSS) and the Two Degree Field Galaxy Redshift Survey (2dFGRS), along with cosmic microwave background (CMB) data from the Wilkinson Microwave Anisotropy Probe (WMAP) and 27 other CMB experiments. We also use the measurement of the Hubble parameter from the Hubble Space Telescope (HST) Key Project. We find the sum of the neutrino masses to be smaller than 0.75 eV at 2\sigma (1.1 eV at 3\sigma).

The current status of the physics of massive neutrinos is reviewed with a forward-looking emphasis. The article begins with the general phenomenology of neutrino oscillations in vacuum and matter and documents the experimental evidence for oscillations of solar, reactor, atmospheric and accelerator neutrinos. Both active and sterile oscillation possibilities are considered. The impact of cosmology (BBN, CMB, leptogenesis) and astrophysics (supernovae, highest energy cosmic rays) on neutrino observables and vice versa, is evaluated. The predictions of grand unified, radiative and other models of neutrino mass are discussed. Ways of determining the unknown parameters of three-neutrino oscillations are assessed, taking into account eight-fold degeneracies in parameters that yield the same oscillation probabilities, as well as ways to determine the absolute neutrino mass scale (from beta-decay, neutrinoless double-beta decay, large scale structure and Z-bursts). Critical unknowns at present are the amplitude of \nu_\mu to \nu_e oscillations and the hierarchy of the neutrino mass spectrum; the detection of CP violation in the neutrino sector depends on these and on an unknown phase. The estimated neutrino parameter sensitivities at future facilities (reactors, superbeams, neutrino factories) are given. The overall agenda of a future neutrino physics program to construct a bottom-up understanding of the lepton sector is presented.

We quantify the extent to which extra relativistic energy density can be concealed by a neutrino asymmetry without conflicting with the baryon asymmetry measured by the Wilkinson Microwave Anisotropy Probe (WMAP). In the presence of a large electron neutrino asymmetry, slightly more than seven effective neutrinos are allowed by Big Bang Nucleosynthesis (BBN) and WMAP at 2\sigma. The same electron neutrino degeneracy that reconciles the BBN prediction for the primordial helium abundance with the observationally inferred value also reconciles the LSND neutrino with BBN by suppressing its thermalization prior to BBN.

We assay how inflationary models whose properties are dominated by the dynamics of a single scalar field are constrained by cosmic microwave background (CMB) data from the Wilkinson Microwave Anisotropy Probe (WMAP). We classify inflationary models in a plane defined by the horizon-flow parameters. Our approach differs from that of the WMAP collaboration in that we analyze only WMAP data and take the spectral shapes from slow-roll inflation rather than power-law parameterizations of the spectra. The only other information we use is the measurement of $h$ from the Hubble Space Telescope (HST) Key Project. We find that the spectral index of primordial density perturbations lies in the 1\sigma range 0.94<n_s<1.04 with no evidence of running. The ratio of the amplitudes of tensor and scalar perturbations is smaller than 0.61 and the inflationary scale is below 2.8\times 10^{16} GeV, both at the 2\sigma C.L. No class of inflation or ekpyrotic/cyclic models is excluded. The \lambda \phi^4 potential is excluded at 3\sigma only if the number of e-folds is assumed to be less than 45.

We place constraints on the number of relativistic degrees of freedom and on the baryon asymmetry at the epoch of Big Bang Nucleosynthesis (BBN) and at recombination, using cosmic background radiation (CBR) data from the Wilkinson Microwave Anisotropy Probe (WMAP), complemented by the Hubble Space Telescope (HST) Key Project measurement of the Hubble constant, along with the latest compilation of deuterium abundances and measurements of the primordial helium abundance. The agreement between the derived values of these key cosmological and particle physics parameters at these widely separated (in time or redshift) epochs is remarkable. From the combination of CBR and BBN data, we find the 2\sigma ranges for the effective number of neutrinos and for the baryon asymmetry (baryon to photon number ratio \eta) to be 1.7-3.0 and 5.53-6.76 \times 10^{-10}, respectively.

The KamLAND reactor antineutrino experiment has detected a 3.4\sigma flux suppression relative to the expectation if no neutrino oscillations occur. We combine KamLAND data with solar neutrino data and show that the LMA solution is the only viable oscillation solution to the solar neutrino problem at the 4.4\sigma C. L.

We show that the temperatures of the emergent non-electron neutrinos and the binding energy released by a galactic Type II supernova are determinable, assuming the Large Mixing Angle (LMA) solution is correct, from observations at the Sudbury Neutrino Observatory (SNO) and at Super-Kamiokande (SK). If the neutrino mass hierarchy is inverted, either a lower or upper bound can be placed on the neutrino mixing angle $\theta_{13}$, and the hierarchy can be deduced for adiabatic transitions. For the normal hierarchy, neither can $\theta_{13}$ be constrained nor can the hierarchy be determined. Our conclusions are qualitatively unchanged for the proposed Hyper-Kamiokande detector.

We perform a global analysis in the framework of two active neutrino oscillations of all solar neutrino data, including the recent SNO day and night spectra (comprised of the charged current (CC), elastic scattering (ES) and neutral current (NC) events), the Super-Kamiokande (SK) day and night spectra (from 1496 days) and the updated SAGE results. We find that the Large Mixing Angle (LMA) solution is selected at the 99% C.L.; the best-fit parameters are \Delta m^2=5.6 \times 10^{-5} eV^2 and \theta=32^{\circ}. No solutions with \theta\geq \pi/4 are allowed at the 5\sigma C.L. Oscillations to a pure sterile state are excluded at 5.3\sigma, but a sizeable sterile neutrino component could still be present in the solar flux.

We examine how constraints can be placed on the neutrino component of dark matter by an accurate measurement of neutrinoless double beta ($0\nu\beta\beta$) decay and the solar oscillation amplitude. We comment on the alleged evidence for $0\nu\beta\beta$ decay.

We dispel the misconception that data from SN 1987A favor the normal neutrino mass hierarchy over the inverted hierarchy for $\sin^2 \theta_{13} \gsim 10^{-4}$. We find comparable fits for the two hierarchies. No bound can be placed on the mixing angle $\theta_{13}$ even at the 1$\sigma$ level.

We show that the combined data from solar neutrino experiments and from the KamLAND reactor neutrino experiment can establish an upper limit on, or detect, potential CPT violation in the neutrino sector of order 10^{-20} GeV to 10^{-21} GeV.

We analyze the 1258-day Super-Kamiokande day and night solar neutrino energy spectra with various $\chi^2$ definitions. The best-fit lies in the LMA region at $(\Delta m^2, tan^2 \theta)=(5.01\times 10^{-5} eV^2, 0.60)$, independently of whether systematic errors are included in the $\chi^2$-definition. We compare the exclusion and allowed regions from the different definitions and choose the most suitable definition to predict the regions from SNO at the end of three years of data accumulation. We first work under the assumption that Super-Kamiokande sees a flux-suppressed flat energy spectrum. Then, we consider the possibility of each one of the three MSW regions being the solution to the solar neutrino problem. We find that the exclusion and allowed regions for the flat spectrum hypothesis and the LMA and LOW solutions are alike. In three years, we expect SNO to find very similar regions to that obtained by Super-Kamiokande. We evaluate whether the zenith angle distribution at SNO with optimum binning will add anything to the analysis of the day and night spectra; for comparison, we show the results of our analysis of the 1258-day zenith angle distribution from Super-Kamiokande, for which the best-fit parameters are $(\Delta m^2, tan^2 \theta)=(5.01\times 10^{-5} eV^2, 0.56)$.

We study how well KamLAND, the first terrestrial neutrino experiment capable of addressing the solar neutrino problem, will perform in ascertaining whether or not the large mixing angle MSW solution (with $10^{-5}<\Delta m_{21}^2< 10^{-4} eV^2$ and oscillation amplitude $sin^2 2 \theta_{12}>0.3$) is correct. We find that in a year of operation KamLAND will provide unequivocal evidence for or against this solution. Furthermore, its sensitivity to the three-neutrino oscillation parameters in this region is sufficiently acute as to determine $\Delta m_{21}^2$ to approximately $\pm 10$ % (for $sin^2 2 \theta_{12}>0.7$) and to fix $sin^2 2 \theta_{12}$ to within $\pm 0.1$ (at the $2\sigma$ level) with three years of accumulated data, independent of the value of $\theta_{13}$.

We consider the efficacy of using luminosity distance measurements of deep redshift supernovae to discriminate between two forms of dark energy, quintessence (a scalar field with canonical kinetic terms rolling down a potential) and k-essence (a scalar field whose cosmic evolution is driven entirely by non-linear kinetic terms). The primary phenomenological distinction between the two types of models that can be quantified by supernova searches (at least in principle) is that the equation of state $w\equiv p/\rho$ of quintessence is falling today while that of k-essence is rising. By simulating $10^5$ possible datasets that SNAP could obtain, we show that even if the mass density $\Omega_m$ is known exactly, an ambiguity remains that may not allow a definitive distinction to be made between the two types of theories.