© ESO, 2017

1. Introduction

Geological evidence suggests that the early Earth had a warmer climate in the first few 100 Myr of its evolution. Such a mild and warm climate on the early Earth 4 Gyr ago was necessary and essential for the evolution and formation of life on our planet (Kasting & Toon 1989; Sackmann & Boothroyd 2003). However, solar standard models predict a lower bolometric luminosity of the Sun at that time, being just 70% the present-day luminosity. The evolution of the Sun’s luminosity had an important effect on the formation of the atmosphere for our Earth and for early Mars. Without an atmosphere on Earth, the average surface temperature would have been 235 K only. Additional present-day greenhouse gases would have raised the temperature to ~253 K, which is still not enough to avoid the completely frozen surfaces on early Earth and Mars (Sagan & Mullen 1972; Kasting & Catling 2003). The discrepancy between the implications from the solar standard models and the geological evidence for a warmer climate on Earth is defined as the “faint young Sun paradox” (FYSP). Apart from a number of proposed solutions of the FYSP (see e.g. Gaidos et al. 2000; Feulner 2012), an astrophysical solution for this problem has been suggested. It assumes that the young main-sequence Sun was brighter than suggested by the standard model, which would be possible if it had been more massive than today and consequently suffered from an increased mass loss during its early main-sequence life through an enhanced solar wind (Graedel et al. 1991; Gaidos et al. 2000).

Winds play an important role in stellar evolution for main-sequence stars like the Sun, especially for the stellar angular momentum. We know that stars spin down with age, because angular momentum is carried away by the magnetized, ionized winds. To understand the mechanism of the interaction between the stellar wind, stellar rotation, and the magnetic field for stars with various ages, information on how winds evolve with time is required. Furthermore, the evolution of stellar winds is important for the evolution of planetary atmospheres and their erosion (Lammer et al. 2010). Most of what we know about stellar winds comes from studies of the solar wind, although the mechanisms for generating, accelerating, and heating the solar wind are still poorly understood (Cranmer 2009; McComas et al. 2003; Schwenn 2006).

Today, the most common way to assess stellar winds and therefore, determine stellar mass loss rates, is to observe the Lyman-α excess of the neutral interstellar hydrogen in high-resolution Hubble Space Telescope spectra of stars, as introduced by Wood et al. (2002, 2005a), Wood (2004). Due to charge exchange interactions between neutral interstellar hydrogen and the ionized wind, a “wall” of hot neutral hydrogen at the edge of the stellar astrosphere is built up. The detected material is not from the fully ionized wind itself, which has no HI, but is interstellar HI instead that is heated within the interaction region between the wind and the local interstellar medium. The amount of astrospheric HI absorption provides diagnostic information on the rate of mass loss of the wind, specifically the momentum in the wind. In Wood et al. (2002) the correlation of the mass loss rates with coronal properties was studied. The coronal X-ray luminosity is a good indicator for the magnetic activity of a star and the scaling relationship between the mass loss rate per unit surface area and the X-ray surface flux is $\mathit{Ṁ}\mathrm{\propto }{\mathit{F}}_{\mathrm{X}}^{\mathrm{1.34}\mathrm{\pm }\mathrm{0.18}}\mathit{,}$(1)which, combined with X-ray luminosity evolution versus time L_X ∝ t^-1.5 (Güdel et al. 1997), suggests that the mass loss rate decreases with time for solar-like stars like Ṁ ∝ t^{− 2.33 ± 0.55} (Wood 2004). The correlation is, however, still not sufficient to solve the FYSP (Wood et al. 2002; Minton & Malhotra 2007). However, as these authors concluded in their study, this correlation fails for the youngest and most active stars for which winds appear to be very weak in Lyman-α observations. Wood et al. (2005b) observed χ¹ Ori, but they were unable to provide any astrospheric detections in their study. They argued that a non-detection does not generally provide a meaningful upper limit to the stellar wind strength, because for non-detections the star could be surrounded by a fully ionized interstellar medium (ISM; Wood et al. 2005b).

Observing and detecting stellar winds similar to the solar case is important to improve our understanding of stellar evolution, such as the correlation between rotation and stellar magnetic activity which provides information on the dynamo and thus magnetic activity. Furthermore, the understanding of acceleration mechanisms of these winds could be improved and the measurements of wind properties of stars with different ages may provide essential information on stellar angular momentum loss. Radio observations of young, solar-type stars are used in our study to test if there was a strong mass loss in the young Sun. A study of the “radio Sun in time”, complementing the “X-ray Sun in time” (Güdel et al. 1997), can explore the range and the long-term evolution of solar and stellar magnetic activity and wind mass loss. First detections of radio emission from low-mass main-sequence stars were reported by Gary & Linsky (1981) and Linsky & Gary (1983). Limits to mass loss have already been established from radio observations of more massive A and F stars (Brown et al. 1990) and active, less massive M stars (Lim & White 1996). Scuderi et al. (1998) observed early type O and B supergiants to make a detailed comparative study of the mass loss evaluated from Hα and radio continuum observations. Güdel et al. (1998) and Gaidos et al. (2000) used the Very Large Array (VLA) to search for radio emission of the active, young, solar-type stars π¹ UMa, κ¹ Cet and β Com at 8.4 GHz. Their observations resulted in 3σ detection limits of 20–30 μJy, which correspond to radio luminosities of ~10^12.5 erg s^-1 Hz^-1 (Villadsen et al. 2014). Early radio observations of EK Dra were recorded in Güdel et al. (1995), where at minimum the 8.4 GHz flux was (34 ± 11) μJy, and at intermediate levels (77 ± 9) μJy.

To derive an estimate or upper limit for the enhanced young solar wind, we observed radio emission of young, solar-like analogues at the main-sequence with the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/Submillimeter Array (ALMA). If we are able to detect free-free radio emission of such winds, their mass loss rates can be calculated. From climate predictions the initial (zero-age main-sequence, ZAMS) solar mass is required to be in the range of 1.03–1.07 M_⊙ if it were to solve the FYSP (Sackmann & Boothroyd 2003; Whitmire et al. 1995), suggesting an enhanced early wind mass loss of the order of 10^-12−10^-10M_⊙ yr^-1. In comparison, the present-day solar wind mass loss amounts to 2 × 10^-14M_⊙ yr^-1 (Feldman et al. 1977).

In this paper we focus on the four young solar analogues EK Dra, π¹ UMa, χ¹ Ori, and κ¹ Cet using the upgraded sensitivity and resolution of the VLA. In Sect. 2, we briefly describe the observations including a description of our targets. Section 3 contains the results of our detections and upper limits of radio emission. The calculation of the mass loss rates of our star sample will be described in Sect. 4, where we will also compare our observational results to results from Lyman-α absorption, presented by Wood et al. (2002, 2005a).

2. Observations

2.1. Target sample

