Kepler K2 long cadence data are used to study white light flares in a sample of 45 L dwarfs. We identified 11 flares on 9 L dwarfs with equivalent durations of (1.3 - 198) hr and total (UV/optical/IR) energies of $\geq$0.9 $\times$ 10$^{32}$ erg. Two superflares with energies of $>$10$^{33}$ erg were detected on an L5 dwarf: this is the coolest object so far on which flares have been identified. The larger superflare on this L5 dwarf has an energy of 4.6$\times$ 10$^{34}$ ergs and an amplitude of $>$300 times the photospheric level: so far, this is the largest amplitude flare detected by the $Kepler/K2$ mission. The next coolest star on which we identified a flare was an L2 dwarf: 2MASS J08585891+1804463. Combining the energies of all the flares which we have identified on 9 L dwarfs with the total observation time which was dedicated by $Kepler$ to all 45 L dwarfs, we construct a composite flare frequency distribution (FFD). The FFD slope is quite shallow (-0.51$\pm$0.17), consistent with earlier results reported by Paudel et al. (2018) for one particular L0 dwarf, for which the FFD slope was found to be -0.34. Using the composite FFD, we predict that, in early and mid-L dwarfs, a superflare of energy 10$^{33}$ erg occurs every 2.4 years and a superflare of energy 10$^{34}$ erg occurs every 7.9 years. Analysis of our L dwarf flares suggests that magnetic fields of $\geq$0.13-1.3 kG are present on the stellar surface: such fields could suppress Type II radio bursts.

Houdebine et al (2017: H17) combined CaII data with projected rotational velocities (v sin i) to construct rotation-activity correlations (RAC) in K-M dwarfs. The RAC slopes were used to argue that a transition between dynamo modes occurs at a spectral type between M2 and M3. H17 suggested that the dynamo transition corresponds to a transition to complete convection (TTCC). An independent study of GAIA data led Jao et al (2018) to suggest that the TTCC sets in near M3.0V, close to the H17 result. However, the changes in a star which cause TTCC signatures in GAIA data might not lead to changes in CaII emission at an identical spectral type: the latter are also affected by magnetic effects which depend on certain properties of convection in the core. Here, we use CaII emission fluxes in a sample of ~600 M dwarfs, and attempt to narrow down the transition from one dynamo mode to another: rather than relying on RAC slopes, we quantify how the CaII emission flux varies with spectral type to identify steps where the flux decreases significantly across a narrow range of spectral types. We suggest that the dynamo mode transition may be narrowed down to between M2.1 and M2.3. This is close to, but earlier than, the TTCC location identified by Jao et al (2018). We suggest that the transition in dynamo mode may be related to the existence of a small convective core which occurs for a finite time interval in certain low mass stars.

The onset of cool massive winds in evolved giants is correlated with an evolutionary feature on the red giant branch known as the bump. Also at the bump, shear instability in the star leads to magnetic fields that occur preferentially on small length scales. Pneuman (1983) has suggested that the emergence of small scale flux tubes in the Sun can give rise to enhanced acceleration of the solar wind as a result of plasmoid acceleration (the melon seed mechanism). In this paper, we examine the Pneuman formalism to determine if it may shed some light on the process that drives mass loss from stars above the bump. Because we do not currently have detailed information for some of the relevant physical parameters, we are not yet able to derive a detailed model. Instead, our goal in this paper is to explore a proof of concept. Using parameters that are known to be plausible in cool giants, we find that the total mass loss rate from such stars can be replicated. Moreover, we find that the radial profile of the wind speed in such stars can be steep or shallow depending on the fraction of the mass loss which is contained in the plasmoids. This is consistent with empirical data which indicate that the velocity profiles of winds from cool giants range from shallow to steep.

Empirical correlations between stellar parameters such as rotation or radius and magnetic activity diagnostics require estimates of the effective temperatures and the stellar radii. The aim of this study is to propose simple methods that can be applied to large samples of stars in order to derive estimates of the stellar parameters. Good empirical correlations between Red/Infra-Red colours (e.g. (R-I)$_C$) and effective temperatures have been well established for a long time (e.g. Veeder 1974, Bessell 1979, Leggett 1992). The more recent (R-I)$_C$ colour-$T_{eff}$ correlation using the data of Mann et al. (2015: hereafter M15) and Boyajian et al. (2012: hereafter B12) shows that this colour can be applied as a temperature estimate for large samples of stars. We find that the mean scatter in $T_{eff}$ relative to the (R-I)$_C$-$T_{eff}$ relationship of B12 and M15 data is only $\pm 3\sigma =$44.6~K for K dwarfs and $\pm 3\sigma =$39.4~K for M dwarfs. These figures are small and show that the (R-I)$_C$ colour can be used as a first guess effective temperature estimator for K and M dwarfs. We derive effective temperatures for about 1910 K and M dwarfs using the calibration of (R-I)$_C$ colour-$T_{eff}$ from B12 and M15 data. We also compiled $T_{eff}$ and metallicity measurements available in the literature using the VizieR database. We determine $T_{eff}$ for 441 stars with previously unknown effective temperatures. We also identified 21 new spectroscopic binaries and 1 triple system from our high resolution spectra.

