Non-gravitational interaction between dark matter and dark energy has been considered in a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) universe. The interaction rate is assumed to be linear in the energy densities of dark matter and dark energy and it is also proportional to the Hubble rate of the FLRW universe. This kind of interaction model leads to an autonomous linear dynamical system, and depending on the coupling parameters, could be solved analytically by calculating the exponential of the matrix, defining the system. We show here that such interaction rate has a very deep connection with some well known cosmological theories. We then investigate the theoretical bounds on the coupling parameters of the interaction rate in order that the energy densities of the dark sector remain positive throughout the evolution of the universe and asymptotically converge to zero at very late times. Our analyses also point out that such linear interacting model may encounter with finite time future singularities depending on the coupling parameters as well as the dark energy state parameter.

We report follow-up observations of transiting exoplanets that have either large uncertainties (>10 minutes) in their transit times or have not been observed for over three years. A fully robotic ground-based telescope network, observations from citizen astronomers and data from TESS have been used to study eight planets, refining their ephemeris and orbital data. Such follow-up observations are key for ensuring accurate transit times for upcoming ground and space-based telescopes which may seek to characterise the atmospheres of these planets. We find deviations from the expected transit time for all planets, with transits occurring outside the 1 sigma uncertainties for seven planets. Using the newly acquired observations, we subsequently refine their periods and reduce the current predicted ephemeris uncertainties to 0.28 - 4.01 minutes. A significant portion of this work has been completed by students at two high schools in London as part of the Original Research By Young Twinkle Students (ORBYTS) programme.

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

In this work we construct a bottom-up reconstruction technique for Loop Quantum Cosmology scalar-tensor theories, from the observational indices. Particularly, the reconstruction technique is based on fixing the functional form of the scalar-to-tensor ratio as a function of the $e$-foldings number. The aim of the technique is to realize viable inflationary scenarios, and the only assumption that must hold true in order for the reconstruction technique to work is that the dynamical evolution of the scalar field obeys the slow-roll conditions. We shall use two functional forms for the scalar-to-tensor ratio, one of which corresponds to a popular inflationary class of models, the $\alpha$-attractors. For the latter, we shall calculate the leading order behavior of the spectral index and we shall demonstrate that the resulting inflationary theory is viable and compatible with the latest Planck and BICEP2/Keck-Array data. In addition, we shall find the classical limit of the theory, and as we demonstrate, the Loop Quantum Cosmology corrected theory and the classical theory are identical at leading order in the perturbative expansion quantified by the parameter $\rho_c$, which is the critical density of the quantum theory. Finally, by using the formalism of slow-roll scalar-tensor Loop Quantum Cosmology, we shall investigate how several inflationary potentials can be realized by the quantum theory, and we shall calculate directly the slow-roll indices and the corresponding observational indices. In addition, the $f(R)$ gravity frame picture is presented.

The PLANCK collaboration has determined values for the spectral parameters of the CMB radiation, namely the spectral index $n_s$, its running $\alpha_s$, the running of the running $\beta_s$, using a growing body of measurements of CMB anisotropies by the Planck satellite and other missions. These values do not follow the hierarchy of sizes predicted by single field, slow roll inflationary theory, and are thus difficult to fit for such inflation models. In this work we present first a study of 49 single field, slow roll inflationary potentials in which we assess the likelyhood of these models fitting the spectral parameters to their currently most accurate determination given by the PLANCK collaboration. We check numerically with a MATLAB program the spectral parameters that each model can yield for a very broad, comprehensive list of possible parameter and field values. The comparison of spectral parameter values supported by the models with their determinations by the PLANCK collaboration leads to the conclusion that the data provided by PLANCK2015 TT+lowP and PLANCK2015 TT,TE,EE+lowP taking into account the running of the running disfavours 40 of the 49 models with confidence level at least 92.8\%. Next, we discuss the reliability of the current computations of these spectral parameters. We identify a bias in the method of determination of the spectral parameters by least residue parameter fitting (using MCMC or any other scheme) currently used to reconstruct the power spectrum of scalar perturbations. This bias can explain the observed contradiction between theory and observations. Its removal is computationally costly, but necessary in order to compare the forecasts of single field, slow roll theories with observations.