Our target sample includes the following objects (see also Table 1 summarizing the properties of our stars):

• EK Dra:

  this is a G1.5 V star that is considered to be among the most active solar analogues in our neighbourhood, with a distance of 34 pc from the Sun. The average rotation period is 2.68 days. Main properties are reviewed by Strassmeier & Rice (1998) and Messina & Guinan (2003). Ribas et al. (2005) adopted an age of about 100 Myr for this near-ZAMS star;

• π¹UMa:

  this is a young, active G1.5 V solar proxy with a rotation period of about 4.9 days (Messina & Guinan 2003) and a distance of 14.3 pc. In the Sun in Time programme, π¹ UMa is reported to have an age of 300 Myr (Ribas et al. 2005);

• χ¹Ori:

  a G1V star with a rotation period of about 5.2 days (Messina et al. 2001), a distance of 8.7 pc and an age of 300 Myr (Ribas et al. 2005). The star χ¹ Ori is classified as a member of the Ursa Major moving group (King et al. 2003);

• κ¹Cet:

  with a spectral type G5 V, it is the coolest star in the sample, with a distance of 9.2 pc from the Sun. Gaidos & Gonzalez (2002) determined spectroscopic parameters. The rotation period is reported by Messina & Guinan (2003) to be about 9.2 days and the age is suggested to be around 650 Myr (Ribas et al. 2005, 2010).

2.2. VLA and ALMA

For the radio measurements we use the Karl G. Jansky VLA, a radio interferometer located in New Mexico near Socorro, operated by the National Radio Astronomy Observatory (NRAO). We use C-band (4–8 GHz, 6 cm) and Ku-band (12–18 GHz, 2 cm) receivers. The Jansky VLA operates with an increased sensitivity relative to the VLA. The observations were performed in C configuration in sessions in spring/summer 2012 and 2013. The Common Astronomy Software Application (CASA) developed by the NRAO has been used for inspecting, editing (including flagging), calibrating, and imaging the data sets. Flux calibrators were observed at the beginning of each observation for several minutes and the phase calibrators were repeatedly observed together with the targets. An overview and summary of the observations is given in Table 2. For the calibration, the raw data needs to be inspected first, which means that bad data due to antenna errors, shadowed antennas, or poor weather conditions need to be flagged, that is removed from the data set. Afterwards, flux, bandpass, and gain calibration steps are applied. We used a pipeline for VLA data¹ that deals with the flagging and calibration. We used this pipeline, but additional flagging was necessary afterwards.

ALMA, located in the Chajnantor plain of the Chilean Andes, was used to observe in band 3 (with a bandwidth of 84–116 GHz) at 100 GHz in December 2013 within Cycle 1. We got observing time for one of our targets, χ¹ Ori. For the ALMA data the NRAO staff provided prefabricated scripts together with our data for flagging and calibration. In the meantime, a calibration pipeline for ALMA has been developed as well². For our data analysis, we used these pipelines, but some extra flagging and a second run through the pipeline calibration were necessary.

From NRAO’s exposure calculator for the VLA, the theoretical noise sensitivity with 2 GHz bandwidth, 27 antennas, and one hour on source is calculated to be around 3.5 μJy rms in C-band and around 3.8 μJy rms in Ku-band. These values represent the expected random noise levels for π¹ UMa and κ¹ Cet, respectively. Except for π¹ UMa in C-band, the achieved noise levels are in good agreement with the expected values (see Table 2). For EK Dra with a bandwidth of 3.5 GHz and 26 antennas, we would expect a noise level of 2.2 μJy in C-band and 2.9 μJy in Ku-band, whereas the achieved values are slightly higher. For χ¹ Ori, 3.5 GHz bandwidth and 26 antennas, the theoretical rms is 1.6 μJy in C-band, which is in good agreement with the observational noise. In Ku-band the expected noise is around 2.2 μJy, whereas the achieved rms is lower, namely 1.6 μJy. The CLEAN procedure is applied to produce images, where the Clark algorithm with natural weighting was chosen for the setting. Depending on the wavelength and the number of antennas, a cell size of 0.7′′ and 0.3′′ for C and Ku-band was used, respectively.

3. Results

3.1. Images

For each target several observation sets are available. To obtain the final images of each target, all data sets in each frequency band are combined to concatenated images which are shown for the detections in Figs. 1 and 2. The crosses mark the expected positions of the sources predicted from the SIMBAD Astronomical Database³ corrected for proper motion to the epoch of observation. Minor offsets in right ascension and declination for EK Dra and χ¹ Ori occur in our analysis. We note that the field of view in the Ku-band images is much smaller than the C-band images. Several sources can be identified in all images, but only two of the four objects of the sample show a radio detection signal at the expected positions. The targets EK Dra and χ¹ Ori, shown in Figs. 1 and 2, are detected in Stokes I (total intensity), with a total flux of around 100 μJy. On the other hand, π¹ UMa and κ¹ Cet (not shown) display non-detections at the expected target positions both in C and Ku-band.

Fig. 1 Contours of EK Dra in a) C-band with contour levels [0.2, 0.4, 0.6, 0.8] · 576 μJy; and b) Ku-band with [0.2, 0.4, 0.6, 0.8] · 74.5 μJy and the beam size in the left corner. The point in the centre gives the coordinates from the Gaussian fit, the cross marks the proper motion corrected Hipparcos position.

Fig. 2 χ1 Ori in a) C-band with contour levels [0.2, 0.4, 0.6, 0.8] · 223 μJy; b) Ku-band with [0.2, 0.4, 0.6, 0.8] · 106.6 μJy; and c) for ALMA at 100 GHz with [0.2, 0.4, 0.6, 0.8] · 108 μJy. The black dot represents the coordinates calculated from the Gaussian fit. The red cross marks the position of the companion χ1 Ori B and the black cross the corrected coordinates of the main target χ1 Ori.

3.2. Radio emission from the stellar wind

The free-free spectrum from thermal bremsstrahlung radiation is characterized by a power-law spectral index α, ranging from −0.1 ≤ α ≤ 2, where the flux density is given by S_ν ∝ ν^α at the frequency ν. The integrated flux densities of the detections are determined by fitting the stellar images by a Gaussian profile. The associated rms values in Stokes I and V (circularly polarized intensity) in both wavelength bands are given in Table 3. For those objects for which no detection was observed, the 3σ upper limit to the flux density from a source-free background region is estimated. Time series of the sources provide information on the variation of each observation interval. For each object, separate images for each time interval of about five minutes in the Stokes I/V plane are created and hence, the time-dependent flux and the related rms are extracted. Because the observing time is much shorter than the stellar rotation periods, no rotational modulation should be seen. On the other hand, short time variations of a few minutes are indicators for flares.

The results for the four targets are summarized in the following sections.

3.2.1. EK Dra

