
We compare analytical computations with numerical simulations for darkmatter
clustering, in general relativity and in the normal branch of DGP gravity
(nDGP). Our analytical frameword is the Effective Field Theory of LargeScale
Structure (EFTofLSS), which we use to compute the oneloop darkmatter power
spectrum, including the resummation of infrared bulk displacement effects. We
compare this to a set of 20 COLA simulations at redshifts $z = 0$, $z=0.5$, and
$z =1$, and fit the free parameter of the EFTofLSS, called the speed of sound,
in both $\Lambda$CDM and nDGP at each redshift. At oneloop at $z = 0$, the
reach of the EFTofLSS is $k_{\rm reach}\approx 0.14 \, h { \rm Mpc^{1}}$ for
both $\Lambda$CDM and nDGP. Along the way, we compare two different infrared
resummation schemes and two different treatments of the time dependence of the
perturbative expansion, concluding that they agree to approximately $1\%$ over
the scales of interest. Finally, we use the ratio of the COLA power spectra to
make a precision measurement of the difference between the speeds of sound in
$\Lambda$CDM and nDGP, and verify that this is proportional to the modification
of the linear coupling constant of the Poisson equation.

We develop a lattice Boltzmann (LB) model for immiscible twophase flow
simulations with central moments (CMs). This successfully combines a
threedimensional nonorthogonal CMbased LB scheme [A. De Rosis, Phys. Rev. E
95, 013310 (2017)] with our previous colorgradient LB model [S. Saito, Y. Abe,
and K. Koyama, Phys. Rev. E 96, 013317 (2017)]. Hydrodynamic meltjet breakup
simulations show that the proposed model is significantly more stable, even for
flow with extremely high Reynolds numbers, up to $O(10^6)$. This enables us to
investigate the phenomena expected under actual reactor conditions.

The recent simultaneous detection of gravitational waves and a gamma ray
burst from a neutron star merger significantly shrank the space of viable
scalartensor theories by demanding that the speed of gravity is equal to that
of light. The survived theories belong to the class of degenerate higher order
scalartensor theories. We study whether these theories are suitable as dark
energy candidates. We find scaling solutions in the matter dominated universe
that lead to de Sitter solutions at late times without the cosmological
constant, realising selfacceleration. We evaluate quasistatic perturbations
around selfaccelerating solutions and show that the stringent constraints
coming from astrophysical objects and gravitational waves can be satisfied,
leaving interesting possibilities to test these theories by cosmological
observations.

Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm
light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize
modifications of the growth of largescale structure in alternative theories of
gravity. We study the values these functions can take in Horndeski theories,
i.e. the class of scalartensor theories with second order equations of motion.
We restrict our attention to models that are in a broad agreement with tests of
gravity and the observed cosmic expansion history. In particular, we require
the speed of gravity to be equal to the speed of light today, as required by
the recent detection of gravitational waves and electromagnetic emission from a
binary neutron star merger. We examine the correlations between the values of
$\Sigma$ and $\mu$ analytically within the quasistatic approximation, and
numerically, by sampling the space of allowed solutions. We confirm that the
conjecture made in [Pogosian:2016pwr], that $(\Sigma1)(\mu 1) \ge 0$ in
viable Horndeski theories, holds very well. Along with that, we check the
validity of the quasistatic approximation within different corners of
Horndeski theory. Our results show that, even with the tight bound on the
present day speed of gravitational waves, there is room within Horndeski
theories for nontrivial signatures of modified gravity at the level of linear
perturbations.

The Vainshtein mechanism, present in many models of gravity, is very
effective at screening dark matter halos such that the fifth force is
negligible and general relativity is recovered within their Vainshtein radii.
Vainshtein screening is independent of halo mass and environment, in contrast
to e.g. chameleon screening, making it difficult to test. However, our previous
studies have found that the dark matter particles in filaments, walls, and
voids are not screened by the Vainshtein mechanism. We therefore investigate
whether cosmic voids, identified as local density minima using a watershed
technique, can be used to test models of gravity that exhibit Vainshtein
screening. We measure density, velocity, and screening profiles of stacked
voids in cosmological $N$body simulations using both dark matter particles and
dark matter halos as tracers of the density field. We find that the voids are
completely unscreened, and the tangential velocity and velocity dispersion
profiles of stacked voids show a clear deviation from $\Lambda$CDM at all
radii. Voids have the potential to provide a powerful test of gravity on
cosmological scales.

