Flare ribbons are bright manifestations of flare energy dissipation in the lower solar atmosphere. For the first time, we report on high-resolution imaging spectroscopy observations of flare ribbons situated off-limb in the H$\beta$ and Ca II 8542 {\AA} lines and make a detailed comparison with radiative hydrodynamic simulations. Observations of the X8.2-class solar flare SOL2017-09-10T16:06 UT obtained with the Swedish Solar Telescope reveal bright horizontal emission layers in H$\beta$ line wing images located near the footpoints of the flare loops. The apparent separation between the ribbon observed in the H$\beta$ wing and the nominal photospheric limb is about 300 - 500 km. The Ca II 8542 {\AA} line wing images show much fainter ribbon emissions located right on the edge of the limb, without clear separation from the limb. RADYN models are used to investigate synthetic spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. The simulations show that, towards the limb, where the line of sight is substantially oblique with respect to the vertical direction, the flaring atmosphere model reproduces the high contrast of the off-limb H$\beta$ ribbons and their significant elevation above the photosphere. The ribbons in the Ca II 8542 {\AA} line wing images are located deeper in the lower solar atmosphere with a lower contrast. A comparison of the height deposition of electron beam energy and the intensity contribution function shows that the H$\beta$ line wing intensities can be an useful tracer of flare energy deposition in the lower solar atmosphere

We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 {\AA}), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 {\AA}), NICER (Neutron Star Interior Composition Explorer; 0.2--12.0 keV), and collaborations of OISTER (Optical and Infrared Synergetic Telescopes for Education and Research) program. Twelve flares are detected in total which include ten H$\alpha$, four X-ray, and four optical-continuum flares; one of them is a superflare with the total energy of $\sim$ 2.0$\times$10$^{33}$ erg. We found that (1) during the superflare, the H$\alpha$ emission line full width at 1/8 maximum dramatically increases to 14 {\AA} from 8 {\AA} in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$\alpha$/X-ray flares are not accompanied with white-light emissions, and (3) the non-flaring emissions show clear rotational modulations in X-ray and H$\alpha$ intensity in the same phase. To understand these observational features, one-dimensional hydrodynamic flare simulations are performed by using the RADYN code. As a result of simulations, we found the simulated H$\alpha$ line profiles with hard and high-energy non-thermal electron beams are consistent with that of the initial phase line profiles of the superflares, while those with more soft- and/or weak-energy beam are consistent with those in decay phases, indicating the changes in the energy fluxes injected to the lower atmosphere. Also, we found that the relation between optical continuum and H$\alpha$ intensity is nonlinear, which can be one cause of the non-white-light flares. The flare energy budget exhibits diversity in the observations and models, and more observations of stellar flares are necessary for constraining the occurrence of various emission line phenomena in stellar flares.

We study the evolution of chromospheric line and continuum emission during the impulsive phase of the X-class SOL2014-09-10T17:45 solar flare. We extend previous analyses of this flare to multiple chromospheric lines of Fe I, Fe II, Mg II, C I, and Si II, observed with IRIS, combined with radiative-hydrodynamical (RHD) modeling. For multiple flaring kernels, the lines all show a rapidly evolving double-component structure: an enhanced, emission component at rest, and a broad, highly red-shifted component of comparable intensity. The red-shifted components migrate from 25-50 km s$^{-1}$ towards the rest wavelength within $\sim$30 seconds. Using Fermi hard X-ray observations, we derive the parameters of an accelerated electron beam impacting the dense chromosphere, using them to drive a RHD simulation with the RADYN code. As in Kowalski et al. 2017a, our simulations show that the most energetic electrons penetrate into the deep chromosphere, heating it to T$\sim$10,000 K, while the bulk of the electrons dissipate their energy higher, driving an explosive evaporation, and its counterpart condensation -- a very dense (n$_e \sim 2 \times 10^{14}$ cm$^{-3}$), thin layer (30--40 km thickness), heated to 8--12,000 K, moving towards the stationary chromosphere at up to 50 km s$^{-1}$. The synthetic Fe II 2814.45A profiles closely resemble the observational data, including a continuum enhancement, and both a stationary and a highly red-shifted component, rapidly moving towards the rest wavelength. Importantly, the absolute continuum intensity, ratio of component intensities, relative time of appearance, and red-shift amplitude, are sensitive to the model input parameters, showing great potential as diagnostics.

Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like stars. In this study, we investigate the possibility of using coronal dimming as a proxy to diagnose stellar CMEs. By simulating a realistic solar CME event and corresponding coronal dimming using a global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we first demonstrate the capability of the model to reproduce solar observations. We then extend the model for simulating stellar CMEs by modifying the input magnetic flux density as well as the initial magnetic energy of the CME flux rope. Our result suggests that with improved instrument sensitivity, it is possible to detect the coronal dimming signals induced by the stellar CMEs.

We present new high fidelity optical coronagraphic imagery of the inner $\sim$50 au of AU Mic's edge-on debris disk using the BAR5 occulter of the Hubble Space Telescope Imaging Spectrograph (HST/STIS) obtained on 26-27 July 2018. This new imagery reveals that "feature A", residing at a projected stellocentric separation of 14.2 au on SE-side of the disk, exhibits an apparent "loop-like" morphology at the time of our observations. The loop has a projected width of 1.5 au and rises 2.3 au above the disk midplane. We also explored TESS photometric observations of AU Mic that are consistent with evidence of two starspot complexes in the system. The likely co-alignment of the stellar and disk rotational axes breaks degeneracies in detailed spot modeling, indicating that AU Mic's projected magnetic field axis is offset from its rotational axis. We speculate that small grains in AU Mic's disk could be sculpted by a time-dependent wind that is influenced by this offset magnetic field axis, analogous to co-rotating Solar interaction regions that sculpt and influence the inner and outer regions of our own Heliosphere. Alternatively, if the observed spot modulation is indicative of a significant mis-alignment of the stellar and disk rotational axes, we suggest the disk could still be sculpted by the differential equatorial versus polar wind that it sees with every stellar rotation.

We explore methods for the identification of stellar flare events in irregularly sampled data of ground-based time domain surveys. In particular, we describe a new technique for identifying flaring stars, which we have implemented in a publicly available Python module called "PyVAN". The approach uses the Differential Evolution algorithm to optimize parameters of empirically derived light-curve templates for different types of stars to fit a candidate light-curve. The difference of the likelihoods that these best-fit templates produced the observed data is then used to delineate targets that are well explained by a flare template but simultaneously poorly explained by templates of common contaminants. By testing on light-curves of known identity and morphology, we show that our technique is capable of recovering flaring status in $69\%$ of all light-curves containing a flare event above thresholds drawn to include $\lt1\%$ of any contaminant population. By applying to Palomar Transient Factory data, we show consistency with prior samples of flaring stars, and identify a small selection of candidate flaring G-type stars for possible follow-up.

The GOES X1 flare SOL2014-10-25T17:08:00 was a three-ribbon solar flare observed with IRIS in the near and far ultraviolet. One of the flare ribbons crossed a sunspot umbra, producing a dramatic, $\sim1000$% increase in the near-ultraviolet (NUV) continuum radiation. We comprehensively analyze the ultraviolet spectral data of the umbral flare brightenings, which provide new challenges for radiative-hydrodynamic modeling of the chromospheric velocity field and the white-light continuum radiation. The emission line profiles in the umbral flare brightenings exhibit redshifts and profile asymmetries, but these are significantly smaller than in another, well-studied X-class solar flare. We present a ratio of the NUV continuum intensity to the Fe II 2814.45 Ang intensity. This continuum-to-line ratio is a new spectral diagnostic of significant heating at high column mass (log $m/$[g cm$^{-2}] >-2$) during solar flares because the continuum and emission line radiation originate from relatively similar temperatures but moderately different optical depths. The full spectral readout of these IRIS data also allow for a comprehensive survey of the flaring NUV landscape: in addition to many lines of Fe II and Cr II, we identify a new solar flare emission line, He I $\lambda2829.91$ (as previously identified in laboratory and early-type stellar spectra). The Fermi/GBM hard X-ray data provide inputs to radiative-hydrodynamic models (which will be presented in Paper II) in order to better understand the large continuum-to-line ratios, the origin of the white-light continuum radiation, and the role of electron beam heating in the low atmosphere.

