# The Active Corona of HD 35850 (F8 V)

## 1 Introduction

By analogy to the Sun, magnetic heating in the chromospheres, transition regions, and coronae of late-type stars is observed as cooling in the Balmer lines, the Ca II lines, and in UV, FUV, EUV, and soft X-ray continuum and line emission. Results from the Mt. Wilson Ca II H and K survey of single lower main-sequence field stars has provided compelling evidence of Solar-like chromospheric active regions, activity cycles, differential rotation, and occasional Maunder minima in mid-F to early-M stars (e.g., Noyes et al. 1984; Baliunas et al. 1996; Gray & Baliunas 1997; Baliunas et al. 1995).

Coronal heating and magnetic dynamos on the Sun and in stars have been discussed extensively in the literature (cf. Haisch & Schmitt 1996). The various flares, microflares, nanoflares, and brightenings observed on the Sun are thought to be manifestations of current sheet reconnection and much of the coronal heating debate is focussed on whether or not these transient events can account for the radiative output of the solar corona (e.g., Oreshina & Somov 1998). While both continuous (e.g., MHD waves) and stochastic (nanoflaring) processes probably contribute to magnetic coronal heating, the dominant process has yet to be identified.

Empirically, coronal emission on active, cool, main-sequence stars is similar to flare emission from solar active regions. For example, Benz & Güdel (1994) find a similar correlation between X-ray luminosity and non-thermal continuum radio luminosity for solar flares and active stars. On the Sun, mildly relativistic electrons accelerated along coronal magnetic field lines produce gyrosynchrotron radio emission. The presence of persistent but variable nonthermal radio emission from cool stars indicates a continual replenishment of the relativistic electron population while a correlation between coronal and radio emissions suggests a causal relationship between electron acceleration in magnetically confined loops and coronal heating. In the 50,000–200,000 K regime, further evidence of continual flaring comes from the broad emission-line components seen in high-resolution UV spectra of active dwarfs and RS CVn binaries. Because the broad-component line widths are 2–4 times the thermal width and because the broad components are blue-shifted relative to the narrow components, Wood et al. (1996) suggest that the broad components are produced by transition-region explosive events.

Two well-studied examples of extreme main-sequence magnetic activity are the Pleiades-age, rapid rotators AB Doradus (K0 V) and EK Draconis (G0 V). Spectro-polarimetric monitoring of AB Dor shows long-lived, cool magnetic spots and latitudinal differential rotation (Donati & Collier Cameron 1997). ROSAT PSPC soft X-ray photometry of AB Dor (Kürster et al. 1997) shows rotationally modulated flare and quiescent emission with the same phase and period (12.4 h) as the photospheric spots. ASCA and EUVE spectra of AB Dor and EK Dra indicate coronal emission-measure distributions, $`\mathrm{Em}(T)`$, with peaks at 5–8 MK and 20–30 MK (Mewe et al. 1996; Güdel et al. 1997).

The 5–8 MK component has been observed in other moderately active coronal sources, most notably Capella (Brickhouse, Raymond, & Smith 1995). As has been pointed out by Gehrels & Williams (1993), optically thin plasma tends to accumulate at temperatures where the cooling curve $`\mathrm{\Omega }(T)`$ has a positive slope. This is a manifestation of the Parker (1953) & Field (1965) instability: when $`\frac{d\mathrm{\Omega }(T)}{dT}>0`$, plasma cools more efficiently at slightly higher temperatures and less efficiently at slightly lower temperatures. At temperatures where $`\frac{d\mathrm{\Omega }(T)}{dT}<0`$, plasma quickly cools to lower temperatures. For solar-abundance plasmas, $`\frac{d\mathrm{\Omega }(T)}{dT}>0`$ in the ranges 0.7–1.2 MK, 4.5–7 MK, and above 21 MK. Thus, the shape of the cooling curve explains the dip from 2–4 MK and from 9–20 MK seen in many EM distributions and provides some physical justification for two- and three-temperature coronal models. The temperature of the hot component and the relative amount of hot emission measure must reflect a balance between magnetic heating and radiative cooling. For stars with significant emission measure above 20 MK, some mechanism must be super-heating the coronal plasma at fairly regular intervals.

