
We perform the covariant canonical quantization of the CPT and
Lorentzsymmetryviolating photon sector of the minimal StandardModel
Extension, which contains a general (timelike, lightlike, or spacelike) fixed
background tensor $k_{AF}^\mu$. Wellknown stability issues, arising from
complexvalued energy states, are solved by introducing a small photon mass,
orders of magnitude below current experimental bounds. We explicitly construct
a covariant basis of polarization vectors, in which the photon field can be
expanded. We proceed to derive the Feynman propagator and show that the theory
is microcausal. Despite the occurrence of negative energies and
vacuumCherenkov radiation, we do not find any runaway stability issues,
because the energy remains bounded from below. An important observation is that
the ordering of the roots of the dispersion relations is the same in any
observer frame, which allows for a frameindependent condition that selects the
correct branch of the dispersion relation. This turns out to be critical for
the consistency of the quantization. To our knowledge, this is the first system
for which quantization has consistently been performed, in spite of the fact
that the theory contains negative energies in some observer frames.

In this followup work to the High Energy Physics Community Summer Study 2013
(HEP CSS 2013, a.k.a. Snowmass), we explore the scientific capabilities of a
future StageIV Cosmic Microwave Background polarization experiment (CMBS4)
under various assumptions on detector count, resolution, and sky coverage. We
use the Fisher matrix technique to calculate the expected uncertainties in
cosmological parameters in $\nu \Lambda$CDM that are especially relevant to the
physics of fundamental interactions, including neutrino masses, effective
number of relativistic species, darkenergy equation of state, darkmatter
annihilation, and inflationary parameters. To further chart the landscape of
future cosmology probes, we include forecasted results from the Baryon Acoustic
Oscillation (BAO) signal as measured by DESI to constrain parameters that would
benefit from low redshift information. We find the following best 1sigma
constraints: $\sigma$ constraints: $\sigma(M_{\nu})= 15$ meV, $\sigma(N_{\rm
eff})= 0.0156$, Dark energy Figure of Merit = 303, $\sigma(p_{ann})=
0.00588\times3\times10^{26}$ cm$^3$/s/GeV, $\sigma(\Omega_K)= 0.00074$,
$\sigma(n_s)= 0.00110$, $\sigma(\alpha_s)= 0.00145$, and $\sigma(r)= 0.00009$.
We also detail the dependences of the parameter constraints on detector count,
resolution, and sky coverage.

This white paper is a condensation of a report by a committee appointed
jointly by the Nuclear Science and Physics Divisions at Lawrence Berkeley
National Laboratory (LBNL). The goal of this study was to identify the most
promising technique(s) for resolving the neutrino mass hierarchy. For the most
part, we have relied on calculations and simulations presented by the
proponents of the various experiments. We have included evaluations of the
opportunities and challenges for these experiments based on what is available
already in the literature.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST,
the Wide Field InfraRed Survey Telescope, recommended by the Astro 2010
Decadal Survey as the highest priority for a large space mission. The SDT was
chartered to work with the WFIRST Project Office at GSFC and the Program Office
at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the
original charge was to produce an interim design reference mission by mid2011.
That document was delivered to NASA and widely circulated within the
astronomical community. In late 2011 the Astrophysics Division augmented its
original charge, asking for two design reference missions. The first of these,
DRM1, was to be a finalized version of the interim DRM, reducing overall
mission costs where possible. The second of these, DRM2, was to identify and
eliminate capabilities that overlapped with those of NASA's James Webb Space
Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's groundbased
Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall
mission cost, while staying faithful to NWNH. This report presents both DRM1
and DRM2.

