The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $\mu$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $\mu$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

The exoplanet WASP-12b is the prototype for the emerging class of ultra-hot, Jupiter-mass exoplanets. Past models have predicted---and near ultra-violet observations have shown---that this planet is losing mass. We present an analysis of two sets of 3.6 $\mu$m and 4.5 $\mu$m $\textit{Spitzer}$ phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5 $\mu$m, not at 3.6 $\mu$m, which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, $\lesssim$ 600 K gas. It is unclear why WASP-12b is the only ultra-hot Jupiter to exhibit this mass loss signature, but perhaps WASP-12b's orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of $6.4\sigma$ ($46.2^{\circ}$) at 3.6 $\mu$m.

In 2010 the first planet was discovered around star HD 34445. Recently, another five planets were announced orbiting the same star. It is a rather dense multi-planet system with some of its planets having separations of fractions of an au and minimum masses ranging from Neptune to sub-Jupiter ones. Given the number of planets and the various uncertainties in their masses and orbital elements, the HD 34445 planetary system is quite interesting as there is the potential for mean motion and secular resonances that could render the outcome of its dynamical evolution and fate an open question. In this paper we investigate the dynamical stability of the six planet system in order to check the validity of the orbital solution acquired. This is achieved by a series of numerical experiments, where the dynamical evolution of the system is tested on different timescales. We vary the orbital elements and masses of the system within the error ranges provided. We find that for a large area of the parameter space we can produce stable configurations and therefore conclude it is very likely that the HD 34445 planetary system is real. Some discussion about the potential habitability of the system is also done.

Heng & Tsai (2016) developed an analytic framework to obtain thermochemical-equilibrium abundances for H$_{2}$O, CO, CO$_2$, CH$_4$, C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, HCN, NH$_3$, and N$_2$ for a system with known temperature, pressure, and elemental abundances (hydrogen, carbon, nitrogen, and oxygen). However, the implementation of their approach can become numerically unstable under certain circumstances, leading to inaccurate solutions (e.g., ${\rm C/O} \ge 1$ atmospheres at low pressures). Building up on their approach, we identified the conditions that prompt inaccurate solutions, and developed a new framework to avoid them, providing a reliable implementation for arbitrary values of temperature (200 to $\sim$2000 K), pressure ($10^{-8}$ to $10^{3}$ bar), and CNO abundances ($10^{-3}$ to $\sim 10^{2}\times$ solar elemental abundances), for hydrogen-dominated atmospheres. The accuracy our analytic framework is better than 10% for the more abundant species that have mixing fractions larger than $10^{-10}$, whereas the accuracy is better than 50% for the less abundant species. Additionally, we added the equilibrium-abundance calculation of atomic and molecular hydrogen into the system, and explored the physical limitations of this approach. Efficient and reliable tools, such as this one, are highly valuable for atmospheric Bayesian studies, which need to evaluate a large number of models. We implemented our analytic framework into the \textsc{rate} Python open-source package, available at https://github.com/pcubillos/rate .

Previous full-lifetime simulations of single-star multi-planet systems across all phases of stellar evolution have predominately assumed coplanar or nearly-coplanar orbits. Here we assess the consequences of this assumption by removing it and exploring the effect of giant branch mass loss on the stability of two-planet systems with small to moderate non-Kozai (<40 degrees) relative inclinations. We run nearly 10^4 simulations over 14 Gyr for F-star, A-star and B-star planet hosts, incorporating main sequence stellar masses of 1.5, 2.0, 2.5, 3.0 and 5.0 solar masses, and initial planetary semimajor axis ratios that straddle their three-dimensional Hill stability limits. We find that the near-coplanar assumption can approximate well the stability frequencies and critical separations found for higher inclinations, except around strong mean-motion commensurabilities. Late instabilities -- after the star has become a white dwarf -- occur throughout the explored mutual inclination range. Consequently, non-Kozai mutual inclination should not be used as a predictive orbital proxy for determining which white dwarf multi-planet systems discovered by Gaia should represent high-priority follow-up targets for the detection of metal pollution and planetary debris discs.