Observations from the \textit{Interface Region Imaging Spectrograph} (\textsl{IRIS}) often reveal significantly broadened and non-reversed profiles of the Mg II h, k and triplet lines at flare ribbons. To understand the formation of these optically thick Mg II lines, we perform plane parallel radiative hydrodynamics modeling with the RADYN code, and then recalculate the Mg II line profiles from RADYN atmosphere snapshots using the radiative transfer code RH. We find that the current RH code significantly underestimates the Mg II h \& k Stark widths. By implementing semi-classical perturbation approximation results of quadratic Stark broadening from the STARK-B database in the RH code, the Stark broadenings are found to be one order of magnitude larger than those calculated from the current RH code. However, the improved Stark widths are still too small, and another factor of 30 has to be multiplied to reproduce the significantly broadened lines and adjacent continuum seen in observations. Non-thermal electrons, magnetic fields, three-dimensional effects or electron density effect may account for this factor. Without modifying the RADYN atmosphere, we have also reproduced non-reversed Mg II h \& k profiles, which appear when the electron beam energy flux is decreasing. These profiles are formed at an electron density of $\sim 8\times10^{14}\ \mathrm{cm}^{-3}$ and a temperature of $\sim1.4\times10^4$ K, where the source function slightly deviates from the Planck function. Our investigation also demonstrates that at flare ribbons the triplet lines are formed in the upper chromosphere, close to the formation heights of the h \& k lines.

We describe how the accurate characterization of exoplanetary atmospheres in the ELT and JWST era will inevitably require taking into consideration the stellar inhomogeneities caused by convection and magnetic fields. The existing evidence that demonstrates the mixture of stellar and planetary signatures in observed transiting spectra is presented. Finally, we discuss how to disentangle these two components through a multipronged approach that includes new solar reference spectra, improved MHD modeling, and synergistic collaborations between the communities involved, from solar to stellar and exoplanet astronomers.

The planetary effective surface temperature alone is insufficient to characterize exoplanet atmospheres and their stability or evolution. Considering the star-planet system as a whole is necessary, and a critical component of the system is the photoionizing stellar extreme ultraviolet emission (EUV; 100-912 {\AA}). EUV photons drive atmospheric mass loss through thermal and nonthermal processes, and an accurate accounting of the EUV energy deposition in a planet's energy budget is essential, especially for terrestrial habitable zone planets and close-in gaseous planets. Direct EUV observations of exoplanet host stars would require a new, dedicated observatory. Archival observations from the $\textit{EUVE}$ satellite, models, and theory alone are insufficient to accurately characterize EUV spectra of the majority of exoplanet host stars, especially for low-mass stars.

The EUV (100-912 {\AA}) is a spectral region notoriously difficult to observe due to attenuation by neutral hydrogen gas in the interstellar medium. Despite this, hundreds to thousands of nearby stars of different spectral types and magnetic activity levels are accessible in the EUV range. The EUV probes interesting and complicated regions in the stellar atmosphere like the lower corona and transition region that are inaccessible from other spectral regions. In this white paper we describe how direct EUV observations, which require a dedicated grazing-incidence observatory, cannot yet be accurately substituted with models and theory. Exploring EUV emission from cool dwarf stars in the time domain can make a major contribution to understanding stellar outer atmospheres and magnetism, and offers the clearest path toward detecting coronal mass ejections on stars other than the Sun.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

As part of the Mega MUSCLES Hubble Space Telescope (HST) Treasury program, we obtained time-series ultraviolet spectroscopy of the M2.5V star, GJ~674. During the FUV monitoring observations, the target exhibited several small flares and one large flare (E_FUV = 10^{30.75} ergs) that persisted over the entirety of a HST orbit and had an equivalent duration >30,000 sec, comparable to the highest relative amplitude event previously recorded in the FUV. The flare spectrum exhibited enhanced line emission from chromospheric, transition region, and coronal transitions and a blue FUV continuum with an unprecedented color temperature of T_c ~ 40,000+/-10,000 K. In this paper, we compare the flare FUV continuum emission with parameterizations of radiative hydrodynamic model atmospheres of M star flares. We find that the observed flare continuum can be reproduced using flare models but only with the ad hoc addition of hot, dense emitting component. This observation demonstrates that flares with hot FUV continuum temperatures and significant EUV/FUV energy deposition will continue to be of importance to exoplanet atmospheric chemistry and heating even as the host M dwarfs age beyond their most active evolutionary phases.

We present NUV flare spectra from the Hubble Space Telescope/Cosmic Origins Spectrograph during two moderate-amplitude U-band flares on the dM4e star GJ 1243. These spectra are some of the first accurately flux-calibrated, NUV flare spectra obtained over the impulsive phase in M dwarf flares. We observed these flares with a fleet of nine ground-based telescopes simultaneously, which provided broadband photometry and low-resolution spectra at the Balmer jump. A broadband continuum increase occurred with a signal-to-noise > 20 in the HST spectra, while numerous Fe II lines and the Mg II lines also increased but with smaller flux enhancements compared to the continuum radiation. These two events produced the most prominent Balmer line radiation and the largest Balmer jumps that have been observed to date in dMe flare spectra. A T=9000 K blackbody under-estimates the NUV continuum flare flux by a factor of two and is a poor approximation to the white-light in these types of flare events. Instead, our data suggest that the peak of the specific continuum flux density is constrained to U-band wavelengths near the Balmer series limit. A radiative-hydrodynamic simulation of a very high energy deposition rate averaged over times of impulsive heating and cooling better explains the lam>2500 Angstrom flare continuum properties. These two events sample only one end of the empirical color-color distribution for dMe flares, and more time-resolved flare spectra in the NUV, U-band, and optical from lam=2000-4200 Angstrom are needed during more impulsive and/or more energetic flares.

The near-ultraviolet and optical (white-light) continuum radiation in M dwarf flares exhibits a range of observed characteristics, suggesting that the amount of heating at large optical depth varies among impulsive-type and gradual-type flares. Specific flux ratios from high-time cadence spectra and narrowband continuum photometry have also shown that these characteristics vary from the peak to the gradual decay phases of flares. In these proceedings, we present the highest-time cadence ($\sim1$~s), highest signal-to-noise ($\sim100$) constraints on the optical color temperature evolution during the rise phase of a large, impulsive-type dMe flare event. The flare exhibits compelling evidence of a hot, color temperature ($T \sim 10,000$ K), but the Balmer jump ratios show that the flare cannot be explained by isothermal slabs or blackbody surfaces at any time in the evolution. The new data analysis establish these properties as critical challenges for any flare model, and we discuss 1D radiative-hydrodynamic modeling that will be compared to the evolution of the flare colors in this intriguing event.

The hydrogen Lyman lines ( 91.2 nm < lambda < 121.6 nm) are significant contributors to the radiative losses of the solar chromosphere, and are enhanced during flares. We have shown previously that the Lyman lines observed by the Extreme Ultraviolet Variability instrument onboard the Solar Dynamics Observatory exhibit Doppler motions equivalent to speeds on the order of 30 km/s. But contrary to expectation, no dominant flow direction was observed, with both redshifts and blueshifts present. To understand the formation of the Lyman lines, particularly their Doppler motions, we have used the radiative hydrodynamic code, RADYN, and the radiative transfer code, RH, to simulate the evolution of the flaring chromosphere and the response of the Lyman lines during solar flares. We find that upflows in the simulated atmospheres lead to blueshifts in the line cores, which exhibit central reversals. We then model the effects of the instrument on the profiles using the EVE instrument's properties. What may be interpreted as downflows (redshifted emission) in the lines after they have been convolved with the instrumental line profile may not necessarily correspond to actual downflows. Dynamic features in the atmosphere can introduce complex features in the line profiles which will not be detected by instruments with the spectral resolution of EVE, but which leave more of a signature at the resolution of the Spectral Investigation of the Coronal Environment (SPICE) instrument on Solar Orbiter.

