EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 614 (June 2018) A&A, 614 (2018) A117 Full HTML
Free Access
Issue
A&A
Volume 614, June 2018
Article Number A117
Number of page(s) 10
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201832891
Published online 22 June 2018

© ESO 2018

1 Introduction

Classical T Tauri stars (cTTS) are pre-main sequence stars (PMSs) of low mass characterized by their rich emission line spectra with signatures of both accreting and outflowing gas. Their complex spectra are understood as a result of magnetospheric accretion of material from circumstellar disks in combination with emission from stellar winds and also winds arising in the surrounding disk. Spectral types can be deduced from the photospheric lines, but in cTTSs the stellar absorption lines are weakened due to overlying excess emission, an effect known as veiling, which has been related to continuous emission from the shocked region below the accretion streams close to the stellar surface. For reviews of properties and models of cTTS, see Petrov (2003), Bouvier et al. (2007), and Hartmann et al. (2016).

Some cTTSs have exceptionally strong and broad absorption components flanking the red wings of certain emission lines, called inverse P Cygni profiles. Such stars are sometimes referred to as YY Orionis stars after Walker’s (1972) spectroscopic survey of objects in two young clusters. The red-shifted components originate in gas falling towards the star with velocities of >100 km s−1.

One PMS system with outstanding YY Ori characteristics is the 1.3″ visual binary S Coronae Australis. At a distance of 130 pc (Reipurth & Zinnecker 1993), this corresponds to a projected separation of 169 AU if physical. Both components, usually designated as NW (north-west) and SE (south-east), have very rich emission line spectra and vary in brightness by up to two magnitudes in the photometric V band on timescales of hours and days and become redder with decreasing brightness (Walter & Miner 2005).

A number of spectroscopic studies of the YY Ori star S CrA followed the survey by Walker (see Appenzeller et al. 1986, and references therein; Krautter et al. 1990; Graham 1992; Hamann & Persson 1992) demonstrating that the red-shifted absorption components are present in different types of lines, and in cases when repeated observations were made, they were found to vary in strength. However, the NW and SE components were not clearly resolved until Prato et al. (2003), Takami et al. (2003), and Carmona et al. (2007) presented individual spectra.

Over the last decade, more detailed analyses of the spectral properties of both components have been made. On the basis of high-resolution spectra obtained of each component at the Very Large Telescope (VLT, ESO), Gahm et al. (2008) focussed on the nature of the veiling, which is variable and on occasion extraordinarily strong. It was demonstrated that when the veiling becomes high, the photospheric lines are filled in by narrow emission lines unrelated to any continuous emission from shocked regions. This phenomenon was related to accretion-powered chromospheric emission in Petrov et al. (2011), and a similar concept was adopted in the models by Dodin & Lamzin (2012) and Dodin et al. (2013).

Emission lines with prominent red-shifted absorption were further studied in Petrov et al. (2014, hereafter Paper I) who found that the radial velocity of the deepest part depends on the excitation and ionization state of the line. On occasion, these components were exceptionally broad and strong in S CrA SE and it was demonstrated that they could be matched with an accretion flow as depicted in the standard magnetospheric models assuming a dipole magnetic field structure. Gas at different density and temperature is found at different distances from the stellar surface, resulting in the observed velocity pattern when viewed along a sight line through the magnetosphere.

The circumstellar disks and associated spectral signatures and disk winds have been discussed in Kessler-Silacci et al. (2006), Schegerer et al. (2009), Pontoppidan et al. (2011), Bast et al (2011), Vural et al. (2012), and Meeus et al. (2013). Recently Garcia et al. (2017, GRAVITY collaboration) resolved the inner gaseous disk in both components from interferometric observationsof their near-infrared (NIR) spectra. They found signatures of both accretion flows and jet streams in the Brγ and He I lines and obtained disk inclinations that are rather similar in the two stars, 28 and 22 degrees, implying that the systems are seen nearly face-on. Finally, the Herbig Haro objects HH 82 and 729 appear to be related to S CrA (Reipurth & Graham 1988; Wang et al. 2004).

In the present paper we revisit the material of high-resolution spectra obtained at VLT of both S CrA NW and SE with the purpose of deriving stellar parameters for the two stars and to disentangle the different spectral signatures of in-flowing and out-flowing matter in an attempt to match these with numerical models of the magnetospheric structures envisioned for the two stars.

The observations and reductions are described in Sect. 2. In Sect. 3 we spot-light different spectral features that are related to stellar winds and accretion streams, and which form the basis for our theoretical models of the objects presented in Sect. 4. The results are further discussed in Sect. 5, and we end with a summary in Sect. 6.

2 Observations

Eleven high-resolution spectra of S CrA were obtained between August 14 and 17, 2005 (MJD 24553597.03 to 24553600.02), with the UVES spectrograph in the dichroic-mode on the 8 m VLT/UT2 at ESO, Chile. The stars were aligned in the slit and were perfectly spatially resolved, therefore the two components were extracted individually. The spectra have a high signal-to-noise ratio, S∕N≳100, a spectral resolution of R ≈ 60 000, and a wavelength coverage from 3500 to 6700 Å. For more information on data reductions see Stempels et al. (2007). In addition, we have compiled a long-term record of the two S CrA components with UBV RI∕JHK photometry from 2003 to 2007, using the ANDICAM dual-channel imager on the SMARTS 1.3 m telescope at Cerro Tololo.

3 Results

3.1 Photometry

Both components in S CrA (referred to as NW and SE below) are strongly variable and become redder when fainter as is evident from our long-term photometry with SMARTS presented in the upper panel of Fig. 1. Most of the time, NW is brighter than SE, on the average by 0.m6 in the V band. The lower panel of Fig. 1 shows how the two stars varied in the I and B bands over the period when the VLT spectra were collected. The variations were moderate in the I band, but in B, changes of almost a magnitude occurred on timescales of one day in both stars.

