This paper explores the stability of an Earth-like planet orbiting a solar-mass star in the presence of a stellar companion using ~ 400,000 numerical integrations. Given the chaotic nature of the systems being considered, we perform a statistical analysis of the ensuing dynamics for ~500 orbital configurations defined by the following set of orbital parameters: the companion mass; the companion eccentricity; the companion periastron; and the planets inclination angle relative to the stellar binary plane. Specifically, we generate a large sample of survival times for each orbital configuration through the numerical integration of N >> 1 equivalent experiments (e.g., with the same orbital parameters but randomly selected initial orbital phases). We then construct distributions of survival time using the variable mu_s = log tau_s (where tau_s is in years) for each orbital configuration. The primary objective of this work is twofold. First, we use the mean of the distributions to gain a better understanding of what orbital configurations, while unstable, have sufficiently long survival times to make them interesting to the study of planet habitability. Second, we calculate the width, skew, and kurtosis of each mu_s distribution and look for general features that may aid further understanding and numerical exploration of these chaotic systems.
Download