For three active, solar-mass stars observed with the ASCA SIS (EK Dra, HN Pegasii, and $`\kappa ^1`$ Ceti in order of decreasing activity), Güdel (1997) has modeled $`\mathrm{Em}(T)`$ by assuming X-rays are produced by an ensemble of flaring and cooling loops. The Güdel (1997) model is a modification of the solar nanoflare model proposed by Kopp & Poletto (1993). While the ASCA data cannot constrain many model parameters (number of loops, loop dimension, and mean magnetic field strength), the models can only reproduce the observed $`\mathrm{Em}(T)`$ provided that the power-law distribution of flare and microflare energies is $`\alpha 2`$. High S/N spectra of other active dwarfs are needed to test the viability of continuous flaring and to establish the importance of this mechanism as a function of rotation rate and effective temperature.

In this paper, we present EUV and X-ray spectroscopy of HD 35850 = HR 1817 (F8 V). HD 35850 was detected serendipitously with EXOSAT (Cutispoto et al. 1991) and follow-up optical spectroscopy has established it as a nearby, single, solar-metallicity Pleiades-age ($`t10^8`$ yr), rapid rotator (Tagliaferri et al. 1994). We are particularly interested in HD 35850 because it is probably in the extreme state of magnetic activity for single, main-sequence F stars. There is considerable interest in probing the magnetic dynamo in F stars because they possess relatively shallow convection zones. Kim & Demarque estimate that, in a $`10^7`$ yr-old $`1.1`$$`1.2M_{}`$ star like HD 35850, the convective turnover time is $`\tau _\mathrm{c}20`$ d. In a young $`1.0M_{}`$ star like EK Dra, $`\tau _\mathrm{c}40`$ d.

## 2 Observations and Data Analysis

HD 35850 was observed by the Extreme Ultraviolet Explorer from 1995 October 23 08:17:07 UT to 1995 October 30 07:20:34 UT. The EUVE Deep Survey/Spectrometer consists of three aligned grazing-incidence telescopes. The telescope beams are intercepted by short-, medium-, and long-wavelength reflection gratings and detected with microchannel-plate detectors (SW, MW, and LW, respectively). The undispersed portions of the three beams are focussed onto the Deep Survey microchannel-plate detector, DS. This way, EUVE obtains simultaneous, time-resolved 75–120 Å broad-band photometry and 70–160 Å, 170–370 Å, and 300–525 Å medium-resolution spectra (cf. Haisch, Bowyer, & Malina 1993). In support of this EUVE observing program, Ca II HK spectra were obtained on 1996 October 2–12. HD 35850 was observed with the ASCA SIS and GIS for $`18`$ ks on 1995 March 12. The reduction and analysis of these data are described below.

### 2.1 EUVE Deep Survey Data

The EUVE Deep Survey/Spectrometer records photon events whenever the source is visible by the satellite, including times when the satellite is passing through the South Atlantic Anomaly (SAA). These SAA passages are characterized by increased background levels which lead to significant and poorly characterized loss of telemetry, nicknamed “primbsching”. Correcting for earth occultations, the on-source time was $`203`$ ks. Correcting for instrument and telemetry dead times, the net exposure time was $`198`$ ks.

The Deep Survey (DS) events were screened at various count-rate and primbsch-correction thresholds. Light curves were generated with bin times ranging from 100 to 2000 s. Source events were extracted from a $`25`$ pixel radius circle centered on the PSF centroid. Background events were extracted from an annulus with inner and outer radii of $`37`$ and $`67`$ pixels, respectively. In Figure 1, we show the background-subtracted, dead-time and primbsch-corrected 75–120 Å DS light curve for HD 35850 using 1000 s bins and including times when primbsching (monitor Det7Q1dpc) was below 25%.

The mean DS background-corrected count rate of HD 35850 is 0.20 counts s<sup>-1</sup>. No large flares are evident in the DS light curve: the variability amplitude is approximately 100% and the most significant deviation from the mean count rate is a $`5\sigma `$ peak near MJD 18.45. The light curve, however, shows small to moderate flares throughout the observation. At least 7 events have a peak count rate $`3\sigma `$ above the apparent “quiescent” level of 0.20 counts s<sup>-1</sup>. To test the hypothesis that the observed counts come from a non-variable source, we have determined the cumulative distribution function of the unbinned photon arrival times for the observed data and for simulated data from a constant, 0.2 counts s<sup>-1</sup> source. A two-sample Kolmogorov-Smirnov test statistic (Feigelson & Babu 1992) was calculated using the observed and simulated distributions; the probability that the two data sets were drawn from the same parent population is low ($`P10^9`$). We conclude that the DS light curve exhibits significant variability.