The veiling changed dramatically over this period, especially in SE, but these changes are not tightly related to the brightness changes, because a contribution of varying line emission dominates when the veiling is high (Gahm et al. 2008). Instead, the main cause of variability is obscuration by passing circumstellar dust clouds in the line-of-sight. The slopes of the observed reddening of the stars in Fig. 1 is consistent with a normal reddening law, so the dust grains must be small.

Both stars show extraordinarily large variations in the optical bands. There is no clear periodicity in these changes judging from our long-term photometric monitoring, and large changes can occur over days as well as over weeks. Variations due to extinction events are not unusual among the TTSs (e.g., Herbst et al. 1994; Guo et al. 2018, and references therein).

thumbnail Fig. 1

Upper panel: colour-magnitude diagram from the long-term monitoring of S CrA (blue: S CrA SE, red: S CrA NW). Lowerpanel: individual I and B magnitudes observed over the period of spectral observations.

3.2 Stellar lines and parameters

The photospheric lines are completely obliterated by superimposed emission in the blue spectral region in both components. The top panel inFig. 2 shows spectra of NW and SE obtained on the four consecutive nights over the region 6000 – 6180 Å, one of the few spectral regions that is not contaminated with line emission. We determined the fundamental stellar parameters of both components by calculating synthetic spectra with the SME software (Valenti & Piskunov 1996; Piskunov & Valenti 2017) extended with routines to compensate for the level of veiling. The radiative transfer calculations in SME were performed using the Kurucz (1993) model atmospheres and atomic data from the VALD database (Piskunov et al. 1995; Kupka et al. 1999). We found that the photospheric spectra of the two stars are extremely similar, with identical v sin i, Teff, and log g.

The photospheric lines are weaker due to veiling compared to normal stars of the same spectral type. The veiling factor (VF) is defined asthe ratio of an assumed excess continuum to the stellar continuum. The lower panel in Fig. 2 shows the excellent agreement between the observed (black line) and synthetic (grey line) spectra. The only difference between the two components is the added level of veiling, here depicted as a continuous excess emission (horizontal dashed lines). However, as described inSect. 1, narrow emission lines entering inside the photospheric lines also contribute to the veiling. This effect of line-dependent veiling in highly accreting stars was studied in more detail by Rei et al. (2018).

In Table 1 we summarize the so-derived stellar parameters, which turn out to be remarkably similar for the two components. From the evolutionary tracks by Siess et al. (2000), assuming Teff = 4250 K and L* = 0.7 L⊙ for each component, we find M* = 1 M⊙ and R* = 1.4 R⊙ for both stars. This is consistent with the derived value of log g = 4.0. Hence, S Cr A is an example of a stellar twin. However, as can be seen in the table, the mean veiling is much higher in NW than inSE.

thumbnail Fig. 2

Upper panel: spectra of the stars obtained during the 4 nights over a spectral region relatively free from emission lines and veiling used to derive stellar parameters (time running upwards). Lower panel: good agreements between observed (black) and synthetic (grey) spectra. Different contributions from a continuous source of veiling is added (horizontal dashed lines).
Table 1

Stellar parameters.

3.3 Spectral signatures of gas flows

Emission lines of H and, He I, and neutral and ionized metals are prominent in the spectra of both stars. Figure 3 shows segments of spectra from August 14. The metal lines are broad in NW (~100 km s−1) on top of which narrow emission lines are present in several lines. In SE, the metal lines are as narrow as the photospheric lines and show velocity changes in anti-phase with the photospheric lines as reported in Petrov et al. (2011) suggesting that the line emission is confined to a certain area on the stellar surface, which is notaligned with the axis of rotation. Strong red-shifted absorption components flank the metal lines in SE.

Spectra taken over the four nights covering the Ca II H & K, Hγ, He I D3, and Na I D lines are shown in Fig. 4. Red- and blue-shifted absorption components are prominent in both stars. The red-shifted components were very stable in NW during the period, while in SE they were highly variable and on occasions extremely strong. The emission components vary on time scales of a day. We note that in the Ca II H profile ofNW the red peak is depressed by the central absorption of Hϵ, and the blue peak ofCa II H is also depressed by absorption in the cool wind. As a result, the H line becomes considerably weaker than the K line.

The presence of both blue- and red-shifted absorption components in the Balmer lines were recognised early by, for example, Rydgren (1977) and Edwards (1979) and at the time it was discussed whether the blue-shifted components could arise in winds as well as in accretion flows. The blue-shifted components in the H & K and D lines reach below the continuum in both stars and must be related to a wind. Both the strength and velocity shifts of these absorption features are larger in NW, indicating a higher mass-loss rate in this star. In NW, neutral gas moves at 120 km s−1 relative to the star as seen in Fig. 5. Obviously, this absorption line originates from gas that has been launched at high speed into the line-of-sight. In addition, narrow interstellar absorption components are present.

Figure 6 shows the observed Balmer line profiles in NW from spectra taken on August 14 and in SE on August 16, the dates of minimum veiling for each star. Common for both stars are the double-peaked profiles with a broad ‘central’ absorption component, strongest in Hα, and a broad red-shifted absorption component in the higher Balmer lines. The strong central absorption component in NW is typical for a disk wind observed at high inclinations of the disk relative to the sky plane (e.g. Kurosawa et al. 2011).

In the case of NW, the deepest part of the wind feature is at about −70 km s−1 in Hα, −20 km s−1 in Hβ, and is close to zero velocityin Hγ and Hδ. This reflects the outward accelerated wind since in the higher Balmer lines the absorption saturates in the inner and densest part of the wind, while in Hα the absorption is saturated already in regions of lower density and at larger expansion velocities.

The emission lines of He and metals are composed of one broad and one narrow component. According to Kwan & Fischer (2011) and Beristain et al. (1998) the broad component can originate in the magnetosphere loaded with accreting gas characterized by Te = 7000–8000 K and log Ne = 10–12, whilst the narrow components are usually attributed to the post-shocked gas with physical conditions close to those in the solar chromosphere, Te ~ 104, log Ne ≈ 12. The stellar wind is characterized by high temperature and low density, Te ≥ 104 K, log Ne = 9–11.

