Search

Neutron stars cool down during their lifetime through the combination of neutrino emission from the interior and photon cooling from the surface. Strongly magnetised neutron stars, called magnetars, are no exception, but the effect of their strong fields adds further complexities to the cooling theory. Besides other factors, modelling the outermost hundred meters (the envelope) plays a crucial role in predicting their surface temperatures. In this letter, we revisit the influence of envelopes on the cooling properties of neutron stars, with special focus on the critical effects of the magnetic field. We explore how our understanding of the relation between the internal and surface temperatures has evolved over the past two decades, and how different assumptions about the neutron star envelope and field topology lead to radically different conclusions on the surface temperature and its cooling with age. In particular, we find that relatively old magnetars with core-threading magnetic fields are actually much cooler than a rotation-powered pulsar of the same age. This is at variance with what is typically observed in crustal-confined models. Our results have important implications for the estimates of the X-ray luminosities of aged magnetars, and the subsequent population study of the different neutron star classes., Comment: 5 pages, 4 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The nearby M dwarf WX UMa has recently been detected at radio wavelengths with LOFAR. The combination of its observed brightness temperature and circular polarisation fraction suggests that the emission is generated via the electron-cyclotron maser instability. Two distinct mechanisms have been proposed to power such emission from low-mass stars: either a sub-Alfv\'enic interaction between the stellar magnetic field and an orbiting planet, or reconnection at the edge of the stellar magnetosphere. In this paper, we investigate the feasibility of both mechanisms, utilising the information about the star's surrounding plasma environment obtained from modelling its stellar wind. Using this information, we show that a Neptune-sized exoplanet with a magnetic field strength of 10-100 G orbiting at ~0.034 au can accurately reproduce the observed radio emission from the star, with corresponding orbital periods of 7.4 days. Due to the stellar inclination, a planet in an equatorial orbit is unlikely to transit the star. While such a planet could induce radial velocity semi-amplitudes from 7 to 396 m s$^{-1}$, it is unlikely that this signal could be detected with current techniques due to the activity of the host star. The application of our planet-induced radio emission model here illustrates its exciting potential as a new tool for identifying planet-hosting candidates from long-term radio monitoring. We also develop a model to investigate the reconnection-powered emission scenario. While this approach produces less favourable results than the planet-induced scenario, it nevertheless serves as a potential alternative emission mechanism which is worth exploring further., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

We present the first 3D fully coupled magneto-thermal simulations of neutron stars (including the most realistic background structure and microphysical ingredients so far) applied to a very complex initial magnetic field topology in the crust, similar to what recently obtained by proto-neutron star dynamo simulations. In such configurations, most of the energy is stored in the toroidal field, while the dipolar component is a few percent of the mean magnetic field. This initial feature is maintained during the long-term evolution (1e6 yr), since the Hall term favours a direct cascade (compensating for Ohmic dissipation) rather than a strong inverse cascade, for such an initial field topology. The surface dipolar component, responsible for the dominant electromagnetic spin-down torque, does not show any increase in time, when starting from this complex initial topology. This is at contrast with the timing properties of young pulsars and magnetars which point to higher values of the surface dipolar fields. A possibility is that the deep-seated magnetic field (currents in the core) is able to self-organize in large scales (during the collapse or in the early life of a neutron star). Alternatively, the dipolar field might be lower than is usually thought, with magnetosphere substantially contributing to the observed high spin-down, via e.g., strong winds or strong coronal magnetic loops, which can also provide a natural explanation to the tiny surface hotspots inferred from X-ray data., 8 pages, 6 figures, submitted for publication in MNRAS, comments are welcome

We present a spectropolarimetric analysis of the hot star V352Peg. We have acquired 18 spectropolarimetric observations of the star with ESPaDOnS at the CFHT between 2018 and 2019 and completed our dataset with one archival ESPaDOnS measurement obtained in 2011. Our analysis of the spectra shows that the star is on the main sequence and chemically peculiar, i.e. it is a Bp star, with overabundances of iron peak elements (Ti, Cr and Fe) and underabundance of He and O. Through a Least-Square Deconvolution of each spectrum, we extracted the mean Zeeman signature and mean line profile of thousands of spectral lines and detected a magnetic field in V352Peg. By modelling the Stokes I and V profiles and using the Oblique Rotator Model, we determined the geometrical configuration of V352Peg. We also performed Zeeman-Doppler Imaging (ZDI) to provide a more detailed characterization of the magnetic field of V352Peg and its surface chemical distributions. We find a magnetic field that is mainly dipolar, dominantly poloidal, and largely non-axisymmetric with a dipole field strength of $\sim$9 kG and a magnetic axis almost perpendicular to the rotation axis. The strong variability of Stokes I profiles also suggests the presence of chemical spots at the stellar surface., 18 pages, 15 figures. Accepted for publication in MNRAS

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M$_\odot$), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly-rotating, nitrogen-enriched ("Group 2") stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of i) revisiting the derived stellar parameters of known magnetic stars, and ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes., Comment: Accepted for publication in MNRAS. A full reproduction package is available on Zenodo at https://doi.org/10.5281/zenodo.7069766

We present the study on continuous high-resolution spectroscopic observations of two long-period single-lined RS Canum Venaticorum (RS CVn) binary stars IM Pegasi (IM Peg) and σ Geminorum (σ Gem), obtained with the Hertzsprung SONG telescope during the 2015–2016 season. Chromospheric activity indicators Hα, $\rm{Na\,\,{\small I}}$ D1, D2 doublet, $\rm{He\,\,{\small I}}$ D3, and Hβ lines have been analysed by using the spectral subtraction technique. The expected chromospheric emission features in the Hα, $\rm{Na\,\,{\small I}}$ D1, D2 doublet, and Hβ lines confirm that both of two stars are very active systems. In the spectra, the $\rm{He\,\,{\small I}}$ D3 line had been always detected in absorption feature. Although the behaviour of chromospheric activity indicators is very similar for both stars, the activity level of IM Peg is much stronger than that of σ Gem. Moreover, the equivalent width variations of the Hα, $\rm{He\,\,{\small I}}$ D3, and Hβ line subtractions correlate well and show different behaviour among different orbital cycles, which indicates the presence and evolution of activity longitudes over the surface of two stars. Furthermore, the subtracted Hα line profile is usually asymmetric. The red-shifted excess absorption features could be interpreted as a strong down-flow of cool absorbing material, while the blue-shifted emission component is probably caused by up-flow of hot materials through microflare events.

Funding: The authors acknowledge support from STFC consolidated grant number ST/R000824/1. Magnetic confinement of dense plasma is found in the magnetospheres of both high and low mass stars. Trapped material traces the magnetic field structure, often at large distances from the star where the magnetic structure is otherwise difficult to observe. This work looks specifically at rapidly rotating, solar-like stars where this behaviour is well observed in the form of “slingshot” prominences. We have produced a model for generating cooled magnetic loops in equilibrium with a range of coronal magnetic fields. These loops can be used to populate model coronae and confine material at a wide range of heights above the stellar surface. We calculate masses for slingshot prominences for the star AB Doradus that are consistent with observational values. The model produces two types of solution: loops with summits at low heights and tall solutions beyond the co-rotation radius. We show that the low-lying solutions are footpoint heavy and generally follow the shape of the background field. We refer to these as solar-like prominences. The tall solutions are summit heavy and are centrifugally supported. These are are the slingshot prominences. These tall solutions can be found within the stellar wind, beyond the closed corona. Hα trails are generated for various coronal field structures with a range of field geometries and coronal extents. Similar Hα trails are produced by a range of global field structures, which implies that magnetic confinement of material should be common in rapidly rotating stars. Publisher PDF