We obtained a clear detection of EK Dra at ${\mathrm{14}}^{\mathrm{h}}{\mathrm{38}}^{\mathrm{m}}\mathrm{59}\stackrel{\mathrm{s}}{.}\mathrm{96}\mathit{,}\mathrm{+}{\mathrm{64}}^{\mathrm{◦}}{\mathrm{17}}^{\mathrm{\prime }}\mathrm{29}.̋\mathrm{49}$. The offsets from the predicted positions in C-band are 0.02s in right ascension and –0.19′′ in declination, and 0.04s and 0.84′′ in Ku-band, well within the beam size of and , respectively. Because EK Dra is a very active star we expect that the radio emission will include coronal emission (Güdel et al. 1995). The Stokes I radio flux was 593 ± 1.7 μJy with an rms of 3.4 μJy in C-band. Judging from the light curve, no flare event seems to be present. In Ku-band the radio emission is observed at 73 ± 2.4 μJy with an rms of 4 μJy. The star EK Dra’s radio emission cannot be only thermal free-free emission as also argued by Güdel et al. (1995), given the variability and the high flux level. The polarization degree r_c = V/I, which ranges from –1 to 1, is found to vary in the range r_c = [–0.088, –0.015] in C-band. In Ku-band our observation does not show any significant non-zero Stokes V flux.

3.2.2. π¹ UMa

The star π¹ UMa, expected at ${\mathrm{08}}^{\mathrm{h}}{\mathrm{39}}^{\mathrm{m}}\mathrm{11}\stackrel{\mathrm{s}}{.}\mathrm{65}\mathit{,}\mathrm{+}{\mathrm{65}}^{\mathrm{◦}}{\mathrm{01}}^{\mathrm{\prime }}\mathrm{16}.̋\mathrm{46}$, is a non-detection and was already studied by other authors (e.g. Gaidos et al. 2000). The 3σ upper limits of the integrated radio intensities are 23.1 μJy in C-band and 6.3 μJy in Ku-band. The C-band intensity limit is high compared to the Ku-band results because, despite heavy flagging and cleaning, the residual of a strong source strongly perturbs our object region and consequently raises the rms. During the observation the fringe pattern directly crossed the expected position of π¹ UMa and caused an increase in the background noise and hence negatively influenced the radio emission estimation for π¹ UMa. Therefore, the radio flux density upper limit in C-band is not as useful as desired. On the other hand, the observations of the Ku-band flux density upper limit of 6.3 μJy are excellent and useful for further analysis and interpretation. The polarization map also shows only noise. Gaidos et al. (2000) reported a non-detection at the location of π¹ UMa as well. They placed a 2σ upper limit of 12 μJy at 3.6 cm (X-band) for the total flux density. Our VLA observations lower these upper limits by a factor of around two.

3.2.3. χ¹ Ori

The star χ¹ Ori is located at ${\mathrm{05}}^{\mathrm{h}}{\mathrm{54}}^{\mathrm{m}}\mathrm{22}\stackrel{\mathrm{s}}{.}\mathrm{78}\mathit{,}\mathrm{+}{\mathrm{20}}^{\mathrm{◦}}{\mathrm{16}}^{\mathrm{\prime }}\mathrm{33}.̋\mathrm{58}$ in our observations. It shows strong radio emission, seen with offsets in C-band of –0.06s in right ascension and 0.62′′ in declination, relative to the expected position (cross in Fig. 2a) using Hipparcos⁴ measurements (van Leeuwen 2007). In Ku-band the offsets to the observational positions are –0.07s in right ascension and 0.61′′ in declination. The integrated radio flux densities in Stokes I as given in Table 3 are 110 ± 0.7μJy with an rms of 1.8 μJy in C-band and 117 ± 2.7 μJy and a corresponding rms of 1.6 μJy in Ku-band. The flux density at 100 GHz measured with ALMA is 103 ± 4.9 μJy in Stokes I. The proper motion corrected offset in right ascension is around –0.03s, in declination it is 0.29′′. The C-band light curve shows a flare that can be clearly identified during an observation interval with a duration of less than 30 min (see Fig. 3). The peak reaches a flux density about three times the quiescent level. The occurrence of the flare and the fact that the slope of the spectrum is slightly negative with increasing frequency, suggest that the radio emission of χ¹ Ori is not exclusively thermal radio bremsstrahlung from a wind but is dominated by gyrosynchrotron emission from accelerated electrons. The third indication supporting this assumption is a ≈10% Stokes V signal (see Table 3). The degree of circular polarization r_c is in the range r_c = [–0.35, 0.63] with a maximum sigma σ = 0.18 in C-band and r_c = [–0.68, 0.72] with σ = 0.49 in Ku-band. For ALMA, no Stokes V measurements were available for Cycle 1 data sets.

Fig. 3 Time series of χ1 Ori in C-band showing an explosive increase in intensity identified as a flare event at the central frequency of 6.15 GHz. The time interval between the data points is around ten minutes.

The images of χ¹ Ori show that the corrected coordinates (black crosses in Fig. 2) do not properly match with the observational positions from the Gaussian fit (black dots). Therefore, we analyzed if the radio signal may come from the M-dwarf companion of χ¹ Ori (Han & Gatewood 2002; König et al. 2002). To derive the position of χ¹ Ori B, the orbit of χ¹ Ori has to be corrected first. The orbital parameters are taken from Han & Gatewood (2002) and König et al. (2002). By correcting the orbit from JD 1991.25 to JD 2012.4 when our VLA observations took place and by including the correction for proper motion from Hipparcos, the expected coordinates for χ¹ Ori are derived. The position of the companion is determined by using the mass ratio between primary and companion, and is displayed by the red cross in the images of Fig. 2 and listed in Table 4 (in C-band only). The two components are separated by 0.49′′ from each other. Some systematic errors occur from Hipparcos itself, especially because Hipparcos did not recognize the binarity of χ¹ Ori, and errors in proper motion and possible position errors of the phase calibrator during the observation may contribute to the residual deviation of the detected coordinates. We checked for new Gaia position measurements⁵ but unfortunately there is no data available for χ¹ Ori. If the companion is responsible for the radio emission, which seems likely, we will still use the observational radio emission of χ¹ Ori A or B for our further analysis and mass loss rate calculations considering it to be an upper limit to the thermal wind emission.

3.2.4. κ¹ Cet

Another non-detection is κ¹ Cet expected at ${\mathrm{03}}^{\mathrm{h}}{\mathrm{19}}^{\mathrm{m}}\mathrm{21}\stackrel{\mathrm{s}}{.}\mathrm{93}\mathit{,}\mathrm{+}{\mathrm{03}}^{\mathrm{◦}}{\mathrm{22}}^{\mathrm{\prime }}\mathrm{13}.̋\mathrm{99}$. We therefore report upper limits for the radio emission. Because of the high sensitivity of the VLA, the surrounding noise can be measured at a very low level although a strong source showing up in the C-band image disturbs the field and contributes to the rms even after careful cleaning. The 3σ rms noise level is used for an upper limit to the radio emission, which is 9 μJy both in C-band and Ku-band.