Regarding the red-shifted absorption components, several examples of line profiles observed in SE were shown in Paper I. In the August 14 spectrum of SE, the veiling is very high, and on this occasion many lines of different elements at different ionization stages show prominent inverse P Cyg profiles. The physical conditions along the accretion funnels could then be derived from a simple dipole model, and mass-loss rates of the order of 10−7 M⊙yr−1 were found. This result agreed with expected rates derived from the equivalent widths of Hydrogen and Helium emission lines, using the empirical relations between line luminosity and mass accretion rate (e.g. Fang et al. 2009). From these relations, one would conclude that the mass-loss rate is twice as large in NW compared to SE. However, the accretion components are stronger in SE for possible reasons discussed in Sect. 5.

thumbnail Fig. 3

Segments of spectra of S CrA NW (upper panel) with broad emission lines and S CrA SE (lower panel) with strong red-shifted absorption components. Narrow emission lines of metals are present in both stars.
thumbnail Fig. 4

Spectral segments of the two stars taken during 4 nights over the Ca II H & K (upper panel) and Hγ and He I plus Na I D lines (lower panels). Intensities are normalized to the continuum level and time is running upwards.
thumbnail Fig. 5

Ca II K and Na I D lines in SCrA NW show strong blue-shifted absorption components in addition to narrow interstellar lines (spectra from August 14). Velocities are astro-centric in this and similar figures.
thumbnail Fig. 6

Balmer line profiles. Left: S CrA NW (August 14). Right: S CrA SE (August 16). From top to bottom Hα, Hβ, Hγ, and Hδ with the continuum level marked with dashed lines.

4 Model calculations

In order to determine the main physical processes in the nearest environment of both components in S CrA and to estimate inclination angles, we have modelled the radiation coming from different regions in the environments of the stars, including one region of magnetospheric accretion (MA) and one with a disk wind (DW). First, we calculated the hydrogen lines of the Balmer series for each region separately, and then we combined these two regions, and calculated line profiles within the hybrid model as depicted in Fig. 7. We took into account all foreground absorption of the radiation from the star and the accretion region by the disk wind. Below, we briefly describe the main assumptions of both models regarding kinematics and temperature regimes. For a full description of these models, see Grinin & Tambovsteva (2011) and Tambovsteva et al. (2016).

thumbnail Fig. 7

Sketch of the hybrid model geometry assumed in the present calculations (not to scale). The inclination angle is measured from the symmetry axis (S).

4.1 Disk wind

The disk wind model is based on the theory of a magneto-centrifugal disk wind developed for the accretion disks surrounding black holes by Blandford & Payne (1982) and was used to study the accretion disks near young stellar objects by Pudritz & Norman (1986). Magneto-centrifugal disk winds in T Tauri stars have been investigated in many works (e.g. Königl & Pudritz 2000; Ferreira 2013, and references therein). Safier (1993) and Garcia et al. (2001) argued that the rising gas is rapidly heated by ambipolar diffusion up to ~10 000 K. Observations of jets from TTSs favour models of extended centrifugally driven winds (Ferreira et al. 2006).

The geometry and kinematics of the disk wind obtained from solutions of the MHD equations can be reconstructed with the method of parametrization (e.g. Kurosowa et al. 2006). In short, the main properties of the disk wind model are as follows. We divide the disk wind region by a number of streamlines extending from one point. For simplicity we assume that the disk wind consists of hydrogen atoms with a constant temperature (10 000 K) along the streamlines. Since the gas temperature just above the disk surface is not high enough to excite the hydrogen line emission, we exclude the base of the disk wind in the calculation and use a scaling factor for the temperature along each streamline.

The main model parameters are as follows. The footpoints of the disk, w1 and wN, refer to the wind launch region for the first and last streamlines, respectively. θ1 is the half openingangle between the first streamline and the vertical axis (Fig. 7). The poloidal velocity component v(l) changes along the streamlines, as given by v(l)=v0+(v∞−v0) (1−li/l)β .\begin{equation*}v(l) = v_0 + (v_{\infty} - v_0)\,(1 - l_i/l)^{\beta}\,. \end{equation*}(1)

Here, v0 and v∞ are the initial and terminal velocities and β is a parameter. We assume that v0 is the sound velocity in the disk wind. v∞ =fuK(wi), where uK(wi) is the Keplerian velocity at the footpoint wi on the disk surface for the ith streamline, and f is a scale factor.

The tangential velocity u(w) changes with a cylindric radius, w, as follows (Grinin et al. 2017):

In the corotation zone (u(w)∕uK(wi)) ≤ fc) u(w)=uK(wi)(w/wi).\begin{equation*}u(w)=u_K(w_i)(w/w_i). \end{equation*}(2)

In the zone of conservation of angular momentum u(w)=uK(wi)fc/(w/wi),\begin{equation*}u(w)=u_K(w_i) f_c/(w/w_i), \end{equation*}(3)

where fc is a parameter.

Two additional parameters enter, namely the mass-loss rate Ṁw and γ. The latter regulates the mass loading among the streamlines. We assumed that the accretion disk is transparent to radiation up to the distance of a sublimation radius. Beyond the sublimation radius, the disk is opaque obscuring the disk wind moving away from the observer.

4.2 Magnetospheric accretion model

We consider models in which the accreting matter falls onto the stellar surface radially within two cones. The cones extend to a radial distance Rout. We neglect the outer regions of MA, which contribute much less to the hydrogen line emission than the inner parts because of the lower gas temperatures. The same scheme of accretion was used by Calvet & Hartmann (1992).