The scenario and efficiency of water transport by icy asteroids and comets are still amongst the most important unresolved questions of planetary systems. A better understanding of cometary dynamics in extrasolar systems shall provide information about cometary reservoirs and give an insight into water transport especially to planets in the habitable zone. The detection of Proxima Centauri-b (PCb), which moves in the habitable zone of this system, triggered a debate whether or not this planet can be habitable. In this work, we focus on the stability of an additional planet in the system and on water transport by minor bodies. We perform numerous N-body simulations with PCb and an outer Oort-cloud like reservoir of comets. We investigate close encounters and collisions with the planet, which are important for the transport of water. Observers found hints for a second planet with a period longer than 60 days. Our dynamical studies show that two planets in this system are stable even for a more massive second planet (~12 Earth masses). Furthermore, we perform simulations including exocomets, a second planet, and the influence of the binary Alpha Centauri. The studies on the dynamics of exocomets reveal that the outer limit for water transport is around 200 au. In addition we show that water transport would be possible from a close-in planetesimal cloud (1-4 au). From our simulations, based on typical M-star protoplanetary disks, we estimate the water mass delivered to the planets up to 51 Earth oceans.

Observations of transiting extra-solar planets provide rich sources of data for probing the in-system environment. In the WASP-12 system, a broad depression in the usually-bright MgII h&k lines has been observed, in addition to atmospheric escape from the extremely hot Jupiter WASP-12b. It has been hypothesized that a translucent circumstellar cloud is formed by the outflow from the planet, causing the observed signatures. We perform 3D hydrodynamic simulations of the full system environment of WASP-12, injecting a planetary wind and stellar wind from their respective surfaces. We find that a torus of density high enough to account for the lack of MgII h&k line core emission in WASP-12 can be formed in approximately 13 years. We also perform synthetic observations of the Lyman-alpha spectrum at different points in the planet's orbit, which demonstrate that significant absorption occurs at all points in the orbit, not just during transits, as suggested by the observations.

The presence of giant planets influences potentially habitable worlds in numerous ways. Massive celestial neighbors can facilitate the formation of planetary cores and modify the influx of asteroids and comets towards Earth-analogs later on. Furthermore, giant planets can indirectly change the climate of terrestrial worlds by gravitationally altering their orbits. Investigating 147 well characterized exoplanetary systems known to date that host a main sequence star and a giant planet we show that the presence of 'giant neighbors' can reduce a terrestrial planet's chances to remain habitable, even if both planets have stable orbits. In a small fraction of systems, however, giant planets slightly increase the extent of habitable zones provided the terrestrial world has a high climate inertia. In providing constraints on where giant planets cease to affect the habitable zone size in a detrimental fashion, we identify prime targets in the search for habitable worlds.

Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source BART code, we investigate the difference between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer, HST, and JWST bandpasses, covering the wavelength range between 1 and 11 um. Using the data from different observing instruments, we present detailed comparisons between the temperature-pressure profiles recovered by BART and those from the 3D simulations. We calculate several averages of the 3D thermal structure and implement two temperature parameterizations to investigate different thermal profile shapes. To assess which part of the thermal structure is best constrained by the data, we generate contribution functions for both our theoretical model and each of our retrieved models. Our conclusions are strongly affected by the spectral resolution of the instruments included, their wavelength coverage, and the number of data points combined. We also see some limitations in each of the temperature parametrizations. The results show that our 1D retrieval is recovering a temperature and pressure profile that most closely matches the arithmetic average of the 3D thermal structure. When we use a higher resolution, more data points, and a parametrized temperature profile that allows more flexibility in the middle part of the atmosphere, we find a better match between the retrieved temperature and pressure profile and the arithmetic average.

Short-period planets exhibit day-night temperature contrasts of hundreds to thousands of degrees K. They also exhibit eastward hotspot offsets whereby the hottest region on the planet is east of the substellar point; this has been widely interpreted as advection of heat due to eastward winds. We present thermal phase observations of the hot Jupiter CoRoT-2b obtained with the IRAC instrument on the Spitzer Space Telescope. These measurements show the most robust detection to date of a westward hotspot offset of 23 $\pm$ 4 degrees, in contrast with the nine other planets with equivalent measurements. The peculiar infrared flux map of CoRoT-2b may result from westward winds due to non-synchronous rotation magnetic effects, or partial cloud coverage, that obscures the emergent flux from the planet's eastern hemisphere. Non-synchronous rotation and magnetic effects may also explain the planet's anomalously large radius. On the other hand, partial cloud coverage could explain the featureless dayside emission spectrum of the planet. If CoRoT-2b is not tidally locked, then it means that our understanding of star-planet tidal interaction is incomplete. If the westward offset is due to magnetic effects, our result represents an opportunity to study an exoplanet's magnetic field. If it has Eastern clouds, then it means that our understanding of large-scale circulation on tidally locked planets is incomplete.