3.3. Chromospheric emission

We expect that the emission from the stellar chromosphere is small compared to the wind emission. Nevertheless, we estimate the emission from the stellar disk of the star to occur in the chromosphere (see e.g. Drake et al. 1993, for Procyon). For an optically thick chromosphere at 10 GHz we assume a temperature of 20 000 K (White 2004). At 100 GHz we expect a lower temperature of typically 10 000 K, although it can be even lower. Furthermore, we assume that the entire surface of the star is covered by chromospheric emission. Using the standard formula for the radio flux from a blackbody with brightness temperature T, the predicted flux density at 100 GHz is ${\mathit{S}}_{\mathit{\nu }}\mathrm{=}\frac{\mathrm{4.94}\mathrm{\times }{\mathrm{10}}^{-26}}{{\mathit{d}}^{\mathrm{2}}}\mathrm{erg}\hspace{0.17em}\mathrm{c}{\mathrm{m}}^{-2}\hspace{0.17em}{\mathrm{s}}^{-1}\hspace{0.17em}\mathrm{H}{\mathrm{z}}^{-1}\mathit{.}$(2)For χ¹ Ori at the distance of d = 8.7 pc the flux density is 65 μJy for ALMA (100 GHz). Hence, part of the emission observed with ALMA can be of chromospheric origin, but it is probably not the only emission source and cannot explain the detected 100 μJy alone. At 10 GHz the expected maximum flux density is 1.3 μJy only, and therefore not significant for the VLA detections. For the non-detections, a chromosphere could probably be detected with deeper observations, as in Villadsen et al. (2014).

4. Mass loss rates

4.1. Spherically symmetric winds

Radio flux density measurements can provide estimates for mass loss rates. The radio free-free flux spectrum for an optically thick, constant velocity, fully ionized isothermal spherical wind is predicted to be of the form (Panagia & Felli 1975; Wright & Barlow 1975; Olnon 1975): ${\mathit{S}}_{\mathit{\nu }}\mathrm{=}\mathrm{0.9}\mathrm{\times }{\mathrm{10}}^{\mathrm{11}}{\left(\frac{\mathit{Ṁ}}{\mathit{v}}\right)}^{\mathrm{4}\mathit{/}\mathrm{3}}{\mathit{T}}^{\mathrm{0.1}}{\mathit{\nu }}^{\mathrm{0.6}}{\mathit{d}}^{-2}\mathrm{mJy}\mathit{,}$(3)where Ṁ is the mass loss rate in M_⊙ yr^-1, T the temperature of the plasma in K, ν the frequency in GHz, v the wind velocity in km s^-1, and d the stellar distance in pc. At any frequency one essentially sees emission from gas down to a level where the gas becomes optically thick. Wright & Barlow (1975) argue that deviations from α_op = 0.6 (which is the exponent of ν in Eq. (3)) may be caused either by variability due to non-uniform mass loss rates or by an increasing fraction of neutral gas with distance responsible for the radio emission. Using this formula and assuming a temperature of T = 10⁶ K and an average wind velocity ofv = 400 km s^-1, the mass loss rate for an (optically thick) wind of π¹ UMa would be Ṁ ≤ 1.1 × 10^-10M_⊙ yr^-1 for C-band and Ṁ ≤ 2.9 × 10^-11M_⊙ yr^-1 for Ku-band. The star κ¹ Cet would show a mass loss rate of Ṁ ≤ 2.8 × 10^-11M_⊙ yr^-1 for C-band and Ṁ ≤ 1.9 × 10^-11M_⊙ yr^-1 for Ku-band with the same assumed temperature and velocity profiles. Apart from spherically symmetric (isotropic) winds we will also discuss the possibility of anisotropic, collimated “jet” flows below.

4.1.1. Radiative transfer equation for non-isothermal winds

A point we have to consider is that the temperature in the solar wind (and presumably in winds from other stars) is not constant but decreases with distance r. Close to the surface the wind is dense and hot but it cools as it expands. This radial temperature can be roughly described by a T ∝ r^-0.5 power law (Richardson et al. 1995). Because of this, we wanted to study the case for variable temperature and therefore re-formulated the general radiation transfer equation. As described in Panagia & Felli (1975), the intensity I_ν(ξ) from any line of sight in local thermodynamic equilibrium (LTE) is given by ${\mathit{I}}_{\mathit{\nu }}\mathrm{\left(}\mathit{\xi }\mathrm{\right)}\mathrm{=}\mathit{B}\mathrm{\left(}\mathit{\nu }\mathrm{\right)}\mathrm{(}\mathrm{1}\mathrm{-}{\mathrm{e}}^{\mathrm{-}\mathit{\tau }\mathrm{\left(}\mathit{\xi }\mathrm{\right)}}\mathrm{)}\mathit{,}$(4)where ξ is the distance from the surface of the star out to a boundary of about 200 stellar radii (to ensure that the entire emission region is contained in the calculation volume) measured in the plan perpendicular to the line of sight. A grid for the temperature and density at each grid point was constructed. Emission and absorption were determined for each grid cell at a given distance from the source to create a ring structure with radius ξ around the source. Moving out to several stellar radii, the contributions from the ring elements are summed up, where the region behind the star is excluded. The optical depth along any line of sight is calculated using: ${\mathit{\tau }}_{\mathit{\nu }}\mathrm{\left(}\mathit{s}\mathrm{\right)}\mathrm{=}{\mathrm{\int }}_{\mathit{s}}^{\mathrm{\infty }}{\mathit{n}}^{\mathrm{2}}{\mathit{\kappa }}_{\mathit{\nu }}\mathrm{\left(}\mathit{T}\mathrm{\right)}\mathrm{d}\mathit{s,}$(5)where κ(ν) is defined as in Mezger & Henderson (1967): $\mathit{\kappa }\mathrm{\left(}\mathit{\nu }\mathrm{\right)}\mathrm{=}\mathrm{8.436}\mathrm{\times }{\mathrm{10}}^{-28}{\left[\frac{\mathit{\nu }}{\mathrm{10}\mathrm{GHz}}\right]}^{-2.1}{\left[\frac{{\mathit{T}}_{\mathrm{e}}}{{\mathrm{10}}^{\mathrm{4}}\mathrm{K}}\right]}^{-1.35}\mathrm{·}$(6)Taking the full geometry into account, we finally obtain a variable temperature transfer equation that can be easily solved numerically. Results are displayed in Fig. 4, shown as the red line. The black spectrum displays the solution for a constant temperature. We see that a variable temperature causes minor changes in the steepness of the spectrum which may lead to a slightly higher flux density and may influence the derived mass loss rate. This is probably because of the n² dependence of Eq. (5) and the strong dependence of the density on distance and hence the temperature, and thus most emission is from very close to the star. Using the equation of mass continuity Ṁ = 4πr²ρv, the mass loss rate will be approximately 1.1–1.6 higher if the temperature is assumed not to be constant. The change in mass loss implied by variations in the temperature is relatively small compared to those due to a change in velocity. The mass loss rate would be enhanced by a factor of about two if the velocity (see Eq. (3)) changed fromv = 400 km s^-1 tov = 800 km s^-1.

