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.

Rather general considerations of the string theory landscape imply a mild statistical draw towards large soft SUSY breaking terms tempered by the requirement of proper electroweak symmetry breaking where SUSY contributions to the weak scale are not too far from m(weak)~ 100 GeV. Such a picture leads to the prediction that m_h~ 125 GeV while most sparticles are beyond current LHC reach. Here we explore the possibility that the magnitude of the Peccei-Quinn (PQ) scale f_a is also set by string landscape considerations within the framework of a compelling SUSY axion model. First, we examine the case where the PQ symmetry arises as an accidental approximate global symmetry from a more fundamental gravity-safe Z(24)^R symmetry and where the SUSY mu parameter arises from a Kim-Nilles operator. The pull towards large soft terms then also pulls the PQ scale as large as possible. Unless this is tempered by rather severe (unknown) cosmological or anthropic bounds on the density of dark matter, then we would expect a far greater abundance of dark matter than is observed. This conclusion cannot be negated by adopting a tiny axion misalignment angle theta_i because WIMPs are also overproduced at large f_a. Hence, we conclude that setting the PQ scale via anthropics is highly unlikely. Instead, requiring soft SUSY breaking terms of order the gravity-mediation scale m_{3/2}~ 10-100 TeV places the mixed axion-neutralino dark matter abundance into the intermediate scale sweet zone where f_a~ 10^{11}-10^{12} GeV. We compare our analysis to the more general case of a generic SUSY DFSZ axion model with uniform selection on theta_i but leading to the measured dark matter abundance: this approach leads to a preference for f_a~ 10^{12} GeV.

The binary black hole merger events observed by the Advanced LIGO (aLIGO) and VIRGO collaboration can shed light on the origins of black holes. Many studies based on black hole stellar origins have shown a maximum mass for stellar black holes, which can be measured or constrained from the observed black hole mass distribution. In this paper, we point out that the redshift distribution of the observed merger events can provide complementary information for studying the black hole mass distribution, because the detectability correlates the event redshift to the black hole masses. Based on the five observed events and using the Kolmogorov-Smirnov test, we have found that the maximum of the stellar black hole masses are constrained to be below $76\,M_\odot$ at 90\% confidence level, for a negative power-law index of 2.3 for the heavier black hole. With the improved sensitivity of aLIGO, a few dozen merger events may be obtained, for which we estimate that the maximum mass will be constrained to $10\,M_\odot$ accuracy.

The string theory landscape of vacua solutions provides physicists with some understanding as to the magnitude of the cosmological constant. Similar reasoning can be applied to the magnitude of the soft SUSY breaking terms in supersymmetric models of particle physics: there appears to be a statistical draw towards large soft terms which is tempered by the anthropic requirement of the weak scale lying not too far from ~100 GeV. For a mild statistical draw of m_{soft}^n with n=1 (as expected from SUSY breaking due to a single F term) then the light Higgs mass is preferred at ~125 GeV while sparticles are all pulled beyond LHC bounds. We confront a variety of LHC and WIMP dark matter search limits with the statistical expectations from a fertile patch of string theory landscape. The end result is that LHC and WIMP dark matter detectors see exactly that which is expected from the string theory landscape: a Standard Model-like Higgs boson of mass 125 GeV but as yet no sign of sparticles or WIMP dark matter. SUSY from the n=1 landscape is most likely to emerge at LHC in the soft opposite-sign dilepton plus jet plus MET channel. Multi-ton noble liquid WIMP detectors should be able to completely explore the n=1 landscape parameter space.

Particle physics models with Peccei-Quinn (PQ) symmetry breaking as a consequence of supersymmetry (SUSY) breaking are attractive in that they solve the strong CP problem with a SUSY DFSZ-like axion, link the SUSY breaking and PQ breaking intermediate mass scales and can resolve the SUSY mu problem with a naturalness-required weak scale mu term whilst soft SUSY breaking terms inhabit the multi-TeV regime as required by LHC sparticle mass limits and the Higgs mass measurement. On the negative ledger, models based on global symmetries suffer a generic gravity spoliation problem. We present a model based on the discrete R-symmetry Z_{24}^R-- which may emerge from compactification of 10-d Lorentzian spacetime in string theory-- where the mu term and dangerous proton decay and R-parity violating operators are either suppressed or forbidden while a gravity-safe PQ symmetry emerges as an accidental approximate global symmetry leading to a solution to the strong CP problem and a weak-scale/natural value for the mu term.

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.

