Recently Aadne and Gr{\o}n have argued that dark energy may follow naturally from a Nash theory of gravity. In this brief note I argue why this cannot be the case.

An infinite class of exact static anisotropic spheres is developed. All members of the class satisfy (i) regularity (meaning no singularities), and in particular at the origin, (ii) positive but monotone decreasing energy density ($\rho(r)$), radial pressure ($p(r)$), and tangential pressure ($P(r)$), (iii) a finite value of $r=R$ such that $p(R)=0$ defining the boundary surface onto vacuum, (iv) $p \leq \rho$, and (v) $p + 2 P=3 \rho$. All standard energy conditions are satisfied except for the dominant energy condition which has an innocuous violation by the tangential stress since $\rho \leq P$ by construction. An infinite number of the solutions violate the Buchdal bound.

We provide an invariant characterization of the physical properties of the Kerr spacetime. We introduce two dimensionless invariants, constructed out of some known curvature invariants, that act as detectors for the event horizon and ergosurface of the Kerr black hole. We also show that the mass and angular momentum can be extracted from local measurements of the curvature invariants, which in the weak field limit could be used to approximate the total angular momentum and mass of a system of merging black holes. Finally, we introduce a dimensionless invariant that gives a local measure of the "Kerrness" of the spacetime.

The Kerr metric is one of the most important solutions to Einstein's field equations, describing the gravitational field outside a rotating black hole. We thoroughly analyze the curvature scalar invariants to study the Kerr spacetime by examining and visualizing their covariant gradient fields. We discover that the part of the Kerr geometry outside the black hole horizon changes qualitatively depending on the spin parameter, a fact previously unknown. The number of observable critical points of the curvature invariants' gradient fields along the axis of rotation changes at several transitional values of the spin parameter. These transitional values are a fundamental property of the Kerr metric. They are physically important since in general relativity these curvature invariants represent the cumulative tidal and frame-dragging effects of rotating black holes in an observer-independent way.

This is the first in a series of papers in which the gradient flows of fundamental curvature invariants are used to formulate a visualization of curvature. We start with the construction of strict Newtonian analogues (not limits) of solutions to Einstein's equations based on the topology of the associated gradient flows. We do not start with any easy case. Rather, we start with the Curzon - Chazy solution, which, as history shows, is one of the most difficult exact solutions to Einstein's equations to interpret physically. We show that the entire field of the Curzon - Chazy solution, up to a region very "close" to the the intrinsic singularity, strictly represents that of a Newtonian ring, as has long been suspected. In this regard, we consider our approach very successful. As regrades the local structure of the singularity of the Curzon - Chazy solution within a fully general relativistic analysis, however, whereas we make some advances, the full structure of this singularity remains incompletely resolved.

Traditional approaches to the study of the dynamics of spacetime curvature in a very real sense hide the intricacies of the nonlinear regime. Whether it be huge formulae, or mountains of numerical data, standard methods of presentation make little use of our remarkable skill, as humans, at pattern recognition. Here we introduce a new approach to the visualization of spacetime curvature. We examine the flows associated with the gradient fields of invariants derived from the spacetime. These flows reveal a remarkably rich structure, and offer fresh insights even for well known analytical solutions to Einstein's equations. This paper serves as an overview and as an introduction to this approach.

Whereas current cosmological observations suggest that the universe is dominated by a positive cosmological constant ($\Lambda > 0$), the AdS/CFT correspondence tells us that the case $\Lambda<0$ is still worthy of consideration. In this paper we study the McVittie solution with $\Lambda<0$. Following a related study, the solution is understood here by way of a systematic construction of conformal diagrams based on detailed numerical integrations of the null geodesic equations. As in the pure Robertson - Walker case, we find that $\Lambda<0$ ensures collapse to a Big Crunch, a feature which completely dominates the global structure.