Multi-band phase variations in principle allow us to infer the longitudinal temperature distributions of planets as a function of height in their atmospheres. For example, 3.6 micron emission originates from deeper layers of the atmosphere than 4.5 micron due to greater water vapor absorption at the longer wavelength. Since heat transport efficiency increases with pressure, we expect thermal phase curves at 3.6 micron to exhibit smaller amplitudes and greater phase offsets than at 4.5 micron; this trend is not observed. Of the seven hot Jupiters with full-orbit phase curves at 3.6 and 4.5 micron, all have greater phase amplitude at 3.6 micron than at 4.5 micron, while four of seven exhibit a greater phase offset at 3.6 micron. We use a 3D radiative-hydrodynamic model to calculate theoretical phase curves of HD 189733b, assuming thermo-chemical equilibrium. The model exhibits temperature, pressure, and wavelength dependent opacity, primarily driven by carbon chemistry: CO is energetically favored on the dayside, while CH4 is favored on the cooler nightside. Infrared opacity therefore changes by orders of magnitude between day and night, producing dramatic vertical shifts in the wavelength-specific photospheres, which would complicate eclipse or phase mapping with spectral data. The model predicts greater relative phase amplitude and greater phase offset at 3.6 micron than at 4.5 micron, in agreement with the data. Our model qualitatively explains the observed phase curves, but is in tension with current thermo-chemical kinetics models that predict zonally uniform atmospheric composition due to transport of CO from the hot regions of the atmosphere.

As the 3D spatial properties of exoplanet atmospheres are being observed in increasing detail by current and new generations of telescopes, the modelling of the 3D scattering effects of cloud forming atmospheres with inhomogeneous opacity structures becomes increasingly important to interpret observational data. We model the scattering and emission properties of a simulated cloud forming, inhomogeneous opacity, hot Jupiter atmosphere of HD 189733b. We compare our results to available HST and Spitzer data and quantify the effects of 3D multiple-scattering on observable properties of the atmosphere. We discuss potential observational properties of HD 189733b for the upcoming TESS and CHEOPS missions. We develop a Monte Carlo radiative transfer code and apply it to post-process output of our 3D radiative-hydrodynamic, cloud formation simulation of HD 189733b. We employ three variance reduction techniques; next event estimation, survival biasing and composite emission biasing to improve signal-to-noise of the output.For cloud particle scattering events, a log-normal area distribution is constructed from the 3D cloud formation RHD results and stochastically sampled in order to model the Rayleigh and Mie scattering behaviour of a mixture of grain sizes. Stellar photon packets incident on the eastern dayside hemisphere show predominantly Rayleigh, single-scattering behaviour, while multiple-scattering occurs on the western hemisphere. Combined scattered and thermal emitted light predictions are consistent with published HST and Spitzer secondary transit observations. Our model predictions are also consistent with geometric albedo constraints from optical wavelength ground based polarimetry and HST B-band measurements. We predict an apparent geometric albedo for HD 189733b of 0.205 and 0.229, in the TESS and CHEOPS photometric bands respectively.

Nearly every star known to host planets will become a white dwarf, and nearly 100 planet-hosts are now known to be accompanied by binary stellar companions. Here, we determine how a binary companion triggers instability in otherwise unconditionally stable single-star two-planet systems during the giant branch and white dwarf phases of the planet host. We perform about 700 full-lifetime (14 Gyr) simulations with A0 and F0 primary stars and secondary K2 companions, and identify the critical binary distance within which instability is triggered at any point during stellar evolution. We estimate this distance to be about seven times the outer planet separation, for circular binaries. Our results help characterize the fates of planetary systems, and in particular which ones might yield architectures that are conducive to generating observable heavy metal pollution in white dwarf atmospheres.