Radio observations are currently the only way to explore accelerated particles in cool stellar environments. We describe how a next generation VLA can contribute to the understanding of the stellar contribution to exo-space weather. This area holds both academic and popular interest, and is expected to grow in the next several decades.

We present new analyses of ALMA 12-m and ACA observations at 233 GHz (1.3 mm) of the Proxima Centauri system with sensitivities of 9.5 and 47 $\mu$Jy beam$^{-1}$, respectively, taken from 2017 January 21 through 2017 April 25. These analyses reveal that the star underwent a significant flaring event during one of the ACA observations on 2017 March 24. The complete event lasted for approximately 1 minute and reached a peak flux density of $100\pm4$ mJy, nearly a factor of $1000\times$ brighter than the star's quiescent emission. At the flare peak, the continuum emission is characterized by a steeply falling spectral index with frequency, $F_\nu \propto \nu^\alpha$ with $\alpha = -1.77\pm0.45$, and a lower limit on the fractional linear polarization of $|Q/I| = 0.19\pm0.02$. Since the ACA observations do not show any quiescent excess emission, we conclude that there is no need to invoke the presence of a dust belt at $1-4$ AU. We also posit that the slight excess flux density of $101\pm9$ $\mu$Jy observed in the 12-m observations compared to the photospheric flux density of $74\pm4$ $\mu$Jy extrapolated from infrared wavelengths may be due to coronal heating from continual smaller flares, as is seen for AU Mic, another nearby, well-studied, M dwarf flare star. If this is true, then the need for warm dust at $\sim0.4$ AU is also removed.

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T~10^4 K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for a time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and the emergent continuum flux spectra of radiative-hydrodynamic model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m_ref) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014-Mar-29 X1 solar flare and that variation in m_ref from 0.001 to 0.02 g/cm2 reproduces most of the observed range of the optical continuum flux ratios at the peaks of M dwarf flares.

We report the first direct detection of a strong, 5 kG magnetic field on the surface of an active brown dwarf. LSR J1835+3259 is an M8.5 dwarf exhibiting transient radio and optical emission bursts modulated by fast rotation. We have detected the surface magnetic field as circularly polarized signatures in the 819 nm sodium lines when an active emission region faced the Earth. Modeling Stokes profiles of these lines reveals the effective temperature of 2800 K and log gravity acceleration of 4.5. These parameters place LSR J1835+3259 on evolutionary tracks as a young brown dwarf with the mass of 55$\pm$4 M$_{\rm J}$ and age of 22$\pm$4 Myr. Its magnetic field is at least 5.1 kG and covers at least 11% of the visible hemisphere. The active region topology recovered using line profile inversions comprises hot plasma loops with a vertical stratification of optical and radio emission sources. These loops rotate with the dwarf in and out of view causing periodic emission bursts. The magnetic field is detected at the base of the loops. This is the first time that we can quantitatively associate brown dwarf non-thermal bursts with a strong, 5 kG surface magnetic field and solve the puzzle of their driving mechanism. This is also the coolest known dwarf with such a strong surface magnetic field. The young age of LSR J1835+3259 implies that it may still maintain a disk, which may facilitate bursts via magnetospheric accretion, like in higher-mass T Tau-type stars. Our results pave a path toward magnetic studies of brown dwarfs and hot Jupiters.

Characterizing the UV spectral energy distribution (SED) of an exoplanet host star is critically important for assessing its planet's potential habitability, particularly for M dwarfs as they are prime targets for current and near-term exoplanet characterization efforts and atmospheric models predict that their UV radiation can produce photochemistry on habitable zone planets different than on Earth. To derive ground-based proxies for UV emission for use when Hubble Space Telescope observations are unavailable, we have assembled a sample of fifteen early-to-mid M dwarfs observed by Hubble, and compared their non-simultaneous UV and optical spectra. We find that the equivalent width of the chromospheric Ca II K line at 3933 Angstroms, when corrected for spectral type, can be used to estimate the stellar surface flux in ultraviolet emission lines, including H I Lyman alpha. In addition, we address another potential driver of habitability: energetic particle fluxes associated with flares. We present a new technique for estimating soft X-ray and >10 MeV proton flux during far-UV emission line flares (Si IV and He II) by assuming solar-like energy partitions. We analyze several flares from the M4 dwarf GJ 876 observed with Hubble and Chandra as part of the MUSCLES Treasury Survey and find that habitable zone planets orbiting GJ 876 are impacted by large Carrington-like flares with peak soft X-ray fluxes >1e-3 W m-2 and possible proton fluxes ~100-1000 pfu, approximately four orders of magnitude more frequently than modern-day Earth.

The broadening of the hydrogen lines during flares is thought to result from increased charge (electron, proton) density in the flare chromosphere. However, disagreements between theory and modeling prescriptions have precluded an accurate diagnostic of the degree of ionization and compression resulting from flare heating in the chromosphere. To resolve this issue, we have incorporated the unified theory of electric pressure broadening of the hydrogen lines into the non-LTE radiative transfer code RH. This broadening prescription produces a much more realistic spectrum of the quiescent, A0 star Vega compared to the analytic approximations used as a damping parameter in the Voigt profiles. We test recent radiative-hydrodynamic (RHD) simulations of the atmospheric response to high nonthermal electron beam fluxes with the new broadening prescription and find that the Balmer lines are over-broadened at the densest times in the simulations. Adding many simultaneously heated and cooling model loops as a "multithread" model improves the agreement with the observations. We revisit the three-component phenomenological flare model of the YZ CMi Megaflare using recent and new RHD models. The evolution of the broadening, line flux ratios, and continuum flux ratios are well-reproduced by a multithread model with high-flux nonthermal electron beam heating, an extended decay phase model, and a "hot spot" atmosphere heated by an ultrarelativistic electron beam with reasonable filling factors: 0.1%, 1%, and 0.1% of the visible stellar hemisphere, respectively. The new modeling motivates future work to understand the origin of the extended gradual phase emission.

We present unique NUV observations of a well-observed X-class flare from NOAA 12087 obtained at Ond\v{r}ejov Observatory. The flare shows a strong white-light continuum but no detectable emission in the higher Balmer and Lyman lines. RHESSI and Fermi observations indicate an extremely hard X-ray spectrum and gamma-ray emission. We use the RADYN radiative hydrodynamic code to perform two type of simulations. One where an energy of 3 x 10^11 erg/cm^2/s is deposited by an electron beam with a spectral index of ~3 and a second where the same energy is applied directly to the photosphere. The combination of observations and simulations allow us to conclude that the white-light emission and the suppression or complete lack of hydrogen emission lines is best explained by a model where the dominant energy deposition layer is located in the lower layers of the solar atmosphere rather than the chromosphere.

Recent observations from the Interface Region Imaging Spectrograph (IRIS) appear to show impulsive brightenings in high temperature lines, which when combined with simultaneous ground based observations in H$\alpha$, appear co-spatial to Ellerman Bombs (EBs). We use the RADYN 1-dimensional radiative transfer code in an attempt to try and reproduce the observed line profiles and simulate the atmospheric conditions of these events. Combined with the MULTI/RH line synthesis codes, we compute the H$\alpha$, Ca II 8542~\AA, and Mg II h \& k lines for these simulated events and compare them to previous observations. Our findings hint that the presence of superheated regions in the photosphere ($>$10,000 K) is not a plausible explanation for the production of EB signatures. While we are able to recreate EB-like line profiles in H$\alpha$, Ca II 8542~\AA, and Mg II h \& k, we cannot achieve agreement with all of these simultaneously.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