The generalized Swiss - cheese model, consisting of a Lema\^itre - Tolman (inhomogeneous dust) region matched, by way of a comoving boundary surface, onto a Robertson-Walker background of homogeneous dust, has become a standard construction in modern cosmology. Here we ask if this construction can be made more realistic by introducing some evolution of the boundary surface. The answer we find is no. To maintain a boundary surface using the Darmois - Israel junction conditions, as opposed to the introduction of a surface layer, the boundary must remain exactly comoving. The options are to drop the assumption of dust or allow the development of surface layers. Either option fundamentally changes the original construction.

Recently Kaloper, Kleban and Martin reexamined the McVittie solution and argued, contrary to a very widely held belief, that the solution contains a black hole in an expanding universe. Here we corroborate their main conclusion but go on to examine, in some detail, a specific solution that asymptotes to the $\Lambda$CDM cosmology. We show that part of the boundary of the solution contains the inner bifurcation two - sphere of the Schwarzschild - de Sitter spacetime and so both the black and white hole horizons together form a partial boundary of this McVittie solution. We go on to show that the null and weak energy conditions are satisfied and that the dominant energy condition is satisfied almost everywhere in the solution. The solution is understood here by way of a systematic construction of a conformal diagram based on detailed numerical integrations of the null geodesic equations. We find that the McVittie solution admits a degenerate limit in which the bifurcation two - sphere disappears. For solutions with zero cosmological constant, we find no evidence for the development of a weak null singularity. Rather, we find that in this case there is either a black hole to the future of an initial singularity or a white hole to its past.

On the basis of the Kerr metric as a model for a spinning black hole accreting test particles from rest at infinity, I show that the center-of-mass energy for a pair of colliding particles is generically divergent at the inner horizon. This shows that not only are classical black holes internally unstable, but also that Planck-scale physics is a characteristic feature within black holes at scales much larger that the Planck length. The novel feature of the divergence discussed here is that the phenomenon is present only for black holes with rotation and in this sense it is distinct from the well known Cauchy horizon instability.

The invariant four-volume $\mathcal{V}$ of a complete black hole (the volume of the spacetime at and interior to the horizon) diverges. However, if one considers the black hole set up by the gravitational collapse of an object, and integrates only a finite time to the future of the collapse, the resultant volume is well defined and finite. In this paper we examine non-degenerate stationary black holes (and cosmological horizons) and find that $\mathcal{V}_{s} \varpropto \ln(\lambda)$ where $s$ is any shell that terminates on the horizon, $\lambda$ is the affine generator of the horizon and the constant of proportionality is the Parikh volume of $s$ divided by the surface gravity. This provides an alternative local and invariant definition of the surface gravity of a stationary black hole.

We generalize the Swiss-cheese cosmologies so as to include nonzero linear momenta of the associated boundary surfaces. The evolution of mass scales in these generalized cosmologies is studied for a variety of models for the background without having to specify any details within the local inhomogeneities. We find that the final effective gravitational mass and size of the evolving inhomogeneities depends on their linear momenta but these properties are essentially unaffected by the details of the background model.

I use the Newtonian equation of hydrostatic equilibrium for an isotropic fluid sphere to generate exact anisotropic solutions of Einstein's equations. The input function is simply the density. An infinite number of regular solutions are constructed, some of which satisfy all the standard energy conditions. Two classes of these solutions generalize the Newtonian polytropes of index 0 and 1.

We generate an explicit four-fold infinity of physically acceptable exact perfect fluid solutions of Einstein's equations by way of conformal transformations of physically unacceptable solutions (one way to view the use of isotropic coordinates). Special cases include the Schwarzschild interior solution and the Einstein static universe. The process we consider involves solving two equations of the Riccati type coupled by a single generating function rather than a specification of one of the two metric functions.

Following a solution generating technique introduced recently by one of us, we transform the Einstein static Universe into a two - fold infinity class of physically acceptable exact perfect fluid solutions of Einstein's equations. Whereas the entire class of solutions can be considered as generalizations of the familiar Tolman IV solution, no member of the class can be written explicitly in isotropic coordinates. Further, except for a set of measure zero, no member of the class can be written explicitly in curvature coordinates either.

