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

Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging) and the stellar radius. We use a sample of 51 stars across a mass range of 0.1-1.34 M$_\odot$ to show how the helicity density relates to stellar mass, Rossby number, magnetic energy and age. We find that the large-scale helicity density increases with decreasing Rossby number $R_o$, peaking at $R_o \simeq 0.1$, with a saturation or decrease below that. For both fully- and partially-convective stars we find that the mean absolute helicity density scales with the mean squared toroidal magnetic flux density according to the power law: $|\langle{h\,}\rangle|$ $\propto$ $\langle{\rm{B_{tor}}^2_{}\,\rangle}^{0.86\,\pm\,0.04}$. The scatter in this relation is consistent with the variation across a solar cycle, which we compute using simulations and observations across solar cycles 23 and 24 respectively. We find a significant decrease in helicity density with age.

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

Although the present-day Sun rotates too slowly to exhibit centrifugally-supported 'slingshot prominences', at some time during its past it may have formed these clouds of cool gas and ejected them into the interplanetary medium. We determine the time period for this behaviour by using a rotation evolution code to derive the properties of the formation and ejection of slingshot prominences during the lifetime of a star similar to our Sun. The mass, mass loss rate and rate of ejection of these prominences are calculated using the analytical expression derived in our previous work. We find that for stars with an initial rotation rate larger than $4.6\, \Omega_\odot$, about half of all solar mass stars, slingshot prominences will be present even after the star reaches the main sequence phase. In a fast rotator, this means that prominences can form until the star reaches $\sim 800$ Myr old. Our results also indicate that the mass and lifetime of this type of prominence have maximum values when the star reaches the ZAMS at an age of $\sim 40$ Myr for a solar mass star.

The behaviour of the large-scale dipolar field, during a star's magnetic cycle, can provide valuable insight into the stellar dynamo and associated magnetic field manifestations such as stellar winds. We investigate the temporal evolution of the dipolar field of the K dwarf 61 Cyg A using spectropolarimetric observations covering nearly one magnetic cycle equivalent to two chromospheric activity cycles. The large-scale magnetic field geometry is reconstructed using Zeeman Doppler imaging, a tomographic inversion technique. Additionally, the chromospheric activity is also monitored. The observations provide an unprecedented sampling of the large-scale field over a single magnetic cycle of a star other than the Sun. Our results show that 61 Cyg A has a dominant dipolar geometry except at chromospheric activity maximum. The dipole axis migrates from the southern to the northern hemisphere during the magnetic cycle. It is located at higher latitudes at chromospheric activity cycle minimum and at middle latitudes during cycle maximum. The dipole is strongest at activity cycle minimum and much weaker at activity cycle maximum. The behaviour of the large-scale dipolar field during the magnetic cycle resembles the solar magnetic cycle. Our results are further confirmation that 61 Cyg A indeed has a large-scale magnetic geometry that is comparable to the Sun's, despite being a slightly older and cooler K dwarf.

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.

One of the aims of the BCool programme is to search for cycles in other stars and to understand how similar they are to the Sun. In this paper we aim to monitor the evolution of $\tau$ Boo's large-scale magnetic field using high-cadence observations covering its chromospheric activity maximum. For the first time, we detect a polarity switch that is in phase with $\tau$ Boo's 120 day chromospheric activity maximum and its inferred X-ray activity cycle maximum. This means that $\tau$ Boo has a very fast magnetic cycle of only 240 days. At activity maximum $\tau$ Boo's large-scale field geometry is very similar to the Sun at activity maximum: it is complex and there is a weak dipolar component. In contrast, we also see the emergence of a strong toroidal component which has not been observed on the Sun, and a potentially overlapping butterfly pattern where the next cycle begins before the previous one has finished.

