Quantifying the evolution of stellar extreme ultraviolet (EUV, 100 -- 1000 $\overset{\circ}{A}$) emission is critical for assessing the evolution of planetary atmospheres and the habitability of M dwarf systems. Previous studies from the HAbitable Zones and M dwarf Activity across Time (HAZMAT) program showed the far- and near-UV (FUV, NUV) emission from M stars at various stages of a stellar lifetime through photometric measurements from the Galaxy Evolution Explorer (GALEX). The results revealed increased levels of short-wavelength emission that remain elevated for hundreds of millions of years. The trend for EUV flux as a function of age could not be determined empirically because absorption by the interstellar medium prevents access to the EUV wavelengths for the vast majority of stars. In this paper, we model the evolution of EUV flux from early M stars to address this observational gap. We present synthetic spectra spanning EUV to infrared wavelengths of 0.4 $\pm$ 0.05 M$_{\odot}$ stars at five distinct ages between 10 and 5000 Myr, computed with the PHOENIX atmosphere code and guided by the GALEX photometry. We model a range of EUV fluxes spanning two orders of magnitude, consistent with the observed spread in X-ray, FUV, and NUV flux at each epoch. Our results show that the stellar EUV emission from young M stars is 100 times stronger than field age M stars, and decreases as t$^{-1}$ after remaining constant for a few hundred million years. This decline stems from changes in the chromospheric temperature structure, which steadily shifts outward with time. Our models reconstruct the full spectrally and temporally resolved history of an M star's UV radiation, including the unobservable EUV radiation, which drives planetary atmospheric escape, directly impacting a planet's potential for habitability.

We search for evidence of the cause of the exoplanet radius gap, i.e. the dearth of planets with radii near $1.8\ R_\oplus$. If the cause was photoevaporation, the radius gap should trend with proxies for the early-life high-energy emission of planet-hosting stars. If, alternatively, the cause was core-powered mass loss, no such trends should exist. Critically, spurious trends between the radius gap and stellar properties arise from an underlying correlation with instellation. After accounting for this underlying correlation, we find no trends remain between the radius gap and stellar mass or present-day stellar activity as measured by near-UV emission. We dismiss the nondetection of a radius gap trend with near-UV emission because present-day near-UV emission is unlikely to trace early-life high-energy emission, but we provide a catalog of GALEX near-UV and far-UV emission measurements for general use. We interpret the nondetection of a radius gap trend with stellar mass by simulating photoevaporation with mass-dependent evolution of stellar high-energy emission. The simulation produces an undetectable trend between the radius gap and stellar mass under realistic sources of error. We conclude that no evidence, from this analysis or others in the literature, currently exists that clearly favors either photoevaporation or core powered mass loss as the primary cause of the exoplanet radius gap. However, repeating this analysis once the body of well-characterized $< 4\ R_\oplus$ planets has roughly doubled could confirm or rule out photoevaporation.

Correct estimates of stellar extreme ultraviolet (EUV; 100 - 1170 \AA) flux are important for studying the photochemistry and stability of exoplanet atmospheres, as EUV radiation ionizes hydrogen and contributes to the heating, expansion, and potential escape of a planet's upper atmosphere. Contamination from interstellar hydrogen makes observing EUV emission from M stars particularly difficult, and impossible past 100 pc, and necessitates other means to predict the flux in this wavelength regime. We present EUV -- infrared (100 \AA - 5.5 $\mu$m) synthetic spectra computed with the PHOENIX atmospheric code of three early M dwarf planet hosts: GJ 832 (M1.5 V), GJ 176 (M2.5 V), and GJ 436 (M3.5 V). These one-dimensional semiempirical nonlocal thermodynamic equilibrium models include simple temperature prescriptions for the stellar chromosphere and transition region, from where ultraviolet (UV; 100 - 3008 \AA) fluxes originate. We guide our models with Hubble Space Telescope far- and near-UV spectra and discuss the ability to constrain these models using Galaxy Evolution Explorer UV photometry. Our models closely reproduce the observations and predict the unobservable EUV spectrum at a wavelength resolution of < 0.1 \AA. The temperature profiles that best reproduce the observations for all three stars are described by nearly the same set of parameters, suggesting that early M type stars may have similar thermal structures in their upper atmospheres. With an impending UV observation gap and the scarcity of observed EUV spectra for stars less luminous and more distant than the Sun, upper-atmosphere models such as these are important for providing realistic spectra across short wavelengths and for advancing our understanding of the effects of radiation on planets orbiting M stars.

