
The Large Synoptic Survey Telescope is designed to provide an unprecedented
optical imaging dataset that will support investigations of our Solar System,
Galaxy and Universe, across half the sky and over ten years of repeated
observation. However, exactly how the LSST observations will be taken (the
observing strategy or "cadence") is not yet finalized. In this
dynamicallyevolving community white paper, we explore how the detailed
performance of the anticipated science investigations is expected to depend on
small changes to the LSST observing strategy. Using realistic simulations of
the LSST schedule and observation properties, we design and compute diagnostic
metrics and Figures of Merit that provide quantitative evaluations of different
observing strategies, analyzing their impact on a wide range of proposed
science projects. This is work in progress: we are using this white paper to
communicate to each other the relative merits of the observing strategy choices
that could be made, in an effort to maximize the scientific value of the
survey. The investigation of some science cases leads to suggestions for new
strategies that could be simulated and potentially adopted. Notably, we find
motivation for exploring departures from a spatially uniform annual tiling of
the sky: focusing instead on different parts of the survey area in different
years in a "rolling cadence" is likely to have significant benefits for a
number of time domain and moving object astronomy projects. The communal
assembly of a suite of quantified and homogeneously coded metrics is the vital
first step towards an automated, systematic, sciencebased assessment of any
given cadence simulation, that will enable the scheduling of the LSST to be as
wellinformed as possible.

The existence of massive ($10^{11}$ solar masses) elliptical galaxies by
redshift z~4 (when the Universe was 1.5 billion years old) necessitates the
presence of galaxies with starformation rates exceeding 100 solar masses per
year at z>6 (corresponding to an age of the Universe of less than 1 billion
years). Surveys have discovered hundreds of galaxies at these early cosmic
epochs, but their starformation rates are more than an order of magnitude
lower. The only known galaxies with very high starformation rates at z>6 are,
with only one exception, the host galaxies of quasars, but these galaxies also
host accreting supermassive (more than $10^9$ solar masses) black holes, which
probably affect the properties of the galaxies. Here we report observations of
an emission line of singly ionized carbon ([CII] at a wavelength of 158
micrometres) in four galaxies at z>6 that are companions of quasars, with
velocity offsets of less than 600 kilometers per second and linear offsets of
less than 600 kiloparsecs. The discovery of these four galaxies was
serendipitous; they are close to their companion quasars and appear bright in
the farinfrared. On the basis of the [CII] measurements, we estimate
starformation rates in the companions of more than 100 solar masses per year.
These sources are similar to the host galaxies of the quasars in [CII]
brightness, linewidth and implied dynamical masses, but do not show evidence
for accreting supermassive black holes. Similar systems have previously been
found at lower redshift. We find such close companions in four out of
twentyfive z>6 quasars surveyed, a fraction that needs to be accounted for in
simulations. If they are representative of the bright end of the [CII]
luminosity function, then they can account for the population of massive
elliptical galaxies at z~4 in terms of cosmic space density.

We study the spectrum of unbounded Jselfadjoint block operator matrices. In
particular, we prove enclosures for the spectrum, provide a sufficient
condition for the spectrum being real and derive variational principles for
certain real eigenvalues even in the presence of nonreal spectrum. The latter
lead to lower and upper bounds and asymptotic estimates for eigenvalues.

For a very general class of unbounded selfadjoint operator function we prove
upper bounds for eigenvalues which lie within arbitrary gaps of the essential
spectrum. These upper bounds are given by triple variations. Furthermore, we
find conditions which imply that a point is in the resolvent set. For norm
resolvent continuous operator functions we show that the variational inequality
becomes an equality.

We consider the problem of how to compute eigenvalues of a selfadjoint
operator when a direct application of the Galerkin (finitesection) method is
unreliable. The last two decades have seen the development of the socalled
quadratic methods for addressing this problem. Recently a new perturbation
approach has emerged, the idea being to perturb eigenvalues off the real line
and, consequently, away from regions where the Galerkin method fails. We
propose a simplified perturbation method which requires no \'a priori
information and for which we provide a rigorous convergence analysis. The
latter shows that, in general, our approach will significantly outperform the
quadratic methods. We also present a new spectral enclosure for operators of
the form $A+iB$ where $A$ is selfadjoint, $B$ is selfadjoint and bounded.
This enables us to control, very precisely, how eigenvalues are perturbed from
the real line. The main results are demonstrated with examples including
magnetohydrodynamics, Schr\"odinger and Dirac operators.

