ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

Debris discs around M stars: non-existence versus non-detection

Motivated by the reported dearth of debris discs around M stars, we use survival models to study the occurrence of planetesimal discs around them. These survival models describe a planetesimal disc with a small number of parameters, determine if it may survive a series of dynamical processes and compute the associated infrared excess. For the Wide-field Infrared Survey Explorer (WISE) satellite, we demonstrate that the dearth of debris discs around M stars may be attributed to the small semimajor axes generally probed if either: (1) the dust grains behave like blackbodies emitting at a peak wavelength coincident with the observed one; (2) or the grains are hotter than predicted by their blackbody temperatures and emit at peak wavelengths that are shorter than the observed one. At these small distances from the M star, planetesimals are unlikely to survive or persist for time-scales of 300 Myr or longer if the disc is too massive. Conversely, our survival models allow for the existence of a large population of low-mass debris discs that are too faint to be detected with current instruments. We gain further confidence in our interpretation by demonstrating the ability to compute infrared excesses for Sun-like stars that are broadly consistent with reported values in the literature. However, our interpretation becomes less clear and large infrared excesses are allowed if only one of these scenarios holds: (3) the dust grains are hotter than blackbody and predominantly emit at the observed wavelength; (4) or are blackbody in nature and emit at peak wavelengths longer than the observed one. Both scenarios imply that the parent planetesimals reside at larger distances from the star than inferred if the dust grains behaved like blackbodies. In all scenarios, we show that the infrared excesses detected at 22 μm (via WISE) and 70 μm (via Spitzer) from AU Mic are easily reconciled with its young age (12 Myr). Conversely, the existence of the old debris disc (2–8 Gyr) from GJ 581 is due to the large semimajor axes probed by the Herschel PACS instrument. We elucidate the conditions under which stellar wind drag may be neglected when considering dust populations around M stars. The WISE satellite should be capable of detecting debris discs around young M stars with ages ∼10 Myr.

ON THE EXISTENCE OF SHOCKS IN IRRADIATED EXOPLANETARY ATMOSPHERES

Supersonic flows are expected to exist in the atmospheres of irradiated exoplanets, but the question of whether shocks develop lingers. Specifically, it reduces to whether continuous flow in a closed loop may become supersonic and if some portions of the supersonic flow steepen into shocks. We first demonstrate that continuous, supersonic flow may exist in two flavors: isentropic and non-isentropic, with shocks being included in the latter class of solutions. Supersonic flow is a necessary but insufficient condition for shocks to develop. The development of a shock requires the characteristics of neighboring points in a flow to intersect. We demonstrate that the intersection of characteristics may be quantified via the knowledge of the Mach number. Finally, we examine three-dimensional simulations of hot Jovian atmospheres and demonstrate that shock formation is expected to occur mostly on the dayside hemisphere, upstream of the substellar point, because the enhanced temperatures near the substellar point provide a natural pressure barrier for the returning flow. Understanding the role of shocks in irradiated exoplanetary atmospheres is relevant to correctly modeling observables such as the peak offsets of infrared phase curves.