In a flat Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) geometry, we consider the expansion of the universe powered by the gravitationally induced `adiabatic' matter creation. To demonstrate how matter creation works well with the expanding universe, we have considered a general creation rate and analyzed this rate in the framework of dynamical analysis. The dynamical analysis hints the presence of a non-singular universe (without the big bang singularity) with two successive accelerated phases, one at the very early phase of the universe (i.e. inflation), and the other one describes the current accelerating universe, where this early, late accelerated phases are associated with an unstable fixed point (i.e. repeller) and a stable fixed (attractor) points, respectively. We have described this phenomena by analytic solutions of the Hubble function and the scale factor of the FLRW universe. Using Jacobi Last multiplier method, we have found a Lagrangian for this matter creation rate describing this scenario of the universe. To match with our early physics results, we introduce an equivalent dynamics driven by a single scalar field and discussed the associated observable parameters compared them with the latest PLANCK data sets. Finally, introducing the teleparallel modified gravity, we have established an equivalent gravitational theory in the framework of matter creation.

A model of a universe without big bang singularity is presented, which displaysanearly inflationary period ending just before a phase transition to a deflationary epoch. The model produces enough heavy particles so as to reheat the universe at temperatures in the MeV regime. After the reheating, it smoothly matches the standard $\Lambda$CDM scenario.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

It is shown that teleparallel $F({\mathcal T})$ theories of gravity combined with Loop Quantum Cosmology support a Matter Bounce Scenario which is an alternative to the inflation scenario in the Big Bang paradigm. It is checked thatthese bouncing models provide theoretical data that fits well with the current observational data, allowing the viability of the Matter Bounce Scenario.

The CMB map provided by the Planck project constrains the value of the ratio of tensor-to-scalar perturbations, namely $r$, to be smaller than $0.11$ (95% CL). This bound rules out the simplest models of inflation. However, recent data from BICEP2 is in strong tension with this constrain, as it finds a value $r=0.20^{+0.07}_{-0.05}$ with $r=0$ disfavored at $7.0 \sigma$, which allows these simplest inflationary models to survive. The remarkable fact is that, even though the BICEP2 experiment was conceived to search for evidence of inflation, its experimental data matches correctly theoretical results coming from the matter bounce scenario (the alternative model to the inflationary paradigm). More precisely, most bouncing cosmologies do not pass Planck's constrains due to the smallness of the value of the tensor/scalar ratio $r\leq 0.11$, but with new BICEP2 data some of them fit well with experimental data. This is the case with the matter bounce scenario in the teleparallel version of Loop Quantum Cosmology.

Several dark energy models on the brane are investigated. They are compared with corresponding theories in the frame of 4d Friedmann-Robertson-Walker cosmology. To constrain the parameters of the models considered, recent observational data, including SNIa apparent magnitude measurements, baryon acoustic oscillation results, Hubble parameter evolution data and matter density perturbations are used. Explicit formulas of the so-called {\it state-finder} parameters in teleparallel theories are obtained that could be useful to test these models and to establish a link between Loop Quantum Cosmology and Brane Cosmology. It is concluded that a joint analysis as the one developed here allows to estimate, in a very convenient way, possible deviation of the real universe cosmology from the standard Friedmann-Robertson-Walker one.

We demonstrate how holonomy corrections in loop quantum cosmology (LQC) prevent the Big Rip singularity by introducing a quadratic modification in terms of the energy density $\rho$ in the Friedmann equation in the Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) space-time in a consistent and useful way. In addition, we investigate whether other kind of singularities like Type II,III and IV singularities survive or are avoided in LQC when the universe is filled by a barotropic fluid with the state equation $P=-\rho-f(\rho)$, where $P$ is the pressure and $f(\rho)$ a function of $\rho$. It is shown that the Little Rip cosmology does not happen in LQC. Nevertheless, the occurrence of the Pseudo-Rip cosmology, in which the phantom universe approaches the de Sitter one asymptotically, is established, and the corresponding example is presented. It is interesting that the disintegration of bound structures in the Pseudo-Rip cosmology in LQC always takes more time than that in Einstein cosmology. Our investigation on future singularities is generalized to that in modified teleparallel gravity, where LQC and Brane Cosmology in the Randall-Sundrum scenario are the best examples. It is remarkable that $F(T)$ gravity may lead to all the kinds of future singularities including Little Rip.