
ProtoDUNESP is the singlephase DUNE Far Detector prototype that is under
construction and will be operated at the CERN Neutrino Platform (NP) starting
in 2018. ProtoDUNESP, a crucial part of the DUNE effort towards the
construction of the first DUNE 10kt fiducial mass far detector module (17 kt
total LAr mass), is a significant experiment in its own right. With a total
liquid argon (LAr) mass of 0.77 kt, it represents the largest monolithic
singlephase LArTPC detector to be built to date. It's technical design is
given in this report.

We derive expressions for the neutrino mixing parameters that result from
complex perturbations on (1) the Majorana neutrino mass matrix (in the basis of
charged lepton mass eigenstates) and on (2) the charged lepton mass matrix, for
arbitrary initial (unperturbed) mixing matrices. In the first case, we find
that the phases of the elements of the perturbation matrix, and the initial
values of the Dirac and Majorana phases, strongly impact the leading order
corrections to the neutrino mixing parameters and phases. For experimentally
compatible scenarios wherein the initial neutrino mass matrix has $\mu\tau$
symmetry, we find that the Dirac phase can take any value under small
perturbations. Similarly, in the second case, perturbations to the charged
lepton mass matrix can generate large corrections to the mixing angles and
phases of the PMNS matrix. As an illustration of our generalized procedure, we
apply it to a situation in which nonstandard scalar and nonstandard vector
interactions simultaneously affect neutrino oscillations.

In the context of the Type I seesaw mechanism, we carry out a systematic
study of the constraints that result from zeros in both the Dirac and
righthanded Majorana neutrino mass matrices. We find that most constraints can
be expressed in the standard form with one or two element/cofactor zeros alone,
while there are 9 classes of nonstandard constraints. We show that all the
constraints are stable under oneloop renormalization group running from the
lightest righthanded neutrino mass scale to the electroweak scale. We study
the predictions of the nonstandard constraints for the lightest neutrino mass,
Dirac CP phase and neutrinoless double beta decay.

We study a model with partial quarklepton universality that can naturally
arise in grand unified theories. We find that constraints on the model can be
reduced to a single condition on the Dirac CP phase $\delta$ in the neutrino
sector. Using our current knowledge of the CKM and PMNS mixing matrices, we
predict $32.4^\circ < \delta < 32.0^\circ$ at $2\sigma$.

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
neutrinonucleon 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
neutrinonucleon cross section, (ii) a new elastic neutral current interaction,
and (iii) a new completely inelastic interaction, individually impact event
rates.

We constrain the possibility of a new pseudoscalar coupling between the muon
and proton using a recent measurement of the 2S hyperfine splitting in muonic
hydrogen.

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

There may be a highenergy cutoff of neutrino events in IceCube data. In
particular, IceCube does not observe either continuum events above 2 PeV, or
the Standard Model Glashowresonance events expected at 6.3 PeV. There are also
no higher energy neutrino signatures in the ANITA and Auger experiments. This
absence of highenergy 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 Lorentzinvariance 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 twobody 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 study Majorana neutrino mass matrices that have two texture zeros, or two
cofactor zeros, or one texture zero and one cofactor zero. The two
texture/cofactor zero conditions give four constraints, which in conjunction
with the five measured oscillation parameters completely determine the nine
independent real parameters of the neutrino mass matrix. We also study the
implications that future measurements of neutrinoless double beta decay and the
Dirac CP phase will have on these cases.

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.

With the recent observation of nonzero \theta_{13}, five neutrino oscillation
parameters are now known. By imposing four zeros in the Yukawa coupling matrix
of the type I seesaw model, the number of parameters in the neutrino mass
matrix is reduced to seven, and we are able to make predictions for the
lightest neutrino mass, Dirac CP phase, and neutrinoless double beta decay.
Four texture zeros in the Yukawa coupling matrix is equivalent to either a
single texture zero or a single cofactor zero for an offdiagonal element of
the light neutrino mass matrix. We find strong similarities between single
texture zero models with one mass ordering and single cofactor zero models with
the opposite mass ordering. In the context of a specific class of singleflavor
leptogenesis models, we find additional constraints on the parameter space.

Many neutrino mixing scenarios that have \mu\tau symmetry with \theta_{13}=0
are in disagreement with recent experimental results that indicate a nonzero
value for \theta_{13}. We investigate the effect of small perturbations on
Majorana mass matrices with \mu\tau symmetry and derive analytic formulae for
the corrections to the mixing angles. We find that since m_1 and m_2 are nearly
degenerate, \mu\tau symmetry mixing scenarios are able to explain the
experimental data with about the same size perturbation for most values of
\theta_{12}. This suggests that the underlying unperturbed mixing need not have
\theta_{12} close to the experimentally preferred value. One consequence of
this is that a new class of models with \mu\tau symmetry is possible, with
unperturbed \theta_{12} equal to zero or 90 degrees for arbitrary unperturbed
\theta_{13}.

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