We have used the binned DS light curve to compute the discrete Fourier transform (Scargle 1989), shown in Figure 2. There is significant power around 0.98 d and its aliases. For HD 35850, $`v\mathrm{sin}i50`$ km s<sup>-1</sup> and $`R_{}1.18R_{}`$ implies $`P_{\mathrm{rot}}/\mathrm{sin}i1.1`$ d, consistent with the 0.98-d period. However, Halpern & Marshall (1996) report that 0.98 d is a beat period associated with the satellite’s passage through the SAA. For example, during the HD 35850 observation, the Det7Q1dpc monitor shows strong periodic signals at 0.984 d and at 0.0656 d. Similar periods have been seen in other primbsch-corrected light curves (Marshall 1998, private communication). We thus conclude that some or all of the 0.98-d signal seen in the DS light curve results from SAA passages. We note that the only period seen in HD 35850’s power spectrum that cannot be attributed to primbsching is at 1.40 d. Until further observations can be carried out, we tentatively identify 1.40 d as the rotation period of HD 35850 (see §2.2).

### 2.2 McMath-Pierce Optical Spectroscopy

Medium-resolution Ca II HK spectra were obtained with the Solar-Stellar Spectrograph at the National Solar Observatory’s McMath-Pierce Telescope on 1996 October 2–12, ten days before the EUVE observation. The spectra around HK are typical of active, late-type dwarfs, with broad absorption troughs and bright, narrow emission cores. To calculate the equivalent width of the H and K emission, we measured the excess emission in the line cores, normalizing by the integrated continuum intensity in a 1 Å box far from line center.

In Figure 3 we plot the equivalent width of the Ca II H (dots) and K (asterisks) emission cores. While the equivalent widths are highly variable on short time scales, HD 35850 was only visible for a few hours each night. The limited sampling makes it difficult to estimate $`P_{\mathrm{rot}}`$ with these data alone. In Figure 4, we have phase-folded the summed H and K equivalent widths using $`P_{\mathrm{rot}}=1.40`$ d as derived from the DS photometry.

### 2.3 EUVE SW and MW Spectrometer Data

For the SW, MW, and LW spectroscopic data, the photon event data were screened to eliminate times of high background using the IRAF EUV package. Specifically, data were excluded when the total detector count rate exceeded 40 counts s<sup>-1</sup>. The SW, MW. and LW photon event lists were used to create wavelength-calibrated FITS images and the SW, MW, and LW images were then analyzed within IDL.

We note that during the HD 35850 observation, the SW microchannel-plate detector developed a hot spot 8–12 detector pixels from the dispersion axis, close to the location of the Fe XXII 114 Å line. A light curve of the hot spot intensity was used to identify the times when the hot spot flickered; those times were removed from the photon event list. Some hot pixels were still visible in the SW image, so a $`9\times 9`$ pixel box around the hot spot was set to the local background level.

Spectra have been extracted from the SW, MW, and LW images using a modified version of the IUE SIPS spectral extraction algorithm (Lenz & Ayres 1992). The optimal extraction routine is specifically tailored to low S/N spectra from photon-counting devices like the EUVE and IUE spectrometers. Briefly, a raw SW, MW, or LW $`2048\times 2048`$ pixel image is trimmed and compressed to $`1024\times 290`$ and 2-pixel Gaussian smoothed. Since the spectra are over-sampled (no less than 7 detector pixels per resolution element), 2-pixel compression and smoothing improves S/N per bin without compromising spectral resolution. Using a predefined background region away from the source spectrum chosen to reduce curved edge effects caused by the EUVE reflection gratings, a background image is generated. A background-subtracted source image and corresponding error image are then created. The source image is used measure the cross-dispersion profile. To perform an optimal extraction, each row in the source image is weighted by the cross-dispersion profile and the source spectrum and corresponding error are extracted using a 7-pixel aperture.

No emission lines longward of 365 Å were detected in the LW count spectrum, presumably because of increased interstellar absorption at long wavelengths. Since the lines shortward of 365 Å are measured in the MW spectrum, we did not consider the LW spectrum in our analysis. In Figures 5 and 6 (bottom panels), we show the flux-calibrated SW and MW spectra.

#### 2.3.1 IDL Line and Continuum Fitting

