Source author record

Kristen Menou

Kristen Menou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

ExoPlaSim: Extending the Planet Simulator for Exoplanets

The discovery of a large number of terrestrial exoplanets in the habitable zones of their stars, many of which are qualitatively different from Earth, has led to a growing need for fast and flexible 3D climate models, which could model such planets and explore multiple possible climate states and surface conditions. We respond to that need by creating ExoPlaSim, a modified version of the Planet Simulator (PlaSim) that is designed to be applicable to synchronously rotating terrestrial planets, planets orbiting stars with non-solar spectra, and planets with non-Earth-like surface pressures. In this paper we describe our modifications, present validation tests of ExoPlaSim's performance against other GCMs, and demonstrate its utility by performing two simple experiments involving hundreds of models. We find that ExoPlaSim agrees qualitatively with more-sophisticated GCMs such as ExoCAM, LMDG, and ROCKE-3D, falling within the ensemble distribution on multiple measures. The model is fast enough that it enables large parameter surveys with hundreds to thousands of models, potentially enabling the efficient use of a 3D climate model in retrievals of future exoplanet observations. We describe our efforts to make ExoPlaSim accessible to non-modellers, including observers, non-computational theorists, students, and educators through a new Python API and streamlined installation through pip, along with online documentation.

preprint2022arXiv

Fundamental Challenges to Remote Sensing of Exo-Earths

Inferring the climate and surface conditions of terrestrial exoplanets in the habitable zone is a major goal for the field of exoplanet science. This pursuit will require both statistical analyses of the population of habitable planets as well as in-depth analyses of the climates of individual planets. Given the close relationship between habitability and surface liquid water, it is important to ask whether the fraction of a planet's surface where water can be a liquid, $χ_\text{hab}$, can be inferred from observations. We have produced a diverse bank of 1,874 3D climate models and computed the full-phase reflectance and emission spectrum for each model to investigate whether surface climate inference is feasible with high-quality direct imaging or secondary eclipse spectroscopy. These models represent the outcome of approximately 200,000 total simulated years of climate and over 50,000 CPU-hours, and the roughly-100 GB model bank and its associated spectra are being made publicly-available for community use. We find that there are correlations between spectra and $χ_\text{hab}$ that will permit statistical approaches. However, spectral degeneracies in the climate observables produced by our model bank indicate that inference of individual climates is likely to be model-dependent, and inference will likely be impossible without exhaustive explorations of the climate parameter space. The diversity of potential climates on habitable planets therefore poses fundamental challenges to remote sensing efforts targeting exo-Earths.

preprint2022arXiv

TOI-1452 b: SPIRou and TESS reveal a super-Earth in a temperate orbit transiting an M4 dwarf

Exploring the properties of exoplanets near or inside the radius valley provides insights on the transition from the rocky super-Earths to the larger, hydrogen-rich atmosphere mini-Neptunes. Here, we report the discovery of TOI-1452 b, a transiting super-Earth ($R_{\rm p} = 1.67 \pm 0.07$ R$_{\oplus}$) in an 11.1--day temperate orbit ($T_{\rm eq} = 326 \pm 7$ K) around the primary member ($H = 10.0$, $T_{\rm eff} = 3185 \pm 50$ K) of a nearby visual binary M dwarf. The transits were first detected by TESS, then successfully isolated between the two $3.2^{\prime\prime}$ companions with ground-based photometry from OMM and MuSCAT3. The planetary nature of TOI-1452 b was established through high-precision velocimetry with the near-infrared SPIRou spectropolarimeter as part of the ongoing SPIRou Legacy Survey. The measured planetary mass ($4.8 \pm 1.3$ M$_{\oplus}$) and inferred bulk density ($5.6^{+1.8}_{-1.6}$ g/cm$^3$) is suggestive of a rocky core surrounded by a volatile-rich envelope. More quantitatively, the mass and radius of TOI-1452 b, combined with the stellar abundance of refractory elements (Fe, Mg and Si) measured by SPIRou, is consistent with a core mass fraction of $18\pm6$ % and a water mass fraction of $22^{+21}_{-13}$%. The water world candidate TOI-1452 b is a prime target for future atmospheric characterization with JWST, featuring a Transmission Spectroscopy Metric similar to other well-known temperate small planets such as LHS 1140 b and K2-18 b. The system is located near Webb's northern Continuous Viewing Zone, implying that is can be followed at almost any moment of the year.

preprint2021arXiv

Climate uncertainties caused by unknown land distribution on habitable M-Earths

A planet's surface conditions can significantly impact its climate and habitability. In this study, we use the 3D general circulation model ExoPlaSim to systematically vary dayside land cover on a synchronously rotating, temperate rocky planet under two extreme and opposite continent configurations, in which either all of the land or all of the ocean is centred at the substellar point. We identify water vapour and sea ice as competing drivers of climate, and we isolate land-dependent regimes under which one or the other dominates. We find that the amount and configuration of land can change the planet's globally averaged surface temperature by up to 20K, and its atmospheric water vapour content by several orders of magnitude. The most discrepant models have partial dayside land cover with opposite continent configurations. Since transit spectroscopy may permit observations of M-dwarf planets' atmospheres, but not their surfaces, these land-related climate differences likely represent a limiting uncertainty in a given planet's climate, even if its atmospheric composition is known. Our results are robust to variations in atmospheric CO2 concentration, stellar temperature, and instellation.

preprint2020arXiv

