-
On the suppression of giant planet formation around low-mass stars in clustered environments
Authors:
Shuo Huang,
Simon Portegies Zwart,
Maite J. C. Wilhelm
Abstract:
Context: Current exoplanet formation studies tend to overlook the birth environment of stars in clustered environments. The effect of this environment on the planet-formation process, however, is important, especially in the earliest stage. Aims: We investigate the differences in planet populations forming in star-cluster environments through pebble accretion and compare these results with the pla…
▽ More
Context: Current exoplanet formation studies tend to overlook the birth environment of stars in clustered environments. The effect of this environment on the planet-formation process, however, is important, especially in the earliest stage. Aims: We investigate the differences in planet populations forming in star-cluster environments through pebble accretion and compare these results with the planet formation around isolated stars. We try to provide potential signatures on the young planetary systems to guide future observation. Methods: We design and present a new planet population synthesis code for clustered environments. The planet formation model is based on pebble accretion and includes migration in the circumstellar disk. The disk's gas and dust are evolved in 1D simulations considering the effects of photo-evaporation of the nearby stars. Results: Planetary systems in a clustered environment are different than those born in isolation; the environmental effects are important for a wide range of observable parameters and the eventual architecture of the planetary systems. Planetary systems born in a clustered environment lack cold Jupiters compared to isolated planetary systems. This effect is more pronounced for low-mass stars ($\lesssim$0.2 $M_\odot$). On the other hand, planetary systems born in clusters show an excess of cold Neptune around these low-mass stars. Conclusions: In future observations, finding an excess of cold Neptunes and a lack of cold Jupiters could be used to constrain the birth environments of these planetary systems. Exploring the dependence of cold Jupiter's intrinsic occurrence rate on stellar mass provides insights into the birth environment of their proto-embryos.
△ Less
Submitted 2 August, 2024; v1 submitted 26 July, 2024;
originally announced July 2024.
-
Massive Star Cluster Formation I. High Star Formation Efficiency While Resolving Feedback of Individual Stars
Authors:
Brooke Polak,
Mordecai-Mark Mac Low,
Ralf S. Klessen,
Jia Wei Teh,
Claude Cournoyer-Cloutier,
Eric P. Andersson,
Sabrina M. Appel,
Aaron Tran,
Sean C. Lewis,
Maite J. C. Wilhelm,
Simon Portegies Zwart,
Simon C. O. Glover,
Long Wang,
Stephen L. W. McMillan
Abstract:
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upg…
▽ More
The mode of star formation that results in the formation of globular clusters and young massive clusters is difficult to constrain through observations. We present models of massive star cluster formation using the Torch framework, which uses AMUSE to couple distinct multi-physics codes that handle star formation, stellar evolution and dynamics, radiative transfer, and magnetohydrodynamics. We upgrade Torch by implementing the N-body code PeTar, thereby enabling Torch to handle massive clusters forming from $10^6\rm\, M_\odot$ clouds with $\ge10^5$ individual stars. We present results from Torch simulations of star clusters forming from $10^4, 10^5$, and $10^6\rm M_\odot$ turbulent, spherical gas clouds (named M4, M5, M6) of radius $R=11.7$ pc. We find that star formation is highly efficient and becomes more so at higher cloud mass and surface density. For M4, M5, and M6 with initial surface densities $2.325\times 10^{1,2,3}\rm\, M_\odot\, pc^{-2}$, after a free-fall time of $t_{ff}=6.7,2.1,0.67$ Myr, we find that $\sim30\%$, 40%, and 60% of the cloud mass has formed into stars, respectively. The final integrated star formation efficiency is $32\%,\, 65\%$, and 85\% for M4, M5, and M6. Observations of nearby clusters similar to M4 have similar integrated star formation efficiencies of $\leq30\%$. The M5 and M6 models represent a different regime of cluster formation that is more appropriate for the conditions in starburst galaxies and gas-rich galaxies at high redshift, and that leads to a significantly higher efficiency of star formation. We argue that young massive clusters build up through short efficient bursts of star formation in regions that are sufficiently dense ($\ge 10^2 \rm\,M_\odot\,pc^{-2}$) and massive ($\ge10^5\rm\, M_\odot$). In such environments, the dynamical time of the cloud becomes short enough that stellar feedback cannot act quickly enough to slow star formation.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Early Evolution and 3D Structure of Embedded Star Clusters
Authors:
Claude Cournoyer-Cloutier,
Alison Sills,
William E. Harris,
Sabrina M. Appel,
Sean C. Lewis,
Brooke Polak,
Aaron Tran,
Martijn J. C. Wilhelm,
Mordecai-Mark Mac Low,
Stephen L. W. McMillan,
Simon Portegies Zwart
Abstract:
We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at…
▽ More
We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the AMUSE framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in FLASH. We simulate a suite of $10^4$ M$_{\odot}$ clouds at 0.0683 pc resolution for $\sim$ 2 Myr after the onset of star formation, with virial parameters $α$ = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial $α$ having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over $\sim$ 0.01 Myr timescales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses $\lesssim$ 1000 M$_{\odot}$ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation.
△ Less
Submitted 16 February, 2023;
originally announced February 2023.
-
Radiation shielding of protoplanetary discs in young star-forming regions
Authors:
Martijn J. C. Wilhelm,
Simon Portegies Zwart,
Claude Cournoyer-Cloutier,
Sean C. Lewis,
Brooke Polak,
Aaron Tran,
Mordecai-Mark Mac Low
Abstract:
Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two en…
▽ More
Protoplanetary discs spend their lives in the dense environment of a star forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the AMUSE framework to couple the Torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Early-Forming Massive Stars Suppress Star Formation and Hierarchical Cluster Assembly
Authors:
Sean C. Lewis,
Stephen L. W. McMillan,
Mordecai-Mark Mac Low,
Claude Cournoyer-Cloutier,
Brooke Polak,
Martijn J. C. Wilhelm,
Aaron Tran,
Alison Sills,
Simon Portegies Zwart,
Ralf S. Klessen,
Joshua E. Wall
Abstract:
Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation…
▽ More
Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suite \texttt{Torch}, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer, $N$-body dynamics, and stellar feedback to model four initially identical $10^4$ M$_\odot$ giant molecular clouds with a Gaussian density profile peaking at $521.5 \mbox{ cm}^{-3}$. Using the \texttt{Torch} software suite through the \texttt{AMUSE} framework we modify three of the models to ensure that the first star that forms is very massive (50, 70, 100 M$_\odot$). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star forming region. The star formation rate is suppressed, reducing the total mass of stars formed. Our fiducial control model without an early massive star has a larger star formation rate and total efficiency by up to a factor of three and a higher average star formation efficiency per free-fall time by up to a factor of seven. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster.
△ Less
Submitted 28 February, 2023; v1 submitted 2 December, 2022;
originally announced December 2022.
-
Expanding shells around young clusters -- S 171/Be 59
Authors:
G. F. Gahm,
M. J. C. Wilhelm,
C. M. Persson,
A. A. Djupvik,
S. F. Portegies Zwart
Abstract:
Some HII regions that surround young stellar clusters are bordered by molecular shells that appear to expand at a rate inconsistent with our current model simulations. In this study we focus on the dynamics of Sharpless 171 (including NGC 7822), which surrounds the cluster Berkeley 59. We aim to compare the velocity pattern over the molecular shell with the mean radial velocity of the cluster for…
▽ More
Some HII regions that surround young stellar clusters are bordered by molecular shells that appear to expand at a rate inconsistent with our current model simulations. In this study we focus on the dynamics of Sharpless 171 (including NGC 7822), which surrounds the cluster Berkeley 59. We aim to compare the velocity pattern over the molecular shell with the mean radial velocity of the cluster for estimates of the expansion velocities of different shell structures, and to match the observed properties with model simulations. Optical spectra of 27 stars located in Berkeley 59 were collected at the Nordic Optical Telescope, and a number of molecular structures scattered over the entire region were mapped in $^{13}$CO(1-0) at Onsala Space Observatory. We obtained radial velocities and MK classes for the cluster's stars. At least four of the O stars are found to be spectroscopic binaries, in addition to one triplet system. From these data we obtain the mean radial velocity of the cluster. From the $^{13}$CO spectra we identify three shell structures, expanding relative to the cluster at moderate velocity (4 km/s), high velocity (12 km/s), and in between. The high-velocity cloudlets extend over a larger radius and are less massive than the low-velocity cloudlets. We performed a model simulation to understand the evolution of this complex. Our simulation of the Sharpless 171 complex and Berkeley 59 cluster demonstrates that the individual components can be explained as a shell driven by stellar winds from the massive cluster members. However, our relatively simple model produces a single component. Modelling of the propagation of shell fragments through a uniform interstellar medium demonstrates that dense cloudlets detached from the shell are decelerated less efficiently than the shell itself. They can reach greater distances and retain higher velocities than the shell.
△ Less
Submitted 11 May, 2022;
originally announced May 2022.
-
The Ice Coverage of Earth-like Planets Orbiting FGK Stars
Authors:
Caitlyn Wilhelm,
Rory Barnes,
Russell Deitrick,
Rachel Mellman
Abstract:
The photometric and spectroscopic signatures of habitable planets orbiting FGK stars may be modulated by surface ice coverage. To estimate its frequency and locations, we simulated the climates of hypothetical planets with a 1D energy balance model and assumed that the planets possess properties similar to modern Earth (mass, geography, atmosphere). We first simulated planets with fixed rotational…
▽ More
The photometric and spectroscopic signatures of habitable planets orbiting FGK stars may be modulated by surface ice coverage. To estimate its frequency and locations, we simulated the climates of hypothetical planets with a 1D energy balance model and assumed that the planets possess properties similar to modern Earth (mass, geography, atmosphere). We first simulated planets with fixed rotational axes and circular orbits, finding that the vast majority (>90%) of planets with habitable surfaces are free of ice. For planets with partial ice coverage, the parameter space for ice caps (interannual ice located at the poles) is about as large as that for "ice belts" (interannual ice located at the equator), but belts only persist on land. We then performed simulations that mimicked perturbations from other planets by forcing sinusoidal orbital and rotational oscillations over a range of frequencies and amplitudes. We assume initially ice-free surfaces and set the initial eccentricity distribution to mirror known exoplanets, while the initial obliquity distribution matches planet formation predictions, ie favoring 90 degrees. For these dynamic cases, we find again that ~90% of habitable planets are free of surface ice for a range of assumptions for ice's albedo. Planets orbiting F dwarfs are three times as likely to have ice caps than belts, but for planets orbiting K and G dwarfs ice belts are twice as likely as caps. In some cases, a planet's surface ice can cycle between the equatorial and polar regions. Future direct imaging surveys of habitable planets may be able to test these predictions.
△ Less
Submitted 6 December, 2021;
originally announced December 2021.
-
Exploring the possibility of Peter Pan discs across stellar mass
Authors:
Martijn J. C. Wilhelm,
Simon Portegies Zwart
Abstract:
Recently, several accreting M dwarf stars have been discovered with ages far exceeding the typical protoplanetary disc lifetime. These `Peter Pan discs' can be explained as primordial discs that evolve in a low-radiation environment. The persistently low masses of the host stars raise the question whether primordial discs can survive up to these ages around stars of higher mass. In this work we ex…
▽ More
Recently, several accreting M dwarf stars have been discovered with ages far exceeding the typical protoplanetary disc lifetime. These `Peter Pan discs' can be explained as primordial discs that evolve in a low-radiation environment. The persistently low masses of the host stars raise the question whether primordial discs can survive up to these ages around stars of higher mass. In this work we explore the way in which different mass loss processes in protoplanetary discs limit their maximum lifetimes, and how this depends on host star mass. We find that stars with masses $\lesssim$ 0.6 M$_\odot$ can retain primordial discs for $\sim$50 Myr. At stellar masses $\gtrsim$ 0.8 M$_\odot$, the maximum disc lifetime decreases strongly to below 50 Myr due to relatively more efficient accretion and photoevaporation by the host star. Lifetimes up to 15 Myr are still possible for all host star masses up to $\sim$2 M$_\odot$. For host star masses between 0.6 and 0.8 M$_\odot$, accretion ceases and an inner gap forms before 50 Myr in our models. Observations suggest that such a configuration is rapidly dispersed. We conclude that Peter Pan discs can only occur around M dwarf stars.
△ Less
Submitted 3 September, 2021;
originally announced September 2021.
-
Evolution of circumstellar discs in young star-forming regions
Authors:
Francisca Concha-Ramírez,
Martijn J. C Wilhelm,
Simon Portegies Zwart
Abstract:
The evolution of circumstellar discs is influenced by their surroundings. The relevant processes include external photoevaporation due to nearby stars, and dynamical truncations. The impact of these processes on disc populations depends on the star-formation history and on the dynamical evolution of the region. Since star formation history and the phase-space characteristics of the stars are impor…
▽ More
The evolution of circumstellar discs is influenced by their surroundings. The relevant processes include external photoevaporation due to nearby stars, and dynamical truncations. The impact of these processes on disc populations depends on the star-formation history and on the dynamical evolution of the region. Since star formation history and the phase-space characteristics of the stars are important for the evolution of the discs, we start simulating the evolution of the star cluster with the results of molecular cloud collapse simulations. In the simulation we form stars with circumstellar discs, which can be affected by different processes. Our models account for the viscous evolution of the discs, internal and external photoevaporation of gas, external photoevaporation of dust, and dynamical truncations. All these processes are resolved together with the dynamical evolution of the cluster, and the evolution of the stars.
An extended period of star formation, lasting for at least 2 Myr, results in some discs being formed late. These late formed discs have a better chance of survival because the cluster gradually expands with time, and a lower local stellar density reduces the effects of photoevaporation and dynamical truncation. Late formed discs can then be present in regions of high UV radiation, solving the proplyd lifetime problem. We also find a considerable fraction of discs that lose their gas content, but remain sufficiently rich in solids to be able to form a rocky planetary system.
△ Less
Submitted 26 May, 2022; v1 submitted 19 January, 2021;
originally announced January 2021.
-
Effects of stellar density on the photoevaporation of circumstellar discs
Authors:
Francisca Concha-Ramírez,
Martijn J. C. Wilhelm,
Simon Portegies Zwart,
Sierk E. van Terwisga,
Alvaro Hacar
Abstract:
Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the…
▽ More
Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the most important ones. Theoretical and observational evidence point to circumstellar discs losing mass quickly when in the vicinity of massive, bright stars. In this work we simulate circumstellar discs in clustered environments in a range of stellar densities, where the photoevaporation mass-loss process is resolved simultaneously with the stellar dynamics, stellar evolution, and the viscous evolution of the discs. Our results indicate that external photoevaporation is efficient in depleting disc masses and that the degree of its effect is related to stellar density. We find that a local stellar density lower than 100 stars pc$^{-2}$ is necessary for discs massive enough to form planets to survive for \SI{2.0}{Myr}. There is an order of magnitude difference in the disc masses in regions of projected density 100 stars pc$^{-2}$ versus $10^4$ stars pc$^{-2}$. We compare our results to observations of the Lupus clouds, the Orion Nebula Cluster, the Orion Molecular Cloud-2, Taurus, and NGC 2024, and find that the trends observed between region density and disc masses are similar to those in our simulations.
△ Less
Submitted 20 November, 2020; v1 submitted 12 June, 2020;
originally announced June 2020.
-
The Milky Way's bar structural properties from gravitational waves
Authors:
Martijn J. C. Wilhelm,
Valeriya Korol,
Elena M. Rossi,
Elena D'Onghia
Abstract:
The Laser Interferometer Space Antenna (LISA) will enable Galactic gravitational wave (GW) astronomy by individually resolving $ > 10^4$ signals from double white dwarf (DWD) binaries throughout the Milky Way. In this work we assess for the first time the potential of LISA data to map the Galactic stellar bar and spiral arms, since GWs are unaffected by stellar crowding and dust extinction unlike…
▽ More
The Laser Interferometer Space Antenna (LISA) will enable Galactic gravitational wave (GW) astronomy by individually resolving $ > 10^4$ signals from double white dwarf (DWD) binaries throughout the Milky Way. In this work we assess for the first time the potential of LISA data to map the Galactic stellar bar and spiral arms, since GWs are unaffected by stellar crowding and dust extinction unlike optical observations of the bulge region. To achieve this goal we combine a realistic population of Galactic DWDs with a high-resolution N-Body simulation a galaxy in good agreement with the Milky Way. We then model GW signals from our synthetic DWD population and reconstruct the structure of the simulated Galaxy from mock LISA observations. Our results show that while the low signal contrast between the background disc and the spiral arms hampers our ability to characterise the spiral structure, the stellar bar will instead clearly appear in the GW map of the bulge. The bar length and bar width derived from these synthetic observations are underestimated, respectively within $1σ$ and at a level greater than $2σ$, but the resulting axis ratio agrees to well within $1σ$, while the viewing angle is recovered to within one degree. These are competitive constraints compared to those from electromagnetic tracers, and they are obtained with a completely independent method. We therefore foresee that the synergistic use of GWs and electromagnetic tracers will be a powerful strategy to map the bar and the bulge of the Milky Way.
△ Less
Submitted 3 November, 2020; v1 submitted 24 March, 2020;
originally announced March 2020.
-
External photoevaporation of circumstellar disks constrains the timescale for planet formation
Authors:
Francisca Concha-Ramírez,
Martijn J. C. Wilhelm,
Simon Portegies Zwart,
Thomas J. Haworth
Abstract:
Planet-forming circumstellar disks are a fundamental part of the star formation process. Since stars form in a hierarchical fashion in groups of up to hundreds or thousands, the UV radiation environment that these disks are exposed to can vary in strength by at least six orders of magnitude. This radiation can limit the masses and sizes of the disks. Diversity in star forming environments can have…
▽ More
Planet-forming circumstellar disks are a fundamental part of the star formation process. Since stars form in a hierarchical fashion in groups of up to hundreds or thousands, the UV radiation environment that these disks are exposed to can vary in strength by at least six orders of magnitude. This radiation can limit the masses and sizes of the disks. Diversity in star forming environments can have long lasting effects in disk evolution and in the resulting planetary populations. We perform simulations to explore the evolution of circumstellar disks in young star clusters. We include viscous evolution, as well as the impact of dynamical encounters and external photoevaporation. We find that photoevaporation is an important process in destroying circumstellar disks: in regions of stellar density $ρ\sim 100 \mathrm{\ M}_\odot \mathrm{\ pc}^{-3}\mathrm{\ }$ around 80% of disks are destroyed before 2 Myr of cluster evolution. Our findings are in agreement with observed disk fractions in young star forming regions and support previous estimations that planet formation must start in timescales < 0.1 - 1 Myr.
△ Less
Submitted 18 October, 2019; v1 submitted 8 July, 2019;
originally announced July 2019.
-
VPLanet: The Virtual Planet Simulator
Authors:
Rory Barnes,
Rodrigo Luger,
Russell Deitrick,
Peter Driscoll,
Thomas R. Quinn,
David P. Fleming,
Hayden Smotherman,
Diego V. McDonald,
Caitlyn Wilhelm,
Rodolfo Garcia,
Patrick Barth,
Benjamin Guyer,
Victoria S. Meadows,
Cecilia M. Bitz,
Pramod Gupta,
Shawn D. Domagal-Goldman,
John Armstrong
Abstract:
We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evo…
▽ More
We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.
△ Less
Submitted 27 August, 2019; v1 submitted 15 May, 2019;
originally announced May 2019.
-
Exo-Milankovitch Cycles II: Climates of G-dwarf Planets in Dynamically Hot Systems
Authors:
Russell Deitrick,
Rory Barnes,
Cecilia Bitz,
David Fleming,
Benjamin Charnay,
Victoria Meadows,
Caitlyn Wilhelm,
John Armstrong,
Thomas R. Quinn
Abstract:
Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball…
▽ More
Using an energy balance model with ice sheets, we examine the climate response of an Earth-like planet orbiting a G dwarf star and experiencing large orbital and obliquity variations. We find that ice caps couple strongly to the orbital forcing, leading to extreme ice ages. In contrast with previous studies, we find that such exo-Milankovitch cycles tend to impair habitability by inducing snowball states within the habitable zone. The large amplitude changes in obliquity and eccentricity cause the ice edge, the lowest latitude extent of the ice caps, to become unstable and grow to the equator. We apply an analytical theory of the ice edge latitude to show that obliquity is the primary driver of the instability. The thermal inertia of the ice sheets and the spectral energy distribution of the G dwarf star increase the sensitivity of the model to triggering runaway glaciation. Finally, we apply a machine learning algorithm to demonstrate how this technique can be used to extend the power of climate models. This work illustrates the importance of orbital evolution for habitability in dynamically rich planetary systems. We emphasize that as potentially habitable planets are discovered around G dwarfs, we need to consider orbital dynamics.
△ Less
Submitted 1 May, 2018;
originally announced May 2018.
-
Exo-Milankovitch Cycles I: Orbits and Rotation States
Authors:
Russell Deitrick,
Rory Barnes,
Thomas R. Quinn,
John Armstrong,
Benjamin Charnay,
Caitlyn Wilhelm
Abstract:
The obliquity of the Earth, which controls our seasons, varies by only ~2.5 degrees over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth's ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities,…
▽ More
The obliquity of the Earth, which controls our seasons, varies by only ~2.5 degrees over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth's ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semi-analytic model that allows us to map broad regions of parameter space. We find that in general, the largest drivers of obliquity variations are secular spin-orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.
△ Less
Submitted 28 December, 2017;
originally announced December 2017.