Stellar dynamos generate magnetic fields that are of fundamental importance to the variability and evolution of Sun-like and low-mass stars, and for the development of their planetary systems. As a key to understanding stellar dynamos, empirical relations between stellar parameters and magnetic fields are required for comparison to ab initio predictions from dynamo models. We report measurements of surface-average magnetic fields in 292 M dwarfs from a comparison with radiative transfer calculations; for 260 of them, this is the first measurement of this kind. Our data were obtained from more than 15 000 high-resolution spectra taken during the CARMENES project. They reveal a relation between average field strength, ⟨B⟩, and Rossby number, Ro, resembling the well-studied rotation–activity relation. Among the slowly rotating stars, we find that magnetic flux, ΦB, is proportional to rotation period, P, and among the rapidly rotating stars that average surface fields do not grow significantly beyond the level set by the available kinetic energy. Furthermore, we find close relations between nonthermal coronal X-ray emission, chromospheric Hα and Ca H&K emission, and magnetic flux. Taken together, these relations demonstrate empirically that the rotation–activity relation can be traced back to a dependence of the magnetic dynamo on rotation. We advocate the picture that the magnetic dynamo generates magnetic flux on the stellar surface proportional to rotation rate with a saturation limit set by the available kinetic energy, and we provide relations for average field strengths and nonthermal emission that are independent of the choice of the convective turnover time. We also find that Ca H&K emission saturates at average field strengths of ⟨B⟩≈800 G while Hα and X-ray emission grow further with stronger fields in the more rapidly rotating stars. This is in conflict with the coronal stripping scenario predicting that in the most rapidly rotating stars coronal plasma would be cooled to chromospheric temperatures. © ESO 2022., CARMENES is an instrument at the Centro Astronómico Hispano-Alemán (CAHA) at Calar Alto (Almería, Spain), operated jointly by the Junta de Andalucía and the Instituto de Astrofísica de Andalucía (CSIC). The authors wish to express their sincere thanks to all members of the Calar Alto staff for their expert support of the instrument and telescope operation. CARMENES was funded by the Max-Planck-Gesellschaft (MPG), the Consejo Superior de Investigaciones Científicas (CSIC), the Ministerio de Economía y Competitividad (MINECO) and the European Regional Development Fund (ERDF) through projects FICTS-2011-02, ICTS-2017-07-CAHA-4, and CAHA16-CE-3978, and the members of the CARMENES Consortium (Max-Planck-Institut für Astronomie, Instituto de Astrofísica de Andalucía, Landessternwarte Königstuhl, Institut de Ciències de l’Espai, Institut für Astrophysik Göttingen, Universidad Complutense de Madrid, Thüringer Landessternwarte Tautenburg, Instituto de Astrofísica de Canarias, Hamburger Sternwarte, Centro de Astrobiología and Centro Astronómico Hispano-Alemán), with additional contributions by the MINECO, the Deutsche Forschungsgemeinschaft through the Major Research Instrumentation Programme and Research Unit FOR2544 “Blue Planets around Red Stars”, the Klaus Tschira Stiftung, the states of Baden-Württemberg and Niedersachsen, by the Junta de Andalucía. We acknowledge financial support from the Agencia Estatal de Investigación of the Ministerio de Ciencia e Innovación and the ERDF “A way of making Europe” through project PID2019-109522GB-C5[1:4]/AEI/10.13039/501100011033 and the Centre of Excellence “Severo Ochoa” and “María de Maeztu” awards to the Instituto de Astrofísica de Canarias (CEX2019-000920-S), Instituto de Astrofísica de Andalucía (SEV-2017-0709), and Centro de Astrobiología (MDM-2017-0737), the Deutsche Forschungsgemeinschaft Heisenberg programme (KA4825/4-1), and the Generalitat de Catalunya/CERCA programme.

Cool main-sequence stars, such as the Sun, have magnetic fields which are generated by an internal dynamo mechanism. In the Sun, the dynamo mechanism produces a balance between the amounts of magnetic flux generated and lost over the Sun's 11-year activity cycle and it is visible in the Sun's different atmospheric layers using multi-wavelength observations. We used the same observational diagnostics, spanning several decades, to probe the emergence of magnetic flux on the two close by, active- and low-mass K dwarfs: 61 Cygni A and Epsilon Eridani. Our results show that 61 Cygni A follows the Solar dynamo with a regular cycle at all wavelengths, while Epsilon Eridani represents a more extreme level of the Solar dynamo, while also showing strong Solar-like characteristics. For the first time we show magnetic butterfly diagrams for stars other than the Sun. For the two K stars and the Sun, the rate at which the toroidal field is generated from surface poloidal field is similar to the rate at which toroidal flux is lost through flux emergence. This suggests that the surface field plays a crucial role in the dynamos of all three stars. Finally, for Epsilon Eridani, we show that the two chromospheric cycle periods, of ~3 and ~13 years, correspond to two superimposed magnetic cycles., Comment: 8 pages, 5 figures: Accepted by A&A

Observations of rapidly-rotating cool stars often show coronal slingshot prominences that remove mass and angular momentum when they are ejected. The derived masses of these prominences show a scatter of some two orders of magnitude. In order to investigate if this scatter could be intrinsic, we use a full magnetic cycle of solar magnetograms to model the coronal structure and prominence distribution in a young Sun, where we scale the field strength in the magnetograms with angular velocity according to $ B \propto \Omega^{-1.32} $. We reproduce both the observed prominence masses and their scatter. We show that both the field strength and the field geometry contribute to the prominence masses that can be supported and to the rate at which they are ejected. Predicted prominence masses follow the magnetic cycle, but with half the period, peaking both at cycle maximum and at cycle minimum. We show that mass loss rates in prominences are less than those predicted for the stellar wind. We also investigate the role of small-scale field that may be unresolved in typical stellar magnetograms. This provides only a small reduction in the predicted total prominence mass, principally by reducing the number of large magnetic loops that can support slingshot prominences. We conclude that the observed scatter in prominence masses can be explained by underlying magnetic cycles., Comment: 10 pages, 7 figures. Accepted for publication in MNRAS