Automated crater shape retrieval using weakly-supervised deep learning

Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while simultaneously producing a mask for each crater that traces its outer rim. Our post-processing pipeline then finds the closest fitting ellipse to these masks, allowing us to retrieve the crater ellipticities. Our model is able to correctly identify 87\% of known craters in the longitude range we hid from the network during training and validation (test set), while predicting thousands of additional craters not present in our training data. Manual validation of a subset of these "new" craters indicates that a majority of them are real, which we take as an indicator of the strength of our model in learning to identify craters, despite incomplete training data. The crater size, ellipticity, and depth distributions predicted by our model are consistent with human-generated results. The model allows us to perform a large scale search for differences in crater diameter and shape distributions between the lunar highlands and maria, and we exclude any such differences with a high statistical significance. The predicted test set catalogue and trained model are available here: https://github.com/malidib/Craters_MaskRCNN/.

preprint2020arXiv

Climate Diversity in the Solar-Like Habitable Zone due to Varying Background Gas Pressure

A large number of studies have responded to the growing body of confirmed terrestrial habitable zone exoplanets by presenting models of various possible climates. However, the impact of the partial pressure of background gases such as N$_2$ has not yet been well-explored, despite the abundance of N$_2$ in Earth's atmosphere and the lack of constraints on its typical abundance in terrestrial planet atmospheres. We use PlaSim, a fast 3D climate model, to simulate many hundreds of climates around Sun-like stars with varying N$_2$ partial pressures, instellations, and surface characteristics to identify the impact of the background gas partial pressure on the climate. We find that the climate's response is nonlinear and highly sensitive to the background gas partial pressure. We identify pressure broadening of greenhouse gas (such as CO$_2$ and H$_2$O) absorption lines, amplification of warming or cooling by the water vapor greenhouse positive feedback, heat transport efficiency, and cooling through Rayleigh scattering as the dominant competing mechanisms that determine the equilibrium climate for a given N$_2$ partial pressure. Finally, we show that different amounts of N$_2$ should have a significant effect on broadband reflected light observations of terrestrial exoplanets around Sun-like stars.

preprint2020arXiv

Predicting the long-term stability of compact multiplanet systems

We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.

preprint2019arXiv

Evolution of the radius valley around low mass stars from $Kepler$ and $K2$

We present calculations of the occurrence rate of small close-in planets around low mass dwarf stars using the known planet populations from the $Kepler$ and $K2$ missions. Applying completeness corrections clearly reveals the radius valley in the maximum a-posteriori occurrence rates as a function of orbital separation and planet radius. We measure the slope of the valley to be $r_{p,\text{valley}} \propto F^{-0.060\pm 0.025}$ which bears the opposite sign from that measured around Sun-like stars thus suggesting that thermally driven atmospheric mass loss may not dominate the evolution of planets in the low stellar mass regime or that we are witnessing the emergence of a separate channel of planet formation. The latter notion is supported by the relative occurrence of rocky to non-rocky planets increasing from $0.5\pm 0.1$ around mid-K dwarfs to $8.5\pm 4.6$ around mid-M dwarfs. Furthermore, the center of the radius valley at $1.54\pm 0.16$ R$_{\oplus}$ is shown to shift to smaller sizes with decreasing stellar mass in agreement with physical models of photoevaporation, core-powered mass loss, and gas-poor formation. Although current measurements are insufficient to robustly identify the dominant formation pathway of the radius valley, such inferences may be obtained by $TESS$ with $\mathcal{O}(85,000)$ mid-to-late M dwarfs observed with 2-minute cadence. The measurements presented herein also precisely designate the subset of planetary orbital periods and radii that should be targeted in radial velocity surveys to resolve the rocky to non-rocky transition around low mass stars.

preprint2019arXiv

Hot Jupiter Atmospheric Flows at High Resolution

Global Circulation Models (GCMs) of atmospheric flows are now routinely used to interpret observational data on Hot Jupiters. Localized "equatorial $β$-plane" simulations by Fromang et al. (2016) have revealed that a barotropic (horizontal shear) instability of the equatorial jet appears at horizontal resolutions beyond those typically achieved in global models; this instability could limit wind speeds and lead to increased atmospheric variability. To address this possibility, we adapt the computationally efficient, pseudo-spectral PlaSim GCM, originally designed for Earth studies, to model Hot Jupiter atmospheric flows and validate it on the Heng et al. (2011) reference benchmark. We then present high resolution global models of HD209458b, with horizontal resolutions of T85 (128x256) and T127 (192x384). The barotropic instability phenomenology found in $β$-plane simulations is not reproduced in these global models, despite comparably high resolutions. Nevertheless, high resolution models do exhibit additional flow variability on long timescales (of order 100 planet days or more), which is absent from the lower resolution models. It manifests as a breakdown of north-south symmetry of the equatorial wind. From post-processing the atmospheric flows at various resolutions (assuming a cloud-free situation), we show that the stronger flow variability achieved at high resolution does not translate into noticeably stronger dayside infrared flux variability. More generally, our results suggest that high horizontal resolutions are not required to capture the key features of hot Jupiter atmospheric flows.

preprint2016arXiv

Disk-fed giant planet formation

Massive giant planets, such as the ones being discovered by direct imaging surveys, likely experience the majority of their growth through a circumplanetary disc. We argue that the entropy of accreted material is determined by boundary layer processes, unlike the "cold-" or "hot-start" hypotheses usually invoked in the core accretion and direct collapse scenarios. A simple planetary evolution model illustrates how a wide range of radius and luminosity tracks become possible, depending on details of the accretion process. Specifically, the proto-planet evolves towards "hot-start" tracks if the scale-height of the boundary layer is $\gtrsim0.24$, a value not much larger than the scale-height of the circumplanetary disc. Understanding the luminosity and radii of young giant planets will thus require detailed models of circumplanetary accretion.

preprint2016arXiv