The underlying assumptions are in short: the accreting gas flows from the outer radius Rout to the stellar surface within two cones with half-opening angles θin and θout and described by one radial (v) and one tangential (u) velocity component. v(r)=v*2+vesc2(R*r−R*Rout),\begin{equation*}v(r) = \sqrt{v_*^2 + v_{esc}^2 \left(\frac{R_*}{r}- \frac {R_*}{R_{\mathrm{out}}}\right)}, \vspace*{-4pt}\end{equation*}(4) u(w)=U*wR*,\begin{equation*}u(w) = U_* \frac{w}{R_*}, \end{equation*}(5)

where r is the distance from the stellar centre, w is the distance from the vertical axis (a cylindric radius), v* is the velocity at the stellarsurface, vesc is the escape velocity, U* is the rotational velocity at the stellar equator, and R* is the stellar radius. The electron temperature Te was chosen following the law Te(r)=Te(R*)(r/R*)−1/α,\begin{equation*}T_e(r) = T_e(R_*) (r/R_*)^{-1/\alpha} ,\end{equation*}(6)

where α is a parameter and Te(R*) the electron temperature at the stellar surface.

Since the gas temperature of accreting matter near the stellar surface can also be large because of heating by stellar radiation and ultraviolet and X-ray radiation from hotspots (Dodin 2018). Therefore, we also took into account broadening of line profiles due to the Stark effect in the MA emitting region.

4.3 Radiative transfer

For all emitting volumes around the star we performed non-LTE modelling of the radiative transfer. The algorithm of the calculations is described in detail by Grinin & Tambovsteva (2011). The solution of the radiative transfer problem was based on the numerical code developed for media with large velocity gradients. The source functions were calculated in the Sobolev approximation (Sobolev 1960), and the intensity of the radiation emergent at frequencies within a spectral line was calculated by exact integration over spatial coordinates in the approximation of full redistribution over frequency in the co-moving coordinate system. We considered 15 atomic levels plus continuum. An inclination angle i is measured from the symmetry axis (for i = 0 the star is seen pole-on).

4.4 Model calculations – results

In our calculations, we considered radiation emerging in the three first lines of the Balmer series, Hα, Hβ, and Hγ. Our aim was not to reproduce exact values of the maximum intensities in both peaks of the profiles because of their strong variability. For a satisfactory agreement with the observed profiles we required a match in line widths, shapes, and asymmetries, in addition to some specific observed features. For the comparison, we normalized both calculated and observed line profiles by dividing the intensities at each radial velocity by the maximum line intensity. Table 2 contains ranges of the main parameters of the disk wind and the magnetospheric accretion models. Parameters of the disk wind and magnetospheric accretion models presented in the paper are shown in Tables 3 and 4 respectively.

As a demonstration of the details of our calculations we present in Fig. 8 all separate components forming the Hα line profiles together with the resulting line profiles, one for a small inclination, 30° and another for 65°. The radiation from the MA region is shown in red (model MA1), and the radiation from the disk wind is shown in blue (model DW1). The total line profiles are depicted in black.

One can see that at small inclinations the disk wind provides an additional emission in spite of the fact that the mass accretion rate is large in comparison with the mass-loss rate. At large inclinations the disk wind absorbs the part of the radiation formed in MA and in the disk wind itself. This is why the lines are double-peaked in this two-component model.

Examples of how the Balmer line profiles depend on inclination angle are presented in Fig. A.1. The left panel shows Hα, Hβ, and Hγ profiles assuming magnetospheric accretion only (model MA1). Since the rotational velocities in both components of S CrA appear to be small (see Table 1) we neglect the tangential velocity component in our calculations. The right panel shows the same lines calculated with the hybrid model (model MA1 plus a disk wind model DW1). It is seen that the disk wind changes the line profiles substantially: at small inclination angles, it gives an emission at low radial velocities, whilst at larger angles (>50°) it absorbs the radiation from the MA region, and also its own radiation, giving rise to a central dip in the line profiles.Nevertheless, the accretion features are still prominent, more clearly seen in the weaker lines (Hβ and Hγ).

Our calculations showed that the best qualitative agreement between observed and calculated line profiles for S CrA SE is reached with the combined models MA1 and DW1. Figure 9 presents the theoretical line profiles of Hα, Hβ and Hγ normalized to the continuum (top panel), and below, these are overlapped on the observed profiles. Best fits were achieved for inclination angles of ~65°. All observedfeatures are accounted for in the model, for instance the ratios between the blue-shifted and red-shifted peaks, which are similar to those observed.

The Balmer line profiles for S CrA NW differ from those in S CrA SE as can be seen in Fig. 10. Firstly, the wind features are revealed more clearly: the blue peak is smaller than the red one in all lines, which is typical for an outflow with rotation. Secondly, the lines are extremely strong, reaching a factor of ~40 relative to the continuum level. One explanation is that the mass-loss rate is exceptionally large. Another possibility is that the stellar continuum has been depressed due to foreground extinction from dust. In this case, the strengths (equivalent widths) of emission lines appear to increase due to contrast. Finally, the central absorption components are very strong and broad. In any case, the presence of a strong disk wind is obvious. We therefore increased the value of the mass-loss rate by an order of magnitude (model DW2) and checked different MA models. Results of these calculations are shown in Fig. 10, where a qualitative agreement between observed and theoretical profiles has been reached.

Finally, we note that models of a polar wind seen nearly pole-on would produce a wide central absorption and double-peaked line profiles. However, in the simulations by Matt & Pudritz (2007) and Cranmer (2008), where the wind is driven by the thermal pressure, there is only a small region close to the star where the gas temperature is suitable for hydrogen line formation. In addition, this wind is optically thin, and the mass-loss rate is small, ~ 10−11 –10−9 M⊙ yr−1. Our calculations showed that the line profiles are symmetric around zero velocity, and the peak intensities are small and practically equal. Thus, this emitting region cannot be considered as an important contributor to the total line radiation.

In conclusion, our calculated line profiles based on the hybrid model reproduce all main features of the Balmer line emission in both S CrA components satisfactorily.

