Observations of stellar rotation show that low-mass stars lose angular momentum during the main sequence. We simulate the winds of Sun-like stars with a range of rotation rates, covering the fast and slow magneto-rotator regimes, including the transition between the two. We generalize an Alfv\'en-wave driven solar wind model that builds on previous works by including the magneto-centrifugal force explicitly. In this model, the surface-averaged open magnetic flux is assumed to scale as $B_\ast f^{\rm open}_\ast \propto {\rm Ro}^{-1.2}$, where $f^{\rm open}_\ast$ and ${\rm Ro}$ are the surface open-flux filling factor and Rossby number, respectively. We find that, 1. the angular momentum loss rate (torque) of the wind is described as $\tau_w \approx 2.59 \times 10^{30} {\rm \ erg} \ \left( \Omega_\ast / \Omega_\odot \right)^{2.82}$, yielding a spin-down law $\Omega_\ast \propto t^{-0.55}$. 2. the mass-loss rate saturates at $\dot{M}_w \sim 3.4 \times 10^{-14} M_\odot {\rm \ yr^{-1}}$, due to the strong reflection and dissipation of Alfv\'en waves in the chromosphere. This indicates that the chromosphere has a strong impact in connecting the stellar surface and stellar wind. Meanwhile, the wind ram pressure scales as $P_w \propto \Omega_\ast^{0.57}$, which is able to explain the lower-envelope of the observed stellar winds by Wood et al. 3. the location of the Alfv\'en radius is shown to scale in a way that is consistent with 1D analytic theory. Additionally, the precise scaling of the Alfv\'en radius matches previous works which used thermally-driven winds. Our results suggest that the Alfv\'en-wave driven magnetic rotator wind plays a dominant role in the stellar spin-down during the main-sequence.

We address the origin of the observed bimodal rotational distribution of stars in massive young and intermediate age stellar clusters. This bimodality is seen as split main sequences at young ages and also has been recently directly observed in the $Vsini$ distribution of stars within massive young and intermediate age clusters. Previous models have invoked binary interactions as the origin of this bimodality, although these models are unable to reproduce all of the observational constraints on the problem. Here we suggest that such a bimodal rotational distribution is set up early within a cluster's life, i.e., within the first few Myr. Observations show that the period distribution of low-mass ($\la 2 M_\odot$) pre-main sequence (PMS) stars is bimodal in many young open clusters and we present a series of models to show that if such a bimodality exists for stars on the PMS that it is expected to manifest as a bimodal rotational velocity (at fixed mass/luminosity) on the main sequence for stars with masses in excess of $\sim1.5$~\msun. Such a bimodal period distribution of PMS stars may be caused by whether stars have lost (rapid rotators) or been able to retain (slow rotators) their circumstellar discs throughout their PMS lifetimes. We conclude with a series of predictions for observables based on our model.

We present a sample of 14 M dwarf stars observed with XMM-Newton and Chandra, for which we also computed rotational periods from Kepler Two-Wheel (K2) Mission light curves. We compiled X-ray and rotation data from the literature and homogenized all data sets to provide the largest uniform sample of M dwarfs (302 stars) for X-ray activity and rotation studies to date. We then fit the relation between $L_{\rm x} - P_{\rm rot}$ using three different mass bins to separate partially and fully convective stars. We found a steeper slope in the unsaturated regime for fully convective stars and a nonconstant $L_{\rm x}$ level in the saturated regime for all masses. In the $L_{\rm x}/L_{\rm bol}-R_{\rm O}$ space we discovered a remarkable double gap that might be related to a discontinuous period evolution. Then we combined the evolution of $P_{\rm rot}$ predicted by angular momentum evolution models with our new results on the empirical $L_{\rm x} - P_{\rm rot}$ relation to provide an estimate for the age decay of X-ray luminosity. We compare predictions of this relationship with the actual X-ray luminosities of M stars with known ages from 100 Myr to a few billion years. We find remarkably good agreement between the predicted $L_{\rm x}$ and the observed values for partially convective stars. However, for fully convective stars at ages of a few billion years, the constructed $L_{\rm x}-$age relation overpredicts the X-ray luminosity because the angular momentum evolution model underpredicts the rotation period of these stars. Finally, we examine the effect of different parameterizations for the Rossby number ($R_{\rm O}$) on the shape of the activity-rotation relation in $L_{\rm x}/L_{\rm bol}-R_{\rm O}$ space, and we find that the slope in the unsaturated regime and the location of the break point of the dual power-law depend sensitively on the choice of $R_{\rm O}$.

Numerous attempts to estimate the rate at which low-mass stars lose angular momentum over their lifetimes exist in the literature. One approach is to use magnetic maps derived from Zeeman-Doppler imaging (ZDI) in conjunction with so-called "braking laws". The use of ZDI maps has advantages over other methods because it allows information about the magnetic field geometry to be incorporated into the estimate. However, ZDI is known to underestimate photospheric field strengths due to flux cancellation effects. Recently, Lehmann et al. (2018) conducted synthetic ZDI reconstructions on a set of flux transport simulations to help quantify the amount by which ZDI underestimates the field strengths of relatively slowly rotating and weak activity solar-like stars. In this paper, we evaluate how underestimated angular momentum-loss rate estimates based on ZDI maps may be. We find that they are relatively accurate for stars with strong magnetic fields but may be underestimated by a factor of up to $\sim$10 for stars with weak magnetic fields. Additionally, we re-evaluate our previous work that used ZDI maps to study the relative contributions of different magnetic field modes to angular momentum-loss. We previously found that the dipole component dominates spin-down for most low-mass stars. This conclusion still holds true even in light of the work of Lehmann et al. (2018).

