An observation of Jupiter's tidal response is anticipated for the on-going Juno spacecraft mission. We combine self-consistent, numerical models of Jupiter's equilibrium tidal response with observed Doppler shifts from the Juno gravity science experiment to test the sensitivity of the spacecraft to tides raised by the Galilean satellites and the Sun. The concentric Maclaurin spheroid (CMS) method finds the equilibrium shape and gravity field of a rotating, liquid planet with the tide raised by a satellite, expanded in Love numbers ($k_{nm}$). We present improvements to CMS theory that eliminate an unphysical center of mass offset and study in detail the convergence behavior of the CMS approach. We demonstrate that the dependence of $k_{nm}$ with orbital distance is important when considering the combined tidal response for Jupiter. Conversely, the details of the interior structure have a negligible influence on $k_{nm}$, for models that match the zonal harmonics $J_2$, $J_4$ and $J_6$, already measured to high precision by Juno. As the mission continues, improved coverage of Jupiter's gravity field at different phases of Io's orbit is expected to yield an observed value for the degree-2 Love number ($k_{22}$) and potentially select higher--degree $k_{nm}$. We present a test of the sensitivity of the Juno Doppler signal to the calculated $k_{nm}$, which suggests the detectability of $k_{33}$, $k_{42}$ and $k_{31}$, in addition to $k_{22}$. A mismatch of a robust Juno observation with the remarkably small range in calculated Io equilibrium $k_{22}=0.58976\pm0.0001$, would indicate a heretofore uncharacterized dynamic contribution to the tides.

The Cassini spacecraft's Grand Finale orbits provided a unique opportunity to probe Saturn's gravity field and interior structure. Doppler measurements yielded unexpectedly large values for the gravity harmonics J_6, J_8, and J_10 that cannot be matched with planetary interior models that assume uniform rotation. Instead we present a suite of models that assume the planet's interior rotates on cylinders, which allows us to match all the observed even gravity harmonics. For every interior model, the gravity field is calculated self-consistently with high precision using the Concentric Maclaurin Spheroid (CMS) method. We present an acceleration technique for this method, which drastically reduces the computational cost, allows us to efficiently optimize model parameters, map out allowed parameter regions with Monte Carlo sampling, and increases the precision of the calculated J_2n gravity harmonics to match the error bars of the observations, which would be difficult without acceleration. Based on our models, Saturn is predicted to have a dense central core of 15-18 Earth masses and an additional 1.5-5 Earth masses of heavy elements in the envelope. Finally, we vary the rotation period in the planet's deep interior and determine the resulting oblateness, which we compare with the value from radio occultation measurements by the Voyager spacecraft. We predict a rotation period of 10:33:34 h +- 55s, which is in agreement with recent estimates derived from ring seismology.

Observations of exoplanet transit spectra are essential to understanding the physics and chemistry of distant worlds. The effects of opacity sources and many physical processes combine to set the shape of a transit spectrum. Two such key processes - refraction and cloud and/or haze forward scattering - have seen substantial recent study. However, models of these processes are typically complex, which prevents their incorporation into observational analyses and standard transit spectrum tools. In this work, we develop analytic expressions that allow for the efficient parameterization of forward scattering and refraction effects in transit spectra. We derive an effective slant optical depth that includes a correction for forward scattered light, and present an analytic form of this correction. We validate our correction against a full-physics transit spectrum model that includes scattering, and we explore the extent to which the omission of forward scattering effects may bias models. Also, we verify a common analytic expression for the location of a refractive boundary, which we express in terms of the maximum pressure probed in a transit spectrum. This expression is designed to be easily incorporated into existing tools, and we discuss how the detection of a refractive boundary could help indicate the background atmospheric composition by constraining the bulk refractivity of the atmosphere. Finally, we show that opacity from Rayleigh scattering and collision induced absorption will outweigh the effects of refraction for Jupiter-like atmospheres whose equilibrium temperatures are above 400-500 K.

The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction of the planet's radius, is helpful in reconciling the calculated $J_n$ with Juno's observations. Although model predictions are strongly affected by the chosen equation of state, the prediction of an enrichment of $Z$ in the deep, metallic envelope over that in the shallow, molecular envelope holds. We estimate Jupiter's core to contain an 7--25 Earth mass of heavy elements. We discuss the current difficulties in reconciling measured $J_n$ with the equations of state, and with theory for formation and evolution of the planet.

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