On the radial velocity detection of additional planets in transiting, slowly rotating M-dwarf systems: the case of GJ 1132

M-dwarfs are known to commonly host high-multiplicity planetary systems. Therefore M-dwarf planetary systems with a known transiting planet are expected to contain additional small planets ($r_p \le 4$ R$_{\oplus}$, $m_p \lesssim 20$ M$_{\oplus}$) that are not seen in transit. In this study we investigate the effort required to detect such planets using precision velocimetry around the sizable subset of M-dwarfs which are slowly rotating ($P_{\mathrm{rot}} \gtrsim 40$ days) and hence more likely to be inactive. We focus on the test case of GJ 1132. Specifically, we perform a suite of Monte-Carlo simulations of the star's radial velocity signal featuring astrophysical contributions from stellar jitter due to rotationally modulated active regions and keplarian signals from the known transiting planet and hypothetical additional planets not seen in transit. We then compute the detection completeness of non-transiting planets around GJ 1132 and consequently estimate the number of RV measurements required to detect those planets. We show that with 1 m s$^{-1}$ RV precision per measurement, only $\sim 50$ measurements are required to achieve a 50\% detection completeness to all non-transiting planets in the system and to planets which are potentially habitable. Throughout we advocate the use of Gaussian process regression as an effective tool for mitigating the effects of stellar jitter including stars with high activity. Given that GJ 1132 is representative of a large population of slowly rotating M-dwarfs, we conclude with a discussion of how our results may be extended to other systems with known transiting planets such as those which will be discovered with TESS.

preprint2015arXiv

Asynchronous rotation of Earth-mass planets in the habitable zone of lower-mass stars

Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than 0.5-0.7Msun. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.

preprint2015arXiv

Dynamical Stability of Imaged Planetary Systems in Formation: Application to HL Tau

A recent ALMA image revealed several concentric gaps in the protoplanetary disk surrounding the young star HL Tau. We consider the hypothesis that these gaps are carved by planets, and present a general framework for understanding the dynamical stability of such systems over typical disk lifetimes, providing estimates for the maximum planetary masses. We collect these easily evaluated constraints into a workflow that can help guide the design and interpretation of new observational campaigns and numerical simulations of gap opening in such systems. We argue that the locations of resonances should be significantly shifted in massive disks like HL Tau, and that theoretical uncertainties in the exact offset, together with observational errors, imply a large uncertainty in the dynamical state and stability in such disks. This presents an important barrier to using systems like HL Tau as a proxy for the initial conditions following planet formation. An important observational avenue to breaking this degeneracy is to search for eccentric gaps, which could implicate resonantly interacting planets. Unfortunately, massive disks like HL Tau should induce swift pericenter precession that would smear out any such eccentric features of planetary origin. This motivates pushing toward more typical, less massive disks. For a nominal non-resonant model of the HL Tau system with five planets, we find a maximum mass for the outer three bodies of approximately 2 Neptune masses. In a resonant configuration, these planets can reach at least the mass of Saturn. The inner two planets' masses are unconstrained by dynamical stability arguments.

preprint2014arXiv

Climate Stability of Habitable Earth-like Planets

The carbon-silicate cycle regulates the atmospheric $CO_2$ content of terrestrial planets on geological timescales through a balance between the rates of $CO_2$ volcanic outgassing and planetary intake from rock weathering. It is thought to act as an efficient climatic thermostat on Earth and, by extension, on other habitable planets. If, however, the weathering rate increases with the atmospheric $CO_2$ content, as expected on planets lacking land vascular plants, the carbon-silicate cycle feedback can become severely limited. Here we show that Earth-like planets receiving less sunlight than current Earth may no longer possess a stable warm climate but instead repeatedly cycle between unstable glaciated and deglaciated climatic states. This has implications for the search for life on exoplanets in the habitable zone of nearby stars.

preprint2013arXiv

Water-Trapped Worlds

Although tidally-locked habitable planets orbiting nearby M-dwarf stars are among the best astronomical targets to search for extrasolar life, they may also be deficient in volatiles and water. Climate models for this class of planets show atmospheric transport of water from the dayside to the nightside, where it is precipitated as snow and trapped as ice. Since ice only slowly flows back to the dayside upon accumulation, the resulting hydrological cycle can trap a large amount of water in the form of nightside ice. Using ice sheet dynamical and thermodynamical constraints, I illustrate how planets with less than about a quarter the Earth's oceans could trap most of their surface water on the nightside. This would leave their dayside, where habitable conditions are met, potentially dry. The amount and distribution of residual liquid water on the dayside depend on a variety of geophysical factors, including the efficiency of rock weathering at regulating atmospheric CO2 as dayside ocean basins dry-up. Water-trapped worlds with dry daysides may offer similar advantages as land planets for habitability, by contrast with worlds where more abundant water freely flows around the globe.

preprint2012arXiv

Planetesimal and Protoplanet Dynamics in a Turbulent Protoplanetary Disk: Ideal Stratified Disks

Due to the gravitational influence of density fluctuations driven by magneto-rotational instability in the gas disk, planetesimals and protoplanets undergo diffusive radial migration as well as changes in other orbital properties. The magnitude of the effect on particle orbits can have important consequences for planet formation scenarios. We use the local-shearing-box approximation to simulate an ideal, isothermal, magnetized gas disk with vertical density stratification and simultaneously evolve numerous massless particles moving under the gravitational field of the gas and the host star. We measure the evolution of the particle orbital properties, including mean radius, eccentricity, inclination, and velocity dispersion, and its dependence on the disk properties and the particle initial conditions. Although the results converge with resolution for fixed box dimensions, we find the response of the particles to the gravity of the turbulent gas correlates with the horizontal box size, up to 16 disk scale heights. This correlation indicates that caution should be exercised when interpreting local-shearing-box models involving gravitational physics of magneto-rotational turbulence. Based on heuristic arguments, nevertheless, the criterion L_h / R ~ O(1), where L_h is the horizontal box size and R is the distance to the host star, is proposed to possibly circumvent this conundrum. If this criterion holds, we can still conclude that magneto-rotational turbulence seems likely to be ineffective at driving either diffusive migration or collisional erosion under most circumstances.