The requirement of electroweak naturalness in supersymmetric (SUSY) models of particle physics necessitates light higgsinos not too far from the weak scale characterized by m(weak)~ m(W,Z,h)~100 GeV. On the other hand, LHC Higgs mass measurements and sparticle mass limits point to a SUSY breaking scale in the multi-TeV regime. Under such conditions, the lightest SUSY particle is expected to be a mainly higgsino-like neutralino with non-negligible gaugino components (required by naturalness). The computed thermal WIMP abundance in natural SUSY models is then found to be typically a factor 5-20 below its measured value. To gain concordance with observations, either an additional DM particle (the axion is a well-motivated possibility) must be present or additional non-thermal mechanisms must augment the neutralino abundance. We compare present direct and indirect WIMP detection limits to three natural SUSY models based on gravity-, anomaly- and mirage-mediation. We show that the case of natural higgsino-only dark matter where non-thermal production mechanisms augment its relic density, is essentially excluded by a combination of direct detection constraints from PandaX-II, LUX and Xenon-1t experiments, and by bounds from Fermi-LAT/MAGIC observations of gamma rays from dwarf spheroidal galaxies.

We argue that if ultrahigh-energy (E > 10^10 GeV) cosmic rays are heavy nuclei (as indicated by existing data), then the pointing of cosmic rays to their nearest extragalactic sources is expected for 10^10.6 < E/GeV < 10^11. This is because for a nucleus of charge Ze and baryon number A, the bending of the cosmic ray decreases as Z/E with rising energy, so that pointing to nearby sources becomes possible at highest energies. In addition, the maximum energy of acceleration capability of the sources grows linearly in Z, while the energy loss per distance traveled decreases with increasing A. Each of these two points tend to favor heavy nuclei at the highest energies. The traditional bi-dimensional analyses, which simultaneously reproduce Auger data on the spectrum and nuclear composition, may not be capable of incorporating the relative importance of these phenomena. In this paper we propose a multi-dimensional reconstruction of the individual emission spectra (in E, direction, and cross-correlation with nearby putative sources) to study the hypothesis that primaries are heavy nuclei subject to GZK photo-disintegration, and to determine the nature of the extragalactic sources. More specifically, we propose to combine information on nuclear composition and arrival direction to associate a potential clustering of events with a 3-dimensional position in the sky. Actually, both the source distance and maximum emission energy can be obtained through a multi-parameter likelihood analysis to accommodate the observed nuclear composition of each individual event in the cluster. We show that nucleus-emitting-sources exhibit a cepa stratis structure on Earth which could be pealed off by future space-missions, such as POEMMA. We also show that metal-rich starburst galaxies are highly-plausible candidate sources, and we use them as an explicit example of our proposed multi-dimensional analysis.

Deep inelastic scattering of very high-energy neutrinos can potentially be enhanced by the production of a single top quark or charm quark via the interaction of a virtual $W$-boson exchange with a $b$-quark or $s$-quark parton in the nucleon. The single top contribution shows a sharp rise at neutrino energies above 0.5 PeV and gives a cross-section contribution of order 5 percent at 10 PeV, while single charm has a low energy threshold and contributes about 25 percent. Semi-leptonic decays of top and charm give di-muon events whose kinematic characteristics are shown. The angular separation of the di-muons from heavy quark production in the IceCube detector can reach up to one degree. Top quark production has a unique, but rare, three muon signal.

As a cold dark matter candidate, the QCD axion may form Bose-Einstein condensates, called axion stars, with masses around $10^{-11}\,M_{\odot}$. In this paper, we point out that a brand new astrophysical object, a Hydrogen Axion Star (HAS), may well be formed by ordinary baryonic matter becoming gravitationally bound to an axion star. We study the properties of the HAS and find that the hydrogen cloud has a high pressure and temperature in the center and is likely in the liquid metallic hydrogen state. Because of the high particle number densities for both the axion star and the hydrogen cloud, the feeble interaction between axion and hydrogen can still generate enough internal power, around $10^{13}~\mbox{W}\times(m_a/5~\mbox{meV})^4$, to make these objects luminous point sources. High resolution ultraviolet, optical and infrared telescopes can discover HAS via black-body radiation.

