
The fourth generation of the Sloan Digital Sky Survey (SDSSIV) has been in
operation since July 2014. This paper describes the second data release from
this phase, and the fourteenth from SDSS overall (making this, Data Release
Fourteen or DR14). This release makes public data taken by SDSSIV in its first
two years of operation (July 20142016). Like all previous SDSS releases, DR14
is cumulative, including the most recent reductions and calibrations of all
data taken by SDSS since the first phase began operations in 2000. New in DR14
is the first public release of data from the extended Baryon Oscillation
Spectroscopic Survey (eBOSS); the first data from the second phase of the
Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE2),
including stellar parameter estimates from an innovative data driven machine
learning algorithm known as "The Cannon"; and almost twice as many data cubes
from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous
release (N = 2812 in total). This paper describes the location and format of
the publicly available data from SDSSIV surveys. We provide references to the
important technical papers describing how these data have been taken (both
targeting and observation details) and processed for scientific use. The SDSS
website (www.sdss.org) has been updated for this release, and provides links to
data downloads, as well as tutorials and examples of data use. SDSSIV is
planning to continue to collect astronomical data until 2020, and will be
followed by SDSSV.

Finite temperature corrections to the effective potential and the
energymomentum tensor of a scalar field are computed in a perturbed Minkoswki
spacetime. We consider the explicit mode decomposition of the field in the
perturbed geometry and obtain analytical expressions in the nonrelativistic
and ultrarelativistic limits to first order in scalar metric perturbations. In
the static case, our results are in agreement with previous calculations based
on the SchwingerDe Witt expansion which indicate that thermal effects in a
curved spacetime can be encoded in the local Tolman temperature at leading
order in perturbations and in the adiabatic expansion. We also study the shift
of the effective potential minima produced by thermal corrections in the
presence of static gravitational fields. Finally we discuss the dependence on
the initial conditions set for the mode solutions.

