Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.sai.msu.ru/dept/adm/lamzin/aa2787.pdf Äàòà èçìåíåíèÿ: Fri Sep 28 16:47:17 2012 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 23:03:38 2012 Êîäèðîâêà:

A&A 377, 557­565 (2001) DOI: 10.1051/0004-6361:20011108
c ESO 2001

Astronomy & Astrophysics

UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau
L. Errico1 , S. A. Lamzin2 , and A. A. Vittone1
1 2

Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy Sternberg Astronomical Institute, Moscow V-234, 119899, Russia

Received 9 April 2001 / Accepted 19 June 2001 Abstract. Ultraviolet spectra of BP Tau observed with HST/GHRS and IUE satellites were analysed. We found that BP Tau activity can be explained in the frame of a disk accretion model if we assume that the stellar magnetic axis is strongly inclined to the disk plane. The following set of accretion process parameters were derived: relative surface area of the accretion zone f 0.25, accretion rate Mac 3.6 â 10-8 M yr-1 , accretion energy flux 10 -1 -2 F = 2.3 â 10 erg s cm and accretion luminosity Lac 0.37 L . The relevance of these parameters is discussed. We argue that the Calvet & Gullbring (1998) accretion shock model is too crude to believe that the accretion spot surface area is indeed proportional to the square of the accretion rate, as Ardila & Basri (2000) found through this model. A strong flare in the Fe ii 2811.8, 2812.1 ° lines was detected, it was probably produced A by an increase of the accretion rate. During the flare, the accretion luminosity was comparable to or even larger than the stellar bolometric luminosity. Key words. stars: pre-main sequence ­ stars: individual: BP Tau ­ ultraviolet: stars ­ X-rays: stars

1. Intro duction
BP Tau is a classical T Tauri star (CTTS): its H equivalent width varies from 47 ° (Cabrit et al. 1990) up to A 202 ° (Muzerolle et al. 1998). These variations are accomA panied by large continuum flux variations: for example the stellar B -magnitudes vary from 10.7m to 13.6m (Herbig & Bell 1988). BP Tau is a single star at least down to 0.01 (Bernacca et al. 1995). According to Johns-Krull (1999b) the effective temperature of the star is 4060 K (spectral type K7) and v sin i = 10.2 ± 1.8 km s-1 . It is usually assumed that the distance to BP Tau is D = 140 pc, i.e. equal to the average distance to the Taurus star forming region. With this value and AV 0.5m Gullbring et al. (1998) found that the luminosity and the radius of the star are L 0.93 L and R 2.0 R respectively. However according to HIPPARCOS measurements the distance to BP Tau is equal to 53+17 pc (ESA -11 1997). Simon et al. (2000) found from IRAM interferometer observations that BP Tau is surrounded by a circumstelSend offprint requests to : L. Errico, e-mail: errico@na.astro.it Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 526555.

lar disk with external radius Rd = 108 ± 4 AU and inclination id = 30+4 degrees to the line of sight. They -2 also found from their data that the mass of BP Tau is M = 1.24+0.2 M if D = 140 pc. Within the approach -0.3 2 they used to derive stellar mass, M D, so the L /M ratio should not depend on the distance. From a comparison of this value with different theoretical track calculations M should be close to 0.8 M . Thus Simon et al. (2000) concluded that "the distance to the star may be closer than 140 pc but not as extreme as the HIPPARCOS value". In this study we accept the following set of BP Tau parameters: M = 0.8 M , R = 2 R and AV = 0.5m. The stellar radius is the only value in this set which depends on the distance: L D2 so R D for a fixed Teff . The main quantitative results of our paper are not sensitive to the current uncertainty in the value of D. For example the escape velocity V = (2GM /R)1/2 1/2 390 (140/D) km s-1 and the relative surface area of the accretion zone derived from Eq. (1) do not depend on the distance at all. A clear signature of circumstellar matter accretion was found in the optical spectrum of BP Tau ­ see Mundt & Giampapa (1982), Edwards et al. (1994), Fernandez et al. (1995), Johns & Basri (1995), Muzerolle et al. (1998), Gullbring et al. (1996, 1998), Alencar & Basri (2000). Johns-Krull et al. (1999b) measured the strength of the


558

L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau

Table 1. Log of HST Program 3845 on 30 July 1993. t0 , UT 17:33 17:34 17:46 17:47 17:49 19:08 Dataset z18e0103t z18e0104t z18e0105t z18e0106t z18e0107t z18e0108t Target W BP W W BP BP Cal Tau Cal Cal Tau Tau Grating G270M G270M G160M G160M G160M G160M , °A A° 2878­2924 2777­2823 1503­1539 1384­1402 1383­1419 1532­1568 N 1 2 1 1 6 5

40

20

0

