
We suggest an approach to perturbative calculations of largescale clustering
in the Universe that includes from the start the stream crossing (multiple
velocities for mass elements at a single position) that is lost in traditional
calculations. Starting from a functional integral over displacement, the
perturbative series expansion is in deviations from (truncated) Zel'dovich
evolution, with terms that can be computed exactly even for streamcrossed
displacements. We evaluate the oneloop formulas for displacement and density
power spectra numerically in 1D, finding dramatic improvement in agreement with
Nbody simulations compared to the Zel'dovich power spectrum (which is exact in
1D up to stream crossing). Beyond 1D, our approach could represent an
improvement over previous expansions even aside from the inclusion of stream
crossing, but we have not investigated this numerically. In the process we show
how to achieve effectivetheorylike regulation of smallscale fluctuations
without free parameters.

We derive and numerically implement an algorithm for estimating the 3D power
spectrum of the Lyman$\alpha$ (Ly$\alpha$) forest flux fluctuations. The
algorithm exploits the unique geometry of Ly$\alpha$ forest data to
efficiently measure the crossspectrum between lines of sight as a function of
parallel wavenumber, transverse separation and redshift. The key to fast
evaluation is to approximate the global covariance matrix as blockdiagonal,
where only pixels from the same spectrum are correlated. We then compute the
eigenvectors of the derivative of the signal covariance with respect to
crossspectrum parameters, and project the inversecovarianceweighted spectra
onto them. This acts much like a radial Fourier transform over redshift
windows. The resulting crossspectrum inference is then converted into our
final product, an approximation of the likelihood for the 3D power spectrum
expressed as second order Taylor expansion around a fiducial model. We
demonstrate the accuracy and scalability of the algorithm and comment on
possible extensions. Our algorithm will allow efficient analysis of the
upcoming Dark Energy Spectroscopic Instrument dataset.

We introduce FastPM, a highlyscalable approximated particle mesh Nbody
solver, which implements the particle mesh (PM) scheme enforcing correct linear
displacement (1LPT) evolution via modified kick and drift factors. Employing a
2dimensional domain decomposing scheme, FastPM scales extremely well with a
very large number of CPUs. In contrast to COmovingLAgrangian (COLA) approach,
we do not require to split the force or track separately the 2LPT solution,
reducing the code complexity and memory requirements. We compare FastPM with
different number of steps ($N_s$) and force resolution factor ($B$) against 3
benchmarks: halo mass function from Friends of Friends halo finder, halo and
dark matter power spectrum, and cross correlation coefficient (or
stochasticity), relative to a high resolution TreePM simulation. We show that
the modified time stepping scheme reduces the halo stochasticity when compared
to COLA with the same number of steps and force resolution. While increasing
$N_s$ and $B$ improves the transfer function and cross correlation coefficient,
for many applications FastPM achieves sufficient accuracy at low $N_s$ and $B$.
For example, $N_s=10$ and $B=2$ simulation provides a substantial saving (a
factor of 10) of computing time relative to $N_s=40$, $B=3$ simulation, yet the
halo benchmarks are very similar at $z=0$. We find that for abundance matched
halos the stochasticity remains low even for $N_s=5$. FastPM compares well
against less expensive schemes, being only 7 (4) times more expensive than 2LPT
initial condition generator for $N_s=10$ ($N_s=5$). Some of the applications
where FastPM can be useful are generating a large number of mocks, producing
nonlinear statistics where one varies a large number of nuisance or
cosmological parameters, or serving as part of an initial conditions solver.

Cosmic voids found in galaxy surveys are defined based on the galaxy
distribution in redshift space. We show that the large scale distribution of
voids in redshift space traces the fluctuations in the dark matter density
field \delta(k) (in Fourier space with \mu being the line of sight projected
kvector): \delta_v^s(k) = (1 + \beta_v \mu^2) b^s_v \delta(k), with a beta
factor that will be in general different than the one describing the
distribution of galaxies. Only in case voids could be assumed to be quasilocal
transformations of the linear (Gaussian) galaxy redshift space field, one gets
equal beta factors \beta_v=\beta_g=f/b_g with f being the growth rate, and b_g,
b^s_v being the galaxy and void bias on large scales defined in redshift space.
Indeed, in our mock void catalogs we measure void beta factors being in good
agreement with the galaxy one. Further work needs to be done to confirm the
level of accuracy of the beta factor equality between voids and galaxies, but
in general the void beta factor needs to be considered as a free parameter for
RSD studies.