preprint2012arXiv

Thermo-Resistive Instability of Hot Planetary Atmospheres

The atmospheres of hot Jupiters and other strongly-forced exoplanets are susceptible to a thermal instability in the presence of ohmic dissipation, weak magnetic drag and strong winds. The instability occurs in radiatively-dominated atmospheric regions when the ohmic dissipation rate increases with temperature faster than the radiative (cooling) rate. The instability domain covers a specific range of atmospheric pressures and temperatures, typically P ~ 3-300 mbar and T ~ 1500-2500K for hot Jupiters, which makes it a candidate mechanism to explain the dayside thermal "inversions" inferred for a number of such exoplanets. The instability is suppressed by high levels of non-thermal photoionization, in possible agreement with a recently established observational trend. We highlight several shortcomings of the instability treatment presented here. Understanding the emergence and outcome of the instability, which should result in locally hotter atmospheres with stronger levels of drag, will require global non-linear atmospheric models with adequate MHD prescriptions.

preprint2011arXiv

A General Circulation Model for Gaseous Exoplanets with Double-Gray Radiative Transfer

We present a new version of our code for modeling the atmospheric circulation on gaseous exoplanets, now employing a "double-gray" radiative transfer scheme, which self-consistently solves for fluxes and heating throughout the atmosphere, including the emerging (observable) infrared flux. We separate the radiation into infrared and optical components, each with its own absorption coefficient, and solve standard two-stream radiative transfer equations. We use a constant optical absorption coefficient, while the infrared coefficient can scale as a powerlaw with pressure. Here we describe our new code in detail and demonstrate its utility by presenting a generic hot Jupiter model. We discuss issues related to modeling the deepest pressures of the atmosphere and describe our use of the diffusion approximation for radiative fluxes at high optical depths. In addition, we present new models using a simple form for magnetic drag on the atmosphere. We calculate emitted thermal phase curves and find that our drag-free model has the brightest region of the atmosphere offset by ~12 degrees from the substellar point and a minimum flux that is 17% of the maximum, while the model with the strongest magnetic drag has an offset of only ~2 degrees and a ratio of 13%. Finally, we calculate rates of numerical loss of kinetic energy at ~15% for every model except for our strong-drag model, where there is no measurable loss; we speculate that this is due to the much decreased wind speeds in that model.

preprint2011arXiv

Atmospheres of Hot Super-Earths

Hot super-Earths likely possess minimal atmospheres established through vapor saturation equilibrium with the ground. We solve the hydrodynamics of these tenuous atmospheres at the surface of Corot-7b, Kepler 10b and 55 Cnc-e, including idealized treatments of magnetic drag and ohmic dissipation. We find that atmospheric pressures remain close to their local saturation values in all cases. Despite the emergence of strongly supersonic winds which carry sublimating mass away from the substellar point, the atmospheres do not extend much beyond the day-night terminators. Ground temperatures, which determine the planetary thermal (infrared) signature, are largely unaffected by exchanges with the atmosphere and thus follow the effective irradiation pattern. Atmospheric temperatures, however, which control cloud condensation and thus albedo properties, can deviate substantially from the irradiation pattern. Magnetic drag and ohmic dissipation can also strongly impact the atmospheric behavior, depending on atmospheric composition and the planetary magnetic field strength. We conclude that hot super-Earths could exhibit interesting signatures in reflection (and possibly in emission) which would trace a combination of their ground, atmospheric and magnetic properties.

preprint2011arXiv

Atmospheric Circulation and Composition of GJ1214b

The exoplanet GJ1214b presents an interesting example of compositional degeneracy for low-mass planets. Its atmosphere may be composed of water, super-solar or solar metallicity material. We present atmospheric circulation models of GJ1214b for these three compositions, with explicit grey radiative transfer and an optional treatment of MHD bottom drag. All models develop strong, superrotating zonal winds (~ 1-2 km/s). The degree of eastward heat advection, which can be inferred from secondary eclipse and thermal phase curve measurements, varies greatly between the models. These differences are understood as resulting from variations in the radiative times at the thermal photosphere, caused by separate molecular weight and opacity effects. Our GJ1214b models illustrate how atmospheric circulation can be used as a probe of composition for similar tidally-locked exoplanets in the mini-Neptune/waterworld class.

preprint2011arXiv

Atmospheric circulation of tidally locked exoplanets: a suite of benchmark tests for dynamical solvers

The complexity of atmospheric modelling and its inherent non-linearity, together with the limited amount of data of exoplanets available, motivate model intercomparisons and benchmark tests. In the geophysical community, the Held-Suarez test is a standard benchmark for comparing dynamical core simulations of the Earth's atmosphere with different solvers, based on statistically-averaged flow quantities. In the present study, we perform analogues of the Held-Suarez test for tidally-locked exoplanets with the GFDL-Princeton Flexible Modeling System (FMS) by subjecting both the spectral and finite difference dynamical cores to a suite of tests, including the standard benchmark for Earth, a hypothetical tidally-locked Earth, a "shallow" hot Jupiter model and a "deep" model of HD 209458b. We find qualitative and quantitative agreement between the solvers for the Earth, tidally-locked Earth and shallow hot Jupiter benchmarks, but the agreement is less than satisfactory for the deep model of HD 209458b. Further investigation reveals that closer agreement may be attained by arbitrarily adjusting the values of the horizontal dissipation parameters in the two solvers, but it remains the case that the magnitude of the horizontal dissipation is not easily specified from first principles. Irrespective of radiative transfer or chemical composition considerations, our study points to limitations in our ability to accurately model hot Jupiter atmospheres with meteorological solvers at the level of ten percent for the temperature field and several tens of percent for the velocity field. Direct wind measurements should thus be particularly constraining for the models. Our suite of benchmark tests also provides a reference point for researchers wishing to adapt their codes to study the atmospheric circulation regimes of tidally-locked Earths/Neptunes/Jupiters.