average magnetic field finding B = 2.1 ± 0.3 kG. Furthermore Johns-Krull et al. (1999a) found strong circular polarization in the He i 5876 emission line, indicating a mean longitudinal magnetic field of 2.5 kG in the line formation region. They concluded that "accretion occurs preferentially along large-scale magnetic loops that occupy a small fraction of the stellar surface". Ardila & Basri (2000) (AB00) analyzed long wavelength low resolution IUE spectra and found, in the frame of the accretion shock model of Calvet & Gullbring (1998) (CG98), an average BP Tau accretion rate Mac = 1.6 â 10-8 M yr-1 such that the accretion zone (spot) occupies 0.4% of the stellar surface. They also found that the spot area increases as the square of Mac and claimed that "current models of the accretion process fail to reproduce such an effect". The evidence for a BP Tau stellar wind is relatively poor in the optical band: no blueshifted absorption features have been observed in the high resolution profiles of strong permitted optical lines (e.g. Balmer, Na i D, Ca ii, O i or He i). At the same time Gullbring et al. (1996) noted "the blue asymmetric appearances in the profiles of the higher Balmer lines". Usually the centroids of [O i] 6300, [O i] 5577 and [S ii] 6731 lines are blueshifted in CTTS spectra indicating a matter outflow with terminal velocity up to 300 km s-1 . But in the case of BP Tau the centroids of [O i] lines are redshifted and the [S ii] 6731 line profile is entirely situated between +100 and +200 km s-1 (Hartigan et al. 1995). In the frame of our detailed analysis of CTTS UV spectra (Errico et al. 2000; Lamzin 2000a, 2000b; Lamzin et al. 2001), in this paper we interpret the spectra of BP Tau obtained with HST and IUE satellites in order to derive the main parameters of the accretion shock and outflow/inflow geometry of this star.

2790

2795

2800

2805

Fig. 1. BP Tau spectrum in the vicinity of Mg II 2800 ° douA blet. The radial velocity of main features is shown in km s-1 relative to the rest wavelength.

20

10

0 1545 1550 1555 1560

Fig. 2. BP Tau spectrum in the vicinity of C iv 1550 ° doublet. A The line profile of the C iv 1548.20 ° is shown in the insert at A top right corner. X -axis of the insert is labelled in km s-1 .

2. Observations
HST observed BP Tau on 30 July 1993 (Program 3845) in 3 spectral bands with the Goddard High Resolution Spectrograph in the medium resolution mode and Detector 2. The wavelength calibration lamp was used in nearly the same spectral bands. The starting times of each observation, Archive Dataset names, grating designation, observed spectral bands and number of independent exposures (RPTOBS+1 parameter) are presented in Table 1. The spectra were adopted from the HST Archive,