We examine updated prospects for detecting WIMPs in supersymmetric models via direct and indirect dark matter search experiments. We examine several historical and also still viable scenarios: projections for well-tempered neutralinos (WTN), projections from the MasterCode (MC), BayesFits (BF) and Fittino (FO) collaborations, non-thermal wino dark matter (NThW) and finally mixed axion-higgsino dark matter from SUSY with radiatively-driven naturalness (RNS). The WTN is ruled out by recent limits from XENON and LUX collaborations. The NThW scenario, previously on tenuous ground due to gamma-line searches, appears also ruled out by recent combined Fermi-LAT/MAGIC limits combined with new HESS results from continuum gamma rays. Substantial portions of MC parameter space and 1 TeV higgsino parameter space from BF group are ruled out. The 100-300 GeV higgsino-like WIMP from RNS survives due to its possible depleted local abundance (where the axion may make up the bulk of dark matter). Projections from ton-scale noble liquid detectors should discover or rule out WIMPs from the remaining parameter space of these surviving models.

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 study the contribution of Galactic sources to the flux of astrophysical neutrinos recently observed by the IceCube Collaboration. We show that in the simplest model of homogeneous and isotropic cosmic ray diffusion in the Milky Way the Galactic diffuse neutrino emission consistent with $\gamma$-ray (Fermi-LAT) and cosmic ray data (KASCADE, KASCADE-Grande and CREAM) is expected to account for only $4\%-8\%$ of the IceCube flux above 60 TeV. Direct neutrino emission from cosmic ray-gas ($pp$) interactions in the sources would require an unusually large average opacity above 0.01. On the other hand, we find that the IceCube events already probe Galactic neutrino scenarios via the distribution of event arrival directions. Based on the latter, we show that most Galactic scenarios can only have a limited contribution to the astrophysical signal: diffuse Galactic emission ($\lesssim50\%$), quasi-diffuse emission of neutrino sources ($\lesssim65\%$), extended diffuse emission from the Fermi Bubbles ($\lesssim25\%$) or unidentified TeV $\gamma$-ray sources ($\lesssim25\%$). The arguments discussed here leave, at present, dark matter decay unconstrained.

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.

By insisting on naturalness in both the electroweak and QCD sectors of the MSSM, the portrait for dark matter production is seriously modified from the usual WIMP miracle picture. In SUSY models with radiatively-driven naturalness (radiative natural SUSY or RNS) which include a DFSZ-like solution to the strong CP and SUSY mu problems, dark matter is expected to be an admixture of both axions and higgsino-like WIMPs. The WIMP/axion abundance calculation requires simultaneous solution of a set of coupled Boltzmann equations which describe quasi-stable axinos and saxions. In most of parameter space, axions make up the dominant contribution of dark matter although regions of WIMP dominance also occur. We show the allowed range of PQ scale f_a and compare to the values expected to be probed by the ADMX axion detector in the near future. We also show WIMP detection rates which are suppressed from usual expectations because now WIMPs comprise only a fraction of the total dark matter. Nonetheless, ton-scale noble liquid detectors should be able to probe the entirety of RNS parameter space. Indirect WIMP detection rates are less propitious since they are reduced by the square of the depleted WIMP abundance.

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 investigate whether a subset of high-energy events observed by IceCube may be due to neutrinos from Sagittarius A*. We check both spatial and temporal coincidences of IceCube events with other transient activities of Sagittarius A*. Among the seven IceCube shower events nearest to the galactic center, we have found that event 25 has a time very close to (around three hours after) the brightest X-ray flare of Sagittarius A* observed by the Chandra X-ray Observatory with a p-value of 0.9%. Furthermore, two of the seven events occurred within one day of each other (there is a 1.6% probability that this would occur for a random distribution in time). Thus, the determination that some IceCube events occur at similar times as X-ray flares and others occur in a burst could be the smoking gun that Sagittarius A* is a point source of very high energy neutrinos. We point out that if IceCube galactic center neutrino events originate from charged pion decays, then TeV gamma rays should come from neutral pion decays at a similar rate. We show that the CTA, HAWC, H.E.S.S. and VERITAS experiments should be sensitive enough to test this.

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.

Recent null results from LHC8 SUSY searches along with the discovery of a SM-like Higgs boson with mass m(h) ~ 125.5 GeV indicates sparticle masses in the TeV range, causing tension with conventional measures of electroweak fine-tuning. We propose a simple Fine-tuning Rule which should be followed under any credible evaluation of fine-tuning. We believe that overestimates of electroweak fine-tuning by conventional measures all arise from violations of this rule. We show that to gain accord with the Fine-tuning Rule, then both Higgs mass and the traditional \Delta_{BG} fine-tuning measures reduce to the model-independent electroweak fine-tuning measure \Delta_{EW}. This occurs by combining dependent contributions to m(Z) or m(h) into independent units. Then, using \Delta_{EW}, we evaluate EW fine-tuning for a variety of SUSY models including mSUGRA, NUHM1, NUHM2, mGMSB, mAMSB, hyper-charged AMSB and nine cases of mixed moduli-anomaly (mirage) mediated SUSY breaking models (MMAMSB) whilst respecting LHC Higgs mass and B-decay constraints (we do not impose LHC8 sparticle mass constraints due to the possibility of compressed spectra within many of these models). We find mSUGRA, mGMSB, mAMSB and MMAMSB models all to be highly fine-tuned. The NUHM1 model is moderately fine-tuned while NUHM2 which allows for radiatively-driven naturalness (RNS) allows for fine-tuning at a meager 10% level in the case where m(higgsinos) ~ 100-200 GeV and the TeV-scale top squarks are well-mixed. Models with RNS may or may not be detect at LHC14. A \sqrt{s} ~ 500 GeV e^+e^- collider will be required to make a definitive search for the requisite light higgsinos.

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.

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