The staid subject of exact static spherically symmetric perfect fluid solutions of Einstein's equations has been reinvigorated in the last decade. We now have several solution generating techniques which give rise to new exact solutions. Here the Einstein static Universe is transformed into a physically acceptable solution the properties of which are examined in detail. The emphasis here is on the importance of the integration constants that these generating techniques introduce.

Recently, Rindler and Ishak have argued that the bending of light is, in principle, changed by the presence of a cosmological constant since one must consider not only the null geodesic equation, but also the process of measurement. I agree with the fact that both must be considered. Here, on the basis of the mathematically exact solution to the classical bending problem, and independent of the cosmological constant, I clarify the approximate argument found in the vast majority of texts (new and old) for the measured value of the bending of light for a single source and show that the result is in part due to the almost perfect cancelation of two terms, one of which is seldom considered. When one considers two sources, this cancelation is of no consequence, and, for example, if the sources are opposite with the same associated apsidal distance, the approximate argument gives the rigorously correct answer (up to numerical evaluation), an answer which is unaffected by the presence of a cosmological constant.

With the many ambitious proposals afoot for new generations of very large telescopes, along with spectrographs of unprecedented resolution, there arises the real possibility that the time evolution of the cosmological redshift may, in the not too distant future, prove to be a useful tool rather than merely a theoretical curiosity. Here I contrast this approach with the standard cosmological procedure based on the luminosity (or any other well-defined) distance. I then show that such observations would not only provide a direct measure of all the associated cosmological parameters of the LCDM model, but would also provide wide-ranging internal consistency checks. Further, in a more general context, I show that without introducing further time derivatives of the redshift one could in fact map out the dark energy equation of state should the LCDM model fail. A consideration of brane-world scenarios and interacting dark energy models serves to emphasize the fact that the usefulness of such observations would not be restricted to high redshifts.

An explicit and complete set of constants of the motion are constructed algorithmically for Friedmann-Lema\^{i}tre-Robertson-Walker (FLRW) models consisting of an arbitrary number of non-interacting species. The inheritance of constants of the motion from simpler models as more species are added is stressed. It is then argued that all FLRW models admit what amounts to a unique candidate for a gravitational epoch function (a dimensionless scalar invariant derivable from the Riemann tensor without differentiation which is monotone throughout the evolution of the universe). The same relations that lead to the construction of constants of the motion allow an explicit evaluation of this function. In the simplest of all models, the $\Lambda$CDM model, it is shown that the epoch function exists for all models with $\Lambda > 0$, but for almost no models with $\Lambda \leq 0$.

It is shown that an Einstein cluster of WIMPs, WIMPs on stable circular geodesic orbits generating the static spherically symmetric gravitational field of a galactic halo, can exactly reproduce the rotation curve of any galaxy simply by adjusting the local angular momentum distribution and consequent number distribution of the WIMPs. No new physics is involved (assuming the forthcoming discovery of WIMPs). Further, stability of the orbits can require an inner truncation of the halo and an explicit example of this is given. There is no exact Newtonian counterpart to the model presented since pure Newtonian gravity fails to account for the contribution made by the angular momentum distribution of the particles in creating the gravitational field. In effect, a galactic halo is supported by the hoop stresses created by the orbiting WIMPs.

GRTensorJ - Books is an active interface to a small part of the computer algebra systems GRTensorII (for Maple) and GRTensorM (for Mathematica) with the specific intent of providing students of General Relativity with an advanced programmable calculator-style tool. All standard functions associated with a classical tensor approach to the subject are available thus reducing these to "elementary functions". This is not a traditional database. The database entries are spacetimes and calculations are done in real time. All spacetimes are referenced directly by equation number in ten current (and classic) texts in notation as close as possible to the original text. The tool is now available free of charge from \texttt{grtensor.org/teaching} .