Cosmological simulations are fundamental tools to study structure formation
and the astrophysics of evolving structures, in particular clusters of
galaxies. While hydrodynamical simulations cannot sample efficiently large
volumes and explore different cosmologies at the same time, Nbody simulations
lack the baryonic physics that is crucial to determine the observed properties
of clusters. One solution is to use (semi)analytical models to implement the
needed baryonic physics. In this way, we can generate the many mock universes
that will be required to fully exploit future large sky surveys, such as that
from the upcoming eROSITA Xray telescope. We developed a phenomenological
model based on observations of clusters to implement gas density and
temperature information on the darkmatteronly halos of the MultiDark
simulations. We generate several fullsky mock lightcones of clusters for the
WMAP and Planck cosmologies, adopting different parameters in our
phenomenological model of the intracluster medium. For one of these
simulations and models, we also generate 100 lightcones corresponding to 100
random observers and explore the variance among them in several quantities. In
this first paper on MultiDark mock galaxy cluster lightcones, we focus on
presenting our methodology and discuss predictions for eROSITA, in particular,
exploring the potential of angular power spectrum analyses of its detected (and
undetected) cluster population to study Xray scaling relations, the
intracluster medium, and the composition of the cosmic Xray background. We
make publicly available online more than 400 GB of lightcones, which include
the expected eROSITA count rate, on Skies & Universes
(http://www.skiesanduniverses.org).

We further develop the description of redshift space distortions within the
Effective Field Theory of Large Scale Structures. First, we generalize the
counterterms to include the effect of baryonic physics and primordial
nonGaussianity. Second, we evaluate the IRresummation of the dark matter
power spectrum in redshift space. This requires us to identify a controlled
approximation that makes the numerical evaluation straightforward and
efficient. Third, we compare the predictions of the theory at one loop with the
power spectrum from numerical simulations up to $\ell=6$. We find that the
IRresummation allows us to correctly reproduce the BAO peak. The $k$reach, or
equivalently the precision for a given $k$, depends on additional counterterms
that need to be matched to simulations. Since the nonlinear scale for the
velocity is expected to be longer than the one for the overdensity, we consider
a minimal and a nonminimal set of counterterms. The quality of our numerical
data makes it hard to firmly establish the performance of the theory at high
wavenumbers. Within this limitation, we find that the theory at redshift
$z=0.56$ and up to $\ell=2$ matches the data to percent level approximately up
to $k \sim 0.13 \, h { \rm Mpc^{1}}$ or $k \sim 0.18 \, h { \rm Mpc^{1}}$,
depending on the number of counterterms used, with potentially large
improvement over former analytical techniques.

We investigate the potential sources of theoretical systematics in the
anisotropic Baryon Acoustic Oscillation (BAO) distance scale measurements from
the clustering of galaxies in configuration space using the final Data Release
(DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS). We perform a
detailed study of the impact on BAO measurements from choices in the
methodology such as fiducial cosmology, clustering estimators, random
catalogues, fitting templates, and covariance matrices.
The theoretical systematic uncertainties in BAO parameters are found to be
0.002 in the isotropic dilation $\alpha$ and 0.003 in the quadrupolar dilation
$\epsilon$. The leading source of systematic uncertainty is related to the
reconstruction techniques. Theoretical uncertainties are subdominant compared
with the statistical uncertainties for BOSS survey, accounting
$0.2\sigma_{stat}$ for $\alpha$ and $0.25\sigma_{stat}$ for $\epsilon$
($\sigma_{\alpha,stat} \sim$0.010 and $\sigma_{\epsilon,stat}\sim$ 0.012
respectively). We also present BAOonly distance scale constraints from the
anisotropic analysis of the correlation function. Our constraints on the
angular diameter distance $D_A(z)$ and the Hubble parameter $H(z)$, including
both statistical and theoretical systematic uncertainties, are 1.5\% and 2.8\%
at $z_{\rm eff}=0.38$, 1.4\% and 2.4\% at $z_{\rm eff}=0.51$, and 1.7\% and
2.6\% at $z_{\rm eff}=0.61$. This paper is part of a set that analyzes the
final galaxy clustering dataset from BOSS. The measurements and likelihoods
presented here are crosschecked with other BAO analysis in \citet{Acacia16}.
The systematic error budget concerning the methodology on postreconstruction
BAO analysis presented here is used in \citet{Acacia16} to produce the final
cosmological constraints from BOSS.

We use the DECaLS DR3 survey photometry matched to the SDSSIII/BOSS DR12
spectroscopic catalog to investigate the morphology and stellar masssize
relation of luminous red galaxies (LRGs) within the CMASS and LOWZ galaxy
samples in the redshift range $0.2<z<0.7$. The large majority of both samples
is composed of earlytype galaxies with De Vaucouleurs profiles, while only
less than 20% are latetype exponentials. We calibrate DECaLS effective radii
using the higher resolution CFHT/MegaCam observations and optimise the
correction for each morphological type. By crossmatching the photometric
properties of the earlytype population with the Portsmouth stellar mass
catalog, we are able to explore the highmass end of the distribution using a
large sample of 313,026 galaxies over 4380 deg$^{2}$. We find a clear
correlation between the sizes and the stellar masses of these galaxies, which
appears flatter than previous estimates at lower masses. The sizes of these
earlytype galaxies do not exhibit significant evolution within the BOSS
redshift range, but a slightly declining redshift trend is found when these
results are combined with $z\sim0.1$ SDSS measurements at the highmass end.
The synergy between BOSS and DECaLS has important applications in other fields,
including galaxy clustering and weak lensing.

We present the clustering measurements of quasars in configuration space
based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended
Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars
spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square
degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach
with a Gaussian Streaming (GS) model for the redshift space distortions of the
correlation function and demonstrate its applicability for dark matter halos
hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$,
we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426
\pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{13}(r_{s}^{\rm
fid}/r_s){\rm km.s}^{1}.{\rm Mpc}^{1}$, and the angular diameter distance
$D_{A}(z_{\rm eff})=1850^{+90}_{115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where
$r_{s}$ is the sound horizon at the end of the baryon drag epoch and
$r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors
include both systematic and statistical contributions. The results on the
evolution of distances are consistent with the predictions of flat
$\Lambda$Cold Dark Matter ($\Lambda$CDM) cosmology with Planck parameters,
and the measurement of $f\sigma_{8}$ extends the validity of General Relativity
(GR) to higher redshifts($z>1$) This paper is released with companion papers
using the same sample. The results on the cosmological parameters of the
studies are found to be in very good agreement, providing clear evidence of the
complementarity and of the robustness of the first fullshape clustering
measurements with the eBOSS DR14 quasar sample.

We present the public release of the MultiDarkGalaxies: three distinct
galaxy catalogues derived from one of the Planck cosmology MultiDark
simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution
of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semianalytic models
GALACTICUS, SAG, and SAGE to it. We compare the three models and their
conformity with observational data for a selection of fundamental properties of
galaxies like stellar mass function, star formation rate, cold gas fractions,
and metallicities  noting that they sometimes perform differently reflecting
model designs and calibrations. We have further selected galaxy subsamples of
the catalogues by number densities in stellar mass, cold gas mass, and star
formation rate in order to study the clustering statistics of galaxies. We show
that despite different treatment of orphan galaxies, i.e. galaxies that lost
their darkmatter host halo due to the finite mass resolution of the Nbody
simulation or tidal stripping, the clustering signal is comparable, and
reproduces the observations in all three models  in particular when selecting
samples based upon stellar mass. Our catalogues provide a powerful tool to
study galaxy formation within a volume comparable to those probed by ongoing
and future photometric and redshift surveys. All model data consisting of a
range of galaxy properties  including broadband SDSS magnitudes  are
publicly available.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in
redshiftspace using the clustering of quasars. We consider a sample of 147,000
quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS)
distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure
their sphericallyaveraged clustering in both configuration and Fourier space.
Our observational dataset and the 1400 simulated realizations of the dataset
allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We
determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent
precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc.
This is the first time the location of the BAO feature has been measured
between redshifts 1 and 2. Our result is fully consistent with the prediction
obtained by extrapolating the Planck flat $\Lambda$CDM bestfit cosmology. All
of our results are consistent with basic largescale structure (LSS) theory,
confirming quasars to be a reliable tracer of LSS, and provide a starting point
for numerous cosmological tests to be performed with eBOSS quasar samples. We
combine our result with previous, independent, BAO distance measurements to
construct an updated BAO distanceladder. Using these BAO data alone and
marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda}
> 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free
curvature.