The Sloan Digital Sky Survey has surveyed 14,555 square degrees of the sky,
and delivered over a trillion pixels of imaging data. We present the
largescale clustering of 1.6 million quasars between z = 0.5 and z = 2.5 that
have been classified from this imaging, representing the highest density of
quasars ever studied for clustering measurements. This data set spans ~11,000
square degrees and probes a volume of 80(Gpc/h)^3. In principle, such a large
volume and medium density of tracers should facilitate highprecision
cosmological constraints. We measure the angular clustering of photometrically
classified quasars using an optimal quadratic estimator in four redshift slices
with an accuracy of ~25% over a bin width of l ~10  15 on scales corresponding
to matterradiation equality and larger (l ~ 2  30). Observational systematics
can strongly bias clustering measurements on large scales, which can mimic
cosmologically relevant signals such as deviations from Gaussianity in the
spectrum of primordial perturbations. We account for systematics by employing a
new method recently proposed by Agarwal et al. (2014) to the clustering of
photometrically classified quasars. We carefully apply our methodology to
mitigate known observational systematics and further remove angular bins that
are contaminated by unknown systematics. Combining quasar data with the
photometric luminous red galaxy (LRG) sample of Ross et al. (2011) and Ho et
al. (2012), and marginalizing over all bias and shot noiselike parameters, we
obtain a constraint on local primordial nonGaussianity of fNL = 113+/154
(1\sigma error). [Abridged]

A new technique for approximating eigenvalues and eigenvectors of a
selfadjoint operator is presented. The method does not incur spectral
pollution, uses trial spaces from the form domain, has a selfadjoint
algorithm, and exhibits superconvergence.

In the first part of this manuscript a relationship between the spectrum of
selfadjoint operator matrices and the spectra of their diagonal entries is
found. This leads to enclosures for spectral points and in particular,
enclosures for eigenvalues. We also consider graph invariant subspaces, and
their corresponding angular operators. The existence of a bounded angular
operator leads to basis properties of the first component of eigenvectors of
operator matrices for which the corresponding eigenvalues lie in a half line.
The results are applied to an example from magnetohydrodynamics.

The method of second order relative spectra has been shown to reliably
approximate the discrete spectrum for a selfadjoint operator. We extend the
method to normal operators and find optimal convergence rates for eigenvalues
and eigenspaces. The convergence to eigenspaces is new, while the convergence
rate for eigenvalues improves on the previous estimate by an order of
magnitude.

We consider the Galerkin method for approximating the spectrum of an operator
$T+A$ where $T$ is semibounded selfadjoint and $A$ satisfies a relative
compactness condition. We show that the method is reliable in all regions where
it is reliable for the unperturbed problem  which always contains
$\mathbb{C}\backslash\mathbb{R}$. The results lead to a new technique for
identifying eigenvalues of $T$, and for identifying spectral pollution which
arises from applying the Galerkin method directly to $T$. The new technique
benefits from being applicable on the form domain.

We present the design and performance of the multiobject fiber spectrographs
for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon
Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999
on the 2.5m aperture Sloan Telescope at Apache Point Observatory, the
spectrographs produced more than 1.5 million spectra for the SDSS and SDSSII
surveys, enabling a wide variety of Galactic and extragalactic science
including the first observation of baryon acoustic oscillations in 2005. The
spectrographs were upgraded in 2009 and are currently in use for BOSS, the
flagship survey of the thirdgeneration SDSSIII project. BOSS will measure
redshifts of 1.35 million massive galaxies to redshift 0.7 and Lymanalpha
absorption of 160,000 high redshift quasars over 10,000 square degrees of sky,
making percent level measurements of the absolute cosmic distance scale of the
Universe and placing tight constraints on the equation of state of dark energy.
The twin multiobject fiber spectrographs utilize a simple optical layout
with reflective collimators, gratings, allrefractive cameras, and
stateoftheart CCD detectors to produce hundreds of spectra simultaneously in
two channels over a bandpass covering the near ultraviolet to the near
infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven
heritage, the spectrographs were upgraded for BOSS with volumephase
holographic gratings and modern CCD detectors, improving the peak throughput by
nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000
nm, and increasing the number of fibers from 640 to 1000 per exposure. In this
paper we describe the original SDSS spectrograph design and the upgrades
implemented for BOSS, and document the predicted and measured performances.

