
THESEUS is a space mission concept aimed at exploiting GammaRay Bursts for
investigating the early Universe and at providing a substantial advancement of
multimessenger and timedomain astrophysics. These goals will be achieved
through a unique combination of instruments allowing GRB and Xray transient
detection over a broad field of view (more than 1sr) with 0.51 arcmin
localization, an energy band extending from several MeV down to 0.3 keV and
high sensitivity to transient sources in the soft Xray domain, as well as
onboard prompt (few minutes) followup with a 0.7 m class IR telescope with
both imaging and spectroscopic capabilities. THESEUS will be perfectly suited
for addressing the main open issues in cosmology such as, e.g., star formation
rate and metallicity evolution of the interstellar and intragalactic medium
up to redshift $\sim$10, signatures of Pop III stars, sources and physics of
reionization, and the faint end of the galaxy luminosity function. In
addition, it will provide unprecedented capability to monitor the Xray
variable sky, thus detecting, localizing, and identifying the electromagnetic
counterparts to sources of gravitational radiation, which may be routinely
detected in the late '20s / early '30s by next generation facilities like
aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide
powerful synergies with the next generation of multiwavelength observatories
(e.g., LSST, ELT, SKA, CTA, ATHENA).

In the framework of phantom quintessence cosmology, we use the Noether
Symmetry Approach to obtain general exact solutions for the cosmological
equations. This result is achieved by the quintessential (phantom) potential
determined by the existence of the symmetry itself. A comparison between the
theoretical model and observations is worked out. In particular, we use type Ia
supernovae and large scale structure parameters determined from the 2degree
Field Galaxy Redshift Survey (2dFGRS)and from the Wide part of the VIMOSVLT
Deep Survey (VVDS). It turns out that the model is compatible with the
presently available observational data. Moreover we extend the approach to
include radiation. We show that it is compatible with data derived from
recombination and it seems that quintessence do not affect nucleosynthesis
results.