The observational signatures of prominences have been detected in single and binary G and K type stars for many years now, but recently this has been extended to the M dwarf regime. Prominences carry away both mass and angular momentum when they are ejected and the impact of this mass on any orbiting planets may be important for the evolution of exoplanetary atmospheres. By means of the classification used in the massive star community, that involves knowledge of two parameters (the co-rotation and Alfv\'en radii, $r_K$ and $r_A$), we have determined which cool stars could support prominences. From a model of mechanical support, we have determined that the prominence mass $m_p/M_\star=(E_M/E_G)(r_\star/r_K)^2 F$ where $E_MB_\star^2r_\star^3$ and $E_G = GM_\star^2/r_\star$ are magnetic and gravitational energies and $F$ is a geometric factor. Our calculated masses and ejection frequencies (typically $10^{16}-10^{17}$g and 0.4 d, respectively) are consistent with observations and are sufficient to ensure that an exoplanet orbiting in the habitable zone of an M dwarf could suffer frequent impacts.

Magnetic activity is known to be correlated to the rotation period for moderately active main sequence solar-like stars. In turn, the stellar rotation period evolves as a result of magnetised stellar winds that carry away angular momentum. Understanding the interplay between magnetic activity and stellar rotation is therefore a central task for stellar astrophysics. Angular momentum evolution models typically employ spin-down torques that are formulated in terms of the surface magnetic field strength. However, these formulations fail to account for the magnetic field geometry, unlike those that are expressed in terms of the open flux, i.e. the magnetic flux along which stellar winds flow. In this work, we model the angular momentum evolution of main sequence solar-mass stars using a torque law formulated in terms of the open flux. This is done using a potential field source surface model in conjunction with the Zeeman-Doppler magnetograms of a sample of roughly solar-mass stars. We explore how the open flux of these stars varies with stellar rotation and choice of source surface radii. We also explore the effect of field geometry by using two methods of determining the open flux. The first method only accounts for the dipole component while the second accounts for the full set of spherical harmonics available in the Zeeman-Doppler magnetogram. We find only a small difference between the two methods, demonstrating that the open flux, and indeed the spin-down, of main sequence solar-mass stars is likely dominated by the dipolar component of the magnetic field.

Magnetic activity and rotation are known to be intimately linked for low-mass stars. Understanding rotation evolution over the stellar lifetime is therefore an important goal within stellar astrophysics. In recent years, there has been increased focus on how the complexity of the stellar magnetic field affects the rate of angular momentum-loss from a star. This is a topic that Zeeman-Doppler imaging (ZDI), a technique that is capable of reconstructing the large-scale magnetic field topology of a star, can uniquely address. Using a potential field source surface model, we estimate the open flux, mass loss-rate and angular momentum-loss rates for a sample of 66 stars that have been mapped with ZDI. We show that the open flux of a star is predominantly determined by the dipolar component of its magnetic field for our choice of source surface radius. We also show that, on the main sequence, the open flux, mass- and angular momentum-loss rates increase with decreasing Rossby number. The exception to this rule is stars less massive than $0.3M_{\odot}$. Previous work suggests that low mass M dwarfs may possess either strong, ordered and dipolar fields or weak and complex fields. This range of field strengths results in a large spread of angular momentum-loss rates for these stars and has important consequences for their spin down behaviour. Additionally, our models do not predict a transition in the mass-loss rates at the so called wind dividing line noted from Ly$\alpha$ studies.

Stellar winds govern the angular momentum evolution of solar-like stars throughout their main-sequence lifetime. The efficiency of this process depends on the geometry of the star's magnetic field. There has been a rapid increase recently in the number of stars for which this geometry can be determined through spectropolarimetry. We present a computationally efficient method to determine the 3D geometry of the stellar wind and to estimate the mass loss rate and angular momentum loss rate based on these observations. Using solar magnetograms as examples, we quantify the extent to which the values obtained are affected by the limited spatial resolution of stellar observations. We find that for a typical stellar surface resolution of 20$^{\rm o}$-30$^{\rm o}$, predicted wind speeds are within 5$\%$ of the value at full resolution. Mass loss rates and angular momentum loss rates are within 5-20$\%$. In contrast, the predicted X-ray emission measures can be under-estimated by 1-2 orders of magnitude, and their rotational modulations by 10-20$\%$.

The magnetic geometry of the surface magnetic fields of more than 55 cool stars have now been mapped using spectropolarimetry. In order to better understand these observations, we compare the magnetic field topology at different surface scale sizes of observed and simulated cool stars. For ease of comparison between the high-resolution non-potential magnetofrictional simulations and the relatively low-resolution observations, we filter out the small-scale field in the simulations using a spherical harmonics decomposition. We show that the large-scale field topologies of the solar-based simulations produce values of poloidal/toroidal fields and fractions of energy in axisymmetric modes that are similar to the observations. These global non-potential evolution model simulations capture key magnetic features of the observed solar-like stars through the processes of surface flux transport and magnetic flux emergence. They do not, however, reproduce the magnetic field of M-dwarfs or stars with dominantly toroidal field. Furthermore, we analyse the magnetic field topologies of individual spherical harmonics for the simulations and discover that the dipole is predominately poloidal, while the quadrupole shows the highest fraction of toroidal fields. Magnetic field structures smaller than a quadrupole display a fixed ratio between the poloidal and toroidal magnetic energies.

Zeeman Doppler imaging has successfully mapped the large-scale magnetic fields of stars over a large range of spectral types, rotation periods and ages. When observed over multiple epochs, some stars show polarity reversals in their global magnetic fields. On the Sun, polarity reversals are a feature of its activity cycle. In this paper, we examine the magnetic properties of stars with existing chromospherically determined cycle periods. Previous authors have suggested that cycle periods lie on multiple branches, either in the cycle period-Rossby number plane or the cycle period-rotation period plane. We find some evidence that stars along the active branch show significant average toroidal fields that exhibit large temporal variations while stars exclusively on the inactive branch remain dominantly poloidal throughout their entire cycle. This lends credence to the idea that different shear layers are in operation along each branch. There is also evidence that the short magnetic polarity switches observed on some stars are characteristic of the inactive branch while the longer chromospherically determined periods are characteristic of the active branch. This may explain the discrepancy between the magnetic and chromospheric cycle periods found on some stars. These results represent a first attempt at linking global magnetic field properties obtained form ZDI and activity cycles.

Wood et al suggested that mass-loss rate is a function of X-ray flux ($\dot{M} \propto F_x^{1.34}$) for dwarf stars with $F_x \lesssim F_{x,6} \equiv 10^6$ erg cm$^{-2}$ s$^{-1}$. However, more active stars do not obey this relation. These authors suggested that the break at $F_{x,6}$ could be caused by significant changes in magnetic field topology that would inhibit stellar wind generation. Here, we investigate this hypothesis by analysing the stars in Wood et al's sample that had their surface magnetic fields reconstructed through Zeeman-Doppler Imaging (ZDI). Although the solar-like outliers in the $\dot{M}$ -- $F_x$ relation have higher fractional toroidal magnetic energy, we do not find evidence of a sharp transition in magnetic topology at $F_{x,6}$. To confirm this, further wind measurements and ZDI observations at both sides of the break are required. As active stars can jump between states with highly toroidal to highly poloidal fields, we expect significant scatter in magnetic field topology to exist for stars with $F_x \gtrsim F_{x,6}$. This strengthens the importance of multi-epoch ZDI observations. Finally, we show that there is a correlation between $F_x$ and magnetic energy, which implies that $\dot{M}$ -- magnetic energy relation has the same qualitative behaviour as the original $\dot{M}$ -- $F_x$ relation. No break is seen in any of the $F_x$ -- magnetic energy

Spectropolarimetric observations have been used to map stellar magnetic fields, many of which display strong bands of azimuthal fields that are toroidal. A number of explanations have been proposed to explain how such fields might be generated though none are definitive. In this paper, we examine the toroidal fields of a sample of 55 stars with magnetic maps, with masses in the range 0.1-1.5$\,{\rm M}_\odot$. We find that the energy contained in toroidal fields has a power law dependence on the energy contained in poloidal fields. However the power index is not constant across our sample, with stars less and more massive than 0.5$\,{\rm M}_\odot$ having power indices of 0.72$\pm$0.08 and 1.25$\pm$0.06 respectively. There is some evidence that these two power laws correspond to stars in the saturated and unsaturated regimes of the rotation-activity relation. Additionally, our sample shows that strong toroidal fields must be generated axisymmetrically. The latitudes at which these bands appear depend on the stellar rotation period with fast rotators displaying higher latitude bands than slow rotators. The results in this paper present new constraints for future dynamo studies.

We investigate the variability of exoplanetary radio emission using stellar magnetic maps and 3D field extrapolation techniques. We use a sample of hot Jupiter hosting stars, focusing on the HD 179949, HD 189733 and tau Boo systems. Our results indicate two time-scales over which radio emission variability may occur at magnetised hot Jupiters. The first is the synodic period of the star-planet system. The origin of variability on this time-scale is the relative motion between the planet and the interplanetary plasma that is co-rotating with the host star. The second time-scale is the length of the magnetic cycle. Variability on this time-scale is caused by evolution of the stellar field. At these systems, the magnitude of planetary radio emission is anticorrelated with the angular separation between the subplanetary point and the nearest magnetic pole. For the special case of tau Boo b, whose orbital period is tidally locked to the rotation period of its host star, variability only occurs on the time-scale of the magnetic cycle. The lack of radio variability on the synodic period at tau Boo b is not predicted by previous radio emission models, which do not account for the co-rotation of the interplanetary plasma at small distances from the star.

Context: The principle definition of habitability for exoplanets is whether they can sustain liquid water on their surfaces, i.e. that they orbit within the habitable zone. However, the planet's magnetosphere should also be considered, since without it, an exoplanet's atmosphere may be eroded away by stellar winds. Aims: The aim of this paper is to investigate magnetospheric protection of a planet from the effects of stellar winds from solar-mass stars. Methods: We study hypothetical Earth-like exoplanets orbiting in the host star's habitable zone for a sample of 124 solar-mass stars. These are targets that have been observed by the Bcool collaboration. Using two wind models, we calculate the magnetospheric extent of each exoplanet. These wind models are computationally inexpensive and allow the community to quickly estimate the magnetospheric size of magnetised Earth-analogues orbiting cool stars. Results: Most of the simulated planets in our sample can maintain a magnetosphere of ~5 Earth radii or larger. This suggests that magnetised Earth analogues in the habitable zones of solar analogues are able to protect their atmospheres and is in contrast to planets around young active M dwarfs. In general, we find that Earth-analogues around solar-type stars, of age 1.5 Gyr or older, can maintain at least a Paleoarchean Earth sized magnetosphere. Our results indicate that planets around 0.6 - 0.8 solar-mass stars on the low activity side of the Vaughan-Preston gap are the optimum observing targets for habitable Earth analogues.

Under sufficiently high electric field gradients, electron behaviour within exactly perpendicular shocks is unstable to the so-called trajectory instability. We extend previous work paying special attention to shortiscale, high amplitude structures as observed within the electric field profile. Via test particle simulations, we show that such structures can cause the electron distribution to heat in a manner that violates conservation of the first adiabatic invariant. This is the case even if the overall shock width is larger than the upstream electron gyroradius. The spatial distance over which these structures occur therefore constitutes a new scale length relevant to the shock heating problem. Furthermore, we find that the spatial location of the short-scale structure is important in determining the total effect of non-adiabatic behaviour - a result that has not been previously noted.