We discuss how to compute certified enclosures for the eigenvalues of
benchmark linear magnetohydrodynamics operators in the plane slab and
cylindrical pinch configurations. For the plane slab, our method relies upon
the formulation of an eigenvalue problem associated to the Schur complement,
leading to highly accurate upper bounds for the eigenvalue. For the cylindrical
configuration, a direct application of this formulation is possible, however,
it cannot be rigourously justified. Therefore in this case we rely on a
specialized technique based on a method proposed by Zimmermann and Mertins. In
turns this technique is also applicable for finding accurate complementary
bounds in the case of the plane slab. We establish convergence rates for both
approaches.

Let A be a selfadjoint operator acting on a Hilbert space. The notion of
second order spectrum of A relative to a given finitedimensional subspace L
has been studied recently in connection with the phenomenon of spectral
pollution in the Galerkin method. We establish in this paper a general
framework allowing us to determine how the second order spectrum encodes
precise information about the multiplicity of the isolated eigenvalues of A.
Our theoretical findings are supported by various numerical experiments on the
computation of inclusions for eigenvalues of benchmark differential operators
via finite element bases.

Promoting racial and ethnic diversity is critically important to the future
success and growth of the field of astronomy. The raw ability, drive and
interest required to excel in the field is distributed without regard to race,
gender, or socioeconomic background. By not actively promoting diversity in our
field we risk losing talented people to other professions (or losing them
entirely), which means that there will be astronomical discoveries that simply
won't get made. There is demonstrated evidence that STEM fields benefit from
diverse perspectives on problems that require more complex thought processes.
This is especially relevant to a field like astronomy where more and more work
is being done collaboratively. The lack of notable growth in African American,
Hispanic, and Native American representation in astronomy indicates that the
'pipeline' for these individuals is systemically leaky at critical junctures.
Substantially more effort must be directed toward improving the educational and
career development of minorities to insure that these potential colleagues are
supported through the process. However, simply recognizing that the pipeline is
faulty is woefully inadequate. There must be very specific, targeted solutions
to help improve the situation. With this in mind, we offer two position papers
addressing specific areas of improvement that we identify as (a) essential for
any foreseeable progress in the field, and (b) attainable in the 20102020
decade. These position papers focus primarily on African Americans, Hispanics,
and Native Americans. Although we do not directly address issues of Asian
Americans, Pacific Islanders, and other groups, many of the recommendations
made here can be adapted to address issues faced by these groups as well.

If the ethnic makeup of the astronomy profession is to achieve parity with
the general population within one generation (~30 years), the number of
underrepresented minorities earning graduate degrees in astronomy and
astrophysics must increase in the coming decade by a factor of 5 to 10. To
accomplish this, the profession must develop and invest in mechanisms to more
effectively move individuals across critical educational junctures to the PhD
and beyond. Early and continuous research engagement starting in the
undergraduate years is critical to this vision, in which the federally funded
research internship programs (e.g. NSF REU, NASA GSRP) and national
centers/observatories play a vital role. Regionally based partnerships with
minorityserving institutions (MSIs) are crucial for tapping extant pools of
minority talent, as are postbaccalaurate and/or masters degree "bridging"
programs that provide critical stepping stones to the PhD. Because of the
strong undergraduate physics, engineering, and computer science backgrounds of
many students from MSIs, we suggest that instrument development and large scale
computing/datamining are particularly promising avenues for engagement in the
coming decade.

In order to attract, recruit and retain underrepresented minority students to
pursue Astronomy and related fields, we must ensure that there continues to be
a well qualified pool of graduate and undergraduate students from which to
recruit. This required pool of people are today's elementary, middle and high
school students. The Astronomy community must be proactive in demonstrating the
importance of pursing scientific study and careers to these students and their
parents. Only by actively engaging these communities can U.S Astronomy hope to
increase the numbers of minority PhDs and continue to be a leader in
Astronomical discovery and knowledge.

