We present a simple choice of integration variables that can be used to exploit the near-integrable character of problems in celestial mechanics. The approach is based on the well-known principle of variation of parameters: instead of orbital elements, we use the phase-space coordinates the object would have at a given point in its (Keplerian) orbit if the perturbing forces were removed. This formulation is suitable for almost any numerical integrator; thus, multistep schemes are easy to build, stepsize can be adjusted, and dissipative forces are allowed. Compared with traditional non-symplectic N-body integrators, the approach often offers increase in speed or accuracy if perturbations are small.

We seek to characterize giant-planet systems by their gravitational scattering properties. We do this to a given system by integrating it numerically along with a large number of hypothetical small bodies that are initially in eccentric habitable zone (HZ)-crossing orbits. Our analysis produces a single number, the escape rate, which represents the rate at which the small-body flux is perturbed away by the giant planets into orbits that no longer pose a threat to terrestrial planets inside the HZ. Obtaining the escape rate this way is similar to computing the largest Liapunov exponent as the exponential rate of divergence of two nearby orbits. For a terrestrial planet inside the HZ, the escape rate value quantifies the "protective" effect that the studied giant-planet system offers. Therefore, escape rates could provide information on whether certain giant-planet configurations produce a more desirable environment for life than the others. We present some computed escape rates on selected planetary systems, focusing on effects of varying the masses and semi-major axes of the giant planets. In the case of our Solar System we find rather surprisingly that Jupiter, in its current orbit, may provide a minimal amount of protection to the Earth.

We have undertaken a thorough dynamical investigation of five extrasolar planetary systems using extensive numerical experiments. The systems Gl 777 A, HD 72659, Gl 614, 47 Uma and HD 4208 were examined concerning the question of whether they could host terrestrial like planets in their habitable zones (=HZ). First we investigated the mean motion resonances between fictitious terrestrial planets and the existing gas giants in these five extrasolar systems. Then a fine grid of initial conditions for a potential terrestrial planet within the HZ was chosen for each system, from which the stability of orbits was then assessed by direct integrations over a time interval of 1 million years. The computations were carried out using a Lie-series integration method with an adaptive step size control. This integration method achieves machine precision accuracy in a highly efficient and robust way, requiring no special adjustments when the orbits have large eccentricities. The stability of orbits was examined with a determination of the Renyi entropy, estimated from recurrence plots, and with a more straight forward method based on the maximum eccentricity achieved by the planet over the 1 million year integration. Additionally, the eccentricity is an indication of the habitability of a terrestrial planet in the HZ; any value of e>0.2 produces a significant temperature difference on a planet's surface between apoapse and periapse. The results for possible stable orbits for terrestrial planets in habitable zones for the five systems are summarized as follows: for Gl 777 A nearly the entire HZ is stable, for 47 Uma, HD 72659 and HD 4208 terrestrial planets can survive for a sufficiently long time, while for Gl 614 our results exclude terrestrial planets moving in stable orbits within the HZ.