No Arabic abstract
We present a model for the rotational evolution of a young, solar-mass star interacting magnetically with an accretion disk. As in a previous paper (Paper I), the model includes changes in the stars mass and radius as it descends the Hayashi track, a decreasing accretion rate, and a prescription for the angular momentum transfer between the star and disk. Paper I concluded that, for the relatively strong magnetic coupling expected in real systems, additional processes are necessary to explain the existence of slowly rotating pre-main-sequence stars. In the present paper, we extend the stellar spin model to include the effect of a spin-down torque that arises from an accretion-powered stellar wind. For a range of magnetic field strengths, accretion rates, initial spin rates, and mass outflow rates, the modeled stars exhibit rotation periods within the range of 1--10 days in the age range of 1--3 Myr. This range coincides with the bulk of the observed rotation periods, with the slow rotators corresponding to stars with the lowest accretion rates, strongest magnetic fields, and/or highest stellar wind mass outflow rates. We also make a direct, quantitative comparison between the accretion-powered stellar wind scenario and the two types of disk-locking models (namely the X-wind and Ghosh & Lamb type models) and identify some remaining theoretical issues for understanding young star spins.
Supermassive stars (SMS; ~ 10^5 M_sun) formed from metal-free gas in the early Universe attract attention as progenitors of supermassive black holes observed at high redshifts. To form SMSs by accretion, central protostars must accrete at as high rates as ~ 0.1-1 M_sun/yr. Such protostars have very extended structures with bloated envelopes, like super-giant stars, and are called super-giant protostars (SGPSs). Under the assumption of hydrostatic equilibrium, SGPSs have density inverted layers, where the luminosity becomes locally super-Eddington, near the surface. If the envelope matter is allowed to flow out, however, a stellar wind could be launched and hinder the accretion growth of SGPSs before reaching the supermassive regime. We examine whether radiation-driven winds are launched from SGPSs by constructing steady and spherically symmetric wind solutions. We find that the wind velocity does not reach the escape velocity in any case considered. This is because once the temperature falls below ~ 10^4 K, the opacity plummet drastically owing to the recombination of hydrogen and the acceleration ceases suddenly. This indicates that, in realistic non-steady cases, even if outflows are launched from the surface of SGPSs, they would fall back again. Such a wind does not result in net mass loss and does not prevent the growth of SGPSs. In conclusion, SGPSs will grow to SMSs and eventually collapse to massive BHs of ~ 10^5 M_sun, as long as the rapid accretion is maintained.
Stellar evolution models are a cornerstone of young star astrophysics, which necessitates that they yield accurate and reliable predictions of stellar properties. Here, I review the current performance of stellar evolution models against young astrophysical benchmarks and highlight recent progress incorporating non-standard physics, such as magnetic field and starspots, to explain observed deficiencies. While addition of these physical processes leads to improved agreement between models and observations, there are several fundamental limitations in our understanding about how these physical processes operate. These limitations inhibit our ability to form a coherent picture of the essential physics needed to accurately compute young stellar models, but provide rich avenues for further exploration.
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).
We investigated stellar winds from zero/low-metallicity low-mass stars by magnetohydrodynamical simulations for stellar winds driven by Alfven waves from stars with mass $M_{star}=(0.6-0.8)M_{odot}$ and metallicity $Z=(0-1)Z_{odot}$, where $M_{odot}$ and $Z_{odot}$ are the solar mass and metallicity, respectively. Alfvenic waves, which are excited by the surface convection, travel upward from the photosphere and heat up the corona by their dissipation. For lower $Z$, denser gas can be heated up to the coronal temperature because of the inefficient radiation cooling. The coronal density of Pop.II/III stars with $Zle 0.01Z_{odot}$ is 1-2 orders of magnitude larger than that of the solar-metallicity star with the same mass, and as a result, the mass loss rate, $dot{M}$, is $(4.5-20)$ times larger. This indicates that metal accretion on low-mass Pop.III stars is negligible. The soft X-ray flux of the Pop.II/III stars is also expected to be $approx (1-30)$ times larger than that of the solar-metallicity counterpart owing to the larger coronal density, even though the radiation cooling efficiency is smaller. A larger fraction of the input Alfvenic wave energy is transmitted to the corona in low $Z$ stars because they avoid severe reflection owing to the smaller density difference between the photosphere and the corona. Therefore, a larger fraction is converted to the thermal energy of the corona and the kinetic energy of the stellar wind. From this energetics argument, we finally derived a scaling of $dot{M}$ as $dot{M}propto L R_{star}^{11/9}M_{star}^{-10/9}T_{rm eff}^{11/2}left[max (Z/Z_{odot},0.01)right]^{-1/5}$, where $L$, $R_{star}$, and $T_{rm eff}$ are stellar luminosity, radius, and effective temperature, respectively.
This work is the first attempt to describe the multi-dimensional structure of accreting young stars based on fully compressible time implicit multi-dimensional hydrodynamics simulations. One major motivation is to analyse the validity of accretion treatment used in previous 1D stellar evolution studies. We analyse the effect of accretion on the structure of a realistic stellar model of the young Sun. Our work is inspired by the numerical work of Kley & Lin (1996, ApJ, 461, 933) devoted to the structure of the boundary layer in accretion disks. We analyse the redistribution of accreted material with a range of values of specific entropy relative to the bulk specific entropy of the material in the accreting objects convective envelope. A primary goal is to understand whether and how accreted energy deposited onto a stellar surface is redistributed in the interior. This study focusses on the high accretion rates characteristic of FU Ori systems. We find that the highest entropy cases produce a distinctive behaviour in the mass redistribution, rms velocities, and enthalpy flux in the convective envelope. This change in behaviour is characterised by the formation of a hot layer on the surface of the accreting object, which tends to suppress convection in the envelope. We analyse the long-term effect of such a hot buffer zone on the structure and evolution of the accreting object with 1D stellar evolution calculations. We study the relevance of the assumption of redistribution of accreted energy into the stellar interior used in the literature. One conclusion is that, for a given amount of accreted energy transferred to the accreting object, a treatment assuming accretion energy redistribution throughout the stellar interior could significantly overestimate the effects on the stellar structure, in particular, on the resulting expansion.