Jurcak et al (2018) have reported that, in a sample of more than 100 umbral cores in sunspots, the umbral-penumbral boundary (UPB) is characterized by a remarkably narrowly-defined numerical value (1867 G) of the vertical component of the magnetic field. Gough and Tayler (1966), in their study of magneto-convection, showed that the onset of convection in the presence of a magnetic field is controlled by a parameter {\delta} which also depends on the vertical component of the field. Combining the Jurcak et al result with various empirical models of sunspots leads us to propose the following hypothesis: the UPB occurs where the vertical field is strong enough to increase the effective adiabatic temperature gradient by at least 100% above its non-magnetic value.

Type II radio bursts are observed in the Sun in association with many coronal mass ejections (CME's. In view of this association, there has been an expectation that, by scaling from solar flares to the flares which are observed on M dwarfs, radio emission analogous to solar Type II bursts should be detectable in association with M dwarf flares. However, several surveys have revealed that this expectation does not seem to be fulfilled. Here we hypothesize that the presence of larger global field strengths in low-mass stars, suggested by recent magneto-convective modeling, gives rise to such large Alfven speeds in the corona that it becomes difficult to satisfy the conditions for the generation of Type II radio bursts. As a result, CME's propagating in the corona/wind of a flare stars are expected to be "radio-quiet" as regards Type II bursts. In view of this, we suggest that, in the context of Type II bursts, scaling from solar to stellar flares is of limited effectiveness.

We observed strong superflares (defined as flares with energy in excess of 10^33 erg) on three late-M dwarfs: 2MASS J08315742+2042213 (hereafter 2M0831+2042; M7 V), 2MASS J08371832+2050349 (hereafter 2M0837+2050; M8 V) and 2MASS J08312608+2244586 (hereafter 2M0831+2244; M9 V). 2M0831+2042 and 2M0837+2050 are members of the young (~700 Myr) open cluster Praesepe. The strong superflare on 2M0831+2042 has an equivalent duration (ED) of 13.7 hr and an estimated energy of 1.3 X 10^35 erg. We observed five superflares on 2M0837+2050, on which the strongest superflare has an ED of 46.4 hr and an estimated energy of 3.5 X 10^35 erg. This energy is larger by 2.7 orders of magnitude than the largest flare observed on the older (7.6 Gyr) planet-hosting M8 dwarf TRAPPIST-1. Furthermore, we also observed five superflares on 2M0831+2244 which is probably a field star. The estimated energy of the strongest superflare on 2M0831+2244 is 6.1 X 10^34 erg. 2M0831+2042, 2M0837+2050 and 2MASS J0831+2244 have rotation periods of 0.556\pm0.002, 0.193\pm0.000 and 0.292\pm0.001 d respectively, which are measured by using K2 light curves. We compare the flares of younger targets with those of TRAPPIST-1 and discuss the possible impacts of such flares on planets in the habitable zone of late-M dwarfs.

We construct evolutionary models of Trappist-1 in which magnetic fields impede the onset of convection according to a physics-based criterion. In the models that best fit all observational constraints, the photospheric fields in Tr-1 are found to be in the range 1450-1700 G. These are weaker by a factor of about 2 than the fields we obtained in previous magnetic models of two other cool dwarfs (GJ65A/B). Our results suggest that Tr-1 possesses a global poloidal field which is some one hundred times stronger than in the Sun. In the context of exoplanets in orbit around Tr-1, the strong poloidal fields on the star may help to protect the planets from the potentially destructive effects of coronal mass ejections. This, in combination with previous arguments about beneficial effects of flare photons in ultraviolet and visible portions of the spectrum, suggests that conditions on Tr-1 are not necessarily harmful to life on a planet in the habitable zone of Tr-1.

On planets near M dwarfs, photosynthesis (PS) will occur with an effectiveness which depends on the supply of visible photons with wavelengths between 400 and 700 nm. In this paper, we quantify the effectiveness of PS in two contexts which are relevant for M dwarfs. First, using photons from an M dwarf in its quiescent non-flaring state, we find that PS on an M dwarf planet in the HZ of its parent star is less effective than on Earth by a factor of 10 for a flare star with mid-M spectral type. For a flare star with late-M spectral type, PS effectiveness is smaller than on Earth by a factor of 100 or more. Second, using photons which are incident on the HZ planet during flares, we find that PS effectiveness can increase by factors of 5-20 above the quiescent values. In the case of a flare star with mid-M spectral type, we find that the PS effectiveness during a flare can increase up to as much as 50-60 percent of the values on Earth. However, for a late-M flare star, even during flares, the PS effectiveness remains almost one order of magnitude smaller than on Earth. We suggest that for biological processes on M dwarf planets, the stellar activity cycle may replace the orbital period as the year.

We present photometric measurements of two superflares observed on a very young brown dwarf CFHT-BD-Tau 4, observed during Campaign 13 of the \textit{Kepler K2} mission. The stronger of the two superflares brightened by a factor of $\sim$48 relative to the quiescent photospheric level, with an increase in \textit{Kepler} magnitude $\Delta \tilde{K_{p}}$ = -4.20. It has an equivalent duration of $\sim$107 hour, a flare duration of 1.7 day, and an estimated total bolometric (ultraviolet/optical/infrared) energy up to 2.1 $\times$ 10$^{38}$ erg. The weaker of the two superflares is a complex (multipeaked) flare with an estimated total bolometric (UV/optical/IR) energy up to 4.7 $\times$ 10$^{36}$ erg. They are the strongest flares observed on any brown dwarf so far. The flare energies are strongly dependent on the value of visual extinction parameter $A_{V}$ used for extinction correction. If we apply a solar flare-model to interpret the two superflares, we find that the magnetic fields are required to be stronger by as much as an order of magnitude than previous reports of field measurements in CFHT-BD-Tau 4 by Reiners et al. (2009b). On the other hand, if we interpret our data in terms of accretion, we find that the requisite rate of accretion for the stronger superflare exceeds the rates which have been reported for other young brown dwarfs.

GJ 65 is a M dwarf binary system consisting of the two flare stars BL Cet (GJ 65A) and UV Cet (GJ 65B). Two teams of investigators have recently reported total magnetic fluxes corresponding to fields of 4.5 and 5.2 kG for GJ65A, and 5.8 and 6.7 kG for GJ65B: for each component, the magnetic results obtained by the two teams agree with each other within 1{\sigma}. For the first time, we can directly compare the predictions of our magneto-convective models, based on fitting observed stellar parameters, with measured field strengths. We find that our models agree with the observed field strengths provided the effects of finite conductivity are accounted for. Thus, GJ65 provides us an opportunity to use observations of field strengths to distinguish between the predictions of our models that assume perfect electrical conductivity and those that allow for finite conductivity.

We report the white light flare rates for 10 ultracool dwarfs (UCDs) using \textit{Kepler K2} short cadence data. Among our sample stars, two have spectral type M6, three are M7, three are M8 and two are L0. Most of our targets are old low mass stars. We identify a total of 283 flares in all of the stars in our sample, with \textit{Kepler} energies in the range log \textit{E$_{Kp}$} $\sim$(29 - 33.5) erg. Using the maximum-likelihood method of line fitting, we find that the flare frequency distribution (FFD) for each star in our sample follows a power law with slope -$\alpha$ in range -(1.3-2.0). We find that cooler objects tend to have shallower slopes. For some of our targets, the FFD follows either a broken power law, or a power law with an exponential cutoff. For the L0 dwarf 2MASS J12321827-0951502, we find a very shallow slope (-$\alpha$ $=$ -1.3) in the \textit{Kepler} energy range (0.82-130)$\times$10$^{30}$ erg: this L0 dwarf has flare rates which are comparable to the rates of high energy flares in stars of earlier spectral types. In addition, we report photometry of two superflares: one on the L0 dwarf 2MASS J12321827-0951502 and another on the M7 dwarf 2MASS J08352366+1029318. In case of 2MASS J12321827-0951502, we report a flare brightening by a factor of $\sim$144 relative to the quiescent photospheric level. Likewise, for 2MASS J08352366+1029318, we report a flare brightening by a factor of $\sim$60 relative to the quiescent photospheric level. These two superflares have bolometric (UV/optical/infrared) energies 3.6 $\times$ 10$^{33}$ erg and 8.9 $\times$ 10$^{33}$ erg respectively, while the FWHM time scales are very short, $\sim$2 minutes. We find that the M8 star TRAPPIST-1 is more active than the M8.5 dwarf: 2M03264453+1919309, but less active than another M8 dwarf (2M12215066-0843197).

Observations of solar and stellar flares have revealed the presence of power law dependences between the flare energy and the time interval between flares. Various models have been proposed to explain these dependences, and to explain the numerical value of the power law indices. Here, we propose a model in which convective flows in granules force the foot-points of coronal magnetic loops, which are frozen-in to photospheric gas, to undergo a random walk. In certain conditions, this can lead to a twist in the loop, which drives the loop unstable if the twist exceeds a critical value. The possibility that a solar flare is caused by such a twist-induced instability in a loop has been in the literature for decades. Here, we quantify the process in an approximate way with a view to replicating the power-law index. We find that, for relatively small flares, the random walk twisting model leads to a rather steep power law slope which agrees very well with the index derived from a sample of 56,000+ solar X-ray flares reported by the GOES satellites. For relatively large flares, we find that the slope of the power law is shallower. The empirical power law slopes reported for flare stars also have a range which overlaps with the slopes obtained here. We suggest that in the coolest stars, a significant change in slope should occur when the frozen-flux assumption breaks down due to low electrical conductivity.

We have previously reported on models of low-mass stars in which the presence of inflated radii is ascribed to magnetic fields which impede the onset of convection (e.g. MacDonald & Mullan [2017a] and citations therein). Some of our magneto-convection models have been criticized because, when they were first reported by Mullan & MacDonald (2001), the deep interior fields were found to be very large (50-100 MG). Such large fields are now known to be untenable. For example, Browning et al. (2016) used stability arguments to suggest that interior fields in low-mass stars cannot be larger than ~1 MG. Moreover, 3D models of turbulent stellar dynamos suggest that fields generated in low-mass interiors may be not much stronger than 10-20 kG (Browning 2008). In the present paper, we present magneto-convective models of inflated low-mass stars in which the interior fields are not permitted to be stronger than 10 kG. These models are used to fit empirical data for 15 low-mass stars for which precise masses and radii have been measured. We show that our 10 kG magneto-convective models can replicate the empirical radii and effective temperatures for 14 of the stars. In the case of the remaining star (in the Praesepe cluster), two different solutions have been reported in the literature. We find that one of these solutions (by Gillen et al. 2017) can be fitted well with our model using the nominal age of Praesepe (800 Myr). However, the second solution (by Kraus et al. 2017) cannot be fitted unless the star's age is assumed to be much younger (~ 150 Myr).

We use Kepler K2 Campaign 4 short-cadence (one-minute) photometry to measure white light flares in the young, moving group brown dwarfs 2MASS J03350208+2342356 (2M0335+23) and 2MASS J03552337+1133437 (2M0355+11), and report on long-cadence (thirty-minute) photometry of a superflare in the Pleiades M8 brown dwarf CFHT-PL-17. The rotation period (5.24 hr) and projected rotational velocity ($45$ km s$^{-1}$) confirm 2M0335+23 is inflated ($R \ge 0.20 R_\odot$) as predicted for a $0.06M_\odot$, 26-Myr old brown dwarf $\beta $Pic moving group member. We detect 22 white light flares on 2M0335+23. The flare frequency distribution follows a power-law distribution with slope $-\alpha = -1.8 \pm 0.2$ over the range $10^{31}$ to $10^{33}$ erg. This slope is similar to that observed in the Sun and warmer flare stars, and is consistent with lower energy flares in previous work on M6-M8 very-low-mass stars; taken the two datasets together, the flare frequency distribution for ultracool dwarfs is a power law over 4.3 orders of magnitude. The superflare ($2.6\times10^{34}$ erg) on CFHT-PL-17 shows higher energy flares are possible. We detect no flares down to a limit of $2 \times 10^{30}$ erg in the nearby L$5\gamma$ AB Dor Moving Group brown dwarf 2M0355+11, consistent with the view that fast magnetic reconnection is suppressed in cool atmospheres. We discuss two multi-peaked flares observed in 2M0335+23, and argue that these complex flares can be understood as sympathetic flares, in which a fast-mode MHD waves similar to EUV waves in the Sun trigger magnetic reconnection in different active regions.

LSPM J1314+1320 (=NLTT 33370) is a binary star system consisting of two nearly identical pre-main sequence stars of spectral type M7. The system is remarkable among ultracool dwarfs for being the most luminous radio emitter over the widest frequency range. Masses and luminosities are at first sight consistent with the system being coeval at age ~80 Myr according to standard (non-magnetic) evolutionary models. However, these models predict an average effective temperature, 2950 +/- 5 K, which is 180 K hotter than the empirical value. Thus, the empirical radii are oversized relative to the standard models by ~13%. We demonstrate that magnetic stellar models can account quantitatively for the oversizing. As a check on our models, we note that the radio emission limits the surface magnetic field strengths: the limits depend on identifying the radio emission mechanism. We find that the field strengths required by our magnetic models are too strong to be consistent with gyrosynchrotron emission, but are consistent with electron cyclotron maser emission.

We study the rotation-activity correlations (RACs) in a sample stars from spectral type dK4 to dM4. We study RACs using chromospheric data and coronal data. We study the Ca\,{\sc ii} line surface fluxes-$P/\sin i$ RACs. We fit the RACs with linear homoscedastic and heteroscedastic regression models. We find that these RACs differ substantially from one spectral sub-type to another. For dM3 and dM4 stars, we find that the RACs cannot be described by a simple model, but instead that there may exist two distinct RAC behaviors for the low activity and the high activity stellar sub-samples respectively. Although these results are preliminary and will need confirmation, the data suggest that these distinct RACs may be associated with different dynamo regimes. We also study R'$_{HK}$ as a function of the Rossby number R$_{0}$. We find that: (i) For dK4 stars, we confirm R'$_{HK}$ as a function of R$_{0}$ agrees well with previous results for F-G-K stars. (ii) In dK6, dM2, dM3 and dM4 stars, we find that, at a given R$_{0}$, the values of R'$_{HK}$ lie a factor of 3, 10, 20 and 90 respectively below the F-G-K RAC. Our results suggest a significant decrease in the efficiency of the dynamo mechanism(s) as regards chromospheric heating before and at dM3, i.e. before and at the TTCC. We also show that the ratio of coronal heating to chromospheric heating $L_{X}/L_{HK}$ increases by a factor of 100 between dK4 and dM4 stars.

Two eclipsing binaries in the USco association have recently yielded precise values of masses and radii for 4 low-mass members of the association. Standard evolution models would require these dM4.5 - dM5 stars to have ages which are younger than the ages of more massive stars in the association by factors which appear (in extreme cases) to be as large as ~3. Are the stars in the association therefore non-coeval? We suggest that the answer is No: by incorporating the effects of magnetic inhibition of convective onset, we show that the stars in USco can be restored to coevality provided that the 4 low-mass member stars have vertical surface fields in the range 200 - 700 G. Fields of such magnitude have already been measured on the surface of certain solar-type stars in other young clusters.

KIC 7177553 is a quadruple system containing two binaries of orbital periods 16.5 and 18 d. All components have comparable masses and are slowly rotating, non-evolved stars of spectral type near G2V. The longer period binary is eclipsing with component masses and radii, M1 = 1.043 +/- 0.014 M_sun, R1 = 0.940 +/- 0.005 R_sun, and M2 = 0.986 +/- 0.015 M_sun, R2 = 0.941 +/- 0.005 R_sun. The essentially equal radii measurements are inconsistent with the two stars being on the main sequence at the same age using standard stellar evolution models. Instead a consistent scenario is found if the stars are in their pre-main sequence phase of evolution and have age 33 - 36 Myr. Such a young age has important implications for the detectability of the massive planet indicated by eclipse time variations.

Reliable determination of Rotation-Activity Correlation (RAC's) depends on precise measurements of the following stellar parameters: $T_{eff}$, parallax, radius, metallicity, and rotational speed $v\sin i$. In this paper, our goal is to focus on the determination of these parameters for a sample of K and M dwarfs. In a future paper (Paper II), we will combine our rotational data with activity data in order to construct RAC's. Here, we report on a determination of effective temperatures based on the (R-I)$_C$ color from the calibrations of Mann et al. (2015) and Kenyon \& Hartmann (1995) for four samples of late-K, dM2, dM3 and dM4 stars. We also determine stellar parameters ($T_{eff}$, $log(g)$ and [M/H]) using the PCA-based inversion technique for a sample of 105 late-K dwarfs. We compile all effective temperatures from the literature for this sample. We determine empirical radius-[M/H] correlations in our stellar samples. This allows us to propose new effective temperatures, stellar radii, and metallicities for a large sample of 612 late-K and M dwarfs. Our mean radii agree well with those of Boyajian et al. (2012). We analyze HARPS and SOPHIE spectra of 105 late-K dwarfs, and have detected $v\sin i$ in 92 stars. In combination with our previous $v\sin i$ measurements in M and K dwarfs, we now derive $P/sin i$ measures for a sample of 418 K and M dwarfs. We investigate the distributions of $P/sin i$ and we show that they are different from one spectral sub-type to another at a 99.9\% confidence level.

Empirically, the X-ray luminosity LX from M dwarfs has been found to have an upper limit of about 0.2% of the bolometric flux Lbol. In the limit where magnetic fields in M dwarfs are generated in equipartition with convective motions, we use stellar models to calculate the energy flux of Alfven waves FA as a function of depth in the sub-surface convection zone. Since Alfven waves have the optimal opportunity for wave modes to reach the corona, we suggest that FA sets an upper limit on the mechanical flux Fmech which causes coronal heating. This suggestion accounts quantitatively for the saturated values of LX/Lbol which have been reported empirically for M dwarfs.

We have measured v sin i with high precision for a sample of dM3 stars (86 targets). We detected rotation in 82 stars (73 dM3 stars, and 9 dM3e stars). We compare our measurements of v sin i for all the stars in our dM0, dM2, dM3 and dM4 samples to those from other authors. We find a good agreement down to v sin i values of less than 1 km/s. The mean of the differences between measurements is only 0.42 km/s. We find that the distribution of P/sini for our dM3 stars is different from the distribution of P/sini among our samples of dM2 and dM4 stars. The mean rotation rate for the dM3 stars (excluding dM3e and sdM3 stars) is significantly slower (25.8 days) than for dM2 (14.4 days) and dM4 stars (11.4 days). Analogous behavior also emerges among the faster rotators (dMe stars): we find that a longer rotation period also occurs at spectral sub-type dM3e. Our data suggest that, as regards the rotational properties of lower main sequence stars, spectral sub-type dM3 stands out as exhibiting unusual slow rotation compared to that of adjoining sub-types. Our data lead us to suggest that the unusual rotational properties of M3 dwarfs may represent a signature of the transition to complete convection (TTCC).

Magnetic fields impede the onset of convection, thereby altering the thermal structure of a convective envelope in a low mass star: this has an effect on the amount of lithium depletion in a magnetized star. In order to quantify this effect, we have applied a magneto-convective model to two low mass stars for which lithium abundances and precise structural parameters are known: YY Gem and CU Cnc. For both stars, we have obtained models which satisfy empirical constraints on the following parameters: R, L, surface magnetic field strength, and Li abundance. In the case of YY Gem, we have obtained a model which satisfies the empirical constraints with an internal magnetic field of several megagauss: such a field strength is within the range of a dynamo where the field energy is in equipartition with rotational energy deep inside the convection zone. However, in the case of CU Cnc, the Li requires an internal magnetic field which is probably too strong for a dynamo origin: we suggest possible alternatives which might account for the reported Li abundance in CU Cnc.

The overstability of the fundamental radial mode in M dwarf models was theoretically predicted by Rodr\'iguez-L\'opez et al. (2012). The periods were found to be in the ranges ~25-40 min and ~4-8 h, depending on stellar age and excitation mechanism. We have extended our initial M dwarf model grid in mass, metallicity, and mixing length parameter. We have also considered models with boundary conditions from PHOENIX NextGen atmospheres to test their influence on the pulsation spectra. We find instability of non-radial modes with radial orders up to k=3, degree l=0-3, including p and g modes, with the period range extending from 20 min up to 11 h. Furthermore, we find theoretical evidence of the potential of M dwarfs as solar-like oscillators.

Precision modeling of M dwarfs has become worthwhile in recent years due to the increasingly precise values of masses and radii which can be obtained from eclipsing binary studies. Torres (2013) has identified 4 prime M dwarf pairs with the most precise empirical determinations of masses and radii. The measured radii are consistently larger than standard stellar models predict. We have previously modeled M dwarfs in the context of a criterion due to Gough & Tayler in which magnetic fields inhibit the onset of convection according to a physics-based prescription. New constraints on the models of M dwarfs are now provided by measurements of lithium abundances. The key aspect of Li in terms of setting constraints on magnetic modeling is that Li burning starts at T = 2.5 MK, and temperatures of just such magnitude are associated with the base of the convection zone: magnetic inhibition of convective onset can shift this base slightly closer to the surface, i.e. to slightly lower temperatures, thereby reducing the amount of Li depletion compared to a non-magnetic model. In the present paper, we consider how our magneto-convection models handle the new test of stellar structure provided by Li measurements. Among the prime systems listed by Torres, we find that plausible magnetic models work well for CM Dra and YY Gem but not for CU Cnc. (The fourth system in Torres's list does not yet have enough information to warrant magnetic modeling.) For CU Cnc, we suggest that the observed lithium may have been accreted from a circumstellar disk. We find that our magneto-convection models of CM Dra, YY Gem and CU Cnc yield results which are consistent with the observed correlation between magnetic flux and X-ray luminosity.

By combining rotation periods with spectroscopic determinations of projected rotation velocity, Jackson, Jeffries & Maxted (2009) have found that the mean radii for low-mass M-dwarfs in the young, open cluster NGC 2516 are larger than model predictions at a given absolute I magnitude or I - K color and also larger than measured radii of magnetically inactive M-dwarfs. The relative radius difference is correlated with magnitude, increasing from a few per cent at MI = 7 to greater than 50 per cent for the lowest luminosity stars in their sample at MI about 9.5. Jackson et al (2009) have suggested that a two-temperature star spot model is capable of explaining the observations, but their model requires spot coverage fractions of at least 50 per cent in rapidly rotating M-dwarfs. Here we examine these results in terms of stellar models that include the inhibiting effects of magnetic fields on convective energy transport, with and without the effects of star spots. We find that a pure spot model is inconsistent with the color - magnitude diagram. The observations of radii versus color and radii versus absolute magnitude in NGC 2516 are consistent with models which include only magnetic inhibition or a combination of magnetic inhibition and spots. At a given mass we find a large dispersion in the strength of the vertical component of the magnetic field in the stellar photosphere but the general trend is that the vertical field increases with decreasing mass from a few hundred Gauss at 0.65 Msun to 600 - 900 Gauss, depending on spot coverage, in the lowest mass stars in the sample at 0.25 Msun.

We use more than 4,500 microflares from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) microflare data set (Christe et al., 2008, Ap. J., 677, 1385) to estimate electron densities and volumetric filling factors of microflare loops using a cooling time analysis. We show that if the filling factor is assumed to be unity, the calculated conductive cooling times are much shorter than the observed flare decay times, which in turn are much shorter than the calculated radiative cooling times. This is likely unphysical, but the contradic- tion can be resolved by assuming the radiative and conductive cooling times are comparable, which is valid when the flare loop temperature is a maximum and when external heating can be ignored. We find that resultant radiative and con- ductive cooling times are comparable to observed decay times, which has been used as an assumption in some previous studies. The inferred electron densities have a mean value of 10^11.6 cm^-3 and filling factors have a mean of 10^-3.7. The filling factors are lower and densities are higher than previous estimates for large flares, but are similar to those found for two microflares by Moore et al. (Ap. J., 526, 505, 1999).

The eclipsing binary CM Dra contains nearly identical red dwarfs of spectral class dM4.5. Their masses and radii have been reported with unprecedentedly small statistical errors. When compared with standard stellar models of appropriate age (\approx4 Gyr), the empirical results indicate that both components are larger in R, and lower in luminosity L, by several standard deviations. Here, we attempt at first to model the components in the context of standard stellar models using differ

The difference in stellar structure above and below spectral type ~M4 is expected to be a very important one, connected directly or indirectly to a variety of observational phenomena in cool stars---such as rotation, activity, magnetic field generation and topology, timescales for evolution of these, and even the basic mass-radius relationship. In this Cool Stars XVI Splinter Session, we aimed to use the M4 transition as an opportunity for discussion about the interiors of low-mass stars and the mechanisms which determine their fundamental properties. By the conclusion of the session, several key points were elucidated. Although M dwarfs exhibit significant changes across the fully convective boundary, this "M4 transition" is not observationally sharp or discrete. Instead, the properties of M dwarfs (i.e. radius, effective temperature, rotation, activity lifetime, magnetic field strength and topology) show smooth changes across M3--M6 spectral types. In addition, a wide range of stellar masses share similar spectral types around the fully convective transition. There appears to be a second transition at M6--M8 spectral types, below which there exists a clear dichotomy of magnetic field topologies. Finally, we used the information and ideas presented in the session to construct a framework for how the structure of an M dwarf star, born with specific mass and chemical composition, responds to the presence of its magnetic field, itself driven by a feedback process that links the star's rotation, interior structure, and field topology.

Although the means of the ages of stars in young groups determined from Li depletion often agree with mean ages determined from Hertzsprung - Russell diagram isochrones, there are often statistically significant differences in the ages of individual stars determined by the two methods. We find that inclusion of the effects of inhibition of convection due to the presence of magnetic fields leads to consistent ages for the individual stars. We illustrate how age consistency arises by applying our results to the \beta Pictoris moving group. We find that, although magnetic inhibition of convection leads to increased ages from the Hertzsprung - Russell diagram isochrones for all stars, Li ages are decreased for fully convective M stars and increased for stars with radiative cores. Our consistent age determination for the \beta Pictoris moving group of 40 Myr is larger than previous determinations by a factor of about two.

Contrary to a common argument that a small increase in the strength of the strong force would lead to destruction of all hydrogen in the big bang due to binding of the diproton and the dineutron with a catastrophic impact on life as we know it, we show that provided the increase in strong force coupling constant is less than about 50% substantial amounts of hydrogen remain. The reason is that an increase in strong force strength leads to tighter binding of the deuteron, permitting nucleosynthesis to occur earlier in the big bang at higher temperature than in the standard big bang. Photodestruction of the less tightly bound diproton and dineutron delays their production to after the bulk of nucleosynthesis is complete. The decay of the diproton can, however, lead to relatively large abundances of deuterium.

There is currently no explanation of why the corona has the temperature and density it has. We present a model which explains how the dynamics of magnetic reconnection regulates the conditions in the corona. A bifurcation in magnetic reconnection at a critical state enforces an upper bound on the coronal temperature for a given density. We present observational evidence from 107 flares in 37 sun-like stars that stellar coronae are near this critical state. The model may be important to self-organized criticality models of the solar corona.

We present Spitzer IRAC and MIPS observations for a sample of eight M dwarfs: six dMe, one dM, and one sdMe star. All of our targets are found to have SEDs which are fitted within the error bars by a purely photospheric spectrum out to 24 micron. We find no evidence for IR excess. None of our targets is detected in the MIPS 70 and 160 micron bands. The estimated ages for all are >10 Myr, suggesting that enough disk dissipation has occurred within the inner several AU of the star. For four of these, Mullan et al. (1989) had reported IRAS detections at 12 micron, although the reported fluxes were below the 5-sigma IRAS detection limit (~0.2 Jy). Mullan et al. also pointed out that V-K colors in dMe stars are larger than those in dM stars, possibly because of the presence of a chromosphere. Here we suggest that metallicity effects provide a better explanation of the V-K data. For reasons of observational selection, our targets are not the most active flare stars known, but being dMe stars indicates the presence of a chromosphere. Scaling from Houdebine's model of the AU Mic chromosphere, we have computed the free-free infrared excesses for a range of densities. Our Spitzer 24 micron data shows that the chromospheres in two of our targets are less dense than in AU Mic by a factor of 10 or more. This is consistent with the fact that our sample includes the less active flare stars. Our models also indicate that the chromospheric contribution to the observed AU Mic emission at submillimeter wavelengths is only about 2%.

The heating of the lower solar corona is examined using numerical simulations and theoretical models of magnetohydrodynamic turbulence in open magnetic regions. A turbulent energy cascade to small length scales perpendicular to the mean magnetic field can be sustained by driving with low-frequency Alfven waves reflected from mean density and magnetic field gradients. This mechanism deposits energy efficiently in the lower corona, and we show that the spatial distribution of the heating is determined by the mean density through the Alfven speed profile. This provides a robust heating mechanism that can explain observed high coronal temperatures and accounts for the significant heating (per unit volume) distribution below two solar radius needed in models of the origin of the solar wind. The obtained heating per unit mass on the other hand is much more extended indicating that the heating on a per particle basis persists throughout all the lower coronal region considered here.

V471 Tauri, an eclipsing system consisting of a hot DA white dwarf (WD) and a dK2 companion in a 12.5-hour orbit, is the prototype of the pre-cataclysmic binaries. The late-type component is magnetically active, due to its being constrained to rotate synchronously with the short orbital period. During a program of UV spectroscopy of V471 Tau, carried out with the Goddard High Resolution Spectrograph (GHRS) onboard the Hubble Space Telescope, we serendipitously detected two episodes in which transient absorptions in the Si III 1206 A resonance line appeared suddenly, on a timescale of <2 min. The observations were taken in a narrow spectral region around Ly-alpha, and were all obtained near the two quadratures of the binary orbit, i.e., at maximum projected separation (~3.3 Rsun) of the WD and K star. We suggest that these transient features arise when coronal mass ejections (CME's) from the K2 dwarf pass across the line of sight to the WD. Estimates of the velocities, densities, and masses of the events in V471 Tau are generally consistent with the properties of solar CME's. Given our detection of 2 events during 6.8 hr of GHRS observing, along with a consideration of the restricted range of latitudes and longitudes on the K star's surface that can give rise to trajectories passing in front of the WD as seen from Earth, we estimate that the active V471 Tau dK star emits some 100-500 CME's per day, as compared to 1-3 per day for the Sun. The K dwarf's mass-loss rate associated with CME's is at least (5-25) x 10^{-14} Msun/yr, but it may well be orders of magnitude higher if most of the silicon is in ionization states other than Si III.