The fourth generation of the Sloan Digital Sky Survey (SDSSIV) began
observations in July 2014. It pursues three core programs: APOGEE2, MaNGA, and
eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS.
This paper describes the first data release from SDSSIV, Data Release 13
(DR13), which contains new data, reanalysis of existing data sets and, like all
SDSS data releases, is inclusive of previously released data. DR13 makes
publicly available 1390 spatially resolved integral field unit observations of
nearby galaxies from MaNGA, the first data released from this survey. It
includes new observations from eBOSS, completing SEQUELS. In addition to
targeting galaxies and quasars, SEQUELS also targeted variabilityselected
objects from TDSS and Xray selected objects from SPIDERS. DR13 includes new
reductions of the SDSSIII BOSS data, improving the spectrophotometric
calibration and redshift classification. DR13 releases new reductions of the
APOGEE1 data from SDSSIII, with abundances of elements not previously
included and improved stellar parameters for dwarf stars and cooler stars. For
the SDSS imaging data, DR13 provides new, more robust and precise photometric
calibrations. Several valueadded catalogs are being released in tandem with
DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and
an updated redclump catalog for APOGEE. This paper describes the location and
format of the data now publicly available, as well as providing references to
the important technical papers that describe the targeting, observing, and data
reduction. The SDSS website, http://www.sdss.org, provides links to the data,
tutorials and examples of data access, and extensive documentation of the
reduction and analysis procedures. DR13 is the first of a scheduled set that
will contain new data and analyses from the planned ~6year operations of
SDSSIV.