recalibrated using the most up-to-date reference files and processed with IRAF v2.11 and STSDAS/TABLES v2.0.2 software as recommended in Chapter 36 of the "HST Data Handbook". The standard "pipeline" wavelength calibration was improved by using the STSDAS waveoff task and respective W CAL observations ­ see Table 1. All wavelengths are corrected for the orbital motion of the HST and the Earth and are presented for vacuum if < 2000 ° and for A air otherwise. The Van Hoof (1999) electronic database was used to identify lines in BP Tau spectra along with an atlas of H2 lines (Roncin & Launay 1994). To improve the signal-to-noise (S/N) ratio we combined all independent exposures for each spectral band and additionally smoothed them via a 4-point running mean, so the resulting spectral resolution is near 15 km s-1 . The five spectral intervals shown by thin lines in Figs. 1­3 indicate the "dead" diodes of Detector 2. High resolution spectra in the 1850-3350 ° spectral A band observed with IUE satellite were adopted from the INES database (http://ines.oat.ts.astro.it). The spectra were not additionally processed. Archive dataset names, date and universal time of the beginning of each


L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau Table 2. High resolution IUE spectra of BP Tau.
8

559

Dataset LWR LWP LWP LWP 11130 06963 09282 09417

Date 19811985198619860710101024 22 10 26

t0 , UT 06:41 02:00 00:51 13:20

t, h 3.5 3.0 4.5 4.5

FESc 218 236 176 162

4

0

1390

1400

Fig. 3. BP Tau spectrum in the vicinity of Si iv 1400 ° douA blet.

observation and exposure time t (in hours) are presented in Table 2. The last column shows the counts of the IUE fine error sensor (FES), which give an estimation of the stellar visual brightness.

3. Results
° The spectrum of BP Tau in the 2787-2810 A band is shown in Fig. 1. The continuum flux was derived from a featureless part of the spectrum between 2807 ° and A 2822 ° Fc = 1.3 ± 0.3 â 10-14 erg s-1 cm-2 °-1 , i.e. A: A Fc > 4c . This value coincides within 1c with the continuum flux derived between 2778 ° and 2788 °. We can see A A strong Mg ii h and k lines with a narrow interstellar (IS) absorption feature slightly redshifted relative to the rest wavelength, in agreement with the BP Tau radial velocity Vr = 15.8 ± 1.0 km s-1 (Hartmann et al. 1986). Two additional absorption features are present in the blue wings of both lines at -60 and -150 km s-1 . The total flux of both Mg ii uv1 lines is about 1.8 â 10-12 erg s-1 cm-2 . The upper levels 3p 2 P0/2,3/2 of Mg ii h and k reso1 nant lines are also the low levels of the Mg ii uv 3 triplet. An emission feature is present near the expected position of two lines of the triplet (2797.99 and 2798.06 ° but A), it is superimposed on the red wing of Mg ii h line, so it is difficult to say anything about these line profiles. The profile of the third ( 2790.84 ° line of the uv 3 mulA) tiplet apparently consists of three emission components. The maximum of the first (more or less symmetric) component is redshifted by +21 km s-1 relative to the rest wavelength. The second component shows a flat top with the center blueshifted by -53 km s-1 , i.e. it corresponds to the position of the first absorption feature in the blue wings of h and k lines. The third emission component of ° the 2790.84 A line seems to be blueshifted by the same value ( -150 km s-1 ) as the "high-velocity" absorption features of the resonant lines. However a few pixels between the second and the third components of the line fall on to the "dead" diodes of Detector 2, so the red wing of the third component could be a bit distorted ­ see Fig. 1.

We did not find any more emission or absorption features in the spectrum with intensity >2c , while relatively strong (emission or absorption) lines of the Fe ii uv 234 multiplet are present in analogous HST/GHRS spectra of RU Lup (Lamzin 2000a), RY Tau (Lamzin 2000b), RW Aur (Errico et al. 2000) and DF Tau (Lamzin et al. 2001). The spectrum of BP Tau in the vicinity of the C iv 1550 ° doublet is shown in Fig. 2. Its S/N ratio is A rather low, so the continuum is underexposed and we can only identify the lines of C+3 and molecular hydrogen. The flux ratio of the C iv 1550 ° doublet compoA nent is close to 2:1, and the total flux of the doublet is 2.4 â 10-13 erg s-1 cm-2 , this is 4% less than the average value found by Valenti et al. (2000) from all IUE spectra. The normalized profile of the C iv 1548.20 ° line A is shown in the insert in the upper right corner of the figure. The line center is redshifted, but it is difficult to judge the redshift value or the reality of the features in the blue wing of the C iv 1550.77 ° line. A ° H2 lines R(3) 1547.34 ° and P(5) 1562.39 A have a A common upper level v = 1,J = 4 of the first excited electronic configuration 2p B1 + of the molecule. If the u lines are optically thin their expected flux ratio is 1.1 (Abgrall et al. 1993), so the R(3) line contributes no more than 10% to C iv 1550 doublet flux. We identified (but with some doubts) the emission feature at 1557 ° A ° line (B1-8 transition) of H2 . with the R(6) 1556.87 A Some useful information can be extracted from the fact that the R(6) 1556.87 ° line is at least twice as weak as the A P(5) 1562.39 ° line. Pumping of these lines occurs from A these levels with a practically identical excitation energy and the absorption coefficients of the pumping transitions are almost the same. Both lines are pumped by quanta from the red wing of L line, but the redshifts for the P(5) and R(6) lines are +100 km s-1 and +17 km s-1 respectively. We conclude therefore that the L line has a deep central depression due to self-absorption as in the case of the solar chromosphere (Mihalas 1978). Unfortunately the BP Tau spectrum in the vicinity of the Si iv 1400 doublet is also very noisy ­ see Fig. 3. Except A Si iv lines, only the H2 B0-5 P(2) 1398.95 ° line can be identified unambiguously. The R(0) 1393.72 ° molecular A hydrogen line also has the same upper level and falls into the spectral bands of Figs. 2 and 3. If both lines are optically thin then the R(0) line flux should be 1.9 times larger than the P(2) line flux, so it contributes significantly to the observed flux of the Si iv 1393.76 ° line. A


560

L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau

0.8

0.4

0

0

0.5

1

Fig. 4. Semiforbidden flux ratio diagram for CTTS accretion shock. Flux ratios of Si iii] 1892, O iii] 1665 and C iii] 1909 lines are plotted along X and Y axis. See text for details.

particle density of the infalling gas in BP Tau does not exceed 1012 cm-3 . The escape velocity on BP Tau is V 390 km s-1 . Probably the angle between the rotational axis and the circumstellar disk is small or even zero, so the equatorial velocity of the star Vrot 20 km s-1 . Thus the coro2/3 tational radius Rc = R (V /Vrot ) 7 R . Bearing in mind the large magnetic field strength of the star we can suppose that the accretion radius is close enough to this value so the expected average infall gas velocity V0 is not far from 350 km s-1 . It then follows from Fig. 4 that N0 5 â 1011 cm-3 . We adopt this value and V0 = 350 km s-1 as the average parameters of the accretion shock. The Si iii] 1892 ° line originates within the H ii preA shock zone of which the extension is much less than the stellar radius for a given set of V0 and N0 parameters (Lamzin 1998). The line formation region can be considered as a plane-parallel slab and the emissivity of the line can be characterized by the specific intensity In (erg s-1 cm-2 ster-1 ) in the direction perpendicular to the shock front. The observed Si iii] 1892 ° line flux F can be A expressed through this value and the total accretion shock surface area Sac as follows (Gomez de Castro & Lamzin 1999): F= In Sac 2D 2 2.0 â 10
-19