In order to measure line fluxes, we have constructed a line list using the emissivity lists of Brickhouse et al. (1995), Monsigniori Fossi & Landini (1993), and Mewe et al. (1985), in order of priority. Lines were grouped to match the SW and MW spectral resolution (0.5 and 1 Å, respectively), allowing us to account for emissivity from blended lines. Based on the ionization states of the brightest identified lines, e.g., Fe XV, Fe XVI, and Fe XX, we made a starting guess at the plasma temperature ($`\mathrm{log}T6.8`$) and assumed a density $`n_\mathrm{e}=10^{11}`$ cm<sup>-3</sup>. The line emissivities, spectrometer effective area curves, and exposure time were used to provide a starting guess for each line strength. The observed spectrum was then fit with multiple Gaussians, varying wavelength ($`\pm `$ one resolution element) and line strength. Lines with fewer than $`2\sigma `$ counts were eliminated from the fit and the procedure was repeated. This way, line counts and upper limits were determined for all high-emissivity Fe IX to Fe XXIV lines.

The SW spectrum shows a small but detectable thermal bremsstrahlung continuum. To determine the SW continuum luminosity from 85–135 Å, the continuum bins were fit with two parameters (emission measure and temperature) using the continuum emissivity model of Mewe et al. (1985). The Hipparcos distance to HD 35850 is 26.8 pc (Perryman et al. 1997). For HD 35850, the 85–135 Å continuum luminosity is $`L_{\mathrm{cont}}4.6\pm 0.4\times 10^{28}`$ ergs s<sup>-1</sup>.

In Figs. 5 and 6 (top and middle panels), we show the observed and fit count spectra with line identifications. Detected lines are indicated with dotted lines and upper limits are indicated with dashed lines. In Tables 1 and 2, we list the ion, possible blends, temperature of peak emissivity, laboratory and measured wavelengths, a detection flag, measured line luminosities, and the signal-to-noise ratio for detected lines. The line luminosities in Tables 1 and 2 use $`N_\mathrm{H}=1.7\times 10^{18}`$ cm<sup>-2</sup> (see §2.3.3). In all, 28 distinct lines are detected above $`2\sigma `$ in the SW and MW spectra.

#### 2.3.2 Electron-Density and Column-Density Estimates

The Fe XV and Fe XVI $`\lambda \lambda 285,335,365`$ line ratios are relatively insensitive to density and provide the best estimate of interstellar $`N_\mathrm{H}`$ in the EUVE bandpass. We find $`N_\mathrm{H}=1.9\pm 0.4\times 10^{18}`$ cm<sup>-2</sup>.

Two density-sensitive line ratios can be measured accurately from the SW spectrum of HD 35850: Fe XXII $`\lambda \lambda 114,117`$ and Fe XXI $`\lambda \lambda 102,129`$. In Figure 7 we plot the observed SW spectrum (histogram) and fitted lines and continuum (solid line). The density-sensitive lines are indicated with dashed lines. The line Fe XXI line ratio is $`0.22\pm 0.05`$. The predicted branching ratio (solid curve) and observed line ratio (cross with 1$`\sigma `$ error bar) are plotted versus density in Figure 8. Fig. 8 suggests that $`\mathrm{log}n_\mathrm{e}<11.6`$ cm<sup>-3</sup>. Consequently, we have chosen $`\mathrm{log}n_\mathrm{e}=11.0`$ cm<sup>-3</sup> for the IDL and SPEX emission-measure analyses.

The Fe XXII line ratio is $`0.47\pm 0.14`$, suggesting $`n_\mathrm{e}>10^{13}`$ cm<sup>-3</sup>. This ratio is unusually high and the implied density does not agree with the more reliable Fe XXI measurement. The detector hot spot (§ 2.3) may have contributed to the anomalous 114 Å flux. We note that line ratio calculations may have systematic uncertainties of up to 50% (Brickhouse et al. 1995) and that the derived density limit is only approximate.

#### 2.3.3 IDL Emission-Measure Analysis

The measured line luminosities and upper limits have been used to estimate the coronal emission-measure distribution, $`\mathrm{Em}(T)`$. Previous analyses of EUVE spectra (e.g., Güdel et al. 1997; Mewe et al. 1996) have generally relied on global fits to estimate $`\mathrm{Em}(T)`$. For low S/N spectra like these, we prefer to use the bright lines in Tables 1 and 2 because (i) the emissivities of brighter lines have been fairly well established and (ii) the $`\chi ^2`$ statistic is not dominated by continuum bins.

We have excluded the He II $`\lambda 304`$ line from our analysis. On the Sun, the He II Ly$`\alpha `$ line is not a reliable temperature indicator: some, perhaps most, of the He II emission occurs in the chromosphere as a result of back-heating from the corona. We note also that around $`\lambda 255`$, the Fe XXIV line is blended with He II. The combined flux from both lines is listed as an upper limit in Table 2. For the remaining lines, our fitting procedure minimizes $`\chi ^2`$ between the measured and predicted line luminosities by varying the amplitude and shape of $`\mathrm{Em}(\mathrm{log}T)`$ (in $`\mathrm{\Delta }\mathrm{log}T=0.1`$ bins).