Fig. 4 Example solution (for arbitrary mass loss rate) of the radiative transfer equation assuming a non-constant temperature, shown in red. The black line represents the result for the constant temperature solution. The initial temperature is set to T = 106 K in both cases. The density is n = 2 × 1010 cm-3. The difference in both spectra is not strongly pronounced. The “peak” at lower frequencies in the red curve results from numerical issues. The arrows mark the upper limits of the observational radio fluxes of π1 UMa in both frequency bands.

4.2. Conical winds

The stars in our sample are very young and active, hence we investigate anisotropic, collimated winds where the magnetic activity is concentrated at the poles (Güdel 2007, and references therein). Reynolds (1986) showed that a well-collimated ionized flow can display a behaviour quite different from that of quasi-spherical flows and they calculated the thermal continuum emission from collimated, ionized winds (“jets”) in the presence of gradients in jet width, velocity, ionization, and temperature. Because the structure in continuum source spectra contains much information about the flow physics, it is important to get a good frequency coverage of the target sample. The total radio flux of a collimated stellar wind given by Reynolds (1986) is: ${\mathit{S}}_{\mathit{\nu }}\mathrm{=}{\mathrm{\int }}_{{\mathit{y}}_{\mathrm{0}}}^{{\mathit{y}}_{\max }}\left[\frac{\mathrm{2}\mathit{w}\mathrm{\left(}\mathit{r}\mathrm{\right)}}{{\mathit{d}}^{\mathrm{2}}}\right]\left(\frac{{\mathit{a}}_{\mathrm{j}}}{{\mathit{a}}_{\mathit{\kappa }}}\mathit{T}{\mathit{\nu }}^{\mathrm{2}}\right)\mathrm{(}\mathrm{1}\mathrm{-}{\mathrm{e}}^{\mathrm{-}\mathit{\tau }}\mathrm{)}\mathrm{d}\mathit{y}\mathit{.}$(7)Here, y is defined as y = rsini with r being the length of the jet and i its inclination (see Fig. 1 in Reynolds 1986). The half-width of the jet is described withw(r), d is the distance to the source, T the temperature, ν the frequency, and τ the optical depth based on the wind density, velocity, and temperature. The constants a_j = 6.50 × 10^-38 and a_κ = 0.212 link the free-free emission and absorption coefficients: j_ν/κ_ν = a_j/a_κTν². The jet half-width, optical depth, temperature, velocity, and density are assumed to vary with r/r₀ like power laws with indices ϵ,q_τ,q_T,q_v, and q_n, respectively. The velocity and density indirectly contribute via their indices to the optical depth in Eq. (7). Different values are assigned to each parameter, depending on the model type (Reynolds 1986), and these quantities are summarized in Table 5. For example for a constant-velocity, fully ionized, adiabatic jet the exponents are chosen to be ϵ = 1, q_n = −2, q_T = −4 / 3, q_v = 0, and q_τ = −1.2 (model B, see Table 5). These variations change the spectral index α_op to 0.83 for a non-isothermal jet instead of α_op = 0.6 for isothermal flows. Calculating Eq. (7) numerically for the properties of π¹ UMa with T = 10⁶ K, n = 2 × 10¹⁰ cm^-3, an opening angle of 40° (centred at the pole) and using the standard spherical model quantities (model A), the total flux spectrum is determined and is displayed in Fig. 5, where it reveals a positive slope of around α_op = 0.6 for the optically thick wind and a change to α = −0.1 at high frequencies for the optically thin regime.

Fig. 5 Upper limits of the total radio flux density for a jet with the properties of π1 UMa theoretically derived from Eq. (7) with an opening angle of around 40°, showing a positive slope of αop = 0.6 for the optically thick wind and a decrease at the turnover frequency for the optically thin part using the quantities for model A (see Table 5). The observational upper limit radio flux density of π1 UMa for Ku-band is shown as a black arrow. The C-band upper limit flux density would lie beyond the plotting range at 23 μJy.

To derive upper limits for mass loss rates, different values for velocity and temperature for a standard spherical jet flow, that is α_op = 0.6, are applied. It is clear that faster but cooler winds lead to stronger mass loss rates. Upper limits for mass loss rates for all three model types, the standard spherical, adiabatic spherical, and adiabatic collimated flow, are calculated using a constant-velocity wind ofv = 400 km s^-1 with the following Reynolds (1986) formula:

$\begin{array}{ccc}\mathit{Ṁ}-6\mathrm{=}\mathrm{0.938}{\mathit{v}}_{\mathrm{8}}{\mathit{x}}_{\mathrm{0}}^{-1}\left(\frac{\mathit{\mu }}{{\mathit{m}}_{\mathrm{p}}}\right)\mathrm{(}{\mathit{S}}_{\mathrm{mJy}}{\mathit{\nu }}_{\mathrm{10}}^{\mathrm{-}{\mathit{\alpha }}_{\mathrm{op}}}{\mathrm{)}}^{\mathrm{3}\mathit{/}\mathrm{4}}{\mathit{d}}_{\mathrm{kpc}}^{\mathrm{3}\mathit{/}\mathrm{2}}{\mathit{\nu }}_{\mathit{m}\mathrm{10}}^{\mathrm{-}\mathrm{0.45}\mathrm{+}\mathrm{3}{\mathit{\alpha }}_{\mathrm{op}}\mathit{/}\mathrm{4}}& & \\ \mathrm{·}{\mathit{\theta }}_{\mathrm{0}}^{\mathrm{3}\mathit{/}\mathrm{4}}{\mathit{T}}_{\mathrm{4}}^{-0.075}\mathrm{(}\sin \mathit{i}{\mathrm{)}}^{\mathrm{-}\mathrm{1}\mathit{/}\mathrm{4}}{\mathit{F}}^{\mathrm{-}\mathrm{3}\mathit{/}\mathrm{4}}\mathit{,}& & \end{array}$(8)