We introduce the Lomonosov suite of highresolution Nbody cosmological
simulations covering a full box of size 32 $h^{1}$ Mpc with lowmass
resolution particles ($2\times10^7$ $h^{1}\,M_\odot$) and three zoomin
simulations of overdense, underdense and mean density regions at much higher
particle resolution ($4\times10^4$ $h^{1}\,M_\odot$). The main purpose of this
simulation suite is to extend the concentrationmass relation of dark matter
halos down to masses below those typically available in large cosmological
simulations. The three different density regions available at higher resolution
provide a better understanding of the effect of the local environment on halo
concentration, known to be potentially important for small simulation boxes and
small halo masses. Yet, we find the correction to be small in comparison with
the scatter of halo concentrations. We conclude that zoom simulations, despite
their limited representativity of the volume of the Universe, can be
effectively used for the measurement of halo concentrations at least at the
halo masses probed by our simulations. In any case, after a precise
characterization of this effect, we develop a robust technique to extrapolate
the concentration values found in zoom simulations to larger volumes with
greater accuracy. All together, Lomonosov provides a measure of the
concentrationmass relation in the halo mass range $10^710^{10}$
$h^{1}\,M_\odot$ with superb halo statistics. This work represents a first
important step to measure halo concentrations at intermediate, yet vastly
unexplored halo mass scales, down to the smallest ones. All Lomonosov data and
files are public for community's use.

We present a new methodology to generate mock halo or galaxy catalogues,
which have accurate clustering properties, nearly indistinguishable from full
$N$body solutions, in terms of the onepoint, twopoint, and threepoint
statistics. In particular, the agreement is remarkable, within $1\%$ up to
$k=0.55$ $h$Mpc$^{1}$ and down to $r=10$ $h^{1}$Mpc, for the power spectrum
and twopoint correlation function respectively, while the bispectrum agrees in
general within $20\%$ for different scales and shapes. Our approach is based on
the Zel'dovich approximation, however, effectively including with the simple
prescriptions the missing physical ingredients, and stochastic scaledependent,
nonlocal and nonlinear biasing contributions. The computing time and memory
required to produce one mock is similar to that using the lognormal model.
With high accuracy and efficiency, the effective Zel'dovich approximation mocks
(EZmocks) provide a reliable and practical method to produce massive mock
galaxy catalogues for the analysis of largescale structure measurements.

We present measurements of the anisotropic galaxy clustering from the Data
Release 9 (DR9) CMASS sample of the SDSSIII Baryon Oscillation Spectroscopic
Survey (BOSS). We analyze the broadrange shape of the monopole and quadrupole
correlation functions to obtain constraints, at the effective redshift $z=0.57$
of the sample, on the Hubble expansion rate $H(z)$, the angulardiameter
distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, the physical
matter density $\Omega_m h^2$, and the biased amplitude of matter fluctuation
b\sigma_8(z). We obtain {$H(0.57)$, $D_A(0.57)$, $f(0.57)\sigma_8(0.57)$,
$\Omega_m h^2$, $b\sigma_8(0.57)$} = {$87.6_{6.8}^{+6.7}$, $1396\pm73$,
$0.126_{0.010}^{+0.008}$, $1.19\pm0.14$, $0.428\pm0.066$} and their covariance
matrix as well. The parameters which are not well constrained by our of galaxy
clustering analysis are marginalized over with wide flat priors. Since no
priors from other data sets (i.e., CMB) are adopted and no dark energy models
are assumed, our results from BOSS CMASS galaxy clustering alone may be
combined with other data sets, i.e. CMB, SNe, lensing or other galaxy
clustering data to constrain the parameters of a given cosmological model. We
show that the major power on constraining dark energy from the anisotropic
galaxy clustering signal, as compared to the angularaveraged one (monopole),
arises from including the normalized growth rate $f(z)\sigma_8(z)$. In the case
of the wCDM cosmological model our singleprobe CMASS constraints, combined
with CMB (WMAP9+SPT), yield a value for the dark energy equation of state
parameter of $w=0.90\pm0.11$. Therefore, it is important to include
$f(z)\sigma_8(z)$ while investigating the nature of dark energy with current
and upcoming largescale galaxy surveys.