Observations of exoplanet atmospheres have revealed the presence of cloud particles in their atmospheres. 3D modelling of cloud formation in atmospheres of extrasolar planets coupled to the atmospheric dynamics has long been a challenge. We investigate the thermo-hydrodynamic properties of cloud formation processes in the atmospheres of hot Jupiter exoplanets. We simulate the dynamic atmosphere of HD 189733b with a 3D model that couples 3D radiative-hydrodynamics with a kinetic, microphysical mineral cloud formation module designed for RHD/GCM exoplanet atmosphere simulations. Our simulation includes the feedback effects of cloud advection and settling, gas phase element advection and depletion/replenishment and the radiative effects of cloud opacity. We model the cloud particles as a mix of mineral materials which change in size and composition as they travel through atmospheric thermo-chemical environments. All local cloud properties such as number density, grain size and material composition are time-dependently calculated. Gas phase element depletion as a result of cloud formation is included in the model. In-situ \textit{effective medium theory} and Mie theory is applied to calculate the wavelength dependent opacity of the cloud component. We present a 3D cloud structure of a chemically complex, gaseous atmosphere of the hot Jupiter HD 189733b. Mean cloud particle sizes are typically sub-micron (0.01-0.5 $\mu$m) at pressures less than 1 bar with hotter equatorial regions containing the smallest grains. Denser cloud structures occur near terminator regions and deeper ($\sim$ 1 bar) atmospheric layers. Silicate materials such as MgSiO$_{3}$[s] are found to be abundant at mid-high latitudes, while TiO$_{2}$[s] and SiO$_{2}$[s] dominate the equatorial regions.

We study the long term orbital evolution of a terrestrial planet under the gravitational perturbations of a giant planet. In particular, we are interested in situations where the two planets are in the same plane and are relatively close. We examine both possible configurations: the giant planet orbit being either outside or inside the orbit of the smaller planet. The perturbing potential is expanded to high orders and an analytical solution of the terrestrial planetary orbit is derived. The analytical estimates are then compared against results from the numerical integration of the full equations of motion and we find that the analytical solution works reasonably well. An interesting finding is that the new analytical estimates improve greatly the predictions for the timescales of the orbital evolution of the terrestrial planet compared to an octupole order expansion. Finally, we briefly discuss possible applications of the analytical estimates in astrophysical problems.

Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-{\alpha} absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption.

3D atmosphere model results are used to comparatively study the kinetic, non-equilibrium cloud formation in the atmospheres of two example planets guided by the giant gas planets HD209458b and HD189733b. Rather independently of hydrodynamic model differences, our cloud modelling suggests that both planets are covered in mineral clouds throughout the entire modelling domain. Both planets harbour chemically complex clouds that are made of mineral particles that have a height-dependent material composition and size. The remaining gas-phase element abundances strongly effects the molecular abundances of the atmosphere in the cloud forming regions. Hydrocarbon and cyanopolyyne molecules can be rather abundant in the inner, dense part of the atmospheres of HD189733b and HD209458b. No one value for metallicity and the C/O ratio can be used to describe an extrasolar planet. Our results concerning the presence and location of water in relation to the clouds explain some of the observed differences between the two planets. In HD189733b, strong water features have been reported while such features appear less strong for HD209458b. By considering the location of the clouds in the two atmospheres, we see that obscuring clouds exist high in the atmosphere of HD209458b, but much deeper in HD189733b. We further conclude that the (self-imposed) degeneracy of cloud parameters in retrieval methods can only be lifted if the cloud formation processes are accurately modelled in contrast to prescribing them by independent parameters.

Detailed observations of the solar system planets reveal a wide variety of local atmospheric conditions. Astronomical observations have revealed a variety of extrasolar planets none of which resembles any of the solar system planets in full. Instead, the most massive amongst the extrasolar planets, the gas giants, appear very similar to the class of (young) Brown Dwarfs which are amongst the oldest objects in the universe. Despite of this diversity, solar system planets, extrasolar planets and Brown Dwarfs have broadly similar global temperatures between 300K and 2500K. In consequence, clouds of different chemical species form in their atmospheres. While the details of these clouds differ, the fundamental physical processes are the same. Further to this, all these objects were observed to produce radio and X-ray emission. While both kinds of radiation are well studied on Earth and to a lesser extent on the solar system planets, the occurrence of emission that potentially originate from accelerated electrons on Brown Dwarfs, extrasolar planets and protoplanetary disks is not well understood yet. This paper offers an interdisciplinary view on electrification processes and their feedback on their hosting environment in meteorology, volcanology, planetology and research on extrasolar planets and planet formation.