One can expect from an analogy with the spectra of RU Lup (Lamzin 2000a), RW Aur (Errico et al. 2000) and DF Tau (Lamzin et al. 2001) that two additional H2 lines will also be superimposed on both lines of Si iv doublet. They are R(1) 1393.96 ° and P(3) 1402.65 ° lines A A with B(0, 2) common upper level and a transition probabilities ratio 1.4 (Abgrall et al. 1993). There are some marginal indications that the lines are indeed present in the BP Tau spectrum. If so the real flux of the Si iv doublet can be almost twice as low as the observed value of 8.7 â 10-14 erg s-1 cm-2 . Beside that a number of bad pixels are at the position of Si iv 1393.76 and 1402.77 ° A lines, so we will not discuss these.

f In

R 2R

2

140 D, pc

2

(1)

4. Discussion 4.1. Accretion flow parameters
The O iii] 1665, Si iii] 1892 and C iii] 1909 semi-forbidden lines originated in the CTTS accretion shock are optically thin and are expected to have comparable intensities (Lamzin 1998). The flux ratio of these lines depends on the infall gas velocity V0 and the infall gas particle density N0 , but does not depend on the geometry of the accretion zone. Thus one can derive V0 and N0 values (Gomez de Castro & Lamzin 1999). According to JohnsKrull et al. (2000) the average flux ratio of the C iii] 1909 and Si iii] 1892 lines derived from low resolution IUE spectra of BP Tau is 0.33 ± 0.04. Unfortunately the flux ratio of the O iii] 1665 and Si iii] 1892 lines is unknown, so in Fig. 4 we plot only the corresponding vertical line. We can see even from this restricted dataset that the average

where = exp(-0.92 A) is the interstellar extinction coefficient at = 1892 ° (it is 0.3 if we assume AV = A 0.5m and the standard extinction law, Seaton 1979); f is the relative fraction of the BP Tau surface taken up by the accretion zone. We assume that the total surface area of the accretion zone is two times larger than the observed area ­ a factor 2 in the denominator of Eq. (1). The average observed Si iii] 1892 line flux is F 4.1 ± 0.2 erg s-1 cm-2 (Valenti et al. 2000) and In 8 â 105 erg s-1 cm-2 ster-1 for our set of V0 and N0 parameters ­ see Fig. 1 of Gomez de Castro & Lamzin (1999). From Eq. (1) follows f 0.25, i.e. the accretion spot takes up 25% of the stellar surface on average. Now we can find the accretion rate, the energy flux F of the accretion flow and the accretion luminosity from the relations: M
ac

= 0 V0 Sac ,

F = 0

V03 , 2

L

ac

2 V = Mac 0 2

(2)

where 0 = µi mH N0 is the infall gas density, mH is the hydrogen atom mass and µi = 1.3 is the average molecular weight of the infalling gas which we assume to have a cosmic abundance (Allen 1973). It follows that Mac 3.6 â 10-8 M yr-1 , F = 2.3 â 1010 erg s-1 cm-2 and Lac 0.37 L i.e. 39% of the BP Tau bolometric luminosity. According to CG98 Mac = 2.9 â 10-8 M yr-1 and Lac 0.18 L . The average values of the accretion rate


L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau

561

and luminosity derived from 45 IUE spectra by AB00 are: Mac = 1.6 â 10-8 M yr-1 , Lac 0.09 L . Our values are about two times larger. This discrepancy can be due to variability and difference in the infall gas velocity adopted by CG98 and AB00, who assumed M = 0.49 M so they adopted V0 = 300 km s-1 instead of our V0 = 350 km s-1 . Thus, our CG98 and AB00 values of Lac and Mac are in a reasonable agreement. The accretion energy flux and filling factor of the accretion spot found in the quoted papers are also similar to each other: f = 7 â 10-3 , log F = 11.5 (CG98) and ¯ ¯ f = 4.4 â 10-3 , log F = 11.7 (AB00). The coincidence of the results is not surprising because AB00 used the accretion shock model calculated by CD98. Our f and F values, derived within the frame of the Lamzin (1998) accretion shock model, differ from the AB00 and CG98 values by more than one order of magnitude. Obviously this discrepancy is the result of the different approach used. As the R/D-ratio does not depend on distance, our accretion shock parameters are also independent on the current distance uncertainty, but they depend on the interstellar extinction value through the coefficient in Eq. (1). For example if AV = 0.3m instead of 0.5m then our f , Mac and Lac values would be 1.6 times less. We suppose however that the main shortcoming of our approach is the assumption on the homogeneous character of the accretion. Equation (1) assumes that In is the same in all points of the accretion zone, otherwise it should be interpreted as an average value. The same is true for the diagnostic diagram in Fig. 4 used to derive the infall gas density. But if the stellar magnetic axis is strongly inclined to the disk plane (see below) one can expect that V0 and N0 values cannot be the same at all points of the accretion zone. Therefore we could overestimate the values for f , Mac and Lac , but not by more than 2­3 times. CG98 derive the accretion shock parameters by comparing the calculated continuum spectral energy distribution with the observed one. The excess of the Balmer continuum emission originates in the upper layers of the heated stellar atmosphere at the bottom of an accretion column. To calculate the spectrum of this region CG98 use a relatively simple approach which gives them the possibility to find only the emergent flux E i.e the specific intensity I () averaged over cosine of the angle between the line of sight and the normal to the plane-parallel gas slab. Then they fit the observed spectrum varying two parameters: N0 and V0 or their equivalent. Sac or f can be derived from the observed absolute flux with a known F value. GC98 and AB00 derive the accretion shock parameters "averaged over cos " ­ see e.g. AB00 and their Table 3. Meanwhile the BP Tau veiling continuum originates in the region where the temperature increases outwards, so the larger the larger the continuum specific intensity should be. This is also why we observe limb brightening effect in the solar chromosphere (Mihalas 1978). Obviously the average over cos results in the preferential contribution of specific intensities with low into the resulting E -value.