preprint2011arXiv

Electromagnetic counterparts of supermassive black hole binaries resolved by pulsar timing arrays

Pulsar timing arrays (PTAs) are expected to detect gravitational waves (GWs) from individual low-redshift (z<1.5) compact supermassive (M>10^9 Msun) black hole (SMBH) binaries with orbital periods of approx. 0.1 - 10 yrs. Identifying the electromagnetic (EM) counterparts of these sources would provide confirmation of putative direct detections of GWs, present a rare opportunity to study the environments of compact SMBH binaries, and could enable the use of these sources as standard sirens for cosmology. Here we consider the feasibility of such an EM identification. We show that because the host galaxies of resolved PTA sources are expected to be exceptionally massive and rare, it should be possible to find unique hosts of resolved sources out to redshift z=0.2. At higher redshifts, the PTA error boxes are larger, and may contain as many as 100 massive-galaxy interlopers. The number of candidates, however, remains tractable for follow-up searches in upcoming wide-field EM surveys. We develop a toy model to characterize the dynamics and the thermal emission from a geometrically thin, gaseous disc accreting onto a PTA-source SMBH binary. Our model predicts that at optical and infrared frequencies, the source should appear similar to a typical luminous active galactic nucleus (AGN). However, owing to the evacuation of the accretion flow by the binary's tidal torques, the source might have an unusually low soft X-ray luminosity and weak UV and broad optical emission lines, as compared to an AGN powered by a single SMBH with the same total mass. For sources near z=1, the decrement in the rest-frame UV should be observable as an extremely red optical color. These properties would make the PTA sources stand out among optically luminous AGN, and could allow their unique identification.

preprint2011arXiv

Magnetic Scaling Laws for the Atmospheres of Hot Giant Exoplanets

We present scaling laws for advection, radiation, magnetic drag and ohmic dissipation in the atmospheres of hot giant exoplanets. In the limit of weak thermal ionization, ohmic dissipation increases with the planetary equilibrium temperature (T_eq >~ 1000 K) faster than the insolation power does, eventually reaching values >~ 1% of the insolation power, which may be sufficient to inflate the radii of hot Jupiters. At higher T_eq values still, magnetic drag rapidly brakes the atmospheric winds, which reduces the associated ohmic dissipation power. For example, for a planetary field strength B=10G, the fiducial scaling laws indicate that ohmic dissipation exceeds 1% of the insolation power over the equilibrium temperature range T_eq ~ 1300-2000 K, with a peak contribution at T_eq ~ 1600 K. Evidence for magnetically dragged winds at the planetary thermal photosphere could emerge in the form of reduced longitudinal offsets for the dayside infrared hotspot. This suggests the possibility of an anticorrelation between the amount of hotspot offset and the degree of radius inflation, linking the atmospheric and interior properties of hot giant exoplanets in an observationally testable way. While providing a useful framework to explore the magnetic scenario, the scaling laws also reveal strong parameter dependencies, in particular with respect to the unknown planetary magnetic field strength.

preprint2011arXiv

The Role of Drag in the Energetics of Strongly Forced Exoplanet Atmospheres

In contrast to the Earth, where frictional heating is typically negligible, we show that drag mechanisms could act as an important heat source in the strongly-forced atmospheres of some exoplanets, with the potential to alter the circulation. We modify the standard formalism of the atmospheric energy cycle to explicitly track the loss of kinetic energy and the associated frictional (re)heating, for application to exoplanets such as the asymmetrically heated "hot Jupiters" and gas giants on highly eccentric orbits. We establish that an understanding of the dominant drag mechanisms and their dependence on local atmospheric conditions is critical for accurate modeling, not just in their ability to limit wind speeds, but also because they could possibly change the energetics of the circulation enough to alter the nature of the flow. We discuss possible sources of drag and estimate the strength necessary to significantly influence the atmospheric energetics. As we show, the frictional heating depends on the magnitude of kinetic energy dissipation as well as its spatial variation, so that the more localized a drag mechanism is, the weaker it can be and still affect the circulation. We also use the derived formalism to estimate the rate of numerical loss of kinetic energy in a few previously published hot Jupiter models with and without magnetic drag and find it to be surprisingly large, at 5-10% of the incident stellar irradiation.

preprint2010arXiv

A New 24 micron Phase Curve for upsilon Andromedae b

We report the detection of 24 micron variations from the planet-hosting upsilon Andromedae system consistent with the orbital periodicity of the system's innermost planet, upsilon And b. We find a peak-to-valley phase curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of constant surface brightness and assuming a planetary radius of 1.3 Jupiter radii gives a planetary temperature contrast of >900 K and an orbital inclination of >28 degrees. We further report the largest phase offset yet observed for an extrasolar planet: the flux maximum occurs ~80 degrees before phase 0.5. Such a large phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on earlier observations of this system in several important ways: (1) observations of a flux calibrator star demonstrate the MIPS detector is stable to 10^-4 on long timescales, (2) we note that the background light varies systematically due to spacecraft operations, precluding use of this background as a flux calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier observations of this system is consistent with our new result.