We use exponential Chebyshev polynomials (Lemen et al. 1989) to describe $`\mathrm{Em}(\mathrm{log}T)`$. In addition to detected emission lines, the undetected Fe X to Fe XIV lines have been used as additional constraints in the non-linear least-squares fitting. The SW continuum luminosity serves to constrain Fe abundance.

Best-fit emission-measure distributions have been derived for a grid of column densities, Fe abundances, and Chebyshev-polynomial orders: $`0.0N_\mathrm{H}4.0\times 10^{18}`$ cm<sup>-2</sup>, $`0.1Z2.0`$, and $`2n11`$, where $`Z`$ is the coronal Fe abundance relative to the solar photospheric value of Anders & Grevesse (1989). We obtain acceptable $`\chi _\nu ^2`$ values for all $`n6`$ and the resulting $`\mathrm{Em}(T)`$ are characterized by peaks at $`\mathrm{log}T6.8`$ and $`\mathrm{log}T7.4`$. The best-fit column densities, abundances, and emission-measure distributions are not a strong function of $`n`$. As a result, we choose $`n=6`$, the lowest polynomial order for which we obtain good fits.

Fits using different order polynomials orders, column densities, and abundances yield qualitatively similar results: prominent peaks in the emission-measure distribution at $`\mathrm{log}T6.8`$ and $`\mathrm{log}T7.4`$, although the sharpness of the high-T peak is not well constrained. The relative strength of the high-T component increases with increasing abundance. The best-fit emission-measure distributions all go to zero at very high temperatures ($`\mathrm{log}T>7.7`$).

In Figure 9, we plot $`\mathrm{log}\mathrm{\Delta }\chi ^2`$ as a function of $`N_\mathrm{H}`$ and $`Z`$. The best fit ($`\chi _\nu ^21.07`$) is obtained for $`N_\mathrm{H}=1.7\times 10^{18}`$ cm<sup>-2</sup> and $`Z=1.15`$. The 68% and 90% confidence contours are plotted for 9 free parameters ($`Z`$, $`N_\mathrm{H}`$, and $`n=0,\mathrm{},6`$). The EUVE line-to-continuum ratio yields $`Z=1.15_{0.35}^{+0.75}`$. To obtain substantially sub-solar abundances ($`Z<0.5`$), would require more than doubling the continuum level (dark solid line in Fig. 7). This is not compatible with the background-subtracted SW spectrum. We have experimented with various background-subtraction algorithms (smoothed, unsmoothed, polynomial fit): all produce bright lines and a weak continuum.

We find adequate fits for a range of column densities and abundances, and, in particular, we find $`\chi _\nu ^21.15`$ for $`N_\mathrm{H}=1.4\times 10^{18}`$ cm<sup>-2</sup> and $`Z=1`$. We adopt these values in Figs. 10–14. In Figure 10, the observed and predicted line luminosities for detected (filled circles) and undetected (open triangles) lines are compared, with a factor of two deviation indicated with dashed lines. Note the anomalously high Fe XXII line flux (see §2.3.2). In Figure 11, we plot $`\mathrm{Em}(\mathrm{log}T)`$ per $`\mathrm{log}T`$ bin of 0.1 dex as derived from the IDL line analysis (dash-dotted line).

#### 2.3.4 SPEX Differential Emission-Measure Analysis

In order to verify the results of the IDL line analysis, we performed a differential emission measure analysis of the optimally-extracted EUVE SW and MW spectra using the SPEX code (see Kaastra et al. 1992 for a detailed description of the DEM method and the SPEX Collisional Ionization Equilibrium model). For fixed values of $`N_\mathrm{H}`$ and $`Z`$, we derived $`\mathrm{Em}(T)`$ using the SPEX regularization method. In Fig. 11, we plot $`\mathrm{Em}(T)`$ for $`N_\mathrm{H}=1.4\times 10^{18}`$ cm<sup>-2</sup> and $`Z=1`$ (solid line). The IDL and SPEX analyses yield similar emission-measure distributions, with 55–60% of the emission measure in the hotter component above $`\mathrm{log}T=7`$.

### 2.4 ASCA SIS data