The weak field limit for a pointlike source of a $f(R) \propto
R^{3/2}$gravity model is studied. We aim to show the viability of such a model
as a valid alternative to GR + dark matter at Galactic and local scales.
Without considering dark matter, within the weak field approximation, we find
general exact solutions for gravity with standard matter, and apply them to
some astrophysical scales, recovering the consistency of the same
$f(R)$gravity model with cosmological results.}{In particular, we show that it
is possible to obtain flat rotation curves for galaxies, [and consistency with]
Solar System tests, as in the socalled "Chameleon Approach". In fact, the
peripheral velocity $ v_\infty $ is shown to be expressed as $ v_\infty =
\lambda \sqrt{M}$, so that the TullyFisher relation is recovered. The results
point out the possibility of achieving alternative theories of gravity in which
exotic ingredients like dark matter and dark energy are not necessary, while
their coarsegrained astrophysical and cosmological effects can be related to a
geometric origin.

We consider cosmological models in scalar tensor theories of gravity that
describe an accelerating universe, and we study a family of inverse power law
potentials, for which exact solutions of the Einstein equations are known. We
also compare theoretical predictions of our models with observations. For this
we use the following data: the publicly available catalogs of type Ia
supernovae and high redshift Gamma Ray Bursts, the parameters of large scale
structure determined by the 2degree Field Galaxy Redshift Survey (2dFGRS), and
measurements of cosmological distances based on the SunyaevZel'dovich effect,
among others.

Scalartensor quintessence models can be constrained by identifying suitable
cosmic clocks which allow to select confidence regions for cosmological
parameters. In particular, we constrain the characterizing parameters of
nonminimally coupled scalartensor cosmological models which admit exact
solutions of the Einstein field equations. Lookback time to galaxy clusters at
low intermediate, and high redshifts is considered. The high redshift
timescale problem is also discussed in order to select other cosmic clocks
such as quasars.

Recent observations of the line of sight velocity profile of elliptical
galaxies have furnished controversial results with some works favouring the
presence of a large amount of dark matter in the outer regions and others
arguing in favour of no dark matter at all. In order to shed new light on this
controversy, we propose here a new phenomenological description of the total
mass profile of galaxies. Under the hypothesis of spherical symmetry, we assume
a double powerlaw expression for the global M/L ratio Upsilon(r)=
Upsilon_0(r/r_0) ^{alpha}(1+r/r_0)^{beta}. In particular, Upsilon propto
r^{alpha} for r/r_0<<1 so that alpha=0 mimics a constant M/L ratio in the inner
regions, while, for (r/r_0>>1), Upsilon propto r^{alpha+beta} thus showing that
models with alpha+beta=0 have an asymptotically constant M/L ratio. A wide
range of possibilities is obtained by varying the slope parameters in the range
we determine on the basis of physical considerations. Choosing a general
expression for the luminosity density profile j(r), we work out an effective
galaxy model that accounts for all the phenomenology observed in real
elliptical galaxies. We derive the main dynamics and lensing properties of such
an effective model. We analyze a general class of models, able to take into
account different dynamical trends. We are able to obtain analytical
expressions for the main dynamical and lensing quantities. We show that
constraining the values of alpha+beta makes it possible to analyze the problem
of the dark matter in elliptical galaxies. Indeed, positive values of
alpha+beta would be a strong evidence for dark matter. Finally we indicate
possible future approaches in order to face the observational data, in
particular using velocity dispersion profiles and lensed quasar events.

We consider scalar tensor theories of gravity assuming that the scalar field
is non minimally coupled with gravity. We use this theory to study evolution of
a flat homogeneous and isotropic universe. In this case the dynamical equations
can be derived form a point like Lagrangian. We study the general properties of
dynamics of this system and show that for a wide range of initial conditions
such models lead in a natural way to an accelerated phase of expansion of the
universe. Assuming that the point like Lagrangian admits a Noether symmetry we
are able to explicitly solve the dynamical equations. We study one particular
model and show that its predictions are compatible with observational data,
namely the publicly available data on type Ia supernovae, the parameters of
large scale structure determined by the 2degree Field Galaxy Redshift Survey
(2dFGRS), the measurements of cosmological distances with the
SunyaevZel'dovich effect and the rate of growth of density perturbations}{It
turns out that this model have a very interesting feature of producing in a
natural way an epoch of accelerated expansion. With an appropriate choice of
parameters our model is fully compatible with several observed characteristics
of the universe

It has been recently shown that, under some general conditions, it is always
possible to find a fourth order gravity theory capable of reproducing the same
dynamics of a given dark energy model. Here, we discuss this approach for a
dark energy model with a scalar field evolving under the action of an
exponential potential. In absence of matter, such a potential can be recovered
from a fourth order theory via a conformal transformation. Including the matter
term, the function f(R) entering the generalized gravity Lagrangian can be
reconstructed according to the dark energy model.

We present a new family of spherically symmetric models for the luminous
components of elliptical and spiral galaxies and their dark matter haloes. Our
starting point is a general expression for the logarithmic slope $\alpha(r) =
d\log{\rho}/d\log{r}$ from which most of the cuspy models yet available in
literature may be derived. We then dedicate our attention to a particular set
of models whose logarithmic slope is a power law function of the radius $r$
investigating in detail their dynamics assuming isotropy in the velocity space.
While the basic properties (such as the density profile and the gravitational
potential) may be expressed analytically, both the distribution function and
the observable quantities (surface brightness and line of sight velocity
dispersion) have to be evaluated numerically. We also consider the extension to
anisotropic models trying two different parameterization. Since the model
recently proposed by Navarro et al. (2004) as the best fit to their sample of
numerically simulated haloes belongs to the family we present here, analytical
approximations are given for the most useful quantities.

We use some of the recently released observational data to test the viability
of two classes of minimally coupled scalar field models of quintessence with
exponential potentials for which exact solutions of the Einstein equations are
known. These models are very sturdy, depending on only one parameter  the
Hubble constant. To compare predictions of our models with observations we
concentrate on the following data: the power spectrum of the CMBR anisotropy as
measured by WMAP, the publicly available data on type Ia supernovae, and the
parameters of large scale structure determined by the 2degree Field Galaxy
Redshift Survey (2dFGRS). We use the WMAP data on the age of the universe and
the Hubble constant to fix the free parameters in our models. We then show that
the predictions of our models are consistent with the observed positions and
relative heights of the first 3 peaks in the CMB power spectrum, with the
energy density of dark energy as deduced from observations of distant type Ia
supernovae, and with parameters of the large scale structure as determined by
2dFGRS, in particular with the average density of dark matter. Our models are
also consistent with the results of the Sloan Digital Sky Survey (SDSS).
Moreover, we investigate the evolution of matter density perturbations in our
quintessential models, solve exactly the evolution equation for the density
perturbations, and obtain an analytical expression for the growth index $f$. We
verify that the approximate relation f ~ Omega_M^(alpha) also holds in our
models.

The structure of the dark matter and the thermodynamical status of the hot
gas in galaxy clusters is an interesting and widely discussed topic in modern
astrophysics. Recently, Rasia et al. (2004) have proposed a new dynamical model
for the mass density profile of clusters of galaxies as a result of a set of
high resolution hydrodynamical simulations of structure formation. We
investigate the lensing properties of this model evaluating the deflection
angle, the lensing potential and the amplification of the images. We reserve
particular attention to the structure and position of the critical curves in
order to see whether this model is able to produce radial and tangential arcs.
To this aim, we also investigate the effect of taking into account the
brightest cluster galaxy in the lensing potential and the deviations from
spherical symmetry mimicked by an external shear. We also analyze the
implication of the gas density and temperature profiles of the Rasia et
al.(2004) model on the properties of the X  ray emission and the
comptonization parameter that determines the CMBR temperature decrement due to
the Sunyaev  Zel'dovich effect.

We present a numerical method to estimate the lensing parameters and the
Hubble constant H_0 from quadruply imaged gravitational lens systems. The lens
galaxy is modeled using both separable deflection potentials and constant
masstolight ratio profiles, while possible external perturbations have been
taken into account introducing an external shear. The model parameters are
recovered inverting the lens and the time delay ratio equations and imposing a
set of physically motivated selection criteria. We investigate correlations
among the model parameters and the Hubble constant. Finally, we apply the codes
to the real lensed quasars PG 1115+080 and RX J0911+0551, and combine the
results from these two systems to get H_0 = 56 +/ 23 km/(s Mpc). In addition,
we are able to fit to the single systems a general elliptical potential with a
non fixed angular part, and then we model the two lens systems with the same
potential and a shared H_0: in this last case we estimate H_0=49_(11)^(+6)
km/(s Mpc).

We show that a general, exact cosmological solution, where dynamics of scalar
field is assigned by an exponential potential, fulfils all the issues of dark
energy approach, both from a theoretical point of view and in comparison with
available observational data. Moreover, tracking conditions are discussed, with
a new treatment of the well known condition $\Gamma>1$. We prove that the
currently used expression for $\Gamma$ is wrong.

We present a dark energy model with a double exponential potential for a
minimally coupled scalar field, which allows general exact integration of the
cosmological equations. The solution can perfectly emulate a $\Lambda$term
model at lowmedium redshift, exhibits tracking behavior, and does not show
eternal acceleration.

We discuss the general and approximate angular diameter distance in the
FriedmanRobertsonWalker cosmological models with nonzero cosmological
constant. We modify the equation for the angular diameter distance by taking
into account the fact that locally the distribution of matter is non
homogeneous. We present exact solutions of this equation in a few special
cases. We propose an approximate analytic solution of this equation which is
simple enough and sufficiently accurate to be useful in practical applications.

Noether symmetry for higher order gravity theory has been explored, with the
introduction of an auxiliary variable which gives the only correct quantum
desccription of the theory, as shown in a series of earlier papers. The
application of Noether theorem in higher order theory of gravity turned out to
be a powerful tool to find the solution of the field equations. A few such
physically reasonable solutions like power law inflation are presented.

We explore the conditions for the existence of Noether symmetries in the
dynamics of FRW metric, non minimally coupled with a scalar field, in the most
general situation, and with nonzero spatial curvature. When such symmetries are
present we find general exact solution for the Einstein equations. We also show
that non Noether symmetries can be found.
Finally,we present an extension of the procedure to the Kantowski Sachs
metric which is particularly interesting in the case of degenerate Lagrangian.

Taking into account a torsion field gives rise to a negative pressure
contribution in cosmological dynamics and then to an accelerated behaviour of
Hubble fluid. The presence of torsion has the same effect of a Lambda  term.
We obtain a general exact solution which well fits the data coming from high 
redshift supernovae and Sunyaev  Zeldovich/X  ray method for the
determination of the cosmological parameters. On the other hand, it is possible
to obtain observational constraints on the amount of torsion density. A dust
dominated Friedmann behaviour is recovered as soon as torsion effects are not
relevant.

We discuss the amplification dispersion in the observed luminosity of
standard candles, like supernovae (SNe) of type Ia, induced by gravitational
lensing in a Universe with dark energy (quintessence). We derive the main
features of the magnification probability distribution function (pdf) of SNe in
the framework of on average FriedmannLemaitreRobertsonWalker (FLRW) models
for both lensing by largescale structures and compact objects. The
magnification pdf is strongly dependent on the equation of state, $w_Q$, of the
quintessence. The dispersion increases with the redshift of the source and is
maximum for dark energy with very large negative pressure; the effects of
gravitational lensing on the magnification pdf, i.e. the mode biased towards
deamplified values and the long tail towards large magnifications, are reduced
for both microscopic DM and quintessence with an intermediate $w_Q$. Different
equations of state of the dark energy can deeply change the dispersion in
amplification for the projected observed samples of SNe Ia by future spaceborn
missions. The "noise" in the Hubble diagram due to gravitational lensing
strongly affects the determination of the cosmological parameters from SNe
data. The errors on the pressureless matter density parameter, $\Omega_M$, and
on $w_Q$ are maximum for quintessence with not very negative pressure. The
effect of the gravitational lensing is of the same order of the other
systematics affecting observations of SNe Ia. Due to the lensing by largescale
structures, in a flat Universe with $\Omega_M =0.4$, at $z=1$ a cosmological
constant ($w_Q=1$) can be interpreted as dark energy with $w_Q <0.84$ (at
2$\sigma$ confidence limit).

We present a new method to estimate the Hubble constant H_0 from the measured
time delays in quadruply imaged gravitational lens systems. We show how it is
possible to get an estimate of H_0 without the need to completely reconstruct
the lensing potential thus avoiding any a priori hypotheses on the expression
of the galaxy lens model. Our method only needs to assume that the lens
potential may be expressed as r^{\alpha} F(\theta), whatever the shape function
F(\theta) is, and it is thus able to fully explore the degeneracy in the mass
models taking also into account the presence of an external shear. We test the
method on simulated cases and show that it does work well in recovering the
correct value of the slope \alpha of the radial profile and of the Hubble
constant H_0. Then, we apply the same method to the real quadruple lenses
PG1115+080 and B1422+231 obtaining H_0 = 58_{15}^{+17} km/s/Mpc (68% CL).

We develop a semi  analytical method to reconstruct the lensing potential in
quadruply imaged gravitational lens systems. Assuming that the potential
belongs to a broad class of boxy non  elliptical models, we show how it is
possible to write down a system of equations which can be numerically solved to
recover the potential parameters directly from image positions and using
physical constraints. We also describe a code developed to search for solutions
of the system previously found and test it on simulated cases. Finally, we
apply the method to the quadruple lens PG1115+080 which allows us to get also
an estimate of the Hubble constant H_0 from the measured time delay as H_0 =
56_{11}^{+17} km/s/Mpc.

We present the first results of the analysis of data collected during the
199899 observational campaign at the 1.3 meter McGrawHill Telescope, towards
the Andromeda galaxy (M31), aimed to the detection of gravitational
microlensing effects as a probe of the presence of dark matter in our and in
M31 halo. The analysis is performed using the pixel lensing technique, which
consists in the study of flux variations of unresolved sources and has been
proposed and implemented by the AGAPE collaboration. We carry out a shape
analysis by demanding that the detected flux variations be achromatic and
compatible with a Paczynski light curve. We apply the DurbinWatson hypothesis
test to the residuals. Furthermore, we consider the background of variables
sources. Finally five candidate microlensing events emerge from our selection.
Comparing with the predictions of a Monte Carlo simulation, assuming a standard
spherical model for the M31 and Galactic haloes, and typical values for the
MACHO mass, we find that our events are only marginally consistent with the
distribution of observable parameters predicted by the simulation.

We investigate the properties of cosmological distances in locally
inhomogeneous universes with pressureless matter and dark energy
(quintessence), with constant equation of state. We give exact solutions for
angular diameter distances in theempty beam approximation. In this hypothesis,
the distanceredshift equation is derived fron the multiple lensplane theory.
The case of a flat universe is considered with particular attention. We show
how this general scheme makes distances degenerate with respect to w_X and the
smoothness parameters, alpha, accounting for the homogeneously distributed
fraction of energy of the icomponents. We analyse how this degeneracy
influences the critical redshift where the angular diameter distance takes its
maximum, and put in evidence future prospects for measuring the smoothness
parameter of the pressureless matter, alpha_M.