We generalise Einstein's formulation of the traceless Einstein equations to $f(R)$ gravity theories. As in the traceless version of general relativity, we find that in traceless $f(R)$ gravity there is no vacuum energy cosmological constant problem. The cosmological constant is a mere integration constant that is unconnected to any vacuum energy density or constant term in the gravitational lagrangian. We obtain the integrability condition for the cosmological constant to appear as an integration constant without a connection to any vacuum energy density. We show that, in $D$% -dimensional spacetime, traceless higher-order gravity is conformally equivalent to general relativity and a scalar field $\phi $ with a potential given by the scale-invariant form: $V(\phi )=\frac{D-2}{4D}Re^{-\phi }$, where $\phi =[2/(D-2)]\ln f^{\prime }(R)$. Unlike in the conformal equivalent of full general relativity, flat potentials are found to be possible in all spacetime dimensions for polynomial lagrangians of all orders. Hence, we can solve the cosmological constant vacuum energy problem and have accelerated inflationary expansion in the very early universe with a very small cosmological constant at late times for a wide range of traceless theories. Fine-tunings required in traceless general relativity or standard non-traceless $f(R)$ theories of gravity are avoided. We show that the predictions of the scale-invariant conformal potential are completely consistent with microwave background observational data.

We show that a general solution of the Einstein equations that describes approach to an inhomogeneous and anisotropic sudden spacetime singularity does not experience geodesic incompleteness. This generalises the result established for isotropic and homogeneous universes. Further discussion of the weakness of the singularity is also included.

We provide an extension of the method of asymptotic decompositions of vector fields with finite-time singularities by applying the central extension technique of Poincar\'e to the dominant part of the vector field on approach to the singularity. This leads to a bundle of fan-out asymptotic systems whose equilibria at infinity govern the dynamics of the asymptotic solutions of the original system. We show how this method can be useful to describe a single-fluid isotropic universe at the time of maximum expansion, and discuss possible relations of our results to structural stability and non-compact phase spaces.

We discuss recent developments related to certain blow up methods suitable for the analysis of cosmological singularities and asymptotics. We review results obtained in a variety of currently popular themes and describe ongoing research about universes with various kinds of extreme states, higher order gravity, and certain models of braneworlds.

We study the dynamics near finite-time singularities of flat isotropic universes filled with two interacting but otherwise arbitrary perfect fluids. The overall dynamical picture reveals a variety of asymptotic solutions valid locally around the spacetime singularity. We find the attractor of all solutions with standard decay, and for `phantom' matter asymptotically at early times. We give a number of special asymptotic solutions describing universes collapsing to zero size and others ending at a big rip singularity. We also find a very complicated singularity corresponding with a logarithmic branch point that resembles of cyclic universe, and give an asymptotic local series representation of the general solution in the neighborhood of infinity.

We construct an asymptotic series for a general solution of the Einstein equations near a sudden singularity. The solution is quasi isotropic and contains nine independent arbitrary functions of the space coordinates as required by the structure of the initial value problem.

We present an elementary account of mathematical cosmology through a series of important unsolved problems. We introduce the fundamental notion of `a cosmology' and focus on the issue of singularities as a theme unifying many current, seemingly unrelated trends of this subject. We discuss problems associated with the definition and asymptotic structure of the notion of cosmological solution and also problems related to the qualification of approximations and to the ranges of validity of given cosmologies.

An introductory guide to mathematical cosmology is given focusing on the issue of the genericity of various important results which have been obtained during the last thirty or so years. Some of the unsolved problems along with certain new and potentially powerful methods which may be used for future progress are also given from a unified perspective.