By way of a complete integration of the Friedmann equations, in terms of observables, it is shown that for the cosmological constant $\Lambda > 0$ there exist non-flat FLRW models for which the total density parameter $\Omega$ remains $\sim 1$ throughout the entire history of the universe. Further, it is shown that in a precise quantitative sense these models are not finely tuned. When observations are brought to bear on the theory, and in particular the WMAP observations, they confirm that we live in just such a universe. The conclusion holds when the classical notion of $\Lambda$ is extended to dark energy.

Exploiting the existence of two "cosmological" constants, one associated with the classical Lovelock theorem and one with the vacuum energy density, we argue, in a model independent way, that in spatially closed FLRW cosmologies with a positive definite effective cosmological constant there exists a range in this constant that serves as a sufficient condition for the satisfaction of the null, weak, strong and dominant energy conditions at a bounce. The application of energy conditions is not unambiguous and we show how the bounce can be considered classically and how, we believe more reasonably, it can be considered a matter of quantum cosmology.

The standard energy conditions of classical general relativity are applied to FLRW cosmologies containing sudden future singularities. Here we show, in a model independent way, that although such cosmologies can satisfy the null, weak and strong energy conditions, they always fail to satisfy the dominant energy condition. They require a divergent spacelike energy flux in all but the comoving frame.

The information contained in galactic rotation curves is examined under a minimal set of assumptions. If emission occurs from stable circular geodesic orbits of a static spherically symmetric field, with information propagated to us along null geodesics, observed rotation curves determine galactic potentials without specific reference to any metric theory of gravity. Given the potential, the gravitational mass can be obtained by way of an anisotropy function of this field. The gravitational mass and anisotropy function can be solved for simultaneously in a Newtonian limit without specifying any specific source. This procedure, based on a minimal set of assumptions, puts very strong constraints on any model of the "dark matter".

We show that a flow (timelike congruence) in any type $B_{1}$ warped product spacetime is uniquely and algorithmically determined by the condition of zero flux. (Though restricted, these spaces include many cases of interest.) The flow is written out explicitly for canonical representations of the spacetimes. With the flow determined, we explore an inverse approach to Einstein's equations where a phenomenological fluid interpretation of a spacetime follows directly from the metric irrespective of the choice of coordinates. This approach is pursued for fluids with anisotropic pressure and shear viscosity. In certain degenerate cases this interpretation is shown to be generically not unique. The framework developed allows the study of exact solutions in any frame without transformations. We provide a number of examples, in various coordinates, including spacetimes with and without unique interpretations. The results and algorithmic procedure developed are implemented as a computer algebra program called GRSource.

An algorithm based on the choice of a single monotone function (subject to boundary conditions) is presented which generates all regular static spherically symmetric perfect fluid solutions of Einstein's equations. For physically relevant solutions the generating functions must be restricted by non-trivial integral-differential inequalities. Nonetheless, the algorithm is demonstrated here by the construction of an infinite number of previously unknown physically interesting exact solutions.

A class of spherically symmetric spacetimes invariantly defined by a zero flux condition is examined first from a purely geometrical point of view and then physically by way of Einstein's equations for a general fluid decomposition of the energy-momentum tensor. The approach, which allows a formal inversion of Einstein's equations, explains, for example, why spherically symmetric perfect fluids with spatially homogeneous energy density must be shearfree.

The r-mode frequencies of the Tolman VII solution for the slowly rotating non-barotropic approximation within the low frequency regime are estimated. The relativistic correction to Newtonian r-mode calculations is shown as function of the tenuity $\frac{R}{M}$ and is shown to be significant only for very compact neutron stars.

The Tolman VII solution is an exact static spherically symmetric perfect fluid solution of Einstein's equations that exhibits a surprisingly good approximation to a neutron star. We show that this solution exhibits trapped null orbits in a causal region even for a tenuity (total radius to mass ratio) $> 3$. In this region the dynamical part of the potential for axial w - modes dominates over the centrifugal part.

