-
Blowin' in the non-isothermal wind: core-powered mass loss with hydrodynamic radiative transfer
Authors:
William Misener,
Matthäus Schulik,
Hilke E. Schlichting,
James E. Owen
Abstract:
The mass loss rates of planets undergoing core-powered escape are usually modeled using an isothermal Parker-type wind at the equilibrium temperature, $T_\mathrm{eq}$. However, the upper atmospheres of sub-Neptunes may not be isothermal if there are significant differences between the opacity to incident visible and outgoing infrared radiation. We model bolometrically-driven escape using aiolos, a…
▽ More
The mass loss rates of planets undergoing core-powered escape are usually modeled using an isothermal Parker-type wind at the equilibrium temperature, $T_\mathrm{eq}$. However, the upper atmospheres of sub-Neptunes may not be isothermal if there are significant differences between the opacity to incident visible and outgoing infrared radiation. We model bolometrically-driven escape using aiolos, a hydrodynamic radiative-transfer code that incorporates double-gray opacities, to investigate the process's dependence on the visible-to-infrared opacity ratio, $γ$. For a value of $γ\approx 1$, we find that the resulting mass loss rates are well-approximated by a Parker-type wind with an isothermal temperature $T = T_\mathrm{eq}/2^{1/4}$. However, we show that over a range of physically plausible values of $γ$, the mass loss rates can vary by orders of magnitude, ranging from $10^{-5} \times$ the isothermal rate for low $γ$ to $10^5 \times$ the isothermal rate for high $γ$. The differences in mass loss rates are largest for small planet radii, while for large planet radii, mass loss rates become nearly independent of $γ$ and approach the isothermal approximation. We incorporate these opacity-dependent mass loss rates into a self-consistent planetary mass and energy evolution model and show that lower/higher $γ$ values lead to more/less hydrogen being retained after core-powered mass loss. In some cases, the choice of opacities determines whether or not a planet can retain a significant primordial hydrogen atmosphere. The dependence of escape rate on the opacity ratio may allow atmospheric escape observations to directly constrain a planet's opacities and therefore its atmospheric composition.
△ Less
Submitted 24 May, 2024;
originally announced May 2024.
-
Aiolos- A multi-purpose 1-D hydrodynamics code for planetary atmospheres
Authors:
Matthäus Schulik,
Richard Booth
Abstract:
We present a new 1-D multi-physics simulation code with use cases intended for, but not limited to, hydrodynamic escapeproblems of planetary atmospheres and planetary accretion models. Our formulation treats an arbitrary number of species asseparated hydrodynamic fields, couples them via friction laws, allows for a multi-band flux-limited radiation transport, and tracksionization fronts in high-en…
▽ More
We present a new 1-D multi-physics simulation code with use cases intended for, but not limited to, hydrodynamic escapeproblems of planetary atmospheres and planetary accretion models. Our formulation treats an arbitrary number of species asseparated hydrodynamic fields, couples them via friction laws, allows for a multi-band flux-limited radiation transport, and tracksionization fronts in high-energy irradiation bands. Besides coupling various known numerical solution techniques together, weimprove on the numerical stability of deep hydrostatic atmospheres by using a well-balanced scheme, hence preventing unphysicaldriving of atmospheric in- or outflow. We demonstrate the correct physical behaviour of the individual code modules and presenta few simple, new applications, such as a proof-of-concept simulations of combined core-powered mass-loss and UV-drivenatmospheric escape, along with a fully time-dependent core-collapse giant planet simulation. The multi-species nature of thecode opens up the area of exploring simulations that are agnostic towards the dominant atmospheric species and can lead toimplementations of advanced planetary evolution schemes.
△ Less
Submitted 14 July, 2022;
originally announced July 2022.
-
Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions
Authors:
Paul Mollière,
Tamara Molyarova,
Bertram Bitsch,
Thomas Henning,
Aaron Schneider,
Laura Kreidberg,
Christian Eistrup,
Remo Burn,
Evert Nasedkin,
Dmitry Semenov,
Christoph Mordasini,
Martin Schlecker,
Kamber R. Schwarz,
Sylvestre Lacour,
Mathias Nowak,
Matthäus Schulik
Abstract:
Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric comp…
▽ More
Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations.
△ Less
Submitted 28 April, 2022;
originally announced April 2022.
-
Dust formation in the outflows of catastrophically evaporating planets
Authors:
Richard A. Booth,
James E. Owen,
Matthäus Schulik
Abstract:
Ultra-short period planets offer a window into the poorly understood interior composition of exoplanets through material evaporated from their rocky interiors. Among these objects are a class of disintegrating planets, observed when their dusty tails transit in front of their host stars. These dusty tails are thought to originate from dust condensation in thermally-driven winds emanating from the…
▽ More
Ultra-short period planets offer a window into the poorly understood interior composition of exoplanets through material evaporated from their rocky interiors. Among these objects are a class of disintegrating planets, observed when their dusty tails transit in front of their host stars. These dusty tails are thought to originate from dust condensation in thermally-driven winds emanating from the sublimating surfaces of these planets. Existing models of these winds have been unable to explain their highly variable nature and have not explicitly modelled how dust forms in the wind. Here we present new radiation-hydrodynamic simulations of the winds from these planets, including a minimal model for the formation and destruction of dust, assuming that nucleation can readily take place. We find that dust forms readily in the winds, a consequence of large dust grains obtaining lower temperatures than the planet's surface. As hypothesised previously, we find that the coupling of the planet's surface temperature to the outflow properties via the dust's opacity can drive time-variable flows when dust condensation is sufficiently fast. In agreement with previous work, our models suggest that these dusty tails are a signature of catastrophically evaporating planets that are close to the end of their lives. Finally, we discuss the implications of our results for the dust's composition. More detailed hydrodynamic models that self-consistently compute the nucleation and composition of the dust and gas are warranted in order to use these models to study the planet's interior composition.
△ Less
Submitted 26 October, 2022; v1 submitted 27 April, 2022;
originally announced April 2022.
-
On the structure and mass delivery towards circumplanetary discs
Authors:
Matthäus Schulik,
Anders Johansen,
Bertram Bitsch,
Elena Lega,
Michiel Lambrechts
Abstract:
Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FA…
▽ More
Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FARGOCA, with a grid structure that allows resolving the planetary gravitational potential sufficiently well for a CPD to form. We then study the gas flows and density/temperature structures that emerge as a function of planet mass, opacity and potential depth. Our results indicate interesting structure formation for Jupiter-mass planets at low opacities, which we subsequently analyse in detail. Using an opacity level that is 100 times lower than that of ISM dust, our Jupiter-mass protoplanet features an envelope that is sufficiently cold for a CPD to form, and a free-fall region separating the CPD and the circumstellar disc emerges. Interestingly, this free-fall region appears to be a result of supersonic erosion of outer envelope material, as opposed to the static structure formation that one would expect at low opacities. Our analysis reveals that the planetary spiral arms seem to pose a significant pressure barrier that needs to be overcome through radiative cooling in order for gas to free-fall onto the CPD. The circulation inside the CPD is near-keplerian and modified by the presence of CPD spiral arms. For high opacities we recover results from the literature, finding an essentially featureless hot envelope. With this work, we demonstrate the first simulation and analysis of a complete detachment process of a protoplanet from its parent disc in a 3D radiation hydrodynamics setting.
△ Less
Submitted 30 March, 2020;
originally announced March 2020.
-
Global 3D radiation-hydrodynamic simulations of gas accretion: The opacity dependent growth of Saturn-mass planets
Authors:
Matthäus Schulik,
Anders Johansen,
Bertram Bitsch,
Elena Lega
Abstract:
The full spatial structure and temporal evolution of the accretion flow into the envelopes of growing gas giants in their nascent discs is only accessible in simulations. Such simulations are constrained in their approach of computing the formation of gas giants by dimensionality, resolution, consideration of self-gravity, energy treatment and the adopted opacity law. Our study explores how a numb…
▽ More
The full spatial structure and temporal evolution of the accretion flow into the envelopes of growing gas giants in their nascent discs is only accessible in simulations. Such simulations are constrained in their approach of computing the formation of gas giants by dimensionality, resolution, consideration of self-gravity, energy treatment and the adopted opacity law. Our study explores how a number of these parameters affect that measured accretion rate of a Saturn-mass planet. We present a global 3D radiative hydrodynamics framework using the FARGOCA-code. The planet is represented by a gravitational potential with a smoothing length at the location of the planet. No mass or energy sink is used, instead luminosity and gas accretion rates are self-consistently computed. We find that the gravitational smoothing length must be resolved by at least 10 grid cells to obtain converged measurements of the gas accretion rates. Secondly, we find gas accretion rates into planetary envelopes compatible with previous studies, and continue to explain those via the structure of our planetary envelopes and their luminosities. Our measured gas accretion rates are formally in the stage of Kelvin-Helmholtz contraction due to the modest entropy loss that can be obtained over the simulation time-scale, but our accretion rates are compatible with those expected during late run-away accretion. Our detailed simulations of the gas flow into the envelope of a Saturn-mass planet provides a framework for understanding the general problem of gas accretion during planet formation and highlights circulation features that develop inside the planetary envelopes. Those circulation features feedback into the envelope energetics and can have further implications for transporting dust into the inner regions of the envelope.
△ Less
Submitted 22 September, 2019; v1 submitted 18 September, 2019;
originally announced September 2019.