We study the temporal evolution of the Na I D1 line profiles in the M3.9 flare SOL2014-06-11T21:03 UT, using high spectral resolution observations obtained with the IBIS instrument on the Dunn Solar Telescope combined with radiative hydrodynamic simulations. Our results show a significant increase in line core and wing intensities during the flare. The analysis of the line profiles from the flare ribbons reveal that the Na I D1 line has a central reversal with excess emission in the blue wing (blue asymmetry). We combine RADYN and RH simulations to synthesise Na I D1 line profiles of the flaring atmosphere and find good agreement with the observations. Heating with a beam of electrons modifies the radiation field in the flaring atmosphere and excites electrons from the ground state $\mathrm{3s~^2S}$ to the first excited state $\mathrm{3p~^2P}$, which in turn modifies relative population of the two states. The change in temperature and the population density of the energy states make the sodium line profile revert from absorption into emission. Analysis of the simulated spectra also reveals that the Na I D1 flare profile asymmetries are produced by the velocity gradients generated %and opacity effects in the lower solar atmosphere.

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from RHESSI. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

On April 23, 2014, the Swift satellite responded to a hard X-ray transient detected by its Burst Alert Telescope, which turned out to be a stellar flare from a nearby, young M dwarf binary DG~CVn. We utilize observations at X-ray, UV, optical, and radio wavelengths to infer the properties of two large flares. The X-ray spectrum of the primary outburst can be described over the 0.3-100 keV bandpass by either a single very high temperature plasma or a nonthermal thick-target bremsstrahlung model, and we rule out the nonthermal model based on energetic grounds. The temperatures were the highest seen spectroscopically in a stellar flare, at T$_{X}$ of 290 MK. The first event was followed by a comparably energetic event almost a day later. We constrain the photospheric area involved in each of the two flares to be $>$10$^{20}$ cm$^{2}$, and find evidence from flux ratios in the second event of contributions to the white light flare emission in addition to the usual hot, T$\sim$10$^{4}$K blackbody emission seen in the impulsive phase of flares. The radiated energy in X-rays and white light reveal these events to be the two most energetic X-ray flares observed from an M dwarf, with X-ray radiated energies in the 0.3-10 keV bandpass of 4$\times$10$^{35}$ and 9$\times$10$^{35}$ erg, and optical flare energies at E$_{V}$ of 2.8$\times$10$^{34}$ and 5.2$\times$10$^{34}$ erg, respectively. The results presented here should be integrated into updated modelling of the astrophysical impact of large stellar flares on close-in exoplanetary atmospheres.

We present a comprehensive study of the active dM4e star GJ 1243. We use previous observations and ground-based echelle spectroscopy to determine that GJ 1243 is a member of the Argus association of field stars, suggesting it is $\sim 30-50$ Myr old. We analyze eleven months of 1-minute cadence data from Kepler, presenting Kepler flare frequency distributions, as well as determining correlations between flare energy, amplitude, duration, and decay time. We find that the exponent $\alpha$ of the power-law flare energy distribution varies in time, primarily due to completeness of sample and the low frequency of high-energy flares. We also find a deviation from a single power law at high energy. We use ground-based spectroscopic observations simultaneous with the Kepler data to provide simultaneous photometric and spectroscopic analysis of three low-energy flares, the lowest-energy dMe flares with detailed spectral analysis to date on any star. The spectroscopic data from these flares extend constraints for radiative hydrodynamic (RHD) flare models to a lower energy regime than has previously been studied. We use this simultaneous spectroscopy and Kepler photometry to develop approximate conversions from the \textit{Kepler} bandpass to the traditional $U$ and $B$ bands. This conversion will be a critical factor in comparing any Kepler flare analyses to the canon of previous ground-based flare studies.

Ground- and space-based planet searches employing radial velocity techniques and transit photometry have detected thousands of planet-hosting stars in the Milky Way. The chemistry of these atmospheres is controlled by the shape and absolute flux of the stellar spectral energy distribution, however, flux distributions of relatively inactive low-mass stars are poorly known at present. To better understand exoplanets orbiting low-mass stars, we have executed a panchromatic (X-ray to mid-IR) study of the spectral energy distributions of 11 nearby planet hosting stars, the {\it Measurements of the Ultraviolet Spectral Characteristics of Low-mass Exoplanetary Systems} (MUSCLES) Treasury Survey. The MUSCLES program consists of contemporaneous observations at X-ray, UV, and optical wavelengths. We show that energetic radiation (X-ray and ultraviolet) is present from magnetically active stellar atmospheres at all times for stars as late as M5. Emission line luminosities of \ion{C}{4} and \ion{Mg}{2} are strongly correlated with band-integrated luminosities. We find that while the slope of the spectral energy distribution, FUV/NUV, increases by approximately two orders of magnitude form early K to late M dwarfs ($\approx$~0.01~to~1), the absolute FUV and XUV flux levels at their corresponding habitable zone distances are constant to within factors of a few, spanning the range 10~--~70 erg cm$^{-2}$ s$^{-1}$ in the habitable zone. Despite the lack of strong stellar activity indicators in their optical spectra, several of the M dwarfs in our sample show spectacular flare emission in their UV light curves. Finally, we interpret enhanced $L(line)$/$L_{Bol}$ ratios for \ion{C}{4} and \ion{N}{5} as tentative observational evidence for the interaction of planets with large planetary mass-to-orbital distance ratios ($M_{plan}$/$a_{plan}$) with the transition regions of their host stars.

We present a large dataset of high cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the NUV and optical continuum spectral evolution on timescales of ~1 second. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously-obtained low-resolution spectra during two moderate-sized flares in the dM4.5e star YZ CMi. By avoiding the spectral complexity within the broadband Johnson filters, the ULTRACAM filters are shown to characterize bona-fide continuum emission in the NUV, blue, and red wavelength regimes. The NUV/blue flux ratio in flares is equivalent to a Balmer jump ratio, and the blue/red flux ratio provides an estimate for the color temperature of the optical continuum emission. We present a new "color-color" relationship for these continuum flux ratios at the peaks of the flares. Using the RADYN and RH codes, we interpret the ULTRACAM filter emission using the dominant emission processes from a radiative-hydrodynamic flare model with a high nonthermal electron beam flux, which explains a hot, T~10,000 K, color temperature at blue-to-red optical wavelengths and a small Balmer jump ratio as are observed in moderate-sized and large flares alike. We also discuss the high time-resolution, high signal-to-noise continuum color variations observed in YZ CMi during a giant flare, which increased the NUV flux from this star by over a factor of 100.

In this talk, we discuss the formation of the near-ultraviolet and optical continuum emission in M dwarf flares through the formation of a dense, heated chromospheric condensation. Results are used from a recent radiative-hydrodynamic model of the response of an M dwarf atmosphere to a high energy flux of nonthermal electrons. These models are used to infer the charge density and optical depth in continuum emitting flare layers from spectra covering the Balmer jump and optical wavelength regimes. Future modeling and observational directions are discussed.

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Halpha and Ca II 8542 {\AA} lines are studied using high spatial, temporal and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1-m Solar Telescope. The temporal evolution of the Halpha line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum, and excess in the blue wing (blue asymmetry) after maximum. However, the Ca II 8542 {\AA} line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesise spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Halpha is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modifies the wavelength of the central reversal in the Halpha line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.

We present a unified computational framework which can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of stellar flares. For a summary, we refer to the paper.