Table 2

Ranges of model parameters.

Table 3

Parameters of the disk wind model.

Table 4

Parameters of the magnetospheric accretion model.

thumbnail Fig. 8

Components of the Hα line profiles calculated in model MA1 + DW1 at 30° (left) and 65° (right): magnetospheric accretion (red), disk wind (blue), and the resulting line profile (black).
thumbnail Fig. 9

Balmer line profiles in S CrA SE. Upper panels: Hα, Hβ and Hγ line profiles as calculated from the hybrid model MA1 plus DW1 and with an inclination angle of 65°. Lower panels: same profiles normalized to peak intensity (solid curves) superimposed on the observed profiles (dashed curves).
thumbnail Fig. 10

As in Fig. 3 but S CrA NW. The line profiles are calculated in the hybrid model MA2 plus DW2 and with an inclination angle of 60°.

5 Discussion

Earlier photometric and spectroscopic studies of S CrA were related to the unresolved binary. A stellar luminosity of 1.5 L⊙ was derived by Hamann & Persson (1992) and of L* = 1.3 L⊙ by Prato et al. (2003). Carmona et al. (2007) resolved the NE and SE components and assigned spectral types G6 and K5, respectively. Vural et al. (2012) found the following parameters for NW and SE: spectral types K3 and M0, stellar luminosities 2.3 ± 0.7 and 0.76 ± 0.24 L⊙. This contradicts our results indicating that the two stars are very similar. It is likely that the earlier spectral type estimated for NW as compared to SE is a consequence of the highly veiled spectrum of NW. Our high-quality and high-resolution spectra enable more accurate determination of the spectral types.

The Balmer line profiles are different in NW and SE, and this fact may help us to understand which component was dominant in the spectra taken of the unresolved star. Although NW is usually referred to as the primary component, the high-resolution spectrum of S CrA in Hamann & Persson (1992) is actually that of the SE component. The Hα line profile reported by Reipurth et al. (1996) in their atlas of pre-main sequence stars also corresponds to the SE component.

From our theoretical modelling of the two S CrA stars in Sect. 4 we could find a good match between calculated and observed Balmer line profiles when considering a combination of a magnetospheric accretion flows and a disk wind, that is, the hybrid model. Best fits were obtained for inclination angles of ~65° for SE and ~60° for NW (for i = 0 the star is seen pole-on). Hence, not only are the two stars very similar in radius, luminosity, and surface temperature (Sect. 3.2), but the calculations indicate that the disks also have similar tilts relative the observer.

In this case, one expects from Fig. 7 that the lines-of-sight to both stars pass the trajectories of accreting gas falling at velocities close to the free-fall velocity just above the stellar surface. This is exactly what we observe in SE, where the red-shifted absorption components extend to velocities of +380 km s−1, which is the free-fall velocity at the stellar surface given the stellar parameters. We note that the red-shifted edge of the absorption line has the same velocity in all our spectra taken over the period studied. Furthermore, with the v sin i:s presented in Table 1 and inclinations of ~60°, the equatorial rotational velocities of the stars would be ~13 km s−1. With stellar radii of R* = 1.4 R⊙ the periods of rotation are 5.4 days, which is typical for young TTSs. This is also congruent with the result presented in Petrov et al. (2011), where it was demonstrated that narrow chromospheric emission lines and photospheric lines show anti-phase variations in velocity caused by an offset between the emitting region and the pole of rotation. From the observed fluctuations, one would infer a rotational period of 5 to 6 days, although the observations cover only part of the cycle.

An additional argument in favour of high inclinations comes from the periodic radial velocity variations observed for the narrow He I emission lines at 4471, 5875, and 6678 Å. This emission originates in hot regions which are off-set from the pole. In Fig. 11, we have depicted the most probable periods with sinusoids for each star. The resulting periods, 4.2 ± 1 days for NW and P = 5.5 ± 1 days for SE, are congruent with the results obtained above. However, the observations span only part of the periods in question, and further observations are warranted. In the case of low inclinations, for example, i = 30° deg, the expected period of rotation would be less than 3 days.

There are, however, some puzzling circumstances related to this scenario. Firstly, our models suggest that the inclinations of the stars are ~65°, and that therefore the disks are inclined to the line-of-sight by ~25°. It follows that the occulting circumstellar grains must be located either very close to the star or high above the disk plane if at larger distances. In the first case, one may consider that dust fragments are stirred up at the innermost part of the disk, something that appears to occur in RW Aur A according to Petrov et al. (2015). This star undergoes large dips in brightness due to extinction changes and it has very strong red-shifted absorption components like those observed in SE. However, in RW Aur A the drops in brightness are “grey”, that is, the dust grains are large, as expected from their origin in the disk plane. Another possibility, in the case of SE and NW, is that dust is carried fromthe inner disk edge along the accreting gas, and that dust grains sputter and become smaller on their way.

Secondly, from the sub-au resolution NIR maps of SE and NW, Garcia et al. (2017) concluded that the disks of both stars are seen nearly face on; the corresponding inclinations were estimated as 22° ± 6° and 28° ± 3°. In line with this, Pontoppidan et al. (2011), Schegerer et al. (2009), and Vural et al. (2012) also reported relatively low inclination angles. At these inclinations we cannot find any reasonable fit between modelled and observed line profiles. The estimates are subject to uncertainties. In our calculations, a change of ±10° would still give reasonable fits to the observed profiles. Furthermore, inclination angles defined from modelling of hydrogen lines and molecular CO lines can differ (e.g. Grady et al. 2010). Estimates of disk inclinations based on different methods and extracted from different wavelength regions can be substantially different (see e.g. Eisner et al. 2004 and Doucet et al. 2006).