This article provides an overview of how models of giant planet interiors are constructed. We review measurements from past space missions that provide constraints for the interior structure of Jupiter. We discuss typical three-layer interior models that consist of a dense central core and an inner metallic and an outer molecular hydrogen-helium layer. These models rely heavily on experiments, analytical theory, and first-principle computer simulations of hydrogen and helium to understand their behavior up to the extreme pressures ~10 Mbar and temperatures ~10,000 K. We review the various equations of state used in Jupiter models and compare them with shock wave experiments. We discuss the possibility of helium rain, core erosion and double diffusive convection may have important consequences for the structure and evolution of giant planets. In July 2016 the Juno spacecraft entered orbit around Jupiter, promising high-precision measurements of the gravitational field that will allow us to test our understanding of gas giant interiors better than ever before.

In anticipation of improved observational data for Jupiter's gravitational field from the Juno spacecraft, we predict the static tidal response for a variety of Jupiter interior models based on ab initio computer simulations of hydrogen-helium mixtures. We calculate hydrostatic-equilibrium gravity terms using the non-perturbative concentric Maclaurin Spheroid (CMS) method that eliminates lengthy expansions used in the theory of figures. Our method captures terms arising from the coupled tidal and rotational perturbations, which we find to be important for a rapidly-rotating planet like Jupiter. Our predicted static tidal Love number $k_2 = 0.5900$ is $\sim$10\% larger than previous estimates. The value is, as expected, highly correlated with the zonal harmonic coefficient $J_2$, and is thus nearly constant when plausible changes are made to interior structure while holding $J_2$ fixed at the observed value. We note that the predicted static $k_2$ might change due to Jupiter's dynamical response to the Galilean moons, and find reasons to argue that the change may be detectable, although we do not present here a theory of dynamical tides for highly oblate Jovian planets. An accurate model of Jupiter's tidal response will be essential for interpreting Juno observations and identifying tidal signals from effects of other interior dynamics in Jupiter's gravitational field.

We extend to three dimensions the Concentric Maclaurin Spheroid method for obtaining the self-consistent shape and gravitational field of a rotating liquid planet, to include a tidal potential from a satellite. We exhibit, for the first time, the important effect of the planetary rotation rate on tidal response of gas giants. Simulations of planets with fast rotation rates like those of Jupiter and Saturn, exhibit significant changes in calculated tidal love numbers $k_{nm}$ when compared with non-rotating bodies. A test model of Saturn fitted to observed zonal gravitational multipole harmonics yields $k_2=0.413$, consistent with a recent observational determination from Cassini astrometry data (Lainey et al., 2016). The calculated love number is robust under reasonable assumptions of interior rotation rate, satellite parameters, and details of Saturn's interior structure. The method is benchmarked against several published test cases.

In anticipation of new observational results for Jupiter's axial moment of inertia and gravitational zonal harmonic coefficients from the forthcoming Juno orbiter, we present a number of preliminary Jupiter interior models. We combine results from ab initio computer simulations of hydrogen-helium mixtures, including immiscibility calculations, with a new nonperturbative calculation of Jupiter's zonal harmonic coefficients, to derive a self-consistent model for the planet's external gravity and moment of inertia. We assume helium rain modified the interior temperature and composition profiles. Our calculation predicts zonal harmonic values to which measurements can be compared. Although some models fit the observed (pre-Juno) second- and fourth-order zonal harmonics to within their error bars, our preferred reference model predicts a fourth-order zonal harmonic whose absolute value lies above the pre-Juno error bars. This model has a dense core of about 12 Earth masses, and a hydrogen-helium-rich envelope with approximately 3 times solar metallicity.

I present exact expressions for the interior gravitational potential V of a system of N concentric constant-density (Maclaurin) spheroids. I demonstrate an iteration procedure to find a self-consistent solution for the shapes of the interfaces between spheroids, and for the interior gravitational potential. The external free-space potential, expressed as a multipole expansion, emerges as part of the self-consistent solution. The procedure is both simpler and more precise than perturbation methods. One can choose the distribution and mass densities of the concentric spheroids so as to reproduce a prescribed barotrope to a specified accuracy. I demonstrate the method's efficacy by comparing its results with several published test cases.