The high energy radiation environment around M dwarf stars strongly impacts the characteristics of close-in exoplanet atmospheres, but these wavelengths are difficult to observe due to geocoronal and interstellar contamination. On account of these observational restrictions, a stellar atmosphere model may be used to compute the stellar extreme ultraviolet (EUV; 100 - 912 \AA) spectrum. We present a case study of the ultra-cool M8 dwarf star, TRAPPIST-1, which hosts seven transiting short-period terrestrial sized planets whose atmospheres will be probed by the James Webb Space Telescope. We construct semi-empirical non-LTE model spectra of TRAPPIST-1 that span EUV to infrared wavelengths (100 \AA\ - 2.5 $\mu$m) using the atmosphere code PHOENIX. These upper-atmosphere models contain prescriptions for the chromosphere and transition region and include newly added partial frequency redistribution capabilities. In the absence of broadband UV spectral observations, we constrain our models using HST Ly$\alpha$ observations from TRAPPIST-1 and GALEX FUV and NUV photometric detections from a set of old M8 stars ($>$1 Gyr). We find that calibrating the models using both data sets separately yield similar FUV and NUV fluxes, and EUV fluxes that range from (1.32 - 17.4) $\times$ 10$^{-14}$ ergs s$^{-1}$ cm$^{-2}$. The results from these models demonstrate that the EUV emission is very sensitive to the temperature structure in the transition region. Our lower activity models predict EUV fluxes similar to previously published estimates derived from semi-empirical scaling relationships, while the highest activity model predicts EUV fluxes a factor of ten higher. Results from this study support the idea that the TRAPPIST-1 habitable zone planets likely do not have much liquid water on their surfaces due to the elevated levels of high energy radiation emitted by the host star.

M stars are powerful emitters of far-ultraviolet light. Over long timescales, a significant, possibly dominant, fraction of this emission is produced by stellar flares. Characterizing this emission is critical to understanding the atmospheres of the stars producing it and the atmospheric evolution of the orbiting planets subjected to it. Ultraviolet emission is known to be elevated for several hundred million years after M stars form. Whether the same is true of ultraviolet flare activity is a key concern for the evolution of exoplanet atmospheres. Hubble Space Telescope (HST) observations by the HAZMAT program (HAbitable Zones and M dwarf Activity across Time) detected 18 flares on young (40 Myr) early M stars in the Tucana-Horologium association over 10 h of observations, ten having energy $>10^{30}$ erg. These imply flares on young M stars are 100--1000$\times$ more energetic than those occurring at the same rate on "inactive," field-age M dwarfs. However, when energies are are normalized by quiescent emission there is no statistical difference between the young and field-age samples. The most energetic flare observed, dubbed the "Hazflare," emitted an energy of $10^{32.1}$ erg in the FUV, 30$\times$ more energetic than any stellar flare previously observed in the FUV with HST's COS or STIS spectrographs. It was accompanied by $15,500\pm400$ K blackbody emission bright enough to designate it a superflare ($E>10^{33}$ erg), with an estimated bolometric energy of $10^{33.6_{-0.2}^{+0.1}}$ erg. This blackbody emitted 18$^{+2}_{-1}$% of its flux in the FUV (912--1700 \AA) where molecules are generally most sensitive to photolysis. Such hot superflares in young, early M stars could play an important role in the evolution of nascent planetary atmospheres.

The NASA K2 mission uses photometry to find planets transiting stars of various types. M dwarfs are of high interest since they host more short period planets than any other type of main sequence stars and transiting planets around M dwarfs have deeper transits compared to other main sequence stars. In this paper, we present stellar parameters from K and M dwarfs hosting transiting planet candidates discovered by our team. Using the SOFI spectrograph on the European Southern Observatory's New Technology Telescope, we obtained R ~ 1000 J-, H-, and K-band (0.95 - 2.52 microns) spectra of 34 late-type K2 planet and candidate planet host systems and 12 bright K4-M5 dwarfs with interferometrically measured radii and effective temperatures. Out of our 34 late-type K2 targets, we identify 27 of these stars as M dwarfs. We measure equivalent widths of spectral features, derive calibration relations using stars with interferometric measurements, and estimate stellar radii, effective temperatures, masses, and luminosities for the K2 planet hosts. Our calibrations provide radii and temperatures with median uncertainties of 0.059 R_sun (16.09%) and 160 K (4.33%), respectively. We then reassess the radii and equilibrium temperatures of known and candidate planets based on our spectroscopically derived stellar parameters. Since a planet's radius and equilibrium temperature depends on the parameters of its host star, our study provides more precise planetary parameters for planets and candidates orbiting late-type stars observed with K2. We find a median planet radius and an equilibrium temperature of approximately 3R_earth and 500 K, with several systems (K2-18b and K2-72e) receiving near-Earth like levels of incident irradiation.