preprint2010arXiv

Habitable Climates: The Influence of Eccentricity

In the outer regions of the habitable zone, the risk of transitioning into a globally frozen "snowball" state poses a threat to the habitability of planets with the capacity to host water-based life. We use a one-dimensional energy balance climate model (EBM) to examine how obliquity, spin rate, orbital eccentricity, and ocean coverage might influence the onset of such a snowball state. For an exoplanet, these parameters may be strikingly different from the values observed for Earth. Since, for constant semimajor axis, the annual mean stellar irradiation scales with (1-e^2)^(-1/2), one might expect the greatest habitable semimajor axis (for fixed atmospheric composition) to scale as (1-e^2)^(-1/4). We find that this standard ansatz provides a reasonable lower bound on the outer boundary of the habitable zone, but the influence of obliquity and ocean fraction can be profound in the context of planets on eccentric orbits. For planets with eccentricity 0.5, our EBM suggests that the greatest habitable semimajor axis can vary by more than 0.8 AU (78%!) depending on obliquity, with higher obliquity worlds generally more stable against snowball transitions. One might also expect that the long winter at an eccentric planet's apoastron would render it more susceptible to global freezing. Our models suggest that this is not a significant risk for Earth-like planets around Sun-like stars since such planets are buffered by the thermal inertia provided by oceans covering at least 10% of their surface. Since planets on eccentric orbits spend much of their year particularly far from the star, such worlds might turn out to be especially good targets for direct observations with missions such as TPF-Darwin. Nevertheless, the extreme temperature variations achieved on highly eccentric exo-Earths raise questions about the adaptability of life to marginally or transiently habitable conditions.

preprint2010arXiv

Magnetic Drag on Hot Jupiter Atmospheric Winds

Hot Jupiters, with atmospheric temperatures T ~ 1000 K, have residual thermal ionization levels sufficient for the interaction of the ions with the planetary magnetic field to result in a sizable magnetic drag on the (neutral) atmospheric winds. We evaluate the magnitude of magnetic drag in a representative three-dimensional atmospheric model of the hot Jupiter HD 209458b and find that it is a plausible mechanism to limit wind speeds in this class of atmospheres. Magnetic drag has a strong geometrical dependence, both meridionally and from the day to the night side (in the upper atmosphere), which could have interesting consequences for the atmospheric flow pattern. By extension, close-in eccentric planets with transiently heated atmospheres will experience time-variable levels of magnetic drag. A robust treatment of magnetic drag in circulation models for hot atmospheres may require iterated solutions to the magnetic induction and Saha equations as the hydrodynamical flow is evolved.

preprint2010arXiv

Ohmic Dissipation in the Atmospheres of Hot Jupiters

Hot Jupiter atmospheres exhibit fast, weakly-ionized winds. The interaction of these winds with the planetary magnetic field generates drag on the winds and leads to ohmic dissipation of the induced electric currents. We study the magnitude of ohmic dissipation in representative, three-dimensional atmospheric circulation models of the hot Jupiter HD 209458b. We find that ohmic dissipation can reach or exceed 1% of the stellar insolation power in the deepest atmospheric layers, in models with and without dragged winds. Such power, dissipated in the deep atmosphere, appears sufficient to slow down planetary contraction and explain the typically inflated radii of hot Jupiters. This atmospheric scenario does not require a top insulating layer or radial currents that penetrate deep in the planetary interior. Circulation in the deepest atmospheric layers may actually be driven by spatially non-uniform ohmic dissipation. A consistent treatment of magnetic drag and ohmic dissipation is required to further elucidate the consequences of magnetic effects for the atmospheres and the contracting interiors of hot Jupiters.

preprint2010arXiv

Photometric and Spectral Signatures of 3D Models of Transiting Giant Exoplanets

Using a 3D GCM, we create dynamical model atmospheres of a representative transiting giant exoplanet, HD 209458b. We post-process these atmospheres with an opacity code to obtain transit radius spectra during the primary transit. Using a spectral atmosphere code, we integrate over the face of the planet seen by an observer at various orbital phases and calculate light curves as a function of wavelength and for different photometric bands. The products of this study are generic predictions for the phase variations of a zero-eccentricity giant planet's transit spectrum and of its light curves. We find that for these models the temporal variations in all quantities and the ingress/egress contrasts in the transit radii are small ($< 1.0$\%). Moreover, we determine that the day/night contrasts and phase shifts of the brightness peaks relative to the ephemeris are functions of photometric band. The $J$, $H$, and $K$ bands are shifted most, while the IRAC bands are shifted least. Therefore, we verify that the magnitude of the downwind shift in the planetary ``hot spot" due to equatorial winds is strongly wavelength-dependent. The phase and wavelength dependence of light curves, and the associated day/night contrasts, can be used to constrain the circulation regime of irradiated giant planets and to probe different pressure levels of a hot Jupiter atmosphere. We posit that though our calculations focus on models of HD 209458b similar calculations for other transiting hot Jupiters in low-eccentricity orbits should yield transit spectra and light curves of a similar character.

preprint2010arXiv

Radiation-Hydrodynamics of Hot Jupiter Atmospheres