We interpret the PeV shower events observed by the IceCube collaboration as an s-channel enhancement of neutrino-quark scattering by a leptoquark that couples to the tau-flavor and light quarks. With a leptoquark mass of ~ 0.6 TeV and a steep 1/E^{2.3} neutrino flux, charged-current scattering gives shower events at ~ 1 PeV and neutral-current scattering gives shower events at ~ 0.5 PeV. This mechanism is also consistent with the paucity of muon-track events above 100 TeV.

Supersymmetric models with low electroweak fine-tuning contain light higgsinos with mass not too far from m_h\simeq 125 GeV, while other sparticles can be much heavier. In the R-parity conserving MSSM, the lightest neutralino is then a higgsino-like WIMP (albeit with non-negligible gaugino components), with thermal relic density well below measured values. This leaves room for axions (or other, perhaps not as well motivated, stable particles) to function as co-dark matter particles. The local WIMP abundance is then expected to be below standard estimates, and direct and indirect detection rates must be accordingly rescaled. We calculate rescaled direct and indirect higgsino-like WIMP detection rates in SUSY models that fulfil the electroweak naturalness condition. In spite of the rescaling, we find that ton-scale noble liquid detectors can probe the entire higgsino-like WIMP parameter space, so that these experiments should either discover WIMPs or exclude the concept of electroweak naturalness in R-parity conserving natural SUSY models. Prospects for spin-dependent or indirect detection are more limited due in part to the rescaling effect.

Supersymmetric models which fulfill the conditions of electroweak naturalness generally contain light higgsinos with mass not too far from m_h ~125 GeV, while other sparticles can be much heavier. In R-parity conserving models, the lightest neutralino is then a higgsino-like WIMP (albeit with non-negligible gaugino components), with thermal relic density well below measured values. This leaves room for axions to function as co-dark matter particles. The local WIMP abundance is then expected to be below standard estimates, and direct and indirect detection rates must be accordingly rescaled. We calculate rescaled direct and indirect higgsino-like WIMP detection rates in SUSY models that fulfil the electroweak naturalness condition. In spite of the rescaling, we find that ton-scale noble liquid detectors can probe the entire higgsino-like WIMP parameter space, so that these experiments should either discover WIMPs or exclude the concept of electroweak naturalness in R-parity conserving natural SUSY models. Prospects for spin-dependent or indirect detection are more limited due in part to the rescaling effect.

We propose an interpretation of the two neutrino initiated cascade events with PeV energies observed by IceCube: Ultra-high energy cosmic ray protons (or Fe nuclei) scatter on CMB photons through the Delta-resonance (the Berezinsky-Zatsepin process) yielding charged pions and neutrons. The neutron decays give electron-antineutrinos which undergo neutrino oscillations to populate all antineutrino flavors, but the electron-antineutrino flux remains dominant. At 6.3 PeV electron-antineutrino energy their annihilation on electrons in the IceCube detector is enhanced by the Glashow resonance (the W-boson) whose decays can give the PeV showers observed in the IceCube detector. The two observed showers with ~1 PeV energies would need to be from W leptonic decays to electrons and taus. An order of magnitude higher event rate of showers at 6.3 PeV is predicted from W to hadron decays. This interpretation can be tested in the near term. It has significant physics implications on the origin of the highest energy cosmic rays, since neutrino events and cosmic ray events likely share a common origin.

We study a two Higgs doublet model augmented by a scalar dark matter particle that provides an excellent fit to the LHC Higgs data and the Fermi-LAT 135 GeV line. The heavy CP-even Higgs boson, which predominantly mediates annihilation and scattering, must have a coupling to weak gauge bosons at or below percent level to suppress the continuum gamma-ray spectrum below the limit from the Fermi-LAT data and the anti-proton spectrum constrained by the PAMELA data. Discovering or excluding this CP-even Higgs boson at the LHC with a mass between 265 and 280 GeV and an enhanced diphoton branching ratio is crucial to test this scenario.

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.

