The orbital eccentricity distribution of planets orbiting M dwarfs.

We investigate the underlying distribution of orbital eccentricities for planets around early-to-mid M dwarf host stars. We employ a sample of 163 planets around early-to mid-M dwarfs across 101 systems detected by NASA's Kepler Mission. We constrain the orbital eccentricity for each planet by leveraging the Kepler lightcurve together with a stellar density prior, constructed using metallicity from spectroscopy, Ks magnitude from 2MASS, and stellar parallax from Gaia. Within a Bayesian hierarchical framework, we extract the underlying eccentricity distribution, assuming alternately Rayleigh, half-Gaussian, and Beta functions for both single- and multi-transit systems. We described the eccentricity distribution for apparently single-transiting planetary systems with a Rayleigh distribution with σ = 0:19+0:04-0:03, and for multitransit systems with σ = 0:03+0:02-0:01. The data suggest the possibility of distinct dynamically warmer and cooler subpopulations within the single-transit distribution: The single-transit data prefer a mixture model composed of two distinct Rayleigh distributions with σ1 = 0:02+0:11-0:00 and σ2 = 0:24+0:20-0:03 over a single Rayleigh distribution, with 7:1 odds. We contextualize our findings within a planet formation framework, by comparing them to analogous results in the literature for planets orbiting FGK stars. By combining our derived eccentricity distribution with other M dwarf demographic constraints, we estimate the underlying eccentricity distribution for the population of early- to mid-M dwarf planets in the local neighborhood. [ABSTRACT FROM AUTHOR]