Motivated by the tentative observation of superluminal neutrinos by the OPERA
experiment, we present a model of activesterile neutrino oscillations in which
sterile neutrinos are superluminal and active neutrinos appear superluminal by
virtue of neutrino mixing. The model demonstrates some interesting
possibilities and challenges that apply to a large class of models aiming to
explain the OPERA result.

Data from the LSND and MiniBooNE experiments, and revised expectations of the
antineutrino flux from nuclear reactors suggest the existence of eVmass
sterile neutrinos. 3+2 and 1+3+1 scenarios accommodate all relevant
shortbaseline neutrino data except for the lowenergy MiniBooNE anomaly. We
analyze the angular distribution of upward going atmospheric neutrino events in
the IceCube40 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.

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 nittygritty 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
helicitydependence 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 zenithangle 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.

We describe a characteristic signature of dark matter (DM) annihilation or
decay into gammarays. 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 3body 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 selfconjugate 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 gammaray data that portends a striking signal at atmospheric
Cherenkov telescopes (ACTs).

We analyze new diffuse gammaray data from the Fermi Gammaray Space
Telescope, which do not confirm an excess in the EGRET data at galactic
midlatitudes, 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 gammaray signal of DM annihilation
from the galactic center while seeing no excess at midlatitudes.

Assuming that the positron excess in PAMELA satellite data is a consequence
of annihilations of cold dark matter, we consider from a modelindependent
perspective if the data show a preference for the spin of dark matter. We then
perform a general analysis of annihilations into twobody 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 AstroParticle Physics
(CCAPP) at The Ohio State University (OSU) organized a 3day 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 NNbar 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 cosigning the white paper.

We show that the rate for upward showers from an isotropic cosmic neutrino
flux at neutrino telescopes like IceCube is independent of the neutrinonucleon
cross section. For bins that span a relatively narrow range in energy, neither
scaling the cross section, nor changing its powerlaw 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
energydependence of the cross section, independent of the powerlaw 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 massvarying 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 discuss the physics of superbeam upgrades, where we focus on T2KK, a NuMI
beam line based experiment NOvA*, and a wide band beam (WBB) experiment
independent of the NuMI beam line. For T2KK, we find that the JapanKorea
baseline helps resolve parameter degeneracies, but the improvement due to
correlated systematics between the two detectors (using identical detectors) is
only moderate. For an upgrade of NOvA with a liquid argon detector, we
demonstrate that the Ash River site is preferred compared to alternatives, such
as at the second oscillation maximum, and is the optimal site within the U.S.
For a WBB experiment, we find that high proton energies and long decay tunnels
are preferable. We compare water Cherenkov and liquid argon technologies, and
find the breakeven point in detector cost at about 4:1. In order to compare
the physics potential of the different experimental configurations, we use the
concept of exposure to normalize the performance. We find that experiments with
WBBs are the best experimental concept. NOvA* could be competitive with
sufficient luminosity. If $\sin^2 2\theta_{13}$ > 0.01, a WBB experiment can
perform better than a neutrino factory.

We show that the combined data from solar, longbaseline and reactor neutrino
experiments can exclude the generalized bicycle model of Lorentz noninvariant
directiondependent and/or directionindependent oscillations of massless
neutrinos. This model has five parameters, which is more than is needed in
standard oscillation phenomenology with neutrino masses. Solar data alone are
sufficient to exclude the pure directiondependent case. The combination of
solar and longbaseline data rules out the pure directionindependent case.
With the addition of KamLAND data, a mixture of directiondependent and
directionindependent terms in the effective Hamiltonian is also excluded.

This report provides the results of an extensive and important study of the
potential for a U.S. scientific program that will extend our knowledge of
neutrino oscillations well beyond what can be anticipated from ongoing and
planned experiments worldwide. The program examined here has the potential to
provide the U.S. particle physics community with world leading experimental
capability in this intensely interesting and active field of fundamental
research. Furthermore, this capability could be unique compared to anywhere
else in the world because of the available beam intensity and baseline
distances. The present study was initially commissioned in April 2006 by top
research officers of Brookhaven National Laboratory and Fermi National
Accelerator Laboratory and, as the study evolved, it also provided responses to
questions formulated and addressed to the study group by the Neutrino
Scientific Advisory Committee (NuSAG) of the U.S. DOE and NSF. The participants
in the study, its Charge and history, plus the study results and conclusions
are provided in this report and its appendices. A summary of the conclusions is
provided in the Executive Summary.