The ATLAS and CMS collaborations have reported an excess of events in the \gamma\gamma, ZZ^*\to 4\ell and WW^* search channels at an invariant mass m \simeq 125 GeV, which could be the first evidence for the long-awaited Higgs boson. We investigate the consequences of requiring m_h\simeq 125 GeV in both the mSUGRA and NUHM2 SUSY models. In mSUGRA, large values of trilinear soft breaking parameter |A_0| are required, and universal scalar m_0\agt 0.8 TeV is favored so that we expect squark and slepton masses typically in the multi-TeV range. This typically gives rise to an "effective SUSY" type of sparticle mass spectrum. In this case, we expect gluino pair production as the dominant sparticle creation reaction at LHC. For m_0< 5 TeV, the superpotential parameter \mu > 2 TeV and m_A> 0.8 TeV, greatly restricting neutralino annihilation mechanisms. These latter conclusions are softened if m_0\sim 10-20 TeV or if one proceeds to the NUHM2 model. The standard neutralino abundance tends to be far above WMAP-measured values unless the neutralino is higgsino-like. We remark upon possible non-standard (but perhaps more attractive) cosmological scenarios which can bring the predicted dark matter abundance into accord with the measured value, and discuss the implications for direct and indirect detection of neutralino cold dark matter.

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

The minimal supergravity (mSUGRA or CMSSM) model is an oft-used framework for exhibiting the properties of neutralino (WIMP) cold dark matter (CDM). However, the recent evidence from Atlas and CMS on a light Higgs scalar with mass m_h\simeq 125 GeV highly constrains the superparticle mass spectrum, which in turn constrains the neutralino annihilation mechanisms in the early universe. We find that stau and stop co-annihilation mechanisms -- already highly stressed by the latest Atlas/CMS results on SUSY searches -- are nearly eliminated if indeed the light Higgs scalar has mass m_h\simeq 125 GeV. Furthermore, neutralino annihilation via the A-resonance is essentially ruled out in mSUGRA so that it is exceedingly difficult to generate thermally-produced neutralino-only dark matter at the measured abundance. The remaining possibility lies in the focus-point region which now moves out to m_0\sim 10-20 TeV range due to the required large trilinear soft SUSY breaking term A_0. The remaining HB/FP region is more fine-tuned than before owing to the typically large top squark masses. We present updated direct and indirect detection rates for neutralino dark matter, and show that ton scale noble liquid detectors will either discover mixed higgsino CDM or essentially rule out thermally-produced neutralino-only CDM in the mSUGRA model.

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$.

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.

While the SUSY flavor, CP and gravitino problems seem to favor a very heavy spectrum of matter scalars, fine-tuning in the electroweak sector prefers low values of superpotential mass \mu. In the limit of low \mu, the two lightest neutralinos and light chargino are higgsino-like. The light charginos and neutralinos may have large production cross sections at LHC, but since they are nearly mass degenerate, there is only small energy release in three-body sparticle decays. Possible dilepton and trilepton signatures are difficult to observe after mild cuts due to the very soft p_T spectrum of the final state isolated leptons. Thus, the higgsino-world scenario can easily elude standard SUSY searches at the LHC. It should motivate experimental searches to focus on dimuon and trimuon production at the very lowest p_T(\mu) values possible. If the neutralino relic abundance is enhanced via non-standard cosmological dark matter production, then there exist excellent prospects for direct or indirect detection of higgsino-like WIMPs. While the higgsino-world scenario may easily hide from LHC SUSY searches, a linear e^+e^- collider or a muon collider operating in the \sqrt{s}\sim 0.5-1 TeV range would be able to easily access the chargino and neutralino pair production reactions.

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.

In this letter we point out that in a class of models for spontaneous R-parity breaking based on gauged B-L, the spectrum for neutrinos is quite peculiar. We find that those models generally predict three layers of neutrinos: one heavy sterile neutrino, two massive active neutrinos, and three near massless (one active and two sterile) neutrinos.

A light Higgs boson is preferred by $M_W$ and $m_t$ measurements. A complex scalar singlet addition to the Standard Model allows a better fit to these measurements through a new light singlet dominated state. It then predicts a light Dark Matter (DM) particle that can explain the signals of DM scattering from nuclei in the CoGeNT and DAMA/LIBRA experiments. Annihilations of this DM in the galactic halo, $AA\rightarrow b\bar{b}, c\bar{c}, \tau^+\tau^-$, lead to gamma rays that naturally improve a fit to the Fermi Large Area Telescope data in the central galactic regions. The associated light neutral Higgs boson may also be discovered at the Large Hadron Collider.

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.

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.

New experiments designed to discover a weakly interacting dark matter (DM) particle via spin dependent scattering can distinguish models of electroweak symmetry breaking. The plane of spin dependent versus spin independent DM scattering cross sections is a powerful model diagnostic. We detail representative predictions of mSUGRA, singlet extended SM and MSSM, a new Dirac neutrino, Littlest Higgs with T-parity (LHT) and Minimal Universal Extra Dimensions (mUED) models. Of these models, the nMSSM has the largest spin dependent (SD) cross section. It has a very light neutralino which would give lower energy nuclear recoils. The Focus Point region of mSUGRA, mUED and the right handed neutrino also predict a very large SD cross section and predict a large signal of high energy neutrinos in the IceCube experiment from annihilations of dark matter in the Sun. We also describe a model independent treatment of the scattering of DM particles of different intrinsic spins.

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.

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.