Using density functional molecular dynamics simulations, we determine the equation of state for hydrogen-helium mixtures spanning density-temperature conditions typical of giant planet interiors, ~0.2-9 g/cc and 1000-80000 K 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 equation of state (EOS) table with 391 density-temperature points is constructed and the results are presented in form of two-dimensional free energy fit for interpolation. Deviations between our ab initio EOS and the semi-analytical EOS model by Saumon and Chabrier 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 which 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".

The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.

Several processes can cause the shape of an extrasolar giant planet's shadow, as viewed in transit, to depart from circular. In addition to rotational effects, cloud formation, non-homogenous haze production and movement, and dynamical effects (winds) could also be important. When such a planet transits its host star as seen from Earth, the asphericity will introduce a deviation in the transit lightcurve relative to the transit of a perfectly spherical (or perfectly oblate) planet. We develop a theoretical framework to interpret planet shapes. We then generate predictions for transiting planet shapes based on a published theoretical dynamical model of HD189733b. Using these shape models we show that planet shapes are unlikely to introduce detectable lightcurve deviations (those >~1e-5 of the host star), but that the shapes may lead to astrophysical sources of systematic error when measuring planetary oblateness, transit time, and impact parameter.

We apply scintillation theory to stellar signal fluctuations in the high-resolution, high signal/noise, dual-wavelength data from the MMT observation of the 2007 March 18 occultation of P445.3 by Pluto. A well-defined high wavenumber cutoff in the fluctuations is consistent with viscous-thermal dissipation of buoyancy waves (internal gravity waves) in Pluto's high atmosphere, and provides strong evidence that the underlying density fluctuations are governed by the gravity-wave dispersion relation.

Recently two groups used first-principles computer simulations to model Jupiter's interior. While both studies relied on the same simulation technique, density functional molecular dynamics, the groups derived very different conclusions. In particular estimates for the size of Jupiter's core and the metallicity of its hydrogen-helium mantle differed substantially. In this paper, we discuss the differences of the approaches and give an explanation for the differing conclusions.

Hydrogen-helium mixtures at conditions of Jupiter's interior are studied with first-principles computer simulations. The resulting equation of state (EOS) implies that Jupiter possesses a central core of 14-18 Earth masses of heavier elements, a result that supports core accretion as standard model for the formation of hydrogen-rich giant planets. Our nominal model has about 2 Earth masses of planetary ices in the H-He-rich mantle, a result that is, within modeling errors, consistent with abundances measured by the 1995 Galileo Entry Probe mission (equivalent to about 5 Earth masses of planetary ices when extrapolated to the mantle), suggesting that the composition found by the probe may be representative of the entire planet. Interior models derived from this first-principles EOS do not give a match to Jupiter's gravity moment J4 unless one invokes interior differential rotation, implying that jovian interior dynamics has an observable effect on the measured gravity field.

We calculate the theoretical evolution of the radii of all fourteen of the known transiting extrasolar giant planets (EGPs) for a variety of assumptions concerning atmospheric opacity, dense inner core masses, and possible internal power sources. We incorporate the effects of stellar irradiation and customize such effects for each EGP and star. Looking collectively at the family as a whole, we find that there are in fact two radius anomalies to be explained. Not only are the radii of a subset of the known transiting EGPs larger than expected from previous theory, but many of the other objects are smaller than the default theory would allow. We suggest that the larger EGPs can be explained by invoking enhanced atmospheric opacities that naturally retain internal heat. This explanation might obviate the necessity for an extra internal power source. We explain the smaller radii by the presence in perhaps all the known transiting EGPs of dense cores, such as have been inferred for Saturn and Jupiter. Importantly, we derive a rough correlation between the masses of our "best-fit" cores and the stellar metallicity that seems to buttress the core-accretion model of their formation. Though many caveats and uncertainties remain, the resulting comprehensive theory that incorporates enhanced-opacity atmospheres and dense cores is in reasonable accord with all the current structural data for the known transiting giant planets.

The observed mass function for all known extrasolar giant planets (EGPs) varies approximately as M^{-1} for mass M between 0.2 Jupiter masses (M_J) and 5 M_J. In order to study evaporation effects for highly-irradiated EGPs in this mass range, we have constructed an observational mass function for a subset of EGPs in the same mass range but with orbital radii <0.07 AU. Surprisingly, the mass function for such highly-irradiated EGPs agrees quantitatively with the M^{-1} law, implying that the mass function for EGPs is preserved despite migration to small orbital radii. Unless there is a remarkable compensation of mass-dependent orbital migration for mass-dependent evaporation, this result places a constraint on orbital migration models and rules out the most extreme mass loss rates in the literature. A theory that predicts more moderate mass loss gives a mass function that is closer to observed statistics but still disagrees for M < 1 M_J.