BigBOSS is a Stage IV groundbased dark energy experiment to study baryon
acoustic oscillations (BAO) and the growth of structure with a widearea galaxy
and quasar redshift survey over 14,000 square degrees. It has been
conditionally accepted by NOAO in response to a call for major new
instrumentation and a highimpact science program for the 4m Mayall telescope
at Kitt Peak. The BigBOSS instrument is a roboticallyactuated, fiberfed
spectrograph capable of taking 5000 simultaneous spectra over a wavelength
range from 340 nm to 1060 nm, with a resolution R = 30004800.
Using data from imaging surveys that are already underway, spectroscopic
targets are selected that trace the underlying dark matter distribution. In
particular, targets include luminous red galaxies (LRGs) up to z = 1.0,
extending the BOSS LRG survey in both redshift and survey area. To probe the
universe out to even higher redshift, BigBOSS will target bright [OII] emission
line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are
obtained to measure the BAO feature, trace the matter power spectrum at smaller
scales, and detect redshift space distortions. BigBOSS will provide additional
constraints on early dark energy and on the curvature of the universe by
measuring the Lyalpha forest in the spectra of over 600,000 2.2 < z < 3.5
quasars.
BigBOSS galaxy BAO measurements combined with an analysis of the broadband
power, including the Lyalpha forest in BigBOSS quasar spectra, achieves a FOM
of 395 with Planck plus Stage III priors. This FOM is based on conservative
assumptions for the analysis of broad band power (kmax = 0.15), and could grow
to over 600 if current work allows us to push the analysis to higher wave
numbers (kmax = 0.3). BigBOSS will also place constraints on theories of
modified gravity and inflation, and will measure the sum of neutrino masses to
0.024 eV accuracy.

We analyze the SDSS Lyalpha forest P_F(k,z) measurement to determine the
linear theory power spectrum. Our analysis is based on fully hydrodynamic
simulations, extended using hydroPM simulations. We account for the effect of
absorbers with damping wings, which leads to an increase in the slope of the
linear power spectrum. We break the degeneracy between the mean level of
absorption and the linear power spectrum without significant use of external
constraints. We infer linear theory power spectrum amplitude
Delta^2_L(k_p=0.009s/km,z_p=3.0)=0.452_{0.0570.116}^{+0.069+0.141} and slope
n_eff=2.321_{0.0470.102}^{+0.055+0.131} (possible systematic errors are
included through nuisance parameters in the fit  a factor >~5 smaller errors
would be obtained on both parameters if we ignored modeling uncertainties). The
errors are correlated and not perfectly Gaussian, so we provide a chi^2 table
to accurately describe the results. The result corresponds to sigma_8=0.85,
n=0.94, for a LCDM model with Omega_m=0.3, Omega_b=0.04, and h=0.7, but is most
useful in a combined fit with the CMB. The inferred curvature of the linear
power spectrum and the evolution of its amplitude and slope with redshift are
consistent with expectations for LCDM models, with the evolution of the slope,
in particular, being tightly constrained. We use this information to constrain
systematic contamination, e.g., fluctuations in the UV background. This paper
should serve as a starting point for more work to refine the analysis,
including technical improvements such as increasing the size and number of the
hydrodynamic simulations, and improvements in the treatment of the various
forms of feedback from galaxies and quasars.

We combine the measurements of luminosity dependence of bias with the
luminosity dependent weak lensing analysis of dark matter around galaxies to
derive the galaxy bias and constrain nonlinear mass and cosmological
parameters. We take advantage of theoretical and simulation predictions that
predict that while halo bias is rapidly increasing with mass for high masses,
it is nearly constant in low mass halos. We use a new weak lensing analysis
around the same SDSS galaxies to determine their halo mass probability
distribution. These halo mass probability distributions are used to predict the
bias for each luminosity subsample and we find an excellent agreement with
observed values. The required nonlinear mass suggests slightly lower matter
density than usually assumed, Omegam=0.25+/ 0.03 for the simplest models. We
combine the bias constraints with those from the WMAP and the SDSS power
spectrum analysis to derive new constraints on bias and sigma_8. For the most
general parameter space we find sigma_8=0.88+/ 0.06 and b_*=0.99+/ 0.07. In
the context of spatially flat models we improve the limit on the neutrino mass
for the case of 3 degenerate families from m_nu<0.6eV without bias to
m_nu<0.18eV with bias (95% c.l.), which is weakened to m_nu<0.24eV if running
is allowed. The corresponding limit for 3 massless + 1 massive neutrino is
1.37eV.