where Ṁ_-6 ≡ Ṁ/ 10⁶M_⊙ yr^-1, v₈ ≡ v/ 10⁸ cm s^-1, ν₁₀ ≡ ν/ 10¹⁰ Hz, T₄ ≡ T/ 10⁴ K and $\mathit{F}\mathrm{\equiv }\frac{{\mathrm{2.1}}^{\mathrm{2}}}{{\mathit{q}}_{\mathit{\tau }}\mathrm{(}{\mathit{\alpha }}_{\mathrm{op}}\mathrm{-}\mathrm{2}\mathrm{)}\mathrm{(}{\mathit{\alpha }}_{\mathrm{op}}\mathrm{+}\mathrm{0.1}\mathrm{)}}\mathrm{·}$(9)The frequency ν_m10 is defined as the turnover frequency where the source becomes completely transparent (its definition can be found in Eq. (13) in Reynolds 1986) and S_mJy is the observed radio flux of our objects in mJy. Table 6 summarizes the maximum mass loss rates for π¹ UMa and κ¹ Cet at 6 GHz and 14 GHz for all three models by changing the model parameters, a temperature of T = 10⁶ K, a velocity ofv = 400 km s^-1, and an opening angle of 40°.

If we assume that the mass loss rates are a function of the opening angle, they increase with increasing opening angle. For example, the mass loss rate with an opening angle of 20° is Ṁ ≤ 3.0 × 10^-12M_⊙ yr^-1 for π¹ UMa for Ku-band. Enlarging the angle to 60° the mass loss rate increases to Ṁ ≤ 6.7 × 10^-12M_⊙ yr^-1. A higher velocity ofv = 800 km s^-1 would raise the mass loss rates by a factor of two. Although we are not able to detect any radio emission signal for π¹ UMa and κ¹ Cet, we can thus give meaningful upper limits to the mass loss rates of these young stars within a range of reasonable wind opening angles and wind temperatures.

As already mentioned, additional coronal, partly flaring radio emission for EK Dra and χ¹ Ori was detected and identified as non-thermal emission, but we can nevertheless provide meaningful upper limits by adopting the detected non-thermal flux densities as upper limits to the thermal emission. We calculate the maximum mass loss of both stars for a spherically symmetric and a conical wind, as done for the non-detections. These mass loss rates for EK Dra and χ¹ Ori in both frequency bands are summarized in Table 7.

4.3. Absorption of the wind due to flares

The presence of flares and polarized emission on EK Dra and χ¹ Ori imply that any radio contribution from winds must be significantly lower than the detected radiation. The non-thermal and flare emission originate close to the surface of the star. The fact that it is detectable implies that the stellar wind is optically thin to this radiation. An assessment for the maximum mass loss possible for an optically thin wind was suggested in Lim & White (1996). A stronger wind would completely absorb the observed radiation from coronal radio flares. The radius at which a spherically symmetric wind becomes optically thick at a given frequency ν can be derived from the expression (Lim & White 1996): $\begin{array}{ccc}\frac{\mathit{R}\mathrm{\left(}\mathit{\nu }\mathrm{\right)}}{{\mathit{R}}_{\mathrm{\odot }}}\mathrm{\approx }\mathrm{6}{\left(\frac{\mathit{\nu }}{\mathrm{10}\mathrm{GHz}}\right)}^{\mathrm{-}\mathrm{2}\mathit{/}\mathrm{3}}{\left(\frac{\mathit{T}}{{\mathrm{10}}^{\mathrm{4}}\mathrm{K}}\right)}^{\mathrm{-}\mathrm{1}\mathit{/}\mathrm{2}}& & \\ \mathrm{\times }{\left(\frac{\mathit{Ṁ}}{{\mathrm{10}}^{-10}{\mathit{M}}_{\mathrm{\odot }}\hspace{0.17em}{\mathrm{yr}}^{-1}}\right)}^{\mathrm{2}\mathit{/}\mathrm{3}}{\left(\frac{{\mathit{v}}_{\mathrm{\infty }}}{\mathrm{300}\mathrm{km}\hspace{0.17em}{\mathrm{s}}^{-1}}\right)}^{\mathrm{-}\mathrm{2}\mathit{/}\mathrm{3}}\mathrm{·}& & \end{array}$(10)Because the non-thermal emission from the star must originate above the optically thick surface at the observing frequency to be detectable, we set R(ν) to R_∗. Assuming the terminal velocity to bev_∞ = 400 km s^-1 and the temperature T = 10⁶ K, and solving Eq. (10) for Ṁ, we find a maximum wind mass loss rate of Ṁ ≤ 1.3 × 10^-10M_⊙ yr^-1 for C-band and Ṁ ≤ 6.9 × 10^-10M_⊙ yr^-1 for Ku-band for EK Dra. With the same velocity and temperature profiles for χ¹ Ori, a wind with Ṁ ≤ 1.3 × 10^-10M_⊙ yr^-1 for C-band at 6 GHz and Ṁ ≤ 7.2 × 10^-10M_⊙ yr^-1 for Ku-band at 14 GHz is derived. Comparing these mass loss rates to those derived for a spherically symmetric and a conical wind, respectively, as given in Table 7, we see that the estimates are similar. We keep the conical wind mass loss results as upper limits for EK Dra and χ¹ Ori.

4.4. Rotational evolution

As magnetized stellar winds remove angular momentum from their host stars and therefore force stars to spin down (Weber & Davis 1967; Skumanich 1972; Kraft 1967) and cause a decrease in rotation rate and magnetic activity as they age (Güdel et al. 1997; Vidotto et al. 2014), it is essential to consider rotational evolution when determining mass loss rates of young, active stars. Several solar wind models (e.g. van der Holst et al. 2007; Zieger & Hansen 2008; Jacobs & Poedts 2011) and rotational evolution models (e.g. Cranmer & Saar 2011; Gallet & Bouvier 2015) have been developed. Johnstone et al. (2015a) developed a solar wind model to estimate the properties of stellar winds for low-mass main-sequence stars between masses of 0.4 M_⊙ and 1.1 M_⊙ at a range of distances from the star based on stellar spin-down and angular momentum loss in a magnetized wind. They used 1D thermal pressure-driven hydrodynamic wind models using the Verstile Advection Code (Tóth 1996) and in-situ measurements of the solar wind. The stellar mass loss rate can then be calculated with$\mathit{Ṁ}\mathrm{=}\dot{{\mathit{M}}_{\mathrm{\odot }}}{\mathit{R}}^{\mathrm{2}}{\mathrm{\Omega }}^{\mathrm{1.33}}{\mathit{M}}^{-3.36}\mathit{,}$(11)where all quantities are in solar units with the Carrington rotation rate of Ω_⊙ = 2.67 × 10^-6 rad s^-1. Graphically, this relation is shown in Fig. 10 in Johnstone et al. (2015b). Applying this formula to our four objects we are able to calculate their mass loss rates considering their rotational evolution, shown as red filled circles in Fig. 6. These values follow a Ṁ ∝ t^-0.75 relation (Johnstone et al. 2015b) and are about two orders of magnitude lower than our upper limits displayed as arrows, but we emphasize that these results are indirect inferences from models.

4.5. Early mass loss of the Sun

