{'fraction_non_alphanumeric': 0.03400871850968996, 'fraction_numerical': 0.029518849287335305, 'mean_word_length': 4.441542288557214, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 8, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Influence of Atmospheric Dynamics on the Infrared Spectra and Light Curves of Hot Jupiters

We explore the infrared spectrum of a three-dimensional dynamical model of planet HD209458b as a function of orbital phase. The dynamical model predicts day-side atmospheric pressure-temperature profiles that are much more isothermal at pressures less than 1 bar than one-dimensional radiative-convective models have found. The resulting day-side thermal spectra are very similar to a blackbody, and only weak water absorption features are seen at short wavelengths. The dayside emission is consequently in significantly better agreement with ground-based and space-based secondary eclipse data than any previous models, which predict strong flux peaks and deep absorption features. At other orbital phases, absorption due to carbon monoxide and methane is also predicted. We compute the spectra under two treatments of atmospheric chemistry: one uses the predictions of equilibrium chemistry, and the other uses non-equilibrium chemistry, which ties the timescales of methane and carbon monoxide chemistry to dynamical timescales. As a function of orbital phase, we predict planet-to-star flux ratios for standard infrared bands and all Spitzer Space Telescope bands. In Spitzer bands, we predict 2-fold to 15-fold variation in planetary flux as a function of orbital phase with equilibrium chemistry, and 2-fold to 4-fold variation with non-equilibrium chemistry. Variation is generally more pronounced in bands from 3-10 $\mu$m than at longer wavelengths. The orbital phase of maximum thermal emission in infrared bands is 15--45 orbital degrees before the time of secondary eclipse. Changes in flux as a function of orbital phase for HD209458b should be observable with Spitzer, given the previously acheived observational error bars.


INTRODUCTION
Astronomers and planetary scientists are just beginning to understand the atmospheres of the short period giant planets known as "hot Jupiters" or "Pegasi planets." A key subset of these planets are those that transit the disk of their parent stars, which make them well-suited for follow-up studies. The most well studied of these transiting hot Jupiters is the first to be discovered, HD 209458b (Charbonneau et al. 2000;Henry et al. 2000), a 0.69 M J planet that orbits its Sun-like parent star at a distance of 0.045 AU.
Currently, considerable work on hot Jupiters is occuring on both the observational and theoretical fronts. In the past few years, several groups have computed dynamical atmosphere models for HD 209458b in an effort to understand the structure, winds, and temperature contrasts of the planet's atmosphere (Showman & Guillot 2002;Cho et al. 2003;Cooper & Showman 2005;Burkert et al. 2005;. If the planet has been tidally de-spun and become locked to its parent star, dynamical models are surely needed to understand the extent to which absorbed stellar energy is transported onto the planet's permanent night side. With the launch of the Spitzer Space Telescope, there is now a platform that is wellsuited for observations of the thermal emission from hot Jupiter planets. Spitzer observations spanning the time of the planet's secondary eclipse (when the planet passes behind its parent star) have been published for HD 209458b at 24 µm , TrES-1 at 4.5 and 8.0 µm (Charbonneau et al. 2005), and HD 189733b at 16µm (Deming et al. 2006a). In all cases, the observed quantity is the planet-to-star flux ratio in Spitzer InfaRed Array Camera (IRAC), InfraRed Spectrograph (IRS), or Multiband Imaging Photometer for Spitzer (MIPS) bands. The dual TrES-1 observations are especially interesting because they allow for a determination of the planet's mid-infrared spectral slope. Several efforts have also been made from the ground to observe the secondary eclipse of HD 209458b. Although no detections have been made, some important, occasionallly overlooked upper limits at K and L band have been obtained. These include Richardson et al. (2003a) around 2.3 µm, Snellen (2005) in K band, and Deming et al. (2006b) at L band. These ground-based observations constrain the flux emitted by the planet in spectral bands where water vapor opacity is expected to be minimal; therefore, emitted flux should be high.
This influx of data has spurred a new generation of radiativeconvective equilibrium models, whose resulting infrared spectra can be compared with data (Fortney et al. 2005;Burrows et al. 2005;Seager et al. 2005;Barman et al. 2005;Fortney et al. 2006;Burrows et al. 2006). See Marley et al. (2006) and Charbonneau et al. (2006) for reviews. The majority of these models are one-dimensional. Authors weight the incident stellar flux by 1/4, to simulate planet-wide average conditions, or by 1/2, to simulate day-side average conditions (with a cold night side). Barman et al. (2005) have investigated twodimensional models with axial symmetry around the planet's substellar-antistellar axis and computed infrared spectra as a function of orbital phase. Iro et al. (2005) have extended one-dimensional models by adding heat transport due to a simple parametrization of winds to generate longitude-dependent temperature maps, but they did not compute disk averaged spectra for these models. Very recently Burrows et al. (2006) have also investigated spectra and light curves of planets with various day-night effective temperature differences, assuming 1D profiles for each hemisphere. These one-and two-dimensional radiative-convective equilibrium models have had some success in matching Spitzer observations, but Seager et al. (2005) and Deming et al. (2006b) have shown that ground-based data for HD 209458b do not indicate promiment flux peaks at 2.3 and 3.8 µm, which solar composition models predict.
The various dynamical models for HD 209458b (Showman & Guillot 2002;Cho et al. 2003;Cooper & Showman 2005;Burkert et al. 2005; are quite varied in their treatment of the planet's atmosphere. We will not review them here, as that is not the focus of this paper. In general, temperature contrasts in the visible atmosphere are expected to be somewhere between 300 and 1000 K. What has been somewhat lacking for these dynamical models are clear observational signatures, which would in principle be testable with Spitzer or other telescopes. The purpose of this paper is to remedy that situation. Here we generate infrared spectra and light curves as a function of orbital phase for the  (hereafter: CS06) dynamical simulation. We present the first spectra generated for three-dimensional models of the atmosphere of HD 209458b.

METHODS
Here we take the first step towards understanding the effects of atmospheric dynamics on the infrared spectra of hot Jupiters. A consistent treatment of coupled atmospheric dynamics, nonequilibrium chemistry, and radiative transfer would be a considerable task. In a coupled scheme, given a three-dimensional P-T grid at a given time step, with corresponding chemical mixing ratios, the radiative transfer scheme would solve for the upward and downward fluxes in each layer. These fluxes would be wavelength dependent and would differ from layer to layer. The thermodynamical heating/cooling rate, which is the vertical divergence of the net flux, would then be calculated. The dynamics scheme would then use this heating rate, together with the velocities and P-T profiles in the grid at the previous time step, to a calculate the chemical abundances, velocities, and P-T profiles on the grid at the new time step. The process steps forward in time, as the radiative transfer solver again finds the new heating/cooling rates. The emergent spectrum of the planet could be found at any stage. In our work presented here, we performed a simplified calculation that contains some aspects of what will eventually be included in a fully consistent treatment.

Dynamical Model
Our input pressure-temperature (P-T) map is from CS06. As the dynamical simulations are described in depth in Cooper & Showman (2005) and CS06, we will only give an overview here. The CS06 model employs the ARIES/GEOS Dynamical Core, version 2 (AGDC2; Suarez & Takacs 1995). The AGDC2 solves the primitive equations of dynamical meteorology, which are the foundation of numerous climate and numerical weather prediction models (Holton 1992;Kalnay 2002). The primitive equations simplify the Navier-Stokes equations of fluid mechanics by assuming hydrostatic balance of each ver-tical column of atmosphere. Forcing is due to incident flux from the parent star, through a Newtonian radiative process described in CS06. The CS06 model is forced from the one-dimensional radiative-convective equilibrium profile of Iro et al. (2005), which assumes globally averaged planetary conditions. For their simulations, CS06 take the top layer of the model to be 1 mbar. The model atmosphere spans ∼15 pressure scale heights between the input top layer and the bottom boundary at 3 kbar. CS06 use 40 layers evenly spaced in log pressure. A P-T profile is generated at locations evenly spaced in longitude (in 5 degree increments) and latitude (in 4 degree increments). The 72 longitude and 44 latitude points create 3168 P-T profiles. In Figure 1 we show the CS06 grid at 3 pressure levels, near ∼2, 20, and 200 mbar. Previous work has shown these levels likely bracket the pressures of interest for forming mid-infrared spectra (Fortney et al. 2005). All three panels of Figure 1 use the same brightness scale for easy comparison between pressure levels.
At the 2 mbar level, where radiative time constants are short (Iro et al. 2005;Seager et al. 2005), the atmosphere responds quickly to incident radiation. Winds, though reaching a speed of up to 8 km s −1 , are not fast enough to lead to significant deviations from a static atmosphere, which implies that the hottest regions remain at the substellar point. The arrows indicate a wind pattern that attempts to carry energy radially away from the hot spot. At this pressure, the atmosphere somewhat resembles the two-dimensional radially symmetric static atmosphere of Barman et al. (2005), who found a hot spot at the substellar point and a uniformly decreasing temperature as radial distance from this point increased. The night side appears nearly uniform, and colder.
At the 25 mbar level, it is clear that a west-to-east circulation pattern has emerged at the equator, and the center of the planet's warm region has been blown downstream by ∼35 degrees. The wind from the west dominates over the predominantly radially outward wind seen at the 2 mbar level. Day-night temperature contrasts are not as large at this pressure as they are at 2 mbar.
At the 220 mbar level, the center of the hot spot has been blown downstream by ∼60 degrees. This jet extends from the equator to the mid-latitudes; the gas in the jet is warmer than gas to the north or south. The radiative time constants become longer the deeper one goes into the atmosphere (Iro et al. 2005). Hence, winds are better able to redistribute energy, leading to weaker temperature contrasts, which cannot simply be characterized as "day-night."

Logistics & Radiative Transfer
Unlike other published models, we stress that the spectra generated here are from a dynamical atmosphere model that is not in radiative-convective equilibrium. Each of the 3168 P-T profiles from CS06 have, without modification (aside from interpolation onto a different pressure grid), been run through our radiative transfer solver. No iteration is done to achieve radiative equilibrium. The equation of radiative transfer is solved with the two-stream source function technique described in Toon et al. (1989). This is the same infrared radiative transfer scheme used in Fortney et al. (2005), Fortney et al. (2006) and M. S. Marley et al. (in prep.). We ignore contributions due to scattered stellar photons, as discussed below.
At a given orbital position, the CS06 map in longitude is remapped into an apparent longitude (as seen from Earth), while the latitude remains unchanged. See Figure 2 for a diagram.
Here we ignore the 3.4 degrees that the orbit differs from being exactly edge-on (Brown et al. 2001). At a given orbital angle ϕ-for each patch of the planet-the cosine of the angle θ from the sub-observer point, µ, is calculated. This µ is consistently included when solving the radiative transfer, which means that the effects of limb darkening (or brightening) are automatically incorporated. We interpolate in a pressure-temperatureabundance grid (Lodders & Fegley 2002), such that any given point in the three-dimensional model has the atomic and molecular abundances that are consistent with that point's pressure and temperature.
At any given time, the one-half of the planet that is not visible is not included in the radiative transfer. The 1584 visible points, at which the emergent specific intensity (erg s −1 cm −2 Hz −1 sr −1 ) is calculated, are then weighted by the apparent visible area of their respective patches. These intensities are summed up to give the total emergent flux density (erg s −1 cm −2 Hz −1 ) from the planet. The spectra generated by the patch-by-patch version of the code were tested against spectra from one-dimensional gray atmospheres and our previously published hot Jupiter profiles.
Emergent spectra are calculated from 0.26 to 325 µm, but since we ignore the contribution due to scattered stellar flux, the spectra at the very shortest wavelengths have little meaning. However, for the radiative-equilibrium HD 209458b model published in Fortney et al. (2005), they found that scattered stellar flux is greater than thermal emission only at wavelengths less than ∼0.68 µm. We expect that a considerable amount of "visible" light that may eventually be seen from hot Jupiters is due to thermal emission. We note that by 1 µm thermal emission is 100 times greater than scattered flux. Here we will present spectra for wavelengths from 1 to 30 µm.
We note that the radiative transfer at every point is solved in the plane-parallel approximation. While this treatment is sufficient for our purposes, we wish to point out two drawbacks. The first is that, near the limb of the planet, we will tend to overestimate the path lengths of photons emerging from the atmosphere, as the curvature of the atmosphere is neglected. The second issue also occurs near the limb. We cannot treat photons whose path, in a completely correct treatment, would start in one column but emerge from an adjoining column. However, atmospheric properties in any two adjoining columns are in general quite similar. The former issue, that of the planeparallel approximation, is likely more important, and should be addressed at a later time, when data precision warrants it. Here we note that in our tests 95% of planetary flux emerges from within 75 degrees of the sub-observer point, such that these limb effects will have little effect on the disk-summed spectra and light curves that we present. When generating our model spectra, we use the elemental abundance data of Lodders (2003) and chemical equilibrium compositions computed with the CONDOR code, as described in Fegley & Lodders (1994), Lodders & Fegley (2002), and Lodders (2002). In §2.3 we will discuss deviations from equilibrium chemistry in the atmosphere of HD 209458b, as calculated by CS06. At this time, we ignore photochemistry, which has been shown by Visscher & Fegley (2005a) to be reasonable at P > 10 mbar. For the most part, the infrared spectra of hot Jupiters are sensitive to opacity at 10 P 200 mbar (Fortney et al. 2005), so equilibrium chemistry calculations are probably sufficient. As discussed in Fortney et al. (2006), we maintain a large and constantly updated opacity database, which is de-scribed in detail in R. S. Freedman & K. Lodders (in prep.). The pressure, temperature, composition, and wavelength-dependent opacity is tabulated beforehand using the correlated-k method (Goody et al. 1989) in 196 wavelength interval bins. The resulting spectra are therefore of low resolution. However, low resolution is suitable for the task at hand, as we are interested in band-integrated fluxes and the radiative transfer must be solved at 1584 locations on the planet at many (here, 36) orbital phases.
In Williams et al. (2006), which focused on examining asymmetrical secondary eclipse light curves caused by dynamical redistribution of flux, we investigated the effects of limb darkening for HD 209458b for the CS06 map. These effects are hard to disentangle from general brightness variations due to temperature differences generated by dynamics. Limb darkening in a particular wavelength would be manifested as a brightness temperature that decreases towards the limb relative to a brightness temperature map computed assuming normal incidence at every point. On the planet's day side, where the P-T profiles are somewhat isothermal, limb darkening is not expected to be significant. Indeed, only at angles greater than ∼80 degrees from the subsolar point was day-side limb darkening as large as 100-200 K calculated. As noted, due to our plane-parallel approximation, this is likely to be somewhat of an overestimation.

Clouds and Chemistry
Before calculating spectra for this dynamic atmosphere, it is worthwhile to step back and look at the atmospheric P-T profiles with an eye towards understanding the effects that clouds and chemistry may have on the emergent spectra. In Figure 3, we plot a random sampling of the 3168 P-T profiles and compare them to cloud condensation and chemical equilibrium boundaries. It should be noted that there is a greater density of profiles on the left side of the plot. These are profiles from the relatively uniform and cool night hemisphere. Also of note is that at pressures less than ∼200 mbar, the warmer profiles are fairly isothermal, with quite shallow temperature gradients.
The CS06 profiles are nearly everywhere cooler than the condensation curves of iron and Mg-silicates. These clouds will form, but for these profiles, cloud bases would lie deep in the atmosphere at pressures greater than 1 kbar, as this is the highest pressure at which the profiles cross the condensation curves of these species. The opacity from such deep clouds would have no effect on the spectrum of HD 209458b.
The temperatures of the CS06 simulation are computed as departures from the equilibrium temperature profile of Iro et al. (2005), as described in CS06. If a warmer base profile had been selected, such as Fortney et al. (2005) or Barman et al. (2005) (which are 100-300 K warmer at our pressures of interest), these curves could be shifted to the right by 100-300 degrees. Marley et al. (2006) provide a graphical comparison of profiles computed by Iro et al. (2005), Fortney et al. (2005), and Barman et al. (2005) under similar assumptions concerning redistribution of stellar flux and atmospheric abundances; they find differences of up to 300 K at 100 mbar. These differences can probably be attributed to different opacity databases, molecular abundances, radiative transfer methods, and perhaps incident stellar fluxes. For the Fortney et al. (2005) profile, the major heating species are H 2 O bands from 1 to 3 µm and neutral atomic Na and K, which absorb strongly in the optical. The major cooling species are CO, at 5 µm, and H 2 O in bands at 3 and 4-10 µm, although for a colder profile, such as CS06 night side profiles, CH 4 would also cool across these wavelengths. CS06 chose the Iro et al. (2005) profile because Iro et al. also computed atmospheric radiative time constants, which other authors have not done to date. Here, the computed night side profiles cross the condensation curve of Na 2 S at low pressure. Showman & Guillot (2002) and Iro et al. (2005) pointed out that the transit detection of weak neutral atomic Na absorption by Charbonneau et al. (2002) could be explained if a large fraction of Na was tied up in this condensate.
If the predictions of equilibrium chemistry hold for this atmosphere, then it would appear that at pressures up to 100 mbar, CH 4 would be the main carbon carrier on the planet's night-side, while CO would be the main carrier on the day side. This would lead to dramatically different spectra for these hemispheres, especially at near-infrared and mid-infrared wavelengths.
However, equilibrium chemical abundances would be expected only if chemical timescales are faster than any mixing timescales. The relative abundances of CO and CH 4 can be driven out of equilibrium if mixing timescales due to dynamical winds are faster than the timescale for relaxation to equilibrium of the (net) reaction (1) The time constant for chemical relaxation toward equilibrium is a strong function of temperature and pressure; it is short in the deep interior and extremely long in the observable atmosphere. Deviations from CH 4 /CO equilibrium are due to the long timescale for conversion of CO to CH 4 , and are well known in the atmospheres of the Jovian planets (Prinn & Barshay 1977;Fegley & Prinn 1985;Fegley & Lodders 1994;Bézard et al. 2002;Visscher & Fegley 2005b) and brown dwarfs (Fegley & Lodders 1996;Noll et al. 1997;Griffith & Yelle 1999;Saumon et al. 2000Saumon et al. , 2003Saumon et al. , 2006, where vertical mixing can be due to convection and/or eddy diffusion. CS06 follow timescales for CH 4 /CO chemistry taken from Yung et al. (1988) and find that vertical winds of 5-10 m sec −1 push this system out of equilibrium at pressures of interest for the formation of infrared spectra. They find that the quench level, where the mixing timescale and chemistry timescale are equal, is ∼1-10 bar, which is below the visible atmosphere. Above this quench level, the mole fractions of CO and CH 4 are essentially constant. They find that the CH 4 /CO ratio becomes homogenized at pressures less than 1 bar everywhere on the planet. For their nomical models this CH 4 /CO ratio is ∼0.014, meaning that the majority of carbon is indeed in CO. However, a non-negligible fraction of carbon remains in CH 4 on both the day and the night hemispheres. We note that CS06 ignore possible effects due to photochemistry on the CO and CH 4 mixing ratios. Given the few explorations into hot Jupiter carbon and oxygen photochemistry to date (Liang et al. 2003(Liang et al. , 2004Visscher & Fegley 2005a), it is unclear how important photochemistry will be in determining the abundances of these species at pressures of tens to hundreds of millibar. CS06 did explore other effects such at atmospheric metallicity and temperature. If the atmosphere is greater than [M/H]=0.0, or if the atmosphere is hotter, the CH 4 /CO ratio would be even smaller. See CS06 for additional discussion.
In our spectral calculations, we find that significant differences arise depending on our treatment of chemistry. We will therefore investigate the effects of a few chemistry cases, as explained below, and shown in Table 1. We label our equi-librium chemistry trial "Case 0," as it is the standard case. For "Case 1," we fix the CH 4 /CO ratio at 0.014 (as found by CS06) at all temperatures and pressures along the profiles. Consistently incorporating the increasing CH 4 /CO ratio at depth (P > 1 bar) would be difficult with previously tabulated opacities, and would have little to no effect on the emergent spectra. Recall that 23% of available oxygen is lost to the formation of Mg-silicate clouds (Lodders 2003), which in this model have cloud bases near 1 kbar. The remaining oxygen is almost entirely found in CO and H 2 O. As the CH 4 and CO abundances are fixed, for Case 1 we will also fix the abundance of H 2 O; the mixing ratios of CH 4 , CO and H 2 O are consistent with the amount of available oxygen at T∼ 1200 K and pressures of tens of millibars. We also briefly consider a "Case 2," in which the CH 4 /CO ratio has been further reduced, in an ad-hoc manner, which could be due to quenching of the abundances at a hotter temperature or the photochemical destruction of CH 4 . The CH 4 /CO ratio is reduced by nearly a factor of 500, to 3 × 10 −5 . The mixing ratios of CO and H 2 O are nearly the same as in Case 1, although the slightly increased CO abundance (due to the drop in the abundance of CH 4 and conservation of carbon atoms) uses up some oxygen at the expense of H 2 O. This case is valuable for comparison purposes because it shows the spectral effects of a negligible CH 4 abundance. In all cases, the mixing ratios of all other chemical species are given by equilibrium values. However, only CH 4 , CO, and H 2 O have a discernable impact on the spectra.
In Figure 4 we plot mixing ratios of CH 4 , CO, and H 2 O predicted from equilibrium chemistry along three P-T profiles in the atmosphere of HD 209458b. In gray are the mixing ratios for the one-dimensional profile of Fortney et al. (2005). The abundance of CH 4 is negligible. In thin black and thick black are the mixing ratios at the sub-stellar point and the anti-stellar point, respectively, of the CS06 simulation. The dominant carbon carrier is clearly CO, rather than CH 4 , except at P <100 mbar for the anti-stellar point. Most of the nightside has a chemistry profile similar to the anti-stellar point, and hence, strong CH 4 absorption will be seen. At the top of the plot, arrows indicate the assigned mixing ratios of Cases 1 and 2 for these three molecules.
In §3 we make quantitative comparisons to current groundbased and space-based infrared data. While we find that the agreement is good for the model presented here, we wish to stress that this study is exploratory. This is the first study that quantitatively explores the influence of atmospheric dynamics on the emergent spectra of a hot Jupiter atmosphere. The greatest uncertainty likely lies in the calculation of heating/cooling rates with a simple Newtonian cooling scheme in the dynamical model, as discussed in Cooper & Showman (2005) and CS06. Another issue is that CS06 calculate temperature deviations from a one dimensional radiative-convective profile published in Iro et al. (2005). However, our emergent spectrum is calculated using our radiative transfer solver, which is a different code, and we use different abundances (Lodders 2003) than these authors (Anders & Grevesse 1989). However, we find that when solving the radiative transfer for the Iro et al. profile we obtain a T eff that differs by only 1% from their value.
In their non-equilibrium chemistry work, CS06 chose the abundances of Lodders & Fegley (1998), while here we use (Lodders 2003), in order for the most natural comparison with our previous work (Fortney et al. 2005. Obviously no choice would be fully self consistent. However, the choice of el-emental abundances is certainly a smaller concern than the current debate concerning the correct chemical timescale for conversion of CO to CH 4 in planetary atmospheres (see Yung et al. 1988;Fegley & Lodders 1994;Griffith & Yelle 1999;Bézard et al. 2002;Visscher & Fegley 2005b). In addion, CS06 have shown that a 300 K decrease in temperature, for instance, leads to a 20-fold increase in CH 4 /CO (see §5.3), so temperature uncertainties likely swamp any abundance issues. We stress that our focus in the following sections is on highlighting the spectral differences between a model that accounts for dynamical redistribution of energy around the planet and one that does not. Indeed the magnitude of the spectral effect suggests that additional, more internally self consistent work, is clearly appropriate.

Spectra as a Function of Orbital Phase
We now turn to the predicted infrared spectrum of our dynamical model atmosphere. As has been shown by many authors since Seager & Sasselov (1998), the infrared spectra of hot Jupiters are believed to be carved predominantly by absorption by H 2 O, CO, and, if temperatures are cool enough, CH 4 . In general, absorption features of hot Jupiters are predicted to be shallower than brown dwarfs of similar effective temperature and abundances. This is predicated on hot Jupiters having a significantly shallower atmospheric temperature gradient, which is due to the intense external irradiation by the parent star. Figure 5 shows absorption cross sections per molecule for CH 4 , CO, and H 2 O from 1 to 30 µm. To avoid clutter on our spectral plots, we will not label absorption features, so referring to Figure 5 will be helpful.
Our first spectral calculation is for Case 0 chemistry, as a function of orbital phase. This is shown in Figure 6. We show the emitted spectrum at six orbital locations: zero (during the transit), 60, 120, 180 (during the secondary eclipse), 240, and 300 degrees, where the degrees are orbital degrees after transit. At zero degrees (the black spectrum), we see the full night side of the planet. The temperature gradient is fairly steep, leading to deep absorption features. Most carbon is in the form of CH 4 , leading to strong methane absorption in bands centered on 1.6, 2.2, 3.3, and 7.8 µm. Absorption due to H 2 O at many wavelengths is also strong. Absorption due to CO at 2.3 and 4.5 µm is quite weak, due to its small mixing ratio.
Sixty degrees later (the red spectrum), the hemisphere we see is 2/3 night and 1/3 day. The part of the day hemisphere we see has the hot west-to-east jet coming towards us. This will later be an interesting point of comparison with the 300 orbital degree spectrum (in magenta), which is also 2/3 night. Here all absorption features are muted because, due to the jet shown in Figure 1, more than ∼1/3 of the planet is showing the nearly isothermal profiles representative of the day side. These profiles lead to a blackbody-like spectrum. At 120 orbital degrees (in green), we are seeing 2/3 day side, and that part of the day side is the hottest, due to the wind from the west blowing the hot spot downstream. The absorption features have become very weak.
The peak in total infrared flux actually occurs 153 orbital degrees after transit. During the secondary eclipse, 27 degrees later, the planet is less luminous. During the secondary eclipse (in dark blue) the planetary spectrum is essentially featureless long-ward of 1.8 µm. The fully illuminated hemisphere possesses profiles that are nearly isothermal in our pressure range of interest. Intriguingly, this leads to a spectrum that is very similar to a blackbody. For comparison, we plot a dashed black curve that is the thermal emission of a 1330 K blackbody. It plots behind the 180 degree spectrum for all λ > 4 µm.
As we move to orbital phases after the secondary eclipse, we see that the spectra are not symmetric with their reflections on the other half of the orbit. The light blue spectrum is significantly lower than the green one, even though at both points we are seeing 2/3 day and 1/3 night. This is again due to the westto-east jet. At 240 orbital degrees, the hottest part of the day side is not in view. At 300 orbital degrees (in magenta) 1/3 of our view is the coolest part of the day side, along with 2/3 of the night side. The spectrum is quite similar to what was seen during the transit. The planet appears least luminous 31 orbital degrees before transit.

Spectra with Alternate Chemistry
Since CS06 find significant deviations from equilibrium carbon chemistry, it is useful to examine spectra with nonequilibrium chemistry. We now investigate the spectra using Cases 1 and 2. Figure 7 shows comparison spectra at two orbital phases: zero and 180 orbital degrees, which is during transit and secondary eclipse, respectively. On the day side, we see that the emergent spectrum is essentially exactly the same, independent of chemistry. Because the P-T profiles are nearly isothermal, essentially no information about the chemistry can be gleaned from the spectrum. For comparison with our previous work, in Figure 7, we plot with a dotted line the cloudless planet-wide average spectrum of HD 209458b from Fortney et al. (2005), which shows absorption due to H 2 O and CO and strong flux peaks between these absorption features.
We will now examine the small inset in Figure 7. This is a zoom-in of the secondary-eclipse spectra near 2.2 µm. The three thick horizontal bars connected by dashes make up an interesting observational one-sigma upper limit that was first reported in Richardson et al. (2003a). These data were analyzed again and presented in Seager et al. (2005). The measurement was a relative one; only the vertical distance between the central band and two side bands is important: they can be moved up or down as a group. Specifically, the flux in the central band cannot exceed the combined flux in the two adjacent bands by an amount greater than vertical distance shown (which is 0.45 × 10 −6 erg s −1 cm −2 Hz −1 ). Our dynamical day-side spectra, which are nearly featureless due to the nearly isothermal atmosphere, easily meet this constraint. We note that an upper limit obtained by Snellen (2005) in K band, which covers a similar wavelength range, has an error bar that is too large to distinguish among published models.
In contrast to the models presented here that include dynamics, published 1D models to date all predict large flux differences between 2.2 µm and the surrounding spectral regions when abundances are solar. The 1D models with solar abundances presented in Seager et al. (2005), for example, were not able to meet the 2.2-µm flux-difference constraint, because the H 2 O bands adjoining the 2.2-µm flux peak were too deep. Seager et al. (2005) suggested that if C/O > 1, then little H 2 O would exist in the atmosphere, and the observational constraint could be met. Here, we propose instead that dynamics produces a relatively isothermal dayside temperature, in which case the constraint can easily be met even with solar C/O. Nevertheless, the observational constraint is not yet firm enough to fully rule out the 1D radiative-equilibrium models at solar abundances. The Fortney et al. (2005) solar-abundance models, for example, marginally meet the current observational flux-difference constraint. (Fortney et al. 2005 find a flux difference of 0.43 × 10 −6 erg s −1 cm −2 Hz −1 .) The cause for the discrepancy between the 1D models in Fortney et al. (2005) and Seager et al. (2005) remains unclear but could occur if the P-T profiles in Seager et al. were steeper than those in Fortney et al., leading to deeper bands in the former study. More precise flux-difference observations like that of Richardson et al. (2003a) during secondary eclipse would help confirm or rule out the 1D solar-abundance models. Even a factor of two decrease in the observational upper limit between the "in" and "out" bands could rule out the solar-abundance Fortney et al. (2005) model, while supporting the isothermal-dayside models (e.g., CS06) and the C/O > 1 models with weak water bands. Alternatively, the detection of a small flux peak at 2.2 µm would lend support to the 1D models with C/O < 1, while constraining the temperature gradient on the planet's dayside.
The fact that both the solar-abundance CS06 models (presented here) and the radiative-equilibrium 1D models with C/O > 1 produce almost no flux difference between 2.2 µm and the surrounding continuum implies that the 2.2-µm band cannot be used to distinguish between these alternatives. Instead, additional constraints at other wavelengths will be necessary to discriminate between them. Measurements of flux differences surrounding CO (rather than H 2 O) bands could provide such a test. The radiative-equilibrium C/O > 1 models, although lacking water bands, would presumably have strong CO bands and would hence predict a strong flux difference between the center of a CO band and the surrounding continuum. On the other hand, because the dayside atmosphere is nearly isothermal in the CS06 models, these circulation-altered models would predict little flux difference between the CO band and the surrounding regions even in the presence of large quantities of CO. An observation of minimal flux difference across CO bands as well as H 2 O bands would support isothermal-dayside models like CS06 while arguing against the radiative-equilibrium 1D models with steep temperature gradients.
Although observations of TrES-1 may not be directly applicable to HD 209458b, as TrES-1 may have a T eff 300 K cooler, there is an important issue to note for our discussion here. The models of Fortney et al. (2005), Seager et al. (2005), and Barman et al. (2005) show a mid-infrared spectral slope that is bluer than observed from 4.5 to 8.0 µm by Charbonneau et al. (2005) for TrES-1, although the model of Burrows et al. (2005) appears to be consistent with this slope.  and Fortney et al. (2006) have discussed that models with an atmospheric temperature inversion would give infrared spectra with a redder spectral slope, due to molecular emission features, in better agreement with observations. The CS06 model of the day side of HD 209458b, which lacks the strong negative temperature gradient of these radiative equilibrium models, also leads to a redder mid-infrared spectral slope.
On the planet's night side, we see significant spectral differences between the three chemistry cases. In Cases 1 and 2, the CH 4 mixing ratio is constrained to a small abundance, weakening these absorption features. Absorption due to CH 4 , CO, and H 2 O is readily seen in Case 1. Case 2 shows essentially the same CO and H 2 O absorption, but CH 4 absorption is no longer seen. The transition, as a function of orbital phase from deep to essentially nonexistent absorption features in Cases 1 and 2, are similar to what was seen in Case 0. In the interest of conciseness, and because Case 2 is ad-hoc, hereafter we concentrate on Cases 0 and 1. The spectra for Cases 1 and 2 are essentially the same, except in regions of CH 4 absorption; we will highlight these differences when necessary. It is important to remember that in Case 1 and Case 2 chemistry, the mixing ratios of our principal absorbers, CH 4 , CO, and H 2 O are fixed. Therefore, changes in the spectra with orbital phase are only due to changes in the planetary P-T profiles on the visible disk, due to the rotation of the planet through its orbit.
One can integrate the spectrum of the planet's visible hemisphere over all wavelengths, as a function of orbital phase, to determine the apparent luminosity of the planet at all phases. Here we divide out 4πR 2 σ to calculate the apparent effective temperature (T eff ) of the visible hemisphere. This is plotted in Figure 8 for Cases 0 and 1. We can clearly see that in both cases, the time of maximum flux precedes the time of secondary eclipse by ∼27 degrees, or 6.3 hours. Since the spectra of the two cases overlap around the time of secondary eclipse, so do the plots of T eff . At other orbital phases, the Case 0 curve always plots lower. The largest effect is before the time of transit. The effect is tied to the CH 4 abundance. If CH 4 is able to attain a large mixing ratio, it leads to an atmosphere that has a higher opacity, meaning one can not see as deeply into the atmosphere. One then reaches an optical depth of ∼1 higher in the atmosphere, which is significantly cooler for night-side P-T profiles, leading to a lower T eff . The point of minimum planetary flux precedes the transit by 31 to 37 degrees, depending on the chemistry.
In their (very similar) previous dynamical model, Cooper & Showman (2005) predicted a time of maximum planetary flux of 60 degrees (or 14 hours) before secondary eclipse. The large timing difference between that work and this one is due almost entirely to the choice of "photospheric pressure" made in Cooper & Showman (2005). They chose a pressure of 220 mbar in their work, which is deeper than the "mean" photosphere that we find here. At higher pressure, the radiative timescales are longer, such that winds are able to blow the atmosphere's hottest point farther downstream. Cooper & Showman (2005) and  previously discussed how their prediction varied as a function of the chosen photospheric pressure.

Spitzer Bands
For the foreseeable future, the most precise data for understanding the atmospheres of hot Jupiters will come from the Spitzer Space Telescope. For HD 209458b, only an observation at 24 µm (the shortest wavelength MIPS band), has been published. It seems likely that observations in all four IRAC bands (3.6, 4.5, 5.8, and 8.0 µm), as well as IRS at 16 µm, will be obtained within the next year or two. As such, we have integrated the spectra of our planet models and a Kurucz (1993) model of star HD 209458 over the Spitzer bands in order to generate planet-to-star flux ratios as a function of orbital phase. These are plotted for Cases 0 and 1 in Figure 9. The stellar model fits the stellar parameters derived in Brown et al. (2001) and is the same model used in Fortney et al. (2005).
The behavior of the planet-to-star flux ratios is quite interesting. While the T eff of the planet was found to reach a maximum 27 degrees before secondary eclipse, the behavior in individual bands is more varied. For instance, in both chemistry cases, the planetary flux in the 24 µm band peaks only 15 degrees be-fore secondary eclipse. This is because the "photospheric pressure" is at a lower pressure in this band than the planet's "mean photospheric pressure." As previously discussed, at lower pressures, the radiative time constants are shorter, and the atmosphere is able to more quickly adjust to changes in incident flux. At higher pressure, winds are better able to blow the planet's hot spot downstream. One should keep in mind that the light curves generated are a function of the dynamical calculation and the radiative transfer. It is the radiative transfer calculation that determines how deeply into the atmosphere (and therefore, to what temperature) we see.
The 3.6 µm band peaks earliest, ∼27 degrees before transit, as the T eff does as well. This band shows a 15-fold variation in flux (peak to trough) as a function of orbital phase for Case 0 because it encompasses a significant CH 4 absorption feature that waxes and wanes. Since the abundances of CH 4 and CO are not free to vary in Case 1, the flux ratios in this case do not show the large amplitudes found in some bands in Case 0. The dotted lines in the Case 1 panels are for Case 2 in the 3.6 and 8.0 µm bands, where CH 4 absorption occurs. The flux variation in these bands is even further reduced as CH 4 absorption is not seen. (See Figure 7.) At 24 µm, our model is 1.6 sigma higher than the secondary eclipse data point published by Deming et al. (2005), indicating that the planet may be dimmer at 180 orbital degrees than we predict. In all bands, differences in chemistry between the two cases have essentially no effect on the timing of the maxima in planetary flux. However, since the night side is much more sensitive to chemistry, the minima can vary by as much as 20 degrees in bands that are sensitive to CH 4 absorption.

Standard Infrared Bands
The results of ground-based observations of flux from hot Jupiters have been mixed. All searches for visible light have yielded only upper limits (Charbonneau et al. 1999;Collier Cameron et al. 2002;Leigh et al. 2003), which have ruled out some models with bright reflecting clouds. In the near infrared, specifically for HD 209458b, Richardson et al. (2003b,a) have constrained molecular bands of CH 4 and H 2 O. The constraint on emission at 2.3 µm between H 2 O absorption features was shown in Figure 7. The predicted planet-to-star flux ratio really does not become favorable until wavelengths longer than 3 µm, which may continue to challenge observers.
In Figure 10, we plot planet-to-star flux ratios in the Mauna Kea Observatory (MKO) H (∼1.6 µm), K (∼2.2 µm), L ′ (∼3.8 µm), and M ′ (∼4.7 µm) bands. Since the wavelength ranges of the L ′ and M ′ bands have significant overlap with the IRAC 3.6 and 4.5 µm bands, the predicted ratios are quite similar. In dark blue, we show the L ′ band upper limit of -0.0007±0.0014 from Deming et al. (2006b). The model is just above the 1 sigma error bar. 5 This is a better match than is attained with one dimensional radiative equilibrium models (Deming et al. 2006b), which predict a flux peak just short of 4 µm, as shown in Figure 7. The Fortney et al. (2005) HD 209458b radiative equilibrium model, which has somewhat muted flux peaks compared with similar models by other authors, has a planet-to-star flux ratio of 0.00114 in L ′ band. In dotted blue in the Case 1 panels is our L ′ band (which encompasses CH 4 absorption) prediction for Case 2 chemistry. Looking at shorter wavelengths, in the H and K bands, the planet is predicted to be dimmer, while the star is brighter, leading to low flux ratios. The peak emission in these bands does occur earlier than in the Spitzer bands. For instance, the H band peak is 42 orbital degrees before secondary eclipse.

Effective Temperature
For the CS06 dynamical model of the atmosphere of HD 209458b, we find that the apparent T eff of the visible hemisphere is strongly variable, with a maximum of 1390 K and a minimum of 915 K for equilibrium chemistry and 1025 K for (the probably more realistic case of) disequilibrium chemistry. This leads to a luminosity of the visible planetary hemisphere that varies by factors of 5.3 and 3.4, respectively for these two cases. For the one dimensional Iro et al. (2005) planet-wide profile, we derive a T eff = 1325 K, 1% less than found by the authors. The day-side T eff for the dynamical model is not significantly larger than this planet-wide T eff . As an additional comparison, we calculate the mean of the planetary luminosity over the entire orbit and then convert to T eff to find a mean T eff for Cases 0 and 1. This gives T eff s of 1195 K and 1227 K, respectively, showing that the CS06 model, at the pressures levels that radiate to space, is as a whole somewhat colder than the one dimensial profile of Iro et al. (2005). This is likely a consequence of the radiative forcing scheme employed in CS06, which will be reinvestigated when models that consistently couple radiative transfer and dynamics are developed.
The ∼400 K T eff contrasts found are predominantly due to changes in the temperature structure of the hemisphere that is visible as function of orbital phase. In addition, atmospheric opacity makes an important contribution when the CH 4 /CO ratio is free to vary with pressure and temperature. When the CH 4 mixing ratio is below that predicted by equilibrium chemistry, this leads to a lower opacity atmosphere, for a given P-T profile. While the few individual bands of CO are somewhat stronger than those of CH 4 , CH 4 absorption across the planet's broad the 2-10 µm flux peak dominates over the two bands of CO. One is able to see more deeply into a CH 4 -depleted atmosphere, leading to a higher T eff .

HD 209458b Infrared Data
Currently, the only published secondary eclipse datum from Spitzer for HD 209458b is the 24 µm detection of Deming et al. (2005). The models presented here have a planet-to-star flux ratio during secondary eclipse that is 1.6 sigma larger than this observational data point. Together with our excellent agreement with the ground-based data at 2.3 µm from Richardson et al. (2003a) and Seager et al. (2005), and our 1.1 sigma difference with the 3.8 µm data from Deming et al. (2006b), we regard this as excellent agreement-significantly better than has been previously obtained with one dimensional radiative equilibrium models.
It is interesting to discuss a few issues that arise if the flux ratios are indeed 10-25% less than we calculate here, as indicated by the L ′ and 24 µm band data. For instance, perhaps day-night temperature contrasts in the atmosphere are not as large as predicted by CS06, leading to smaller deviations from a "planet-wide" T eff ∼ 1200K. This may involve radiative time constants that are significantly longer than predicted by Iro et al. (2005), winds faster than predicted by CS06, or both. It would be worthwhile for other groups that possess hot Jupiter radiative transfer codes to compute radiative time constants at these temperatures and pressures. This is an area that we will pursue in the near future.
Another possibility for a smaller planet-to-star flux ratio during secondary eclipse would be if the planet had a larger Bond albedo than calculations currently indicate. This would mean less absorbed stellar flux and correspondingly lower temperatures everywhere on the planet. Cooler temperatures everywhere on the planet would lead to chemical abundances that differ from our previous cases. More CH 4 would form, at the expense of CO, which would also lead to a slightly higher mixing ratio for H 2 O, which shares oxygen with CO.
Hot Jupiters are believed to have very low Bond albedos-on the order of 90% or more of incident stellar flux is expected to be absorbed. In Fortney et al. (2005), we found that our onedimensional model atmosphere for HD 209458b scattered only 8% of incident flux. For TrES-1, this was 6%. To date, there is at least a hint that the Bond albedo of TrES-1 may have been underestimated. Charbonneau et al. (2005), under the assumption that the planet emits as a blackbody, determined a Bond albedo of 0.31 ± 0.14 from their Spitzer IRAC observations at 4.5 and 8.0 µm. Perhaps hot Jupiters are not quite as hot as had been previously thought. However, this determination should be considered very preliminary at this time. Richardson et al. (2006) have recently observed the transit of HD 209458b at 24 µm, as well. From an observed transit depth of 0.0149 ± 0.0003, they determined the radius of the planet in this band to be 1.26 ±0.08 R J , which includes uncertainties in the stellar radius. Our model predicts a change in the apparent planet-to-star flux ratio of ∼0.00008 during the 20 orbital degrees of the Richardson et al. (2006) transit observations, ∼4 times smaller than their uncertainty, and hence too small to have been detected. As was seen in Figures 9 and 10, the change in planetary flux as a function of orbital phase is not as pronounced near the time of transit (and secondary eclipse) as it is at other phases.

Temperature Sensitivity
To illustrate the sensitivity of our results to temperature changes, we have computed spectra as a function of orbital phase, using equilibrium (Case 0) chemistry for two additional models. These are additional dynamical models described in CS06, in which the base P-T profile of Iro et al. (2005) has been increased or decreased by 300 K, with correspondingly warmer or colder night sides. The full dynamical simulations have been run again with these parameters. The resulting light curves, in Spitzer bands, are shown in Figures 11a, and c. For the "cold" (-300 K) model CH 4 is the dominant carbon carrier on the night hemisphere, and CO on the day hemisphere, leading to large flux variation in excess of that found for our nominal CS06 simulation, as shown in Figure 11. For the "hot" (+300 K) model, on both the day and night hemispheres, CO is the dominant carbon carrier. Methane absorption features are very weak on the night side, leading to flux variations in every band no larger than a factor of 2.7 from peak to trough. This model is somewhat similar to our earlier Case 2, but at warmer temperatures, as at all phases CO is the dominant carbon carrier. Fluxes are everywhere greater in the hot model than the nominal model, and everywhere less in the cold model than in the nominal model. This is is due to the differing atmospheric temperatures. The cold model best fits the Deming et al. (2005) datum at 24 µm. In addition, the phase of maximum and minimum flux in a given band can change by up to ∼10 orbital degrees between these simulations, due both to differing chemistry and atmospheric dynamics. As discussed in Cooper & Showman (2005) and CS06, while these models all possess similarly strong east-to-west jets, the dynamical atmospheres differ slightly in detail.
CS06 also examined non-equilibrium CH 4 /CO chemistry for these models. For the "cold" (-300 K) model, non-equilibrium chemistry was important, and the CH 4 /CO ratio at P < 1 bar became homogenized at 0.20 around the entire planet. This ratio was 0.014 for the nominal Case 1 decribed earlier. For the cold model equilibrium chemistry would predict a night side dominated by CH 4 and a day side dominated by CO. For the "hot" model (+300 K) both equilibrium and non-equilibrium chemistry predicts that CO is the dominant carbon carrier on both the day and side hemispheres. Figure 11d shows our computed flux ratios for non-equilibrium chemistry for the cold case. As was shown previously for our nominal case, flux variation is smaller with non-equilibrium chemistry, because the mixing ratios of CH 4 , CO, and H 2 O are the same on the night and day hemispheres. Again, because of the nearly isothermal temperature structure of the day-side, chemical abundances have little effect on the day-side planet-to-star flux ratios, which are nearly equal for these two chemistry cases.
These additional cases further highlight the importance of CH 4 /CO chemistry in the computation of infrared light curves. Infrared fluxes as a function of orbital phase are sensitive to the temperatures of the hemisphere facing the observer, as well as the abundances of important molecular absorbers. For planetary P-T profiles that cross important CH 4 /CO chemical boundaries, as most hot Jupiters surely do, whether or not these species are found in their equilibrium mixing ratios has a major impact on resulting infrared flux in Spitzer bands, especially on the night side.

The Future & Conclusions
From the size of the error bar from the Deming et al. (2005) 24 µm observations, it is clear that, could this instrument stability be sustained over the course of tens of hours of observations, the change in flux over time presented here could be detected. One-half of an orbital period for HD 209458b is 42 hours. If flux in the 8.0 µm band for HD 209458b is higher than models predict, as was the case for TrES-1 (Fortney et al. 2005;Barman et al. 2005), then this would be an attractive band as well, as the error bars should be smaller.
What might one hope to see with sustained observations? If the peak in infrared flux does occur only ∼25-30 orbital degrees before secondary eclipse, this would be difficulty to discern. However, the detection of any sort of ramp up in flux from the time of transit to secondary eclipse would give us important information on the day-night temperature contrast. The recently discovered transiting planet HD 189733b (Bouchy et al. 2005) would likely be an even more attractive target, as the planet-star flux ratios are likely twice as large Deming et al. 2006a), and the orbital period is ∼40% shorter. We predict secondary eclipse planet-to-star flux ratios for this system in Fortney et al. (2006); our calculation at 16 µm is a good match to the published observation of Deming et al. (2006a).
A clear prediction from our calculations here is that when one uses realistic non-equilibrium chemistry calculations, the change in planetary flux as a function of orbital phase is greatly reduced, relative to equilibrium chemistry calculations, because the atmosphere's composition is fixed. For the Spitzer bands, for Case 0, the maximum flux variation is in the 3.6 µm band, which varies by a factor of 15 from peak to trough. The minimum variation is a factor of 2.2, in the 16 µm band. For Case 1, this variation drops significantly, and the maximum variation is a factor of 3.7, in the 5.8 µm band, and the minimum is 2.0, again in the 16 µm band. We suggest that the 5.8 and 8.0 µm bands may be the best Spitzer bands in which to search for flux variations as a function of orbital phase, as these bands combine a high planet-to-star flux ratio and the sensitivity and stability of the IRAC detectors.
If the day-side thermal emission of hot Jupiters is similar to a blackbody, as we find here, problems arise with the notion that the emission can be used to characterize the atmospheric chemistry from secondary eclipse observations. Absorption features due to CH 4 , CO, and H 2 O would be nonexistent or extremely weak. Observations at other orbital phases would then take on additional importance.
Given the significant spectral differences between our model and radiative-equilibrium models, it is clear that more work in this area is certainly warranted. We note that blackbody-like hot Jupiter emission can simultaneously explain all observations to date. First is the secondary eclipse observation at 24 µm by Deming et al. (2005), a clear detection. Second is the very low flux ratio upper limit in L ′ band by Deming et al. (2006b). Third is the 2.3 µm relative flux observation of Richardson et al. (2003a), which most one-dimensional models cannot fit . Fourth is the set of observations of TrES-1 by Charbonneau et al. (2005), who found an infrared spectral slope from 4.5 to 8 µm that was redder than that found by most one-dimensional models (Fortney et al. 2005;Barman et al. 2005). Additional observations of these planets, especially in the Spitzer 3.6 µm band, which catches much of the predicted 4 µm flux peak, would strengthen or refute this argument and provide a critical test for the CS06 dynamical simulation. Deming et al. (2006b) posited, perhaps with a wink, that a blackbodylike spectrum would be more consistent with observations to date than any published hot Jupiter model. Due to the dynamically altered temperature structure of the atmosphere of HD 209458b, we find that this could be the reality. Note. -Mixing ratios are given for each chemistry case. Abundances from equilibrium chemistry calculations are used when "equil." is specified. The CS06 latitude/longitude grid covers the entire planet. The planet's night side is shaded gray. At every angle in the orbit (ϕ), the longitude is remapped to give the apparent longitude and latitude, as viewed by the distant observer. The cosine of the angle θ made between a patch of the planet and the sub-observer point is µ .  FIG. 4.-The mixing ratios of CH 4 , CO, and H 2 O, as function of pressure, along three P-T profiles in the atmosphere of HD 209458b. These profiles are the one-dimensional profile from Fortney et al. (2005), in gray, the subsolar point of the CS06 simulation, in thin black, and the anti-stellar point of the CS06 simulation, in thick black. Arrows at the top of the plot indicate the fixed mixing ratios of Case 1 and 2 chemistry, given in Table 1. FIG. 6.-Planetary emergent flux density (erg s −1 cm −2 Hz −1 ) vs. wavelength as a function of orbital phase for Case 0 equilibrium chemistry. The spectra are color-coded with the diagram in the upper right corner of the figure. All spectra are calculated at intervals of 60 orbital degrees. The black spectrum is for the night side of the planet, which is seen during transit. The red spectrum is 60 orbital degrees later. The dark blue spectrum shows the planet during secondary eclipse (180 degrees), when the fully illuminated hemisphere is visible. The magenta spectrum is 300 orbital degrees after transit. The dashed black curve is the flux of a 1330 K blackbody, which plots behind the dark blue curve at λ > 4µm. Normalized Spitzer band passes are shown in dotted lines at the bottom and standard H, K, L, and M bands are shown at the top.
FIG. 7.-Planetary emergent flux density (erg s −1 cm −2 Hz −1 ) vs. wavelength at two orbital phases for all three of our chemistry cases. The spectrum at 180 degrees, during the secondary eclipse, is essentially the same for all three cases. At zero degrees, during the transit, CH 4 is very abundant in the Case 0 trial, leading to deep absorption features centered at 1.6, 2.2, 3.3, and 7.8 µm. Absorption due to CO is muted in Case 0, but is strong in Cases 1 and 2. Water vapor absorption is strong in all cases. Inset: Zoom near 2.2 µm. The three thick horizontal lines connected by dashes are a constraint on the maximum relative height of the flux peak at 2.2 µm (middle band), compared to the average of the flux in the two adjacent bands, from Seager et al. (2005). The bars can be shifted together vertically. The dotted lines indicate predictions for Case 2 chemistry. At 180 orbital degrees is the Deming et al. (2005) secondary eclipse observation at 24 µm, with one-sigma error bars (0.00260±0.00046). This is shown in magenta, which is the color used for the 24 µm curve.  Deming et al. (2006b) secondary eclipse observation at L ′ band, with one-sigma error bars. The observation is -0.0007±0.0014, so the actual data point is well off of the bottom of the plot. This is shown in dark blue, which is the color used for the L ′ curve.  Figure 9a. Panel C shows the CS06 dynamical simulation with a 300 K temperature decrease. The right-most panel (D) is also the simulation with the 300 K temperature decrease, but uses non-equilibrium chemistry with a fixed CH 4 /CO ratio of 0.20. (See text.) For all panels, the Deming et al. (2005) datum at 24 µm is shown.