These are the smallest intensities, so if the accretion zone is observed within a large range of then the comparison of an observed spectrum with E results in an overestimation of F and a corresponding underestimation of f ­ one can derive the same accretion luminosity with a large energy flux density and a small f or vice versa. The extent of the C iv 1548 ° line wings does not A exceed 200 km s-1 (Fig. 2) while the infall gas velocity is at least 1.5 times larger. It means that the "average" angle between the line of sight and the accretion zone is large enough, so CG98 and AB00 have overestimated F and underestimated f . What is more, CG98 also assume that the accretion is homogeneous, as we do. It cannot be said without special calculations and a specification of the accretion zone geometry how large the resulting error is. But the CG98 approach seems too crude to believe 2 that the f Mac dependence found by AB00, applying this model, is real. CG98 found that highly veiled CTTSs have similar energy fluxes F as less veiled stars but a larger surface coverage f . On the other hand as far as it becomes possible to fit the observed continuum spectral energy distribution more or less well, the value of Lac appears to be automatically correct. The spread of CTTS infall gas velocity is relatively small (250 < V0 < 400 km s-1 ), so according to Eq. (2) the accretion rate derived by CG98 is likely. In fact our values of Lac and Mac are in reasonable agreement with CG98's values taking into account the BP Tau variability and the uncertainty in the V0 value. Both our methods and those of GC98 to derive accretion shock parameters of BP Tau are not perfect because they are based on the assumption that accretion is homogeneous. But an additional averaging of the specific intensity over makes it impossible, in principle, to derive more accurate accretion shock parameters within the frame of the GC98 model. In contrast our approach gives the possibility to derive N0 and V0 values in each point of the accretion zone as well as its geometry if we succeed in observing the star at different phases of the rotational period (Lamzin 2000c).

4.2. The outflow
Two absorption features are present in the blue wings of the Mg ii h and k resonant lines in the HST spectrum. There are also two emission components in the profile of the Mg ii 2790.84 ° line blueshifted to the same values: A V -60 and -150 km s-1 . Thus at the moment of observations there were two separate regions of gas moving outwards between the star and the Earth with significantly different radial velocities. We call them low velocity (LV) and high velocity (HV) region. There are no Fe ii absorption lines in the spectrum shown in Fig. 1 including uv234 multiplet lines which are expected to be the strongest (Errico et al. 2000). The excitation energy of the uv234 multiplet low level (Ei 3.2 eV) is less than the energy of the Mg ii uv3 subordinate multiplet emission lines: h 4.4 eV.


562

L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau

accretion disk along magnetic field lines and later falls onto the stellar surface. In other words the LV region is the part of the accretion column that pro jects onto the stellar limb. We did not observe redshifted Fe ii lines nor Mg ii redshifted absorption features because at large distance from the stellar surface the inflowing matter pro jects out of the stellar limb. At the moment of the observation and at distance <0.1 R these ions disappear due to photoionization (Lamzin 1998). The HV region is a real outflow, i.e. the base of the magnetospheric wind. In this picture (Fig. 5) blueshifted HV and LV Mg ii absorption features should appear and disappear periodically, unless the stellar rotation axis is also strongly inclined to the disk rotational axis. A similar behavior can be probably seen in the resonant lines of Na i and Ca ii as well.

4.3. Detection of a flare and Mg II line profile variability
AB00 found that the BP Tau Mg ii line flux varies in the range 0.8-8.2 10-12 erg s-1 cm-2 ­ see also Simon et al. (1990) and Gomez de Castro & Franqueira (1997). We found a flux 1.8 â 10-12 erg s-1 cm-2 , so the HST spectrum represents a low level of stellar activity. To investigate the reason of the doublet flux variability we plot in Fig. 6 all IUE high resolution spectra in the vicinity of ° 2800 A along with the HST spectrum from Fig. 1 for reference. Gomez de Castro & Franqueira (1997) compared the same IUE spectra with the HST/GHRS one ­ see their Fig. 4 ­ but they did not recalculate vacuum wavelengths of the HST spectrum to air wavelengths. For this reason the IS absorption features in Mg ii uv1 lines appear redshifted in their plot and were misinterpreted as a signature of a "warm infalling gas". Unfortunately the quality of IUE spectra is not good enough to definitely detect either HV or LV absorption features if they have the same strength as in the HST spectrum. Nevertheless it is possible to say something about the variability of Mg ii line profiles. In LWR11130 and LWP09282 spectra the intensity of the lines was 1.5-2 times as big as one of the HST spectrum, the red wings of h and k lines increased more significantly than the blue ones. In the LWP06963 spectrum the situation is opposite: the blue part of the line profiles increased much stronger than the red ones. According to AB00 the timescale of the Mg ii 2800 ° A doublet flux and equivalent width variability is much less than the stellar rotation period, which means that the variability is the result of some non-stationary process rather than rotational modulation or occultation by dust clouds. The intensity of the h and k lines in the LWP09417 spectrum is close to that of the HST spectrum but the extent of their blue wings is significantly less (Fig. 6). By the way, the doublet line blue wings also look less extended in other IUE spectra in comparison with the HST spectrum. It means that the observed geometry of the gas flow is also time variable in agreement with the predictions of

