We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

Jupiter's radio emission has been linked to its planetary-scale magnetic field, and spacecraft investigations have revealed that most planets, and some moons, have or had a global magnetic field. Generated by internal dynamos, magnetic fields are one of the few remote sensing means of constraining the properties of planetary interiors. For the Earth, its magnetic field has been speculated to be partially responsible for its habitability, and knowledge of an extrasolar planet's magnetic field may be necessary to assess its habitability. The radio emission from Jupiter and other solar system planets is produced by an electron cyclotron maser, and detections of extrasolar planetary electron cyclotron masers will enable measurements of extrasolar planetary magnetic fields. This white paper draws heavily on the W. M. Keck Institute for Space Studies report Planetary Magnetic Fields: Planetary Interiors and Habitability (Lazio, Shkolnik, Hallinan, et al.), it incorporates topics discussed at the American Astronomical Society Topical Conference "Radio Exploration of Planetary Habitability," it complements the Astrobiology Science Strategy white paper "Life Beyond the Solar System: Space Weather and Its Impact on Habitable Worlds" (Airapetian et al.), and it addresses aspects of planetary magnetic fields discussed in the NASA Astrobiology Strategy.

We analyze the evolution of the potentially habitable planet Proxima Centauri b to identify environmental factors that affect its long-term habitability. We consider physical processes acting on size scales ranging from the galactic to the stellar system to the planet's core. We find that there is a significant probability that Proxima Centauri has had encounters with its companion stars, Alpha Centauri A and B, that are close enough to destabilize an extended planetary system. If the system has an additional planet, as suggested by the discovery data, then it may perturb planet b's eccentricity and inclination, possibly driving those parameters to non-zero values, even in the presence of strong tidal damping. We also model the internal evolution of the planet, evaluating the roles of different radiogenic abundances and tidal heating and find that magnetic field generation is likely for billions of years. We find that if planet b formed in situ, then it experienced 169 +/- 13 million years in a runaway greenhouse as the star contracted during its formation. This early phase could remove up to 5 times as much water as in the modern Earth's oceans, possibly producing a large abiotic oxygen atmosphere. On the other hand, if Proxima Centauri b formed with a substantial hydrogen atmosphere (0.01 - 1% of the planet's mass), then this envelope could have shielded the water long enough for it to be retained before being blown off itself. After modeling this wide range of processes we conclude that water retention during the host star's pre-main sequence phase is the biggest obstacle for Proxima b's habitability. These results are all obtained with a new software package called VPLANET.

How can scientists conclude with high confidence that an exoplanet hosts life? As telescopes come on line over the next 20 years that can directly observe photons from terrestrial exoplanets, this question will dictate the activities of many scientists across many fields. The expected data will be sparse and with low signal-to-noise, which will make disentangling biosignatures from abiotic features challenging. Our Earth is not just unique in that it hosts life, it is also the only terrestrial planet with direct observations of its interior through seismic waves, and compositional evolution through field and laboratory measurements. This extensive research reveals a planet born from collisions between worlds, followed by a complicated biogeochemical evolution. Exoplanet interiors, on the other hand, can only be constrained by the following observations: 1) photometric and spectroscopic analysis of the planet's atmosphere, 2) spectroscopic and photometric analysis of the host star, and 3) companion planet properties. From these (future) data, astrobiologists must generate plausible compositional and evolutionary models that constrain a potentially habitable exoplanet's internal properties and history, provide environmental context, and rule out geochemical explanations for any putative biosignatures. The goal of this white paper is to frame the role of geophysical and geochemical processes relevant to the search for life beyond the Solar System and to identify critical, but understudied, areas of future research.