Dark matter annihilations in the Sun to neutrino-antineutrino pairs have striking signatures in neutrino detectors such as IceCube and KM3. We make a model independent study of the signals after propagation of the neutrinos from the center of the Sun to the Earth. A large spin-dependent DM capture cross section in the Sun make the discovery prospects robust.

Neutralino annihilations in the Sun to weak boson and top quark pairs lead to high-energy neutrinos that can be detected by the IceCube and KM3 experiments in the search for neutralino dark matter. We calculate the neutrino signals from real and virtual WW, ZZ, Zh, and $t \bar t$ production and decays, accounting for the spin-dependences of the matrix elements, which can have important influences on the neutrino energy spectra. We take into account neutrino propagation including neutrino oscillations, matter-resonance, absorption, and nu_tau regeneration effects in the Sun and evaluate the neutrino flux at the Earth. We concentrate on the compelling Focus Point (FP) region of the supergravity model that reproduces the observed dark matter relic density. For the FP region, the lightest neutralino has a large bino-higgsino mixture that leads to a high neutrino flux and the spin-dependent neutralino capture rate in the Sun is enhanced by 10^3 over the spin-independent rate. For the standard estimate of neutralino captures, the muon signal rates in IceCube are identifiable over the atmospheric neutrino background for neutralino masses above M_Z up to 400 GeV.

We investigate the correlated predictions of singlet extended MSSM models for direct detection of the lightest neutralino with its cosmological relic density. To illustrate the general effects of the singlet, we take heavy sleptons and squarks. We apply LEP, $(g-2)_\mu$ and perturbativity constraints. We find that the WMAP upper bound on the cold dark matter density limits much of the parameter space to regions where the lightest neutralino can be discovered in recoil experiments. The results for the NMSSM and UMSSM are typically similar to the MSSM since their light neutralinos have similar compositions and masses. In the n/sMSSM the neutralino is often very light and its recoil detection is within the reach of the CDMS II experiment. In general, most points in the parameter spaces of the singlet models we consider are accessible to the WARP experiment.

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 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.

The cosmological matter-antimatter asymmetry can originate from CP-violating interactions of seesaw Majorana neutrinos via leptogenesis in the thermal phase of the early universe. Having the cosmological CP-phase for leptogenesis requires at least two right-handed Majorana neutrinos. Using only the low energy neutrino observables we quantitatively reconstruct a minimal neutrino seesaw. We establish a general criterion for minimal seesaw schemes in which the cosmological CP-phase is {\it completely} reconstructed from the low energy CP-phases measured by neutrino oscillation and neutrinoless double-beta decay experiments. We reveal and analyze two distinct classes of such minimal schemes that are shown to be highly predictive. Extension of our reconstruction formalism to a three-heavy-neutrino seesaw is discussed.

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.

In models involving new TeV-scale Z' gauge bosons, the new U(1)' symmetry often prevents the generation of Majorana masses needed for a conventional neutrino seesaw, leading to three superweakly interacting ``right-handed'' neutrinos nu_R, the Dirac partners of the ordinary neutrinos. These can be produced prior to big bang nucleosynthesis by the Z' interactions, leading to a faster expansion rate and too much ^4He. We quantify the constraints on the Z' properties from nucleosynthesis for Z' couplings motivated by a class of E_6 models parametrized by an angle theta_E6. The rate for the annihilation of three approximately massless right-handed neutrinos into other particle pairs through the Z' channel is calculated. The decoupling temperature, which is higher than that of ordinary left-handed neutrinos due to the large Z' mass, is evaluated, and the equivalent number of new doublet neutrinos Delta N_nu is obtained numerically as a function of the Z' mass and couplings for a variety of assumptions concerning the Z-Z' mixing angle and the quark-hadron transition temperature T_c. Except near the values of theta_E6 for which the Z' decouples from the right-handed neutrinos, the Z' mass and mixing constraints from nucleosynthesis are much more stringent than the existing laboratory limits from searches for direct production or from precision electroweak data, and are comparable to the ranges that may ultimately be probed at proposed colliders. For the case T_c = 150 MeV with the theoretically favored range of Z-Z' mixings, Delta N_nu < 0.3 for M_Z' > 4.3 TeV for any value of theta_E6. Larger mixing or larger T_c often lead to unacceptably large Delta N_nu except near the nu_R decoupling limit.

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.