Some close-in gaseous exoplanets are nearly in Roche-lobe contact, and previous studies show tidal decay can drive hot Jupiters into contact during the main sequence of their host stars. Improving upon a previous model, we present a revised model for mass transfer in a semi-detached binary system that incorporates an extended atmosphere around the donor and allows for an arbitrary mass ratio. We apply this new formalism to hypothetical, confirmed, and candidate planetary systems to estimate mass loss rates and compare with models of evaporative mass loss. Overflow may be significant for hot Neptunes out to periods of $\sim$ 2 days, while for hot Jupiters, it may only be important inward of 0.5 days. We find that CoRoT-24 b may be losing mass at a rate of more than an Earth mass in a Gyr. The hot Jupiter WASP-12 b may lose an Earth mass in a Myr, while the putative planet orbiting a T-Tauri star PTFO8-8695 might shed its atmosphere in a few Myrs. We point out that the orbital expansion that can accompany mass transfer may be less effective than previously considered because the gas accreted by the host star removes some of the system's angular momentum from the orbit, but simple scaling arguments suggest that the Roche-lobe overflow might remain stable. Consequently, the recently discovered small planets in ultra-short-periods ($<$ 1 day) may not be the remnants of hot Jupiters/Neptunes. The new model presented here has been incorporated into Modules for Experiments in Stellar Astrophysics (MESA).

Many gaseous exoplanets in short-period orbits are on the verge or are in the process of Roche-lobe overflow (RLO). Moreover, orbital stability analysis shows tides can drive many hot Jupiters to spiral inevitably toward their host stars. Thus, the coupled processes of orbital evolution and RLO likely shape the observed distribution of close-in exoplanets and may even be responsible for producing some of the short-period rocky planets. However, the exact outcome for an overflowing planet depends on its internal response to mass loss, and the accompanying orbital evolution can act to enhance or inhibit RLO. In this study, we apply the fully-featured and robust Modules for Experiments in Stellar Astrophysics (MESA) suite to model RLO of short-period gaseous planets. We show that, although the detailed evolution may depend on several properties of the planetary system, it is largely determined by the core mass of the overflowing gas giant. In particular, we find that the orbital expansion that accompanies RLO often stops and reverses at a specific maximum period that depends on the core mass. We suggest that RLO may often strand the remnant of a gas giant near this orbital period, which provides an observational prediction that can corroborate the hypothesis that short-period gas giants undergo RLO. We conduct a preliminary comparison of this prediction to the observed population of small, short-period planets and find some planets in orbits that may be consistent with this picture. To the extent that we can establish some short-period planets are indeed the remnants of gas giants, that population can elucidate the properties of gas giant cores, the properties of which remain largely unconstrained.

A star's UV emission can greatly affect the atmospheric chemistry and physical properties of closely orbiting planets with the potential for severe mass loss. In particular, the Lyman-alpha emission line at 1216 Angstroms, which dominates the far-ultraviolet spectrum, is a major source of photodissociation of important atmospheric molecules such as water and methane. The intrinsic flux of Lyman-alpha, however, cannot be directly measured due to the absorption of neutral hydrogen in the interstellar medium and contamination by geocoronal emission. To date, reconstruction of the intrinsic Lyman-alpha line based on Hubble Space Telescope spectra has been accomplished for 46 FGKM nearby stars, 28 of which have also been observed by the Galaxy Evolution Explorer (GALEX). Our investigation provides a correlation between published intrinsic Lyman-alpha and GALEX far- and near-ultraviolet chromospheric fluxes for K and M stars. The negative correlations between the ratio of the Lyman-alpha to the GALEX fluxes reveal how the relative strength of Lyman-alpha compared to the broadband fluxes weakens as the FUV and NUV excess flux increase. We also correlate GALEX fluxes with the strong near-ultraviolet Mg II h+k spectral emission lines formed at lower chromospheric temperatures than Lyman-alpha. The reported correlations provide estimates of intrinsic Lyman-alpha and Mg II fluxes for the thousands of K and M stars in the archived GALEX all-sky surveys. These will constrain new stellar upper-atmosphere models for cool stars and provide realistic inputs to models describing exoplanetary photochemistry and atmospheric evolution in the absence of ultraviolet spectroscopy.