[Paper Review] BinHab: A Numerical Tool for the Calculation of S/P-Type Habitable Zones in Binary Systems

M. Cuntz, R. Bruntz|arXiv (Cornell University)|Sep 11, 2014

Astro and Planetary Science3 citations

TL;DR

BinHab is a publicly accessible numerical tool that calculates S-type and P-type habitable zones in binary star systems by combining radiative habitable zone constraints with orbital stability criteria based on Holman & Wiegert (1999). It enables the identification of habitable zones—classified as S, ST, P, or PT-type—across diverse binary configurations, including low-mass stars, with results showing that habitable zone width decreases linearly with increasing binary eccentricity, and that systems with 0.50 M⊙ stars sustain habitability up to eccentricities of 0.65 for GHZ.

ABSTRACT

The aim of this contribution is to introduce the numerical tool BinHab, a publicly accessible code, available at The University of Texas at Arlington, that allows the calculation of S-type and P-type habitable zones of general binary systems.

Motivation & Objective

• To develop a flexible, publicly accessible numerical tool for calculating habitable zones in binary stellar systems.
• To integrate both radiative habitable zone constraints and orbital stability requirements for putative Earth-like planets.
• To enable the classification of habitable zones as S-type, ST-type, P-type, or PT-type based on system parameters.
• To support studies of habitability in systems with low-mass main-sequence stars, which constitute ~90% of Milky Way stars.
• To extend applicability to evolved stars (e.g., subgiants, giants) and future multi-stellar system modeling.

Proposed method

• The tool computes habitable zones using the radiative habitable zone (RHZ) formalism based on Kasting et al. (1993), with conservative (CHZ), general (GHZ), and extended (EHZ) zone definitions.
• It applies the orbital stability criterion from Holman & Wiegert (1999) as an upper bound for S-type orbits and a lower bound for P-type orbits.
• A fourth-order polynomial solution is used to analytically determine the boundaries of habitable zones in both S-type and P-type configurations.
• The system uses a web-based interface with HTML for input/output and PHP for input validation and server-side processing.
• The core calculations are performed by a Fortran binary, ensuring high-speed execution and security, with inputs and outputs handled via server-side scripts.
• The tool supports both circular and eccentric binary orbits and allows for stellar parameters (mass, temperature, luminosity, radius) to be specified independently for each component.

Experimental results

Research questions

• RQ1 What are the conditions under which S-type habitable zones exist in binary systems with low-mass main-sequence stars?
• RQ2 How does binary eccentricity affect the width and existence of S/ST-type habitable zones in such systems?
• RQ3 To what extent does orbital stability truncate the radiative habitable zone in S-type configurations?
• RQ4 How do different habitable zone definitions (CHZ, GHZ, EHZ) influence the classification and extent of habitable regions?
• RQ5 What is the maximum binary eccentricity for which habitable zones remain viable in systems with 0.50 M⊙ stars?

Key findings

• For systems with 2a_b = 5.0 AU, S/ST-type habitable zones exist only for binary eccentricities below 0.20 when using the GHZ definition.
• Systems with 0.75 M⊙ stars exhibit the broadest radiative habitable zones but experience significant truncation due to orbital stability constraints, resulting in ST-type classification.
• For 0.50 M⊙ binary pairs, habitable zones persist up to eccentricities of 0.65 (GHZ) and 0.62 (CHZ), indicating greater resilience to eccentricity.
• The width of S/ST-type habitable zones decreases linearly as a fraction of binary eccentricity across all mass combinations, due to the Holman & Wiegert stability criterion.
• ST-type habitability is the dominant classification for most models, especially in systems with higher-mass stars, due to strong orbital stability constraints.
• The width of habitable zones is consistently smaller in CHZ models than in GHZ models, as expected from the more restrictive greenhouse gas and stellar flux assumptions.