We study the galaxy clustering dependence on the
$\left[\mathrm{O\scriptsize{II}}\right]$ emission line luminosity in the SDSS
DR7 Main galaxy sample at mean redshift $z\sim0.1$. We select volumelimited
samples of galaxies with different $\left[\mathrm{O\scriptsize{II}}\right]$
luminosity thresholds and measure their projected, monopole and quadrupole
twopoint correlation functions. We model these observations using the
1$h^{1}\rm{Gpc}$ MultiDark Planck cosmological simulation and generate
lightcones with the SUrvey GenerAtoR algorithm. To interpret our results, we
adopt a modified (Sub)Halo Abundance Matching scheme, accounting for the
stellar mass incompleteness of the emission line galaxies. The satellite
fraction constitutes an extra parameter in this model and allows to optimize
the clustering fit on both small and intermediate scales (i.e. $r_p\lesssim
30h^{1}\rm{Mpc})$, with no need of any velocity bias correction. We find that,
in the local Universe, the $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity
correlates with all the clustering statistics explored and with the galaxy
bias. This latter quantity correlates more strongly with the SDSS $r$band
magnitude than $\left[\mathrm{O\scriptsize{II}}\right]$ luminosity. In
conclusion, we propose a straightforward method to produce reliable clustering
models, entirely built on the simulation products, which provides robust
predictions of the typical ELG host halo masses and satellite fraction values.
The SDSS galaxy data, MultiDark mock catalogues and clustering results are made
publicly available.

We describe the Sloan Digital Sky Survey IV (SDSSIV), a project encompassing
three major spectroscopic programs. The Apache Point Observatory Galactic
Evolution Experiment 2 (APOGEE2) is observing hundreds of thousands of Milky
Way stars at high resolution and high signaltonoise ratio in the
nearinfrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA)
survey is obtaining spatiallyresolved spectroscopy for thousands of nearby
galaxies (median redshift of z = 0.03). The extended Baryon Oscillation
Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas
distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using
baryon acoustic oscillations, redshift space distortions, and the shape of the
power spectrum. Within eBOSS, we are conducting two major subprograms: the
SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating Xray
AGN and galaxies in Xray clusters, and the Time Domain Spectroscopic Survey
(TDSS), obtaining spectra of variable sources. All programs use the 2.5meter
Sloan Foundation Telescope at Apache Point Observatory; observations there
began in Summer 2014. APOGEE2 also operates a second nearinfrared
spectrograph at the 2.5meter du Pont Telescope at Las Campanas Observatory,
with observations beginning in early 2017. Observations at both facilities are
scheduled to continue through 2020. In keeping with previous SDSS policy,
SDSSIV provides regularly scheduled public data releases; the first one, Data
Release 13, was made available in July 2016.