We utilize multi-dimensional simulations of varying equatorial jet strength to predict wavelength dependent variations in the eclipse times of gas-giant planets. A displaced hot-spot introduces an asymmetry in the secondary eclipse light curve that manifests itself as a measured offset in the timing of the center of eclipse. A multi-wavelength observation of secondary eclipse, one probing the timing of barycentric eclipse at short wavelengths and another probing at longer wavelengths, will reveal the longitudinal displacement of the hot-spot and break the degeneracy between this effect and that associated with the asymmetry due to an eccentric orbit. The effect of time offsets was first explored in the IRAC wavebands by Williams et. al (2006). Here we improve upon their methodology, extend to a broad ranges of wavelengths, and demonstrate our technique on a series of multi-dimensional radiative-hydrodynamical simulations of HD 209458b with varying equatorial jet strength and hot-spot displacement. Simulations with the largest hot-spot displacement result in timing offsets of up to 100 seconds in the infrared. Though we utilize a particular radiative hydrodynamical model to demonstrate this effect, the technique is model independent. This technique should allow a much larger survey of hot-spot displacements with JWST then currently accessible with time-intensive phase curves, hopefully shedding light on the physical mechanisms associated with thermal energy advection in irradiated gas-giants.

We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind.

Context. Observations suggest that exoplanets such as HD 189733b form clouds in their atmospheres which have a strong feedback onto their thermodynamical and chemical structure, and overall appearance. Aims. Inspired by mineral cloud modelling efforts for Brown Dwarf atmospheres, we present the first spatially varying kinetic cloud model structures for HD 189733b. Methods. We apply a 2-model approach using results from a 3D global radiation-hydrodynamic simulation of the atmosphere as input for a detailed, kinetic cloud formation model. Sampling the 3D global atmosphere structure with 1D trajectories allows us to model the spatially varying cloud structure on HD 189733b. The resulting cloud properties enable the calculation of the scattering and absorption properties of the clouds. Results. We present local and global cloud structure and property maps for HD 189733b. The calculated cloud properties show variations in composition, size and number density of cloud particles which are strongest between the dayside and nightside. Cloud particles are mainly composed of a mix of materials with silicates being the main component. Cloud properties, and hence the local gas composition, change dramatically where temperature inversions occur locally. The cloud opacity is dominated by absorption in the upper atmosphere and scattering at higher pressures in the model. The calculated 8{\mu}m single scattering Albedo of the cloud particles are consistent with Spitzer bright regions. The cloud particles scattering properties suggest that they would sparkle/reflect a midnight blue colour at optical wavelengths.

Three-dimensional equatorial trapped waves excited by stellar isolation and the resulting equatorial superrotating jet in a vertical stratified atmosphere of a tidally-locked hot Jupiter are investigated. Taking the hot Jupiter HD 189733b as a fiducial example, we analytically solve a set of linear equations subject to stationary stellar heating with a uniform zonal-mean flow included. We also extract wave information in the final equilibrium state of the atmosphere from the radiative hydrodynamical simulation for HD 189733b by Dobbs-Dixon & Agol (2013). We find that the analytic wave solutions are able to qualitatively explain the three-dimensional simulation results. Studying the vertical structure of waves allows us to explore new wave features such as the westward tilt of wavefronts related to the Rossby-wave resonance as well as double gyres of dispersive Rossby waves. We also make an attempt to apply our linear wave analysis to explain some numerical features associated with the equatorial jet development seen in the GCM by Showman & Polvani (2011). During the spin-up phase of the equatorial jet, the acceleration of the jet as a result of the divergence of the wave momentum flux can be in principle boosted by the Rossby-wave resonance as the superrotating jet speed is close to the phase speed of the Rossby waves. However, we also find that as the jet speed increases, the Rossby-wave structure shifts eastward, while the Kelvin-wave structure remains approximately stationary, leading to the decline of the acceleration rate. Our analytic model of jet evolution implies that there exists only one stable equilibrium state of the atmosphere, possibly implying that the final state of the atmosphere is independent of initial conditions in the linear regime. Limitations and future improvements are also discussed.

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