The combination of new cosmic microwave background data with cosmological priors has determined the physical cold dark matter (cdm) density to be $\Omega_{\rm cdm} h^2 = 0.13\pm0.01$. We find the corresponding regions of parameters in the minimal supergravity model for which the lightest neutralino is the cdm particle. We then compare with the muon anomalous magnetic moment ($a_\mu$) measurement, the mass of the lighter CP-even Higgs boson ($h^0$), the trilepton search at the upgraded Fermilab Tevatron collider, and a direct search for relic neutralinos. The intersection of $\Omega_{\rm cdm} h^2$ and $a_\mu$ constraints selects $\tan\beta\approx 40-45$.

We investigate the prospects of indirect searches for supersymmetric neutralino dark matter. Relic neutralinos gravitationally accumulate in the Sun and their annihilations produce high energy neutrinos. Muon neutrinos of this origin can be seen in large detectors like AMANDA, IceCube and ANTARES. We evaluate the relic density and the detection rate in several models---the minimal supersymmetric model, minimal supergravity, and supergravity with non-universal Higgs boson masses at the grand unification scale. We make realistic estimates for the indirect detection rates including effects of the muon detection threshold, quark hadronization, and solar absorption. We find good prospects for detection of neutralinos with mass above 200 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.

A revolution in our understanding of the neutrino sector is underway, driven by observations that are interpreted in terms of changes in neutrino flavors as they propagate. Since neutrino oscillations occur only if neutrinos are massive, these phenomena indicate physics beyond the Standard Model. With the present evidence for oscillations from atmospheric, solar and accelerator data, we are already able to begin to make strong inferences about the mass spectrum and mixings of neutrinos. Theoretical efforts to achieve a synthesis have produced a variety of models with differing testable consequences. A combination of particle physics, nuclear physics and astrophysics is needed for a full determination of the fundamental properties of neutrinos. This article reviews what has been achieved thus far and the future prospects for understanding the nature of neutrino masses and mixing.

Recent evidence for neutrino oscillations has revolutionized the study of neutrino masses and mixing. This report gives an overview of what we are learning from the neutrino oscillation experiments, the prospects for the near term, and the bright future of neutrino mass studies.

We reconsider neutrino decay as an explanation for atmospheric neutrino observations. We show that if the mass-difference relevant to the two mixed states \nu_\mu and \nu_\tau is very small (< 10^{-4} eV^2), then a very good fit to the observations can be obtained with decay of a component of \nu_\mu to a sterile neutrino and a Majoron. We discuss how the K2K and MINOS long-baseline experiments can distinguish the decay and oscillation scenarios.

In the Kaluza-Klein (KK) scenario with n large extra dimensions where gravity propagates in the 4+n dimensional bulk of spacetime while gauge and matter fields are confined to a four dimensional subspace, the light graviton KK modes can be produced in the Sun, red giants and supernovae. We study the energy-loss rates through photon-photon annihilation, electron-positron annihilation, gravi-Compton-Primakoff scattering, gravi-bremsstrahlung and nucleon-nucleon bremsstrahlung, and derive lower limits to the string scale M_S. The most stringent lower limit obtained from SN1987A leads to $M_S> 30 - 130$ TeV (2.1-9.2 TeV) for the case of two (three) large extra dimensions.

When three Majorana neutrinos describe the solar and atmospheric neutrino data via oscillations, a nonzero measurement of neutrinoless double beta ($0\nu\beta\beta$) decay can determine the sum of neutrino masses $\sum m_\nu$ if the solar solution has small-angle mixing, and place a lower bound on $\sum m_\nu$ for large-angle solar mixing. If in addition a nonzero $\sum m_\nu$ is deduced from cosmology, the neutrino mass spectrum may be uniquely specified for some ranges of neutrino parameters. For $\sum m_\nu > 0.75$ eV, the small-angle solar solution is excluded by the current upper limit on neutrinoless double beta decay. In models with maximal solar mixing the $CP$ phases of the neutrinos may be strongly constrained by stringent upper bounds on $0\nu\beta\beta$ decay.

We determine the three-neutrino mixing and mass parameters that are allowed by the solar and atmospheric neutrino data when vacuum oscillations are responsible for both phenomena. The global fit does not appreciably change the allowed regions for the parameters obtained from effective two-neutrino fits. We discuss how measurements of the solar electron energy spectrum below 6.5 GeV in Super-Kamiokande and seasonal variations in the Super-Kamiokande, $^{71}$Ga, and BOREXINO experiments can distinguish the different solar vacuum solutions.

