Exact solutions of the Dirac general relativistic equation (DE) that describe the dynamics of a massive, electrically charged particle with half-integer spin in the curved spacetime geometry of an electrically charged, rotating Kerr-Newman-(anti) de Sitter black hole (BH) are investigated. We first, derive the DE in the Kerr-Newman-de Sitter (KNdS) BH background using a generalised Kinnersley null tetrad in the Newman-Penrose formalism. In this frame, we prove the separation of the DE into ordinary differential equations for the radial and angular parts. Under specific transformations of the independent and dependent variables we prove that the transformed radial equation for a massive charged spin $\frac{1}{2}$ fermion in the background KNdS BH constitutes a highly non-trivial generalisation of Heun's equation. Using a Regge-Wheeler-like independent variable we transform the radial equation in the KNdS background into a Schr\"{o}dinger-like differential equation and investigate its asymptotic behaviour near the event and cosmological horizons. For a massive fermion (MF) in the background of a Kerr-Newman (KN) BH we first prove that the radial and angular equations that result from the separation of DE reduce to the generalised Heun differential equation (GHE). The local solutions of such GHE are derived and can be described by holomorphic functions whose power series coefficients are determined by a four-term recurrence relation. Using asymptotic analysis we derive the solutions for the MF far away from the KN BH and the solutions near the event horizon . The determination of the separation constant as an eigenvalue problem in the KN background is investigated. Using the aforementioned four-term recursion formula we prove that in the non-extreme KN geometry there are no bound states with $\omega^2<\mu^2$, where $\omega$ and $\mu$ are the energy and mass of the fermion respectively.

The null geodesics that describe photon orbits in the spacetime of a rotating electrically charged black hole (Kerr-Newman) are solved exactly including the contribution from the cosmological constant. We derive elegant closed form solutions for relativistic observables such as the deflection angle and frame dragging effect that a light ray experiences in the gravitational fields (i) of a Kerr-Newman black hole and (ii) of a Kerr-Newman-de Sitter black hole. We then solve the more involved problem of treating a Kerr-Newman black hole as a gravitational lens, i.e. a KN black hole along with a static source of light and a static observer both located far away but otherwise at arbitrary positions in space. For this model, we derive the analytic solutions of the lens equations in terms of Appell and Lauricella hypergeometric functions and the Weierstra\ss modular form. The exact solutions derived for null, spherical polar and non-polar orbits, are applied for the calculation of frame dragging for the orbit of a photon around the galactic centre, assuming that the latter is a Kerr-Newman black hole. We also apply the exact solution for the deflection angle of an equatorial light ray in the gravitational field of a Kerr-Newman black hole for the calculation of bending of light from the gravitational field of the galactic centre for various values of the Kerr parameter, electric charge and impact factor. In addition, we derive analytic expressions for the Maxwell tensor components for a Zero-Angular-Momentum-Observer (ZAMO) in the Kerr-Newman-de Sitter spacetime.

