-
FargoCPT: A 2D Multi-Physics Code for Simulating the Interaction of Disks with Stars, Planets and Particles
Authors:
Thomas Rometsch,
Lucas M. Jordan,
Tobias W. Moldenhauer,
Dennis Wehner,
Steven Rendon Restrepo,
Tobias W. A. Müller,
Giovanni Picogna,
Wilhelm Kley,
Cornelis P. Dullemond
Abstract:
Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, a…
▽ More
Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, an updated version of FARGO incorporating other previous enhancements to the code, to provide a simulation environment tailored to study interactions between stars, planets, and disks, ensuring accurate representation of planet systems, hydrodynamics, and dust dynamics with a focus on usability. Methods: The radiation-hydrodynamics part of FargoCPT uses a second-order upwind scheme in 2D polar coordinates supporting multiple equations of state, radiation transport, heating and cooling, and self-gravity. Shocks are considered using artificial viscosity. Integration of the N-body system is achieved by leveraging the REBOUND code. The dust module utilizes massless tracer particles, adapting to drag laws for the Stokes and Epstein regimes. Moreover, FargoCPT provides mechanisms to simulate accretion onto the stars and planets. Results: The code has been tested in practice by its use in various publications. Additionally, it comes with an automated test suite to test the physics modules. Conclusions: FargoCPT offers a unique set of simulation capabilities within the current landscape of publicly available planet-disk interaction simulation tools. Its structured interface and underlying technical updates are intended to assist researchers in the ongoing exploration of planet formation.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
Recycling of the first atmospheres of embedded planets: Dependence on core mass and optical depth
Authors:
T. W. Moldenhauer,
R. Kuiper,
W. Kley,
C. W. Ormel
Abstract:
Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypo…
▽ More
Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypothesis and showed that atmosphere-disk recycling is able to fully compensate for radiative cooling and thereby halt Kelvin-Helmholtz contraction to prevent runaway gas accretion.
To understand the parameters that determine the efficiency of atmospheric recycling, we extend our earlier studies by exploring the effects of the core mass, the effect of circumstellar gas on sub-Keplerian orbits (headwind), and the optical depth of the surrounding gas on the recycling timescale. Additionally, we analyze their effects on the size and mass of the forming atmosphere.
For the explored parameter space, all simulations eventually reach an equilibrium where heating due to hydrodynamic recycling fully compensates radiative cooling. In this equilibrium, the atmosphere-to-core mass ratio stays well below $10 \, \%$, preventing the atmosphere from becoming self-gravitating and entering runaway gas accretion. Higher core masses cause the atmosphere to become turbulent, which further enhances recycling.
Even for our highest core mass of $10 \, M_\mathrm{Earth}$, atmosphere-disk recycling is efficient enough to fully compensate for radiative cooling and prevent the atmosphere from becoming self-gravitating. Hence, in-situ formation of hot Jupiters is very unlikely, and migration of gas giants is a key process required to explain their existence. Our findings imply that atmosphere-disk recycling is the most natural explanation for the prevalence of close-in super-Earths and mini-Neptunes.
△ Less
Submitted 24 February, 2022; v1 submitted 23 February, 2022;
originally announced February 2022.
-
Steady State by Recycling prevents Premature Collapse of Protoplanetary Atmospheres
Authors:
T. W. Moldenhauer,
R. Kuiper,
W. Kley,
C. W. Ormel
Abstract:
In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or…
▽ More
In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or cooling) timescale, but it has also been suggested that thermodynamical cooling can be prevented by hydrodynamical planet atmosphere-disk recycling. We investigate the efficacy of the recycling hypothesis in preventing the collapse of the atmosphere, check for the existence of a steady state configuration, and determine the final atmospheric mass to core mass ratio. We use three-dimensional radiation-hydrodynamic simulations to model the formation of planetary proto-atmospheres. Equations are solved in a local frame centered on the planet. Ignoring small oscillations that average to zero over time, the simulations converge to a steady state where the velocity field of the gas becomes constant in time. In a steady state, the energy loss by radiative cooling is fully compensated by the recycling of the low entropy gas in the planetary atmosphere with high entropy gas from the circumstellar disk. For close-in planets, recycling naturally halts the cooling of planetary proto-atmospheres, preventing them from contracting toward the runaway regime and collapsing into gas giants.
△ Less
Submitted 1 February, 2021;
originally announced February 2021.
-
Spirals inside the millimeter cavity of transition disk SR 21
Authors:
G. A. Muro-Arena,
C. Ginski,
C. Dominik,
M. Benisty,
P. Pinilla,
A. J. Bohn,
T. Moldenhauer,
W. Kley,
D. Harsono,
T. Henning,
R. G. van Holstein,
M. Janson,
M. Keppler,
F. Ménard,
L. M. Pérez,
T. Stolker,
M. Tazzari,
M. Villenave,
A. Zurlo,
C. Petit,
F. Rigal,
O. Möller-Nilsson,
M. Llored,
T. Moulin,
P. Rabou
Abstract:
Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In thi…
▽ More
Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11°. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk.
△ Less
Submitted 18 March, 2020;
originally announced March 2020.