We combine the constraints from the recent Lyalpha forest and bias analysis
of the SDSS with previous constraints from SDSS galaxy clustering, the latest
supernovae, and WMAP . Combining WMAP and the lya we find for the primordial
slope n_s=0.98\pm 0.02. We see no evidence of running, dn/d\ln k=0.003\pm
0.010, a factor of 3 improvement over previous constraints. We also find no
evidence of tensors, r<0.36 (95% c.l.). A positive correlation between tensors
and primordial slope disfavors chaotic inflation type models with steep slopes:
V \propto \phi^4 is outside the 3sigma contour. For the amplitude we find
sigma_8=0.90\pm 0.03 from the lyaf and WMAP alone. We find no evidence of
neutrino mass: for the case of 3 massive neutrino families with an inflationary
prior, \sum m_{\nu}<0.42eV and the mass of lightest neutrino is m_1<0.13eV at
95% c.l. For the 3 massless + 1 massive neutrino case we find m_{\nu}<0.79eV
for the massive neutrino, excluding at 95% c.l. all neutrino mass solutions
compatible with the LSND results. We explore dark energy constraints in models
with a fairly general time dependence of dark energy equation of state, finding
Omega_lambda=0.72\pm 0.02, w(z=0.3)=0.98^{+0.10}_{0.12}, the latter changing
to w(z=0.3)=0.92^{+0.09}_{0.10} if tensors are allowed. We find no evidence
for variation of the equation of state with redshift,
w(z=1)=1.03^{+0.21}_{0.28}. These results rely on the current understanding
of the lyaf and other probes, which need to be explored further both
observationally and theoretically, but extensive tests reveal no evidence of
inconsistency among different data sets used here.

We explore several physical effects on the power spectrum of the Lymanalpha
forest transmitted flux. The effects we investigate here are usually not part
of hydrodynamic simulations and so need to be estimated separately. The most
important effect is that of high column density absorbers with damping wings,
which add power on large scales. We compute their effect using the
observational constraints on their abundance as a function of column density.
Ignoring their effect leads to an underestimation of the slope of the linear
theory power spectrum. The second effect we investigate is that of fluctuations
in the ionizing radiation field. For this purpose we use a very large high
resolution Nbody simulation, which allows us to simulate both the fluctuations
in the ionizing radiation and the small scale LyaF within the same simulation.
We find an enhancement of power on large scales for quasars and a suppression
for galaxies. The strength of the effect rapidly increases with increasing
redshift, allowing it to be uniquely identified in cases where it is
significant. We develop templates which can be used to search for this effect
as a function of quasar lifetime, quasar luminosity function, and attenuation
length. Finally, we explore the effects of galactic winds using hydrodynamic
simulations. We find the wind effects on the LyaF power spectrum to be be
degenerate with parameters related to the temperature of the gas that are
already marginalized over in cosmological fits. While more work is needed to
conclusively exclude all possible systematic errors, our results suggest that,
in the context of data analysis procedures where parameters of the LyaF model
are properly marginalized over, the flux power spectrum is a reliable tracer of
cosmological information.

We investigate the promise of the Lyalpha forest for high precision
cosmology in the era of the Sloan Digital Sky Survey using low order Npoint
statistics. We show that with the existing data one can determine the
amplitude, slope and curvature of the slope of the matter power spectrum with a
few percent precision. Higher order statistics such as the bispectrum provide
independent information that can confirm and improve upon the statistical
precision from the power spectrum alone. The achievable precision is comparable
to that from the cosmic microwave background with upcoming satellites, and
complements it by measuring the power spectrum amplitude and shape at smaller
scales. Since the data cover the redshift range 2<z<4, one can also extract the
evolution of the growth factor and Hubble parameter over this range, and
provide useful constraints on the presence of dark energy at z>2.

We consider the possibility that massive primordial black holes are the
dominant form of dark matter. Black hole formation generates entropy
fluctuations that adds a Poisson noise to the matter power spectrum. We use
Lymanalpha forest observations to constrain this Poisson term in matter power
spectrum, then we constrain the mass of black holes to be less than few times
10^4 solar mass. We also find that structures with less than ~ 10^3 primordial
black holes evaporate by now.