The strength and direction of migration of low mass planets depends on the disc's thermodynamics. In discs where the viscous heating is balanced by radiative transport, the migration can be directed outwards, a process which extends the lifetime of growing planetary embryos. We investigate the influence of opacity and stellar irradiation on the disc thermodynamics. Utilizing the resulting disc structure, we determine the regions of outward migration. We perform two-dimensional numerical simulations of equilibrium discs with viscous heating, radiative cooling and stellar irradiation. We use the hydrodynamical code NIRVANA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration is studied by using torque formulae. In the constant opacity case, we reproduce the analytical results of a black-body disc: the stellar irradiation dominates in the outer regions -- leading to flaring -- while the viscous heating dominates close to the star. We find that the inner edge of the disc should not be significantly puffed-up by the stellar irradiation. If the opacity depends on the local density and temperature, the structure of the disc is different, and several bumps in the aspect ratio H/r appear, due to transitions between different opacity regimes. The bumps in the disc can shield the outer disc from stellar irradiation. Stellar irradiation is an important factor for determining the disc structure and has dramatic consequences for the migration of embedded planets. Compared to discs with only viscous heating, a stellar irradiated disc features a much smaller region of outward migration for a range of planetary masses. This suggests that the region where the formation of giant planet cores takes place is smaller, which in turn might lead to a shorter growth phase.

We present a detailed three-dimensional radiative-hydrodynamical simulation of the well known irradiated exoplanet HD189733b. Our model solves the fully compressible Navier-Stokes equations coupled to wavelength-dependent radiative transfer throughout the entire planetary envelope. We provide detailed comparisons between the extensive observations of this system and predictions calculated directly from the numerical models. The atmospheric dynamics is characterized by supersonic winds that fairly efficiently advect energy from the dayside to the nightside. Super-rotating equatorial jets form for a wide range of pressures from 10^-5 to 10 bars while counter rotating jets form at higher latitudes. Calculated transit spectrum agree well with the data from the infrared to the UV including the strong Rayleigh scattering seen at short wavelength, though we slightly under-predict the observations at wavelengths shorter then ~0.6 microns. Our predicted emission spectrum agrees remarkably well at 5.8 and 8 microns, but slightly over-predicts the emission at 3.6 and 4.5 microns when compared to the latest analysis by Knutson et. al (2012). Our simulated IRAC phasecurves agree fairly well with the amplitudes of variations, shape, and phases of minimum and maximum flux. However, we over-predict the peak amplitude at 3.6 and 4.5 microns, and slightly under-predict the location of the phasecurve maximum and minimum. These simulations include, for the first time in a multi-dimensional simulation, a strong Rayleigh scattering component to the absorption opacity, necessary to explain observations in the optical and UV. The agreement between our models and observations suggest that including the effects of condensates in simulations as the dominant form of opacity will be very important in future models.

We present theoretical wavelength-dependent transit light curves for the giant planet HD209458b based on a number of state of the art 3D radiative hydrodynamical models. By varying the kinematic viscosity in the model we calculate observable signatures associated with the emergence of a super-rotating circumplanetary jet that strengthens with decreased viscosity. We obtain excellent agreement between our mid-transit transit spectra and existing data from Hubble and Spitzer, finding the best fit for intermediate values of viscosity. We further exploit dynamically driven differences between eastern and western hemispheres to extract the spectral signal imparted by a circumplanetary jet. We predict that: (i) the transit depth should decrease as the jet becomes stronger; (ii) the measured transit times should show timing offsets of up to 6 seconds at wavelengths with higher opacity, which increases with jet strength; (iii) wavelength-dependent differences between ingress and egress spectra increase with jet strength; (iv) the color-dependent transit shape should vary more strongly with wavelength for stronger jets. These techniques and trends should be valid for other hot Jupiters as well. Observations of transit timing offsets may be accessible with current instrumentation, though the other predictions may require the capabilities JWST and other future missions. Hydrodynamical models utilized solve the 3D Navier-Stokes equations together with decoupled thermal and radiative energy equations and wavelength dependent stellar heating.

We present a new approach for simulating the atmospheric dynamics of the close-in giant planet HD209458b that allows for the decoupling of radiative and thermal energies, direct stellar heating of the interior, and the solution of the full 3D Navier Stokes equations. Simulations reveal two distinct temperature inversions (increasing temperature with decreasing pressure) at the sub-stellar point due to the combined effects of opacity and dynamical flow structure and exhibit instabilities leading to changing velocities and temperatures on the nightside for a range of viscosities. Imposed on the quasi-static background, temperature variations of up to 15% are seen near the terminators and the location of the coldest spot is seen to vary by more than 20 degrees, occasionally appearing west of the anti-solar point. Our new approach introduces four major improvements to our previous methods including simultaneously solving both the thermal energy and radiative equations in both the optical and infrared, incorporating updated opacities, including a more accurate treatment of stellar energy deposition that incorporates the opacity relevant for higher energy stellar photons, and the addition of explicit turbulent viscosity.

