Saturation of Stellar Winds from Young Suns

We investigate mass losses via stellar winds from sun-like main sequence stars with a wide range of activity levels. We perform forward-type magnetohydrodynamical numerical experiments for Alfven wave-driven stellar winds with a wide range of the input Poynting flux from the photosphere. Increasing the magnetic field strength and the turbulent velocity at the stellar photosphere from the current solar level, the mass loss rate rapidly increases at first owing to the suppression of the reflection of the Alfven waves. The surface materials are lifted up by the magnetic pressure associated with the Alfven waves, and the cool dense chromosphere is intermittently extended to 10 -- 20 % of the stellar radius. The dense atmospheres enhance the radiative losses and eventually most of the input Poynting energy from the stellar surface escapes by the radiation. As a result, there is no more sufficient energy remained for the kinetic energy of the wind; the stellar wind saturates in very active stars, as observed in Wood et al. The saturation level is positively correlated with B_{r,0}f_0, where B_{r,0} and f_0 are the magnetic field strength and the filling factor of open flux tubes at the photosphere. If B_{r,0}f_0 is relatively large >~ 5 G, the mass loss rate could be as high as 1000 times. If such a strong mass loss lasts for ~ 1 billion years, the stellar mass itself is affected, which could be a solution to the faint young sun paradox. We derive a Reimers-type scaling relation that estimates the mass loss rate from the energetics consideration of our simulations. Finally, we derive the evolution of the mass loss rates, dot{M} t^{-1.23}, of our simulations, combining with an observed time evolution of X-ray flux from sun-like stars, which is shallower than dot{M} t^{-2.33+/-0.55} in Wood et al.(2005).