How did the solar wind evolve over time? For a simple evaluation of the total early solar mass loss, power laws are placed through the sample of young, solar-type stars observed in this study. The upper limits of the mass loss rate for the conical wind (α = 0.6) of the two non-detections of π¹ UMa and κ¹ Cet in Ku-band and the solar wind mass loss rate are used to define a piecewise power law through the sample. These relationships are shown in Fig. 6. Additionally, the mass loss rates from rotational evolution are marked as red circles. Although EK Dra and χ¹ Ori are marked in the plot, they are not used for the evaluation of the power laws, since these Ṁ are estimated based on the detected non-thermal radiation. The power laws are extrapolated from 0.3 Gyr down to the age of 0.1 Gyr. First, we apply our results to spherically symmetric winds with the corresponding power laws, which give an upper limit to the solar mass loss of 2.02% after the integration from 100 Myr to 4.5 Gyr, resulting in an initial solar mass of 1.02 M_⊙. Conical winds (using α_op = 0.6) follow similar power laws as shown in Fig. 6: $\begin{array}{ccc}\mathit{Ṁ}& \mathrm{\propto }& {\mathit{t}}^{-0.46}\mathrm{from}\mathrm{0.1}\mathrm{to}\mathrm{0.65}\mathrm{Gyr}\\ \mathit{Ṁ}& \mathrm{\propto }& {\mathit{t}}^{-2.66}\mathrm{from}\mathrm{0.65}\mathrm{to}\mathrm{4.5}\mathrm{Gyr}\mathit{.}\end{array}$Cranmer & Saar (2011) also estimate a mass loss rate versus time resulting in a power law index of –1.1, lying below our mass loss upper limits but above the model calculations for rotational evolution, shown as the blue solid line in Fig. 6. The relation of Wood (2004) as given in Eq. (1) is shown as the red line in the figure. The dashed part of the line marks the age region where the relation fails for most of young and active stars. Furthermore, Alvarado-Gómez et al. (2016) find a weaker power law relation of mass loss rate versus age based on magnetohydrodynamics (MHD) simulations, resulting in Ṁ ∝ t^-1.37.

Fig. 6 Mass loss evolution described by the non-detections π1 UMa and κ1 Cet at 0.3 and 0.65 Gyr, respectively, and the present Sun at 4.5 Gyr, shown as a black solid line assuming a conical wind with αop = 0.6 (Model A, solid line). The mass loss rates of all sources including EK Dra and χ1 Ori are determined from their observational flux densities in Ku-band and opening angles of 40°. For Models B and C, this evolution would lie below the solid line implying a lower upper limit of the solar mass. The black dashed line shows the evolution assuming a spherically symmetric wind. For EK Dra and χ1 Ori only the mass loss rates of spherically symmetric winds are shown, whereas for π1 UMa and κ1 Cet, mass loss rates of symmetric and conical winds, respectively, are shown. The red circles are the mass loss rate estimates using rotational evolution model calculations described in Johnstone et al. (2015b). The red line shows the result of Wood (2004), where the dashed line indicates the age region where the power law fails. Cranmer & Saar (2011) also estimate a mass loss rate versus time; this fit is shown as the blue solid line.

After the integration in time from 100 Myr to 4.5 Gyr, the total mass is in our case for the above given power laws at most 0.4% (α_op = 0.6) higher than at present, depending on the model for the spectral index α_op, resulting in a solar mass of 1.004 M_⊙ only. Considering the theoretical model calculation for rotational evolution, the solar mass would be even lower at 1.0002 M_⊙. The boundaries necessary for solving the FYSP are at 3% to 7% total mass loss, required to keep liquid water on early Mars and to control and avoid the runaway greenhouse effect on Earth at early stages up to a few 100 Myr (Whitmire et al. 1995). Our limits for the spherically symmetric and conical wind models are definitely below the 3% boundary and therefore imply that the faint young Sun problem cannot be solved by assuming increased wind mass loss rates and therefore a higher mass for the young Sun.

5. Summary and discussion

In this study, we analyzed four young, solar-type stars on the main-sequence of different ages, which are part of the Sun in Time programme to study the decline of magnetic activity and wind mass loss in solar analogues. For the analysis, observations of the VLA at 2 cm and 6 cm wavelength and ALMA at 100 GHz are used, aiming to detect thermal radio emission, that is free-free radio bremsstrahlung, which is indicative of the existence of a stellar wind. The well-studied analogues of the Sun cover the young evolutionary stages on the main-sequence. Our sample of four stars results in two detections: EK Dra and χ¹ Ori; and two non-detections: π¹ UMa and κ¹ Cet. For both detections we can conclude that the radio emission is not thermal bremsstrahlung alone, but consists of additional coronal radio emission in the form of non-thermal, partly flaring emission. Indicators for that assumption are a negative slope of the radio spectrum, the presence of flares seen in the light curves, and the presence of circular polarization. Furthermore, we have argued that instead of χ¹ Ori, we have actually detected its M-dwarf companion. For the non-detections π¹ UMa and κ¹ Cet, we can estimate their maximum wind radio emission flux densities by placing the 3σ rms value as upper limits. We have to clearly state that we cannot rule out other contributing but also undetected emission processes in these sources.

The estimated radio emissions are used to derive upper limits to mass loss rates for the observed targets. Mass loss rates are important quantities for the study of the evolution of young stars including the Sun. They could possibly result in an explanation and solution for the problem of the famous FYSP. Furthermore, the evolution of mass loss rates of the young Sun leads to essential information for the formation and evolution of the atmospheres of the early Earth and other planetary atmospheres. We estimate mass loss rates for all targets for spherically symmetric and anisotropic collimated winds. We note that any additional neutral wind component would increase the mass loss rate, but such winds would not be detected by our methods. However, the solar wind is essentially fully ionized, so we assume the same for solar analogues.

We applied three different model types (standard spherical, adiabatic spherical, and adiabatic collimated) for the mass loss rate calculation by changing the different parameter quantities. If we vary the velocity and the temperature, we see that the mass loss rate increases for a cooler and faster wind.

The star EK Dra’s mass loss is estimated to be Ṁ ≤ 1.3 × 10^-10M_⊙ yr^-1 for C-band and Ṁ ≤ 6.9 × 10^-10M_⊙ yr^-1 for Ku-band. The star χ¹ Ori shows a mass loss rate of Ṁ ≤ 1.3 × 10^-10M_⊙ yr^-1 for C-band at 6 GHz and Ṁ ≤ 7.2 × 10^-10M_⊙ yr^-1 for Ku-band at 14 GHz. Here, we assume that non-thermal emission from coronal radio flares contributes to the total mass loss following Lim & White (1996), implying that these features originate close to the stellar surface propagating through an optically thin wind. Mass loss rates from the spherically symmetric wind and conical wind calculations are similar to these results.