We present the flare occurrence rates and starspot evolution for GJ 1245 A and B, two active M5 stars, based on nine months of Kepler short cadence observations, and four years of nearly continuous long cadence observations. The A component is separated from the B component by 7 arcseconds, and the stars are not resolved in the Kepler pipeline processing due to Kepler's large plate scale of 4 arcseconds/pixel. Analyzing the target pixel data, we have generated separate light curves for components A and B using the PyKE pixel response function modeling procedures, and note the effects of CCD saturation and non-linear response to high energy flares. In our sample, GJ 1245A and B exhibit an average of 3.0 and 2.6 flares per day, respectively. We introduce a new metric, $L_{fl}/L_{\mathrm{Kp}}$, to compare the flare rates between stars, and discuss this in the context of GJ 1245 A and B. Both stars exhibit starspot features that evolve on long time scales, with the slower rotating B component showing evidence of differential rotation. Intriguingly, the angular separation between the A and B component photocenters decreases during the four years of observations in a manner consistent with a shift in the position of the A photocenter due to the orbit of its unseen M8 companion (GJ 1245C), which is $\sim$94% less bright. Among the most detailed photometric studies of fully convective M dwarfs in a multiple system, these results provide an important constraint on stellar age-rotation-activity models.

We present the largest sample of flares ever compiled for a single M dwarf, the active M4 star GJ 1243. Over 6100 individual flare events, with energies ranging from $10^{29}$ to $10^{33}$ erg, are found in 11 months of 1-minute cadence data from Kepler. This sample is unique for its completeness and dynamic range. We have developed automated tools for finding flares in short-cadence Kepler light curves, and performed extensive validation and classification of the sample by eye. From this pristine sample of flares we generate a median flare template. This template shows that two exponential cooling phases are present during the white-light flare decay, providing fundamental constraints for models of flare physics. The template is also used as a basis function to decompose complex multi-peaked flares, allowing us to study the energy distribution of these events. Only a small number of flare events are not well fit by our template. We find that complex, multi-peaked flares occur in over 80% of flares with a duration of 50 minutes or greater. The underlying distribution of flare durations for events 10 minutes and longer appears to follow a broken power law. Our results support the idea that sympathetic flaring may be responsible for some complex flare events.

We analyze optical spectra of a two-ribbon, long duration C1.1 flare that occurred on 18 Aug 2011 within AR 11271 (SOL2011-08-18T15:15). The impulsive phase of the flare was observed with a comprehensive set of space-borne and ground-based instruments, which provide a range of unique diagnostics of the lower flaring atmosphere. Here we report the detection of enhanced continuum emission, observed in low-resolution spectra from 3600 \AA\ to 4550 \AA\ acquired with the Horizontal Spectrograph at the Dunn Solar Telescope. A small, $\le$0''.5 ($10^{15}$ cm$^2$) penumbral/umbral kernel brightens repeatedly in the optical continuum and chromospheric emission lines, similar to the temporal characteristics of the hard X-ray variation as detected by the Gamma-ray Burst Monitor (GBM) on the Fermi spacecraft. Radiative-hydrodynamic flare models that employ a nonthermal electron beam energy flux high enough to produce the optical contrast in our flare spectra would predict a large Balmer jump in emission, indicative of hydrogen recombination radiation from the upper flare chromosphere. However, we find no evidence of such a Balmer jump in the bluemost spectral region of the continuum excess. Just redward of the expected Balmer jump, we find evidence of a "blue continuum bump" in the excess emission which may be indicative of the merging of the higher order Balmer lines. The large number of observational constraints provides a springboard for modeling the blue/optical emission for this particular flare with radiative-hydrodynamic codes, which are necessary to understand the opacity effects for the continuum and emission line radiation at these wavelengths.

We analyzed Kepler short-cadence M dwarf observations. Spectra from the ARC 3.5m telescope identify magnetically active (H$\alpha$ in emission) stars. The active stars are of mid-M spectral type, have numerous flares, and well-defined rotational modulation due to starspots. The inactive stars are of early-M type, exhibit less starspot signature, and have fewer flares. A Kepler to U-band energy scaling allows comparison of the Kepler flare frequency distributions with previous ground-based data. M dwarfs span a large range of flare frequency and energy, blurring the distinction between active and inactive stars designated solely by the presence of H$\alpha$. We analyzed classical and complex (multiple peak) flares on GJ 1243, finding strong correlations between flare energy, amplitude, duration and decay time, with only a weak dependence on rise time. Complex flares last longer and have higher energy at the same amplitude, and higher energy flares are more likely to be complex. A power law fits the energy distribution for flares with log $E_{K_p} >$ 31 ergs, but the predicted number of low energy flares far exceeds the number observed, at energies where flares are still easily detectable, indicating that the power law distribution may flatten at low energy. There is no correlation of flare occurrence or energy with starspot phase; the flare waiting time distribution is consistent with flares occurring randomly in time; and the energies of consecutive flares are uncorrelated. These observations support a scenario where many independent active regions on the stellar surface are contributing to the observed flare rate.

The eruptive variable V838 Monocerotis gained notoriety in 2002 when it brightened nine magnitudes in a series of three outbursts and then rapidly evolved into an extremely cool supergiant. We present optical, near-IR, and mid-IR spectroscopic and photometric observations of V838 Monocerotis obtained between 2008 and 2012 at the Apache Point Observatory 3.5m, NASA IRTF 3m, and Gemini South 8m telescopes. We contemporaneously analyze the optical & IR spectroscopic properties of V838 Monocerotis to arrive at a revised spectral type L3 supergiant and effective temperature Teff~2000--2200 K. Because there are no existing optical observational data for L supergiants in the optical, we speculate that V838 Monocerotis may represent the prototype for L supergiants in this wavelength regime. We find a low level of Halpha emission present in the system, consistent with interaction between V838 Monocerotis and its B3V binary; however, we cannot rule out a stellar collision as the genesis event, which could result in the observed Halpha activity. Based upon a two-component blackbody fit to all wavelengths of our data, we conclude that, as of 2009, a shell of ejecta surrounded V838 Monocerotis at a radius of R=263+/-10 AU with a temperature of T=285+/-2 K. This result is consistent with IR interferometric observations from the same era and predictions from the Lynch et al. model of the expanding system, which provides a simple framework for understanding this complicated system.

In this splinter session, ten speakers presented results on solar and stellar activity and how the two fields are connected. This was followed by a lively discussion and supplemented by short, one-minute highlight talks. The talks presented new theoretical and observational results on mass accretion on the Sun, the activity rate of flare stars, the evolution of the stellar magnetic field on time scales of a single cycle and over the lifetime of a star, and two different approaches to model the radial-velocity jitter in cool stars that is due to the granulation on the surface. Talks and discussion showed how much the interpretation of stellar activity data relies on the sun and how the large number of objects available in stellar studies can extend the parameter range of activity models.

We present a homogeneous survey of line and continuum emission from near-ultraviolet (NUV) to optical wavelengths during twenty M dwarf flares with simultaneous, high cadence photometry and spectra. These data were obtained to study the white-light continuum components at bluer and redder wavelengths than the Balmer jump. Our goals were to break the degeneracy between emission mechanisms that have been fit to broadband colors of flares and to provide constraints for radiative-hydrodynamic (RHD) flare models that seek to reproduce the white-light flare emission. The main results from the continuum analysis are the following: 1) the detection of Balmer continuum (in emission) that is present during all flares and with a wide range of relative contributions to the continuum flux at bluer wavelengths than the Balmer jump; 2) a blue continuum at flare maximum that is linearly decreasing with wavelength from \lambda = 4000-4800\AA, matched by the spectral shape of hot, blackbody emission with typical temperatures of T_{BB}~9000-14,000 K; 3) a redder continuum apparent at wavelengths longer than H\beta\ (\lambda > 4900\AA) which becomes relatively more important to the energy budget during the late gradual phase. We calculate Balmer jump flux ratios and compare to RHD model spectra. The model ratios are too large and the blue-optical (\lambda = 4000-4800\AA) slopes are too red in both the impulsive and gradual decay phases of all twenty flares. This discrepancy implies that further work is needed to understand the heating at high column mass during dMe flares. (Abridged)

