-
Empirical stability criteria for 3D hierarchical triple systems I: Circumbinary planets
Authors:
Nikolaos Georgakarakos,
Siegfried Eggl,
Mohamad Ali-Dib,
Ian Dobbs-Dixon
Abstract:
In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We carry out more than 3 10^8 numerical simulations of planets between the size of Mercury and the lower fusion boundary (13 Jupiter masses) which revolve around the center of mass of a stellar binary over long timescales. For the first time, three dimen…
▽ More
In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We carry out more than 3 10^8 numerical simulations of planets between the size of Mercury and the lower fusion boundary (13 Jupiter masses) which revolve around the center of mass of a stellar binary over long timescales. For the first time, three dimensional and eccentric planetary orbits are considered. We explore systems with a variety of binary and planetary mass ratios, binary and planetary eccentricities from 0 to 0.9 and orbital mutual inclinations ranging from 0 to 180 degrees. The simulation time is set to 10^6 planetary orbital periods. We classify the results of those long term numerical integrations into three categories: stable, unstable and mixed. We provide empirical expressions in the form of multidimensional, parameterized fits for the two borders that separate the three dynamical domains . In addition, we train a machine learning model on our data set in order to have an alternative tool of predicting the stability of circumbinary planets. Both the empirical fits and the machine learning model are tested against randomly generated circumbinary systems with very good results regarding the predictions of orbital stability. The empirical formulae are also applied to the Kepler and TESS circumbinary systems, confirming the stability of the planets in these systems. Finally, we present a REST API with a web based application for convenient access of our simulation data set.
△ Less
Submitted 21 April, 2024;
originally announced April 2024.
-
Protoplanet collisions: new scaling laws from SPH simulations
Authors:
Samuele Crespi,
Mohamad Ali-Dib,
Ian Dobbs-Dixon
Abstract:
One common approach for solving collisions between protoplanets in simulations of planet formation is to employ analytical scaling laws. The most widely used one was developed by Leinhardt & Stewart (2012) from a catalog of ~ 180 N-body simulations of rubble-pile collisions. In this work, we use a new catalogue of more than 20,000 SPH simulations to test the validity and the prediction capability…
▽ More
One common approach for solving collisions between protoplanets in simulations of planet formation is to employ analytical scaling laws. The most widely used one was developed by Leinhardt & Stewart (2012) from a catalog of ~ 180 N-body simulations of rubble-pile collisions. In this work, we use a new catalogue of more than 20,000 SPH simulations to test the validity and the prediction capability of Leinhardt & Stewart (2012) scaling laws. We find that these laws overestimate the fragmentation efficiency in the merging regime and they are not able to properly reproduce the collision outcomes in the super-catastrophic regime. In the merging regime, we also notice a significant dependence between the collision outcome, in terms of the largest remnant mass, and the relative mass of the colliding protoplanets. Here, we present a new set of scaling laws that are able to better predict the collision outcome in all regimes and it is also able to reproduce the observed dependence on the mass ratio. We compare our new scaling laws against a machine learning approach and obtain similar prediction efficiency.
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b
Authors:
Taylor J. Bell,
Nicolas Crouzet,
Patricio E. Cubillos,
Laura Kreidberg,
Anjali A. A. Piette,
Michael T. Roman,
Joanna K. Barstow,
Jasmina Blecic,
Ludmila Carone,
Louis-Philippe Coulombe,
Elsa Ducrot,
Mark Hammond,
João M. Mendonça,
Julianne I. Moses,
Vivien Parmentier,
Kevin B. Stevenson,
Lucas Teinturier,
Michael Zhang,
Natalie M. Batalha,
Jacob L. Bean,
Björn Benneke,
Benjamin Charnay,
Katy L. Chubb,
Brice-Olivier Demory,
Peter Gao
, et al. (58 additional authors not shown)
Abstract:
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5…
▽ More
Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions).
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
High-resolution resonant portraits of a single-planet white dwarf system
Authors:
Dimitri Veras,
Nikolaos Georgakarakos,
Ian Dobbs-Dixon
Abstract:
The dynamical excitation of asteroids due to mean motion resonant interactions with planets is enhanced when their parent star leaves the main sequence. However, numerical investigation of resonant outcomes within post-main-sequence simulations is computationally expensive, limiting the extent to which detailed resonant analyses have been performed. Here, we combine the use of a high-performance c…
▽ More
The dynamical excitation of asteroids due to mean motion resonant interactions with planets is enhanced when their parent star leaves the main sequence. However, numerical investigation of resonant outcomes within post-main-sequence simulations is computationally expensive, limiting the extent to which detailed resonant analyses have been performed. Here, we combine the use of a high-performance computer cluster and the general semianalytical libration width formulation of Gallardo et al. (2021) in order to quantify resonant stability, strength and variation instigated by stellar evolution for a single-planet system containing asteroids on both crossing and non-crossing orbits. We find that resonant instability can be accurately bound with only main-sequence values by computing a maximum libration width as a function of asteroid longitude of pericentre. We also quantify the relative efficiency of mean motion resonances of different orders to stabilize versus destabilize asteroid orbits during both the giant branch and white dwarf phases. The 4:1, 3:1 and 2:1 resonances represent efficient polluters of white dwarfs, and even when in the orbit-crossing regime, both the 4:3 and 3:2 resonances can retain small reservoirs of asteroids in stable orbits throughout giant branch and white dwarf evolution. This investigation represents a preliminary step in characterising how simplified extrasolar Kirkwood gap structures evolve beyond the main-sequence.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
GCM-Motivated Multidimensional Temperature Parametrization Scheme for Phasecurve Retrieval
Authors:
Ian Dobbs-Dixon,
Jasmina Blecic
Abstract:
We present a novel physically motivated, parametrized temperature model for phase-curve retrieval, able to self-consistently assess the variation in thermal structure in multidimensions. To develop this approach, we drew motivation from both full three-dimensional general circulation models and analytic formulations, accounting for the dominant dynamical feature of tidally locked planets, the plan…
▽ More
We present a novel physically motivated, parametrized temperature model for phase-curve retrieval, able to self-consistently assess the variation in thermal structure in multidimensions. To develop this approach, we drew motivation from both full three-dimensional general circulation models and analytic formulations, accounting for the dominant dynamical feature of tidally locked planets, the planetary jet. Our formulation shows notable flexibility. It can generate planetary jets of various characteristics and redistribution efficiencies seen in the literature, including both standard eastward and unusual westward offset hotspots, as well as more exotic configurations for potential future observations. In our modeling scheme we utilize a tractable set of parameters efficient enough to enable future Bayesian analysis and, in addition to the resolved temperature structure, we return physical insights not yet derived from retrievals: the amplitude and the phase offset, and the location and the extent of the equatorial jet.
△ Less
Submitted 8 September, 2022; v1 submitted 23 February, 2022;
originally announced February 2022.
-
Protoplanet Collisions: Statistical Properties of Ejecta
Authors:
Samuele Crespi,
Ian Dobbs-Dixon,
Nikolaos Georgakarakos,
Nader Haghighipour,
Thomas I. Maindl,
Christoph M. Schäfer,
Philip Matthias Winter
Abstract:
The last phase of the formation of rocky planets is dominated by collisions among Moon- to Mars-sized planetary embryos. Simulations of this phase need to handle the difficulty of including the post-impact material without saturating the numerical integrator. A common approach is to include the collision-generated material by clustering it into few bodies with the same mass and uniformly scatterin…
▽ More
The last phase of the formation of rocky planets is dominated by collisions among Moon- to Mars-sized planetary embryos. Simulations of this phase need to handle the difficulty of including the post-impact material without saturating the numerical integrator. A common approach is to include the collision-generated material by clustering it into few bodies with the same mass and uniformly scattering them around the collision point. However, this approach oversimplifies the properties of the collision material by neglecting features that can play important roles in the final structure and composition of the system. In this study, we present a statistical analysis of the orbital architecture, mass, and size distributions of the material generated through embryo-embryo collisions and show how they can be used to develop a model that can be directly incorporated into the numerical integrations. For instance, results of our analysis indicate that the masses of the fragments follow an exponential distribution with an exponent of $-2.21\pm0.17$ over the range of $10^{-7}$ to $2\times 10^{-2}$ Earth-masses. The distribution of the post-impact velocities show that a large number of fragments are scattered toward the central star. The latter is a new finding that may be quite relevant to the delivery of material from the outer regions of the asteroid belt to the accretion zones of terrestrial planets. Finally, we present an analytical model for the 2D distribution of fragments that can be directly incorporated into numerical integrations.
△ Less
Submitted 6 October, 2021;
originally announced October 2021.
-
The entry geometry and velocity of planetary debris into the Roche sphere of a white dwarf
Authors:
Dimitri Veras,
Nikolaos Georgakarakos,
Alexander J. Mustill,
Uri Malamud,
Tim Cunningham,
Ian Dobbs-Dixon
Abstract:
Our knowledge of white dwarf planetary systems predominately arises from the region within a few Solar radii of the white dwarfs, where minor planets break up, form rings and discs, and accrete onto the star. The entry location, angle and speed into this Roche sphere has rarely been explored but crucially determines the initial geometry of the debris, accretion rates onto the photosphere, and ulti…
▽ More
Our knowledge of white dwarf planetary systems predominately arises from the region within a few Solar radii of the white dwarfs, where minor planets break up, form rings and discs, and accrete onto the star. The entry location, angle and speed into this Roche sphere has rarely been explored but crucially determines the initial geometry of the debris, accretion rates onto the photosphere, and ultimately the composition of the minor planet. Here we evolve a total of over 10^5 asteroids with single-planet N-body simulations across the giant branch and white dwarf stellar evolution phases to quantify the geometry of asteroid injection into the white dwarf Roche sphere as a function of planetary mass and eccentricity. We find that lower planetary masses increase the extent of anisotropic injection and decrease the probability of head-on (normal to the Roche sphere) encounters. Our results suggest that one can use dynamical activity within the Roche sphere to make inferences about the hidden architectures of these planetary systems.
△ Less
Submitted 8 June, 2021;
originally announced June 2021.
-
An open-source Bayesian atmospheric radiative transfer (BART) code: III. Initialization, atmospheric profile generator, post-processing routines, and application to exoplanet WASP-43b
Authors:
Jasmina Blecic,
Joseph Harrington,
Patricio E. Cubillos,
M. Oliver Bowman,
Patricio Rojo,
Madison Stemm,
Ryan C. Challener,
Michael D. Himes,
Austin J. Foster,
Ian Dobbs-Dixon,
Andrew S. D. Foster,
Nathaniel B. Lust,
Sarah D. Blumenthal,
Dylan Bruce,
Thomas J. Loredo
Abstract:
This and companion papers by Harrington et al. 2021, submitted and Cubillos et al. 2021, submitted describe an open-source retrieval framework, Bayesian Atmospheric Radiative Transfer (BART), available to the community under the reproducible-research license via https://github.com/exosports/BART . BART is a radiative-transfer code (transit, https://github.com/exosports/transit , Rojo 2009, 2009ASP…
▽ More
This and companion papers by Harrington et al. 2021, submitted and Cubillos et al. 2021, submitted describe an open-source retrieval framework, Bayesian Atmospheric Radiative Transfer (BART), available to the community under the reproducible-research license via https://github.com/exosports/BART . BART is a radiative-transfer code (transit, https://github.com/exosports/transit , Rojo 2009, 2009ASPC..420..321R), initialized by the Thermochemical Equilibrium Abundances (TEA, https://github.com/dzesmin/TEA , Blecic et al. 2016, arXiv:1505.06392) code, and driven through the parameter phase space by a differential-evolution Markov-chain Monte Carlo (MC3, https://github.com/pcubillos/mc3 , Cubillos et al. 2017, arXiv:1610.01336) sampler. In this paper we give a brief description of the framework, and its modules that can be used separately for other scientific purposes; outline the retrieval analysis flow; present the initialization routines, describing in detail the atmospheric profile generator and the temperature and species parameterizations; and specify the post-processing routines and outputs, concentrating on the spectrum band integrator, the best-fit model selection, and the contribution functions. We also present an atmospheric analysis of WASP-43b secondary eclipse data obtained from space- and ground-based observations. We compare our results with the results from the literature, and investigate how the inclusion of additional opacity sources influence the best-fit model.
△ Less
Submitted 26 April, 2021;
originally announced April 2021.
-
An Open-Source Bayesian Atmospheric Radiative Transfer (BART) Code: I. Design, Tests, and Application to Exoplanet HD 189733 b
Authors:
Joseph Harrington,
Michael D. Himes,
Patricio E. Cubillos,
Jasmina Blecic,
Patricio M. Rojo,
Ryan C. Challener,
Nate B. Lust,
M. Oliver Bowman,
Sarah D. Blumenthal,
Ian Dobbs-Dixon,
Andrew S. D. Foster,
Austin J. Foster,
M. R. Green,
Thomas J. Loredo,
Kathleen J. McIntyre,
Madison M. Stemm,
David C. Wright
Abstract:
We present the open-source Bayesian Atmospheric Radiative Transfer (BART) retrieval package, which produces estimates and uncertainties for an atmosphere's thermal profile and chemical abundances from observations. Several BART components are also stand-alone packages, including the parallel Multi-Core Markov chain Monte Carlo (MC3), which implements several Bayesian samplers; a line-by-line radia…
▽ More
We present the open-source Bayesian Atmospheric Radiative Transfer (BART) retrieval package, which produces estimates and uncertainties for an atmosphere's thermal profile and chemical abundances from observations. Several BART components are also stand-alone packages, including the parallel Multi-Core Markov chain Monte Carlo (MC3), which implements several Bayesian samplers; a line-by-line radiative-transfer model, transit; a code that calculates Thermochemical Equilibrium Abundances, TEA; and a test suite for verifying radiative-transfer and retrieval codes, BARTTest. The codes are in Python and C. BART and TEA are under a Reproducible Research (RR) license, which requires reviewed-paper authors to publish a compendium of all inputs, codes, and outputs supporting the paper's scientific claims. BART and TEA produce the compendium's content. Otherwise, these codes are under permissive open-source terms, as are MC3 and BARTTest, for any purpose. This paper presents an overview of the code, BARTTest, and an application to eclipse data for exoplanet HD 189733 b. Appendices address RR methodology for accelerating science, a reporting checklist for retrieval papers, the spectral resolution required for synthetic tests, and a derivation of the effective sample size required to estimate any Bayesian posterior distribution to a given precision, which determines how many iterations to run. Paper II, by Cubillos et al., presents the underlying radiative-transfer scheme and an application to transit data for exoplanet HAT-P-11b. Paper III, by Blecic et al., discusses the initialization and post-processing routines, with an application to eclipse data for exoplanet WASP-43b. We invite the community to use and improve BART and its components at http://GitHub.com/ExOSPORTS/BART/.
△ Less
Submitted 26 April, 2021;
originally announced April 2021.
-
Modeling Solar Proton Event-induced Martian Surface Radiation Dose
Authors:
Dimitra Atri,
Caitlin MacArthur,
Ian Dobbs-Dixon
Abstract:
Solar Proton Events (SPEs) can cause abrupt and significant enhancements to the Martian surface radiation dose. Observations of the impact of SPEs on the Martian surface are available from satellites and surface detectors, but the data set is very limited in time, and the energy range is limited in scope, which makes it insufficient to estimate the impact of major events on the Martian surface. On…
▽ More
Solar Proton Events (SPEs) can cause abrupt and significant enhancements to the Martian surface radiation dose. Observations of the impact of SPEs on the Martian surface are available from satellites and surface detectors, but the data set is very limited in time, and the energy range is limited in scope, which makes it insufficient to estimate the impact of major events on the Martian surface. On the other hand, long-term data of SPEs impacting the Earth spanning a large energy range is widely available, and can be used to estimate the impact of major events on Mars on long timescales. Herein, we take major SPEs observed during the past several decades on Earth (1956 - 2014), along with PAMELA observations (2006 - 2014) and use the GEANT4 Monte Carlo code to calculate the Martian surface radiation dose. We study the contribution of proton fluence and spectral shape of events on the surface radiation dose and estimated the impact of possible major SPEs on the Martian surface in the future. These results have major implications for the planned human exploration of Mars. Overall we find that the radiation dose from extreme events can have a significant impact on astronaut health, and in rare, worst case scenarios, the estimated dose can even reach lethal levels.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST
Authors:
Olivia Venot,
Vivien Parmentier,
Jasmina Blecic,
Patricio E. Cubillos,
Ingo P. Waldmann,
Quentin Changeat,
Julianne I. Moses,
Pascal Tremblin,
Nicolas Crouzet,
Peter Gao,
Diana Powell,
Pierre-Olivier Lagage,
Ian Dobbs-Dixon,
Maria E. Steinrueck,
Laura Kreidberg,
Natalie Batalha,
Jacob L. Bean,
Kevin B. Stevenson,
Sarah Casewell,
Ludmila Carone
Abstract:
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, in…
▽ More
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 $μ$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 $μ$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.
△ Less
Submitted 14 January, 2020;
originally announced January 2020.
-
Do the planets in the HD 34445 system really exist?
Authors:
Nikolaos Georgakarakos,
Ian Dobbs-Dixon
Abstract:
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…
▽ More
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.
△ Less
Submitted 11 July, 2019;
originally announced July 2019.
-
Mass Loss from the Exoplanet WASP-12b Inferred from $\textit{Spitzer}$ Phase Curves
Authors:
Taylor J. Bell,
Michael Zhang,
Patricio E. Cubillos,
Lisa Dang,
Luca Fossati,
Kamen O. Todorov,
Nicolas B. Cowan,
Drake Deming,
Robert T. Zellem,
Kevin B. Stevenson,
Ian J. M. Crossfield,
Ian Dobbs-Dixon,
Jonathan J. Fortney,
Heather A. Knutson,
Michael R. Line
Abstract:
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 $μ$m and 4.5 $μ$m $\textit{Spitzer}$ phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from…
▽ More
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 $μ$m and 4.5 $μ$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 $μ$m, not at 3.6 $μ$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σ$ ($46.2^{\circ}$) at 3.6 $μ$m.
△ Less
Submitted 14 September, 2019; v1 submitted 11 June, 2019;
originally announced June 2019.
-
Toward More Reliable Analytic Thermochemical-equilibrium Abundances
Authors:
Patricio Cubillos,
Jasmina Blecic,
Ian Dobbs-Dixon
Abstract:
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 circum…
▽ More
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 .
△ Less
Submitted 31 January, 2019; v1 submitted 11 January, 2019;
originally announced January 2019.
-
Effects of non-Kozai mutual inclinations on two-planet system stability through all phases of stellar evolution
Authors:
Dimitri Veras,
Nikolaos Georgakarakos,
Boris T. Gänsicke,
Ian Dobbs-Dixon
Abstract:
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. W…
▽ More
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.
△ Less
Submitted 4 September, 2018;
originally announced September 2018.
-
Generation of a Circumstellar Gas Disk by Hot Jupiter WASP-12b
Authors:
Alex Debrecht,
Jonathan Carroll-Nellenback,
Adam Frank,
Luca Fossati,
Eric G. Blackman,
Ian Dobbs-Dixon
Abstract:
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, causin…
▽ More
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.
△ Less
Submitted 1 May, 2018;
originally announced May 2018.
-
Giant planets: good neighbors for habitable worlds?
Authors:
Nikolaos Georgakarakos,
Siegfried Eggl,
Ian Dobbs-Dixon
Abstract:
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 characteri…
▽ More
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.
△ Less
Submitted 6 April, 2018;
originally announced April 2018.
-
The Implications of 3D Thermal Structure on 1D Atmospheric Retrieval
Authors:
Jasmina Blecic,
Ian Dobbs-Dixon,
Thomas Greene
Abstract:
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 waveleng…
▽ More
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.
△ Less
Submitted 18 March, 2018;
originally announced March 2018.
-
Detection of a Westward Hotspot Offset in the Atmosphere of a Hot Gas Giant CoRoT-2b
Authors:
Lisa Dang,
Nicolas B. Cowan,
Joel C. Schwartz,
Emily Rauscher,
Michael Zhang,
Heather A. Knutson,
Michael Line,
Ian Dobbs-Dixon,
Drake Deming,
Sudarsan Sundararajan,
Jonathan J. Fortney,
Ming Zhao
Abstract:
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…
▽ More
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.
△ Less
Submitted 19 January, 2018;
originally announced January 2018.
-
Wavelength Does Not Equal Pressure: Vertical Contribution Functions and their Implications for Mapping Hot Jupiters
Authors:
Ian Dobbs-Dixon,
Nicolas B. Cowan
Abstract:
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…
▽ More
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.
△ Less
Submitted 22 November, 2017;
originally announced November 2017.
-
Exocomets in the Proxima Centauri system and their importance for water transport
Authors:
Richard Schwarz,
Akos Bazso,
Nikolaos Georgakarakos,
Birgit Loibnegger,
David Bancelin,
Elke Pilat-Lohinger,
Kristina Kislyakova,
Rudolf Dvorak,
Thomas Maindl,
Ian Dobbs-Dixon
Abstract:
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)…
▽ More
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.
△ Less
Submitted 7 August, 2018; v1 submitted 13 November, 2017;
originally announced November 2017.
-
Dynamic mineral clouds on HD 189733b II. Monte Carlo radiative transfer for 3D cloudy exoplanet atmospheres: combining scattering and emission spectra
Authors:
Elspeth Lee,
Kenneth Wood,
Ian Dobbs-Dixon,
Anna Rice,
Christiane Helling
Abstract:
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, inhomo…
▽ More
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.
△ Less
Submitted 18 January, 2017; v1 submitted 4 January, 2017;
originally announced January 2017.
-
Binary star influence on post-main-sequence multi-planet stability
Authors:
Dimitri Veras,
Nikolaos Georgakarakos,
Ian Dobbs-Dixon,
Boris T. Gaensicke
Abstract:
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) si…
▽ More
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.
△ Less
Submitted 17 October, 2016;
originally announced October 2016.
-
Long term evolution of planetary systems with a terrestrial planet and a giant planet
Authors:
Nikolaos Georgakarakos,
Ian Dobbs-Dixon,
Michael J. Way
Abstract:
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 expan…
▽ More
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.
△ Less
Submitted 9 June, 2016;
originally announced June 2016.
-
Hot planetary winds near a star: dynamics, wind-wind interactions, and observational signatures
Authors:
Jonathan Carroll-Nellenback,
Adam Frank,
Baowei Liu,
Alice C. Quillen,
Eric G. Blackman,
Ian Dobbs-Dixon
Abstract:
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 l…
▽ More
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-α 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.
△ Less
Submitted 28 April, 2016;
originally announced April 2016.
-
Dynamic mineral clouds on HD 189733b I. 3D RHD with kinetic, non-equilibrium cloud formation
Authors:
E. K. H. Lee,
I. Dobbs-Dixon,
Ch. Helling,
K. Bognar,
P. Woitke
Abstract:
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…
▽ More
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 $μ$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.
△ Less
Submitted 18 July, 2016; v1 submitted 30 March, 2016;
originally announced March 2016.
-
The mineral clouds on HD 209458b and HD189733b
Authors:
Ch. Helling,
G. Lee,
I. Dobbs-Dixon,
N. Mayne,
D. S. Amundsen,
J. Khaimova,
A. A. Unger,
J. Manners,
D. Acreman,
C. Smith
Abstract:
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 ch…
▽ More
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.
△ Less
Submitted 23 March, 2016; v1 submitted 13 March, 2016;
originally announced March 2016.
-
Atmospheric electrification in dusty, reactive gases in the solar system and beyond
Authors:
Ch. Helling,
R. G. Harrison,
F. Honary,
D. A. Diver,
K. Aplin,
I. Dobbs-Dixon,
U. Ebert,
S. Inutsuka,
F. J. Gordillo-Vazquez,
S. Littlefair
Abstract:
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 oldes…
▽ More
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.
△ Less
Submitted 29 January, 2016; v1 submitted 18 January, 2016;
originally announced January 2016.
-
Spectral Eclipse Timing
Authors:
Ian Dobbs-Dixon,
Eric Agol,
Drake Deming
Abstract:
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 timi…
▽ More
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.
△ Less
Submitted 19 November, 2015;
originally announced November 2015.
-
Planetary Evaporation and the Dynamics of Planet Wind/Stellar Wind Bow Shocks
Authors:
A. Frank,
B. Liu,
J. Carroll-Nellenback,
A. C. Quillen,
J. F. Kasting,
I. Dobbs-Dixon,
E. G. Blackman
Abstract:
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 stren…
▽ More
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.
△ Less
Submitted 4 August, 2015;
originally announced August 2015.
-
Modelling the local and global cloud formation on HD 189733b
Authors:
E. K. H. Lee,
Ch. Helling,
I. Dobbs-Dixon,
D. Juncher
Abstract:
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 app…
▽ More
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μ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.
△ Less
Submitted 25 May, 2015;
originally announced May 2015.
-
3D Structures of equatorial waves and the resulting superrotation in the atmosphere of a tidally locked hot Jupiter
Authors:
Shang-Min Tsai,
Ian Dobbs-Dixon,
Pin-Gao Gu
Abstract:
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 ext…
▽ More
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.
△ Less
Submitted 20 August, 2014; v1 submitted 30 April, 2014;
originally announced May 2014.
-
Infrared Transmission Spectroscopy of the Exoplanets HD209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope
Authors:
Drake Deming,
Ashlee Wilkins,
Peter McCullough,
Adam Burrows,
Jonathan J. Fortney,
Eric Agol,
Ian Dobbs-Dixon,
Nikku Madhusudhan,
Nicolas Crouzet,
Jean-Michel Desert,
Ronald L. Gilliland,
Korey Haynes,
Heather A. Knutson,
Michael Line,
Zazralt Magic,
Avi M. Mandell,
Sukrit Ranjan,
David Charbonneau,
Mark Clampin,
Sara Seager,
Adam P. Showman
Abstract:
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,…
▽ More
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.
△ Less
Submitted 12 July, 2013; v1 submitted 5 February, 2013;
originally announced February 2013.
-
Stellar irradiated discs and implications on migration of embedded planets I: equilibrium discs
Authors:
Bertram Bitsch,
Aurélien Crida,
Alessandro Morbidelli,
Willy Kley,
Ian Dobbs-Dixon
Abstract:
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 str…
▽ More
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.
△ Less
Submitted 2 December, 2012; v1 submitted 27 November, 2012;
originally announced November 2012.
-
Three Dimensional Radiative Hydrodynamical Simulations of the Highly Irradiated Short Period Exoplanet HD189733b
Authors:
Ian Dobbs-Dixon,
Eric Agol
Abstract:
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…
▽ More
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.
△ Less
Submitted 7 November, 2012;
originally announced November 2012.
-
The Impact of Circumplantary Jets on Transit Spectra and Timing Offsets for Hot-Jupiters
Authors:
Ian Dobbs-Dixon,
Eric Agol,
Adam Burrows
Abstract:
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 betw…
▽ More
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.
△ Less
Submitted 19 October, 2011;
originally announced October 2011.
-
Radiative Hydrodynamic Simulations of HD209458b: Temporal Variability
Authors:
Ian Dobbs-Dixon,
Andrew Cumming,
D. N. C Lin
Abstract:
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…
▽ More
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.
△ Less
Submitted 6 January, 2010;
originally announced January 2010.
-
Radiative Hydrodynamical Studies of Irradiated Atmospheres
Authors:
Ian Dobbs-Dixon
Abstract:
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 radiati…
▽ More
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.
△ Less
Submitted 28 July, 2008;
originally announced July 2008.
-
Atmospheric Dynamics of Short-period Extra Solar Gas Giant Planets I: Dependence of Night-Side Temperature on Opacity
Authors:
Ian Dobbs-Dixon,
D. N. C. Lin
Abstract:
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 contr…
▽ More
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.
△ Less
Submitted 24 April, 2007;
originally announced April 2007.
-
Tidal Barrier and the Asymptotic Mass of Proto Gas-Giant Planets
Authors:
Ian Dobbs-Dixon,
Shu Lin Li,
D. N. C. Lin
Abstract:
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 alon…
▽ More
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.
△ Less
Submitted 9 January, 2007;
originally announced January 2007.
-
Spin-Orbit Evolution Of Short-Period Planets
Authors:
Ian Dobbs-Dixon,
D. N. C. Lin,
Rosemary A. Mardling
Abstract:
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 discrepanc…
▽ More
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.
△ Less
Submitted 10 August, 2004;
originally announced August 2004.