We present calculations for the evolution and surviving mass of highly-irradiated extrasolar giant planets (EGPs) at orbital semimajor axes ranging from 0.023 to 0.057 AU using a generalized scaled theory for mass loss, together with new surface-condition grids for hot EGPs and a consistent treatment of tidal truncation. Theoretical estimates for the rate of energy-limited hydrogen escape from giant-planet atmospheres differ by two orders of magnitude, when one holds planetary mass, composition, and irradiation constant. Baraffe et al. (2004, A&A 419, L13-L16) predict the highest rate, based on the theory of Lammer et al. (2003, Astrophys. J. 598, L121-L124). Scaling the theory of Watson et al. (1981, Icarus 48, 150-166) to parameters for a highly-irradiated exoplanet, we find an escape rate ~102 lower than Baraffe's. With the scaled Watson theory we find modest mass loss, occurring early in the history of a hot EGP. In this theory, mass loss including the effect of Roche-lobe overflow becomes significant primarily for masses below a Saturn mass, for semimajor axes = 0.023 AU. This contrasts with the Baraffe model, where hot EGPs are claimed to be remnants of much more massive bodies, originally several times Jupiter and still losing substantial mass fractions at present.

We calculate radius versus age trajectories for the photometrically-selected transiting extrasolar giant planet, OGLE-TR-56b, and find agreement between theory and observation, without introducing an ad hoc extra source of heat in its core. The fact that the radius of HD209458b seems larger than the radii of the recently discovered OGLE family of extremely close-in transiting planets suggests that HD209458b is anomalous. Nevertheless, our good fit to OGLE-TR-56b bolsters the notion that the generic dependence of transit radii on stellar irradiation, mass, and age is, to within error bars, now quantitatively understood.

We build on recent new evolutionary models of Jupiter and Saturn and here extend our calculations to investigate the evolution of extrasolar giant planets of mass 0.15 to 3.0 M_J. Our inhomogeneous thermal history models show that the possible phase separation of helium from liquid metallic hydrogen in the deep interiors of these planets can lead to luminosities ~2 times greater than have been predicted by homogeneous models. For our chosen phase diagram this phase separation will begin to affect the planets' evolution at ~700 Myr for a 0.15 M_J object and ~10 Gyr for a 3.0 M_J object. We show how phase separation affects the luminosity, effective temperature, radii, and atmospheric helium mass fraction as a function of age for planets of various masses, with and without heavy element cores, and with and without the effect of modest stellar irradiation. This phase separation process will likely not affect giant planets within a few AU of their parent star, as these planets will cool to their equilibrium temperatures, determined by stellar heating, before the onset of phase separation. We discuss the detectability of these objects and the likelihood that the energy provided by helium phase separation can change the timescales for formation and settling of ammonia clouds by several Gyr. We discuss how correctly incorporating stellar irradiation into giant planet atmosphere and albedo modeling may lead to a consistent evolutionary history for Jupiter and Saturn.

We present the first models of Saturn and Jupiter to couple their evolution to both a radiative-atmosphere grid and to high-pressure phase diagrams of hydrogen with helium. The purpose of these models is to quantify the evolutionary effects of helium phase separation in Saturn's deep interior. We find that prior calculated phase diagrams in which Saturn's interior reaches a region of predicted helium immiscibility do not allow enough energy release to prolong Saturn's cooling to its known age and effective temperature. We explore modifications to published phase diagrams that would lead to greater energy release, and find a modified H-He phase diagram that is physically reasonable, leads to the correct extension of Saturn's cooling, and predicts an atmospheric helium mass fraction Y_atmos in agreement with recent estimates. We then expand our inhomogeneous evolutionary models to show that hypothetical extrasolar giant planets in the 0.15 to 3.0 Jupiter mass range may have T_effs 10-15 K greater than one would predict with models that do not incorporate helium phase separation.