An 18-ks ASCA exposure of HD 35850 obtained in 1995 has been analyzed by Tagliaferri et al. (1997). We use their reduced ASCA SIS0 spectrum to further constrain the EUVE results. Because the SIS cannot resolve individual emission lines from coronal sources, we cannot use the IDL line analysis method described in § 2.1.4. Consequently, we fitted the ASCA SIS0 spectrum using the SPEX DEM code. For the SIS0 fits, $`N_\mathrm{H}`$ was fixed at $`1.4\times 10^{18}`$ cm<sup>-2</sup>. The best-fit abundance using the DEM collisional ionization equilibrium (CIE) model is $`Z0.5`$ with acceptable fits in the range 0.34–0.81. The ASCA upper bound on $`Z`$ is plotted in Fig. 9 as a dashed line. We note that $`Z0.8`$ is marginally compatible with EUVE and ASCA.

To illustrate, the ASCA SIS0 data are shown in Figure 12 (points) with the DEM model spectrum (histogram) for $`Z=1`$. Discrepancies between the SIS data and the $`Z=1`$ model are seen around 1.2 keV and 2.4 keV. These are the spectral features which normally drive the fit towards lower abundances. Brickhouse et al. (1997) have identified a complex of lines from 1.0–1.3 keV from highly excited ($`n>5`$), highly ionized (Fe XVII to Fe XXV) states that are missing from the plasma codes in XSPEC and SPEX. The missing lines may partly explain the large discrepancy around 1.2 keV.

In Figure 13, we show the corresponding ASCA SPEX $`\mathrm{Em}(\mathrm{log}T)`$ for $`Z=1`$. Like the EUVE emission-measure distributions, the ASCA $`\mathrm{Em}(T)`$ peaks at $`\mathrm{log}T`$ of 6.7 and 7.4. Finally, we note that the EUVE count rates are higher than expected from the SIS spectrum. Based on the SIS emission-measure analysis, the expected EUVE DS count rate is approximately 0.16 counts s<sup>-1</sup>, indicated in Fig. 1 as a dashed line. HD 35850 appears to have been more active (by $`25\%`$) during the EUVE observation.

### 2.5 Comparison with Previous Results

For HD 35850, we find that the EUVE line-to-continuum ratio indicates approximately solar coronal Fe abundance while a SPEX DEM analysis of the ASCA SIS0 spectrum is consistent with moderately sub-solar abundances. We note that HD 35850’s photospheric Fe abundance is close to solar (Tagliaferri et al. 1994). Analyses of other active coronal sources often find coronal abundances far below the measured photospheric abundances (e.g., AB Dor, Mewe et al. 1996).

The ASCA SIS and GIS spectra of HD 35850 were fit by Tagliaferri et al. (1997) using a number of multi-temperature plasma models in XSPEC (Arnaud 1996). For example, the best-fit SIS0 MEKAL parameters indicate a two-temperature corona ($`kT_10.6`$ keV and $`kT_21.2`$ keV) with sub-solar abundances ($`0.12<Z<0.25`$). Tagliaferri et al. (1997) also performed a combined PSPC/SIS/GIS analysis using a 5-ks ROSAT PSPC observation of a field near HD 35850. The combined analysis involved comparing 2-T and 3-T MEKAL and RS models in XSPEC with solar and non-solar abundances and using different cross-detector normalization constraints (see their Table 5 and Figure 3). The XSPEC fits generally favor sub-solar abundance models; the best-fit MEKAL 3-T model yields $`Z=0.34\pm 0.04`$, $`T_1=0.52\pm 0.07`$, $`T_2=0.78_{0.09}^{+0.15}`$, and $`T_3=1.9_{0.4}^{+1.0}`$, with $`47\%`$ of the total emission measure in the coolest component. This cool component is the $`\mathrm{log}T=6.8`$ component seen in the EUVE DEM. The two hottest components appear to represent the hotter DEM component. Our DEM analysis of the ASCA SIS spectra yields somewhat higher abundances ($`0.34<Z<0.81`$) than Tagliaferri et al. (1997) and somewhat lower Fe abundance than our EUVE analysis.

Systematic differences between ASCA and EUVE analyses have been reported for other bright coronal sources. For example, Brickhouse et al. (1997) obtained simultaneous observations of Capella with ASCA and EUVE in 1996 March. As with previously obtained EUVE spectra of Capella, they find essentially solar photospheric Fe abundance. On the other hand, the plasma codes in XSPEC yield consistently poorer fits to the ASCA SIS spectra of Capella and require sub-solar abundances ($`Z0.7`$).