Over the last decade optical spectroscopic surveys have characterized the low
redshift galaxy population and uncovered populations of starforming galaxies
back to z ~ 7. This work has shown that the primary epoch of galaxy building
and black hole growth occurs at redshifts of 2 to 3. The establishment of the
concordance LCDM cosmology shifted the focus of galaxy population studies from
constraining cosmological parameters to characterizing the processes which
regulate the formation and evolution of galaxies.In the next decade, high
redshift observers will attempt to formulate a coherent evolutionary picture
connecting galaxies in the high redshift Universe to galaxies today. In order
to link galaxy populations at different redshifts, we must not only
characterize their evolution in a systematic way, we must establish which
physical processes are responsible for it. Considerable progress has already
been made in understanding how galaxies evolved from z ~ 1 to the present day.
Large spectroscopic surveys in the near infrared are required to push these
studies back towards the main epoch of galaxy building. Only then will we
understand the full story of the formation of L* galaxies like our own Milky
Way. A large nearIR spectroscopic survey will also provide the calibration
needed to avoid systematics in the large photometric programs proposed to study
the nature of dark matter and dark energy. We provide an outline design for a
multiobject 0.4 to 1.8 micron spectrograph, which could be placed on an
existing telescope, and which would allow a full characterization of the galaxy
population out to z ~ 2. We strongly recommend a serious further study to
design a real instrument, which will be required for galaxy formation studies
to advance to the next frontier.

Much of the science case for the next generation of deep, widefield
optical/infrared surveys has been driven by the further study of dark energy.
This is a laudable goal (and the subject of a companion white paper by Zhan et
al.). However, one of the most important lessons of the current generation of
surveys is that the interesting science questions at the end of the survey are
quite different than they were when the surveys were being planned. The current
surveys succeeded in this evolving terrain by being very general tools that
could be applied to a number of very fundamental measurements. Likewise, the
accessibility of the data enabled the broader cosmological and astronomical
community to generate more science than the survey collaborations could alone.
With that in mind, we should consider some of the basic physical and
cosmological questions that surveys like LSST and JDEMWide will be able to
address.

The coming decade will be an exciting period for dark energy research, during
which astronomers will address the question of what drives the accelerated
cosmic expansion as first revealed by type Ia supernova (SN) distances, and
confirmed by later observations.
The mystery of dark energy poses a challenge of such magnitude that, as
stated by the Dark Energy Task Force (DETF), "nothing short of a revolution in
our understanding of fundamental physics will be required to achieve a full
understanding of the cosmic acceleration." The lack of multiple complementary
precision observations is a major obstacle in developing lines of attack for
dark energy theory. This lack is precisely what nextgeneration surveys will
address via the powerful techniques of weak lensing (WL) and baryon acoustic
oscillations (BAO)  galaxy correlations more generally  in addition to SNe,
cluster counts, and other probes of geometry and growth of structure. Because
of their unprecedented statistical power, these surveys demand an accurate
understanding of the observables and tight control of systematics.
This white paper highlights the opportunities, approaches, prospects, and
challenges relevant to dark energy studies with widedeep multiwavelength
photometric redshift surveys. Quantitative predictions are presented for a
20000 sq. deg. groundbased 6band (ugrizy) survey with 5sigma depth of
r~27.5, i.e., a Stage 4 survey as defined by the DETF.

It is widely expected that the coming decade will witness the first direct
detection of gravitational waves (GWs). The groundbased LIGO and Virgo GW
observatories are being upgraded to advanced sensitivity, and are expected to
observe a significant binary merger rate. The launch of The Laser
Interferometer Space Antenna (LISA) would extend the GW window to low
frequencies, opening new vistas on dynamical processes involving massive (M >~
10^5 M_Sun) black holes. GW events are likely to be accompanied by
electromagnetic (EM) counterparts and, since information carried
electromagnetically is complementary to that carried gravitationally, a great
deal can be learned about an event and its environment if it becomes possible
to measure both forms of radiation in concert. Measurements of this kind will
mark the dawn of transspectral astrophysics, bridging two distinct spectral
bands of information. The aim of this whitepaper is to articulate future
directions in both theory and observation that are likely to impact broad
astrophysical inquiries of general interest. What will EM observations reflect
on the nature and diversity of GW sources? Can GW sources be exploited as
complementary probes of cosmology? What crossfacility coordination will expand
the science returns of gravitational and electromagnetic observations?

In this paper we discuss the stability of an alternative pollutionfree
procedure for computing spectra. The main difference with the Galerkin method
lies in the fact that it gives rise to a weak approximate problem which is
quadratic in the spectral parameter, instead of linear. Previous accounts on
this new procedure can be found in Levitin and Shargorodsky (2002)
[math.SP/0212087] and Boulton (2006) [math.SP/0503126].