The almost simultaneous detection of gravitational waves and a short
gammaray burst from a neutron star merger has put a tight constraint on the
difference between the speed of gravity and light. In the fourdimensional
scalartensor theory with second order equations of motion, the Horndeski
theory, this translates into a significant reduction of the viable parameter
space of the theory. Recently, extensions of Horndeski theory, which are free
from Ostrogradsky ghosts despite the presence of higher order derivatives in
the equations of motion, have been identified and classified exploiting the
degeneracy criterium. In these new theories, the fifth force mediated by the
scalar field must be suppressed in order to evade the stringent Solar System
constraints. We study the Vainshtein mechanism in the most general degenerate
higher order scalartensor theory in which light and gravity propagate at the
same speed. We find that the Vainshtein mechanism generally works outside a
matter source but it is broken inside matter, similarly to beyond Horndeski
theories. This leaves interesting possibilities to test these theories that are
compatible with gravitational wave observations using astrophysical objects.

The effect of massive neutrinos on the growth of cold dark matter
perturbations acts as a scaledependent Newton's constant and leads to
scaledependent growth factors just as we often find in models of gravity
beyond General Relativity. We show how to compute growth factors for
$\Lambda$CDM and general modified gravity cosmologies combined with massive
neutrinos in Lagrangian perturbation theory for use in COLA and extensions
thereof. We implement this together with the gridbased massive neutrino method
of Brandbyge and Hannestad in $\texttt{MGPICOLA}$ and compare COLA simulations
to full $\it N$body simulations of $\Lambda$CDM and $f(R)$ gravity with
massive neutrinos. Our implementation is computationally cheap if the
underlying cosmology already has scaledependent growth factors and it is shown
to be able to produce results that match $\it N$body to percent level accuracy
for both the total and CDM matter powerspectra up to $k\lesssim 1 h/$Mpc.

Bayesian inference techniques are used to investigate situations where an
additional light scalar field is present during inflation and reheating. This
includes (but is not limited to) curvatontype models. We design a numerical
pipeline where $\simeq 200$ inflaton setups $\times\, 10$ reheating scenarios
$= 2000$ models are implemented and we present the results for a few
prototypical potentials. We find that singlefield models are remarkably robust
under the introduction of light scalar degrees of freedom. Models that are
ruled out at the singlefield level are not improved in general, because good
values of the spectral index and the tensortoscalar ratio can only be
obtained for very finetuned values of the extra field parameters and/or when
large nonGaussianities are produced. The only exception is quartic largefield
inflation, so that the best models after Planck are of two kinds: plateau
potentials, regardless of whether an extra field is added or not, and quartic
largefield inflation with an extra light scalar field, in some specific
reheating scenarios. Using Bayesian complexity, we also find that more
parameters are constrained for the models we study than for their singlefield
versions. This is because the added parameters not only contribute to the
reheating kinematics but also to the cosmological perturbations themselves, to
which the added field contributes. The interplay between these two effects lead
to a suppression of degeneracies that is responsible for having more
constrained parameters.

We show how standard Newtonian Nbody simulations can be interpreted in terms
of the weakfield limit of general relativity by employing the recently
developed Newtonian motion gauge. Our framework allows the inclusion of
radiation perturbations and the nonlinear evolution of matter. We show how to
construct the weakfield metric by combining Newtonian simulations with results
from EinsteinBoltzmann codes. We discuss observational effects on weak lensing
and ray tracing, identifying important relativistic corrections.

We extend our previous redshift space power spectrum code to the redshift
space correlation function. Here we focus on the Gaussian Streaming Model
(GSM). Again, the code accommodates a wide range of modified gravity and dark
energy models. For the nonlinear real space correlation function used in the
GSM, we use the Fourier transform of the RegPT 1loop matter power spectrum. We
compare predictions of the GSM for a Vainshtein screened and Chameleon screened
model as well as GR. These predictions are compared to the Fourier transform of
the Taruya, Nishimichi and Saito (TNS) redshift space power spectrum model
which is fit to NBody data. We find very good agreement between the Fourier
transform of the TNS model and the GSM predictions, with $\leq 6\%$ deviations
in the first two correlation function multipoles for all models for separations
in $50$Mpc$/h \leq s \leq 180$Mpc/$h$. Excellent agreement is found in the
differences between the modified gravity and GR multipole predictions for both
approaches to the redshift space correlation function, highlighting their
matched ability in picking up deviations from GR. We elucidate the timeliness
of such nonstandard templates at the dawn of stageIV surveys and discuss
necessary preparations and extensions needed for upcoming high quality data.