Fig. 5. Qualitative scheme of the geometry of BP Tau inflow and outflow with both open and close magnetic field lines (solid lines). The broken lines show the line of sight.

The absence of Fe ii lines in absorption is possible if the gas temperature in both LV and HV regions is relatively ° low and the excitation of the Mg ii 2790.84 A line is due to radiation rather than to electron collisions. Cold external gas absorbs Mg ii h and k line quanta produced by an accretion shock resulting in the pumping of 2 P0/2,3/2 levels. 1 Quanta from the red wing of Mg ii h line excite the 2 D3/2 level and its de-excitation occurs mainly via emission of the Mg ii 2790.84 line (2 D3/2 2 P0/2 transition). 1 The Fe i 2788.17 ° line of the a5 F-y5 G0 multiplet is A also absent in the BP Tau spectrum, while the excitation energy of the a5 F term is only 0.86 eV and the log gf value of the line is -0.05 (Nave et al. 1994), i.e. the line is expected to be strong unless the relative number of iron atoms is small enough. We conclude that the gas in the LV and HV regions is cold but that the iron and magnesium are singly ionized. We suppose that the source of ionization is H i L quanta produced by the accretion shock as in the case of RW Aur (Errico et al. 2000). Indeed the BP Tau H line equivalent width is 100 ° which A, means that 1% of stellar luminosity is radiated in this line and the L line is expected to be even stronger. In this case both LV and HV regions should be close to the star. At the same time the circumstellar disk is seen nearly pole-on: id 30 (Simon et al. 2000). Therefore the only possibility to observe two regions of moving outward gas pro jected onto the stellar limb is shown in Fig. 5. The proposed large inclination angle im of the stellar magnetic axis is in agreement with Muzerolle et al. (2001) who found im 70 modelling H i Balmer and Na i D line profiles. We suggest that the LV absorption component originates in the gas, which initially moves outwards from the


L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau
2790 10 5 0 10 5 0 10 5 0 10 5 0 2790 2800 2810 2800 2810

563

Fig. 6. High resolution IUE spectra of BP Tau in the vicinity of Mg II 2800 doublet (solid lines). Gaps in the curves respect to pixels with non-zero quality flag. HST/GHRS spectrum is plotted for comparison (thin lines). Rest frame positions of the Fe II 2811.8, 2812.1 ° lines are shown ­ see text for details. A

the inclined magnetic rotator model. Furthermore if the inclination of the stellar magnetic axis is large we can expect strong variability due to the non-stationary character of the accretion process. In the LWP06963 spectrum a strong emission feature is present at 2811.2 ° but it is absent in all other A spectra (Fig. 6). There is no reason to doubt the reality of this feature because its pixel fluxes are 4-5 times larger than the flux errors. We identified this feature with a blend of two Fe ii lines of the c4 D-s4 Do multiplet with a common upper level J = 5/2 (2811.81, 2812.14 ° A) blueshifted to 80 km s-1 relative to the rest frame position. The upper level of these lines can be pumped from the a4 P3/2 level (Ei 1 eV) by quanta from the blue wing (V -230 km s-1 ) of the H i L line. Fluorescent Fe ii lines pumped by L line quanta were observed in the RW Aur spectrum as well (Errico et al. 2000). The excitation energy of the Fe ii c4 D term is Ei 7.5 eV so the fluorescent lines are definitely optically thin and their flux ratio should be equal to the ratio of transition probabilities. According to Nahar (1995) the expected 2811.81 ° line is 2.5 times weaker than the 2812.14 ° A A line, in agreement with the observations. We also searched for other Fe ii lines with an s4 Do/2 upper level which fall 5 in the range of LWP06963 spectrum. Unfortunately they are 2-10 times weaker than the 2811.2 ° emission feaA ture and fall in regions with low S/N-ratio, so we could not identify them. We believe that the observed increase of h and k line blue wing intensity is accompanied by the increase of

hydrogen L line intensity, resulting in the appearance of blueshifted Fe ii fluorescent lines. The BP Tau visual brightness was also large; during this flare ­ see Table 2 ­ only two stellar spectra out of 49 from the INES database were observed with larger values of FES counts. In the frame of the inclined rotator model these facts can be interpreted as follows. The increase of the accretion rate results in an increase of the continuum as well as the Mg ii 2800 and L line emission from the accretion shock. We assume that the flare has occurred when the stellar magnetic pole was on the visible hemisphere of the star, with a rotational period phase shifted to relative to the phase shown in Fig. 5. The regions of the accretion column where gas moves away from the observer are pro jected onto the star limb, while the regions with negative radial velocities pro ject out of the limb. Thus the red wings of h and k lines are attenuated due to the absorption inside the accretion column and the intensity of the blue wings is increased. In other words we also observe an "absorption feature" in h and k lines, but it is redshifted and much wider than the LV feature because, apparently, the radial velocity range is larger. It was shown in Sect. 3 that the accretion shock L line has a deep central depression, i.e. its profile is double peaked. During the flare only those regions of the external gas flow could produce observed Fe ii fluorescent lines which had a Doppler shift of -230 km s-1 relative to one of the two peaks of the accretion shock L line. If we suppose it was the red peak, then the Fe ii fluorescent line