Earth's climate, mantle, and core interact over geologic timescales. Climate influences whether plate tectonics can take place on a planet, with cool climates being favorable for plate tectonics because they enhance stresses in the lithosphere, suppress plate boundary annealing, and promote hydration and weakening of the lithosphere. Plate tectonics plays a vital role in the long-term carbon cycle, which helps to maintain a temperate climate. Plate tectonics provides long-term cooling of the core, which is vital for generating a magnetic field, and the magnetic field is capable of shielding atmospheric volatiles from the solar wind. Coupling between climate, mantle, and core can potentially explain the divergent evolution of Earth and Venus. As Venus lies too close to the sun for liquid water to exist, there is no long-term carbon cycle and thus an extremely hot climate. Therefore plate tectonics cannot operate and a long-lived core dynamo cannot be sustained due to insufficient core cooling. On planets within the habitable zone where liquid water is possible, a wide range of evolutionary scenarios can take place depending on initial atmospheric composition, bulk volatile content, or the timing of when plate tectonics initiates, among other factors. Many of these evolutionary trajectories would render the planet uninhabitable. However, there is still significant uncertainty over the nature of the coupling between climate, mantle, and core. Future work is needed to constrain potential evolutionary scenarios and the likelihood of an Earth-like evolution.

Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here we use 1D coupled climate-photochemical models to generate self-consistent atmospheres for evolutionary scenarios predicted in our companion paper (Barnes et al., 2016). These include high-O2, high-CO2, and more Earth-like atmospheres, with either oxidizing or reducing compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b, but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere, and JWST observations longward of 7 microns could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST, but may be within the reach of larger aperture telescopes. Direct imaging spectra may detect O4, which is diagnostic of massive water loss and O2 retention, rather than a photosynthesis. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 {\mu}m would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4 and either photosynthetically produced O2 or a hydrocarbon haze layer.

The internal thermal and magnetic evolution of rocky exoplanets is critical to their habitability. We focus on the thermal-orbital evolution of Earth-mass planets around low mass M stars whose radiative habitable zone overlaps with the "tidal zone". We develop a thermal-orbital evolution model calibrated to Earth that couples tidal dissipation, with a temperature-dependent Maxwell rheology, to orbital circularization and migration. We illustrate thermal-orbital steady states where surface heat flow is balanced by tidal dissipation and cooling can be stalled for billions of years until circularization occurs. Orbital energy dissipated as tidal heat in the interior drives both inward migration and circularization, with a circularization time that is inversely proportional to the dissipation rate. We identify a peak in the internal dissipation rate as the mantle passes through a visco-elastic state at mantle temperatures near 1800 K. Planets orbiting a 0.1 solar-mass star within $0.07$ AU circularize before 10 Gyr, independent of initial eccentricity. Once circular, these planets cool monotonically and maintain dynamos similar to Earth. Planets forced into eccentric orbits can experience a super-cooling of the core and rapid core solidification, inhibiting dynamo action for planets in the habitable zone. We find that tidal heating is insignificant in the habitable zone around $0.45$ (or larger) solar mass stars because tidal dissipation is a stronger function of orbital distance than stellar mass, and the habitable zone is further from larger stars. Suppression of the planetary magnetic field exposes the atmosphere to stellar wind erosion and the surface to harmful radiation. In addition to weak magnetic fields, massive melt eruption rates and prolonged magma oceans may render eccentric planets in the habitable zone of low mass stars inhospitable for life.

We present time series Doppler data from Lick Observatory that reveal the presence of long-period planetary companions orbiting nearby stars. The typical eccentricity of these massive planets are greater than the mean eccentricity of known exoplanets. HD30562b has Msini = 1.29 Mjup, with semi-major axis of 2.3 AU and eccentricity 0.76. The host star has a spectral type F8V and is metal rich. HD86264b has Msini = 7.0 Mjup, arel = 2.86 AU, an eccentricity, e = 0.7 and orbits a metal-rich, F7V star. HD87883b has Msini = 1.78 Mjup, arel = 3.6 AU, e = 0.53 and orbits a metal-rich K0V star. HD89307b has Msini = 1.78 Mjup, arel = 3.3 AU, e = 0.24 and orbits a G0V star with slightly subsolar metallicity. HD148427b has Msini = 0.96 Mjup, arel = 0.93 AU, eccentricity of 0.16 and orbits a metal rich K0 subgiant. We also present velocities for a planet orbiting the F8V metal-rich binary star, HD196885A. The planet has Msini = 2.58 Mjup, arel = 2.37 AU, and orbital eccentricity of 0.48, in agreement with the independent discovery by Correia et al. 2008.

