Planets residing in circumstellar habitable zones offer us the best opportunities to test hypotheses of life’s potential pervasiveness and complexity. Constraining the precise boundaries of habitability and its observational discriminants is critical to maximizing our chances at remote life detection with future instruments. Conventionally, calculations of the inner edge of the habitable zone (IHZ) have been performed using both 1D radiative-convective and 3D general circulation models. However, these models lack interactive 3D chemistry and do not resolve the mesosphere and lower thermosphere region of the upper atmosphere. Here, we employ a 3D high-top chemistry-climate model (CCM) to simulate the atmospheres of synchronously rotating planets orbiting at the inner edge of habitable zones of K- and M-dwarf stars (between T eff = 2600 and 4000 K). While our IHZ climate predictions are in good agreement with general circulation model studies, we find noteworthy departures in simulated ozone and HOx photochemistry. For instance, climates around inactive stars do not typically enter the classical moist greenhouse regime even with high (gsim10?3 mol mol?1) stratospheric water vapor mixing ratios, which suggests that planets around inactive M-stars may only experience minor water-loss over geologically significant timescales. In addition, we find much thinner ozone layers on potentially habitable moist greenhouse atmospheres, as ozone experiences rapid destruction via reaction with hydrogen oxide radicals. Using our CCM results as inputs, our simulated transmission spectra show that both water vapor and ozone features could be detectable by instruments NIRSpec and MIRI LRS on board the James Webb Space Telescope.
Habitability is a measure of an environment’s potential to support life, and a habitable exoplanet supports liquid water on its surface. However, a planet’s success in maintaining liquid water on its surface is the end result of a complex set of interactions between planetary, stellar, planetary system and even Galactic characteristics and processes, operating over the planet’s lifetime. In this chapter, we describe how we can now determine which exoplanets are most likely to be terrestrial, and the research needed to help define the habitable zone under different assumptions and planetary conditions. We then move beyond the habitable zone concept to explore a new framework that looks at far more characteristics and processes, and provide a comprehensive survey of their impacts on a planet’s ability to acquire and maintain habitability over time. We are now entering an exciting era of terrestiral exoplanet atmospheric characterization, where initial observations to characterize planetary composition and constrain atmospheres is already underway, with more powerful observing capabilities planned for the near and far future. Understanding the processes that affect the habitability of a planet will guide us in discovering habitable, and potentially inhabited, planets.
We use a one-dimensional (1D) cloud-free climate model to estimate habitable zone (HZ) boundaries for terrestrial planets of masses 0.1 ME and 5 ME around circumbinary stars of various spectral type combinations. Specifically, we consider binary systems with host spectral types F-F, F-G, F-K, F-M, G-G, G-K, G-M, K-K, K-M and M-M. Scaling the background N2 atmospheric pressure with the radius of the planet, we find that the inner edge of the HZ moves inwards toward the star for 5 ME compared to 0.1 ME planets for all spectral types. This is because the water-vapor column depth is smaller for larger planets and higher temperatures are needed before water vapor completely dominates the outgoing longwave radiation. The outer edge of the HZ changes little due to competing effects of the albedo and greenhouse effect. While these results are broadly consistent with the trend of single star HZ results for different mass planets, there are significant differences between single star and binary star systems for the inner edge of the HZ. Interesting combinations of stellar pairs from our 1D model results can be used to explore for in-depth climate studies with 3D climate models. We identify a common HZ stellar flux domain for all circumbinary spectral types.
The pathways through which incoming energy is distributed between the surface and atmosphere have been analyzed for the Earth. However, the effect of the spectral energy distribution of a host star on the energy budget of an orbiting planet may be significant given the wavelength-dependent absorption properties of atmospheric CO2, water vapor, surface ice, and snow. We have quantified the flow of energy on aqua planets orbiting M-, G-, and F-dwarf stars, using a 3D Global Climate Model with a static ocean. The atmosphere and surface of an M-dwarf planet receiving an instellation equal to 88% of the modern solar constant at the top of the atmosphere absorb 12% more incoming stellar radiation than those of a G-dwarf planet receiving 100% of the modern solar constant, and 17% more radiation than an F-dwarf planet receiving 108% of the modern solar constant, resulting in climates similar to that of modern-day Earth on all three planets, assuming a 24 hr rotation period and fixed CO2. At 100% instellation, a synchronously rotating M-dwarf planet exhibits smaller flux absorption in the atmosphere and on the surface of the dayside, and a dayside mean surface temperature that is 37 K colder than its rapidly rotating counterpart. Energy budget diagrams are included to illustrate the variations in global energy budgets as a function of host star spectral class, and can contribute to habitability assessments of planets as they are discovered.
The macroturbulent atmospheric circulation of Earth-like planets mediates their equator-to-pole heat transport. For fast-rotating terrestrial planets, baroclinic instabilities in the mid-latitudes lead to turbulent eddies that act to transport heat poleward. In this work, we derive a scaling theory for the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanet atmospheres. This theory is built on the work of Jansen & Ferrari and determines how unstable the atmosphere is to baroclinic instability (the baroclinic “criticality”) through a balance between the baroclinic eddy heat flux and radiative heating/cooling. We compare our scaling theory to General Circulation Model (GCM) simulations and find that the theoretical predictions for equator-to-pole temperature contrast and bulk lapse rate broadly agree with GCM experiments with varying rotation rate and surface pressure throughout the baroclincally unstable regime. Our theoretical results show that baroclinic instabilities are a strong control of heat transport in the atmospheres of Earth-like exoplanets, and our scalings can be used to estimate the equator-to-pole temperature contrast and bulk lapse rate of terrestrial exoplanets. These scalings can be tested by spectroscopic retrievals and full-phase light curves of terrestrial exoplanets with future space telescopes.
This article presents the proportional evolutionary time (PET) hypothesis, which posits that the mean time required for the evolution of complex life is a function of stellar mass. The biological available window is defined as the region of a stellar spectrum between 200 and 1200?nm that generates free energy for life. Over the ?4 Gyr history of Earth, the total energy incident at the top of the atmosphere and within the biological available window is ?1034 J. The hypothesis assumes that the rate of evolution from the origin of life to complex life is proportional to this total energy, which would suggest that planets orbiting other stars should not show signs of complex life if the total energy incident on the planet is below this energy threshold. The PET hypothesis predicts that late K- and M-dwarf stars (M??0.7 ) represent the best targets for the next generation of space telescopes to search for spectroscopic biosignatures indicative of complex life.
We investigate planetesimal accretion via a direct N-body simulation of an annulus at 1 au orbiting a 1 M? star. The planetesimal ring, which initially contains N = 106 bodies is evolved into the oligarchic growth phase. Unlike previous lower resolution studies, we find that the mass distribution of planetesimals develops a bump at intermediate mass after the oligarchs form. This feature marks a boundary between growth modes. The smallest planetesimals are packed tightly enough together to populate mean motion resonances with the oligarchs, which heats the small bodies, enhancing their growth. If we depopulate most of the resonances by decreasing the width of the annulus, this effect becomes weaker. To clearly demonstrate the dynamics driving these growth modes, we also examine the evolution of a planetary embryo embedded in an annulus of collisionless planetesimals. In this case, we find that the resonances push planetesimals away from the embryo, decreasing the surface density of the bodies adjacent to the embryo. This effect only occurs when the annulus is wide enough and the mass resolution of the planetesimals is fine enough to populate the resonances. The bump we observe in the mass distribution resembles the 100 km power-law break seen in the size distribution of asteroid belt objects. Although the bump produced in our simulations occurs at a size larger than 100 km, we show that the bump location is sensitive to the initial planetesimal mass, which implies that this feature is potentially useful for constraining planetesimal formation models.
New research from VPL Team members Jacob Lustig-Yaeger, Andrew Lincowski, and their faculty advisor Dr. Victoria Meadows, recently published in…
Thermal history models, historically used to understand Earth’s geologic history, are being coupled to climate models to map conditions that allow planets to maintain life. However, the lack of structural uncertainty assessment has blurred guidelines for how thermal history models can be used toward this end. Structural uncertainty is intrinsic to the modeling process. Model structure refers to the cause and effect relations that define a model and are assumed to adequately represent a particular real world system. Intrinsic/structural uncertainty is different from input and parameter uncertainties (which are often evaluated for thermal history models). A full uncertainty assessment requires that input/parametric and intrinsic/structural uncertainty be evaluated (one is not a substitute for the other). We quantify the intrinsic uncertainty for several parameterized thermal history models (a subclass of planetary models). We use single perturbation analysis to determine the reactance time of different models. This provides a metric for how long it takes low?amplitude, unmodeled effects to decay or grow. Reactance time is shown to scale inversely with the strength of the dominant model feedback (negative or positive). A perturbed physics analysis is then used to determine uncertainty shadows for model outputs. This provides probability distributions for model predictions. It also tests the structural stability of a model (do model predictions remain qualitatively similar, and within assumed model limits, in the face of intrinsic uncertainty?). Once intrinsic uncertainty is accounted for, model outputs/predictions and comparisons to observational data should be treated in a probabilistic way.
The search for an inhabited planet, beyond our own, is a driver of planetary exploration in our solar system and beyond. Using information from our own planet to inform search strategies allows for a targeted search. It is, however, worth considering some span in the strategy and in a priori expectation. An inhabited, Earth-like planet is one that would be similar to Earth in ways that extend beyond having biota. To facilitate a comparative cost/risk/benefit analysis of different potential search strategies, we use a metric akin to the Earth-similarity index. The metric extends from zero, for an inhabited planet that is like Earth in all other regards (i.e., zero differences), toward end-member values for planets that differ from Earth but maintain life potential. The analysis shows how finding inhabited planets that do not share other Earth characteristics could improve our ability to assess galactic life potential without a large increase in time-commitment costs. Search strategies that acknowledge the possibility of such planets can minimize the potential of exploration losses (e.g., searching for long durations to reach conclusions that are biased). Discovering such planets could additionally provide a test of the Gaia hypothesisa test that has proven difficult when using only Earth as a laboratory. Finally, we discuss how an Earth2.0 narrative that has been presented to the public as a search strategy comes with nostalgia-laden baggage that does not best serve exploration.
The discovery of planets orbiting stars other than the Sun has accelerated over the past decade, and this trend will continue as new space- and ground-based observatories employ next-generation instrumentation to search the skies for habitable worlds. However, many factors and processes can affect planetary habitability and must be understood to accurately determine a planet’s habitability potential. While climate models have long been used to understand and predict climate and weather patterns on the Earth, a growing community of researchers has begun to apply these models to extrasolar planets. This work has provided a better understanding of how orbital, surface, and atmospheric properties affect planetary climate and habitability; how these climatic effects might change for different stellar and planetary environments; and how the habitability and observational signatures of newly discovered planets might be influenced by these climatic factors. This review summarizes the origins and evolution of the burgeoning field of exoplanet climatology, discusses recent work using a hierarchy of computer models to identify those planets most capable of supporting life, and offers a glimpse into future directions of this quickly evolving subfield of exoplanet science.
VPL Scientist Dr. Lucianne Walkowicz (Adler Planetarium) was recently featured on the Ask an Astrobiologist program! Dr. Lucianne Walkowicz is…
Super-Earths planets with sizes between the Earth and Neptune are found in tighter orbits than that of the Earth around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that also formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two clearly distinct formation pathways that are regulated by the pebble reservoir of the disc. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of approximately Mars-mass embryos when the gas disc dissipates. Subsequently, without gas being present, the embryos become unstable due to mutual gravitational interactions and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most five Earth masses. Instead, if the pebble mass flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the five to twenty Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. Therefore, the key difference between the two growth modes is whether embryos grow fast enough to undergo significant migration. The terrestrial growth mode produces small rocky planets on wider orbits like those in the solar system whereas the super-Earth growth mode produces planets in short-period orbits inside 1 AU, with masses larger than the Earth that should be surrounded by a primordial H/He atmosphere, unless subsequently lost by stellar irradiation. The pebble flux which controls the transition between the two growth modes may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.
The discovery of the first interstellar object passing through the Solar System, 1I/2017 U1 (Oumuamua), provoked intense and continuing interest from the scientific community and the general public. The faintness of Oumuamua, together with the limited time window within which observations were possible, constrained the information available on its dynamics and physical state. Here we review our knowledge and find that in all cases, the observations are consistent with a purely natural origin for Oumuamua. We discuss how the observed characteristics of Oumuamua are explained by our extensive knowledge of natural minor bodies in our Solar System and our current knowledge of the evolution of planetary systems. We highlight several areas requiring further investigation.
Terrestrial planets orbiting M dwarfs may soon be observed with the James Webb Space Telescope (JWST) to characterize their atmospheric composition and search for signs of habitability or life. These planets may undergo significant atmospheric and ocean loss due to the superluminous pre-main-sequence phase of their host stars, which may leave behind abiotically generated oxygen, a false positive for the detection of life. Determining if ocean loss has occurred will help assess potential habitability and whether or not any O2 detected is biogenic. In the solar system, differences in isotopic abundances have been used to infer the history of ocean loss and atmospheric escape (e.g., Venus, Mars). We find that isotopologue measurements using transit transmission spectra of terrestrial planets around late-type M dwarfs like TRAPPIST-1 may be possible with JWST, if the escape mechanisms and resulting isotopic fractionation were similar to Venus. We present analyses of post-ocean-loss O2- and CO2-dominated atmospheres containing a range of trace gas abundances. Isotopologue bands are likely detectable throughout the near-infrared (18 ?m), especially 34 ?m, although not in CO2-dominated atmospheres. For Venus-like D/H ratios 100 times that of Earth, TRAPPIST-1b transit signals of up to 79 ppm are possible by observing HDO. Similarly, 18O/16O ratios 100 times that of Earth produce signals at up to 94 ppm. Detection at signal-to-noise ratio = 5 may be attained on these bands with as few as four to 11 transits, with optimal use of JWST’s NIRSpec Prism. Consequently, H2O and CO2 isotopologues could be considered as indicators of past ocean loss and atmospheric escape for JWST observations of terrestrial planets around M dwarfs.
The habitable zone (HZ) is commonly defined as the range of distances from a host star within which liquid water, a key requirement for life, may exist on a planet’s surface. Substantially more CO2 than present in Earth’s modern atmosphere is required to maintain clement temperatures for most of the HZ, with several bars required at the outer edge. However, most complex aerobic life on Earth is limited by CO2 concentrations of just fractions of a bar. At the same time, most exoplanets in the traditional HZ reside in proximity to M dwarfs, which are more numerous than Sun-like G dwarfs but are predicted to promote greater abundances of gases that can be toxic in the atmospheres of orbiting planets, such as carbon monoxide (CO). Here we show that the HZ for complex aerobic life is likely limited relative to that for microbial life. We use a 1D radiative-convective climate and photochemical models to circumscribe a Habitable Zone for Complex Life (HZCL) based on known toxicity limits for a range of organisms as a proof of concept. We find that for CO2 tolerances of 0.01, 0.1, and 1 bar, the HZCL is only 21%, 32%, and 50% as wide as the conventional HZ for a Sun-like star, and that CO concentrations may limit some complex life throughout the entire HZ of the coolest M dwarfs. These results cast new light on the likely distribution of complex life in the universe and have important ramifications for the search for exoplanet biosignatures and technosignatures.
Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location of water vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.
Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location of water vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.
A rigorous definition of the habitable zone and its dependence on planetary properties is part of the search for habitable exoplanets. In this work, we use the general circulation model ExoCAM to determine how the inner edge of the habitable zone of tidally locked planets orbiting M dwarf stars depends on planetary radius, surface gravity, and surface pressure. We find that the inner edge of the habitable zone for more massive planets occurs at higher stellar irradiation, as found in previous 1D simulations. We also determine the relative effects of varying planetary radius and surface gravity. Increasing the planetary radius leads to a lower planetary albedo and warmer climate, pushing the inner edge of the habitable zone to lower stellar irradiation. This results from a change in circulation regime that leads to the disruption of the thick, reflective cloud deck around the substellar point. Increasing gravity increases the outgoing longwave radiation, which moves the inner edge of the habitable zone to higher stellar irradiation. This is because the column mass of water vapor decreases with increasing gravity, leading to a reduction in the greenhouse effect. The effect of gravity on the outgoing longwave radiation is stronger than the effect of radius on the planetary albedo, so that increasing gravity and radius together causes the inner edge of the habitable zone to move to higher stellar irradiation. Our results show that the inner edge of the habitable zone for more massive terrestrial planets occurs at a larger stellar irradiation.
Recent observations have found a valley in the size distribution of close-in super-Earths that is interpreted as a signpost that close-in super-Earths are mostly rocky in composition. However, new models predict that planetesimals should first form at the water ice line such that close-in planets are expected to have a significant water ice component. Here we investigate the water contents of super-Earths by studying the interplay between pebble accretion, planet migration and disc evolution. Planets compositions are determined by their position relative to different condensation fronts (ice lines) throughout their growth. Migration plays a key role. Assuming that planetesimals start at or exterior to the water ice line (r > rH2O), inward migration causes planets to leave the source region of icy pebbles and therefore to have lower final water contents than in discs with either outward migration or no migration. The water ice line itself moves inward as the disc evolves, and delivers water as it sweeps across planets that formed dry. The relative speed and direction of planet migration and inward drift of the water ice line is thus central in determining planets water contents. If planet formation starts at the water ice line, this implies that hot close-in super-Earths (r < 0.3 AU) with water contents of a few percent are a signpost of inward planet migration during the early gas phase. Hot super-Earths with larger water ice contents on the other hand, experienced outward migration at the water ice line and only migrated inwards after their formation was complete either because they become too massive to be contained in the region of outward migration or in chains of resonant planets. Measuring the water ice content of hot super-Earths may thus constrain their migration history.
Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an intercomparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (2030 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance.
Giant planets migrate though the protoplanetary disc as they grow their solid core and attract their gaseous envelope. Previously, we have studied the growth and migration of an isolated planet in an evolving disc. Here, we generalise such models to include the mutual gravitational interaction between a high number of growing planetary bodies. We have investigated how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 2040 AU when using nominal type-I and type-II migration rates and a pebble flux of approximately 100200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30 AU migrate into the inner system (rP < 1AU). There they form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explored slower migration rates which allow the formation of gas giants from embryos originating from the 510 AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. These giant planets can also form in discs with lower pebbles fluxes (50100 Earth masses per Myr). We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. In our model, giant planet growth requires a sufficiently high pebble flux to enable growth to out-compete migration. An even higher pebble flux produces systems with multiple gas giants. We show that planetary embryos starting interior to 5 AU do not grow into gas giants, even if migration is slow and the pebble flux is large. These embryos instead grow to just a few Earth masses, the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is therefore independent of the pebble flux. Additionally, we show that the long-term evolution of our formed planetary systems can naturally produce systems with inner super-Earths and outer gas giants as well as systems of giant planets on very eccentric orbits.
Oxygen and methane are considered to be the canonical biosignatures of modern Earth, and the simultaneous detection of these gases in a planetary atmosphere is an especially strong biosignature. However, these gases may be challenging to detect together in the planetary atmospheres because photochemical oxygen radicals destroy methane. Previous work has shown that the photochemical lifetime of methane in oxygenated atmospheres is longer around M dwarfs, but M dwarf planet habitability may be hindered by extreme stellar activity and evolution. Here, we use a 1D photochemical-climate model to show that K dwarf stars also offer a longer photochemical lifetime of methane in the presence of oxygen compared to G dwarfs. For example, we show that a planet orbiting a K6V star can support about an order of magnitude more methane in its atmosphere compared to an equivalent planet orbiting a G2V star. In the reflected-light spectra of worlds orbiting K dwarf stars, strong oxygen and methane features could be observed at visible and near-infrared wavelengths. Because K dwarfs are dimmer than G dwarfs, they offer a better planet-star contrast ratio, enhancing the signal-to-noise ratio (S/N) possible in a given observation. For instance, a 50 hr observation of a planet at 7 pc with a 15 m telescope yields S/N = 9.2 near 1 ?m for a planet orbiting a solar-type G2V star, and S/N = 20 for the same planet orbiting a K6V star. In particular, nearby mid-late K dwarfs such as 61 Cyg A/B, Epsilon Indi, Groombridge 1618, and HD 156026 may be excellent targets for future biosignature searches.
At least 30% of main sequence stars host planets with sizes between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk-planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes (?15M?) and super-massive planetary cores potentially able to become gas giant planets (?15M?). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the disk’s inner edge forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to ?95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the resonant configuration. Our simulations match observations if we combine a dominant fraction (?95%) of unstable systems with a sprinkling (?5%) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (N?2) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet-planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice-rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths (abridged)
Knowing the high-energy radiation environment of a star over a planet’s formation and evolutionary period is critical in determining if that planet is potentially habitable and if any biosignatures could be detected, as UV radiation can severely change or destroy a planet’s atmosphere. Current efforts for finding a potentially habitable planet are focused on M stars, yet K stars may offer more habitable conditions due to decreased stellar activity and more distant and wider habitable zones (HZs). While M star activity evolution has been observed photometrically and spectroscopically, there has been no dedicated investigation of K star UV evolution. We present the first comprehensive study of the near-UV, far-UV, and X-ray evolution of K stars. We used members of young moving groups and clusters ranging in age from 10 to 625 Myr combined with field stars and their archived GALEX UV and ROSAT X-ray data to determine how the UV and X-ray radiation evolve. We find that the UV and X-ray flux incident on an HZ planet is 550 times lower than that of HZ planets around early-M stars and 501000 times lower than those around late-M stars, due to both an intrinsic decrease in K dwarf stellar activity occurring earlier than for M dwarfs and the more distant location of the K dwarf HZ.