The images of many distant galaxies are displaced, distorted and often multiplied by the presence of foreground massive galaxies near the line of sight; the foreground galaxies act as gravitational lenses. Commonly, the lens equation, which relates the placement and distortion of the images to the real source position in the thin-lens scenario, is obtained by extremizing the time of arrival among all the null paths from the source to the observer (Fermat's principle). We show that the construction of envelopes of certain families of null surfaces consitutes an alternative variational principle or version of Fermat's principle that leads naturally to a lens equation in a generic spacetime with any given metric. We illustrate the construction by deriving the lens equation for static asymptotically flat thin lens spacetimes. As an application of the approach, we find the bending angle for moving thin lenses in terms of the bending angle for the same deflector at rest. Finally we apply this construction to cosmological spacetimes (FRW) by using the fact they are all conformally related to Minkowski space.

We present a framework, based on the null-surface formulation of general relativity, for discussing the dynamics of Fermat potentials for gravitational lensing in a generic situation without approximations of any kind. Additionally, we derive two lens equations: one for the case of thick compact lenses and the other one for lensing by gravitational waves. These equations in principle generalize the astrophysical scheme for lensing by removing the thin-lens approximation while retaining the weak fields.

Predictions of the standard thin lens approximation and a new iterative approach to gravitational lensing are compared with an ``exact'' approach in simple test cases involving one or two lenses. We show that the thin lens and iterative approaches are remarkably accurate in predicting time delays, source positions and image magnifications for a single monopole lens and combinations of two monopole lenses. In the cases studied, the iterative method provided greater accuracy than the thin lens method. We also study the accuracy of a ``2 lens, single lens plane model,'' where two monopole lenses colinear with the observer are modeled by a mass distribution in a single lens plane lying between them. We see that this model can lead to large inaccuracies in physically meaningful situations. A previous version of this paper was published as Phys.Rev.D62, 024025, (2000) with errors in the computation of two lens comparisons. This paper corrects these errors and presents new conclusions which differ from the previous version.

We develop an iterative approach to gravitational lensing theory based on approximate solutions of the null geodesic equations. The approach can be employed in any space-time which is ``close'' to a space-time in which the null geodesic equations can be completely integrated, such as Minkowski space-time, Robertson-Walker cosmologies, or Schwarzschild-Kerr geometries. To illustrate the method, we construct the iterative gravitational lens equations and time of arrival equation for a single Schwarzschild lens. This example motivates a discussion of the relationship between the iterative approach, the standard thin lens formulation, and an exact formulation of gravitational lensing.