564

L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau

formation region was located in the stellar wind or in that part of the accretion column where gas moved away from the star (Fig. 5). This is the reason of the blueshift observed in the iron lines. The energy release was very large during the flare. The observed flux of the Fe ii 2811.2 emission feature is 3 â 10-13 erg s-1 cm-2 and its dereddened luminosity is 2 â 1030 erg s-1 , assuming D = 140 pc, AV 0.5m and the normal extinction law (Seaton 1979). Meanwhile the lines of the b4 D-s4 D0 multiplet (1547.24, 1547.80 and 1550.09 °) have a transition probA ability 30 times larger than that of the c4 D-s4 D0 multiplet which is the largest value among Fe ii lines with an s4 D0/2 upper level. The excitation energy of the b4 D term 5 3.9 eV) to expect that the lines is high enough (Ei of this multiplet are optically thin. Thus the luminosity of all iron fluorescent lines originating from the s4 D0/2 5 upper level was near 3% of the BP Tau total luminosity L 0.9 L (Gullbring et al. 1998). Obviously only a small fraction of L quanta was re-radiated in the Fe ii fluorescent lines, so the total excess L luminosity at the moment of the LWP6963 observation was comparable with L , and Lac apparently even exceeded L . In the case of RU Lupi Lamzin et al. (1996) also observed a flare like event with an excess luminosity larger than the stellar bolometric luminosity. A large energy release in Fe ii fluorescent lines means that the optical depth 0 of the iron 1214.73 ° pumping A line (a4 P3/2 -s4 D0/2 transition) exceeds unity ­ we take 5 0 = 3 as a lower limit. Let ni be the relative population of the a4 P3/2 level, nFe II the relative abundance of the Fe+ ion, NH the hydrogen (both ionized & neutral) particle density (cm-3 )and l the radial extent of the line formation region. We have 0 = k0 ni n
Fe I I

The soft X-ray emission of the CTTS accretion shock is practically independent on the infall gas density (Lamzin 1999), i.e. the theory predicts the absence of a correlation between the variations in optical and X-ray bands in agreement with the observations of BP Tau (Gullbring et al. 1997).

5. Conclusions
The main results of our investigation can be summarized as follows. ­ BP Tau activity can be explained in the frame of a disk accretion model if we assume that the stellar magnetic axis is strongly inclined to the disk plane. In the frame of this hypothesis, the nontrivial character of line flux and profile variability is the result of the nonstationary accretion observed at different orientations of the stellar magnetosphere relative to the Earth. ­ We found the following average parameters of BP Tau accretion process: relative surface area of the accre 3.6 â tion zone f 0.25, accretion rate Mac 10-8 M yr-1 , accretion energy flux F = 2.3 â 1010 erg s-1 cm-2 and accretion luminosity Lac 0.37 L . All values were derived assuming that the accretion is homogeneous, i.e. the infall gas density and infall gas velocity are the same in all points of the accretion zone. ­ Ardila & Basri (2000) found that the accretion spot surface area is indeed proportional to the square of the accretion rate, applying the CG98 accretion shock model. We argue that this model is too crude to believe their result. ­ A strong flare in the Fe ii 2811.8, 2812.1 ° lines was A detected. It was probably produced by an increase of the accretion rate. During the flare, the accretion luminosity was comparable to or even larger than the stellar bolometric luminosity. ­ Finally we identified a number of lines in the UV spectra of BP Tau with H2 lines. We were strongly restricted by the quantity and quality of analysed data. It is necessary to observe the star in a wider UV spectral range and at different phases of the rotational period in order to get information on the geometry and parameters of the BP Tau accretion shock. Fe ii fluorescent lines, pumped by accretion shock L radiation, can be used as a tracer of the gas flow in the vicinity of the star.
Acknow ledgements. This work was accomplished during the stage of S. L. at Astronomical Observatory of Capodimonte (OAC) in Naples. S. L. thanks the OAC staff for support as well as the Russian Fund of Fundamental Research for the grant 99­02­17184. The authors thank the anonymous referee for his helpful remarks.

Fe NH l 3

where k0 is the absorption cross section in the center of the line and Fe 4 â 10-5 is the iron relative abundance (Allen 1973). Assuming that the local line broadening is due to thermal motion and adopting the atomic data by Nahar (1995) we can derive from Eq. (9.43) of Mihalas (1978) that k0 4.5 â 10-16/ T/104 cm2 , where T is the gas temperature. As long as nFe II 1 we find that, during the flare, the hydrogen column density NH l is larger than 1.5 â 1020 /ni cm-2 . Adopting statistical weights and excitation energies of the Fe ii levels by the NIST Database (http://physics.nist.gov) we find that with T < 104 K Boltzman (LTE) population of the a4 P3/2 level is less than 0.014, it is the upper limit to the ni -value. Therefore during the flare NH l was 1022 cm-2 . From the Cruddace et al. (1974) data we find that at this moment the BP Tau X-ray emission should be strongly attenuated below 1 keV (1kev > 2) if the iron fluorescent line formation region is pro jected onto the stellar limb. There is direct evidence that in some cases CTTS X-ray variability can be due to circumstellar envelope variable extinction (Walter & Kuhi 1981; Kastner et al. 1999).

References
Abgrall, H., Roueff, E., Launay, F., Roncin, J.-Y., & Subtil, J. L. 1993, A&AS, 101, 273


L. Errico et al.: UV spectra of T Tauri stars from the HST and IUE satellites: BP Tau Allen, C. W. 1973, Astrophysical Quantities (Athlone Press) Alencar, S. H. P., & Basri, G. 2000, AJ, 119, 1881 Ardila, D. R., & Basri, G. 2000, ApJ, 539, 834 Bernacca, P. L., Lattanzi, M. G., Porro, I., Neuhauser, R., & Bucciarelli, B. 1995, A&A, 299, 933 Cabrit, S., Edwards, S., Strom, S. E., & Strom K. M. 1990, ApJ, 354, 687 Calvet, N., & Gullbring, E. 1998, ApJ, 509, 802 Cruddace, R., Paresce, F., Bowyer, S., & Lampton, M. 1974, ApJ, 187, 497 Edwards, S., Hartigan, P., Ghandour, L., & Andrulis C. 1994, AJ, 108, 1056 Errico, L., Lamzin, S. A., & Vittone, A. A. 2000, A&A, 357, 951 ESA, 1997, The Hipparcos and Tycho Catalogues, ESA SP1200 Fernandez, M., Ortiz, E., Eiroa, C., & Miranda, L. F. 1995, A&AS, 114, 439 Gomez de Castro, A. I., & Franqueira, M. 1997, ApJ, 482, 465 Gomez de Castro, A. I., & Lamzin, S. A. 1999, MNRAS, 304, L41 Gullbring, E., Petrov, P. P., Ilyin, I., et al. 1996, A&A, 314, 835 Gullbring, E., Barwig, H., & Schmitt, J. H. M. M. 1997, A&A, 324, 155 Gullbring, E., Hartmann, L., Briceno, C., & Calvet, N. 1998, ApJ, 492, 323 Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736 Hartmann, L., Stahler, S., & Mathieu, R. D. 1986, ApJ, 309, 275 Herbig, G. H., & Bell, K. R. 1988, Lick Observ. Bull., 1111 Johns, C. M., & Basri, G. 1995, AJ, 109, 2800 Johns-Krull, C. M., Valenti, J. A., Hatzes, A. P., & Kanaan, A. 1999a, ApJ, 510, L41

565

Johns-Krull, C. M., Valenti, J. A., & Koresko, C. 1999b, ApJ, 516, 900 Johns-Krull, C. M., Valenti, J. A., & Linsky, J. L. 2000, ApJ, 539, 815 Kastner, J. H., Huenemoerder, D. P., Schultz, N. S., & Weintraub, D. A. 1999, ApJ, 525, 837 Lamzin, S. A., Bisnovatyi-Kogan, G. S., Errico, L., et al. 1996, A&A, 306, 877 Lamzin, S. A. 1998, Astron. Rep., 42, 322 Lamzin, S. A. 1999, Astron. Lett., 25, 430 Lamzin, S. A. 2000a, Astron. Lett., 26, 225 Lamzin, S. A. 2000b, Astron. Lett., 26, 589 Lamzin, S. A. 2000c, Astron. Rep., 44, 323 Lamzin, S. A., Vittone, A. A., & Errico, L. 2001, Astron. Lett., 27, 313 Mihalas, D. 1978, Stellar Atmospheres (W. H. Freeman & Co., San Francisco) Mundt, R., & Giampapa, M. S. 1982, ApJ, 256, 156 Muzerolle, J., Hartmann, L., & Calvet, N. 1998, AJ, 116, 455 Muzerolle, J., Calvet, N., & Hartmann, L. 2001, ApJ, in press Nahar, S. N. 1995, A&A, 293, 967 Nave, G., Johansson, S., Learner, R. C. M., Thorne, A. P., & Brault, J. W. 1994, ApJS, 94, 221 Roncin, J.-Y., & Launay, F. 1994, J. Phys. Chem. Ref. Data, 4, 1 Seaton, M.J. 1979, MNRAS, 187, 73 Simon, T., Vrba, F. J., & Herbst, W. 1990, AJ, 100, 1957 Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034 Van Hoof, P. 1999, Atomic Line List, v 2.01, http://www.pa.uky.edu/peter/atomic Valenti, J. A., Johns-Krull, C. M., & Linsky, J. L. 2000, ApJS, 129, 399 Walter, F. M., & Kuhi, L. V. 1981, ApJ, 250, 254