An explicit necessary condition for the occurrence of resonance scattering of axial gravitational waves, along with the internal trapping of null geodesics, is proposed for static spherically symmetric perfect fluid solutions to Einstein's equations. Some exact inhomogeneous solutions which exhibit this trapping are given with special attention to boundary conditions and the physical acceptability of the space times. In terms of the tenuity ($\alpha = R/M$ at the boundary) all the examples given lie in the narrow range $2.8 < \alpha < 2.9$. The tenuity can be raised to more interesting values by the addition of an envelope without altering the trapping.

This paper explores properties of the instantaneous ergo surface of a Kerr black hole. The surface area is evaluated in closed form. In terms of the mass ($m$) and angular velocity ($a$), to second order in $a$, the area of the ergo surface is given by $16 \pi m^2 + 4 \pi a^2$ (compared to the familiar $16 \pi m^2 - 4 \pi a^2$ for the event horizon). Whereas the total curvature of the instantaneous event horizon is $4 \pi$, on the ergo surface it ranges from $4 \pi$ (for $a=0$) to 0 (for $a=m$) due to conical singularities on the axis ($\theta=0,\pi$) of deficit angle $2 \pi (1-\sqrt{1-(a/m)^2})$. A careful application of the Gauss-Bonnet theorem shows that the ergo surface remains topologically spherical. Isometric embeddings of the ergo surface in Euclidean 3-space are defined for $0 \leq a/m \leq 1$ (compared to $0 \leq a/m \leq \sqrt{3}/2$ for the horizon).

The recent analysis of Markovic and Shapiro on the effect of a cosmological constant on the evolution of a spherically symmetric homogeneous dust ball is extended to include the inhomogeneous and degenerate cases. The histories are shown by way of effective potential and Penrose-Carter diagrams.

GRLite and GRTensorJ are first and second generation graphical user interfaces to the computer algebra system GRTensorII. Current development centers on GRTensorJ, which provides fully customizable symbolic procedures that reduce many complex calculations to "elementary functions". Although still in development, GRTensorJ, which is now available (free of charge) over the internet, is sufficiently advanced to be of interest to researchers in general relativity and related fields.

We ask the following question: Of the exact solutions to Einstein's equations extant in the literature, how many could represent the field associated with an isolated static spherically symmetric perfect fluid source? The candidate solutions were subjected to the following elementary tests: i) isotropy of the pressure, ii) regularity at the origin, iii) positive definiteness of the energy density and pressure at the origin, iv) vanishing of the pressure at some finite radius, v) monotonic decrease of the energy density and pressure with increasing radius, and vi) subluminal sound speed. A total of 127 candidate solutions were found. Only 16 of these passed all the tests. Of these 16, only 9 have a sound speed which monotonically decreases with radius. The analysis was facilitated by use of the computer algebra system GRTensorII.

We present the GRjunction package which allows boundary surfaces and thin-shells in general relativity to be studied with a computer algebra system. Implementing the Darmois-Israel thin shell formalism requires a careful selection of definitions and algorithms to ensure that results are generated in a straight-forward way. We have used the package to correctly reproduce a wide variety of examples from the literature. We present several of these verifications as a means of demonstrating the packages capabilities. We then use GRjunction to perform a new calculation - joining two Kerr solutions with differing masses and angular momenta along a thin shell in the slow rotation limit.

Motivated by the recent determinations of the Hubble constant $(H_0)$ from observations of Cepheid variables in NGC4571 and M100, we plot the Tinsley diagram with level curves of the cosmological constant $(\Lambda)$. Based on current estimates of the absolute ages of globular clusters we conclude that $\Lambda > 0$ and, irrespective of the background spatial curvature, the universe will not recollapse. These conclusions hold for both relativistic and Newtonian models and are {\it independent} of the density parameter.