Description
Using density functional molecular dynamics simulations, we determine the equation of state (EOS) for hydrogen-helium mixtures spanning density-temperature conditions typical of giant-planet interiors, ~0.2-9g/cm3 and 1000-80000K for a typical helium mass fraction of 0.245. In addition to computing internal energy and pressure, we determine the entropy using an ab initio thermodynamic integration technique. A comprehensive EOS table with 391 density-temperature points is constructed and the results are presented in the form of a two-dimensional free energy fit for interpolation. Deviations between our ab initio EOS and the semi-analytical EOS model by Saumon and Chabrier (1992PhRvA..46.2084S, 1995ApJS...99..713S) are analyzed in detail, and we use the results for initial revision of the inferred thermal state of giant planets with known values for mass and radius. Changes are most pronounced for planets in the Jupiter mass range and below. We present a revision to the mass-radius relationship that makes the hottest exoplanets increase in radius by ~0.2 Jupiter radii at fixed entropy and for masses greater than ~0.5 Jupiter mass. This change is large enough to have possible implications for some discrepant "inflated giant exoplanets."
|