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.

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 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.

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.