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.

Energetic feedback from active galactic nuclei (AGNs) is often used in simulations to resolve several outstanding issues in galaxy formation, but its impact is still not fully understood. Here we derive new constraints on AGN feedback by comparing observations and simulations of the thermal Sunyaev-Zel'dovich (tSZ) effect. We draw on observational results presented in Spacek et al. (2016, 2017) who used data from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) to measure the tSZ signal from >= 10^11 M_Sun and >= 1 Gyr galaxies at z=0.5-1.0 (low-z) and z=1.0-1.5 (high-z). Using the large-scale cosmological hydrodynamical simulations Horizon-AGN and Horizon-NoAGN, which include and omit AGN feedback, we extract simulated tSZ measurements around galaxies equivalent to the observational work. We find that the Horizon-AGN results only differ from the SPT measurements at levels of 0.4 sigma (low-z) and 0.6 sigma (high-z), but differ from the ACT measurements by 3.4 sigma (low-z) and 2.3 sigma (high-z). The Horizon-NoAGN results provide a slightly better fit to the SPT measurements by differing by 0.2 sigma (low-z) and 0.4 sigma (high-z), but a significantly better match to the ACT measurements by differing by only 0.5 sigma (low-z) and 1.4 sigma (high-z). We conclude that, while the lower-mass (<~ 5 x 10^11 M_Sun) SPT results allow for the presence AGN feedback energy, the higher-mass (>~ 5 x 10^11 M_Sun) ACT results show significantly less energy than predicted in the simulation including AGN feedback, while more closely matching the simulation without AGN feedback, indicating that AGN feedback may be milder than often predicted in simulations.

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.

We directly measure the thermal energy of the gas surrounding galaxies through the thermal Sunyaev-Zel'dovich (tSZ) effect. We perform a stacking analysis of microwave background images from the Atacama Cosmology Telescope, around 1179 massive quiescent elliptical galaxies at 0.5 <= z <= 1.0 ('low-z') and 3274 galaxies at 1.0 <= z <= 1.5 ('high-z'), selected using data from the Wide-Field Infrared Survey Explorer All-Sky Survey and the Sloan Digital Sky Survey (SDSS) within the SDSS Stripe-82 field. The gas surrounding these galaxies is expected to contain energy from past episodes of active galactic nucleus (AGN) feedback, and after using modeling to subtract undetected contaminants, we detect a tSZ signal at a significance of 0.9-sigma for our low-z galaxies and 1.8-sigma for our high-z galaxies. We then include data from the high-frequency Planck bands for a subset of 227 low-z galaxies and 529 high-z galaxies and find low-z and high-z tSZ detections of 1.0-sigma and 1.5-sigma, respectively. These results indicate an average thermal heating around these galaxies of 5.6(+5.9/-5.6) x 10^60 erg for our low-z galaxies and 7.0(+4.7/-4.4) x 10^60 erg for our high-z galaxies. Based on simple heating models, these results are consistent with gravitational heating without additional heating due to AGN feedback.

Energetic feedback due to active galactic nuclei (AGN) is likely to play an important role in the observed anti-hierarchical trend in the evolution of galaxies, and yet the energy injected into the circumgalactic medium by this process is largely unknown. One promising approach to constrain this feedback is through measurements of CMB spectral distortions due to the thermal Sunyaev-Zel'dovich (tSZ) effect, whose magnitude is directly proportional to the energy input by AGN. Here we co-add South Pole Telescope SZ (SPT-SZ) survey data around a large set of massive quiescent elliptical galaxies at z >= 0.5. We use data from the Blanco Cosmology Survey and VISTA Hemisphere Survey to create a large catalog of galaxies split up into two redshift bins, with 3394 galaxies at 0.5 <= z <= 1.0 and 924 galaxies at 1.0 <= z <= 1.5, with typical stellar masses of 1.5 x 10^11 M_Sun. We then co-add the emission around these galaxies, resulting in a measured tSZ signal at 2.2 sigma significance for the lower redshift bin and a contaminating signal at 1.1 sigma for the higher redshift bin. To remove contamination due to dust emission, we use SPT-SZ source counts to model a contaminant source population in both the SPT-SZ bands and Planck high-frequency bands for a subset of 937 low-redshift galaxies and 240 high-redshift galaxies. This increases our detection to 3.6 sigma for low redshifts and 0.9 sigma for high redshifts. We find the mean angularly-integrated Compton-y values to be 2.2 (-0.7+0.9) x 10^-7 Mpc^2 for low redshifts and 1.7 (-1.8+2.2) x 10^-7 Mpc^2 for high redshifts, corresponding to total thermal energies of 7.6 (-2.3+3.0) x 10^60 ergs and 6.0 (-6.3+7.7) x 10^60 ergs, respectively. These numbers are higher than expected from simple theoretical models that do not include AGN feedback, and serve as constraints that can be applied to current simulations of massive galaxy formation. (abridged)