
(Abridged) We describe here the most ambitious survey currently planned in
the optical, the Large Synoptic Survey Telescope (LSST). A vast array of
science will be enabled by a single widedeepfast sky survey, and LSST will
have unique survey capability in the faint time domain. The LSST design is
driven by four main science themes: probing dark energy and dark matter, taking
an inventory of the Solar System, exploring the transient optical sky, and
mapping the Milky Way. LSST will be a widefield groundbased system sited at
Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m
effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel
camera. The standard observing sequence will consist of pairs of 15second
exposures in a given field, with two such visits in each pointing in a given
night. With these repeats, the LSST system is capable of imaging about 10,000
square degrees of sky in a single filter in three nights. The typical 5$\sigma$
pointsource depth in a single visit in $r$ will be $\sim 24.5$ (AB). The
project is in the construction phase and will begin regular survey operations
by 2022. The survey area will be contained within 30,000 deg$^2$ with
$\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$,
covering the wavelength range 3201050 nm. About 90\% of the observing time
will be devoted to a deepwidefast survey mode which will uniformly observe a
18,000 deg$^2$ region about 800 times (summed over all six bands) during the
anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The
remaining 10\% of the observing time will be allocated to projects such as a
Very Deep and Fast time domain survey. The goal is to make LSST data products,
including a relational database of about 32 trillion observations of 40 billion
objects, available to the public and scientists around the world.

Intrinsic alignments (IA), correlations between the intrinsic shapes and
orientations of galaxies on the sky, are both a significant systematic in weak
lensing and a probe of the effect of largescale structure on galactic
structure and angular momentum. In the era of precision cosmology, it is thus
especially important to model IA with high accuracy. Efforts to use
cosmological perturbation theory to model the dependence of IA on the
largescale structure have thus far been relatively successful; however, extant
models do not consistently account for time evolution. In particular, advection
of galaxies due to peculiar velocities alters the impact of IA, because galaxy
positions when observed are generally different from their positions at the
epoch when IA is believed to be set. In this work, we evolve the galaxy IA from
the time of galaxy formation to the time at which they are observed, including
the effects of this advection, and show how this process naturally leads to a
dependence of IA on the velocity shear. We calculate the galaxygalaxyIA
bispectrum to tree level (in the linear matter density) in terms of the evolved
IA coefficients. We then discuss the implications for weak lensing systematics
as well as for studies of galaxy formation and evolution. We find that
considering advection introduces nonlocality into the bispectrum, and that the
degree of nonlocality represents the memory of a galaxy's path from the time of
its formation to the time of observation. We discuss how this result can be
used to constrain the redshift at which IA is determined and provide Fisher
estimation for the relevant measurements using the example of SDSSBOSS.

Previous studies have shown the filamentary structures in the cosmic web
influence the alignments of nearby galaxies. We study this effect in the LOWZ
sample of the Sloan Digital Sky Survey using the "Cosmic Web Reconstruction"
filament catalogue of Chen et al. (2016). We find that LOWZ galaxies exhibit a
small but statistically significant alignment in the direction parallel to the
orientation of nearby filaments. This effect is detectable even in the absence
of nearby galaxy clusters, which suggests it is an effect from the matter
distribution in the filament. A nonparametric regression model suggests that
the alignment effect with filaments extends over separations of $3040$ Mpc. We
find that galaxies that are bright and earlyforming align more strongly with
the directions of nearby filaments than those that are faint and lateforming;
however, trends with stellar mass are less statistically significant, within
the narrow range of stellar mass of this sample.