The null geodesic equations that describe motion of photons in Kerr spacetime are solved exactly in the presence of the cosmological constant $\Lambda$. The exact solution for the deflection angle for generic light orbits (i.e. non-polar, non-equatorial) is calculated in terms of the generalized hypergeometric functions of Appell and Lauricella. We then consider the more involved issue in which the black hole acts as a `gravitational lens'. The constructed Kerr black hole gravitational lens geometry consists of an observer and a source located far away and placed at arbitrary inclination with respect to the black hole's equatorial plane. The resulting lens equations are solved elegantly in terms of Appell-Lauricella hypergeometric functions and the Weierstrass elliptic function. We then, systematically, apply our closed form solutions for calculating the image and source positions of generic photon orbits that solve the lens equations and reach an observer located at various values of the polar angle for various values of the Kerr parameter and the first integrals of motion. In this framework, the magnification factors for generic orbits are calculated in closed analytic form for the first time. The exercise is repeated with the appropriate modifications for the case of a non-zero cosmological constant.

The exact solution for the motion of a test particle in a non-spherical polar orbit around a Kerr black hole is derived. Exact novel expressions for frame dragging (Lense-Thirring effect), periapsis advance and the orbital period are produced. The resulting formulae, are expressed in terms of Appell's first hypergeometric function $F_1$, Jacobi's amplitude function, and Appell's $F_1$ and Gau$\ss$ hypergeometric function respectively. The exact expression for frame dragging is applied for the calculation of the Lense-Thirring effect for the orbits of S-stars in the central arcsecond of our Galaxy assuming that the galactic centre is a Kerr black hole, for various values of the Kerr parameter including those supported by recent observations. In addition, we apply our solutions for the calculation of frame dragging and periapsis advance for stellar non-spherical polar orbits in regions of strong gravitational field close to the event horizon of the galactic black hole, e.g. for orbits in the central milliarcsecond of our galaxy. Such orbits are the target of the GRAVITY experiment. We provide examples with orbital periods in the range of 100min - 54 days. Detection of such stellar orbits will allow the possibility of measuring the relativistic effect of periapsis advance with high precision at the strong field realm of general relativity. Further, an exact expression for the orbital period of a test particle in a non-circular equatorial motion around a Kerr black hole is produced. We also derive exact expressions for the periapsis advance and the orbital period for a test particle in a non-circular equatorial motion in the Kerr field in the presence of the cosmological constant in terms of Lauricella's fourth hypergeometric function $F_D$.

The timelike geodesic equations resulting from the Kerr gravitational metric element are derived and solved exactly including the contribution from the cosmological constant. The geodesic equations are derived, by solving the Hamilton-Jacobi partial differential equation by separation of variables. The solutions can be applied in the investigation of the motion of a test particle in the Kerr and Kerr-(anti) de Sitter gravitational fields. In particular, we apply the exact solutions of the timelike geodesics i) to the precise calculation of dragging (Lense-Thirring effect) of a satellite's spherical polar orbit in the gravitational field of Earth assuming Kerr geometry, ii) assuming the galactic centre is a rotating black hole we calculate the precise dragging of a stellar polar orbit aroung the galactic centre for various values of the Kerr parameter including those supported by recent observations. The exact solution of non-spherical geodesics in Kerr geometry is obtained by using the transformation theory of elliptic functions. The exact solution of spherical polar geodesics with a nonzero cosmological constant can be expressed in terms of Abelian modular theta functions that solve the corresponding Jacobi's inversion problem.

The geodesic equations resulting from the Schwarzschild gravitational metric element are solved exactly including the contribution from the Cosmological constant. The exact solution is given by genus 2 Siegelsche modular forms. For zero cosmological constant the hyperelliptic curve degenerates into an elliptic curve and the resulting geodesic is solved by the Weierstra$\ss$ Jacobi modular form. The solution is applied to the precise calculation of the perihelion precession of the orbit of planet Mercury around the Sun.

Strong field (exact) solutions of the gravitational field equations of General Relativity in the presence of a Cosmological Constant are investigated. In particular, a full exact solution is derived within the inhomogeneous Szekeres-Szafron family of space-time line element with a nonzero Cosmological Constant. The resulting solution connects, in an intrinsic way, General Relativity with the theory of modular forms and elliptic curves. The homogeneous FLRW limit of the above space-time elements is recovered and we solve exactly the resulting Friedmann Robertson field equation with the appropriate matter density for generic values of the Cosmological Constant %Lambda and curvature constant K. A formal expression for the Hubble constant is derived. The cosmological implications of the resulting non-linear solutions are systematically investigated. Two particularly interesting solutions i) the case of a flat universe K=0, Lambda not= 0 and ii) a case with all three cosmological parameters non-zero, are described by elliptic curves with the property of complex multiplication and absolute modular invariant j=0 and 1728, respectively. The possibility that all non-linear solutions of General Relativity are expressed in terms of theta functions associated with Riemann-surfaces is discussed.

This paper describes how the non-gravitational contribution to Galactic Velocity Rotation Curves can be explained in terms of a negative Cosmological Constant ($\Lambda$). It will be shown that the Cosmological Constant leads to a velocity contribution proportional to the radii, at large radii, and depending on the mass of the galaxy. This explanation contrasts with the usual interpretation that this effect is due to Dark Matter halos. The velocity rotation curve for the galaxy NGC 3198 will be analysed in detail, while several other galaxies will be studied superficially. The Cosmological Constant derived experimentally from the NGC 3198 data was found to be:$|\Lambda|_{Exp}= 5.0\times 10^{-56} cm^{-2}$. This compares favourably with the theoretical value obtained from the Large Number Hypothesis of: $|\Lambda|_{Theory}=2.1\times 10^{-56}cm^{-2}$. The Extended LNH is then used to define other cosmological parameters: gravitational modification constant, energy density, and the Cosmological Constant in terms of a fundamental length. A speculative theory for the evolution of the Universe is outlined where it is shown how the Universe can be defined, in any particular era, by two parameters: the fundamental length and the energy density of the vacuum for that epoch. The theory is applied to the time evolution of the universe where a possible explanation for the $\rho_{Planck}/\rho_{\Lambda}^{QH} \approx 10^{120}$ problem is proposed. The nature of the ''vacuum'' is reviewed along with a speculative approach for calculating the Cosmological Constant via formal M-theory.The experimentally derived results presented in this paper support a decelerating Universe, in contrast with recent indicationsfrom Type Ia Supernovae experiments, for an accelerating Universe.

In this paper a new theory of Dark Matter is proposed. Experimental analysis of several Galaxies show how the non-gravitational contribution to galactic Velocity Rotation Curves can be interpreted as that due to the Cosmological Constant $\Lambda$. The experimentally determined values for $\Lambda$ are found to be consistent with those expected from Cosmological Constraints. The Cosmological Constant is interpreted as leading to a constant energy density which in turn can be used to partly address the energy deficit problem (Dark Energy) of the Universe. The work presented here leads to the conclusion that the Cosmological Constant is negative and that the universe is de-accelerating. This is in clear contradiction to the Type Ia Supernovae results which support an accelerating universe.

Cosmological inflation is studied in the case where the inflaton is the overall modulus $T$ for an orbifold. General forms of the (non-perturbative) superpotential are considered to ensure that $G=K+{\rm ln}|W|^2$ is modular invariant. We find generically that these models do not produce a potential flat enough for slow roll to a supersymmetric minimum, although we do find a model which produces up to 20 e-folds of inflation to a non-supersymmetric minimum.