We compute the power spectrum at oneloop order in standard perturbation
theory for the matter density field to which a standard Lagrangian Baryonic
acoustic oscillation (BAO) reconstruction technique is applied. The BAO
reconstruction method corrects the bulk motion associated with the
gravitational evolution using the inverse Zel'dovich approximation (ZA) for the
smoothed density field. We find that the overall amplitude of oneloop
contributions in the matter power spectrum substantially decrease after
reconstruction. The reconstructed power spectrum thereby approaches the initial
linear spectrum when the smoothed density field is close enough to linear,
i.e., the smoothing scale $R_s$ larger than around 10$h^{1}$Mpc. On smaller
$R_s$,however, the deviation from the linear spectrum becomes significant on
large scales ($k\lt R_s^{1}$) due to the nonlinearity in the smoothed density
field, and the reconstruction is inaccurate. Compared with Nbody simulations,
we show that the reconstructed power spectrum at one loop order agrees with
simulations better than the unreconstructed power spectrum. We also calculate
the treelevel bispectrum in standard perturbation theory to investigate
nonGaussianity in the reconstructed matter density field. We show that the
amplitude of the bispectrum significantly decreases for small $k$ after
reconstruction and that the treelevel bispectrum agrees well with Nbody
results in the weakly nonlinear regime.

We present a general parallelized and easytouse code to perform numerical
simulations of structure formation using the COLA (COmoving Lagrangian
Acceleration) method for cosmological models that exhibit scaledependent
growth at the level of first and second order Lagrangian perturbation theory.
For modified gravity theories we also include screening using a fast
approximate method that covers all the main examples of screening mechanisms in
the literature. We test the code by comparing it to full simulations of two
popular modified gravity models, namely $f(R)$ gravity and nDGP, and find good
agreement in the modified gravity boostfactors relative to $\Lambda$CDM even
when using a fairly small number of COLA time steps.

A flat FriedmanRobersonWalker universe dominated by a cosmological constant
($\Lambda$) and cold dark matter (CDM) has been the working model preferred by
cosmologists since the discovery of cosmic acceleration. However, tensions of
various degrees of significance are known to be present among existing datasets
within the $\Lambda$CDM framework. In particular, the Lyman$\alpha$ forest
measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation
Spectroscopic Survey (BOSS) prefers a smaller value of the matter density
fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave
background (CMB). Also, the recently measured value of the Hubble constant,
$H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{1} \ {\rm Mpc}^{1}$, is $3.4\sigma$
higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{1} \ {\rm Mpc}^{1}$ inferred
from the Planck CMB data. In this work, we investigate if these tensions can be
interpreted as evidence for a nonconstant dynamical dark energy (DE). Using
the KullbackLeibler (KL) divergence to quantify the tension between datasets,
we find that the tensions are relieved by an evolving DE, with the dynamical DE
model preferred at a $3.5\sigma$ significance level based on the improvement in
the fit alone. While, at present, the Bayesian evidence for the dynamical DE is
insufficient to favour it over $\Lambda$CDM, we show that, if the current best
fit DE happened to be the true model, it would be decisively detected by the
upcoming DESI survey.

We present a novel cosmological solution in the framework of extended
quasidilaton theory which underwent scrutiny recently. We only consider terms
that do not generate the BoulwareDeser degree of freedom, hence the
"ghostfree" quasidilaton theory, and show three new branches of cosmological
evolution therein. One of the solutions passes the perturbative stability
tests. This new solution exhibits a late time selfacceleration and all
graviton polarizations acquire masses that converge to a constant in the
asymptotic future. Moreover, all modes propagate at the speed of light. We
propose that this solution can be used as a benchmark model for future
phenomenological studies.

A threedimensional colorfluid lattice Boltzmann model for immiscible
twophase flows is developed in the framework of a threedimensional
27velocity (D3Q27) lattice. The collision operator comprises the D3Q27
versions of three suboperators: a multiplerelaxationtime (MRT) collision
operator, a generalized LiuValocchiKang perturbation operator, and a
LatvaKokkoRothman recoloring operator. A D3Q27 version of an enhanced
equilibrium distribution function is also incorporated into this model to
improve the Galilean invariance. Three types of numerical tests, namely, a
static droplet, an oscillating droplet, and the RayleighTaylor instability,
show a good agreement with analytical solutions and numerical simulations.
Following these numerical tests, this model is applied to liquidjetbreakup
simulations. The simulation conditions are matched to the conditions of the
previous experiments. In this case, numerical stability is maintained
throughout the simulation, although the kinematic viscosity for the continuous
phase is set as low as $1.8\times10^{4}$, in which case the corresponding
Reynolds number is $3.4\times10^{3}$; the developed lattice Boltzmann model
based on the D3Q27 lattice enables us to perform the simulation with parameters
directly matched to the experiments. The jet's liquid column transitions from
an asymmetrical to an axisymmetrical shape, and entrainment occurs from the
side of the jet. The measured time history of the jet's leadingedge position
shows a good agreement with the experiments. Finally, the reproducibility of
the regime map for liquidliquid systems is assessed. The present lattice
Boltzmann simulations well reproduce the characteristics of predicted regimes,
including varicose breakup, sinuous breakup, and atomization.