We study the relationship between the geometry of smoothly bounded domains in
complete Riemannian manifolds and the associated sequence of $L^1$norms of
exit time moments for Brownian motion. We establish bounds for Dirichlet
eigenvalues and, for closed manifolds, we establish a comparison result for
elements of the moment sequence.

We present a Bayesian phasespace reconstruction of the cosmic largescale
matter density and velocity fields from the SDSSIII Baryon Oscillations
Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering
catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the
range of 610 $h^{1}$ Mpc, and a lognormalPoisson model with a redshift
dependent nonlinear bias. The bias parameters are derived from the data and a
general renormalised perturbation theory approach. We use combined Gibbs and
Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively
reconstruct the dark matter density field and the coherent peculiar velocities
of individual galaxies, correcting hereby for coherent redshift space
distortions (RSD). Our tests relying on accurate $N$body based mock galaxy
catalogues, show unbiased real space power spectra of the nonlinear density
field up to $k\sim0.2\, h$ Mpc$^{1}$, and vanishing quadrupoles down to
$r\sim20\,h^{1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be
obtained from the tidal field tensor based on the Gaussian component of the
reconstructed density field. We find that the reconstructed velocities have a
statistical correlation coefficient compared to the true velocities of each
individual lightcone mock galaxy of $r\sim0.68$ including about 10% of
satellite galaxies with virial motions (about $r=0.75$ without satellites). The
power spectra of the velocity divergence agree well with theoretical
predictions up to $k\sim0.2\,h\,{\rm Mpc}^{1}$. This work will be especially
useful to improve, e.g. BAO reconstructions, kinematic SunyaevZeldovich (kSZ),
integrated SachsWolfe (ISW) measurements, or environmental studies.

We study estimates involving the principal Dirichlet eigenvalue associated to
a smoothly bounded domain in a complete Riemannian manifold and L1norms of
exit time moments of Brownian motion. Our results generalize a classical
inequality of Polya.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

We study torsional rigidity for graph and quantum graph analogs of wellknown
pairs of isospectral nonisometric planar domains. We prove that such
isospectral pairs are distinguished by torsional rigidity.

We prove that heat content determines planar triangles.

We analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon
Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is
performed in Fourierspace, using the power spectrum monopole and quadrupole.
The dataset includes $1\,198\,006$ galaxies over the redshift range $0.2 < z <
0.75$. We divide this dataset into three (overlapping) redshift bins with the
effective redshifts $\zeff = 0.38$, $0.51$ and $0.61$. We demonstrate the
reliability of our analysis pipeline using Nbody simulations as well as $\sim
1000$ MultiDarkPatchy mock catalogues, which mimic the BOSSDR12 target
selection. We apply density field reconstruction to enhance the BAO
signaltonoise ratio. By including the power spectrum quadrupole we can
separate the lineofsight and angular modes, which allows us to constrain the
angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$ separately.
We obtain two independent $1.6\%$ and $1.5\%$ constraints on $D_A(z)$ and
$2.9\%$ and $2.3\%$ constraints on $H(z)$ for the low ($\zeff=0.38$) and high
($\zeff=0.61$) redshift bin, respectively. We obtain two independent $1\%$ and
$0.9\%$ constraints on the angular averaged distance $D_V(z)$, when ignoring
the AlcockPaczynski effect. The detection significance of the BAO signal is of
the order of $8\sigma$ (postreconstruction) for each of the three redshift
bins. Our results are in good agreement with the Planck prediction within
$\Lambda$CDM. This paper is part of a set that analyses the final galaxy
clustering dataset from BOSS. The measurements and likelihoods presented here
are combined with others in~\citet{Alam2016} to produce the final cosmological
constraints from BOSS.

The usual fluid equations describing the largescale evolution of mass
density in the universe can be written as local in the density, velocity
divergence, and velocity potential fields. As a result, the perturbative
expansion in small density fluctuations, usually written in terms of
convolutions in Fourier space, can be written as a series of products of these
fields evaluated at the same location in configuration space. Based on this, we
establish a new method to numerically evaluate the 1loop power spectrum (i.e.,
Fourier transform of the 2point correlation function) with onedimensional
Fast Fourier Transforms. This is exact and a few orders of magnitude faster
than previously used numerical approaches. Numerical results of the new method
are in excellent agreement with the standard quadrature integration method.
This fast model evaluation can in principle be extended to higher loop order
where existing codes become painfully slow. Our approach follows by writing
higher order corrections to the 2point correlation function as, e.g., the
correlation between two secondorder fields or the correlation between a linear
and a thirdorder field. These are then decomposed into products of
correlations of linear fields and derivatives of linear fields. The method can
also be viewed as evaluating threedimensional Fourier space convolutions using
products in configuration space, which may also be useful in other contexts
where similar integrals appear.

The source of CPTviolation in the photon sector of the Standard Model
Extension arises from a ChernSimonslike contribution that involves a coupling
to a fixed background vector field $k_{AF}^\mu$. These Lorentz and
CPTviolating photons have wellknown theoretical issues that arise from
missing states at low momenta when $k_{AF}^\mu$ is timelike. In order to make
the theory consistent, a tiny mass for the photon can be introduced, well below
current experimental bounds. The implementation of canonical quantization can
then be implemented as in the CPTpreserving case by using the St\"uckelberg
mechanism. We explicitly construct a covariant basis of properlynormalized
polarization vectors at fixed threemomentum satisfying the momentum space
field equations, in terms of which the vector field can be expanded. As an
application of the theory, we calculate the Cherenkov radiation rate for the
case of purely timelike $k_{AF}^\mu$, and find a radiation rate at high
energies that has a contribution that does not depend on the mass used to
regulate the photons.

Several lagrangians associated to classical limits of lorenzviolating
fermions in the Standard Model extension (SME) have been shown to yield Finsler
functions when the theory is expressed in Euclidean space. When spincouplings
are present, the lagrangian can develop singularities that obstruct the
construction of a globally defined Legendre transformation, leading to singular
Finsler spaces. A specific sector of the SME where such problems arise is
studied. It is found that the singular behavior can be eliminated by an
appropriate lifting of the problem to an associated algebraic variety. This
provides a smooth classical model for the singular problem. In Euclidean space,
the procedure involves combining two related singular Finsler functions into a
single smooth function with a semipositive definite quadratic form defined on
a desingularized variety.

We develop a general framework for generating estimators of a given quantity
which are unbiased to a given order in the difference between the true value of
the underlying quantity and the fiducial position in theory space around which
we expand the likelihood. We apply this formalism to rederive the optimal
quadratic estimator and show how the replacement of the second derivative
matrix with the Fisher matrix is a generic way of creating an unbiased
estimator (assuming choice of the fiducial model is independent of data). Next
apply the approach to estimation of shear lensing, closely following the work
of Bernstein and Armstrong (2014). Our first order estimator reduces to their
estimator in the limit of zero shear, but it also naturally allows for the case
of nonconstant shear and the easy calculation of correlation functions or
power spectra using standard methods. Both our firstorder estimator and
Bernstein and Armstrong's estimator exhibit a bias which is quadratic in true
shear. Our thirdorder estimator is, at least in the realm of the toy problem
of Bernstein and Armstrong, unbiased to 0.1% in relative shear errors $\Delta
g/g$ for shears up to $g=0.2$.

We present Fisher matrix projections for future cosmological parameter
measurements, including neutrino masses, dark energy, curvature, modified
gravity, the inflationary perturbation spectrum, nonGaussianity, and dark
radiation. We focus on DESI and generally redshift surveys (BOSS, HETDEX,
eBOSS, Euclid, and WFIRST), but also include CMB (Planck) and weak
gravitational lensing (DES and LSST) constraints. The goal is to present a
consistent set of projections, for concrete experiments, which are otherwise
scattered throughout many papers and proposals. We include neutrino mass as a
free parameter in most projections, as it will inevitably be relevant  DESI
and other experiments can measure the sum of neutrino masses to ~0.02 eV or
better, while the minimum possible sum is ~0.06 eV. We note that the BAOonly
use of galaxy clustering is substantially degraded as a dark energy probe in
the presence of neutrino mass uncertainty  using broadband galaxy power is
critical, especially pushing it to as small a scale as possible, and big gains
are achieved by combining lensing survey constraints with redshift survey
constraints. We do not try to be especially innovative, e.g., in careful
treatments of potential systematic errors  these projections are intended as
a straightforward baseline for comparison to more detailed analyses.

McDonald (2007) presented an approach to improving perturbation theory (PT)
calculations of the dark matter power spectrum, with a derivation based on the
idea of renormalization group flow with time. In spite of a questionable
approximation made in deriving it, subsequent comparisons by several groups
between the predictions of the resulting equation and Nbody simulations showed
remarkable improvement relative to "standard" PT (SPT) at similar order. In
this brief note I show that the same final equation can be derived cleanly from
the point of view not of flowing with time but with nonlinear coupling
strength, i.e., gradually dialing the coupling from the trivial value zero to
the physical one. This understanding makes it clear how to extend the approach
to higher order and other statistics. While I do not necessarily think this
approach is best among the many, it may be interesting in that it contains a
unique way of suppressing UV sensitivity. In passing I remind the reader of
references demonstrating that SPT works remarkably well without improvement
(except near redshift zero, where, fortunately for SPT, there is very little
volume in the Universe).

Photon quantization is implemented in the standard model extension (SME)
using the GuptaBleuler method and BRST concepts. The quantization prescription
applies to both the birefringent and nonbirefringent CPTeven couplings. A
curious incompatibility is found between the presence of the Lorentzviolating
terms and the existence of a nontrivial conjugate momentum $\Pi^0$ yielding
problems with covariant quantization procedure. Introduction of a mass
regulator term can avoid the vanishing of $\Pi^0$ and allows for the
implementation of a covariant quantization procedure. Fieldtheoretic
calculations involving the SME photons can then be performed using the mass
regulator, similar to the conventional procedure used in electrodynamics for
infrareddivergence regulation.

The Lyman$\beta$ forest refers to a region in the spectra of distant quasars
that lies between the restframe Lyman$\beta$ and Lyman$\gamma$ emissions.
The forest in this region is dominated by a combination of absorption due to
resonant Ly$\alpha$ and Ly$\beta$ scattering. When considering the 1D Ly$\beta$
forest in addition to the 1D Ly$\alpha$ forest, the full statistical
description of the data requires four 1D power spectra: Ly$\alpha$ and
Ly$\beta$ autopower spectra and the Ly$\alpha$Ly$\beta$ real and imaginary
crosspower spectra. We describe how these can be measured using an optimal
quadratic estimator that naturally disentangles Ly$\alpha$ and Ly$\beta$
contributions. Using a sample of approximately 60,000 quasar sightlines from
the BOSS Data Release 9, we make the measurement of the onedimensional power
spectrum of fluctuations due to the Ly$\beta$ resonant scattering. While we
have not corrected our measurements for resolution damping of the power and
other systematic effects carefully enough to use them for cosmological
constraints, we can robustly conclude the following: i) Ly$\beta$ power
spectrum and Ly$\alpha$Ly$\beta$ cross spectra are detected with high
statistical significance; ii) the crosscorrelation coefficient is $\approx 1$
on large scales; iii) the Ly$\beta$ measurements are contaminated by the
associated OVI absorption, which is analogous to the SiIII contamination of the
Ly$\alpha$ forest. Measurements of the Ly$\beta$ forest will allow extension of
the usable pathlength for the Ly$\alpha$ measurements while allowing a better
understanding of the physics of intergalactic medium and thus more robust
cosmological constraints.