The standard relation between the cosmological redshift and cosmic scale
factor underlies cosmological inference from virtually all kinds of
cosmological observations, leading to the emergence of the LambdaCDM
cosmological model. This relation is not a fundamental theory and thus
observational determination of this function (redshift remapping) should be
regarded as an insightful alternative to holding its standard form in analyses
of cosmological data. Here we present nonparametric reconstructions of
redshift remapping in darkmatterdominated models and constraints on
cosmological parameters from a joint analysis of all primary cosmological
probes including the local measurement of the Hubble constant, Type Ia
supernovae, baryonic acoustic oscillations (BAO), Planck observations of the
cosmic microwave background (CMB) radiation (temperature power spectrum) and
cosmic chronometers. The reconstructed redshift remapping points to an
additional boost of redshift operating in late epoch of cosmic evolution, but
affecting both lowredshift observations and the CMB. The model predicts a
significant difference between the actual Hubble constant, h=0.48+/0.02, and
its local determination, h_obs=0.73+/0.02. The ratio of these two values
coincides closely with the maximum expansion rate inside voids formed in the
corresponding open cosmological model with Omega_m=0.87+/0.03, whereas the
actual value of the Hubble constant implies the age of the Universe that is
compatible with the Planck LambdaCDM cosmology. The new darkmatterdominated
model with redshift remapping provides excellent fits to all data and
eliminates recently reported tensions between the Planck LambdaCDM cosmology,
the local determination of the Hubble constant and the BAO measurements from
the Lyalpha forest of highredshift quasars.

The onepoint probability distribution function (PDF) of the matter density
field in the universe is a fundamental property that plays an essential role in
cosmology for estimates such as gravitational weak lensing, nonlinear
clustering, massive production of mock galaxy catalogs, and testing predictions
of cosmological models. Here we make a comprehensive analysis of the dark
matter PDF using a suite of 7000 Nbody simulations that covers a wide range of
numerical and cosmological parameters. We find that the PDF has a simple shape:
it declines with density as a powerlaw P~rho**(2), which is exponentially
suppressed on both small and large densities. The proposed doubleexponential
approximation provides an accurate fit to all our Nbody results for small
filtering scales R< 5Mpc/h with rms density fluctuations sigma>1. In
combination with the spherical infall model that works well for small
fluctuations sigma<1, the PDF is now approximated with just few percent errors
over the range of twelve orders of magnitude  a remarkable example of
precision cosmology. We find that at 510% level the PDF explicitly depends on
redshift (at fixed sigma) and on cosmological density parameter Omega_m. We
test different existing analytical approximations and find that the often used
lognormal approximation is always 35 times less accurate than either the
doubleexponential approximation or the spherical infall model.

$N$body cosmological simulations are an essential tool to understand the
observed distribution of galaxies. We use the MultiDark simulation suite, run
with the Planck cosmological parameters, to revisit the mass and velocity
functions. At redshift $z=0$, the simulations cover four orders of magnitude in
halo mass from $\sim10^{11}M_\odot$ with 8,783,874 distinct halos and 532,533
subhalos. The total volume used is $\sim$515 Gpc$^3$, more than 8 times larger
than in previous studies. We measure and model the halo mass function, its
covariance matrix w.r.t halo mass and the large scale halo bias. With the
formalism of the excursionset mass function, we explicit the tight
interconnection between the covariance matrix, bias and halo mass function. We
obtain a very accurate ($<2\%$ level) model of the distinct halo mass function.
We also model the subhalo mass function and its relation to the distinct halo
mass function. The set of models obtained provides a complete and precise
framework for the description of halos in the concordance Planck cosmology.
Finally, we provide precise analytical fits of the $V_{max}$ maximum velocity
function up to redshift $z<2.3$ to push for the development of halo occupation
distribution using $V_{max}$. The data and the analysis code are made publicly
available in the \textit{Skies and Universes} database.