The rotation rates and magnetic activity of Sun-like and low-mass (< 1.4 Msun) main-sequence stars are known to decline with time, and there now exist several models for the evolution of rotation and activity. However, the role that chemical composition plays during stellar spin-down has not yet been explored. In this work, we use a structural evolution code to compute the rotational evolution of stars with three different masses (0.7, 1.0, and 1.3Msun and six different metallicities, ranging from [Fe/H]=-1.0 to [Fe/H]=+0.5. We also implement three different wind-braking formulations from the literature (two modern and one classical) and compare their predictions for rotational evolution. The effect that metallicity has on stellar structural properties, and in particular the convective turnover timescale, leads the two modern wind-braking formulations to predict a strong dependence of the torque on metallicity. Consequently, they predict that metal rich stars spin-down more effectively at late ages (> 1 Gyr) than metal poor stars, and the effect is large enough to be detectable with current observing facilities. For example, the formulations predict that a Sun-like (solar-mass and solar-aged) star with [Fe/H]=-0.3 will have a rotation period of less than 20 days. Even though old, metal poor stars are predicted to rotate more rapidly at a given age, they have larger Rossby numbers and are thus expected to have lower magnetic activity levels. Finally, the different wind-braking formulations predict quantitative differences in the metallicity-dependence of stellar rotation, which may be used to test them.

The rate at which the solar wind extracts angular momentum from the Sun has been predicted by theoretical models for many decades, and yet we lack a conclusive measurement from in-situ observations. In this letter we present a new estimate of the time-varying angular momentum flux in the equatorial solar wind, as observed by the \textit{Wind} spacecraft from 1994-2019. We separate the angular momentum flux into contributions from the protons, alpha particles, and magnetic stresses, showing that the mechanical flux in the protons is $\sim$3 times larger than the magnetic field stresses. We observe the tendency for the angular momentum flux of fast wind streams to be oppositely signed to the slow wind streams, as noted by previous authors. From the average total flux, we estimate the global angular momentum loss rate of the Sun to be $3.3\times10^{30}$erg, which lies within the range of various MHD wind models in the literature. This angular momentum loss rate is a factor of $\sim$2 weaker than required for a Skumanich-like rotation period evolution ($\Omega_*\propto $ stellar age$^{-1/2}$), which should be considered in studies of the rotation period evolution of Sun-like stars.

Main sequence low-mass stars are known to spin-down as a consequence of their magnetised stellar winds. However, estimating the precise rate of this spin-down is an open problem. The mass-loss rate, angular momentum-loss rate and the magnetic field properties of low-mass stars are fundamentally linked making this a challenging task. Of particular interest is the stellar magnetic field geometry. In this work, we consider whether non-dipolar field modes contribute significantly to the spin-down of low-mass stars. We do this using a sample of stars that have all been previously mapped with Zeeman-Doppler imaging. For a given star, as long as its mass-loss rate is below some critical mass-loss rate, only the dipolar fields contribute to its spin-down torque. However, if it has a larger mass-loss rate, higher order modes need to be considered. For each star, we calculate this critical mass-loss rate, which is a simple function of the field geometry. Additionally, we use two methods of estimating mass-loss rates for our sample of stars. In the majority of cases, we find that the estimated mass-loss rates do not exceed the critical mass-loss rate and hence, the dipolar magnetic field alone is sufficient to determine the spin-down torque. However, we find some evidence that, at large Rossby numbers, non-dipolar modes may start to contribute.

The Sun and Sun-like stars lose angular momentum to their magnetised stellar winds. This braking torque is coupled to the stellar magnetic field, such that changes in the strength and/or geometry of the field modifies the efficiency of this process. Since the space-age, we have been able to directly measure solar wind properties using in-situ spacecraft. Furthermore, indirect proxies such as sunspot number, geomagnetic indices, and cosmogenic radionuclides, constrain the variation of solar wind properties on centennial, and millennial timescales. We use near-Earth measurements of the solar wind plasma and magnetic field to calculate the torque on the Sun throughout the space-age. Then, reconstructions of the solar open magnetic flux are used to estimate the time-varying braking torque during the last nine millennia. We assume a relationship for the solar mass loss rate based on observations during the space-age which, due to the weak dependence of the torque on mass loss rate, does not strongly affect our predicted torque. The average torque during the last nine millennia is found to be $2.2\times10^{30}$erg, which is comparable to the average value from the last two decades. Our dataset includes grand minima (such as the Maunder Minimum), and maxima in solar activity, where the torque varies from $\sim1-5\times10^{30}$erg (averaged on decadal timescales), respectively. We find no evidence for any secular variation of the torque on timescales of less than $9000$ years.

The magnetic fields of low-mass stars are observed to be variable on decadal timescales, ranging in behaviour from cyclic to stochastic. The changing strength and geometry of the magnetic field should modify the efficiency of angular momentum loss by stellar winds, but this has not been well quantified. In Finley et al. (2018) we investigated the variability of the Sun, and calculated the time-varying angular momentum loss rate in the solar wind. In this work, we focus on four low-mass stars that have all had their surface magnetic fields mapped for multiple epochs. Using mass loss rates determined from astrospheric Lyman-$\alpha$ absorption, in conjunction with scaling relations from the MHD simulations of Finley & Matt (2018), we calculate the torque applied to each star by their magnetised stellar winds. The variability of the braking torque can be significant. For example, the largest torque for $\epsilon$ Eri is twice its decadal averaged value. This variation is comparable to that observed in the solar wind, when sparsely sampled. On average, the torques in our sample range from 0.5-1.5 times their average value. We compare these results to the torques of Matt et al. (2015), which use observed stellar rotation rates to infer the long-time averaged torque on stars. We find that our stellar wind torques are systematically lower than the long-time average values, by a factor of ~3-30. Stellar wind variability appears unable to resolve this discrepancy, implying that there remain some problems with observed wind parameters, stellar wind models, or the long-term evolution models, which have yet to be understood.

Low-mass stars are known to have magnetic fields that are believed to be of dynamo origin. Two complementary techniques are principally used to characterise them. Zeeman-Doppler imaging (ZDI) can determine the geometry of the large-scale magnetic field while Zeeman broadening can assess the total unsigned flux including that associated with small-scale structures such as spots. In this work, we study a sample of stars that have been previously mapped with ZDI. We show that the average unsigned magnetic flux follows an activity-rotation relation separating into saturated and unsaturated regimes. We also compare the average photospheric magnetic flux recovered by ZDI, $\langle B_V\rangle$, with that recovered by Zeeman broadening studies, $\langle B_I\rangle$. In line with previous studies, $\langle B_V\rangle$ ranges from a few % to $\sim$20% of $\langle B_I\rangle$. We show that a power law relationship between $\langle B_V\rangle$ and $\langle B_I\rangle$ exists and that ZDI recovers a larger fraction of the magnetic flux in more active stars. Using this relation, we improve on previous attempts to estimate filling factors, i.e. the fraction of the stellar surface covered with magnetic field, for stars mapped only with ZDI. Our estimated filling factors follow the well-known activity-rotation relation which is in agreement with filling factors obtained directly from Zeeman broadening studies. We discuss the possible implications of these results for flux tube expansion above the stellar surface and stellar wind models.

The rotational evolution of cool stars is governed by magnetised stellar winds which slow the stellar rotation during their main sequence lifetimes. Magnetic variability is commonly observed in Sun-like stars, and the changing strength and topology of the global field is expected to affect the torque exerted by the stellar wind. We present three different methods for computing the angular momentum loss in the solar wind. Two are based on MHD simulations from Finley & Matt (2018), with one using the open flux measured in the solar wind, and the other using remotely-observed surface magnetograms. Both methods agree in the variation of the solar torque seen through the solar cycle and show a 30-40% decrease from cycle 23 to 24. The two methods calculate different average values, $2.9\times10^{30}$erg (open flux) and $0.35\times10^{30}$erg (surface field). This discrepancy results from the already well-known difficulty with reconciling the magnetograms with observed open flux, which is currently not understood, leading to an inability to discriminate between these two calculated torques. The third method is based on the observed spin-rates of Sun-like stars, which decrease with age, directly probing the average angular momentum loss. This method gives $6.2\times10^{30}$erg for the solar torque, larger than the other methods. This may be indicative of further variability in the solar torque on timescales much longer than the magnetic cycle. We discuss the implications for applying the formula to other Sun-like stars, where only surface field measurements are available, and where the magnetic variations are ill-constrained.

During the lifetime of sun-like or low mass stars a significant amount of angular momentum is removed through magnetised stellar winds. This process is often assumed to be governed by the dipolar component of the magnetic field. However, observed magnetic fields can host strong quadrupolar and/or octupolar components, which may influence the resulting spin-down torque on the star. In Paper I, we used the MHD code PLUTO to compute steady state solutions for stellar winds containing a mixture of dipole and quadrupole geometries. We showed the combined winds to be more complex than a simple sum of winds with these individual components. This work follows the same method as Paper I, including the octupole geometry which increases the field complexity but also, more fundamentally, looks for the first time at combining the same symmetry family of fields, with the field polarity of the dipole and octupole geometries reversing over the equator (unlike the symmetric quadrupole). We show, as in Paper I, that the lowest order component typically dominates the spin down torque. Specifically, the dipole component is the most significant in governing the spin down torque for mixed geometries and under most conditions for real stars. We present a general torque formulation that includes the effects of complex, mixed fields, which predicts the torque for all the simulations to within 20% precision, and the majority to within ~5%. This can be used as an input for rotational evolution calculations in cases where the individual magnetic components are known.

Sun-like and low-mass stars possess high temperature coronae and lose mass in the form of stellar winds, driven by thermal pressure and complex magnetohydrodynamic processes. These magnetized outflows probably do not significantly affect the star's structural evolution on the Main Sequence, but they brake the stellar rotation by removing angular momentum, a mechanism known as magnetic braking. Previous studies have shown how the braking torque depends on magnetic field strength and geometry, stellar mass and radius, mass-loss rate, and the rotation rate of the star, assuming a fixed coronal temperature. For this study we explore how different coronal temperatures can influence the stellar torque. We employ 2.5D, axisymmetric, magnetohydrodynamic simulations, computed with the PLUTO code, to obtain steady-state wind solutions from rotating stars with dipolar magnetic fields. Our parameter study includes 30 simulations with variations in coronal temperature and surface-magnetic-field strength. We consider a Parker-like (i.e. thermal-pressure-driven) wind, and therefore coronal temperature is the key parameter determining the velocity and acceleration profile of the flow. Since the mass loss rates for these types of stars are not well constrained, we determine how torque scales for a vast range of stellar mass loss rates. Hotter winds lead to a faster acceleration, and we show that (for a given magnetic field strength and mass-loss rate) a hotter outflow leads to a weaker torque on the star. We derive new predictive torque formulae for each temperature, which quantifies this effect over a range of possible wind acceleration profiles.

Cool stars with outer convective envelopes are observed to have magnetic fields with a variety of geometries, which on large scales are dominated by a combination of the lowest order fields such as the dipole, quadrupole and octupole modes. Magnetised stellar wind outflows are primarily responsible for the loss of angular momentum from these objects during the main sequence. Previous works have shown the reduced effectiveness of the stellar wind braking mechanism with increasingly complex, but singular, magnetic field geometries. In this paper, we quantify the impact of mixed dipolar and quadrupolar fields on the spin-down torque using 50 MHD simulations with mixed field, along with 10 of each pure geometries. The simulated winds include a wide range of magnetic field strength and reside in the slow-rotator regime. We find that the stellar wind braking torque from our combined geometry cases are well described by a broken power law behaviour, where the torque scaling with field strength can be predicted by the dipole component alone or the quadrupolar scaling utilising the total field strength. The simulation results can be scaled and apply to all main-sequence cool stars. For Solar parameters, the lowest order component of the field (dipole in this paper) is the most significant in determining the angular momentum loss.

We seek to characterize how the change of global rotation rate influences the overall dynamics and large scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the latter profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number $R_{of} = \tilde{\omega}/2\Omega_* \sim \tilde{v}/2\Omega_* R_*$, which we calibrate based on our 3-D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map ($R_{of} \gtrsim 1.5-2$) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

We study the relation between stellar rotation and magnetic activity for a sample of 134 bright, nearby M dwarfs observed in the Kepler Two-Wheel (K2) mission during campaigns C0 to C4. The K2 lightcurves yield photometrically derived rotation periods for 97 stars (79 of which without previous period measurement), as well as various measures for activity related to cool spots and flares. We find a clear difference between fast and slow rotators with a dividing line at a period of ~10d at which the activity level changes abruptly. All photometric diagnostics of activity (spot cycle amplitude, flare peak amplitude and residual variability after subtraction of spot and flare variations) display the same dichotomy, pointing to a quick transition between a high-activity mode for fast rotators and a low-activity mode for slow rotators. This unexplained behavior is reminiscent of a dynamo mode-change seen in numerical simulations that separates a dipolar from a multipolar regime. A substantial number of the fast rotators are visual binaries. A tentative explanation is accelerated disk evolution in binaries leading to higher initial rotation rates on the main-sequence and associated longer spin-down and activity lifetimes. We combine the K2 rotation periods with archival X-ray and UV data. X-ray, FUV and NUV detections are found for 26, 41, and 11 stars from our sample, respectively. Separating the fast from the slow rotators, we determine for the first time the X-ray saturation level separately for early- and for mid-M stars.

Magnetic interactions between a star and a close-in planet are postulated to be a source of enhanced emissions and to play a role in the secular evolution of the orbital system. Close-in planets generally orbit in the sub-alfv\'enic region of the stellar wind, which leads to efficient transfers of energy and angular momentum between the star and the planet. We model the magnetic interactions occurring in close-in star-planet systems with three-dimensional, global, compressible magneto-hydrodynamic numerical simulations of a planet orbiting in a self-consistent stellar wind. We focus on the cases of magnetized planets and explore three representative magnetic configurations. The Poynting flux originating from the magnetic interactions is an energy source for enhanced emissions in star-planet systems. Our results suggest a simple geometrical explanation for ubiquitous on/off enhanced emissions associated with close-in planets, and confirm that the Poynting fluxes can reach powers of the order of $10^{19}$ W. Close-in planets are also showed to migrate due to magnetic torques for sufficiently strong stellar wind magnetic fields. The topology of the interaction significantly modifies the shape of the magnetic obstacle that leads to magnetic torques. As a consequence, the torques can vary by at least an order of magnitude as the magnetic topology of the interaction varies.

The results from a spectro-polarimetric study of the planet-hosting Sun-like star, HD 147513 (G5V), are presented here. Robust detections of Zeeman signatures at all observed epochs indicate a surface magnetic field, with longitudinal magnetic field strengths varying between 1.0-3.2 G. Radial velocity variations from night to night modulate on a similar timescale to the longitudinal magnetic field measurements. These variations are therefore likely due to the rotational modulation of stellar active regions rather than the much longer timescale of the planetary orbit (Porb=528 d). Both the longitudinal magnetic field measurements and radial velocity variations are consistent with a rotation period of 10 +/- 2 days, which are also consistent with the measured chromospheric activity level of the star (log R'(HK)=-4.64). Together, these quantities indicate a low inclination angle, i~18 degrees. We present preliminary magnetic field maps of the star based on the above period and find a simple poloidal large-scale field. Chemical analyses of the star have revealed that it is likely to have undergone a barium-enrichment phase in its evolution because of a higher mass companion. Despite this, our study reveals that the star has a fairly typical activity level for its rotation period and spectral type. Future studies will enable us to explore the long-term evolution of the field, as well as to measure the stellar rotation period, with greater accuracy.

Observations of surface magnetic fields are now within reach for many stellar types thanks to the development of Zeeman-Doppler Imaging. These observations are extremely useful for constraining rotational evolution models of stars, as well as for characterizing the generation of magnetic field. We recently demonstrated that the impact of coronal magnetic field topology on the rotational braking of a star can be parametrized with a scalar parameter: the open magnetic flux. However, without running costly numerical simulations of the stellar wind, reconstructing the coronal structure of the large scale magnetic field is not trivial. An alternative -broadly used in solar physics- is to extrapolate the surface magnetic field assuming a potential field in the corona, to describe the opening of the field lines by the magnetized wind. This technique relies on the definition of a so-called source surface radius, which is often fixed to the canonical value of 2.5Rsun. However this value likely varies from star to star. To resolve this issue, we use our extended set of 2.5D wind simulations published in 2015, to provide a criteria for the opening of field lines as well as a simple tool to assess the source surface radius and the open magnetic flux. This allows us to derive the magnetic torque applied to the star by the wind from any spectropolarimetric observation. We conclude by discussing some estimations of spin-down time scales made using our technique, and compare them to observational requirements.

We analyse the magnetic activity characteristics of the planet hosting Sun-like star, HD 1237, using HARPS spectro-polarimetric time-series data. We find evidence of rotational modulation of the magnetic longitudinal field measurements consistent with our ZDI analysis, with a period of 7 days. We investigate the effect of customising the LSD mask to the line depths of the observed spectrum and find that it has a minimal effect on shape of the extracted Stokes V profile but does result in a small increase in the S/N ($\sim$ 7%). We find that using a Milne-Eddington solution to describe the local line profile provides a better fit to the LSD profiles in this slowly rotating star, which also impacts the recovered ZDI field distribution. We also introduce a fit-stopping criterion based on the information content (entropy) of the ZDI maps solution set. The recovered magnetic field maps show a strong (+90 G) ring-like azimuthal field distribution and a complex radial field dominating at mid latitudes ($\sim$45 degrees). Similar magnetic field maps are recovered from data acquired five months apart. Future work will investigate how this surface magnetic field distribution impacts the coronal magnetic field and extended environment around this planet-hosting star.

To better understand the observed distributions of rotation rate and magnetic activity of sun-like and low-mass stars, we derive a physically motivated scaling for the dependence of the stellar-wind torque on Rossby number. The torque also contains an empirically-derived scaling with stellar mass (and radius), which provides new insight into the mass-dependence of stellar magnetic and wind properties. We demonstrate that this new formulation explains why the lowest mass stars are observed to maintain rapid rotation for much longer than solar-mass stars, and simultaneously, why older populations exhibit a sequence of slowly rotating stars, in which the low-mass stars rotate more slowly than solar-mass stars. The model also reproduces some previously unexplained features in the period-mass diagram for the Kepler field, notably: the particular shape of the "upper envelope" of the distribution, suggesting that ~95% of Kepler field stars with measured rotation periods are younger than ~4 Gyr; and the shape of the "lower envelope," corresponding to the location where stars transition between magnetically saturated and unsaturated regimes.

Stellar winds are thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. 60 simulations made with a 2.5D, cylindrical and axisymmetric set-up and computed with the PLUTO code were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phase as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486- 4943-1) whose topology has been obtained by Zeeman-Doppler Imaging (ZDI).

The braking of main sequence stars originates mainly from their stellar wind. The efficiency of this angular momentum extraction depends on the rotation rate of the star, the acceleration profile of the wind and the coronal magnetic field. The derivation of scaling laws parametrizing the stellar wind torque is important for our understanding of gyro-chronology and the evolution of the rotation rates of stars. In order to understand the impact of complex magnetic topologies on the stellar wind torque, we present three- dimensional, dynamical simulations of the corona of HD 189733. Using the observed complex topology of the magnetic field, we estimate how the torque associated with the wind scales with model parameters and compare those trends to previously published scaling laws.

Stellar winds are thought to be the main process responsible for the spin down of main-sequence stars. The extraction of angular momentum by a magnetized wind has been studied for decades, leading to several formulations for the resulting torque. However, previous studies generally consider simple dipole or split monopole stellar magnetic topologies. Here we consider in addition to a dipolar stellar magnetic field, both quadrupolar and octupolar configurations, while also varying the rotation rate and the magnetic field strength. 60 simulations made with a 2.5D, cylindrical and axisymmetric set-up and computed with the PLUTO code were used to find torque formulations for each topology. We further succeed to give a unique law that fits the data for every topology by formulating the torque in terms of the amount of open magnetic flux in the wind. We also show that our formulation can be applied to even more realistic magnetic topologies, with examples of the Sun in its minimum and maximum phase as observed at the Wilcox Solar Observatory, and of a young K-star (TYC-0486-4943-1) whose topology has been obtained by Zeeman-Doppler Imaging (ZDI). Ideas about how to compute the open flux from ZDI Maps are discussed.

This paper explores and compares the pitfalls of modelling the three-dimensional wind of a spherical star with a cartesian grid. Several numerical methods are compared, using either uniform and stretched grid or adaptative mesh refinement (AMR). An additional numerical complication is added, when an orbiting planet is considered. In this case a rotating frame is added to the model such that the orbiting planet is at rest in the frame of work. The three-dimensional simulations are systematically compared to an equivalent two-dimensional, axisymmetric simulation. The comparative study presented here suggests to limit the rotation rate of the rotating frame below the rotating frame of the star and provides guidelines for further three-dimensional modelling of stellar winds in the context of close-in star-planet interactions.

Magnetic interactions between close-in planets and their host star can play an important role in the secular orbital evolution of the planets, as well as the rotational evolution of their host. As long as the planet orbits inside the Alfv\'en surface of the stellar wind, the magnetic interaction between the star and the planet can modify the wind properties and also lead to direct angular momentum transfers between the two. We model these star-planet interactions using compressible magneto-hydrodynamic (MHD) simulations, and quantify the angular momentum transfers between the star, the planet, and the stellar wind. We study the cases of magnetized and non-magnetized planets and vary the orbital radius inside the Alfv\'en surface of the stellar wind. Based on a grid of numerical simulations, we propose general scaling laws for the modification of the stellar wind torque, for the torque between the star and the planet, and for the planet migration associated with the star-planet magnetic interactions. We show that when the coronal magnetic field is large enough and the star is rotating sufficiently slowly, the effect of the magnetic star-planet interaction is comparable to tidal effects and can lead to a rapid orbital decay.

In our own solar system, the necessity of understanding space weather is readily evident. Fortunately for Earth, our nearest stellar neighbor is relatively quiet, exhibiting activity levels several orders of magnitude lower than young, solar-type stars. In protoplanetary systems, stellar magnetic phenomena observed are analogous to the solar case, but dramatically enhanced on all physical scales: bigger, more energetic, more frequent. While coronal mass ejections (CMEs) could play a significant role in the evolution of protoplanets, they could also affect the evolution of the central star itself. To assess the consequences of prominence eruption/CMEs, we have invoked the solar-stellar connection to estimate, for young, solar-type stars, how frequently stellar CMEs may occur and their attendant mass and angular momentum loss rates. We will demonstrate the necessary conditions under which CMEs could slow stellar rotation.

We model the magnetized interaction between a star and a close-in planet (SPMIs), using global, magnetohydrodynamic numerical simulations. In this proceedings, we study the effects of the numerical boundary conditions at the stellar surface, where the stellar wind is driven, and in the planetary interior. We show that is it possible to design boundary conditions that are adequate to obtain physically realistic, steady-state solutions for cases with both magnetized and unmagnetized planets. This encourages further development of numerical studies, in order to better constrain and understand SPMIs, as well as their effects on the star-planet rotational evolution.

This chapter aims at providing the most complete review of both the emerging concepts and the latest observational results regarding the angular momentum evolution of young low-mass stars and brown dwarfs. In the time since Protostars & Planets V, there have been major developments in the availability of rotation period measurements at multiple ages and in different star-forming environments that are essential for testing theory. In parallel, substantial theoretical developments have been carried out in the last few years, including the physics of the star-disk interaction, numerical simulations of stellar winds, and the investigation of angular momentum transport processes in stellar interiors. This chapter reviews both the recent observational and theoretical advances that prompted the development of renewed angular momentum evolution models for cool stars and brown dwarfs. While the main observational trends of the rotational history of low mass objects seem to be accounted for by these new models, a number of critical open issues remain that are outlined in this review.

We present 2D magnetohydrodynamic simulations performed with the PLUTO code to model magnetized star-planet interactions. We study two simple scenarios of magnetized star-planet interactions: the unipolar and dipolar interactions. Despite the simplified hypotheses we consider in the model, the qualitative behavior of the interactions is very well recovered. These encouraging results promote further developments of the model to obtain predictions on the effect and the physical manifestation of magnetized star--close-in planet interactions.

We develop an empirical model to estimate mass-loss rates via coronal mass ejections (CMEs) for solar-type pre-main-sequence (PMS) stars. Our method estimates the CME mass-loss rate from the observed energies of PMS X-ray flares, using our empirically determined relationship between solar X-ray flare energy and CME mass: log(M_CME [g]) = 0.63 x log(E_flare [erg]) - 2.57. Using masses determined for the largest flaring magnetic structures observed on PMS stars, we suggest that this solar-calibrated relationship may hold over 10 orders of magnitude in flare energy and 7 orders of magnitude in CME mass. The total CME mass-loss rate we calculate for typical solar-type PMS stars is in the range 1e-12 to 1e-9 M_sun/yr. We then use these CME mass-loss rate estimates to infer the attendant angular momentum loss leading up to the main sequence. Assuming the CME outflow rate for a typical ~1 M_sun T Tauri star is < 1e-10 M_sun/yr, the resulting spin-down torque is too small during the first ~1 Myr to counteract the stellar spin-up due to contraction and accretion. However, if the CME mass-loss rate is >1e-10 M_sun/yr, as permitted by our calculations, the CME spin-down torque may influence the stellar spin evolution after an age of a few Myr.

It is a well established fact that some YSO jets (e.g. RW Aur) display different propagation speeds between their blue and red shifted parts, a feature possibly associated with the central engine or the environment in which the jet propagates. In order to understand the origin of asymmetric YSO jet velocities, we investigate the efficiency of two candidate mechanisms, one based on the intrinsic properties of the system and one based on the role of the external medium. In particular, a parallel or anti-parallel configuration between the protostellar magnetosphere and the disk magnetic field is considered and the resulting dynamics are examined both in an ideal and a resistive magneto-hydrodynamical (MHD) regime. Moreover, we explore the effects of a potential difference in the pressure of the environment, as a consequence of the non-uniform density distribution of molecular clouds. Ideal and resistive axisymmetric numerical simulations are carried out for a variety of models, all of which are based on a combination of two analytical solutions, a disk wind and a stellar outflow. We find that jet velocity asymmetries can indeed occur both when multipolar magnetic moments are present in the star-disk system as well as when non-uniform environments are considered. The latter case is an external mechanism that can easily explain the large time scale of the phenomenon, whereas the former one naturally relates it to the YSO intrinsic properties. [abridged]

We investigate how the evolution of the stellar spin rate affects, and is affected by, planets in close orbits, via star-planet tidal interactions. To do this, we used a standard equilibrium tidal model to compute the orbital evolution of single planets orbiting both Sun-like stars and 0.1 M\odot M-dwarfs. We tested two stellar spin evolution profiles, one with fast initial rotation (P=1.2 day) and one with slow initial rotation (P=8 day). We tested the effect of varying the stellar and planetary dissipation and the planet's mass and initial orbital radius. Conclusions: Tidal evolution allows to differentiate the early behaviors of extremely close-in planets orbiting either a rapidly rotating star or a slowly rotating star. The early spin-up of the star allows the close-in planets around fast rotators to survive the early evolution. For planets around M-dwarfs, surviving the early evolution means surviving on Gyr timescales whereas for Sun-like stars the spin-down brings about late mergers of Jupiter planets. In light of this study, we can say that differentiating between one spin evolution from another given the present position of planets can be very tricky. Unless we can observe some markers of former evolution it is nearly impossible to distinguish the two very different spin profiles, let alone intermediate spin profiles. Though some conclusions can still be drawn from statistical distributions of planets around fully convective M-dwarfs. However, if the tidal evolution brings about a merger late in its history it can also entail a noticeable acceleration of the star in late ages, so that it is possible to have old stars that spin rapidly. This raises the question of better constraining the age of stars.

We use two-dimensional axisymmetric magnetohydrodynamic simulations to compute steady-state solutions for solar-like stellar winds from rotating stars with dipolar magnetic fields. Our parameter study includes 50 simulations covering a wide range of relative magnetic field strengths and rotation rates, extending from the slow- and approaching the fast-magnetic-rotator regimes. Using the simulations to compute the angular momentum loss, we derive a semi-analytic formulation for the external torque on the star that fits all of the simulations to a precision of a few percents. This formula provides a simple method for computing the magnetic braking of sun-like stars due to magnetized stellar winds, which properly includes the dependence on the strength of the magnetic field, mass loss rate, stellar radius, suface gravity, and spin rate and which is valid for both slow and fast rotators.

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 star's 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.

The stellar luminosity and depth of the convective envelope vary rapidly with mass for G- and K-type main sequence stars. In order to understand how these properties influence the convective turbulence, differential rotation, and meridional circulation, we have carried out 3D dynamical simulations of the interiors of rotating main sequence stars, using the anelastic spherical harmonic (ASH) code. The stars in our simulations have masses of 0.5, 0.7, 0.9, and 1.1 M_sun, corresponding to spectral types K7 through G0, and rotate at the same angular speed as the sun. We identify several trends of convection zone properties with stellar mass, exhibited by the simulations. The convective velocities, temperature contrast between up- and down-flows, and meridional circulation velocities all increase with stellar luminosity. As a consequence of the trend in convective velocity, the Rossby number (at a fixed rotation rate) increases and the convective turnover timescales decrease significantly with increasing stellar mass. The 3 lowest mass cases exhibit solar-like differential rotation, in a sense that they show a maximum rotation at the equator and minimum at higher latitudes, but the 1.1 M_sun case exhibits anti-solar rotation. At low mass, the meridional circulation is multi-cellular and aligned with the rotation axis; as the mass increases, the circulation pattern tends toward a unicellular structure covering each hemisphere in the convection zone.

We report on new X-ray observations of the classical T Tauri star DG Tau. DG Tau drives a collimated bi-polar jet known to be a source of X-ray emission perhaps driven by internal shocks. The rather modest extinction permits study of the jet system to distances very close to the star itself. Our initial results presented here show that the spatially resolved X-ray jet has been moving and fading during the past six years. In contrast, a stationary, very soft source much closer (~ 0.15-0.2") to the star but apparently also related to the jet has brightened during the same period. We report accurate temperatures and absorption column densities toward this source, which is probably associated with the jet base or the jet collimation region.

Classical T Tauri stars (CTTS) differ in their X-ray signatures from older pre-main sequence stars, e.g. weak-lined TTS (WTTS). CTTS show a soft excess and deviations from the low-density coronal limit in the He-like triplets. We test whether these features correlate with accretion or the presence of a disk by observing IM Lup, a disk-bearing object apparently in transition between CTTS and WTTS without obvious accretion. We analyse a Chandra grating spectrum and additional XMM-Newton data of IM Lup and accompanying optical spectra, some of them taken simultaneously to the X-ray observations. We fit the X-ray emission lines and decompose the Ha emission line in different components. In X-rays IM Lup has a bright and hot active corona, where elements of low first-ionisation potential are depleted. The He-like Ne IX triplet is in the low-density state, but due to the small number of counts a high-density scenario cannot be excluded on the 90% confidence level. In all X-ray properties IM Lup resembles a main-sequence star, but it is also compatible with CTTS signatures on the 90% confidence level, thus we cannot decide if the soft excess and deviations from the low-density coronal limit in the He-like triplets in CTTS require accretion or only the presence of a disk. IM Lup is chromospherically active, which explains most of the emission in Ha. Despite its low equivalent width, the complexity of the Ha line profile is reminiscent of CTTS. We present an estimate for the mass accretion rate of 10e-11 Msun/yr.

The Chandra Orion Ultradeep Project observed hundreds of young, low-mass stars undergoing highly energetic X-ray flare events. The 32 most powerful cases have been modeled with the result that the magnetic structures responsible for these flares can be many stellar radii in extent. In this paper, we model the observed spectral energy distributions of these 32 stars in order to determine, in detail for each star, whether there is circumstellar disk material situated in sufficient proximity to the stellar surface for interaction with the large magnetic loops inferred from the observed X-ray flares. Our spectral energy distributions span the wavelength range 0.3-8 um (plus 24 um for some stars), allowing us to constrain the presence of dusty circumstellar material out to >10 AU from the stellar surface in most cases. For 24 of the 32 stars in our sample the available data are sufficient to constrain the location of the inner edge of the dusty disks. Six of these (25%) have spectral energy distributions consistent with inner disks within reach of the observed magnetic loops. Another four stars may have gas disks interior to the dust disk and extending within reach of the magnetic loops, but we cannot confirm this with the available data. The remaining 14 stars (58%) appear to have no significant disk material within reach of the large flaring loops. Thus, up to ~40% of the sample stars exhibit energetic X-ray flares that possibly arise from a magnetic star-disk interaction, and the remainder are evidently associated with extremely large, free-standing magnetic loops anchored only to the stellar surface.

We present a model for the rotational evolution of a young, solar mass star interacting with an accretion disk. The model incorporates a description of the angular momentum transfer between the star and disk due to a magnetic connection, and includes changes in the star's mass and radius and a decreasing accretion rate. The model also includes, for the first time in a spin evolution model, the opening of the stellar magnetic field lines, as expected to arise from twisting via star-disk differential rotation. In order to isolate the effect that this has on the star-disk interaction torques, we neglect the influence of torques that may arise from open field regions connected to the star or disk. For a range of magnetic field strengths, accretion rates, and initial spin rates, we compute the stellar spin rates of pre-main-sequence stars as they evolve on the Hayashi track to an age of 3~Myr. How much the field opening affects the spin depends on the strength of the coupling of the magnetic field to the disk. For the relatively strong coupling (i.e., high magnetic Reynolds number) expected in real systems, all models predict spin periods of less than $\sim3$ days, in the age range of 1--3~Myr. Furthermore, these systems typically do not reach an equilibrium spin rate within 3~Myr, so that the spin at any given time depends upon the choice of initial spin rate. This corroborates earlier suggestions that, in order to explain the full range of observed rotation periods of approximately $1$--$10$ days, additional processes, such as the angular momentum loss from powerful stellar winds, are necessary.

Some T Tauri stars show a peculiar X-ray spectrum that can be modelled by two components with different absorbing column densities. We seek to explain the soft X-ray component in DG Tau, the best studied of these sources, with an outflow model, taking observations at other wavelengths into consideration. We constrain the outflow properties through spectral fitting and employ simple semi-analytical formulae to describe properties of a shock wave that heats up the X-ray emitting region. The X-ray emission is consistent with its arising from the fastest and innermost component of the optically detected outflow. Only a small fraction of the total mass loss is required for this X-ray emitting component. Our favoured model requires shock velocities between 400 and 500 km/s. For a density >10^5 /cm^3 all dimensions of the shock cooling zone are only a few AU, so even in optical observations this cannot be resolved. This X-ray emission mechanism in outflows may also operate in other, less absorbed T Tauri stars, in addition to corona and accretion spots.

Using numerical simulations of magnetized stellar winds, we carry out a parameter study to find the dependence of the stellar wind torque on observable parameters. We find that the power-law dependencies of the torque on parameters is significantly different than what has been used in all spin evolution models to date.

Today, the generation of magnetic fields in solar-type stars and its relation to activity and rotation can coherently be explained, although it is certainly not understood in its entirety. Rotation facilitates the generation of magnetic flux that couples to the stellar wind, slowing down the star. There are still many open questions, particularly at early phases (young age), and at very low mass. It is vexing that rotational braking becomes inefficient at the threshold to fully convective interiors, although no threshold in magnetic activity is seen, and the generation of large scale magnetic fields is still possible for fully convective stars. This article briefly outlines our current understanding of the rotation-magnetic field relation.

Models of magnetospheric accretion on to classical T Tauri stars often assume that the stellar magnetic field is a simple dipole. Recent Zeeman-Doppler imaging studies of V2129 Oph and BP Tau have shown however that their magnetic fields are more complex. V2129 Oph is a high mass T Tauri star and despite its young age is believed to have already developed a radiative core. In contrast to this, the lower mass BP Tau is likely to be completely convective. As the internal structure and therefore the magnetic field generation process is different in both stars, it is of particular interest to compare the structure of their magnetic fields obtained by field extrapolation from magnetic surface maps. We compare both field structures to mulitpole magnetic fields, and calculate the disk truncation radius for both systems. We find that by considering magnetic fields with a realistic degree of complexity, the disk is truncated at, or within, the radius obtained for dipole fields.

We present results from a low resolution (R~300) near-infrared spectroscopic variability survey of actively accreting T Tauri stars (TTS) in the Taurus-Auriga star forming region. Paschen and Brackett series H I recombination lines were detected in 73 spectra of 15 classical T Tauri systems. The values of the Pan/PaB, Brn/BrG, and BrG/Pan H I line ratios for all observations exhibit a scatter of < 20% about the weighted mean, not only from source to source, but also for epoch-to-epoch variations in the same source. A representative or `global' value was determined for each ratio in both the Paschen and Brackett series as well as the BrG/Pan line ratios. A comparison of observed line ratio values was made to those predicted by the temperature and electron density dependent models of Case B hydrogen recombination line theory. The measured line ratios are statistically well-fit by a tightly constrained range of temperatures (T < 2000 K) and electron densities 1e9 < n_e < 1e10 cm^-3. A comparison of the observed line ratio values to the values predicted by the optically thick and thin local thermodynamic equilibrium cases rules out these conditions for the emitting H I gas. Therefore, the emission is consistent with having an origin in a non-LTE recombining gas. While the range of electron densities is consistent with the gas densities predicted by existing magnetospheric accretion models, the temperature range constrained by the Case B comparison is considerably lower than that expected for accreting gas. The cooler gas temperatures will require a non-thermal excitation process (e.g., coronal/accretion-related X-rays and UV photons) to power the observed line emission.

Models of magnetospheric accretion on to classical T Tauri stars often assume that stellar magnetic fields are simple dipoles. Recently published surface magnetograms of BP Tau and V2129 Oph have shown, however, that their fields are more complex. The magnetic field of V2129 Oph was found to be predominantly octupolar. For BP Tau the magnetic energy was shared mainly between the dipole and octupole field components, with the dipole component being almost four times as strong as that of V2129 Oph. From the published surface maps of the photospheric magnetic fields we extrapolate the coronal fields of both stars, and compare the resulting field structures with that of a dipole. We consider different models where the disc is truncated at, or well-within, the Keplerian corotation radius. We find that although the structure of the surface magnetic field is particularly complex for both stars, the geometry of the larger scale field, along which accretion is occurring, is somewhat simpler. However, the larger scale field is distorted close to the star by the stronger field regions, with the net effect being that the fractional open flux through the stellar surface is less than would be expected with a dipole magnetic field model. Finally, we estimate the disc truncation radius, assuming that this occurs where the magnetic torque from the stellar magnetosphere is comparable to the viscous torque in the disc.

We compare the stellar wind torque calculated in a previous work (Paper II) to the spin-up and spin-down torques expected to arise from the magnetic interaction between a slowly rotating ($\sim 10$% of breakup) pre-main-sequence star and its accretion disk. This analysis demonstrates that stellar winds can carry off orders of magnitude more angular momentum than can be transferred to the disk, provided that the mass outflow rates are greater than the solar wind. Thus, the equilibrium spin state is simply characterized by a balance between the angular momentum deposited by accretion and that extracted by a stellar wind. We derive a semi-analytic formula for predicting the equilibrium spin rate as a function only of the ratio of $\dot M_{\rm w} / \dot M_{\rm a}$ and a dimensionless magnetization parameter, $\Psi \equiv B_*^2 R_*^2 (\dot M_{\rm a} v_{\rm esc})^{-1}$, where $\dot M_{\rm w}$ is the stellar wind mass outflow rate, $\dot M_{\rm a}$ the accretion rate, $B_*$ the stellar surface magnetic field strength, $R_*$ the stellar radius, and $v_{\rm esc}$ the surface escape speed. For parameters typical of accreting pre-main-sequence stars, this explains spin rates of $\sim 10$% of breakup speed for $\dot M_{\rm w} / \dot M_{\rm a} \sim 0.1$. Finally, the assumption that the stellar wind is driven by a fraction of the accretion power leads to an upper limit to the mass flow ratio of $\dot M_{\rm w} / \dot M_{\rm a} \la 0.6$.

[Abridged] In order to explain the slow rotation observed in a large fraction of accreting pre-main-sequence stars (CTTSs), we explore the role of stellar winds in torquing down the stars. For this mechanism to be effective, the stellar winds need to have relatively high outflow rates, and thus would likely be powered by the accretion process itself. Here, we use numerical magnetohydrodynamical simulations to compute detailed 2-dimensional (axisymmetric) stellar wind solutions, in order to determine the spin down torque on the star. We explore a range of parameters relevant for CTTSs, including variations in the stellar mass, radius, spin rate, surface magnetic field strength, the mass loss rate, and wind acceleration rate. We also consider both dipole and quadrupole magnetic field geometries. Our simulations indicate that the stellar wind torque is of sufficient magnitude to be important for spinning down a ``typical'' CTTS, for a mass loss rate of $\sim 10^{-9} M_\odot$ yr$^{-1}$. The winds are wide-angle, self-collimated flows, as expected of magnetic rotator winds with moderately fast rotation. The cases with quadrupolar field produce a much weaker torque than for a dipole with the same surface field strength, demonstrating that magnetic geometry plays a fundamental role in determining the torque. Cases with varying wind acceleration rate show much smaller variations in the torque suggesting that the details of the wind driving are less important. We use our computed results to fit a semi-analytic formula for the effective Alfv\'en radius in the wind, as well as the torque. This allows for considerable predictive power, and is an improvement over existing approximations.

Stellar winds may be important for angular momentum transport from accreting T Tauri stars, but the nature of these winds is still not well-constrained. We present some simulation results for hypothetical, hot (~1e6 K) coronal winds from T Tauri stars, and we calculate the expected emission properties. For the high mass loss rates required to solve the angular momentum problem, we find that the radiative losses will be much greater than can be powered by the accretion process. We place an upper limit to the mass loss rate from accretion-powered coronal winds of ~1e-11 solar masses per year. We conclude that accretion powered stellar winds are still a promising scenario for solving the stellar angular momentum problem, but the winds must be cool (~1e4 K) and thus are not driven by thermal pressure.

We review the theoretical efforts to understand why pre-main-sequence stars spin much more slowly than expected. The first idea put forward was that massive stellar winds may remove substantial angular momentum. Since then, it has become clear that the magnetic interaction between the stars and their accretion disks explains many of the observed emission properties. The disk locking scenario, which assumes the magnetic star-disk interaction also solves the stellar spin problem, has received the most attention in the literature to date. However, recent considerations suggest that the torques in the star-disk interaction are insufficient for disk locking to explain the slow rotators. This prompts us to revisit stellar winds, and we conclude that stellar winds, working in conjunction with magnetospheric accretion, are a promising candidate for solving the angular momentum problem. We suggest future directions for both observations and theory, to help shed light on this issue.

We present the results of a numerical magnetohydrodynamic simulation that demonstrates a mechanism by which magnetic fields tap rotational energy of a stellar core and expel the envelope. Our numerical setup, designed to focus on the basic physics of the outflow mechanism, consists of a solid, gravitating sphere, which may represent the compact core of a star, surrounded by an initially hydrostatic envelope of ionized gas. The core is threaded by a dipolar magnetic field that also permeates the envelope. At the start of the simulation, the core begins to rotate at 10% of the escape speed. The magnetic field is sufficiently strong to drive a magneto-rotational explosion, whereby the entire envelope is expelled, confirming the expectation of analytical models. Furthermore, the dipolar nature of the field results in an explosion that is enhanced simultaneously along the rotation axis (a jet) and along the magnetic equator. While the initial condition is simplified, the simulation approximates circumstances that may arise in astrophysical objects such as Type II supernovae, gamma ray bursts, and proto-planetary nebulae.

We compare the angular momentum extracted by a wind from a pre-main-sequence star to the torques arising from the interaction between the star and its Keplerian accretion disk. We find that the wind alone can counteract the spin-up torque from mass accretion, solving the mystery of why accreting pre-main-sequence stars are observed to spin at less than 10% of break-up speed, provided that the mass outflow rate in the stellar winds is ~10% of the accretion rate. We suggest that such massive winds will be driven by some fraction $\epsilon$ of the accretion power. For observationally constrained typical parameters of classical T-Tauri stars, $\epsilon$ needs to be between a few and a few tens of percent. In this scenario, efficient braking of the star will terminate simultaneously with accretion, as is usually assumed to explain the rotation velocities of stars in young clusters.

We formulate a general, steady-state model for the torque on a magnetized star from a surrounding accretion disc. For the first time, we include the opening of dipolar magnetic field lines due to the differential rotation between the star and disc, so the magnetic topology then depends on the strength of the magnetic coupling to the disc. This coupling is determined by the effective slip rate of magnetic field lines that penetrate the diffusive disc. Stronger coupling (i.e., lower slip rate) leads to a more open topology and thus to a weaker magnetic torque on the star from the disc. In the expected strong coupling regime, we find that the spin-down torque on the star is more than an order of magnitude smaller than calculated by previous models. We also use our general approach to examine the equilibrium (`disc-locked') state, in which the net torque on the star is zero. In this state, we show that the stellar spin rate is roughly an order of magnitude faster than predicted by previous models. This challenges the idea that slowly-rotating, accreting protostars are disc locked. Furthermore, when the field is sufficiently open (e.g., for mass accretion rates > 5 x 10^{-9} M_sun / yr, for typical accreting protostars), the star will receive no magnetic spin-down torque from the disc at all. We therefore conclude that protostars must experience a spin-down torque from a source that has not yet been considered in the star-disc torque models--possibly from a stellar wind along the open field lines.

The luminous blue variable Eta Carinae is surrounded by a complex and highly structured nebula of ejected material. The best-studied and axisymmetric components of this outflow consist of bipolar lobes (the ``homunculus'') and an equatorial ``skirt.'' Recent proper motion measurements suggest that the skirt was ejected at the same time as the lobes, contrary to the assumptions of all current theoretical models for the formation of the nebula (which use the skirt to collimate stellar winds into lobes). We present a magnetohydrodynamic (MHD) stellar wind model that produces an outflowing disk and bipolar lobes in a single, steady-state wind. The basic model consists of a wind from a rotating star with a rotation-axis-aligned dipole magnetic field. The azimuthal component of the magnetic field, generated by stellar rotation, compresses the wind toward the equator and also toward the rotation axis, simultaneously producing an outflowing disk and jet. We use numerical MHD simulations to study the wind for various amounts of stellar rotation and to show a range of wind morphologies. In order to produce wide angle lobes similar to the homunculus (which have roughly a 30 degree opening angle), a high-speed polar wind (with enhanced energy density) from the star is also required. In that case, the structure of the wind bears a remarkable resemblance to the skirt plus homunculus morphology of the Eta Car nebulae, and a significant fraction of the stellar angular momentum is carried away by the wind. Although the model assumes a steady-state wind (rather than an eruption) and thermal wind driving (rather than radiation pressure), the structure of the wind is encouraging.

We critically examine the theory of disk locking, which assumes that the angular momentum deposited on an accreting protostar is exactly removed by torques carried along magnetic field lines connecting the star to the disk. In this letter, we consider that the differential rotation between the star and disk naturally leads to an opening (i.e., disconnecting) of the magnetic field between the two. We find that this significantly reduces the spin-down torque on the star by the disk. Thus, disk-locking cannot account for the slow rotation ($\sim$10% of breakup speed) observed in several systems and for which the model was originally developed.

In this talk, we take a new look at the theory of disk locking, which assumes that an accreting protostar rids itself of accreted angular momentum through a magnetic coupling with the accretion disk. We consider that differential rotation between the star and disk twists the field lines. For large enough twist, the magnetic field lines connecting the star and disk open and disconnect. This significantly reduces the spin-down torque on the star by the disk, and so we find that disk-locking theory predicts spin periods that are much too short to account for typical observed systems.

We discuss a mechanism by which a giant star can expel its envelope in an outburst, leaving its core exposed. The outburst is powered by rotational kinetic energy of the core, transferred to the envelope via the twisting of magnetic fields. We show that, if the core is magnetized, and if it has sufficient angular momentum, this mechanism may be triggered at the end of the asymptotic giant branch phase, and drive a proto-planetary nebula (pPN) outflow. This explosion of magnetic energy self-consistently explains some of the asymmetries and dynamics of pPNe.

Studies of jets from young stellar objects (YSO's) suggest that material is launched from a small central region at wide opening angles and collimated by an interaction with the surrounding environment. Using time-dependent, numerical magnetohydrodynamic simulations, we follow the detailed launching of a central wind via the coupling of a stellar dipole field to the inner edge of an accretion disc. Our method employs a series of nested computational grids, which allows the simulations to follow the central wind out to scales of tens of AU, where it may interact with its surroundings. The coupling between the stellar magnetosphere and disc inner edge has been known to produce an outflow containing both a highly collimated jet plus a wide-angle flow. The jet and wide-angle wind flow at roughly the same speed (100--200 km/s), and most of the energy and mass is carried off at relatively wide angles. We show that the addition of a weak disc-associated field (<< 0.1 Gauss) has little effect on the wind launching, but it collimates the entire flow (jet + wide wind) at a distance of several AU. The collimation is inevitable, regardless of the relative polarity of the disc field and stellar dipole, and the result is a more powerful and physically broader collimated flow than from the star-disc interaction alone. Within the collimation region, the morphology of the large-scale flow resembles a pitchfork, in projection. We compare these results with observations of outflows from YSO's and discuss the possible origin of the disc-associated field.

We present the results of time-dependent, numerical magnetohydrodynamic simulations of a realistic young stellar object outflow model with the addition of a disk-associated magnetic field. The outflow produced by the magnetic star-disk interaction consists of an episodic jet plus a wide-angle wind with an outflow speed comparable to that of the jet (100--200 km/s). An initially vertical field of << 0.1 Gauss, embedded in the disk, has little effect on the wind launching mechanism, but we show that it collimates the entire flow (jet + wide wind) at large (several AU) distances. The collimation does not depend on the polarity of the vertical field. We also discuss the possible origin of the disk-associated field.

To gain insight into the nature of the peculiar Herbig-Haro object HH 255 (also called Burnham's nebula), we use previously published observations to derive information about the emission line fluxes as a function of position within HH 255 and compare them with the well-studied, and relatively well-behaved bow shock HH 1. There are some qualitative similarities in the H$\alpha$ and [O III] 5007 lines in both objects. However, in contrast to the expectation of the standard bow shock model, the fluxes of the [O I] 6300, [S II] 6731, and [N II] 6583 lines are essentially constant along the axis of the flow, while the electron density decreases, over a large distance within HH 255. We also explore the possibility that HH 255 represents the emission behind a standing or quasi-stationary shock. The shock faces upwind, and we suggest, using theoretical arguments, that it may be associated with the collimation of the southern outflow from T Tauri. Using a simplified magnetohydrodynamic simulation to illustrate the basic concept, we demonstrate that the existence of such a shock at the north edge of HH 255 could indeed explain its unusual kinematic and ionization properties. Whether or not such a shock can explain the detailed emission line stratification remains an open question.

We examine, parametrically, the interaction between the magnetosphere of a rotating, young stellar object (YSO) and a circumstellar accretion disk using 2.5-D (cylindrically symmetric) numerical magnetoydrodynamic simulations. The interaction drives a collimated outflow, and we find that the jet formation mechanism is robust. For variations in initial disk density of a factor of 16, variations of stellar dipole strength of a factor of 4, and for various initial conditions with respect to the disk truncation radius and the existence of a disk field, outflows with similar morphologies were consistently produced. Secondly, the system is self-regulating, where the outflow properties depend relatively weakly on the parameters above. The large scale magnetic field structure rapidly evolves to a configuration that removes angular momentum from the disk at a rate that depends most strongly on the field and weakly on the rotation rate of the foot-points of the field in the disk and the mass outflow rate. Third, the simulated jets are episodic, with the timescale of jet outbursts identical to the timescale of magnetically induced oscillations of the inner edge of the disk. To better understand the physics controlling these disk oscillations, we present a semi-analytical model and confirm that the oscillation period is set by the spin down rate of the disk inner edge. Finally, our simulations offer strong evidence that it is indeed the interaction of the stellar magnetosphere with the disk, rather than some primordial field in the disk itself, that is responsible for the formation of jets from these systems.

(modified for brevity) It is important to determine whether the observed bow shocks in the working surfaces of Herbig-Haro outflows has lead to a destruction of dust grains and consequently to a change in the gas phase metal abundances (say of Fe) in the cooling regions of HH bow shocks. Detailed studies are currently available for only 5 HH outflows (Beck-Winchatz et al.\ 1996), due to the large observational and theoretical effort required to determine metal abundances in HH objects. Information about metal abundances in more HH objects is badly needed. We therefore use a very approximate method by introducing a "characteristic number", A_me, whose definition is based only on often observed line fluxes. We find a good correlation between A_me and the Fe abundance for the 5 well-studied HH objects. We use this correlation to determine approximate values of the gas phase Fe abundance in 13 additional high excitation and in 3 additional low excitation HH objects. Of the 16 high excitation HH objects studied, there are 6 which approximately agree with normal population I abundance (i.e., no depletion due to dust formation). The remaining 10 show some Fe gas phase depletion which never gets larger than by a factor of 2.5. This result is in agreement with our qualitative expectations that fast shocks efficiently destroy dust grains. Of the 5 low excitation HH objects studied, there are 4 which, unexpectedly, show a normal population I abundance. In low excitation objects, one might expect strong gas phase Fe depletion (showing the unchanged molecular cloud composition), unless the matter has previously gone through shocks of much higher shock velocities.

We present a simple, robust mechanism by which an isolated star can produce an equatorial disk. The mechanism requires that the star have a simple dipole magnetic field on the surface and an isotropic wind acceleration mechanism. The wind couples to the field, stretching it until the field lines become mostly radial and oppositely directed above and below the magnetic equator, as occurs in the solar wind. The interaction between the wind plasma and magnetic field near the star produces a steady outflow in which magnetic forces direct plasma toward the equator, constructing a disk. In the context of a slow (10 km/s) outflow (10^{-5} M_sun/yr) from an AGB star, MHD simulations demonstrate that a dense equatorial disk will be produced for dipole field strengths of only a few Gauss on the surface of the star. A disk formed by this model can be dynamically important for the shaping of Planetary Nebulae.