We study the formation and evolution of a failed filament eruption observed in NOAA active region 11121 near the southeast limb on November 6, 2010. We use a time series of SDO/AIA 304, 171, 131, 193, 335, 94 {\AA} images, SDO/HMI magnetograms, plus ROSA and ISOON H\alpha images, to study the erupting active region. We identify coronal loop arcades associated with a quadrupolar magnetic configuration, and show that the expansion and cancelation of the central loop arcade system over the filament is followed by the eruption of the filament. The erupting filament reveals a clear helical twist and develops a same sign of writhe in the form of inverse \gamma-shape. The observations support the "magnetic breakout" process with the eruption been triggered by quadrupolar reconnection in the corona. We suggest that the formation mechanism of the inverse \gamma-shape flux rope may be the MHD helical kink instability. The eruption has failed due to the large-scale, closed, overlying magnetic loop arcade that encloses the active region.

This contribution summarizes the splinter session "Non-thermal processes in coronae and beyond" held at the Cool Stars 17 workshop in Barcelona in 2012. It covers new developments in high energy non-thermal effects in the Earth's exosphere, solar and stellar flares, the diffuse emission in star forming regions and reviews the state and the challenges of the underlying atomic databases.

We present continuous, high-precision photometric monitoring data with 1 minute cadence of the dM3e flare star AD Leo with the {\it MOST} satellite. We observed 19 flares in 5.8 days, and find a flare frequency distribution that is similar to previous studies. The light curve reveals a sinusoidal modulation with period of $2.23^{+0.36}_{-0.27}$ days that we attribute to the rotation of a stellar spot rotating into and out of view. We see no correlation between the occurrence of flares and rotational phase, indicating that there may be many spots distributed at different longitudes, or possibly that the modulation is caused by varying surface coverage of a large polar spot that is viewed nearly pole-on. The data show no correlation between flare energy and the time since the previous flare. We use these results to reject a simple model in which all magnetic energy is stored in one active region and released only during flares.

This summary reports on papers presented at the Cool Stars-16 meeting in the splinter session "Solar and Stellar flares." Although many topics were discussed, the main themes were the commonality of interests, and of physics, between the solar and stellar flare communities, and the opportunities for important new observations in the near future.

We utilize the Sagittarius Window Eclipsing Extrasolar Planet Search (SWEEPS) HST/ACS dataset for a Deep Rapid Archival Flare Transient Search (DRAFTS) to constrain the flare rate toward the older stellar population in the Galactic bulge. During 7 days of monitoring 229,293 stars brighter than V=29.5, we find evidence for flaring activity in 105 stars between V=20 and V=28. We divided the sample into non-variable stars and variable stars whose light curves contain large-scale variability. The flare rate on variable stars is \sim 700 times that of non-variable stars, with a significant correlation between the amount of underlying stellar variability and peak flare amplitude. The flare energy loss rates are generally higher than those of nearby well-studied single dMe flare stars. The distribution of proper motions is consistent with the flaring stars being at the distance and age of the Galactic bulge. If they are single dwarfs, they span a range of \approx 1.0 - 0.25M\odot. A majority of the flaring stars exhibit periodic photometric modulations with P <3d. If these are tidally locked magnetically active binary systems, their fraction in the bulge is enhanced by a factor of \sim20 compared to the local value. These stars may be useful for placing constraints on the angular momentum evolution of cool close binary stars. Our results expand the type of stars studied for flares in the optical band, and suggest that future sensitive optical time-domain studies will have to contend with a larger sample of flaring stars than the M dwarf flare stars usually considered.

We present the first rates of flares from M dwarf stars in both red optical and near infrared (NIR) filters. We have studied ~50,000 M dwarfs from the SDSS Stripe 82 area, and 1,321 M dwarfs from the 2MASS Calibration Scan Point Source Working Database that overlap SDSS imaging fields. We assign photometric spectral types from M0 to M6 using (r-i) and (i-z) colors for every star in our sample. Stripe 82 stars each have 50-100 epochs of data, while 2MASS Calibration stars have ~1900 epochs. From these data we estimate the observed rates and theoretical detection thresholds for flares in eight photometric bands as a function of spectral type. Optical flare rates are found to be in agreement with previous studies, while the frequency per hour of NIR flare detections is found to be more than two orders of magnitude lower. An excess of small amplitude flux increases in all bands exhibits a power-law distribution, which we interpret as the result of flares below our detection thresholds. In order to investigate the recovery efficiency for flares in each filter, we extend a two-component flare model into the NIR. Quiescent M0-M6 spectral templates were used with the model to predict the photometric response of flares from u to Ks. We determine that red optical filters are sensitive to flares with u-band amplitudes >2 mag, and NIR filters to flares with delta u>4.5 mag. Our model predicts that M0 stars have the best color-contrast for J-band detections, but M4-M6 stars should show the highest rate of NIR flares with amplitudes of delta J ~0.01 mag. Characterizing flare rates and photometric variations at longer wavelengths is important for predicting the signatures of M dwarf variability in next-generation surveys, and we discuss their impact on surveys such as LSST.

Observations of extreme-ultraviolet (EUV) emission from an X-class solar flare that occurred on 2011 February 15 at 01:44 UT are presented, obtained using the EUV Variability Experiment (EVE) onboard the Solar Dynamics Observatory. The complete EVE spectral range covers the free-bound continua of H I (Lyman continuum), He I, and He II, with recombination edges at 91.2, 50.4, and 22.8 nm, respectively. By fitting the wavelength ranges blue-ward of each recombination edge with an exponential function, lightcurves of each of the integrated continua were generated over the course of the flare, as well as emission from the free-free continuum (6.5-37 nm). The He II 30.4 nm and Lyman-alpha 121.6 nm lines, and soft X-ray (0.1-0.8 nm) emission from GOES are also included for comparison. Each free-bound continuum was found to have a rapid rise phase at the flare onset similar to that seen in the 25-50 keV lightcurves from RHESSI, suggesting that they were formed by recombination with free electrons in the chromosphere. However, the free-free emission exhibited a slower rise phase seen also in the soft X-ray emission from GOES, implying a predominantly coronal origin. By integrating over the entire flare the total energy emitted via each process was determined. We find that the flare energy in the EVE spectral range amounts to at most a few per cent of the total flare energy, but EVE gives us a first comprehensive look at these diagnostically important continuum components.

We present the results of a campaign to monitor active M dwarfs using infrared spectroscopy, supplemented with optical photometry and spectroscopy. We detected 16 flares during nearly 50 hours of observations on EV Lac, AD Leo, YZ CMi, and VB8. The three most energetic flares also showed infrared emission, including the first reported detections of P\beta, P\gamma, He I 10830\AA and Br\gamma during an M dwarf flare. The strongest flare (\Delta u = 4.02 on EV Lac) showed emission from H\gamma, H\delta, He I 4471\AA, and Ca II K in the UV/blue and P\beta, P\gamma, P\delta, Br\gamma, and He I 10830\AA in the infrared. The weaker flares (\Delta u = 1.68 on EV Lac and \Delta U = 1.38 on YZ CMi) were only observed with photometry and infrared spectroscopy; both showed emission from P\beta, P\gamma, and He I 10830\AA. The strongest infrared emission line, P\beta, occurred in the active mid-M dwarfs with a duty cycle of ~3-4%. To examine the most energetic flare, we used the static NLTE radiative transfer code RH to produce model spectra based on a suite of one-dimensional model atmospheres. Using a hotter chromosphere than previous one-dimensional atmospheric models, we obtain line ratios that match most of the observed emission lines.