We study the first year of the eBOSS quasar sample in the redshift range
$0.9<z<2.2$ which includes 68,772 homogeneously selected quasars. We show that
the main source of systematics in the evaluation of the correlation function
arises from inhomogeneities in the quasar target selection, particularly
related to the extinction and depth of the imaging data used for targeting. We
propose a weighting scheme that mitigates these systematics. We measure the
quasar correlation function and provide the most accurate measurement to date
of the quasar bias in this redshift range, $b_Q = 2.45 \pm 0.05$ at $\bar
z=1.55$, together with its evolution with redshift. We use this information to
determine the minimum mass of the halo hosting the quasars and the
characteristic halo mass, which we find to be both independent of redshift
within statistical error. Using a recentlymeasured quasarluminosityfunction
we also determine the quasar duty cycle. The size of this first year sample is
insufficient to detect any luminosity dependence to quasar clustering and this
issue should be further studied with the final $\sim$500,000 eBOSS quasar
sample.

We analyze the spectra of 300,000 luminous red galaxies (LRGs) with stellar
masses $M_* \gtrsim 10^{11} M_{\odot}$ from the SDSSIII Baryon Oscillation
Spectroscopic Survey (BOSS). By studying their starformation histories, we
find two main evolutionary paths converging into the same quiescent galaxy
population at $z\sim0.55$. Fastgrowing LRGs assemble $80\%$ of their stellar
mass very early on ($z\sim5$), whereas slowgrowing LRGs reach the same
evolutionary state at $z\sim1.5$. Further investigation reveals that their
clustering properties on scales of $\sim$130 Mpc are, at a high level of
significance, also different. Fastgrowing LRGs are found to be more strongly
clustered and reside in overall denser largescale structure environments than
slowgrowing systems, for a given stellarmass threshold. Our results imply a
dependence of clustering on stellarmass assembly history (naturally connected
to the massformation history of the corresponding halos) for a homogeneous
population of similar mass and color, which constitutes a strong observational
evidence of galaxy assembly bias.

We investigate the cosmological implications of studying galaxy clustering
using a tomographic approach applied to the final BOSS DR12 galaxy sample,
including both auto and crosscorrelation functions between redshift shells.
We model the signal of the full shape of the angular correlation function,
$\omega(\theta)$, in redshift bins using stateoftheart modelling of
nonlinearities, bias and redshiftspace distortions. We present results on the
redshift evolution of the linear bias of BOSS galaxies, which cannot be
obtained with traditional methods for galaxyclustering analysis. We also
obtain constraints on cosmological parameters, combining this tomographic
analysis with measurements of the cosmic microwave background (CMB) and type Ia
supernova (SNIa). We explore a number of cosmological models, including the
standard $\Lambda$CDM model and its most interesting extensions, such as
deviations from $w_\rm{DE} = 1$, nonminimal neutrino masses, spatial
curvature and deviations from general relativity using the growthindex
$\gamma$ parametrisation. These results are, in general, comparable to the most
precise presentday constraints on cosmological parameters, and show very good
agreement with the standard model. In particular, combining CMB,
$\omega(\theta)$ and SNIa, we find a value of $w_\rm{DE}$ consistent with $1$
to a precision better than 5\% when it is assumed to be constant in time, and
better than 6\% when we also allow for a spatiallycurved Universe.

We perform a tomographic baryon acoustic oscillations analysis using the
twopoint galaxy correlation function measured from the combined sample of BOSS
DR12, which covers the redshift range of $0.2<z<0.75$. Splitting the sample
into multiple overlapping redshift slices to extract the redshift information
of galaxy clustering, we obtain a measurement of $D_A(z)/r_d$ and $H(z)r_d$ at
nine effective redshifts with the full covariance matrix calibrated using
MultiDarkPatchy mock catalogues. Using the reconstructed galaxy catalogues, we
obtain the precision of $1.3\%2.2\%$ for $D_A(z)/r_d$ and $2.1\%6.0\%$ for
$H(z)r_d$. To quantify the gain from the tomographic information, we compare
the constraints on the cosmological parameters using our 9bin BAO
measurements, the consensus 3bin BAO and RSD measurements at three effective
redshifts in \citet{Alam2016}, and the nontomographic (1bin) BAO measurement
at a single effective redshift. Comparing the 9bin with 1bin constraint
result, it can improve the dark energy Figure of Merit by a factor of 1.24 for
the ChevallierPolarskiLinder parametrisation for equation of state parameter
$w_{\rm DE}$. The errors of $w_0$ and $w_a$ from 9bin constraints are slightly
improved when compared to the 3bin constraint result.