The reason for this discrepancy remains unclear. Perturbed disks developing open spiral patterns may take forms that differ from the circular case, but it seems unlikely that the disks in both the SE and NW components would be affected in the same way. Interferometric NIR observations are sensitive to the thermal radiation from a disk wind, whichaccording to Bans & Königl (2012) can give a large contribution to the NIR excess of TTSs. This additional source, if present, is not as flat as the disk and would result in lower limits of the disk inclinations in both components. Observations using the sub-millimetre interferometer ALMA can shed more light on the precise forms of the disks.

Alternatively, the concept of standard magnetospheric accretion depicted in Fig. 7 must be revised. To match both the low inclinations and observed accretion signatures, one might consider more complex configurations where the magnetosphere is tilted relative to the disk plane or multipole magnetic structures.

In the case where the disks are oriented nearly face-on, there is a problem in accounting for the results obtained from the rotational periods and the variable circumstellar extinction. In this case, the brightness changes might be caused by dust carried along with the disk wind or caused by dusty cloudlets containing small grains passing the stars, possibly in remote remnants of the parent cloud surrounding the binary, explanations one might consider also in our configuration.

Nevertheless, from the present study we find that with inclinations of ~65° we not only get a satisfactory fit between modelled and observed Balmer line profiles, but also an agreement between the resulting rotational periods and those inferred from the anti-phase changes described above. In addition, we have found an explanation of the extent of the red-shifted absorption components in SE. These features are weak in NW, while the strength of the Hα emission indicates a very high accretion rate. The absence of inverse P Cygni absorptions does not imply a low accretion rate however, since their strength also depends on filling factors, the inner radius of the accretion disk, and the orientation and structure of the magnetosphere.

thumbnail Fig. 11

Time variations of radial velocity of HeI emission lines. Upper panel: S CrA NW, lower panel: S CrA SE. Dashed lines: approximation by sinusoid with periods of 4.2 days (NW) and 5.5 days (SE).

6 Conclusions

We have investigated the properties of the two stars forming the visual T Tauri binary S CrA based on the results from spectral and photometric monitoring and theoretical modelling of their magnetospheric structures and winds. Both stars show large brightness variations (≥2 mag in V band), an unusually high and variable degree of veiling, remarkable emission line spectra with both broad and narrow emission lines, and broad red-shifted absorption lines flanking strong lines, properties that are rare among the T Tauri stars. The stars are labelled NW and SE, and this study leads to the following conclusions:

We have derived a number of stellar properties and found that the two stars are remarkably similar in all relevant parameters: Teff (K) = 4250 K, log g = 4.0, v sin i = 12 km s−1, M* = 1 M⊙, L* = 0.7 L⊙, and R* = 1.4 R⊙, valid for both stars. Hence, S CrA NW and SW represents a T Tauri twin. However, their emission spectra are different. Both stars have narrowemission lines superimposed on broad emission lines, and the latter are strongest in NW. This star also has the highest degree of veiling reaching close to VF = 10 on occasion, in terms of the standard veiling factor. As demonstrated in previous studies the contribution from narrow emission lines filling in the photospheric lines dominates in VF, especially when the veiling is high. Therefore, variations in the continuous excess emission contribute only to a smaller extent in the light changes. The narrow emission lines are formed in surface areas off-set from the pole and vary in radial velocity in anti-phase with the photospheric lines. From these fluctuations we infer rotational periods of ≈4.2 days for NW and ≈5.5 days for SE, and combined with the measured v sin i:s we obtain inclinations of the stars of around 60°.

Both stars have strong emission lines of neutral and ionized elements with absorption components related to both accretion flows and winds. The red-shifted absorption components are strongest in SE, extending to +380 km s−1, the estimated free-fall velocity at the surface.

In our model simulations of the gas emission in NW and SE we took into account magnetospheric accretion, a polar wind, and a disk wind. We compared the resulting Balmer line profiles with those observed and found that very good agreements can be obtained for inclinations of ~65° in a hybrid model, including magnetospheric accretion and a disk wind only. From our models we derive mass-loss rates in winds ranging from 5 × 10−9 to 5 × 10−8 M⊙ yr−1 and mass accretion rates of ~5 × 10−7. The accretion rates appear to be larger in NW by a factor of approximately two judging from spectral signatures.

For inclinations of ~65° the disks are inclined by ~25° to the line-of-sight. At this orientation the trajectories of the accreting funnels just above the stellar surface are parallel to the line-of-sight in agreement with our observation of gas moving towards SE at free-fall velocity. Moreover, the inclinations inferred from the anti-phase velocity variations mentioned above agree with those adopted in the model.

The brightness changes are caused mainly by rapid variations in the amount of dust in the line-of-sight. The stars become redder with decreasing brightness, indicating that the foreground occulting dust grains are small, and it is unlikely that this dust has been stirred up from the inner disk. We discuss whether sputtering of larger grains carried along the accretion funnels can produce the observed fluctuations.

Another open issue is that the inclinations derived from our observations and model calculations are different from those obtained from interferometric NIR spectroscopy, with disks oriented almost face-on, and we discuss possible causes for this puzzling discrepancy. In the case where the disks are seen nearly face-on, it is more likely that the occulting dust is carried by the disk wind or resides in dusty cloudlets further out from the system. Future observations with ALMA can shed more light on the shapes and orientations of the disks in question.

Acknowledgements

This work was supported by the INTAS grant 03-51-6311, the Swedish National Space Board, and the Magnus Bergvall foundation. FMW thanks past Provosts and Vice Presidents for Research of Stony Brook University for providing support for the SMARTS partnership.

Appendix A Hydrogen line profiles in models assuming magnetospheric accretion and a disk wind

Figure A.1 shows results of calculations of the Balmer line profiles in models of the magnetospheric accretion and in the hybrid model combining the magnetospheric accretion and a disk wind (models MA1 and DW1). The line profiles are shown for inclinations from 10 to 80 degrees.

thumbnail Fig. A.1

Left: Hα, Hβ and Hγ line profiles resulting from the magnetospheric accretion region (model MA1); Right: as in the left panel but in the framework of the hybrid model combining magnetospheric accretion (MA1) and a disk wind (DW1). Inclination angles are marked.