Note: After the submission of this paper, a preprint of a paper to appear in Astronomy and Astrophysics by Mathioudakis & Mullan (1998) came to our attention which analyzes the EUVE observations of HR 1817=HD 35850. This paper and the Mathioudakis & Mullan paper arrive at substantially different conclusions regarding the Fe abundance of HD 35850. Mathioudakis & Mullan (1998) visually compare the observed EUVE SW and MW spectra with synthetic spectra based on the 3-T MEKAL model suggested by Tagliaferri et al. (1997). They conclude that the EUVE spectra are consistent with Tagliaferri et al. (1997) and are inconsistent with any solar-abundance model, particularly because the Fe XIII–Fe XV lines in the MW spectrum are so weak. While it is difficult to visually compare spectra, their observed and synthetic MW spectra appear to be quite different: the predicted lines are weak and the the low-Z model appears to over-predict the MW continuum (see Figures 2 and 5 in Mathioudakis & Mullan 1998). Given the limitations of ASCA and EUVE, higher-resolution, higher-S/N AXAF HETG or XMM RGS spectra may be needed to determine HD 35850’s coronal abundance.

## 3 Microflaring

The double-peaked EM distribution observed on three active, solar-type stars has been modeled by Güdel (1997) on the basis of a simplified stochastic flare model derived from the solar nanoflare model of Kopp & Poletto (1993). Rather than approaching the EM distribution problem hydrostatically, as is the case when using loop scaling laws, Kopp & Poletto (1993) use a simplified hydrodynamic model to treat a large number of flares.

The flare model reduces the loop hydrodynamics to a point model with a chromospheric energy sink. A series of energy pulses of finite duration is fed into the loop. This heating energy is lost by conduction into the chromosphere and by radiation into space. Without flaring, all loops are kept at the same equilibrium temperature dictated by hydrostatic loop scaling laws (Rosner, Tucker, & Vaiana 1978). Therefore, our model starts out at a lower threshold temperature and does not consider cooler loops. The salient feature of this model is its phenomenological similarity with more sophisticated hydrodynamic simulations in terms of the emission measure, temperature, and radiation history. The point model cannot treat the more complicated 1-D problem of radiation and conduction occurring between the corona, the transition region, and the chromosphere. The need to simulate a large number of flares requires such a simplified approach.

Flares are ignited in each loop randomly distributed in time but satisfying a statistical number distribution in total flare energy $`E`$ above a pre-defined threshold energy. I.e., $`\mathrm{d}N/\mathrm{d}EE^\alpha `$ where $`N(E)`$ is the number density of flares in the energy interval $`[E,E+\mathrm{d}E]`$. On the basis of solar and stellar optical flare monitoring, $`1.8<\alpha <2.0`$ has been suggested by Hudson (1991). For consistency with Güdel (1997), we chose a threshold energy $`E>10^{27}`$ ergs. In order to reduce the number of free parameters, all loops have the same height and vary only in their thickness between a minimum of $`6.4\times 10^6`$ cm and a maximum of $`4.3\times 10^9`$ cm to accommodate the range of flare energies. Each run produces a large number of small flares and only a small number of very large flares. Because no very large flares were observed by EUVE, the largest flares in the simulations were eliminated.

The total EM distribution is determined by averaging (in time and space) the EM distribution of all loops. Similarly, the light curve is the superposition of coronal radiation leaving the star into $`2\pi `$ steradians, i.e., half the radiative losses are absorbed in the chromosphere, half are radiated into space. We have experimented with various parameters and find that our model reproduces the required coronal luminosity and EM distribution only if we assume compact loops, with a semi-length of $`2.6\times 10^9`$ cm and essentially full surface coverage. The simulations provide fairly tight constraints on the power-lax index of flare energies ($`\alpha 1.8`$). In this case, the equilibrium temperature is $`6\times 10^6`$ K.

Our best-fit microflare model for HD 35850 uses 1100 loops, distributed evenly over the star’s surface, and accommodating $`2.64\times 10^7`$ flares during a typical observation lasting $`5\times 10^5`$ s. The resulting EM distribution is shown in Figure 14. The soft X-ray light curves from these simulations show the same 1–12 hour variability seen in ASCA SIS and EUVE DS light curves. Because the power-law distribution of flares favors small (microflare) events, the simulated light curves also show lower-amplitude, shorter-time scale variability seen in higher-S/N X-ray observations of dMe stars but not discernible for HD 35850.