Radiative transfer in planetary atmospheres is usually treated in the static limit, i.e., neglecting atmospheric motions. We argue that hot Jupiter atmospheres, with possibly fast (sonic) wind speeds, may require a more strongly coupled treatment, formally in the regime of radiation-hydrodynamics. To lowest order in v/c, relativistic Doppler shifts distort line profiles along optical paths with finite wind velocity gradients. This leads to flow-dependent deviations in the effective emission and absorption properties of the atmospheric medium. Evaluating the overall impact of these distortions on the radiative structure of a dynamic atmosphere is non-trivial. We present transmissivity and systematic equivalent width excess calculations which suggest possibly important consequences for radiation transport in hot Jupiter atmospheres. If winds are fast and bulk Doppler shifts are indeed important for the global radiative balance, accurate modeling and reliable data interpretation for hot Jupiter atmospheres may prove challenging: it would involve anisotropic and dynamic radiative transfer in a coupled radiation-hydrodynamical flow. On the bright side, it would also imply that the emergent properties of hot Jupiter atmospheres are more direct tracers of their atmospheric flows than is the case for Solar System planets. Radiation-hydrodynamics may also influence radiative transfer in other classes of hot exoplanetary atmospheres with fast winds.

preprint2010arXiv

Three Dimensional Modeling of Hot Jupiter Atmospheric Flows

We present a three dimensional hot Jupiter model, extending from 200 bar to 1 mbar, using the Intermediate General Circulation Model from the University of Reading. Our horizontal spectral resolution is T31 (equivalent to a grid of 48x96), with 33 logarithmically spaced vertical levels. A simplified (Newtonian) scheme is employed for the radiative forcing. We adopt a physical set up nearly identical to the model of HD 209458b by Cooper & Showman (2005,2006) to facilitate a direct model inter-comparison. Our results are broadly consistent with theirs but significant differences also emerge. The atmospheric flow is characterized by a super-rotating equatorial jet, transonic wind speeds, and eastward advection of heat away from the dayside. We identify a dynamically-induced temperature inversion ("stratosphere") on the planetary dayside and find that temperatures at the planetary limb differ systematically from local radiative equilibrium values, a potential source of bias for transit spectroscopic interpretations. While our model atmosphere is quasi-identical to that of Cooper & Showman (2005,2006) and we solve the same meteorological equations, we use different algorithmic methods, spectral-implicit vs. grid-explicit, which are known to yield fully consistent results in the Earth modeling context. The model discrepancies identified here indicate that one or both numerical methods do not faithfully capture all of the atmospheric dynamics at work in the hot Jupiter context. We highlight the emergence of a shock-like feature in our model, much like that reported recently by Showman et al. (2009), and suggest that improved representations of energy conservation may be needed in hot Jupiter atmospheric models, as emphasized by Goodman (2009).

preprint2010arXiv

Time-Dependent Models for the Afterglows of Massive Black Hole Mergers

The Laser Interferometer Space Antenna (LISA) will detect gravitational wave signals from coalescing pairs of massive black holes in the total mass range (10^5 - 10^7)/Msol out to cosmological distances. Identifying and monitoring the electromagnetic counterparts of these events would enable cosmological studies and offer new probes of gas physics around well-characterized massive black holes. Milosavljevic & Phinney (2005) proposed that a circumbinary disk around a binary of mass ~10^6 Msol will emit an accretion-powered X-ray afterglow approximately one decade after the gravitational wave event. We revisit this scenario by using Green's function solutions to calculate the temporal viscous evolution and the corresponding electromagnetic signature of the circumbinary disk. Our calculations suggest that an electromagnetic counterpart may become observable as a rapidly brightening source soon after the merger, i.e. several years earlier than previously thought. The afterglow can reach super-Eddington luminosities without violating the local Eddington flux limit. It is emitted in the soft X-ray by the innermost circumbinary disk, but it may be partially reprocessed at optical and infrared frequencies. We also find that the spreading disk becomes increasingly geometrically thick close to the central object as it evolves, indicating that the innermost flow could become advective and radiatively inefficient, and generate a powerful outflow. We conclude that the mergers of massive black holes detected by LISA offer unique opportunities for monitoring on humanly tractable timescales the viscous evolution of accretion flows and the emergence of outflows around massive black holes with precisely known masses, spins and orientations.

preprint2010arXiv

White dwarfs stripped by massive black holes: sources of coincident gravitational and electromagnetic radiation

White dwarfs inspiraling into black holes of mass $\MBH\simgt 10^5M_\odot$ are detectable sources of gravitational waves in the LISA band. In many of these events, the white dwarf begins to lose mass during the main observational phase of the inspiral. The mass loss starts gently and can last for thousands of orbits. The white dwarf matter overflows the Roche lobe through the $L_1$ point at each pericenter passage and the mass loss repeats periodically. The process occurs very close to the black hole and the released gas can accrete, creating a bright source of radiation with luminosity close to the Eddington limit, $L\sim 10^{43}$~erg~s$^{-1}$. This class of inspirals offers a promising scenario for dual detections of gravitational waves and electromagnetic radiation.

preprint2010arXiv

Witnessing the Birth of a Quasar

The coalescence of a supermassive black hole binary (SMBHB) is thought to be accompanied by an electromagnetic (EM) afterglow, produced by the viscous infall of the surrounding circumbinary gas disk after the merger. It has been proposed that once the merger has been detected in gravitational waves (GWs) by LISA, follow-up EM searches for this afterglow can help identify the EM counterpart of the LISA source. Here we study whether the afterglows may be sufficiently bright and numerous to be detectable in EM surveys alone. The viscous afterglow, which lasts for years to decades for SMBHBs in LISA's sensitivity window, is characterized by rapid increases in both the bolometric luminosity and in the spectral hardness of the source. If quasar activity is triggered by the same major galaxy mergers that produce SMBHBs, then the afterglow could be interpreted as a signature of the birth of a quasar. Using an idealized model for the post-merger viscous spreading of the circumbinary disk and the resulting light curve, and using the observed luminosity function of quasars as a proxy for the SMBHB merger rate, we delineate the survey requirements for identifying such birthing quasars. If circumbinary disks have a high disk surface density and viscosity, an all-sky soft X-ray survey with a sensitivity of ~<3x10^-14 erg s^-1 cm^-2 and a time resolution of ~months could identify dozens of birthing quasars with sustained brightening rates of >10%/yr. If >1% of the X-ray emission is reprocessed into optical frequencies, birthing quasars could also be identified in optical transient surveys such as the LSST. Distinguishing a birthing quasar from other variable sources may be facilitated by the monotonic hardening of its spectrum, but will likely remain challenging. This reinforces the notion that joint EM-plus-GW observations offer the best prospects for identifying the EM signatures of SMBHB mergers.

preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

preprint2009arXiv

Rayleigh Adjustment of Narrow Barriers in Protoplanetary Discs

Sharp density features in protoplanetary discs, for instance at the edge of a magnetically dead zone, have recently been proposed as effective barriers to slow down or even stop the problematically fast migration of planetary cores into their central star. Density features on a radial scale approaching the disc vertical scale height might not exist, however, since they could be Rayleigh (or more generally Solberg-Hoiland) unstable. Stability must be checked explicitly in one-dimensional viscous accretion disc models because these instabilities are artificially eliminated in the process of reducing the full set of axisymmetric equations. The disc thermodynamics, via the entropy stratification, and its vertical structure also influence stability when sharp density features are present. We propose the concept of Rayleigh adjustment for viscous disc models: any density feature that violates Rayleigh stability (or its generalization) should be diffused radially by hydrodynamical turbulence on a dynamical time-scale, approaching marginal stability in a quasi-continuous manner.

preprint2008arXiv

Habitable Climates: The Influence of Obliquity

Extrasolar terrestrial planets with the potential to host life might have large obliquities or be subject to strong obliquity variations. We revisit the habitability of oblique planets with an energy balance climate model (EBM) allowing for dynamical transitions to ice-covered snowball states as a result of ice-albedo feedback. Despite the great simplicity of our EBM, it captures reasonably well the seasonal cycle of global energetic fluxes at Earth's surface. It also performs satisfactorily against a full-physics climate model of a highly oblique Earth-like planet, in an unusual regime of circulation dominated by heat transport from the poles to the equator. Climates on oblique terrestrial planets can violate global radiative balance through much of their seasonal cycle, which limits the usefulness of simple radiative equilibrium arguments. High obliquity planets have severe climates, with large amplitude seasonal variations, but they are not necessarily more prone to global snowball transitions than low obliquity planets. We find that terrestrial planets with massive CO2 atmospheres, typically expected in the outer regions of habitable zones, can also be subject to such dynamical snowball transitions. Some of the snowball climates investigated for CO2-rich atmospheres experience partial atmospheric collapse. Since long-term CO2 atmospheric build-up acts as a climatic thermostat for habitable planets, partial CO2 collapse could limit the habitability of such planets. A terrestrial planet's habitability may thus depend sensitively on its short-term climatic stability.

preprint2007arXiv

Probing Gravity with Spacetime Sirens

A gravitational observatory such as LISA will detect coalescing pairs of massive black holes, accurately measure their luminosity distance and help identify a host galaxy or an electromagnetic counterpart. If dark energy is a manifestation of modified gravity on large scales, gravitational waves from cosmologically-distant spacetime sirens are direct probes of this new physics. For example, a gravitational Hubble diagram based on black hole pair luminosity distances and host galaxy redshifts could reveal a large distance extra-dimensional leakage of gravity. Various additional signatures may be expected in a gravitational signal propagated over cosmological scales.

preprint2007arXiv

Toward Eclipse Mapping of Hot Jupiters

Recent Spitzer infrared measurements of hot Jupiter eclipses suggest that eclipse mapping techniques could be used to spatially resolve the day-side photospheric emission of these planets using partial occultations. As a first step in this direction, we simulate ingress/egress lightcurves for the three brightest known eclipsing hot Jupiters and evaluate the degree to which parameterized photospheric emission models can be distinguished from each other with repeated, noisy eclipse measurements. We find that the photometric accuracy of Spitzer is insufficient to use this tool effectively. On the other hand, the level of photospheric details that could be probed with a few JWST eclipse measurements could greatly inform hot Jupiter atmospheric modeling efforts. A JWST program focused on non-parametric eclipse map inversions for hot Jupiters should be actively considered.

preprint2003arXiv

Black Hole Accretion Disks On The Edge

The local axisymmetric stability of hydrodynamical and magnetized, nearly-Keplerian gaseous accretion disks around non-rotating black holes is examined in the vicinity of the classical marginally-stable orbit (at radii ~ R_ms). An approximate Paczynski-Wiita pseudo-Newtonian potential is used. Hydrodynamical disks are linearly unstable inside a radius which differs slightly from the classical R_ms value because of finite pressure and radial stratification effects. Linear stresses associated with unstable hydrodynamical modes vanish exactly at the radius of marginal stability and are generally positive inside of that radius. When a magnetic field is introduced, however, the concept of radius of marginal stability becomes largely irrelevant because there are linearly unstable magneto-rotational modes everywhere. Associated linear stresses are positive and continuous across the region of hydrodynamical marginal stability, even for large-scale "hydro-like" modes subject only to weak magnetic tension. This conclusion is valid for arbitrarily thin disks (in ideal MHD) and it does not require a large-scale "radially-connecting" magnetic field. Results on hydrodynamical diskoseismic modes trapped in deep relativistic potential wells should be revised to account for the short-wavelength Alfven-like behavior of inertio-gravity waves in magnetized disks.