References

  1. Appenzeller, I., Jankovics, I., & Jetter, R. 1986, A&AS, 64, 65 [Google Scholar]
  2. Bans, A., & Königl, A. 2012, ApJ, 758, 100 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bast, J. E., Brown, J. M., Herczeg, G. J., van Dishoeck, E. F., & Pontoppidan, K. M. 2011, A&A, 527, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Beristain, G., Edwards, S., & Kwan, J. 1998, ApJ, 499, 828 [NASA ADS] [CrossRef] [Google Scholar]
  5. Blandford, L. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bouvier, J., Alencar, S. H. P., Harries, T. J., et al. 2007, in Protostars and Planets V, eds. B. Reipurth, D. Jewitt, and K. Keil (Tucson: University of Arizona Press), 479 [Google Scholar]
  7. Calvet, N., & Hartmann, L. 1992, ApJ, 386, 239 [NASA ADS] [CrossRef] [Google Scholar]
  8. Carmona, A., van den Ancker, M. E., & Henning, Th. 2007, A&A, 464, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Cranmer, S. R. 2008, ApJ, 689, 316 [NASA ADS] [CrossRef] [Google Scholar]
  10. Dodin, A. V. 2018, MNRAS, 475, 4367 [NASA ADS] [CrossRef] [Google Scholar]
  11. Dodin, A. V., & Lamzin, S. A. 2012, Astron. Lett., 39, 389 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dodin, A.V., Lamzin, S. A., & Sitnova, T. M. 2013, Astron. Lett., 39, 315 [NASA ADS] [CrossRef] [Google Scholar]
  13. Doucet, C., Pantin, E., Lagage, P. O., & Dullemond, C. P. 2006, A&A, 460, 117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Edwards, S. 1979, PASP, 91, 329 [CrossRef] [Google Scholar]
  15. Eisner, J. A., Lane, B. F., Hillenbrand, L. A., Akeson, R. L. & Sargent, A. I. 2004, ApJ, 613, 1049 [NASA ADS] [CrossRef] [Google Scholar]
  16. Fang, M., van Boekel, R., Wang, W., et al. 2009, A&A, 504, 461 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Ferreira, J. 2013, EAS Pub. Ser., 62, 169 [CrossRef] [Google Scholar]
  18. Ferreira, J., Dougados, C., & Cabrit, S. 2006, A&A, 453, 785 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Gahm, G. F., Walter, F. M., Stempels, H. C., Petrov, P. P., & Herczeg, G. J. 2008, A&A, 482, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Garcia, P. J. V., Ferreira, J., Cabrit, S., & Binette, L. 2001, A&A, 377, 589 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Garcia Lopez, R., Perraut, K., Caratti o Garatti, A., et al. 2017, A&A, 608, A78 [Google Scholar]
  22. Grady, C. A., Hamaguchi, K., Schneider, G., et al. 2010, ApJ, 719, 1565 [NASA ADS] [CrossRef] [Google Scholar]
  23. Graham, J. A. 1992, PASP, 104, 479 [NASA ADS] [CrossRef] [Google Scholar]
  24. Grinin, V. P., & Tambovsteva, L. V. 2011, Astron. Rep., 55, 704 [NASA ADS] [CrossRef] [Google Scholar]
  25. Grinin, V., Tambovtseva, L., Potravnov, I., & Mkrtichian, D. 2017, ASP Conf. Ser., 508, 185 [Google Scholar]
  26. Guo, Z., Herczeg, G. J., Jose, J., et al. 2018, ApJ, 852, 56 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hamann, F., & Persson, S.E. 1992, ApJS, 82, 247 [NASA ADS] [CrossRef] [Google Scholar]
  28. Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135 [NASA ADS] [CrossRef] [Google Scholar]
  29. Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 1994, AJ, 108, 1906 [NASA ADS] [CrossRef] [Google Scholar]
  30. Kessler-Silacci, J., Augereau, J.-C., Dullemond, C. P., et al. 2006, ApJ, 639, 275 [NASA ADS] [CrossRef] [Google Scholar]
  31. Königl, A., & Pudritz, R. E. 2000, in Protostars and Planets IV, eds. Mannings, V., Boss, A. P., & Russell, S. S., 759 [Google Scholar]
  32. Krautter, J., Appenzeller, I., & Jankovics, I. 1990, A&A, 236, 416 [NASA ADS] [Google Scholar]
  33. Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W. 1999, A&AS, 138, 119 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
  34. Kurucz, R. L. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid, CD-ROM No. 13 (Cambridge, MA: Smithsonian Astrophysical Observatory) [Google Scholar]
  35. Kurosowa, R., Harries, T. J., & Symington, N. H. 2006, MNRAS, 370, 580 [NASA ADS] [CrossRef] [Google Scholar]
  36. Kurosawa, R., Romanova, M. M., & Harries, T. J. 2011, MNRAS, 416, 2623 [NASA ADS] [CrossRef] [Google Scholar]
  37. Kwan, J., & Fischer, W. 2011, MNRAS, 411, 2383 [NASA ADS] [CrossRef] [Google Scholar]
  38. Matt, S, & Pudritz, R. E. 2007, in Proceedings of the International Astronomical Union, eds. J. Bouvier & I. Appenzeller, IAU Symp., 243, 299 [CrossRef] [Google Scholar]
  39. Meeus, G., Salyk, C., Bruderer, S., et al. 2013, A&A, 559, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Petrov, P. P. 2003, Astrophysics, 46, 4 [CrossRef] [Google Scholar]
  41. Petrov, P. P., Gahm, G. F., Stempels, H. C., Walter, F. M., & Artemenko, S. A. 2011, A&A, 535, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Petrov, P. P., Gahm, G. F., Herczeg, G. J., Stempels, H. C., & Walter, F. M. 2014, A&A, 568, L10 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Petrov, P. P., Gahm, G. F., Djupvik, A. A., et al. 2015, A&A, 577, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Piskunov, N., & Valenti, J. A. 2017, A&A, 597, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525 [NASA ADS] [Google Scholar]
  46. Pontoppidan, K. M., Blake, G. A., & Smette, A. 2011, ApJ, 733, 84 [NASA ADS] [CrossRef] [Google Scholar]
  47. Prato, L., Greene, T. P., & Simon, M. 2003, ApJ 584, 853 [NASA ADS] [CrossRef] [Google Scholar]
  48. Pudritz, R. E. & Norman, C. A. 1986, ApJ, 301, 571 [NASA ADS] [CrossRef] [Google Scholar]
  49. Rei, A. C. S., Petrov, P. P., & Gameiro, J. F. 2018, A&A, 610, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Reipurth, B., & Graham, J. A. 1988, A&A, 202, 219 [NASA ADS] [Google Scholar]
  51. Reipurth, B., & Zinnecker, H. 1993, A&A, 278, 81 [NASA ADS] [Google Scholar]
  52. Reipurth, B., Pedrosa, A., & Lago, M. T. V. T. 1996, A&AS, 120, 229 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Rydgren, A. E. 1977, PASP, 89, 557 [CrossRef] [Google Scholar]
  54. Safier, P. N. 1993, ApJ, 408, 115 [NASA ADS] [CrossRef] [Google Scholar]
  55. Schegerer, A. A., Wolf, S., Hummel, C. A., Quanz, S. P., & Richichi, A. 2009, A&A, 502, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [Google Scholar]
  57. Sobolev, V. V. 1960, in Moving Envelopes of Stars (Cambridge: Harvard University Press) [CrossRef] [Google Scholar]
  58. Stempels, H. C., Gahm, G. F., & Petrov, P. P. 2007, A&A, 461, 253 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Takami, M., Bailey, J., & Chrysostomou, A. 2003, A&A, 397, 675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Tambovsteva, L. V., Grinin, V. P., Rodgers, B., & Kozlova, O. V. 2001, Astron. Rep., 45, 442 [NASA ADS] [CrossRef] [Google Scholar]
  61. Tambovsteva, L. V., Grinin, V. P., & Weigelt, G.. 2016, A&A, 590, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Vural, J., Kreplin, A., Kraus, S., & et al. 2012, A&A, 543, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Walter, F. M., & Miner J. 2005, In Proc. 13th Cool Stars Workshop, Hamburg, July 2004, eds. F. Favata, G. Hussain, & B. Battrick, ESA SP, 560, 1021 [Google Scholar]
  65. Walker, M. F. 1972, ApJ, 175, 89 [NASA ADS] [CrossRef] [Google Scholar]
  66. Wang, H., Mundt, R., Henning, T., & Apai, D. 2004, ApJ, 617, 1191 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Stellar parameters.