We examine the prospects for determining the neutrino mixing matrix and for observing CP and T violation in neutrino oscillations in four-neutrino models. We focus on a general class of four-neutrino models with two pairs of nearly degenerate mass eigenstates separated by approximately 1 eV, which can describe the solar, atmospheric and LSND neutrino data. We present a general parametrization of these models and discuss in detail the determination of the mixing parameters and the mass matrix texture from current and future neutrino data in the case where $\nu_e$ and $\nu_\mu$ each mix primarily with one other neutrino. We find that measurable CP/T violation in long-baseline experiments, with amplitude at the level of the LSND signal, is possible given current experimental constraints. Also, additional oscillation effects in short- and long-baseline experiments may be measurable in many cases. We point out that, given separate scales for the mass-squared differences of the solar and atmospheric oscillations, observable CP/T violation effects in neutrino oscillations signals the existence of a sterile neutrino. We examine several textures of the neutrino mass matrix and determine which textures can have measurable CP/T violation in neutrino oscillations in long-baseline experiments. We also briefly discuss some possible origins of the neutrino mass terms in straightforward extensions of the Standard Model.

We show that the observed zenith angle dependence of the atmospheric neutrinos can be accounted for by neutrino decay. Furthermore, it is possible to account for all neutrino anomalies with just three flavors.

By combining experimental constraints from atmospheric and solar neutrino oscillations and the tritium beta decay endpoint, we infer upper limits of 4.4 eV on the $\nu_\mu$ and $\nu_\tau$ masses, if the universe consists of three neutrinos. For hierarchical mass spectra $m_3 \gg m_1, m_2$ or $m_3 \simeq m_2 \gg m_1$ we infer that $m_{\nu_\alpha} \lsim 0.08$ eV for $\alpha = e, \mu, \tau$. In a four neutrino universe, and assuming neutrino oscillations also account for the LSND experimental results, $m_{\nu_\alpha} \lsim 5.4$ eV. We also obtain lower limits on the masses.

We show that if the solar and atmospheric data are both described by maximal vacuum oscillations at the relevant mass scales then there exists a unique mixing matrix for three neutrino flavors. The solution necessarily conserves CP and automatically implies that there is no disappearance of atmospheric $\nu_e$, consistent with indications from the Super-Kamiokande experiment. We also investigate the consequences for three-neutrino mixing if the solar and atmospheric oscillations exhibit mixing that is large but not maximal. For non-maximal mixing $\nu_e\leftrightarrow\nu_\tau$ and $\nu_e\leftrightarrow\nu_\mu$ oscillations are predicted that may be observable in future long-baseline experiments.

We make a model-independent analysis of all available data that indicate neutrino oscillations. Using probability diagrams, we confirm that a mass spectrum with two nearly degenerate pairs of neutrinos separated by a mass gap of $\simeq1$ eV is preferred over a spectrum with one mass eigenstate separated from the others. We derive some new relations among the four-neutrino mixing matrix elements. We design four-neutrino mass matrices with three active neutrinos and one sterile neutrino that naturally incorporate maximal oscillations of atmospheric $\nu_\mu$ and explain the solar neutrino and LSND results. The models allow either a large or small angle MSW or vacuum oscillation description of the solar neutrino deficit. The models predict (i) oscillations of either $\nu_e \to \nu_\tau$ or $\nu_e \to\nu_s$ in long-baseline experiments at $L/E \gg 1$ km/GeV, with amplitude determined by the LSND oscillation amplitude and argument given by the atmospheric $\delta m^2$, and (ii) the equality of the $\nu_e$ disappearance probability, the $\nu_\mu$ disappearance probability, and the LSND $\nu_\mu\to\nu_e$ appearance probability in short-baseline experiments.

We present a four-neutrino model with three active neutrinos and one sterile neutrino which naturally has maximal $\nu_\mu\to\nu_\tau$ oscillations of atmospheric neutrinos and can explain the solar neutrino and LSND results. The model predicts $\nu_e \to \nu_\tau$ and $\nu_e \to \nu_\mu$ oscillations in long-baseline experiments with $L/E \gg 1$ km/GeV with amplitudes that are determined by the LSND oscillation amplitude and argument controlled by the atmospheric $\delta m^2$.

A recently reported anomaly in the time structure of signals in the KARMEN neutrino detector suggests the decay of a new particle $x$, produced in $\pi^+ \to \mu^+ x$ with mass $m_x=33.9$ MeV. We discuss the constraints and difficulties in interpreting $x$ as a neutrino. We show that a mainly-sterile neutrino scenario is compatible with all laboratory constraints, within narrow limits on the mixing parameters, although there are problems with astrophysical and cosmological constraints. This scenario predicts that appreciable numbers of other $x$-decay events with different origins and time structures should also be observable in the KARMEN detector. Such $x$-decay events should also be found in the LSND experiment and may be relevant to the search for $\bar\nu_\mu\to\bar\nu_e$ oscillations.