We study the possible effects of classical gravitational backgrounds on the
Higgs field through the modifications induced in the oneloop effective
potential and the vacuum expectation value of the energymomentum tensor. We
concentrate our study on the Higgs selfinteraction contribution in a perturbed
FRW metric. For weak and slowly varying gravitational fields, a complete set of
mode solutions for the KleinGordon equation is obtained to leading order in
the adiabatic approximation. Dimensional regularization has been used in the
integral evaluation and a detailed study of the integration of nonrational
functions in this formalism has been presented. As expected, the regularized
effective potential contains the same divergences as in flat spacetime, which
can be renormalized without the need of additional counterterms. We find that,
in contrast with other regularization methods, even though metric perturbations
affect the mode solutions, they do not contribute to the leading adiabatic
order of the potential. We also obtain explicit expressions of the complete
energymomentum tensor for general nonminimal coupling in terms of the
perturbed modes. The corresponding leading adiabatic contributions are also
obtained.

In current and future surveys, quasars play a key role. The new data will
extend our knowledge of the Universe as it will be used to better constrain the
cosmological model at redshift $z>1$ via baryon acoustic oscillation and
redshift space distortion measurements. Here, we present the first clustering
study of quasars observed by the extended Baryon Oscillation Spectroscopic
Survey. We measure the clustering of $\sim 70,000$ quasars located in the
redshift range $0.9<z<2.2$ that cover 1,168 deg$^2$. We model the clustering
and produce highfidelity quasar mock catalogues based on the BigMultiDark
Planck simulation. Thus, we use a modified (Sub)Halo Abundance Matching model
to account for the specificities of the halo population hosting quasars. We
find that quasars are hosted by halos with masses $\sim10^{12.7}M_\odot$ and
their bias evolves from 1.54 ($z=1.06$) to 3.15 ($z=1.98$). Using the current
eBOSS data, we cannot distinguish between models with different fractions of
satellites. The highfidelity mock lightcones, including properties of halos
hosting quasars, are made publicly available.

Upcoming and existing largescale surveys of galaxies require accurate
theoretical predictions of the dark matter clustering statistics for thousands
of mock galaxy catalogs. We demonstrate that this goal can be achieve with our
new Parallel ParticleMesh (PM) Nbody code (PPMGLAM) at a very low
computational cost. We run about 15,000 simulations with ~2 billion particles
that provide ~1% accuracy of the dark matter power spectra P(k) for
wavenumbers up to k~ 1h/Mpc. Using this large dataset we study the power
spectrum covariance matrix, the stepping stone for producing mock catalogs. In
contrast to many previous analytical and numerical results, we find that the
covariance matrix normalised to the power spectrum C(k,k')/P(k)P(k') has a
complex structure of nondiagonal components. It has an upturn at small k,
followed by a minimum at k=0.10.2h/Mpc. It also has a maximum at
k=0.50.6h/Mpc. The normalised covariance matrix strongly evolves with
redshift: C(k,k')~delta(t)^alpha P(k)P(k'), where delta is the linear growth
factor and alpha ~ 11.25, which indicates that the covariance matrix depends
on cosmological parameters. We also show that waves longer than 1Gpc have very
little impact on the power spectrum and covariance matrix. This significantly
reduces the computational costs and complexity of theoretical predictions:
relatively small volume ~ (1Gpc)^3 simulations capture the necessary properties
of dark matter clustering statistics. All the power spectra obtained from many
thousands of our simulations are publicly available.