The Sloan Digital Sky Survey can be used to detect and characterize red and
reddened quasars. In Richards et al. (2003), we showed that 6% of SDSS quasars
have red colors consistent with significant dust reddening by an extinction
curve similar to that of the Small Magellanic Cloud (SMC). We estimate that a
further 10% of the luminous quasar population is missing from our
magnitudelimited SDSS sample. More recent work (Hopkins et al., in
preparation) confirms the conclusion that dust reddening is the primary
explanation for SDSS quasars in the red tail of the color distribution. Fitting
orthogonal first and secondorder polynomials to SDSS quasar photometry
measures the slope and curvature of each object's UV/optical spectrum. The
slope vs. curvature distribution is elongated along the axis predicted for
SMClike reddening, while the axes predicted for LMC or MWlike reddening
provide significantly poorer fits. Extension to longer wavelengths using a
smaller sample of SDSS/2MASS matches confirms this result at high significance.

We compute the angular power spectrum C_l from 1.5 million galaxies in early
SDSS data on large angular scales, l<600. The data set covers about 160 square
degrees, with a characteristic depth of order 1 Gpc/h in the faintest (21<r<22)
of our four magnitude bins. Cosmological interpretations of these results are
presented in a companion paper by Dodelson et al (2001). The data in all four
magnitude bins are consistent with a simple flat ``concordance'' model with
nonlinear evolution and linear bias factors of order unity. Nonlinear evolution
is particularly evident for the brightest galaxies. A series of tests suggest
that systematic errors related to seeing, reddening, etc., are negligible,
which bodes well for the sixtyfold larger sample that the SDSS is currently
collecting. Uncorrelated error bars and wellbehaved window functions make our
measurements a convenient starting point for cosmological model fitting.

Precision measurements of the galaxy power spectrum P(k) require a data
analysis pipeline that is both fast enough to be computationally feasible and
accurate enough to take full advantage of highquality data. We present a
rigorous discussion of different methods of power spectrum estimation, with
emphasis on the traditional Fourier method, the linear (KarhunenLoeve; KL),
and quadratic data compression schemes, showing in what approximations they
give the same result. To improve speed, we show how many of the advantages of
KL data compression and power spectrum estimation may be achieved with a
computationally faster quadratic method. To improve accuracy, we derive
analytic expressions for handling the integral constraint, since it is crucial
that finite volume effects are accurately corrected for on scales comparable to
the depth of the survey. We also show that for the KL and quadratic techniques,
multiple constraints can be included via simple matrix operations, thereby
rendering the results less sensitive to galactic extinction and misestimates
of the radial selection function. We present a data analysis pipeline that we
argue does justice to the increases in both quality and quantity of data that
upcoming redshift surveys will provide. It uses three analysis techniques in
conjunction: a traditional Fourier approach on small scales, a pixelized
quadratic matrix method on large scales and a pixelized KL eigenmode analysis
to probe anisotropic effects such as redshiftspace distortions.

We employed the Very Large Array (VLA) of the National Radio Astronomy
Observatory in C configuration to map 39 ultraluminous IRAS galaxies at 6~cm
and 20~cm, at resolutions of ~ 4" and 15", respectively, and 24 sources at 6~cm
with in the A configuration with a resolution of ~0.5". Most of the sources
have radio spectral indices indicative of synchrotron emission (alpha ~ 0.65).
There is one source, however, that shows an inverted spectrum with alpha =
+2.1; observations at higher frequencies show that the spectrum peaks between 5
and 8 GHz, as high as any of the ``gigahertz peaked spectrum'' sources studied
by O'Dea etal. We discuss the implications of this source for observations of
fluctuations in the CMB.
Two of the sources show multiple unresolved components, another four are
doubles with at least one resolved component, 14 show extended emission which
could arise from a disk, and two show arcsecond long jets. Our data fit the
tight correlation found by Helou etal (1985) between farinfrared and microwave
luminosity; this correlation extends to the highest infrared luminosities. The
correlation is weaker if only the extended or the nuclear components of the
radio luminosity are used. Therefore the farinfrared emission in the majority
of these higher luminosity galaxies is due to the same mechanism as the lower
luminosity FIR sources, which is believed to be star formation, rather than
nonthermal activity in the nucleus. Moreover, the star formation is not
confined to the extended disk in these sources, but is important in the nucleus
as well.