In the text
Table 2

Ranges of model parameters.

In the text
Table 3

Parameters of the disk wind model.

In the text
Table 4

Parameters of the magnetospheric accretion model.

In the text

All Figures

thumbnail Fig. 1

Upper panel: colour-magnitude diagram from the long-term monitoring of S CrA (blue: S CrA SE, red: S CrA NW). Lowerpanel: individual I and B magnitudes observed over the period of spectral observations.
In the text
thumbnail Fig. 2

Upper panel: spectra of the stars obtained during the 4 nights over a spectral region relatively free from emission lines and veiling used to derive stellar parameters (time running upwards). Lower panel: good agreements between observed (black) and synthetic (grey) spectra. Different contributions from a continuous source of veiling is added (horizontal dashed lines).
In the text
thumbnail Fig. 3

Segments of spectra of S CrA NW (upper panel) with broad emission lines and S CrA SE (lower panel) with strong red-shifted absorption components. Narrow emission lines of metals are present in both stars.
In the text
thumbnail Fig. 4

Spectral segments of the two stars taken during 4 nights over the Ca II H & K (upper panel) and Hγ and He I plus Na I D lines (lower panels). Intensities are normalized to the continuum level and time is running upwards.
In the text
thumbnail Fig. 5

Ca II K and Na I D lines in SCrA NW show strong blue-shifted absorption components in addition to narrow interstellar lines (spectra from August 14). Velocities are astro-centric in this and similar figures.
In the text
thumbnail Fig. 6

Balmer line profiles. Left: S CrA NW (August 14). Right: S CrA SE (August 16). From top to bottom Hα, Hβ, Hγ, and Hδ with the continuum level marked with dashed lines.
In the text
thumbnail Fig. 7

Sketch of the hybrid model geometry assumed in the present calculations (not to scale). The inclination angle is measured from the symmetry axis (S).
In the text
thumbnail Fig. 8

Components of the Hα line profiles calculated in model MA1 + DW1 at 30° (left) and 65° (right): magnetospheric accretion (red), disk wind (blue), and the resulting line profile (black).
In the text
thumbnail Fig. 9

Balmer line profiles in S CrA SE. Upper panels: Hα, Hβ and Hγ line profiles as calculated from the hybrid model MA1 plus DW1 and with an inclination angle of 65°. Lower panels: same profiles normalized to peak intensity (solid curves) superimposed on the observed profiles (dashed curves).
In the text
thumbnail Fig. 10

As in Fig. 3 but S CrA NW. The line profiles are calculated in the hybrid model MA2 plus DW2 and with an inclination angle of 60°.
In the text
thumbnail Fig. 11

Time variations of radial velocity of HeI emission lines. Upper panel: S CrA NW, lower panel: S CrA SE. Dashed lines: approximation by sinusoid with periods of 4.2 days (NW) and 5.5 days (SE).
In the text
thumbnail Fig. A.1

Left: Hα, Hβ and Hγ line profiles resulting from the magnetospheric accretion region (model MA1); Right: as in the left panel but in the framework of the hybrid model combining magnetospheric accretion (MA1) and a disk wind (DW1). Inclination angles are marked.
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.