We present the results of an observational campaign which obtained high time cadence, high precision, simultaneous optical and IR photometric observations of three M dwarf flare stars for 47 hours. The campaign was designed to characterize the behavior of energetic flare events, which routinely occur on M dwarfs, at IR wavelengths to milli-magnitude precision, and quantify to what extent such events might influence current and future efforts to detect and characterize extrasolar planets surrounding these stars. We detected and characterized four highly energetic optical flares having U-band total energies of ~7.8x10^30 to ~1.3x10^32 ergs, and found no corresponding response in the J, H, or Ks bandpasses at the precision of our data. For active dM3e stars, we find that a ~1.3x10^32 erg U-band flare (delta Umax ~1.5 mag) will induce <8.3 (J), <8.5 (H), and <11.7 (Ks) milli-mags of a response. A flare of this energy or greater should occur less than once per 18 hours. For active dM4.5e stars, we find that a ~5.1x10^31 erg U-band flare (delta Umax ~1.6 mag) will induce <7.8 (J), <8.8 (H), and <5.1 (Ks) milli-mags of a response. A flare of this energy or greater should occur less than once per 10 hours. No evidence of stellar variability not associated with discrete flare events was observed at the level of ~3.9 milli-mags over 1 hour time-scales and at the level of ~5.6 milli-mags over 7.5 hour time-scales. We therefore demonstrate that most M dwarf stellar activity and flares will not influence IR detection and characterization studies of M dwarf exoplanets above the level of ~5-11 milli-mags, depending on the filter and spectral type. We speculate that the most energetic megaflares on M dwarfs, which occur at rates of once per month, are likely to be easily detected in IR observations with sensitivity of tens of milli-mags.

The white light during M dwarf flares has long been known to exhibit the broadband shape of a T~10,000 K blackbody, and the white light in solar flares is thought to arise primarily from Hydrogen recombination. Yet, a current lack of broad wavelength coverage solar-flare spectra in the optical/near-UV prohibits a direct comparison of the continuum properties to determine if they are indeed so different. New spectroscopic observations of a secondary flare during the decay of a megaflare on the dM4.5e star YZ CMi have revealed multiple components in the white-light continuum of stellar flares, including both a blackbody-like spectrum and a hydrogen recombination spectrum. One of the most surprising findings is that these two components are anti-correlated in their temporal evolution. We combine initial phenomenological modeling of the continuum components with spectra from radiative-hydrodynamic models to show that continuum veiling gives rise to the measured anti-correlation. This modeling allows us to use the components' inferred properties to predict how a similar spatially resolved, multiple-component white-light continuum might appear using analogies to several solar flare phenomena. We also compare the properties of the optical stellar flare white light to Ellerman bombs on the Sun.

We present sub-second, continuous-coverage photometry of three flares on the dM3.5e star, EQ Peg A, using custom continuum filters with WHT/ULTRACAM. These data provide a new view of flare continuum emission, with each flare exhibiting a very distinct light curve morphology. The spectral shape of flare emission for the two large-amplitude flares is compared with synthetic ULTRACAM measurements taken from the spectra during the large 'megaflare' event on a similar type flare star. The white light shape during the impulsive phase of the EQ Peg flares is consistent with the range of colors derived from the megaflare continuum, which is known to contain a Hydrogen recombination component and compact, blackbody-like components. Tentative evidence in the ULTRACAM photometry is found for an anti-correlation between the emission of these components.

We present a catalog of periodic stellar variability in the "Stripe 82" region of the Sloan Digital Sky Survey (SDSS). After aggregating and recalibrating catalog-level data from the survey, we ran a period-finding algorithm (Supersmoother) on all point-source lightcurves. We used color selection to identify systems that are likely to contain low-mass stars, in particular M dwarfs and white dwarfs. In total, we found 207 candidates, the vast majority of which appear to be in eclipsing binary systems. The catalog described in this paper includes 42 candidate M dwarf / white dwarf pairs, 4 white-dwarf pairs, 59 systems whose colors indicate they are composed of 2 M dwarfs and whose lightcurve shapes suggest they are in detached eclipsing binaries, and 28 M dwarf systems whose lightcurve shapes suggest they are in contact binaries. We find no detached systems with periods longer than 3 days, thus the majority of our sources are likely to have experienced orbital spin-up and enhanced magnetic activity. Indeed, twenty-six of twenty-seven M dwarf systems that we have spectra for show signs of chromospheric magnetic activity, far higher than the 24% seen in field stars of the same spectral type. We also find binaries composed of stars that bracket the expected boundary between partially and fully convective interiors, which will allow the measurement of the stellar mass-radius relationship across this transition. The majority of our contact systems have short orbital periods, with small variance (0.02 days) in the sample near the observed cutoff of 0.22 days. The accumulation of these stars at short orbital period suggests that the process of angular momentum loss, leading to period evolution, becomes less efficient at short periods. (Abridged)

We present a spectroscopic catalog of 70,841 visually inspected M dwarfs from the seventh data release of the Sloan Digital Sky Survey (SDSS). For each spectrum, we provide measurements of the spectral type, a number of molecular bandheads, and the H-alpha, H-beta, H-gamma, H-delta and Ca II K emission lines. In addition, we calculate the metallicity-sensitive parameter zeta and identify a relationship between zeta and the g-r and r-z colors of M dwarfs. We assess the precision of our spectral types (which were assigned by individual examination), review the bulk attributes of the sample, and examine the magnetic activity properties of M dwarfs, in particular those traced by the higher order Balmer transitions. Our catalog is cross-matched to Two Micron All Sky Survey (2MASS) infrared data, and contains photometric distances for each star. Lastly, we identify eight new late-type M dwarfs that are possibly within 25 pc of the Sun. Future studies will use these data to thoroughly examine magnetic activity and kinematics in late-type M dwarfs and examine the chemical and dynamical history of the local Milky Way.

Time-resolved databases with large spatial coverage are quickly becoming a standard tool for all types of astronomical studies. We report preliminary results from our search for stellar flares in the 2MASS calibration fields. A sample of 4343 M dwarfs, spatially matched between the SDSS and the 2MASS calibration fields, each with hundreds to thousands of epochs in near infrared bandpasses, is analyzed using a modified Welch-Stetson index to characterize the variability. A Monte Carlo model was used to assess the noise of the variability index. We find significnat residuals above the noise with power-law slopes of -3.37 and -4.05 for our JH and HKs distributions respectively. This is evidence for flares being observed from M dwarfs in infrared photometry.

We present a spectroscopic catalog of 70,841 visually inspected M dwarfs from the seventh data release (DR7) of the Sloan Digital Sky Survey (SDSS). For each spectrum, we provide measurements of the spectral type, a number of molecular bandheads, and the H-alpha, H-beta, H-gamma, H-delta and Ca II K emission lines. In addition, we calculate the metallicity-sensitive parameter zeta and 3D space motions for most of the stars in the sample. Our catalog is cross-matched to Two Micron All Sky Survey (2MASS) infrared data, and contains photometric distances for each star. Future studies will use these data to thoroughly examine magnetic activity and kinematics in late-type M dwarfs and examine the chemical and dynamical history of the local Milky Way.

M dwarfs are known to flare on timescales from minutes to hours, with flux increases of several magnitudes in the blue/near-UV. These frequent, powerful events, which are caused by magnetic reconnection, will have a strong observational signature in large, time-domain surveys. The radiation and particle fluxes from flares may also exert a significant influence on the atmospheres of orbiting planets, and affect their habitability. We present a statistical model of flaring M dwarfs in the Galaxy that allows us to predict the observed flare rate along a given line of sight for a particular survey depth and cadence. The parameters that enter the model are the Galactic structure, the distribution of magnetically active and inactive M dwarfs, and the flare frequency distribution (FFD) of both populations. The FFD is a function of spectral type, activity, and Galactic height. Although inactive M dwarfs make up the majority of stars in a magnitude-limited survey, the FFD of inactive stars is very poorly constrained. We have organized a flare monitoring campaign comprising hundreds of hours of new observations from both the ground and space to better constrain flare rates. Incorporating the new observations into our model provides more accurate predictions of stellar variability caused by flares on M dwarfs. We pay particular attention to the likelihood of flares appearing as optical transients (i.e., host star not seen in quiescent data).