Using a full frequency-dependent atmosphere code that can incorporate irradiation by a central primary star, we calculate self-consistent boundary conditions for the evolution of the radius of the transiting planet HD 209458b. Using a well-tested extrasolar giant planet evolutionary code, we then calculate the behavior of this planet's radius with age. The measured radius is in fact a transit radius that resides high in HD 209458b's inflated atmosphere. Using our derived atmospheric and interior structures, we find that irradiation plus the proper interpretation of the transit radius can yield a theoretical radius that is within the measured error bars. We conclude that if HD 209458b's true transit radius is at the lower end of the measured range, an extra source of core heating power is not necessary to explain the transit observations.

We present the first models of Jupiter and Saturn to couple their evolution to both a radiative-atmosphere grid and to high-pressure phase diagrams of hydrogen with helium and other admixtures. We find that prior calculated phase diagrams in which Saturn's interior reaches a region of predicted helium immiscibility do not allow enough energy release to prolong Saturn's cooling to its known age and effective temperature. We explore modifications to published phase diagrams that would lead to greater energy release, and propose a modified H-He phase diagram that is physically reasonable, leads to the correct extension of Saturn's cooling, and predicts an atmospheric helium mass fraction Y_atmos = 0.185, in agreement with recent estimates. We also explore the possibility of internal separation of elements heavier than helium, and find that, alternatively, such separation could prolong Saturn's cooling to its known age and effective temperature under a realistic phase diagram and heavy element abundance (in which case Saturn's Y_atmos would be solar but heavier elements would be depleted). In none of these scenarios does Jupiter's interior evolve to any region of helium or heavy-element immiscibility: Jupiter evolves homogeneously to the present day. We discuss the implications of our calculations for Saturn's primordial core mass.

Using a self-consistent atmosphere code, we construct a new model of the atmosphere of the transiting extrasolar giant planet HD 209458b to investigate the disparity between the observed strength of the sodium absorption feature at 589 nm and the predictions of previous models. For the atmospheric temperature-pressure profile we derive, silicate and iron clouds reside at a pressure of several mbar in the planet's atmosphere. These clouds have significant vertical extent and optical depth due to our slant viewing geometry and lead to increased absorption in bands directly adjacent to the sodium line core. Using a non-LTE sodium ionization model that includes photoionization by stellar UV flux, collisional processes with H_2, and radiative recombination, we show that the ionization depth in the planet's atmosphere reaches ~1/2 mbar at the day/night terminator. Ionization leads to a slight weakening of the sodium feature. We present our baseline model, including ionization and clouds, which falls near the observational error bars. The sensitivity of our conclusions to the derived atmospheric temperature-pressure profile is discussed.

We present a synthesis of physical effects influencing the observed lightcurve of an extrasolar giant planet (EGP) transiting its host star. The synthesis includes a treatment of Rayleigh scattering, cloud scattering, refraction, and molecular absorption of starlight in the EGP atmosphere. Of these effects, molecular absorption dominates in determining the transit-derived radius $R$. Using a generic model for the atmosphere of EGP HD209458b, we perform a fit to the best available transit lightcurve data, and infer that this planet has a radius at a pressure of 1 bar, $R_1$, equal to 94430 km. We predict that $R$ will be a function of wavelength of observation, with a robust prediction of variations of $\pm 1$% at infrared wavelengths where H$_2$O opacity in the high EGP atmosphere dominates.

Straddling the traditional realms of the planets and the stars, objects below the edge of the main sequence have such unique properties, and are being discovered in such quantities, that one can rightly claim that a new field at the interface of planetary science and and astronomy is being born. In this review, we explore the essential elements of the theory of brown dwarfs and giant planets, as well as of the new spectroscopic classes L and T. To this end, we describe their evolution, spectra, atmospheric compositions, chemistry, physics, and nuclear phases and explain the basic systematics of substellar-mass objects across three orders of magnitude in both mass and age and a factor of 30 in effective temperature. Moreover, we discuss the distinctive features of those extrasolar giant planets that are irradiated by a central primary, in particular their reflection spectra, albedos, and transits. Aspects of the latest theory of Jupiter and Saturn are also presented. Throughout, we highlight the effects of condensates, clouds, molecular abundances, and molecular/atomic opacities in brown dwarf and giant planet atmospheres and summarize the resulting spectral diagnostics. Where possible, the theory is put in its current observational context.

The recent discovery that the close-in extrasolar giant planet, HD209458b, transits its star has provided a first-of-its-kind measurement of the planet's radius and mass. In addition, there is a provocative detection of the light reflected off of the giant planet, $\tau$ Boo b. Including the effects of stellar irradiation, we estimate the general behavior of radius/age trajectories for such planets and interpret the large measured radii of HD209458b and $\tau$ Boo b in that context. We find that HD209458b must be a hydrogen-rich gas giant. Furthermore, the large radius of close-in gas giant is not due to the thermal expansion of its atmosphere, but to the high residual entropy that remains throughout its bulk by dint of its early proximity to a luminous primary. The large stellar flux does not inflate the planet, but retards its otherwise inexorable contraction from a more extended configuration at birth. This implies either that such a planet was formed near its current orbital distance or that it migrated in from larger distances ($\geq$0.5 A.U.), no later than a few times $10^7$ years of birth.

We present evolutionary sequences for Jupiter and Saturn, based on new nongray model atmospheres, which take into account the evolution of the solar luminosity and partitioning of dense components to deeper layers. The results are used to set limits on the extent to which possible interior phase separation of hydrogen and helium may have progressed in the two planets. When combined with static models constrained by the gravity field, our evolutionary calculations constrain the helium mass fraction in Jupiter to be between 0.20 and 0.27, relative to total hydrogen and helium. This is in agreement with the Galileo determination. The helium mass fraction in Saturn's atmosphere lies between 0.11 and 0.25, higher than the Voyager determination. Based on the discrepancy between the Galileo and Voyager results for Jupiter, and our models, we predict that Cassini measurements will yield a higher atmospheric helium mass fraction for Saturn relative to the Voyager value.

We summarize and extend recent work on the theory of extrasolar giant planets (EGPs) and brown dwarfs, paying particular attention to Gliese 229 B, the albedos of EGPs, the compositions of substellar atmospheres, the connections with the giant planets in the solar system, cloud physics, and non-gray spectral synthesis. The role of condensates in altering the optical spectrum of Gliese 229 B is explored, as are the systematics of the reflection spectra from extrasolar giant planets near their primaries. In addition, we discuss the role of convection and disequilibrium chemistry in explaining the anomalous detection of CO in Gliese 229 B. Throughout, we highlight the distinctive chemistry that defines this new class of objects and set goals for future study.

Giant planets in circumstellar disks can migrate inward from their initial (formation) positions. Radial migration is caused by inward torques between the planet and the disk; by outward torques between the planet and the spinning star; and by outward torques due to Roche lobe overflow and consequent mass loss from the planet. We present self-consistent numerical considerations of the problem of migrating giant planets. Summing torques on planets for various physical parameters, we find that Jupiter-mass planets can stably arrive and survive at small heliocentric distances, thus reproducing observed properties of some of the recently discovered extra-solar planets. Inward migration timescales can be approximately equal to or less than disk lifetimes and star spindown timescales. Therefore, the range of fates of massive planets is broad, and generally comprises three classes: (I) planets which migrate inward too rapidly and lose all their mass; (II) planets which migrate inward, lose some but not all of their mass, and survive in very small orbits; and (III) planets which do not lose any mass. Some planets in Class III do not migrate very far from their formation locations. Our results show that there is a wide range of possible fates for Jupiter-mass planets for both final heliocentric distance and final mass.

Since October of 1995, employing Doppler spectroscopy astronomers have discovered as many as 20 giant planets and brown dwarfs, including companions to $\tau$ Boo, 51 Peg, $\upsilon$ And, 55 Cnc, $\rho$ CrB, 70 Vir, 16 Cyg B, and 47 UMa. These discoveries have excited the planetary science community, astronomers, and the public at large. Within hours of the announcement of the planet 51 Peg b, the first direct detection of an unimpeachable brown dwarf, Gl229 B, was also announced. Gl229 B is a watershed since it has methane spectral features and a surface temperature below 1000 Kelvin, characteristics unique to objects with substellar masses. During the last two years, building upon our previous experience in the modeling of brown dwarfs and M stars, we published theoretical studies of the evolution and spectra of extrasolar giant planets. We have recently upgraded our capabilities and now generate non-gray spectral and color models of both giant planets and brown dwarfs. This theory will soon encompass objects whose effective temperatures range from 100 K to 4000 K and whose masses span three orders of magnitude. The evolutionary, spectral, and color calculations upon which we have embarked are in direct support of the searches now being planned in earnest with the HST (WFPC2, NICMOS), the IRTF, the MMT 6.5-meter upgrade, the LBT, Keck's I and II, ISO, UKIRT, NGST, the VLT, COROT, DENIS, 2MASS, and SIRTF.

Using the helium abundance measured by Galileo in the atmosphere of Jupiter and interior models reproducing the observed external gravitational field, we derive new constraints on the composition and structure of the planet. We conclude that, except for helium which must be more abundant in the metallic interior than in the molecular envelope, Jupiter could be homogeneous (no core) or could have a central dense core up to 12 Earth masses. The mass fraction of heavy elements is less than 7.5 times the solar value in the metallic envelope and between 1 and 7.2 times solar in the molecular envelope. The total amount of elements other than hydrogen and helium in the planet is between 11 and 45 Earth masses.

We have developed a new non-gray theory of the evolution, spectra, and colors of extrasolar giant planets (EGPs) and brown dwarfs that reveals their exotic nature and uniqueness. We have discovered that the fluxes of such objects for T$_{eff}$s from 1300 K to 100 K can be spectacularly higher in the near infrared bands than black body values and that their infrared colors are anomalously blue. As a consequence, EGPs and brown dwarfs reside in hitherto unoccupied realms of the H--R diagram and may be more easily found with current and planned telescopes than previously imagined.

We present the results of a new series of non-gray calculations of the atmospheres, spectra, colors, and evolution of extrasolar giant planets (EGPs) and brown dwarfs for effective temperatures below 1300 K. This theory encompasses most of the mass/age parameter space occupied by substellar objects and is the first spectral study down to 100 K. These calculations are in aid of the multitude of searches being conducted or planned around the world for giant planets and brown dwarfs and reveal the exotic nature of the class. Generically, absorption by H_2 at longer wavelengths and H_2O opacity windows at shorter wavelengths conspire to redistribute flux blueward. Below 1200 K, methane is the dominant carbon bearing molecule and is a universal diagnostic feature of EGP and brown dwarf spectra. We find that the primary bands in which to search are $Z$ (\sim 1.05 \mic), $J$ (\sim 1.2 \mic), $H$ (\sim 1.6 \mic), $K$ (\sim 2.2 \mic), $M$ (\sim 5 \mic), and $N$ (\sim 10 \mic), that enhancements of the emergent flux over blackbody values, in particular in the near infrared, can be by many orders of magnitude, and that the infrared colors of EGPs and brown dwarfs are much bluer than previously believed. In particular, relative to $J$ and $H$, the $K$ band flux is reduced by CH_4 and H_2 absorption. Furthermore, we derive that for T_{eff}s below 1200 K most or all true metals are sequestered below the photosphere, that an interior radiative zone is a generic feature of substellar objects, and that clouds of H_2O and NH_3 are formed for T_{eff}s below \sim 400 K and \sim 200 K, respectively. This study is done for solar-metallicity objects in isolation and does not include the effects of stellar insolation. Nevertheless, it is an attempt to bridge the gap between the planetary and stellar realms and to develop a non-gray theory of objects from 0.3 \mj (``saturn'') to 70 \mj ($\sim$0.07 \mo). We find that the detection ranges for brown dwarf/EGP discovery of both ground-- and space-based telescopes are larger than previously estimated.

With the discovery of the companions of 51 Peg, 55 Cnc, $\tau$ Boo, $\upsilon$ And, 70 Vir, 47 UMa, and Gl229, evolutionary and spectral models of gas giants and/or brown dwarfs with masses from 0.3 through 60 times that of Jupiter assume a new and central role in the emerging field of extrasolar planetary studies. In this contribution, we describe the structural, spectral, and evolutionary characteristics of such exotic objects, as determined by our recent theoretical calculations. These calculations can be used to establish direct search strategies via SIRTF, ISO, and HST (NICMOS), and via various ground-based adaptive optics and interferometric platforms planned for the near future.

Electron-degenerate, pressure-ionized hydrogen (usually referred to as metallic hydrogen) is the principal constituent of brown dwarfs, the long-sought objects which lie in the mass range between the lowest-mass stars (about eighty times the mass of Jupiter) and the giant planets. The thermodynamics and transport properties of metallic hydrogen are important for understanding the properties of these objects, which, unlike stars, continually and slowly cool from initial nondegenerate (gaseous) states. Within the last year, a brown dwarf (Gliese 229 B) has been detected and its spectrum observed and analyzed, and several examples of extrasolar giant planets have been discovered. The brown dwarf appears to have a mass of about forty to fifty Jupiter masses and is now too cool to be fusing hydrogen or deuterium, although we predict that it will have consumed all of its primordial deuterium. This paper reviews the current understanding of the interrelationship between its interior properties and its observed spectrum, and also discusses the current status of research on the structure of giant planets, both in our solar system and elsewhere.

The recent discoveries of giant planets around nearby stars have galvanized the planetary science community, astronomers, and the public at large. Since {\it direct} detection is now feasible, and is suggested by the recent acquisition of Gl229 B, it is crucial for the future of extrasolar planet searches that the fluxes, evolution, and physical structure of objects from Saturn's mass to 15 Juipter masses be theoretically investigated. We discuss our first attempts to explore the characteristics of extrasolar giant planets (EGPs), in aid of both NASA's and ESA's recent plans to search for such planets around nearby stars.

Theoretical spectra and evolutionary models that span the giant planet--brown dwarf continuum have been computed based on the recent discovery of the brown dwarf, Gliese 229 B. A flux enhancement in the 4--5 micron window is a universal feature from Jovian planets to brown dwarfs. We confirm the existence of methane and water in Gl 229 B's spectrum and find its mass to be 30 to 55 Jovian masses. Although these calculations focus on Gliese 229 B, they are also meant to guide future searches for extra-solar giant planets and brown dwarfs.

Using Doppler spectroscopy to detect the reflex motion of the nearby star, 51 Pegasi, Mayor \& Queloz (1995) claim to have discovered a giant planet in a 0.05-A.U., 4.23-day orbit. They estimate its mass to be in the range 0.5 \mj to 2 \mj, but are not able to determine its nature or origin. Including the effects of the severe stellar insolation implied, we extend the theory of giant planets we have recently developed to encompass those at very small orbital distances. Our calculations can be used to help formulate search strategies for luminous planets in tight orbits around other nearby stars. We calculate the radii and luminosities of such giants for a variety of compositions (H/He, He, H$_2$O, and olivine), evolutionary tracks for solar- composition gas giants, and the geometry of the Hayashi forbidden zone in the gas-giant mass regime. We show that such planets are stable and estimate the magnitude of classical Jeans evaporation and of photodissociation and loss due to EUV radiation. Even over the lifetime of the primary, the present companion would not have lost a large fraction of its mass. In addition, we demonstrate that for the mass range quoted, such planets are well within their Roche lobes. We show that the strong composition-dependence of the model radii and distinctive spectral signatures provide clear diagnostics that might reveal 51 Peg B's nature, should interferometric or adaptive-optics techniques ever succeed in photometrically separating planet from star.

We present a broad suite of models of extrasolar giant planets (EGP's), ranging in mass from 0.3 to 15 Jupiter masses. The models predict luminosity (both reflected and emitted) as a function of age, mass, deuterium abundance and distance from parent stars of various spectral type. We also explore the effects of helium mass fraction, rotation rate and the presence of a rock-ice core. The models incorporate the most accurate available equation of state for the interior, including a new theory for the enhancement of deuterium fusion by electron screening which is potentially important in these low mass objects. The results of our calculations reveal the enormous sensitivity of EGP's to the presence of the parent star, particularly for G and earlier spectral types. They also show a strong sensitivity of the flux contrast in the mid-infrared between parent star and EGP to the mass and age of the EGP's. We interpret our results in terms of search strategies for ground- and space-based observatories in place or anticipated in the near future.

Interest among astronomers in the detection of extra-solar planets is accelerating with the growing realization that it may soon be technically feasible. The ongoing renaissance in telescope construction and the anticipated launches of new space platforms are encouraging many scientists to review and improve the means by which planets can be discovered. The direct detection of the light from a distant planet would be the most compelling means of discovery and to gauge the feasibility of various search strategies, astronomers have traditionally used the current Jupiter as a benchmark planet. However, in principle, extra-solar giant planets (EGPs) can have a wide range of masses and, hence, can be significantly brighter than Jupiter. Furthermore, the maximum mass a planet can have is not known a priori, and observations will be needed to determine it. We predict the optical and infrared fluxes of EGPs with masses from 0.3 through 15 Jupiter masses and ages from 10$^7$ through $5 \times 10^{9}$ years that searches in the next few years may reveal.