Spectropolarimetric data allow for simultaneous monitoring of stellar chromospheric $\log{R^{\prime}_{\rm{HK}}}$ activity and the surface-averaged longitudinal magnetic field, $B_l$, giving the opportunity to probe the relationship between large-scale stellar magnetic fields and chromospheric manifestations of magnetism. We present $\log{R^{\prime}_{\rm{HK}}}$ and/or $B_l$ measurements for 954 mid-F to mid-M stars derived from spectropolarimetric observations contained within the PolarBase database. Our magnetically active sample complements previous stellar activity surveys that focus on inactive planet-search targets. We find a positive correlation between mean $\log{R^{\prime}_{\rm{HK}}}$ and mean $\log|B_l|$, but for G stars the relationship may undergo a change between $\log{R'_{\rm{HK}}}\sim-4.4$ and $-4.8$. The mean $\log{R^{\prime}_{\rm{HK}}}$ shows a similar change with respect to the $\log{R^{\prime}_{\rm{HK}}}$ variability amplitude for intermediately-active G stars. We also combine our results with archival chromospheric activity data and published observations of large-scale magnetic field geometries derived using Zeeman Doppler Imaging. The chromospheric activity data indicate a slight under-density of late-F to early-K stars with $-4.75\leq\log{R'_{\rm HK}}\leq-4.5$. This is not as prominent as the original Vaughan-Preston gap, and we do not detect similar under-populated regions in the distributions of the mean $|B_l|$, or the $B_l$ and $\log{R'_{\rm HK}}$ variability amplitudes. Chromospheric activity, activity variability and toroidal field strength decrease on the main sequence as rotation slows. For G stars, the disappearance of dominant toroidal fields occurs at a similar chromospheric activity level as the change in the relationships between chromospheric activity, activity variability and mean field strength., 20 pages, 15 figures, 3 tables

To push the radial velocity (RV) exoplanet detection threshold, it is crucial to find more reliable radial velocity extraction methods. The Least-Squares Deconvolution (LSD) technique has been used to infer the stellar magnetic flux from spectropolarimetric data for the past two decades. It relies on the assumption that stellar absorption lines are similar in shape. Although this assumption is simplistic, LSD provides a good model for intensity spectra and likewise an estimate for their Doppler shift. We present the Multi-Mask Least-Squares Deconvolution (MM-LSD) RV extraction pipeline which extracts the radial velocity from two-dimensional echelle-order spectra using LSD with multiple tailored masks after continuum normalisation and telluric absorption line correction. The flexibility of LSD allows to exclude spectral lines or pixels at will, providing a means to exclude variable lines or pixels affected by instrumental problems. The MM-LSD pipeline was tested on HARPS-N data for the Sun and selected well-observed stars with 5.7 < Vmag < 12.6. For FGK-type stars with median signal-to-noise above 100, the pipeline delivered RV time series with on average 12 per cent lower scatter as compared to the HARPS-N RV extraction pipeline based on the Cross-Correlation Function technique. The MM-LSD pipeline may be used as a standalone RV code, or modified and extended to extract a proxy for the magnetic field strength., Accepted for publication in MNRAS. Code available on github: https://github.com/florian-lienhard/MM-LSD. 16 pages, 15 figures

The long-term evolution of the internal, strong magnetic fields of neutron stars needs a specific numerical modelling. The diversity of the observed phenomenology of neutron stars indicates that their magnetic topology is rather complex and three-dimensional simulations are required, for example, to explain the observed bursting mechanisms and the creation of surface hotspots. We present MATINS, a new three dimensions numerical code for magneto-thermal evolution in neutron stars, based on a finite-volume scheme that employs the cubed-sphere system of coordinates. In this first work, we focus on the crustal magnetic evolution, with the inclusion of realistic calculations for the neutron star structure, composition and electrical conductivity assuming a simple temperature evolution profile. MATINS follows the evolution of strong fields (1e14-1e15 Gauss) with complex non-axisymmetric topologies and dominant Hall-drift terms, and it is suitable for handling sharp current sheets. After introducing the technical description of our approach and some tests, we present long-term simulations of the non-linear field evolution in realistic neutron star crusts. The results show how the non-axisymmetric Hall cascade redistributes the energy over different spatial scales. Following the exploration of different initial topologies, we conclude that during a few tens of kyr, an equipartition of energy between the poloidal and toroidal components happens at small-scales. However, the magnetic field keeps a strong memory of the initial large-scales, which are much harder to be restructured or created. This indicates that large-scale configuration attained during the neutron star formation is crucial to determine the field topology at any evolution stage., 21 pages, 17 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The primary component of the multiple star phi Dra is one of the brightest magnetic chemically peculiar stars in the northern sky. Here we report results of a comprehensive study of the rotational photometric variability, binarity, magnetic field geometry, and surface chemical spot structure for this star. We derived a precise photometric rotational period of 1.71650213(21) d based on one year of TESS nearly continuous space observations and discovered modulation of the stellar light curve with the phase of the 127.9-d binary orbit due to the light time travel effect. We revised parameters of the binary orbit and detected spectroscopic contribution of the secondary. A tomographic mapping technique was applied to the average intensity and circular polarisation profiles derived from Narval high-resolution spectropolarimetric observations. This analysis yielded a detailed map of the global magnetic field topology together with the surface distributions of Si, Cr, and Fe abundances. Magnetic mapping demonstrates that the surface field structure of phi Dra is dominated by a distorted dipolar component with a peak field strength of 1.4 kG and a large asymmetry between the poles. Chemical maps show an enhancement of Cr, Fe and, to a lesser extent, Si in a series of spots encircling intersections of the magnetic and rotational equators. These chemical spot geometries do not directly correlate with either the local field strength or the field inclination., Comment: 13 pages, 13 figures; accepted for publication in MNRAS

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star V410 Tau based on data collected mostly with SPIRou, the near-infrared (NIR) spectropolarimeter recently installed at the Canada-France-Hawaii Telescope, as part of the SPIRou Legacy Survey large programme, and with TESS between October and December 2019. Using Zeeman-Doppler Imaging (ZDI), we obtained the first maps of photospheric brightness and large-scale magnetic field at the surface of this young star derived from NIR spectropolarimetric data. For the first time, ZDI is also simultaneously applied to high-resolution spectropolarimetric data and very-high-precision photometry. V410 Tau hosts both dark and bright surface features and magnetic regions similar to those previously imaged with ZDI from optical data, except for the absence of a prominent dark polar spot. The brightness distribution is significantly less contrasted than its optical equivalent, as expected from the difference in wavelength. The large-scale magnetic field (~410 G), found to be mainly poloidal, features a dipole of ~390 G, again compatible with previous studies at optical wavelengths. NIR data yield a surface differential rotation slightly weaker than that estimated in the optical at previous epochs. Finally, we measured the radial velocity of the star and filtered out the stellar activity jitter using both ZDI and Gaussian Process Regression down to a precision of ~0.15 and 0.08 $\mathrm{km\,s^{-1}}$ RMS, respectively, confirming the previously published upper limit on the mass of a potential close-in massive planet around V410 Tau., 19 pages, 27 figures, accepted for publication in MNRAS

Neutron star models with maximum mass close to $2 \ M_{\odot}$ reach high central densities, which may activate nucleonic and hyperon direct Urca neutrino emission. To alleviate the tension between fast theoretical cooling rates and thermal luminosity observations of moderately magnetized, isolated thermally-emitting stars (with $L_{\gamma} \gtrsim 10^{31}$ erg s$^{-1}$ at $t \gtrsim 10^{5.3}$ yr), some internal heating source is required. The power supplied by the internal heater is estimated for both a phenomenological source in the inner crust and Joule heating due to magnetic field decay, assuming different superfluidity models and compositions of the outer stellar envelope. It is found that a thermal power of $W(t) \approx 10^{34}$ erg s$^{-1}$ allows neutron star models to match observations of moderately magnetized, isolated stars with ages $t \gtrsim 10^{5.3}$ yr. The requisite $W(t)$ can be supplied by Joule heating due to crust-confined initial magnetic configurations with (i) mixed poloidal-toroidal fields, with surface strength $B_{\textrm{dip}} = 10^{13}$ G at the pole of the dipolar poloidal component and $\sim 90$ per cent of the magnetic energy stored in the toroidal component; and (ii) poloidal-only configurations with $B_{\textrm{dip}} = 10^{14}$ G., Comment: 15 pages, 9 figures. Accepted for publication in MNRAS

Context. The magnetic activity of the Sun changes with the solar cycle. Similar cycles are found in other stars as well, but their details are not known to a similar degree. Characterising stellar magnetic cycles is important for the understanding of the stellar and solar dynamos that are driving the magnetic activity. Aims: We present spectropolarimetric observations of five young, solar-type stars and compare them to previous observations, with the aim to identify and characterise stellar equivalents of the solar cycle. Methods: We use Zeeman-Doppler imaging (ZDI) to map the surface magnetic field and brightness of our targets. The magnetic field is decomposed into spherical harmonic expansions, from which we report the strengths of the axisymmetric versus non-axisymmetric and poloidal versus toroidal components, and we compare them to the Rossby numbers of the stars. Results: We present five new ZDI maps of young, solar-type stars from December 2017. Of special interest is the case of V1358 Ori, which had gone through a polarity reversal between our observations and earlier ones. A less evident polarity reversal might also have occurred in HD 35296. There is a preference for a more axisymmetric field, and possibly a more toroidal field, for the more active stars with lower Rossby number, but a larger sample should be studied to draw any strong conclusions from this. For most of the individual stars, the amounts of toroidal and poloidal field have stayed on levels similar to those in earlier observations. Conclusions: We find evidence for a magnetic polarity reversal having occurred in V1358 Ori. An interesting target for future observations is χ1 Ori, which may have a short magnetic cycle of a few years. The correlation between the brightness maps and the magnetic field is mostly poor, which could indicate the presence of small-scale magnetic features of different polarities that cancel one another out and are thus not resolved in our maps. The data for the magnetic field and brightness maps of the stars are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/659/A71 Based on observations made with the HARPSpol instrument on the ESO 3.6 m telescope in La Silla (Chile), under programme IDs 091.D-0836 and 0100.D-0176.

Recent large-scale spectropolarimetric surveys have established that a small but significant percentage of massive stars host stable, surface dipolar magnetic fields with strengths on the order of kG. These fields channel the dense, radiatively driven stellar wind into circumstellar magnetospheres, whose density and velocity structure can be probed using ultraviolet (UV) spectroscopy of wind-sensitive resonance lines. Coupled with appropriate magnetosphere models, UV spectroscopy provides a valuable way to investigate the wind-field interaction, and can yield quantitative estimates of the wind parameters of magnetic massive stars. We report a systematic investigation of the formation of UV resonance lines in slowly rotating magnetic massive stars with dynamical magnetospheres. We pair the Analytic Dynamical Magnetosphere (ADM) formalism with a simplified radiative transfer technique to produce synthetic UV line profiles. Using a grid of models, we examine the effect of magnetosphere size, the line strength parameter, and the cooling parameter on the structure and modulation of the line profile. We find that magnetic massive stars uniquely exhibit redshifted absorption at most viewing angles and magnetosphere sizes, and that significant changes to the shape and variation of the line profile with varying line strengths can be explained by examining the individual wind components described in the ADM formalism. Finally, we show that the cooling parameter has a negligible effect on the line profiles., 16 pages, 15 figures, accepted to MNRAS

$\tau$ Sco, a well-studied magnetic B-type star in the Upper Sco association, has a number of surprising characteristics. It rotates very slowly and shows nitrogen excess. Its surface magnetic field is much more complex than a purely dipolar configuration which is unusual for a magnetic massive star. We employ the CMFGEN radiative transfer code to determine the fundamental parameters and surface CNO and helium abundances. Then, we employ MESA and GENEC stellar evolution models accounting for the effects of surface magnetic fields. To reconcile $\tau$ Sco's properties with single-star models, an increase is necessary in the efficiency of rotational mixing by a factor of 3 to 10 and in the efficiency of magnetic braking by a factor of 10. The spin down could be explained by assuming a magnetic field decay scenario. However, the simultaneous chemical enrichment challenges the single-star scenario. Previous works indeed suggested a stellar merger origin for $\tau$ Sco. However, the merger scenario also faces similar challenges as our magnetic single-star models to explain $\tau$ Sco's simultaneous slow rotation and nitrogen excess. In conclusion, the single-star channel seems less likely and versatile to explain these discrepancies, while the merger scenario and other potential binary-evolution channels still require further assessment as to whether they may self-consistently explain the observables of $\tau$ Sco., Comment: Accepted for publication in MNRAS. A full reproduction package is shared on zenodo in accordance with the Research Data Management plan of the Anton Pannekoek Institute for Astronomy at the University of Amsterdam: https://doi.org/10.5281/zenodo.4633408

There have recently been detections of radio emission from low-mass stars, some of which are indicative of star-planet interactions. Motivated by these exciting new results, in this paper we present Alfv\'en wave-driven stellar wind models of the two active planet-hosting M dwarfs Prox Cen and AU Mic. Our models incorporate large-scale photospheric magnetic field maps reconstructed using the Zeeman-Doppler Imaging method. We obtain a mass-loss rate of $0.25~\dot{M}_{\odot}$ for the wind of Prox Cen. For the young dwarf AU Mic, we explore two cases: a low and high mass-loss rate. Depending on the properties of the Alfv\'en waves which heat the corona in our wind models, we obtain mass-loss rates of $27$ and $590~\dot{M}_{\odot}$ for AU Mic. We use our stellar wind models to assess the generation of electron cyclotron maser instability emission in both systems, through a mechanism analogous to the sub-Alfv\'enic Jupiter-Io interaction. For Prox Cen we do not find any feasible scenario where the planet can induce radio emission in the star's corona, as the planet orbits too far from the star in the super-Alfv\'enic regime. However, in the case that AU Mic has a stellar wind mass-loss rate of $27~\dot{M}_{\odot}$, we find that both planets b and c in the system can induce radio emission from $\sim10$ MHz to 3 GHz in the corona of the host star for the majority of their orbits, with peak flux densities of $\sim10$ mJy. Detection of such radio emission would allow us to place an upper limit on the mass-loss rate of the star., Comment: 8 pages, 5 figures. Accepted for publication in MNRAS

Energetic particles, such as stellar cosmic rays, produced at a heightened rate by active stars (like the young Sun) may have been important for the origin of life on Earth and other exoplanets. Here we compare, as a function of stellar rotation rate ($\Omega$), contributions from two distinct populations of energetic particles: stellar cosmic rays accelerated by impulsive flare events and Galactic cosmic rays. We use a 1.5D stellar wind model combined with a spatially 1D cosmic ray transport model. We formulate the evolution of the stellar cosmic ray spectrum as a function of stellar rotation. The maximum stellar cosmic ray energy increases with increasing rotation i.e., towards more active/younger stars. We find that stellar cosmic rays dominate over Galactic cosmic rays in the habitable zone at the pion threshold energy for all stellar ages considered ($t_*=0.6-2.9\,$Gyr). However, even at the youngest age, $t_*=0.6\,$Gyr, we estimate that $\gtrsim\,80$MeV stellar cosmic ray fluxes may still be transient in time. At $\sim1\,$Gyr when life is thought to have emerged on Earth, we demonstrate that stellar cosmic rays dominate over Galactic cosmic rays up to $\sim$4$\,$GeV energies during flare events. Our results for $t_*=0.6\,$Gyr ($\Omega = 4\Omega_\odot$) indicate that $\lesssim$GeV stellar cosmic rays are advected from the star to 1$\,$au and are impacted by adiabatic losses in this region. The properties of the inner solar wind, currently being investigated by the Parker Solar Probe and Solar Orbiter, are thus important for accurate calculations of stellar cosmic rays around young Sun-like stars., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

The helicity of a magnetic field is a fundamental property that is conserved in ideal MHD. It can be explored in the stellar context by mapping large-scale magnetic fields across stellar surfaces using Zeeman-Doppler imaging. A recent study of 51 stars in the mass range 0.1-1.34 M$_\odot$ showed that the photospheric magnetic helicity density follows a single power law when plotted against the toroidal field energy, but splits into two branches when plotted against the poloidal field energy. These two branches divide stars above and below $\sim$ 0.5 M$_\odot$. We present here a novel method of visualising the helicity density in terms of the linkage of the toroidal and poloidal fields that are mapped across the stellar surface. This approach allows us to classify the field linkages that provide the helicity density for stars of different masses and rotation rates. We find that stars on the lower-mass branch tend to have toroidal fields that are non-axisymmetric and so link through regions of positive and negative poloidal field. A lower-mass star may have the same helicity density as a higher-mass star, despite having a stronger poloidal field. Lower-mass stars are therefore less efficient at generating large-scale helicity., Comment: 8 pages, 7 figures, 1 table

The discovery of the moderate differential rotation between the core and the envelope of evolved solar-like stars could be the signature of a strong magnetic field trapped inside the radiative interior. The population of intermediate-mass red giants presenting a surprisingly low-amplitude of their mixed modes could also arise from the effect of an internal magnetic field. Indeed, stars more massive than about 1.1Ms are known to develop a convective core during their main sequence, which could relax into a strong fossil magnetic field trapped inside the core of the star for the rest of its evolution. The observations of mixed modes can constitute an excellent probe of the deepest layers of evolved solar-like stars. The magnetic perturbation on mixed modes may thus be visible in asteroseismic data. To unravel which constraints can be obtained from observations, we theoretically investigate the effects of a plausible mixed axisymmetric magnetic field with various amplitudes on the mixed-mode frequencies of evolved solar-like stars. The first-order frequency perturbations are computed for dipolar and quadrupolar mixed modes. These computations are carried out for a range of stellar ages, masses, and metallicities. We show that typical fossil-field strengths of 0.1-1 MG, consistent with the presence of a dynamo in the convective core during the main sequence, provoke significant asymmetries on mixed-mode frequency multiplets during the red-giant branch. We show that these signatures may be detectable in asteroseismic data for field amplitudes small enough for the amplitude of the modes not to be affected by the conversion of gravity into Alfven waves inside the magnetised interior. Finally, we infer an upper limit for the strength of the field, and the associated lower limit for the timescale of its action, to redistribute angular momentum in stellar interiors., 30 pages, 19 figures, Submitted to A&A, Revised version

Studying cool star magnetic activity gives an important insight into the stellar dynamo and its relationship with stellar properties, as well as allowing us to place the Sun's magnetism in the context of other stars. Only 61 Cyg A (K5V) and $\tau$ Boo (F8V) are currently known to have magnetic cycles like the Sun's, where the large-scale magnetic field polarity reverses in phase with the star's chromospheric activity cycles. ${\tau}$ Boo has a rapid $\sim$240 d magnetic cycle, and it is not yet clear whether this is related to the star's thin convection zone or if the dynamo is accelerated by interactions between ${\tau}$ Boo and its hot Jupiter. To shed light on this, we studied the magnetic activity of HD75332 (F7V) which has similar physical properties to ${\tau}$ Boo and does not appear to host a hot Jupiter. We characterized its long term chromospheric activity variability over 53 yrs and used Zeeman Doppler Imaging to reconstruct the large-scale surface magnetic field for 12 epochs between 2007 and 2019. Although we observe only one reversal of the large-scale magnetic dipole, our results suggest that HD75332 has a rapid $\sim$1.06 yr solar-like magnetic cycle where the magnetic field evolves in phase with its chromospheric activity. If a solar-like cycle is present, reversals of the large-scale radial field polarity are expected to occur at around activity cycle maxima. This would be similar to the rapid magnetic cycle observed for ${\tau}$ Boo, suggesting that rapid magnetic cycles may be intrinsic to late-F stars and related to their shallow convection zones., Comment: 23 pages, 14 figures, accepted for publication in MNRAS

Whilst ``slingshot'' prominences have been observed on M-dwarfs, most if not all theoretical studies have focused on solar-like stars. We present an investigation into stellar prominences around rapidly rotating young M-dwarfs. We have extrapolated the magnetic field in the corona from Zeeman-Doppler maps and determined the sites of mechanical stability where prominences may form. We analyse the prominence mass that could be supported and the latitude range over which this material is distributed. We find that for these maps, much of this prominence mass may be invisible to observation - typically, 13 pages, 10 figures

We thank the anonymous referee for a careful reading of the manuscript. DLH acknowledges financial support from the Science and Technology Facilities Council (STFC) via grant ST/M000877/1. MSC acknowledges the support by FCT/MCTES through the research grants UIDB/04434/2020, UIDP/04434/2020, and PTDC/FIS-AST/30389/2017, and by FEDER - Fundo Europeu de Desenvolvimento Regional through COMPETE2020 - Programa Operacional Competitividade e Internacionalizacao (grant: POCI01-0145-FEDER-030389). MSC is supported by national funds through FCT in the form of a work contract. VA was supported by a research grant (00028173) fromVILLUMFONDEN. Funding for the Stellar Astrophysics Centre is provided by The Danish National Research Foundation (Grant agreement no.: DNRF106). The research leading to these results has received funding from the Research Foundation Flanders (FWO) by means of a senior postdoctoral fellowship to DMB with grant agreement No. 1286521N. This project has been supported by the Lendulet Program of the Hungarian Academy of Sciences, project No. LP2018-7/2020 and by the EU's MW-Gaia COST Action (CA18104). ASG acknowledges financial support from the Max Planck Society under grant 'Preparations for PLATO science' and from ALMA-CONICYT under grant #31170029. ZsB acknowledges the support by the Janos Bolyai Research Scholarship of the Hungarian Academy of Sciences. LFM acknowledges the financial support from the UNAM under grant PAPIIT IN100918. CCL acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC). DM acknowledges the support of the National Astronomical Research Institute of Thailand (NARIT). JPG acknowledges funding support from Spanish public funds for research from project PID2019-107061GB-C63 from the 'Programas Estatales de Generacion de Conocimiento y Fortalecimiento Cientifico y Tecnologico del Sistema de I+D+i y de I+D+i Orientada a los Retos de la Sociedad', and from the State Agency for Research through the Center of Excellence Severo Ochoa' award to the Instituto de Astrof ' isica de Andaluc ' ia (SEV2017-0709), all from the Spanish Ministry of Science, Innovation and Universities (MCIU). TRY acknowledges support from the NSF REU program, grant number PHY-1359195. The material is based upon work supported by NASA under award number 80GSFC21M0002. AD was supported by the UNKP-20-5 New National Excellence Program of the Ministry of Human Capacities and the J ' anos Bolyai Research Scholarship of the Hungarian Academy of Sciences. AD would like to thank the City of Szombathely for support under Agreement No. 67.177-21/2016. AGH acknowledges support from 'FEDER/Junta de Andalucia-Consejeria de Economia y Conocimiento' under project E-FQM-041-UGR18 by Universidad de Granada and from Spanish public funds (including FEDER funds) for research under project ESP2017-87676-C5-2-R. ARB acknowledges funding support from Spanish public funds for research through the research grant PRE2018-084322, through projects ESP2017-87676-C5-5-R and PID2019-107061GB-C63, and from the State Agency for Research through the `Center of Excellence Severo Ochoa' award to the Instituto de Astrof ' isica de Andaluc ' ia (SEV-2017-0709), all from the Spanish Ministry of Science, Innovation and Universities (MCIU). We are grateful to Sowgata Chowdhury, supported by Polish NCN grant 2015/18/A/ST9/00578, for conducting some spectroscopic observations that are reported in this work. This paper includes data collected by the TESS mission, which are publicly available from the Mikulski Archive for Space Telescopes (MAST). Funding for the TESS mission is provided by NASA's Science Mission directorate. Funding for the TESS Asteroseismic Science Operations Centre is provided by the Danish National Research Foundation (Grant agreement no.: DNRF106), ESA PRODEX (PEA 4000119301) and Stellar Astrophysics Centre (SAC) at Aarhus University. We thank the TESS team and staff and TASC/TASOC for their support of this work. This research has made use of the SIMBAD data base, operated at CDS, Strasbourg, France. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement., We present the results of a systematic search for new rapidly oscillating Ap (roAp) stars using the 2-min cadence data collected by the Transiting Exoplanet Survey Satellite (TESS) during its Cycle 1 observations. We identify 12 new roAp stars. Amongst these stars we discover the roAp star with the longest pulsation period, another with the shortest rotation period, and six with multiperiodic variability. In addition to these new roAp stars, we present an analysis of 44 known roAp stars observed by TESS during Cycle 1, providing the first high-precision and homogeneous sample of a significant fraction of the known roAp stars. The TESS observations have shown that almost 60 per cent (33) of our sample of stars are multiperiodic, providing excellent cases to test models of roAp pulsations, and from which the most rewarding asteroseismic results can be gleaned. We report four cases of the occurrence of rotationally split frequency multiplets that imply different mode geometries for the same degree modes in the same star. This provides a conundrum in applying the oblique pulsator model to the roAp stars. Finally, we report the discovery of non-linear mode interactions in α Cir (TIC 402546736, HD128898) around the harmonic of the principal mode – this is only the second case of such a phenomenon., UK Research & Innovation (UKRI), Science & Technology Facilities Council (STFC), Science and Technology Development Fund (STDF) ST/M000877/1, Portuguese Foundation for Science and Technology, European Commission UIDB/04434/2020 UIDP/04434/2020 PTDC/FIS-AST/30389/2017, FEDER - Fundo Europeu de Desenvolvimento Regional through COMPETE2020 - Programa Operacional Competitividade e Internacionalizacao POCI-01-0145-FEDER-030389, European Commission, VILLUMFONDEN 00028173, Danmarks Grundforskningsfond DNRF106 FWO 1286521N, Lendulet Program of the Hungarian Academy of Sciences LP2018-7/2020, EU's MW-Gaia COST Action CA18104, Max Planck Society, Foundation CELLEX, ALMA-CONICYT 31170029, Hungarian Academy of Sciences, UNAM under grant PAPIIT IN100918, Natural Sciences and Engineering Research Council of Canada (NSERC), National Astronomical Research Institute of Thailand (NARIT), Spanish public funds the Spanish Ministry of Science, Innovation and Universities (MCIU) PID2019-107061GB-C63 PRE2018-084322 ESP2017-87676-C5-5-R, 'Programas Estatales de Generacion de Conocimiento y Fortalecimiento Cientifico y Tecnologico del Sistema de I+D+i y de I+D+i Orientada a los Retos de la Sociedad' from the Spanish Ministry of Science, Innovation and Universities (MCIU), State Agency for Research through the 'Center of Excellence Severo Ochoa' award from the Spanish Ministry of Science, Innovation and Universities (MCIU) SEV-2017-0709, National Science Foundation (NSF), NSF - Office of the Director (OD) PHY-1359195, National Aeronautics & Space Administration (NASA) 80GSFC21M0002, New National Excellence Program of the Ministry of Human Capacities UNKP-20-5, City of Szombathely 67.177-21/2016, 'FEDER/Junta de Andalucia-Consejeria de Economia y Conocimiento' by Universidad de Granada E-FQM-041-UGR18, Spanish public funds (FEDER) ESP2017-87676-C5-2-R, Polish NCN 2015/18/A/ST9/00578, NASA's Science Mission directorate, European Space Agency PEA 4000119301, Stellar Astrophysics Centre (SAC) at Aarhus University

Context. Recent observational progress has challenged the dust grain-alignment theories used to explain the polarized dust emission routinely observed in star-forming cores. Aims. In an effort to improve our understanding of the dust grain alignment mechanism(s), we have gathered a dozen ALMA maps of (sub)millimeter-wavelength polarized dust emission from Class 0 protostars and carried out a comprehensive statistical analysis of dust polarization quantities. Methods. We analyze the statistical properties of the polarization fraction Pfrac and the dispersion of polarization position angles S. More specifically, we investigate the relationship between S and Pfrac as well as the evolution of the product S × Pfrac as a function of the column density of the gas in the protostellar envelopes. We compare the observed trends with those found in polarization observations of dust in the interstellar medium and in synthetic observations of non-ideal magneto-hydrodynamic (MHD) simulations of protostellar cores. Results. We find a significant S Pfrac-0.79 correlation in the polarized dust emission from protostellar envelopes seen with ALMA; the power-law index significantly differs from the one observed by Planck in star-forming clouds. The product S × Pfrac, which is sensitive to the dust grain alignment efficiency, is approximately constant across three orders of magnitude in envelope column density (from NH2 = 1022 cm-2 to NH2 = 1025 cm-2), with a mean value of 0.36-0.17+0.10. This suggests that the grain alignment mechanism producing the bulk of the polarized dust emission in star-forming cores may not systematically depend on the local conditions such as the local gas density. However, in the lowest-luminosity sources in our sample, we find a hint of less efficient dust grain alignment with increasing column density. Our observations and their comparison with synthetic observations of MHD models suggest that the total intensity versus the polarized dust are distributed at different intrinsic spatial scales, which can affect the statistics from the ALMA observations, for example, by producing artificially high Pfrac. Finally, synthetic observations of MHD models implementing radiative alignment torques (RATs) show that the statistical estimator S × Pfrac is sensitive to the strength of the radiation field in the core. Moreover, we find that the simulations with a uniform perfect alignment (PA) of dust grains yield, on average, much higher S × Pfrac values than those implementing RATs; the ALMA values lie among those predicted by PA, and they are significantly higher than the ones obtained with RATs, especially at large column densities. Conclusions. Ultimately, our results suggest that dust alignment mechanism(s) are efficient at producing dust polarized emission in the various local conditions typical of Class 0 protostars. The grain alignment efficiency found in these objects seems to be higher than the efficiency produced by the standard RAT alignment of paramagnetic grains. Further studies will be needed to understand how more efficient grain alignment via, for example, different irradiation conditions, dust grain characteristics, or additional grain alignment mechanisms can reproduce the observations., V.J.M.L.G. acknowledges the support of the ESO Studentship Program. A.J.M. acknowledges support from the Joint ALMA Observatory Visitor Program, and ESO Visitor program. J.M.G. is supported by the Spanish grant AYA2017-84390-C2-R (AEI/FEDER, UE). C.L.H.H. acknowledges the support of both the NAOJ Fellowship as well as JSPS KAKENHI grants 18K13586 and 20K14527. This work has benefited from the support of the European Research Council under the Horizon 2020 Framework Programme (Starting Grant MagneticYSOs with grant agreement no. 679937)

Strong magnetic fields in chemically peculiar A-type (Ap) stars typically suppress low-overtone pressure modes (p modes) but allow high-overtone p modes to be driven. KIC 11296437 is the first star to show both. We obtained and analysed a Subaru spectrum, from which we show that KIC 11296437 has abundances similar to other magnetic Ap stars, and we estimate a mean magnetic field modulus of $2.8\pm0.5$ kG. The same spectrum rules out a double-lined spectroscopic binary, and we use other techniques to rule out binarity over a wide parameter space, so the two pulsation types originate in one $\delta$ Sct--roAp hybrid pulsator. We construct stellar models depleted in helium and demonstrate that helium settling is second to magnetic damping in suppressing low-overtone p modes in Ap stars. We compute the magnetic damping effect for selected p and g modes, and find that modes with frequencies similar to the fundamental mode are driven for polar field strengths $\lesssim4$ kG, while other low-overtone p modes are driven for polar field strengths up to $\sim$1.5 kG. We find that the high-order g modes commonly observed in $\gamma$ Dor stars are heavily damped by polar fields stronger than 1--4 kG, with the damping being stronger for higher radial orders. We therefore explain the observation that no magnetic Ap stars have been observed as $\gamma$ Dor stars. We use our helium-depleted models to calculate the $\delta$ Sct instability strip for metallic-lined A (Am) stars, and find that driving from a Rosseland mean opacity bump at $\sim$$5\times10^4$ K caused by the discontinuous H-ionization edge in bound-free opacity explains the observation of $\delta$ Sct pulsations in Am stars., Comment: Accepted for publication in MNRAS

For moderate and slow rotation, magnetic activity of solar-like stars is observed to strongly depend on rotation. These observations do not yet have a solid explanation in terms of dynamo theory. We aim to find such an explanation by numerically investigated the rotational dependency of dynamo drivers in solar-like stars. We ran semi-global convection simulations of stars with rotation rates from 0 to 30 times the solar value, corresponding to Coriolis numbers, Co, of 0 to 110. We measured the turbulent transport coefficients describing the magnetic field evolution with the help of the test-field method, and compared with the dynamo effect arising from the differential rotation. The trace of the $\alpha$ tensor increases for moderate rotation rates with Co$^{0.5}$ and levels off for rapid rotation. This behavior agrees with the kinetic $\alpha$, if one considers the decrease of the convective scale with increasing rotation. The $\alpha$ tensor becomes highly anisotropic for Co$\gtrsim 1$. Furthermore, $\alpha_{rr}$ dominates for moderate rotation (1, Comment: 16 pages with 13 figures, plus 3 pages with 4 figures in the appendix, accepted for publication in A&A

We report time-resolved, high-resolution optical spectropolarimetric observations of the young double-lined spectroscopic binary V1878 Ori. Our observations were collected with the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope through the BinaMIcS large programme. V1878 Ori A and B are partially convective intermediate mass weak-line T Tauri stars on an eccentric and asynchronous orbit. We also acquired X-ray observations at periastron and outside periastron. Using the least-squares deconvolution technique (LSD) to combine information from many spectral lines, we clearly detected circular polarization signals in both components throughout the orbit. We refined the orbital solution for the system and obtained disentangled spectra for the primary and secondary components. The disentangled spectra were then employed to determine atmospheric parameters of the two components using spectrum synthesis. Applying our Zeeman Doppler imaging code to composite Stokes $IV$ LSD profiles, we reconstructed brightness maps and the global magnetic field topologies of the two components. We find that V1878 Ori A and B have strikingly different global magnetic field topologies and mean field strengths. The global magnetic field of the primary is predominantly poloidal and non-axisymmetric (with a mean field strength of 180 G). while the secondary has a mostly toroidal and axisymmetric global field (mean strength of 310 G). These findings confirm that stars with very similar parameters can exhibit radically different global magnetic field characteristics. The analysis of the X-ray data shows no sign of enhanced activity at periastron, suggesting the lack of strong magnetospheric interaction at this epoch., 13 pages, 9 figures, accepted for publication in MNRAS

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}$., Accepted for publication by A&A

In this paper, we present results from the spectropolarimetric follow-up of photometrically selected candidate magnetic B stars from the MOBSTER project. Out of four observed targets, one (HD 38170) is found to host a detectable surface magnetic field, with a maximum longitudinal field measurement of 105$\pm$14 G. This star is chemically peculiar and classified as an $\alpha^2$ CVn variable. Its detection validates the use of TESS to perform a photometric selection of magnetic candidates. Furthermore, upper limits on the strength of a putative dipolar magnetic field are derived for the remaining three stars, and we report the discovery of a previously unknown spectroscopic binary system, HD 25709. Finally, we use our non-detections as case studies to further inform the criteria to be used for the selection of a larger sample of stars to be followed up using high-resolution spectropolarimetry., Comment: 9 pages, 6 figures, 4 tables, accepted by MNRAS (2021 March 23)

Global magnetic fields of active solar-like stars are nowadays routinely detected with spectropolarimetric measurements and are mapped with Zeeman-Doppler imaging (ZDI). However, due to the cancellation of opposite field polarities, polarimetry captures only a tiny fraction of the magnetic flux and cannot assess the overall stellar surface magnetic field if it is dominated by a small-scale component. Analysis of Zeeman broadening in high-resolution intensity spectra can reveal these hidden complex magnetic fields. Historically, there were very few attempts to obtain such measurements for G dwarf stars due to the difficulty of disentangling Zeeman effect from other broadening mechanisms affecting spectral lines. Here we developed a new magnetic field diagnostic method based on relative Zeeman intensification of optical atomic lines with different magnetic sensitivity. Using this technique we obtained 78 field strength measurements for 15 Sun-like stars, including some of the best-studied young solar twins. We find that the average magnetic field strength $Bf$ drops from 1.3-2.0 kG in stars younger than about 120 Myr to 0.2-0.8 kG in older stars. The mean field strength shows a clear correlation with the Rossby number and with the coronal and chromospheric emission indicators. Our results suggest that magnetic regions have roughly the same local field strength $B\approx3.2$ kG in all stars, with the filling factor $f$ of these regions systematically increasing with stellar activity. Comparing our results with the spectropolarimetric analyses of global magnetic fields in the same stars, we find that ZDI recovers about 1% of the total magnetic field energy in the most active stars. This figure drops to just 0.01% for the least active targets., Comment: 25 pages, 26 figures; accepted for publication in A&A

We present the evolution of the X-ray emission properties of the magnetar 1E 1547.0-5408 since February 2004 over a time period covering three outbursts. We analyzed new and archival observations taken with the Swift, NuSTAR, Chandra, and XMM-Newton X-ray satellites. The source has been observed at a relatively steady soft X-ray flux of ≈10-11 erg cm-2 s-1 (0.3-10 keV) over the last 9 years, which is about an order of magnitude fainter than the flux at the peak of the last outburst in 2009, but a factor of ∼30 larger than the level in 2006. The broad-band spectrum extracted from two recent NuSTAR observations in April 2016 and February 2019 showed a faint hard X-ray emission up to ∼70 keV. Its spectrum is adequately described by a flat power law component, and its flux is ∼7 × 10-12 erg cm-2 s-1 (10-70 keV), that is a factor of ∼20 smaller than at the peak of the 2009 outburst. The hard X-ray spectral shape has flattened significantly in time, which is at variance with the overall cooling trend of the soft X-ray component. The pulse profile extracted from these NuSTAR pointings displays variability in shape and amplitude with energy (up to ≈25 keV). Our analysis shows that the flux of 1E 1547.0-5408 is not yet decaying to the 2006 level and that the source has been lingering in a stable, high-intensity state for several years. This might suggest that magnetars can hop among distinct persistent states that are probably connected to outburst episodes and that their persistent thermal emission can be almost entirely powered by the dissipation of currents in the corona., FCZ acknowledges Matthew Baring for helpful discussions. FCZ, AB, NR and DV acknowledge support from grants SGR2017-1383 and PGC2018-095512-B-I00 and the support of the PHAROS COST Action (CA16214). FCZ is also supported by a Juan de la Cierva fellowship. NR, AB and DV are also supported by the ERC Consolidator Grant “MAGNESIA” (nr. 817661). DV acknowledges support from grants AYA2016-80289-P and AYA2017-82089-ERC. TE is supported by JSPS KAKENHI grant numbers 16H02198, 18H01246 and 18H04584. JAP acknowledges support from the Spanish Agencia Estatal de Investigación (grant PGC2018-095984-B-I00) and the Generalitat Valenciana (grant PROMETEO/2019/071). The scientific results reported in this article are based on data obtained with Swift, the Chandra X-ray Observatory, XMM–Newton and the NuSTAR mission. XMM–Newton is an ESA science mission with instruments and contributions directly funded by ESA Member States and the National Aeronautics and Space Administration (NASA). The NuSTAR mission is a project led by the California Insitute of Technology, managed by the Jet Propulsion Laboratory, and funded by NASA. We made use of data supplied by the UK Swift Science Data Centre at the University of Leicester and of the XRT Data Analysis Software (XRTDAS) developed under the responsibility of the ASI Science Data Center (ASDC), Italy. We also made use of the NuSTAR Data Analysis Software (NUSTARDAS jointly developed by the ASI Science Data Center (ASDC, Italy) and the California Institute of Technology (USA), and of softwares and tools provided by the High Energy Astrophysics Science Archive Research Center (HEASARC), which is a service of the Astrophysics Science Division at NASA/GSFC and the High Energy Astrophysics Division of the Smithsonian Astrophysical Observatory.

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., 11 pages, 7 figures

Magnetic fields rule many physical processes in and around stars throughout their lifetime. All cool stars possess a magnetic field, likely generated by dynamo processes. In order to properly understand the evolution of cool stars, we need to understand their magnetism. Stellar magnetic fields can be directly observed through the imprint of the Zeeman effect in intensity and polarized spectra. In intensity spectra (Stokes I), spectral lines are broadened or split into several components by the magnetic field. Modelling this effect in high-resolution spectra allows us to determine the average unsigned magnetic field strength over the stellar surface. The magnetic field also induces circular (Stokes V) and linear polarization (Stokes QU) in spectral lines, according to its orientation. These polarization signals can be used to map the large-scale magnetic field at the surface of the star using tomographic techniques such as Zeeman Doppler imaging (ZDI). In this thesis, we investigated pre-main-sequence T Tauri stars and the active M dwarf AD Leo with the goal to understand their magnetic fields. We modelled the Zeeman broadening in high-resolution near-infrared spectra of low-mass and intermediate-mass T Tauri stars and derived their mean magnetic field strengths. In intermediate-mass T Tauri stars, we only found fields weaker than 2-3 kG. However, we found that low-mass T Tauri stars can have a wide range of magnetic field strength from relatively weak fields of 1.5 kG to fields as strong as 4.4 kG, and that their field strengths do not correlate with stellar parameters. Our observations of the M dwarf AD Leo led to the first detection of linear polarization in the spectral lines of an M dwarf. We also discovered that its Stokes V profiles, which were constant over many years, had changed in our observations. We mapped its global magnetic field using ZDI and found that it became concentrated into smaller areas on the stellar surface. Finally, we analyzed Stokes IV observations of the spectroscopic binary V1878 Ori. Both components of this system are intermediate-mass T Tauri stars with very similar properties. We determined stellar parameters by studying orbital motion of the components and comparing their disentangled spectra to theoretical models. We then mapped the global magnetic fields of the two stars simultaneously using ZDI. We found that their magnetic fields have radically different geometries and different strengths.

The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by the 6.7 hour period of the Galactic magnetar candidate, 1E 161348--5055. We explore three physical mechanisms that could lead to the creation of ultra long period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares and (iii) fallback leading to long lasting accretion disks. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultra long periods, sufficiently long to accommodate FRB 180916 or 1E 161348--5055. If confirmed, such periods implicate magnetars in relatively mature states (ages $1-10$ kyr) and which possessed large internal magnetic fields at birth $B_{\rm int}\gtrsim 10^{16}$ G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars which terminate their magnetic activity at shorter periods $P \lesssim 10$ s., Comment: 12 pages, 5 figures. Accepted for publication in MNRAS

Our knowledge of the dynamics of stars has undergone a revolution thanks to the simultaneous large amount of high-quality photometric observations collected by space-based asteroseismology and ground-based high-precision spectropolarimetry. They allowed us to probe the internal rotation of stars and their surface magnetism in the whole Hertzsprung-Russell diagram. However, new methods should still be developed to probe the deep magnetic fields in those stars. Our goal is to provide seismic diagnoses that allow us to sound the internal magnetism of stars. Here, we focus on asymptotic low-frequency gravity modes and high-frequency acoustic modes. Using a first-order perturbative theory, we derive magnetic splittings of their frequencies as explicit functions of stellar parameters. As in the case of rotation, we show how asymptotic gravity and acoustic modes can allow us to probe the different components of the magnetic field in the cavities where they propagate. This demonstrates again the high potential of using mixed-modes when this is possible., Comment: 12 pages, 4 figures, accepted for publication in Astronomy and Astrophysics

The authors acknowledge support from STFC. We present models for the magnetic support of the ‘slingshot prominences’ observed in the coronae of rapidly rotating stars. We calculate mechanical equilibria of loops in a spherical geometry. Prominence-forming loops are found first for dipolar and quadrupolar stellar fields that are fully closed. Equilibria are then found within the stellar wind for a dipolar field that becomes open beyond a given radius. We identify two physical processes that may produce gaps in the distribution of prominence heights: the location of this opening radius, and the behaviour of the buoyancy force. The buoyancy may differ from one prominence-bearing loop to another if they are at different temperatures, thus potentially smearing out any gap in observed height distributions. We produce synthetic prominence distributions and compare to the observations of two well-observed stars: AB Doradus and Speedy Mic. The model recovers the more compact prominence distribution observed for Speedy Mic and reproduces better the overall shape of the height distributions for both stars when the opening radius is beyond the co-rotation radius. Publisher PDF

HD 149277 is a rare SB2 system with a slowly rotating magnetic He-rich primary with $P_{\rm rot}=25.4$ d. The CFHT/ESPaDOnS archive spectra revealed $P_{\rm orb}=11.5192 \pm 0.0005$ d indicating strong subsynchronous rotation of the primary component. Such a strong subsynchronous rotation was not detected in any other SB2 system with a magnetic chemically peculiar component. Our inspection of the spectra revealed the presence of resolved Zeeman split spectral lines allowing us to determine the variability of the mean magnetic field modulus over the rotation period. The maximum of the magnetic field modulus concides roughly with the positive extremum of the longitudinal field, whereas the minimum of the modulus with the negative extremum of the longitudinal field. No evidence for a longitudinal magnetic field was seen in the circularly polarized spectra of the secondary component. Using archival data from the ASAS3 survey, we find in the frequency spectrum only one significant peak, corresponding to the period $P_{\rm phot}=25.390 \pm 0.014$ d. This value is in good agreements with the previous determination of the rotation period, $P_{\rm rot}=25.380\pm0.007$ d, which was based on longitudinal magnetic field measurements., Comment: Accepted as MNRAS Letter

We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wide parametric search in order to understand the properties of the solutions. We confirm our previous results that the lower energy solutions have up to $\sim 25\%$ more energy than the vacuum case, helicity of the order of $\sim 5$ (in some defined units), maximum twist of $\sim 1.5$ rad, and a dipole strength that is up to $\sim 40\%$ larger than the vacuum dipole. Including the degenerate higher energy solutions allows for larger theoretical limits of up to $\sim 80\%$ more energy with respect to the vacuum case, helicity of the order of $\sim 8$, and a dipole strength that can now be up to four times that of the vacuum dipole, while the twist can be significantly larger and even diverge for configurations with disconnected domains. The higher energy solutions are probably unstable, therefore, it is unlikely that such magnetospheres exist under normal conditions in magnetars and high magnetic field pulsars., Comment: 12 pages, 9 figures; added missing reference; accepted for publication in MNRAS

Recently, high-precision optical 2~min cadence light curves obtained with \emph{TESS} for targets located in the mission's defined first four sectors have been released. The majority of these high-cadence and high-precision measurements currently span $\sim28\,{\rm d}$, thereby allowing periodic variability occurring on timescales $\lesssim14\,{\rm d}$ to potentially be detected. Magnetic chemically peculiar (mCP) A-type stars are well known to exhibit rotationally modulated photometric variability that is produced by inhomogeneous chemical abundance distributions in their atmospheres. While mCP stars typically exhibit rotation periods that are significantly longer than those of non-mCP stars, both populations exhibit typical periods $\lesssim10\,{\rm d}$; therefore, the early \emph{TESS} releases are suitable for searching for rotational modulation of the light curves of both mCP and non-mCP stars. We present the results of our search for A-type stars that exhibit variability in their \emph{TESS} light curves that is consistent with rotational modulation based on the first two data releases obtained from sectors 1 to 4. Our search yielded $134$ high-probability candidate rotational variables -- $60$ of which have not been previously reported. Approximately half of these stars are identified in the literature as Ap (mCP) stars. Comparisons between the subsample of high-probability candidate rotationally variable Ap stars and the subsample of stars that are not identified as Ap reveal that the latter subsample exhibits statistically (i) shorter rotation periods and (ii) significantly lower photometric amplitudes., Comment: 24 pages, 11 figures, 4 tables. Accepted for publication in MNRAS

Magnetic cataclysmic variables are close binaries which consist of a compact object - a white dwarf - and a red dwarf filling its Roche Lobe. Such systems are physical laboratories which enable study of the influence of magnetic fields on matter flows. They often exhibit spin-up or spin-down of the white dwarf, while some systems exhibit more complex behaviour of the spin period change. We monitor changes of the spin periods of white dwarfs in a sample of close binary systems to study interaction of the magnetic field and accretion processes as well as evolution of intermediate polars. Within the framework of our intermediate polar monitoring program, we obtained photometric CCD observations at several observatories. Two-period trigonometric polynomial fitting was used for determination of extrema timings. The (O-C) analysis was performed to study the variability of the orbital and spin periods of the systems. Using data taken during 9 years of observations of the magnetic cataclysmic variable V2306 Cygni (formerly known as 1WGA J1958.2+3232), we detected the spin period variability which shows a spin-up of the white dwarf with a characteristic time of $(53\pm5)\cdot10^4$ years. The value of the spin period was $733.33976$ seconds with the formal accuracy of $0.00015$ seconds. We derived an improved value of the orbital period of the system to be $4.371523\pm0.000009$ hours., Accepted to MNRAS 2019 July 23. 5 pages, 3 figures, supporting data online

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS., Comment: 15 pages, 3 figures, accepted for publication by MNRAS (2019 April 22)

Context. Stellar magnetic fields are one of the candidates often invoked to explain the missing transport of angular momentum observed in the models of stellar interiors. However, the properties of an internal magnetic field and the consequences of its presence on stellar evolution are largely unknown. Aims. We study the effect of an axisymmetric internal magnetic field on the frequency of gravity modes in rapidly rotating stars to check whether gravity modes can be used to detect and probe such a field. Methods. Rotation is taken into account using the traditional approximation of rotation and the effect of the magnetic field is computed using a perturbative approach. As a proof of concept, we compute frequency shifts due to a mixed (i.e. with both poloidal and toroidal components) fossil magnetic field for a representative model of a known magnetic, rapidly rotating, slowly pulsating B-type star, namely HD 43317. Results. We find that frequency shifts induced by the magnetic field scale with the square of its amplitude. A magnetic field with a near-core strength of order 150 kG (which is consistent with the observed surface field strength of order 1 kG) leads to signatures that are detectable in period spacings for high-radial-order gravity modes. Conclusions. The predicted frequency shifts can be used to constrain internal magnetic fields and offers the potential for a significant step forward in our interpretation of the observed structure of gravity-mode period spacing patterns in rapidly rotating stars., Comment: 9 pages, 5 figures, 1 table, accepted in A&A

We analyze the statistics of pulse arrival times in fast radio burst (FRB) 121102 and demonstrate that they are remarkably similar to statistics of magnetar high-energy short bursts. Motivated by this correspondence, we propose that repeating FRBs are generated during short bursts in the closed field line zone of magnetar magnetospheres via a pulsar-like emission mechanism. Crustal slippage events dislocate field line foot points, initiating intense particle acceleration and pair production, giving rise to coherent radio emission similar to that generated near pulsar polar caps. We argue that the energetics of FRB 121102 can be readily accounted for if the efficiency of the conversion of Poynting flux into coherent radio emission is $\sim10^{-4}-10^{-2}$, values consistent with empirical efficiencies of radio emission in pulsars and radio-loud magnetars. Such a mechanism could operate only in magnetars with preexisting low twist of the magnetosphere, so that the charge density in the closed zone is initially insufficient to screen the electric field provoked by the wiggling of magnetic field lines and is low enough to let $\sim 1$ GHz radio emission escape the magnetosphere, which can explain the absence of FRBs from known magnetars. The pair cascades crowd the closed flux tubes with plasma, screening the accelerating electric field, thus limiting the radio pulse duration to $\sim1$ ms. Within the framework of our model, the current dataset of the polarization angle variation in FRB 121102 suggests a magnetic obliquity $\alpha\lesssim40^\circ$ and viewing angle $\zeta$ with respect to the spin axis $\alpha, Comment: 11 pages, 4 figures, accepted to ApJ, proofs corrected

In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached. Twisting beyond that point triggers instabilities that results in the rapid expansion of magnetic field lines, followed by reconnection, as observed in previous axi-symmetric simulations. If the injection of magnetic helicity goes on, the process is recurrent, periodically releasing a similar amount of energy, of the order of a few % of the total magnetic energy. From our current distribution, we estimate the local temperature assuming that dissipation occurs mainly in the highly resistive outermost layer of the neutron star. We find that the temperature smoothly increases with injected twist, being larger for spots located in the tropical regions than in polar regions, and rather independent of their sizes. After the injection of helicity ceases, the magnetosphere relaxes to a new stable state, in which the persistent currents maintain the footprints area slightly hotter than before the onset of the instability., 6 pages, 5 figures

HD 62658 (B9p V) is a little-studied chemically peculiar star. Light curves obtained by the Kilodegree Extremely Little Telescope (KELT) and Transiting Exoplanet Survey Satellite (TESS) show clear eclipses with a period of about 4.75 d, as well as out-of-eclipse brightness modulation with the same 4.75 d period, consistent with synchronized rotational modulation of surface chemical spots. High-resolution ESPaDOnS circular spectropolarimetry shows a clear Zeeman signature in the line profile of the primary; there is no indication of a magnetic field in the secondary. PHOEBE modelling of the light curve and radial velocities indicates that the two components have almost identical masses of about 3 M$_\odot$. The primary's longitudinal magnetic field $\langle B_z \rangle$ varies between about $+100$ and $-250$ G, suggesting a surface magnetic dipole strength $B_{\rm d} = 850$~G. Bayesian analysis of the Stokes $V$ profiles indicates $B_{\rm d} = 650$~G for the primary and $B_{\rm d} < 110$ G for the secondary. The primary's line profiles are highly variable, consistent with the hypothesis that the out-of-eclipse brightness modulation is a consequence of rotational modulation of that star's chemical spots. We also detect a residual signal in the light curve after removal of the orbital and rotational modulations, which might be pulsational in origin; this could be consistent with the weak line profile variability of the secondary. This system represents an excellent opportunity to examine the consequences of magnetic fields for stellar structure via comparison of two stars that are essentially identical with the exception that one is magnetic. The existence of such a system furthermore suggests that purely environmental explanations for the origin of fossil magnetic fields are incomplete., 14 pages, 10 figures, accepted for publication in MNRAS

$\varepsilon$ Lupi A is a binary system consisting of two main sequence early B-type stars Aa and Ab in a short period, moderately eccentric orbit. The close binary pair is the only doubly-magnetic massive binary currently known. Using photometric data from the BRITE-Constellation we identify a modest heartbeat variation. Combining the photometry with radial velocities of both components we determine a full orbital solution including empirical masses and radii. These results are compared with stellar evolution models as well as interferometry and the differences discussed. We also find additional photometric variability at several frequencies, finding it unlikely these frequencies can be caused by tidally excited oscillations. We do, however, determine that these signals are consistent with gravity mode pulsations typical for slowly pulsating B stars. Finally we discuss how the evolution of this system will be affected by magnetism, determining that tidal interactions will still be dominant., Comment: 18 pages, 7 figures

Context. M dwarfs are known to generate the strongest magnetic fields among main-sequence stars with convective envelopes, but we are still lacking a consistent picture of the link between the magnetic fields and underlying dynamo mechanisms, rotation, and activity. Aims. In this work we aim to measure magnetic fields from the high-resolution near-infrared spectra taken with the CARMENES radial-velocity planet survey in a sample of 29 active M dwarfs and compare our results against stellar parameters. Methods. We used the state-of-The-Art radiative transfer code to measure total magnetic flux densities from the Zeeman broadening of spectral lines and filling factors. Results. We detect strong kG magnetic fields in all our targets. In 16 stars the magnetic fields were measured for the first time. Our measurements are consistent with the magnetic field saturation in stars with rotation periods P < 4 d. The analysis of the magnetic filling factors reveal two different patterns of either very smooth distribution or a more patchy one, which can be connected to the dynamo state of the stars and/or stellar mass. Conclusions. Our measurements extend the list of M dwarfs with strong surface magnetic fields. They also allow us to better constrain the interplay between the magnetic energy, stellar rotation, and underlying dynamo action. The high spectral resolution and observations at near-infrared wavelengths are the beneficial capabilities of the CARMENES instrument that allow us to address important questions about the stellar magnetism.© D. Shulyak et al. 2019., CARMENES is an instrument for the Centro Astronomico Hispano-Aleman de Calar Alto (CAHA, Almeria, Spain). CARMENES is funded by the German Max-Planck-Gesellschaft (MPG), the Spanish Consejo Superior de Investigaciones Cientificas (CSIC), the European Union through FEDER/ERF FICTS-2011-02 funds, and the members of the CARMENES Consortium (Max-Planck-Institut fur Astronomie, Instituto de Astrofisica de Andalucia, Landessternwarte Konigstuhl, Institut de Ciencies de l'Espai, Insitut fur Astrophysik Gottingen, Universidad Complutense de Madrid, Thuringer Landessternwarte Tautenburg, Instituto de Astrofisica de Canarias, Hamburger Sternwarte, Centro de Astrobiologia and Centro Astronomico Hispano-Aleman), with additional contributions by the Spanish Ministry of Economy, the German Science Foundation through the Major Research Instrumentation Programme and DFG Research Unit FOR2544 >Blue Planets around Red Stars>, the Klaus Tschira Stiftung, the states of Baden-Wurttemberg and Niedersachsen, and by the Junta de Andalucia. L.T.-O. acknowledges support from the Israel Science Foundation (grant No. 848/16). This work has made use of the VALD database, operated at Uppsala University, the Institute of Astronomy RAS in Moscow, and the University of Vienna. We also acknowledge the use of electronic databases SIMBAD and NASA's ADS.