Initial conditions for (Newtonian) cosmological Nbody simulations are
usually set by rescaling the presentday power spectrum obtained from linear
(relativistic) Boltzmann codes to the desired initial redshift of the
simulation. This backscaling method can account for the effect of
inhomogeneous residual thermal radiation at early times, which is absent in the
Newtonian simulations. We analyse this procedure from a fully relativistic
perspective, employing the recentlyproposed Newtonian motion gauge framework.
We find that Nbody simulations for LambdaCDM cosmology starting from
backscaled initial conditions can be selfconsistently embedded in a
relativistic spacetime with firstorder metric potentials calculated using a
linear Boltzmann code. This spacetime coincides with a simple "Nbody gauge"
for z<50 for all observable modes. Care must be taken, however, when simulating
nonstandard cosmologies. As an example, we analyse the backscaling method in
a cosmology with decaying dark matter, and show that metric perturbations
become large at early times in the backscaling approach, indicating a
breakdown of the perturbative description. We suggest a suitable "forwards
approach" for such cases.

We elucidate the importance of the consistent treatment of gravitymodel
specific nonlinearities when estimating the growth of cosmological structures
from redshift space distortions (RSD). Within the context of standard
perturbation theory (SPT), we compare the predictions of two theoretical
templates with redshift space data from COLA (COmoving Lagrangian Acceleration)
simulations in the normal branch of DGP gravity (nDGP) and General Relativity
(GR). Using COLA for these comparisons is validated using a suite of full
Nbody simulations for the same theories. The two theoretical templates
correspond to the standard general relativistic perturbation equations and
those same equations modelled within nDGP. Gravitational clustering nonlinear
effects are accounted for by modelling the power spectrum up to one loop order
and redshift space clustering anisotropy is modelled using the Taruya,
Nishimichi and Saito (TNS) RSD model. Using this approach, we attempt to
recover the simulation's fiducial logarithmic growth parameter $f$. By
assigning the simulation data with errors representing an idealised survey with
a volume of $10\mbox{Gpc}^3/h^3$, we find the GR template is unable to recover
fiducial $f$ to within 1$\sigma$ at $z=1$ when we match the data up to $k_{\rm
max}=0.195h$/Mpc. On the other hand, the DGP template recovers the fiducial
value within $1\sigma$. Further, we conduct the same analysis for sets of mock
data generated for generalised models of modified gravity using SPT, where
again we analyse the GR template's ability to recover the fiducial value. We
find that for models with enhanced gravitational nonlinearity, the theoretical
bias of the GR template becomes significant for stage IV surveys. Thus, we show
that for the future large data volume galaxy surveys, the selfconsistent
modelling of nonGR gravity scenarios will be crucial in constraining theory
parameters.

We describe and demonstrate the potential of a new and very efficient method
for simulating certain classes of modified gravity theories, such as the widely
studied $f(R)$ gravity models. High resolution simulations for such models are
currently very slow due to the highly nonlinear partial differential equation
that needs to be solved exactly to predict the modified gravitational force.
This nonlinearity is partly inherent, but is also exacerbated by the specific
numerical algorithm used, which employs a variable redefinition to prevent
numerical instabilities. The standard NewtonGaussSeidel iterative method used
to tackle this problem has a poor convergence rate. Our new method not only
avoids this, but also allows the discretised equation to be written in a form
that is analytically solvable. We show that this new method greatly improves
the performance and efficiency of $f(R)$ simulations. For example, a test
simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is
now 5 times faster than before, while a Millenniumresolution simulation for
$f(R)$ gravity is estimated to be more than 20 times faster than with the old
method. Our new implementation will be particularly useful for running very
high resolution, largesized simulations which, to date, are only possible for
the standard model, and also makes it feasible to run large numbers of lower
resolution simulations for covariance analyses. We hope that the method will
bring us to a new era for precision cosmological tests of gravity.

