The catalog of stellar evolution tracks discussed in our previous work is meant to help characterize exoplanet host-stars of interest for follow-up observations with future missions like JWST. However, the utility of the catalog has been predicated on the assumption that we would precisely know the age of the particular host-star in question; in reality, it is unlikely that we will be able to accurately estimate the age of a given system. Stellar age is relatively straightforward to calculate for stellar clusters, but it is difficult to accurately measure the age of an individual star to high precision. Unfortunately, this is the kind of information we should consider as we attempt to constrain the long-term habitability potential of a given planetary system of interest. This is ultimately why we must rely on predictions of accurate stellar evolution models, as well a consideration of what we can observably measure (stellar mass, composition, orbital radius of an exoplanet) in order to create a statistical framework wherein we can identify the best candidate systems for follow-up characterization. In this paper we discuss a statistical approach to constrain long-term planetary habitability by evaluating the likelihood that at a given time of observation, a star would have a planet in the 2 Gy continuously habitable zone (CHZ2). Additionally, we will discuss how we can use existing observational data (i.e. data assembled in the Hypatia catalog and the Kepler exoplanet host star database) for a robust comparison to the catalog of theoretical stellar models.

We present stellar evolution models for 0.5 - 1.2 \Msol at scaled metallicities of 0.1 - 1.5 Z\sol and O/Fe values of 0.44 - 2.28 O/Fe\sol. The time dependent evolution of habitable zone boundaries are calculated for each stellar evolution track based on stellar mass, effective temperature, and luminosity parameterizations. The rate of change of stellar surface quantities and the surrounding habitable zone position are strong functions of all three quantities explored. The range of orbits that remain continuously habitable, or habitable for at least 2 Gyr, are provided. The results show that the detailed chemical characterization of exoplanet host stars and a consideration of their evolutionary history are necessary to assess the likelihood that a planet found in the instantaneous habitable zone has had sufficient time to develop a biosphere capable of producing detectable biosignatures. This model grid is designed for use by the astrobiology and exoplanet communities to efficiently characterize the time evolution of host stars and their habitable zones for planetary candidates of interest.

Building on previous work, we have expanded our catalog of evolutionary models for stars with variable composition; here we present models for stars of mass 0.5 - 1.2 Msol, at scaled metallicities of 0.1 - 1.5 Zsol, and specific C/Fe, Mg/Fe, and Ne/Fe values of 0.58 - 1.72 C/Fe_sol, 0.54 - 1.84 Mg/Fe_sol and 0.5 - 2.0 Ne/Fe_sol, respectively. We include a spread in abundance values for carbon and magnesium based on observations of their variability in nearby stars; we choose an arbitrary spread in neon abundance values commensurate with the range seen in other low Z elements due to the difficult nature of obtaining precise measurements of neon abundances in stars. As indicated by the results of Truitt et al. (2015), it is essential that we understand how differences in individual elemental abundances, and not just the total scaled metallicity, can measurably impact a star's evolutionary lifetime and other physical characteristics. In that work we found that oxygen abundances significantly impacted the stellar evolution; carbon, magnesium, and neon are potentially important elements to individually consider due to their relatively high (but also variable) abundances in stars. We present 528 new stellar main sequence models, and we calculate the time-dependent evolution of the associated habitable zone boundaries for each based on mass, temperature, and luminosity. We also reintroduce the 2 Gyr "Continuously Habitable Zone" as a useful tool to help gauge the habitability potential for a given planetary system.

{\tau} Ceti (HD10700), a G8 dwarf with solar mass of 0.78, is a close (3.65 pc) sun-like star where 5 possibly terrestrial planet candidates (minimum masses of 2, 3.1, 3.5, 4.3, and 6.7 Earth masses) have recently been discovered. We report abundances of 23 elements using spectra from the MIKE spectrograph on Magellan. Using stellar models with the abundances determined here, we calculate the position of the classical habitable zone with time. At the current best fit age, 7.63 Gy, up to two planets (e and f) may be in the habitable zone, depending on atmospheric properties. The Mg/Si ratio of the star is found to be 1.78, which is much greater than for Earth (about 1.2). With a system that has such an excess of Mg to Si ratio it is possible that the mineralogical make-up of planets around {\tau} Ceti could be significantly different from that of Earth, with possible oversaturation of MgO, resulting in an increase in the content of olivine and ferropericlase compared with Earth. The increase in MgO would have a drastic impact on the rheology of the mantles of the planets around {\tau} Ceti.