We report the detection of five Jovian mass planets orbiting high metallicity stars. Four of these stars were first observed as part of the N2K program and exhibited low RMS velocity scatter after three consecutive observations. However, follow-up observations over the last three years now reveal the presence of longer period planets with orbital periods ranging from 21 days to a few years. HD 11506 is a G0V star with a planet of \msini = 4.74 \mjup in a 3.85 year orbit. HD 17156 is a G0V star with a 3.12 \mjup planet in a 21.2 day orbit. The eccentricity of this orbit is 0.67, one of the highest known for a planet with a relatively short period. The orbital period for this planet places it in a region of parameter space where relatively few planets have been detected. HD 125612 is a G3V star with a planet of \msini = 3.5 \mjup in a 1.4 year orbit. HD 170469 is a G5IV star with a planet of \msini = 0.67 \mjup in a 3.13 year orbit. HD 231701 is an F8V star with planet of 1.08 \mjup in a 142 day orbit. All of these stars have supersolar metallicity. Three of the five stars were observed photometrically but showed no evidence of brightness variability. A transit search conducted for HD 17156 was negative but covered only 25% of the search space and so is not conclusive.

We report the discovery of two new planets: a 1.94 M_Jup planet in a 1.8-year orbit of HD 5319, and a 2.51 M_Jup planet in a 1.1-year orbit of HD 75898. The measured eccentricities are 0.12 for HD 5319 b and 0.10 for HD 75898 b, and Markov Chain Monte Carlo simulations based on derived orbital parameters indicate that the radial velocities of both stars are consistent with circular planet orbits. With low eccentricity and 1 < a < 2 AU, our new planets have orbits similar to terrestrial planets in the solar system. The radial velocity residuals of both stars have significant trends, likely arising from substellar or low-mass stellar companions.

We report precision Doppler measurements of three intermediate-mass subgiants from Lick and Keck Observatories. All three stars show variability in their radial velocities consistent with planet-mass companions in Keplerian orbits. We find a planet with a minimum mass of 2.5 Mjup in a 351.5 day orbit around HD 192699, a planet with a minimum mass of 2.0 Mjup in a 341.1 day orbit around HD 210702, and a planet with a minimum mass of 0.61 Mjup in a 297.3 day orbit around HD 175541. Stellar mass estimates from evolutionary models indicate that all of these stars were formerly A-type dwarfs with masses ranging from 1.65 to 1.85 Msun. These three long-period planets would not have been detectable during their stars' main-sequence phases due to the large rotational velocities and stellar jitter exhibited by early-type dwarfs. There are now 9 "retired" (evolved) A-type stars (Mstar > 1.6 Msun) with known planets. All 9 planets orbit at distances a \geq 0.78 AU, which is significantly different than the semimajor axis distribution of planets around lower-mass stars. We examine the possibility that the observed lack of close-in planets is due to engulfment by their expanding host stars, but we find that this explanation is inadequate given the relatively small stellar radii of K giants (Rstar < 32 Rsun = 0.15 AU) and subgiants (Rstar < 7 Rsun = 0.03 AU). Instead, we conclude that planets around intermediate-mass stars reside preferentially beyond ~0.8 AU, which may be a reflection of different formation and migration histories of planets around A-type stars.

The N2K consortium is carrying out a distributed observing campaign with the Keck, Magellan and Subaru telescopes, as well as the automatic photometric telescopes of Fairborn Observatory, in order to search for short-period gas giant planets around metal-rich stars. We have established a reservoir of more than 14,000 main sequence and subgiant stars, closer than 110 pc, brighter than V=10.5 and with $0.4 < B-V < 1.2$. Because the fraction of stars with planets is a sensitive function of stellar metallicity, a broadband photometric calibration has been developed to identify a subset of 2000 stars with [Fe/H] $> 0.1$ dex for this survey. We outline the strategy and report the detection of a planet orbiting the metal-rich G5IV star HD 88133 with a period of 3.41 days, semi-velocity amplitude, K=35.7 \mse and \msini = 0.29 \mjupe. Photometric observations reveal that HD 88133 is constant on the 3.415-day radial velocity period to a limit of 0.0005 mag. Despite a transit probability of 15.6%, our photometry rules out the shallow transits predicted by the large stellar radius.