Galaxy intrinsic alignments (IA) are a critical uncertainty for current and
future weak lensing measurements. We describe a perturbative expansion of IA,
analogous to the treatment of galaxy biasing. From an astrophysical
perspective, this model includes the expected largescale alignment mechanisms
for galaxies that are pressuresupported (tidal alignment) and
rotationsupported (tidal torquing) as well as the crosscorrelation between
the two. Alternatively, this expansion can be viewed as an effective model
capturing all relevant effects up to the given order. We include terms up to
second order in the density and tidal fields and calculate the resulting IA
contributions to twopoint statistics at oneloop order. For fiducial
amplitudes of the IA parameters, we find the quadratic alignment and
linearquadratic cross terms can contribute orderunity corrections to the
total intrinsic alignment signal at $k\sim0.1\,h^{1}{\rm Mpc}$, depending on
the source redshift distribution. These contributions can lead to significant
biases on inferred cosmological parameters in Stage IV photometric weak lensing
surveys. We perform forecasts for an LSSTlike survey, finding that use of the
standard "NLA" model for intrinsic alignments cannot remove these large
parameter biases, even when allowing for a more general redshift dependence.
The model presented here will allow for more accurate and flexible IA treatment
in weak lensing and combined probes analyses, and an implementation is made
available as part of the public FASTPT code. The model also provides a more
advanced framework for understanding the underlying IA processes and their
relationship to fundamental physics.

We examine the cosmological information available from the 1point
probability distribution (PDF) of the weaklensing convergence field, utilizing
fast LPICOLA simulations and a Fisher analysis. We find competitive
constraints in the $\Omega_m$$\sigma_8$ plane from the convergence PDF with
$188\ arcmin^2$ pixels compared to the cosmic shear power spectrum with an
equivalent number of modes ($\ell < 886$). The convergence PDF also partially
breaks the degeneracy cosmic shear exhibits in that parameter space. A joint
analysis of the convergence PDF and shear 2point function also reduces the
impact of shape measurement systematics, to which the PDF is less susceptible,
and improves the total figure of merit by a factor of $23$, depending on the
level of systematics. Finally, we present a correction factor necessary for
calculating the unbiased Fisher information from finite differences using a
limited number of cosmological simulations.

We present cosmological constraints from the Dark Energy Survey (DES) using a
combined analysis of angular clustering of red galaxies and their
crosscorrelation with weak gravitational lensing of background galaxies. We
use a 139 square degree contiguous patch of DES data from the Science
Verification (SV) period of observations. Using large scale measurements, we
constrain the matter density of the Universe as Omega_m = 0.31 +/ 0.09 and the
clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/ 0.13
after marginalizing over seven nuisance parameters and three additional
cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16}
= 0.74 +/ 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 =
0.78 +/ 0.09 using two bins over the range 0.2 <z< 0.5. We study the
robustness of the results under changes in the data vectors, modelling and
systematics treatment, including photometric redshift and shear calibration
uncertainties, and find consistency in the derived cosmological parameters. We
show that our results are consistent with previous cosmological analyses from
DES and other data sets and conclude with a joint analysis of DES angular
clustering and galaxygalaxy lensing with Planck CMB data, Baryon Accoustic
Oscillations and Supernova type Ia measurements.

At the epoch of decoupling, cosmic baryons had supersonic velocities relative
to the dark matter that were coherent on large scales. These velocities
subsequently slow the growth of smallscale structure and, via feedback
processes, can influence the formation of larger galaxies. We examine the
effect of streaming velocities on the galaxy correlation function, including
all leadingorder contributions for the first time. We find that the impact on
the BAO peak is dramatically enhanced (by a factor of ~5) over the results of
previous investigations, with the primary new effect due to advection: if a
galaxy retains memory of the primordial streaming velocity, it does so at its
Lagrangian, rather than Eulerian, position. Since correlations in the streaming
velocity change rapidly at the BAO scale, this advection term can cause a
significant shift in the observed BAO position. If streaming velocities impact
tracer density at the 1% level, compared to the linear bias, the recovered BAO
scale is shifted by approximately 0.5%. This new effect, which is required to
preserve Galilean invariance, greatly increases the importance of including
streaming velocities in the analysis of upcoming BAO measurements and opens a
new window to the astrophysics of galaxy formation.