The model EM distribution exhibits the characteristic two-temperature structure with a broad intermediate minimum. However, in this model the EM structure is not just a result of radiative cooling as suggested by Gehrels & Williams (1993). It is a consequence of the balance between flare heating, cooling by conduction, and cooling by radiation. Given our simplistic model, it is noteworthy that the two-temperature structure (or rather the broad minimum around 10 MK) observed on many active normal stars can be explained, in part, by a hydrodynamic property of frequently-flaring loops.

## 4 The Active Corona of HD 35850

We summarize the combined EUVE and ASCA results for HD 35850 in Table 3. The ASCA SIS spectral analysis suggests sub-solar abundances and the EUVE line-to-continuum ratio indicates approximately solar photospheric abundances. Although $`Z0.8`$ is marginally consistent with the EUVE and ASCA spectra, all the ASCA analyses yield systematically lower Fe abundance than the EUVE Fe line-to-continuum ratio.

Assuming photospheric abundances, the corona is characterized by emission measure in a warm component from 5–8 MK and a hot component from 21–30 MK. Like other rapidly rotating, Pleiades-age dwarfs, HD 35850 appears to represent an activity extremum for main-sequence solar-type stars. HD 35850’s inverse Rossby number is $`\tau _\mathrm{c}/P_{\mathrm{rot}}20.7/1.4015`$ and its X-ray surface flux is $`F_\mathrm{X}1.8\times 10^7`$ ergs s<sup>-1</sup> cm<sup>-2</sup>. For AB Dor (K0 V), $`F_\mathrm{X}1.7\times 10^7`$ ergs s<sup>-1</sup> cm<sup>-2</sup> (Hempelmann et al. 1995) and $`\tau _\mathrm{c}/P_{\mathrm{rot}}>120`$. For EK Dra (G0 V), $`F_\mathrm{X}1.5\times 10^7`$ ergs s<sup>-1</sup> cm<sup>-2</sup> and $`\tau _\mathrm{c}/P_{\mathrm{rot}}15`$. AB Dor, HD 35850, and EK Dra appear to have saturated or nearly saturated X-ray activity and, despite having more rapid rotation and a deeper convection zone, AB Dor’s X-ray surface flux is the same as HD 35850’s.

The EUVE DS light curve of HD 35850 shows about one moderate-amplitude (40% increase) flare per day. The X-ray and EUV variability and the presence of substantial emission measure above 20 MK suggests that some flare-like mechanism must be heating the corona. To test this hypothesis, we have modeled HD 35850’s light curves and $`\mathrm{Em}(T)`$ using the hydrodynamic microflare point model of Güdel (1997). The simulations suggest a power-lax index of flare energies $`\alpha 1.8`$.

On the Sun, transient brightenings seen in Yohkoh Soft X-ray Telescope images of solar active regions provide a direct measure of the microflaring amplitude and power-law index. Microflaring can provide, at most, 20% of the heating rate required to power the active-region corona (Shimizu 1995). Ofman, Davila, & Shimizu (1996) have shown that the transient brightenings seen by Yohkoh may be a consequence of resonant absorption of global-mode Alfvén waves in coronal loops, excited by random footpoint motions of these loops. A more recent analysis of Yohkoh data suggests that a steeper distribution of smaller flares, the so-called nanoflares, may also be present (Shimizu & Tsuneta 1997). Whether the microflares occur as a result of Alfvén waves or magnetic reconnection, they cannot account for the radiative output of the Solar corona. For the rapidly rotating F and G dwarfs, however, the microflare models are able to reproduce the observed luminosity, EM distribution, and variability.

To better test microflaring on stars, accurate coronal density and abundance measurements as a function of temperature are needed to better constrain the emission-measure distribution. Photospheric magnetic field strength and filling fractions can further constrain loop models. Also, higher S/N, higher cadence hard and soft X-ray light curves will provide a statistically robust estimate of the distribution of flares. Planned observations of AB Dor with the Advanced X-ray Astrophysics Facility’s High-Energy Transmission Grating Spectrometer will help resolve some of the density and abundance issues discussed in this paper and may lead to more quantitative tests of coronal heating models.

MG would like to thank Tom Ayres for providing some of the IDL routines to perform optimal extraction and line fitting. MG made extensive use of the NASA Astrophysics Data System abstract service, the HEASARC database at NASA/GSFC, and the Simbad database at the Centre de Données astronomiques de Strasbourg. The authors would like to thank an anonymous referee for many helpful suggestions. This research was supported under NASA grant NAG-2891 to the University of Colorado.