Transiting planets provide a unique opportunity to study the atmospheres of extra-solar planets. Radiative hydrodynamical models of the atmosphere provide a crucial link between the physical characteristics of the atmosphere and the observed properties. Here I present results from 3D simulations which solve the full Navier-Stokes equations coupled to a flux-limited diffusion treatment of radiation transfer for planets with 1, 3, and 7 day periods. Variations in opacity amongst models leads to a variation in the temperature differential across the planet, while atmospheric dynamics becomes much more variable at longer orbital periods. I also present 3D radiative simulations illustrating the importance of distinguishing between optical and infrared opacities.

More than two dozen short-period Jupiter-mass gas giant planets have been discovered around nearby solar-type stars in recent years, several of which undergo transits, making them ideal for the detection and characterization of their atmospheres. Here we adopt a three-dimensional radiative hydrodynamical numerical scheme to simulate atmospheric circulation on close-in gas giant planets. In contrast to the conventional GCM and shallow water algorithms, this method does not assume quasi hydrostatic equilibrium and it approximates radiation transfer from optically thin to thick regions with flux-limited diffusion. In the first paper of this series, we consider synchronously-spinning gas giants. We show that a full three-dimensional treatment, coupled with rotationally modified flows and an accurate treatment of radiation, yields a clear temperature transition at the terminator. Based on a series of numerical simulations with varying opacities, we show that the night-side temperature is a strong indicator of the opacity of the planetary atmosphere. Planetary atmospheres that maintain large, interstellar opacities will exhibit large day-night temperature differences, while planets with reduced atmospheric opacities due to extensive grain growth and sedimentation will exhibit much more uniform temperatures throughout their photosphere's. In addition to numerical results, we present a four-zone analytic approximation to explain this dependence.

Extrasolar planets found with radial velocity surveys have masses ranging from several Earth to several Jupiter masses. While mass accretion onto protoplanetary cores in weak-line T-Tauri disks may eventually be quenched by a global depletion of gas, such a mechanism is unlikely to have stalled the growth of some known planetary systems which contain relatively low-mass and close-in planets along with more massive and longer period companions. Here, we suggest a potential solution for this conundrum. In general, supersonic infall of surrounding gas onto a protoplanet is only possible interior to both of its Bondi and Roche radii. At a critical mass, a protoplanet's Bondi and Roche radii are equal to the disk thickness. Above this mass, the protoplanets' tidal perturbation induces the formation of a gap. Although the disk gas may continue to diffuse into the gap, the azimuthal flux across the protoplanets' Roche lobe is quenched. Using two different schemes, we present the results of numerical simulations and analysis to show that the accretion rate increases rapidly with the ratio of the protoplanet's Roche to Bondi radii or equivalently to the disk thickness. In regions with low geometric aspect ratios, gas accretion is quenched with relatively low protoplanetary masses. This effect is important for determining the gas-giant planets' mass function, the distribution of their masses within multiple planet systems around solar type stars, and for suppressing the emergence of gas-giants around low mass stars.

The negligible eccentricity of all extra solar planets with periods less than six days can be accounted for by dissipation of tidal disturbances within their envelopes which are induced by their host stars. In the period range of 7-21 days, planets with circular orbits coexist with planets with eccentric orbits. These will be referred to as the borderline planets. We propose that this discrepancy can be attributed to the variation in spin-down rates of young stars. In particular, prior to spin-down, dissipation of a planet's tidal disturbance within the envelope of a sufficiently rapidly spinning star can excite eccentricity growth, and for a more slowly spinning star, at least reduce the eccentricity damping rate. In contrast, tidal dissipation within the envelope of a slowly spinning low-mass mature star can enhance the eccentricity damping process. Based on these results, we suggest that short-period planets around relatively young stars may have a much larger dispersion in eccentricity than those around mature stars. We also suggest that because the rate of angular momentum loss from G and K dwarfs via stellar winds is much faster than the tidal transfer of angular momentum between themselves and their very-short (3-4 days) period planets, they cannot establish a dynamical configuration in which the stellar and planetary spins are approximately parallel and synchronous with the orbital frequency. In principle, however, such configurations may be established for planets (around G and K dwarfs) with orbital periods of up to several weeks. In contrast to G and K dwarfs, the angular momentum loss due to stellar winds is much weaker in F dwarfs. It is therefore possible for synchronized short-period planets to exist around such stars. The planet around Tau Boo is one such example.