Intrinsic alignment (IA) of source galaxies is one of the major astrophysical
systematics for ongoing and future weak lensing surveys. This paper presents
the first forecasts of the impact of IA on cosmic shear measurements for
current and future surveys (DES, Euclid, LSST, WFIRST) using simulated
likelihood analyses and realistic covariances that include higherorder moments
of the density field in the computation. We consider a range of possible IA
scenarios and test mitigation schemes, which parameterize IA by the fraction of
red galaxies, normalization, luminosity and redshift dependence of the IA
signal (for a subset we consider joint IA and photoz uncertainties). Compared
to previous studies we find smaller biases in timedependent dark energy models
if IA is ignored in the analysis; the amplitude and significance of these
biases vary as a function of survey properties (depth, statistical
uncertainties), luminosity function, and IA scenario: Due to its small
statistical errors and relatively shallow observing strategy Euclid is
significantly impacted by IA. LSST and WFIRST benefit from their increased
survey depth, while the larger statistical errors for DES decrease IA's
relative impact on cosmological parameters. The proposed IA mitigation scheme
removes parameter biases due to IA for DES, LSST, and WFIRST even if the shape
of the IA power spectrum is only poorly known; successful IA mitigation for
Euclid requires more prior information. We explore several alternative IA
mitigation strategies for Euclid; in the absence of alignment of blue galaxies
we recommend the exclusion of red (IA contaminated) galaxies in cosmic shear
analyses. We find that even a reduction of 20% in the number density of
galaxies only leads to a 410% loss in cosmological constraining power.

Coherent alignments of galaxy shapes, often called "intrinsic alignments"
(IA), are the most significant source of astrophysical uncertainty in weak
lensing measurements. We develop the tidal alignment model of IA and
demonstrate its success in describing observational data. We also describe a
technique to separate IA from galaxygalaxy lensing measurements. Applying this
technique to luminous red galaxy lenses in the Sloan Digital Sky Survey, we
constrain potential IA contamination from associated sources to be below a few
percent.

We develop an analytic model for galaxy intrinsic alignments (IA) based on
the theory of tidal alignment. We calculate all relevant nonlinear corrections
at oneloop order, including effects from nonlinear density evolution, galaxy
biasing, and source density weighting. Contributions from density weighting are
found to be particularly important and lead to bias dependence of the IA
amplitude, even on large scales. This effect may be responsible for much of the
luminosity dependence in IA observations. The increase in IA amplitude for more
highly biased galaxies reflects their locations in regions with large tidal
fields. We also consider the impact of smoothing the tidal field on halo
scales. We compare the performance of this consistent nonlinear model in
describing the observed alignment of luminous red galaxies with the linear
model as well as the frequently used "nonlinear alignment model," finding a
significant improvement on small and intermediate scales. We also show that the
crosscorrelation between density and IA (the "GI" term) can be effectively
separated into source alignment and source clustering, and we accurately model
the observed alignment down to the onehalo regime using the tidal field from
the fully nonlinear halomatter cross correlation. Inside the onehalo regime,
the average alignment of galaxies with density tracers no longer follows the
tidal alignment prediction, likely reflecting nonlinear processes that must be
considered when modeling IA on these scales. Finally, we discuss tidal
alignment in the context of cosmic shear measurements.