A white paper prepared for the Space Studies Board, National Academy of Sciences (USA), for its Decadal Survey of Solar and Space Physics (Heliophysics), reviewing and encouraging studies of flare physics in the chromosphere.

M dwarfs produce explosive flare emission in the near-UV and optical continuum, and the mechanism responsible for this phenomenon is not well-understood. We present a near-UV/optical flare spectrum from the rise phase of a secondary flare, which occurred during the decay of a much larger flare. The newly formed flare emission resembles the spectrum of an early-type star, with the Balmer lines and continuum in absorption. We model this observation phenomonologically as a temperature bump (hot spot) near the photosphere of the M dwarf. The amount of heating implied by our model (\Delta T_phot ~ 16,000K) is far more than predicted by chromospheric backwarming in current 1D RHD flare models (\Delta T_phot ~ 1200K).

We have identified 63 flares on M dwarfs from the individual component spectra in the Sloan Digital Sky Survey using a novel measurement of emission line strength called the Flare Line Index. Each of the ~38,000 M dwarfs in the SDSS low mass star spectroscopic sample of West et al. was observed several times (usually 3-5) in exposures that were typically 9-25 minutes in duration. Our criteria allowed us to identify flares that exhibit very strong H-alpha and H-beta emission line strength and/or significant variability in those lines throughout the course of the exposures. The flares we identified have characteristics consistent with flares observed by classical spectroscopic monitoring. The flare duty cycle for the objects in our sample is found to increase from 0.02% for early M dwarfs to 3% for late M dwarfs. We find that the flare duty cycle is larger in the population near the Galactic plane and that the flare stars are more spatially restricted than the magnetically active but non-flaring stars. This suggests that flare frequency may be related to stellar age (younger stars are more likely to flare) and that the flare stars are younger than the mean active population.

On UT 2009 January 16, we observed a white light megaflare on the dM4.5e star YZ CMi as part of a long-term spectroscopic flare-monitoring campaign to constrain the spectral shape of optical flare continuum emission. Simultaneous U-band photometric and 3350A-9260A spectroscopic observations were obtained during 1.3 hours of the flare decay. The event persisted for more than 7 hours and at flare peak, the U-band flux was almost 6 magnitudes brighter than in the quiescent state. The properties of this flare mark it as one of the most energetic and longest-lasting white light flares ever to be observed on an isolated low-mass star. We present the U-band flare energetics and a flare continuum analysis. For the first time, we show convincingly with spectra that the shape of the blue continuum from 3350A to 4800A can be represented as a sum of two components: a Balmer continuum as predicted by the Allred et al radiative hydrodynamic flare models and a T ~ 10,000K blackbody emission component as suggested by many previous studies of the broadband colors and spectral distributions of flares. The areal coverage of the Balmer continuum and blackbody emission regions vary during the flare decay, with the Balmer continuum emitting region always being significantly (~3-16 times) larger. These data will provide critical constraints for understanding the physics underlying the mysterious blue continuum radiation in stellar flares.

(Abridged) Classical Be stars occasionally transition from having a gaseous circumstellar disk (''Be phase'') to a state in which all observational evidence for the presence of these disks disappears (''normal B-star phase''). We present one of the most comprehensive spectropolarimetric views to date of such a transition for two Be stars, pi Aquarii and 60 Cygni. 60 Cyg's disk loss episode was characterized by a monotonic decrease in emission strength over a time-scale of 1000 days, consistent with the viscous time-scale of the disk, assuming alpha is 0.14. pi Aqr's disk loss was episodic in nature and occurred over a time-scale of 2440 days. An observed time lag between the behavior of the polarization and H-alpha in both stars indicates the disk clearing proceeded in an ''inside-out'' manner. We determine the position angle of the intrinsic polarization to be 166.7 +/- 0.1 degrees for pi Aqr and 107.7 +/- 0.4 degrees for 60 Cyg, and model the observed polarization during the quiescent diskless phase of each star to determine the interstellar polarization along the line of sight. Minor outbursts observed during the quiescent phase of each star shared similar lifetimes as those previously reported for mu Cen, suggesting that the outbursts represent the injection and subsequent viscous dissipation of individual blobs of material into the inner circumstellar environments of these stars. We also observe deviations from the mean intrinsic polarization position angle during polarization outbursts in each star, indicating deviations from axisymmetry. We propose that these deviations might be indicative of the injection (and subsequent circularization) of new blobs into the inner disk, either in the plane of the bulk of the disk material or in a slightly inclined (non-coplanar) orbit.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present a flare rate analysis of 50,130 M dwarf light curves in SDSS Stripe 82. We identified 271 flares using a customized variability index to search ~2.5 million photometric observations for flux increases in the u- and g-bands. Every image of a flaring observation was examined by eye and with a PSF-matching and image subtraction tool to guard against false positives. Flaring is found to be strongly correlated with the appearance of H-alpha in emission in the quiet spectrum. Of the 99 flare stars that have spectra, we classify 8 as relatively inactive. The flaring fraction is found to increase strongly in stars with redder colors during quiescence, which can be attributed to the increasing flare visibility and increasing active fraction for redder stars. The flaring fraction is strongly correlated with |Z| distance such that most stars that flare are within 300 pc of the Galactic plane. We derive flare u-band luminosities and find that the most luminous flares occur on the earlier-type M dwarfs. Our best estimate of the lower limit on the flaring rate (averaged over Stripe 82) for flares with \Delta u \ge 0.7 magnitudes on stars with u < 22 is 1.3 flares hour^-1 square degree^-1 but can vary significantly with the line-of-sight.

Asteroid families, traditionally defined as clusters of objects in orbital parameter space, often have distinctive optical colors. We show that the separation of family members from background interlopers can be improved with the aid of SDSS colors as a qualifier for family membership. Based on an ~88,000 object subset of the Sloan Digital Sky Survey Moving Object Catalog 4 with available proper orbital elements, we define 37 statistically robust asteroid families with at least 100 members using a simple Gaussian distribution model in both orbital and color space. The interloper rejection rate based on colors is typically ~10% for a given orbital family definition, with four families that can be reliably isolated only with the aid of colors. About 50% of all objects in this data set belong to families, and this fraction varies from about 35% for objects brighter than an H magnitude of 13 and rises to 60% for objects fainter than this. The fraction of C-type objects in families decreases with increasing H magnitude for H > 13, while the fraction of S-type objects above this limit remains effectively constant. This suggests that S-type objects require a shorter timescale for equilibrating the background and family size distributions via collisional processing. The size distributions for 15 families display a well-defined change of slope and can be modeled as a "broken" double power-law. Such "broken" size distributions are twice as likely for S-type familes than for C-type families, and are dominated by dynamically old families. The remaining families with size distributions that can be modeled as a single power law are dominated by young families. When size distribution requires a double power-law model, the two slopes are correlated and are steeper for S-type families.

The first results from a near-contemporaneous optical and infrared spectroscopic observing program designed to probe the detailed density structure of classical Be circumstellar disks are presented. We report the discovery of asymmetrical infrared emission lines of He I, O I, Fe II, and the Brackett, Paschen, and Pfund series lines of H I which exhibit an opposite V/R orientation (V $>$ R) to that observed for the optical Balmer H$\alpha$ line (V $<$ R) in the classical Be star $\zeta$ Tau. We interpret these data as evidence that the density wave which characterizes $\zeta$ Tau's disk has a significantly different average azimuthal morphology in the inner disk region as compared to the outer disk region. A follow-up multi-wavelength observational campaign to trace the temporal evolution of these line profile morphologies, along with detailed theoretical modeling, is suggested to test this hypothesis.