For π¹ UMa the mass loss rate is derived to be Ṁ ≤ 1.9 × 10^-11M_⊙ yr^-1 for C-band and Ṁ ≤ 5 × 10^-12M_⊙ yr^-1 for Ku-band for a jet-like wind with an opening angle of 40°. For κ¹ Cet the determined mass loss rate for a collimated wind is Ṁ ≤ 5.1 × 10^-11M_⊙ yr^-1 for C-band and Ṁ ≤ 3.5 × 10^-12M_⊙ yr^-1 for Ku-band.

The resulting maximum mass loss rate of π¹ UMa in Ku-band (with α_op = 0.6) is about 250 times stronger than the present day solar mass loss rate of Ṁ = 2 × 10^-14M_⊙ yr^-1. Wood et al. (2014) studied the stellar wind and mass loss of π¹ UMa using Ly-α observations. With hydrodynamic models for the astrosphere to infer the stellar wind strength, the study of Wood et al. (2014) results in a wind for π¹ UMa only half as strong as the solar wind. From their research, Wood et al. (2014) concluded that the Sun and solar-like stars do not experience particularly strong coronal winds in their past. Drake et al. (2013) studied coronal mass ejections in connection to stellar winds, where the authors found that coronal mass ejection (CME) induced mass loss rates can amount to several percent of the steady wind rate. Their estimation for a CME mass loss rate for π¹ UMa implies Ṁ ~ 3 × 10^-12M_⊙ yr^-1, comparable with our upper limits. For κ¹ Cet the mass loss rate in the Drake et al. (2013) study is similar. We see that the measurements of Wood et al. (2014) of Ṁ = 0.5 Ṁ_⊙ are much lower than the Ṁ = 150 Ṁ_⊙ predictions of Drake et al. (2013) and our observational upper limit. The rotation-wind model by Johnstone et al. (2015b) in fact also requires a wind mass loss rate significantly above the value suggested by Wood et al. (2014) to explain the observed spin-down rate for solar-like stars in this age range.

Finally, the maximum total solar mass for the young Sun was derived for three cases: spherically symmetric winds, conical jet flows, and rotational evolution models. The results are quite different: 1) the mass loss rate of spherically symmetric winds indicates a total maximum mass of 1.02 M_⊙; 2) conical winds lead to a total mass of 1.004 M_⊙; and 3) the rotational evolution model suggests an initial solar mass of only 1.0002 M_⊙ at an age of about 100 Myr.

6. Conclusion

If the FYSP is to be solved with a larger initial solar mass, Earth and Mars climate constraints require the solar mass to be in the range of 1.03–1.07 M_⊙ near the zero-age main-sequence, requiring an enhanced early wind mass loss rate of order 10^-12−10^-10M_⊙ yr^-1. Our results for mass loss rates derived with radio observations of solar analogues indicate an early solar mass of at most 1.02 M_⊙ assuming spherically symmetric winds. This is not sufficient to solve the faint young Sun paradox. It appears that other explanations such as higher concentrations of greenhouse gases and aerosols (e.g. Sagan & Mullen 1972; Kasting 1993), a lower global albedo, either through less cloud coverage (e.g. Shaviv 2003), and/or a smaller continental land mass (Rosing et al. 2010) are required.

Acknowledgments

We thank the referee, Jeffrey Linsky, for very helpful comments that improved the paper. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. B.F and M.G. acknowledge the support of the FWF “Nationales Forschungsnetzwerk” project S116601-N16 “Pathways to Habitability: From Disks to Active Stars, Planets and Life” and the related FWF NFN subproject S116604-N16 “Radiation and Wind Evolution from the T Tauri Phase to ZAMS and Beyond”. Financial support of this project by the University of Vienna is also acknowledged. This publication is supported by the Austrian Science Fund (FWF). The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This paper makes use of the following ALMA data: ADS / JAO.ALMA#2011.0.01234.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ.

References

All Tables

All Figures

	Fig. 1 Contours of EK Dra in a) C-band with contour levels [0.2, 0.4, 0.6, 0.8] · 576 μJy; and b) Ku-band with [0.2, 0.4, 0.6, 0.8] · 74.5 μJy and the beam size in the left corner. The point in the centre gives the coordinates from the Gaussian fit, the cross marks the proper motion corrected Hipparcos position.
In the text

	Fig. 2 χ1 Ori in a) C-band with contour levels [0.2, 0.4, 0.6, 0.8] · 223 μJy; b) Ku-band with [0.2, 0.4, 0.6, 0.8] · 106.6 μJy; and c) for ALMA at 100 GHz with [0.2, 0.4, 0.6, 0.8] · 108 μJy. The black dot represents the coordinates calculated from the Gaussian fit. The red cross marks the position of the companion χ1 Ori B and the black cross the corrected coordinates of the main target χ1 Ori.
In the text

	Fig. 3 Time series of χ1 Ori in C-band showing an explosive increase in intensity identified as a flare event at the central frequency of 6.15 GHz. The time interval between the data points is around ten minutes.
In the text

Fig. 4 Example solution (for arbitrary mass loss rate) of the radiative transfer equation assuming a non-constant temperature, shown in red. The black line represents the result for the constant temperature solution. The initial temperature is set to T = 106 K in both cases. The density is n = 2 × 1010 cm-3. The difference in both spectra is not strongly pronounced. The “peak” at lower frequencies in the red curve results from numerical issues. The arrows mark the upper limits of the observational radio fluxes of π1 UMa in both frequency bands.

In the text

Fig. 5 Upper limits of the total radio flux density for a jet with the properties of π1 UMa theoretically derived from Eq. (7) with an opening angle of around 40°, showing a positive slope of αop = 0.6 for the optically thick wind and a decrease at the turnover frequency for the optically thin part using the quantities for model A (see Table 5). The observational upper limit radio flux density of π1 UMa for Ku-band is shown as a black arrow. The C-band upper limit flux density would lie beyond the plotting range at 23 μJy.

In the text

Fig. 6 Mass loss evolution described by the non-detections π1 UMa and κ1 Cet at 0.3 and 0.65 Gyr, respectively, and the present Sun at 4.5 Gyr, shown as a black solid line assuming a conical wind with αop = 0.6 (Model A, solid line). The mass loss rates of all sources including EK Dra and χ1 Ori are determined from their observational flux densities in Ku-band and opening angles of 40°. For Models B and C, this evolution would lie below the solid line implying a lower upper limit of the solar mass. The black dashed line shows the evolution assuming a spherically symmetric wind. For EK Dra and χ1 Ori only the mass loss rates of spherically symmetric winds are shown, whereas for π1 UMa and κ1 Cet, mass loss rates of symmetric and conical winds, respectively, are shown. The red circles are the mass loss rate estimates using rotational evolution model calculations described in Johnstone et al. (2015b). The red line shows the result of Wood (2004), where the dashed line indicates the age region where the power law fails. Cranmer & Saar (2011) also estimate a mass loss rate versus time; this fit is shown as the blue solid line.

In the text