We examine the signature of dynamic (redshiftspace) distortions and
geometric distortions (including the AlcockPaczynski effect) in the context of
the galaxy power spectrum measured in upcoming galaxy redshift surveys.
Information comes from both the baryon acoustic oscillation (BAO) feature and
the broadband power spectrum shape. Accurate modeling is required to extract
this information without introducing systematic biases in the result. We
consider an analytic model for the power spectrum of dark matter halos in
redshift space, based on the distribution function expansion, and compare with
halo clustering measured in Nbody simulations. We forecast that the
distribution function model is sufficiently accurate to allow the inclusion of
broadband information on scales down to k~0.2 h/Mpc, with somewhat better
accuracy for higher bias halos. Compared with a BAOonly analysis with
reconstruction, including broadband shape information can improve unbiased
constraints on distance measures H(z) and D_A(z) by ~30% and 20%, respectively,
for a galaxy sample similar to the DESI luminous red galaxies. The gains in
precision are larger in the absence of BAO reconstruction. Furthermore,
including broadband shape information allows the measurement of structure
growth, through redshiftspace distortions. For the same galaxy sample, the
distribution function model is able to constrain f*sigma_8 to ~2%, when
simultaneously fitting for H(z) and D_A(z). We discuss techniques to optimize
the analysis of the power spectrum, including removing modes near the
lineofsight that are particularly challenging to model, and whether these
approaches can improve parameter constraints. We find that such techniques are
unlikely to significantly improve constraints on geometry, although they may
allow higher precision measurements of redshiftspace distortions.

The coherent physical alignment of galaxies is an important systematic for
gravitational lensing studies as well as a probe of the physical mechanisms
involved in galaxy formation and evolution. We develop a formalism for treating
this intrinsic alignment (IA) in the context of galaxygalaxy lensing and
present an improved method for measuring IA contamination, which can arise when
sources physically associated with the lens are placed behind the lens due to
photometric redshift scatter. We apply the technique to recent Sloan Digital
Sky Survey (SDSS) measurements of Luminous Red Galaxy lenses and typical (L*)
source galaxies with photometric redshifts selected from the SDSS imaging data.
Compared to previous measurements, this method has the advantage of being fully
selfconsistent in its treatment of the IA and lensing signals, solving for the
two simultaneously. We find an IA signal consistent with zero, placing tight
constraints on both the magnitude of the IA effect and its potential
contamination to the lensing signal. While these constraints depend on source
selection and redshift quality, the method can be applied to any measurement
that uses photometric redshifts. We obtain a modelindependent upperlimit of
roughly 10% IA contamination for projected separations of approximately 0.1100
Mpc/h. With more stringent photoz cuts and reasonable assumptions about the
physics of intrinsic alignments, this upper limit is reduced to 12%. These
limits are well below the statistical error of the current lensing
measurements. Our results suggest that IA will not present intractable
challenges to the next generation of galaxygalaxy lensing experiments, and the
methods presented here should continue to aid in our understanding of alignment
processes and in the removal of IA from the lensing signal.

Weak gravitational lensing has become a powerful probe of largescale
structure and cosmological parameters. Precision weak lensing measurements
require an understanding of the intrinsic alignment of galaxy ellipticities,
which can in turn inform models of galaxy formation. It is hypothesized that
elliptical galaxies align with the background tidal field and that this
alignment mechanism dominates the correlation between ellipticities on
cosmological scales (in the absence of lensing). We use recent largescale
structure measurements from the Sloan Digital Sky Survey to test this picture
with several statistics: (1) the correlation between ellipticity and galaxy
overdensity, w_{g+}; (2) the intrinsic alignment autocorrelation functions;
(3) the correlation functions of curlfree, E, and divergencefree, B, modes
(the latter of which is zero in the linear tidal alignment theory); (4) the
alignment correlation function, w_g(r_p,theta), a recently developed statistic
that generalizes the galaxy correlation function to account for the angle
between the galaxy separation vector and the principle axis of ellipticity. We
show that recent measurements are largely consistent with the tidal alignment
model and discuss dependence on galaxy luminosity. In addition, we show that at
linear order the tidal alignment model predicts that the angular dependence of
w_g(r_p,theta) is simply w_{g+}*cos(2*theta) and that this dependence is
consistent with recent measurements. We also study how stochastic nonlinear
contributions to galaxy ellipticity impact these statistics. We find that a
significant fraction of the observed LRG ellipticity can be explained by
alignment with the tidal field on scales >~10 h^1 Mpc. These considerations
are relevant to galaxy formation and evolution.