Theories of gravity in the beyond Horndeski class recover the predictions of
general relativity in the solar system whilst admitting novel cosmologies,
including latetime de Sitter solutions in the absence of a cosmological
constant. Deviations from Newton's law are predicted inside astrophysical
bodies, which allow for falsifiable, smokinggun tests of the theory. In this
work we study the pulsations of stars by deriving and solving the wave equation
governing linear adiabatic oscillations to find the modified period of
pulsation. Using both semianalytic and numerical models, we perform a
preliminary survey of the stellar zoo in an attempt to identify the best
candidate objects for testing the theory. Brown dwarfs and Cepheid stars are
found to be particularly sensitive objects and we discuss the possibility of
using both to test the theory.

Theories of gravity in the beyond Horndeski class encompass a wide range of
scalartensor theories that will be tested on cosmological scales over the
coming decade. In this work, we investigate the possibility of testing them in
the strongfield regime by looking at the properties of compact
objectsneutron, hyperon, and quark starsembedded in an asymptotically de
Sitter spacetime, for a specific subclass of theories. We extend previous
works to include slow rotation and find a relation between the dimensionless
moment of intertia, ($\bar{I}=Ic^2/G_{\rm N} M^3$), and the compactness,
$\cal{C}=G_{\rm N} M/Rc^2$ (an $\bar{I}$$\cal{C}$ relation), independent of
the equation of state, that is reminiscent of but distinct from the general
relativity prediction. Several of our equations of state contain hyperons and
free quarks, allowing us to revisit the hyperon puzzle. We find that the
maximum mass of hyperon stars can be larger than $2M_\odot$ for small values of
the beyond Horndeski parameter, thus providing a resolution of the hyperon
puzzle based on modified gravity. Moreover, stable quark stars exist when
hyperonic stars are unstable, which means that the phase transition from
hyperon to quark stars is predicted just as in general relativity, albeit with
larger quark star masses. Two important and potentially observable consequences
of some of the theories we consider are the existence of neutron stars in a
range of masses significantly higher than in GR, and $\bar{I}$$\mathcal{C}$
relations that differ from their GR counterparts. In the former case, we find
objects that, if observed, could not be accounted for in GR because they
violate the usual GR causality condition. We end by discussing several
difficult technical issues that remain to be addressed in order to reach more
realistic predictions that may be tested using gravitational wave searches or
neutron star observations.

We derive the first systematic observational constraints on reheating in
models of inflation where an additional light scalar field contributes to
primordial density perturbations and affects the expansion history during
reheating. This encompasses the original curvaton model but also covers a
larger class of scenarios. We find that, compared to the singlefield case,
lower values of the energy density at the end of inflation and of the reheating
temperature are preferred when an additional scalar field is introduced. For
instance, if inflation is driven by a quartic potential, which is one of the
most favoured models when a light scalar field is added, the upper bound
$T_{\mathrm{reh}}<5\times 10^{4}\,\mathrm{GeV}$ on the reheating temperature
$T_{\mathrm{reh}}$ is derived, and the implications of this value on
postinflationary physics are discussed. The information gained about reheating
is also quantified and it is found that it remains modest in plateau inflation
(though still larger than in the singlefield version of the model) but can
become substantial in quartic inflation. The role played by the vev of the
additional scalar field at the end of inflation is highlighted, and opens
interesting possibilities for exploring stochastic inflation effects that could
determine its distribution.

We present all scalartensor Lagrangians that are cubic in second derivatives
of a scalar field, and that are degenerate, hence avoiding Ostrogradsky
instabilities. Thanks to the existence of constraints, they propagate no more
than three degrees of freedom, despite having higher order equations of motion.
We also determine the viable combinations of previously identified quadratic
degenerate Lagrangians and the newly established cubic ones. Finally, we study
whether the new theories are connected to known scalartensor theories such as
Horndeski and beyond Horndeski, through conformal or disformal transformations.

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.

The extended quasidilaton theory is one of the simplest Lorentzinvariant
massive gravity theories which can accommodate a stable selfaccelerating
vacuum solution. In this paper we revisit this theory and study the effect of
matter fields. For a matter sector that couples minimally to the physical
metric, we find hints of a Jeans type instability in the IR. In the analogue
kessence field setup, this instability manifests itself as an IR ghost for
the scalar field perturbation, but this can be interpreted as a classical
instability that becomes relevant below some momentum scale in terms of matter
density perturbations. We also consider the effect of the background evolution
influenced by matter on the stability of the gravity sector perturbations. In
particular, we address the previous claims of ghost instability in the IR
around the late time attractor. We show that, although the matterinduced
modification of the evolution potentially brings tension to the stability
conditions, one goes beyond the regime of validity of the effective theory well
before the solutions become unstable. We also draw attention to the fact that
the IR stability conditions are also enforced by the existence requirements of
consistent background solutions.