The Dark Energy Spectroscopic Instrument (DESI) is a massively multiplexed
fiberfed spectrograph that will make the next major advance in dark energy in
the timeframe 20182022. On the Mayall telescope, DESI will obtain spectra and
redshifts for at least 18 million emissionline galaxies, 4 million luminous
red galaxies and 3 million quasistellar objects, in order to: probe the
effects of dark energy on the expansion history using baryon acoustic
oscillations (BAO), measure the gravitational growth history through
redshiftspace distortions, measure the sum of neutrino masses, and investigate
the signatures of primordial inflation. The resulting 3D galaxy maps at z<2
and Lymanalpha forest at z>2 will make 1%level measurements of the distance
scale in 35 redshift bins, thus providing unprecedented constraints on
cosmological models.

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9
(DR9) to detect and measure the position of the Baryonic Acoustic Oscillation
(BAO) feature in the threedimensional correlation function in the Lymanalpha
forest flux fluctuations at a redshift z=2.4. The feature is clearly detected
at significance between 3 and 5 sigma (depending on the broadband model and
method of error covariance matrix estimation) and is consistent with
predictions of the standard LCDM model. We assess the biases in our method,
stability of the error covariance matrix and possible systematic effects. We
fit the resulting correlation function with several models that decouple the
broadband and acoustic scale information. For an isotropic dilation factor, we
measure 100x(alpha_iso1) = 1.6 ^{+2.0+4.3+7.4}_{2.04.16.8} (stat.) +/ 1.0
(syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits)
with respect to the acoustic scale in the fiducial cosmological model (flat
LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and
transversal dilation factors we find marginalised constraints 100x(alpha_par1)
= 1.3 ^{+3.5+7.6 +12.3}_{3.36.710.2} (stat.) +/ 2.0 (syst.) and
100x(alpha_perp1) = 2.2 ^{+7.4+17}_{7.115} +/ 3.0 (syst.). The dilation
factor measurements are significantly correlated with crosscorrelation
coefficient of ~ 0.55. Errors become significantly nonGaussian for deviations
over 3 standard deviations from best fit value. Because of the data cuts and
analysis method, these measurements give tighter constraints than a previous
BAO analysis of the BOSS DR9 Lymanalpha forest sample, providing an important
consistency test of the standard cosmological model in a new redshift regime.

This work describes a new instrument optimized for a detection of the neutral
hydrogen 21cm power spectrum between redshifts of 0.51.5: the Baryon Acoustic
Oscillation Broadband and Broadbeam (BAOBAB) Array. BAOBAB will build on the
efforts of a first generation of 21cm experiments which are targeting a
detection of the signal from the Epoch of Reionization at z ~ 10. At z ~ 1, the
emission from neutral hydrogen in selfshielded overdense halos also presents
an accessible signal, since the dominant, synchrotron foreground emission is
considerably fainter than at redshift 10. The principle science driver for
these observations are Baryon Acoustic Oscillations in the matter power
spectrum which have the potential to act as a standard ruler and constrain the
nature of dark energy. BAOBAB will fully correlate dualpolarization antenna
tiles over the 600900MHz band with a frequency resolution of 300 kHz and a
system temperature of 50K. The number of antennas will grow in staged
deployments, and reconfigurations of the array will allow for both traditional
imaging and high power spectrum sensitivity operations. We present calculations
of the power spectrum sensitivity for various array sizes, with a 35element
array measuring the cosmic neutral hydrogen fraction as a function of redshift,
and a 132element system detecting the BAO features in the power spectrum,
yielding a 1.8% error on the z ~ 1 distance scale, and, in turn, significant
improvements to constraints on the dark energy equation of state over an
unprecedented range of redshifts from ~0.51.5.

We develop a perturbative approach to redshift space distortions (RSD) using
the phase space distribution function approach and apply it to the dark matter
redshift space power spectrum and its moments. RSD can be written as a sum over
density weighted velocity moments correlators, with the lowest order being
density, momentum density and stress energy density. We use standard and
extended perturbation theory (PT) to determine their auto and cross
correlators, comparing them to Nbody simulations. We show which of the terms
can be modeled well with the standard PT and which need additional terms that
include higher order corrections which cannot be modeled in PT. Most of these
additional terms are related to the small scale velocity dispersion effects,
the so called finger of god (FoG) effects, which affect some, but not all, of
the terms in this expansion, and which can be approximately modeled using a
simple physically motivated ansatz such as the halo model. We point out that
there are several velocity dispersions that enter into the detailed RSD
analysis with very different amplitudes, which can be approximately predicted
by the halo model. In contrast to previous models our approach systematically
includes all of the terms at a given order in PT and provides a physical
interpretation for the small scale dispersion values. We investigate RSD power
spectrum as a function of \mu, the cosine of the angle between the Fourier mode
and line of sight, focusing on the lowest order powers of \mu and multipole
moments which dominate the observable RSD power spectrum. Overall we find
considerable success in modeling many, but not all, of the terms in this
expansion.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the
scale of baryon acoustic oscillations (BAO) in the clustering of matter over a
larger volume than the combined efforts of all previous spectroscopic surveys
of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as
i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7.
Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000
quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5.
Early results from BOSS include the first detection of the largescale
threedimensional clustering of the Lyman alpha forest and a strong detection
from the Data Release 9 data set of the BAO in the clustering of massive
galaxies at an effective redshift z = 0.57. We project that BOSS will yield
measurements of the angular diameter distance D_A to an accuracy of 1.0% at
redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the
same redshifts. Forecasts for Lyman alpha forest constraints predict a
measurement of an overall dilation factor that scales the highly degenerate
D_A(z) and H^{1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey
is complete. Here, we provide an overview of the selection of spectroscopic
targets, planning of observations, and analysis of data and data quality of
BOSS.