Search

Early-type stars have convective cores due to a steep temperature gradient produced by the CNO cycle. These cores can host dynamos, and the generated magnetic fields can be relevant to explain the magnetism observed in Ap/Bp stars. Our main objective is to characterise the convective core dynamos and differential rotation, and to do the first quantitative analysis of the relation between magnetic activity cycle and rotation period. We use numerical 3D star-in-a-box simulations of a $2.2~M_\odot$ A-type star with a convective core of roughly $20\%$ of the stellar radius surrounded by a radiative envelope. Rotation rates from 8 to 20 days are explored. We use two models of the entire star, and an additional zoom set, where $50\%$ of the radius is retained. The simulations produce hemispheric core dynamos with cycles and typical magnetic field strengths around 60 kG. However, only a very small fraction of the magnetic energy is able to reach the surface. The cores have solar-like differential rotation, and a substantial part of the radiative envelope has quasi-rigid rotation. In the most rapidly rotating cases the magnetic energy in the core is roughly 40\% of the kinetic energy. Finally, we find that the magnetic cycle period $P_\mathrm{cyc}$ increases with decreasing the rotation period $P_\mathrm{rot}$ which is also observed in many simulations of solar-type stars. Our simulations indicate that a strong hemispherical core dynamo arises routinely, but that it is not enough the explain the surface magnetism of Ap/Bp stars. Nevertheless, as the core dynamo produces dynamically relevant magnetic fields it should not be neglected when other mechanisms are explored., Comment: Article accepted for publication in A&A

The canonical undestanding of stellar convection has recently been put under doubt due to helioseismic results and global 3D convection simulations. This "convective conundrum" is manifested by much higher velocity amplitudes in simulations at large scales in comparison to helioseismic results, and the difficulty in reproducing the solar differential rotation and dynamo with global 3D simulations. Here some aspects of this conundrum are discussed from the viewpoint of hydrodynamic Cartesian 3D simulations targeted at testing the rotational influence and surface forcing on deep convection. More specifically, the dominant scale of convection and the depths of the convection zone and the weakly subadiabatic -- yet convecting -- Deardorff zone are discussed in detail., Comment: 10 pages, 4 figures, to appear in IAU Symposium 365 proceedings (A. Getling & L. Kitchatinov, eds)

(abridged) Context: Rotation is thought to influence the size of convective eddies and the efficiency of convective energy transport in the deep convection zones of stars. Rotationally constrained convection has been invoked to explain the lack of large-scale power in observations of solar flows. Aims: The main aims are to quantify the effects of rotation on the scale of convective eddies and velocity, the depths of convective overshoot, and the subadiabatic Deardorff layers. Methods: Three-dimensional hydrodynamic simulations of rotating convection in Cartesian domains were run. The results were compared with theoretical scaling results that assume a balance between Coriolis, inertial, and buoyancy (Archemedean) forces (CIA balance). Results: The scale of convective eddies decreases as rotation increases, and ultimately reaches a rotationally constrained regime consistent with the CIA balance. Using a new measure of the rotational influence on the system, it is shown that even the deep parts of the solar convection zone are not in the rotationally constrained regime. The simulations capture the slowly and rapidly rotating scaling laws predicted by theory, and the Sun appears to be in between these two regimes. Both, the overshooting depth and the extent of the Deardorff layer, decrease as rotation becomes more rapid. For sufficiently rapid rotation the Deardorff layer is absent. Conclusions: Relating the simulations with the Sun suggests that the convective scale even in the deep parts of the Sun is only mildly affected by rotation and that some other mechanism is needed to explain the lack of strong large-scale flows in the Sun. Taking the current results at face value, the overshoot and Deardorff layers are estimated to span roughly five per cent of the pressure scale height at the base of the convection zone in the Sun., Comment: 16 pages, 14 figures, submitted to Astron. Astrophys

We review the state of the art of three dimensional numerical simulations of solar and stellar dynamos. We summarize fundamental constraints of numerical modelling and the techniques to alleviate these restrictions. Brief summary of the relevant observations that the simulations seek to capture is given. We survey the current progress of simulations of solar convection and the resulting large-scale dynamo. We continue to studies that model the Sun at different ages and to studies of stars of different masses and evolutionary stages. Both simulations and observations indicate that rotation, measured by the Rossby number which is the ratio of rotation period and convective turnover time, is a key ingredient in setting the overall level and characteristics of magnetic activity. Finally, efforts to understand global 3D simulations in terms of mean-field dynamo theory are discussed., Comment: 70 pages, 9 figures, submitted to Space Science Reviews

The aim of this study is to explore the magnetic and flow properties of fully convective M dwarfs as a function of rotation period Prot and magnetic Reynolds ReM and Prandlt numbers PrM. We performed three-dimensional simulations of fully convective stars using a star-in-a-box setup. This setup allows global dynamo simulations in a sphere embedded in a Cartesian cube. The equations of non-ideal magnetohydrodynamics were solved with the Pencil Code. We used the stellar parameters of an M5 dwarf with 0.21M_odot at three rotation rates corresponding to rotation periods (Prot): 43, 61 and 90 days, and varied the magnetic Prandtl number in the range from 0.1 to 10. We found systematic differences in the behaviour of the large-scale magnetic field as functions of rotation and PrM. For the simulations with Prot = 43 days and PrM <= 2, we found cyclic large-scale magnetic fields. For PrM > 2 the cycles vanish and field shows irregular reversals. In simulations with Prot = 61 days for PrM <= 2 the cycles are less clear and the reversal are less periodic. In the higher-PrM cases, the axisymmetric mean field shows irregular variations. For the slowest rotation case with Prot = 90 days, the field has an important dipolar component for PrM <= 5. For the highest PrM the large-scale magnetic field is predominantly irregular at mid-latitudes, with quasi-stationary fields near the poles. For the simulations with cycles, the cycle period length slightly increases with increasing ReM., Comment: 11 pages, 12 figures, submitted to A&A

The centrifugal force is often omitted from simulations of stellar convection either for numerical reasons or because it is assumed to be weak compared to the gravitational force. However, it might be an important factor in rapidly rotating stars, such as solar analogs, due to its $\Omega^2$ scaling, where $\Omega$ is the rotation rate of the star. We study the effects of the centrifugal force in a set of 21 semi-global stellar dynamo simulations with varying rotation rates. Included in the set are three control runs aimed at distinguishing the effects of the centrifugal force from the nonlinear evolution of the solutions. We decomposed the magnetic field into spherical harmonics and studied the migration of azimuthal dynamo waves (ADWs), the energy of different large-scale magnetic modes, and differential rotation. In the regime with the lowest rotation rates, $\Omega = 5-10\Omega_\odot$, where $\Omega_\odot$ is the rotation rate of the Sun, we see no marked changes in either the differential rotation or the magnetic field properties. For intermediate rotation, $\Omega = 20-25\Omega_\odot$, we identify an increase in the differential rotation as a function of centrifugal force. The axisymmetric magnetic energy tends to decrease with centrifugal force, while the non-axisymmetric one increases. The ADWs are also affected, especially in the propagation direction. In the most rapidly rotating set with $\Omega=30\Omega_\odot$, these changes are more pronounced, and in one case the propagation direction of the ADW changes from prograde to retrograde. The control runs suggest that the results are a consequence of the centrifugal force and not due to the details of the initial conditions or the history of the run. We find that the differential rotation and the ADWs only change as a function of the centrifugal force when rotation is rapid enough., Comment: 8 pages, 7 figures, accepted for publication in A&A

According to our understanding of stellar evolution, early-type stars have radiative envelopes and convective cores due to a steep temperature gradient produced by the CNO cycle. Some of these stars (mainly, the subclasses Ap and Bp) have strong magnetic fields, enough to be directly observed using the Zeeman effect. Here, we present 3D magnetohydrodynamic simulations of an $2 ~M_{\odot}$ A-type star using the star-in-a-box model. Our goal is to explore if the modeled star is able to maintain a magnetic field as strong as the observed ones, via a dynamo driven by its convective core, or via maintaining a stable fossil field configuration coming from its early evolutionary stages, using different rotation rates. We created two models, a partially radiative and a fully radiative one, which are determined by the value of the heat conductivity. Our model is able to explore both scenarios, including convection-driven dynamos., Comment: Conference proceeding (Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia), 3 pages, 3 figures

Turbulent motions enhance the diffusion of large-scale flows and temperature gradients. Such diffusion is often parameterized by coefficients of turbulent viscosity ($\nu_{\rm t}$) and turbulent thermal diffusivity ($\chi_{\rm t}$) that are analogous to their microscopic counterparts. We compute the turbulent diffusion coefficients by imposing large-scale velocity and temperature gradients on a turbulent flow and measuring the response of the system. We also confirm our results using experiments where the imposed gradients are allowed to decay. To achieve this, we use weakly compressible three-dimensional hydrodynamic simulations of isotropically forced homogeneous turbulence. We find that the turbulent viscosity and thermal diffusion, as well as their ratio the turbulent Prandtl number, ${\rm Pr}_{\rm t} = \nu_{\rm t}/\chi_{\rm t}$, approach asymptotic values at sufficiently high Reynolds and Pecl\'et numbers. We also do not find a significant dependence of ${\rm Pr}_{\rm t}$ on the microscopic Prandtl number ${\rm Pr} = \nu/\chi$. These findings are in stark contrast to results from the $k-\epsilon$ model which suggests that ${\rm Pr}_{\rm t}$ increases monotonically with decreasing ${\rm Pr}$. The current results are relevant for the ongoing debate of, for example, the nature of the turbulent flows in the very low ${\rm Pr}$ regimes of stellar convection zones., Comment: 10 pages, 5 figures, submitted to JFM Rapids

(abridged) Context: Late-type stars such as the Sun rotate differentially due to the interaction of turbulent convection and rotation. Aims: The aim of the study is to investigate the effects of the thermal Prandtl number on the transition from anti-solar (slow equator, fast poles) to solar-like (fast equator, slow poles) differential rotation. Methods: Three-dimensional hydrodynamic and magnetohydrodynamic simulations in semi-global spherical wedge geometry are used to model convection zones of solar-like stars. Results: The overall convective velocity amplitude increases as the Prandtl number decreases in accordance with earlier studies. The transition from anti-solar to solar-like differential rotation is insensitive to the Prandtl number for Prandtl numbers below unity but for Prandtl numbers greater than unity, solar-like differential rotation becomes significantly harder to excite. Magnetic fields and more turbulent regimes with higher fluid and magnetic Reynolds numbers help in achieving solar-like differential rotation in near-transition cases. Solar-like differential rotation occurs only in cases with radially outward angular momentum transport at the equator. The dominant contribution to such outward transport near the equator is due to prograde propagating thermal Rossby waves. Conclusions: The differential rotation is sensitive to the Prandtl number only for large Prandtl numbers in the parameter regime explored in the current study. Magnetic fields have a greater effect on the differential rotation, although the inferred presence of a small-scale dynamo does not lead to drastically different results in the present study. The dominance of the thermal Rossby waves in the simulations is puzzling given the non-detection in the Sun. The current simulations are shown to be incompatible with the currently prevailing mean-field theory of differential rotation., Comment: 11 pages, 8 figures, submitted to Astron. Astrophys

M dwarfs are low-mass main-sequence stars, the most numerous type of stars in the solar neighbourhood, which are known to have significant magnetic activity. The aim of this work is to explore the dynamo solutions and magnetic fields of fully convective M dwarfs with varying magnetic Prandtl numbers ${\rm Pr_M}$, and a rotation period (${\rm P_{rot}}$) of 43 days. ${\rm Pr_M}$ is known to play an important role in the dynamo action; dynamos for low-$\rm Pr_M$ and large-$\rm Pr_M$ have very different properties. We performed three-dimensional magnetohydrodynamical (MHD) numerical simulations with the ``star-in-a-box'' model using stellar parameters for an M5 dwarf with $0.21 M_{\odot}$. We found that the dynamo solutions are sensitive to ${\rm Pr_M}$. The simulations at this rotation period present periodic cycles of the large-scale magnetic field up to ${\rm Pr_M} \leq 2$; for higher values the cycles disappear and irregular solutions to arise. Our results are consistent with previous studies and suggest that the dynamos operating in fully convective stars behave similarly as those in partially convective stars., Comment: 3 pages, 1 figure, Accepted for publication in Bolet\'in de la Asociaci\'on Argentina de Astronom\'ia

Eclipsing time variations (ETVs) in post-common-envelope binaries (PCEBs) were proposed to be due to the time-varying component of the stellar gravitational quadrupole moment. This is suggested to be produced by changes in the stellar structure due to an internal redistribution of angular momentum and the effect of the centrifugal force. We examined this hypothesis and present 3D simulations of compressible magnetohydrodynamics (MHD) performed with the {\sc Pencil Code}. We modeled the stellar dynamo for a solar-mass star with angular velocities of 20 and 30 times solar. We included and varied the strength of the centrifugal force and compared the results with reference simulations without the centrifugal force and with a simulation in which its effect is enhanced. The centrifugal force causes perturbations in the evolution of the numerical model, so that the outcome in the details becomes different as a result of nonlinear evolution. While the average density profile is unaffected by the centrifugal force, a relative change in the density difference between high latitudes and the equator of $\sim10^{-4}$ is found. The power spectrum of the convective velocity is found to be more sensitive to the angular velocity than to the strength of the centrifugal force. The quadrupole moment of the stars includes a fluctuating and a time-independent component, which vary with the rotation rate. As very similar behavior is produced in absence of the centrifugal force, we conclude that it is not the main ingredient for producing the time-averaged and fluctuating quadrupole moment of the star. In a real physical system, we thus expect contributions from both components, that is, from the time-dependent gravitational force from the variation in the quadrupole term and from the spin-orbit coupling that is due to the persistent part of the quadrupole., Comment: 8 pages, 7 figures. Accepted in A&A

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. 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, , and Rossby number, $Ro$, resembling the well-studied rotation-activity relation. Among the slowly rotating stars, we find that magnetic flux, $\Phi_\textrm{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$\alpha$ 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 $\langle B \rangle \approx 800$ G while H$\alpha$ 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., Comment: A&A accepted, abstract abbreviated, 20 pages, Table B.1 included in arXiv version, full MCMC posteriors and linefit plots are available at http://carmenes.cab.inta-csic.es/

Magnetohydrodynamic star-in-a-box simulations of convection and dynamos in a solar-like star with different rotation rates are presented. These simulations produce solar-like differential rotation with a fast equator and slow poles, and magnetic activity that resembles that of the Sun with equatorward migrating activity at the surface. Furthermore, the ratio of rotation to cycle period is almost constant as the rotation period decreases in the limited sample considered here. This is reminiscent of the suggested inactive branch of stars from observations and differs from most earlier simulation results from spherical shell models. While the exact excitation mechanism of the dynamos in the current simulations is not yet clear, it is plausible that the greater freedom that the magnetic field has due to the inclusion of the radiative core and regions exterior to the star are important in shaping the dynamo., Comment: 8 pages, 4 figures, submitted to AAS journals. Animation related to Fig. 2 is available as ancillary material

The radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical mean-field models have been successful in explaining this property of the solar NSSL, while global direct convection models have been unsuccessful. We investigate reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under NSSL conditions. Simulations have minimal ingredients. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. Parameters of the simulations are chosen such they match the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch is varied by changing the rotation rate such that the resulting Coriolis numbers<1. We measure the turbulent transport coefficient relevant for the non-diffusive and diffusive parts of the Reynolds stress and compare them with predictions of current mean-field theories. A negative radial gradient of mean flow similar to the solar NSSL is generated only at the equator where meridional flows are absent. At other latitudes the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50\% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional stresses are neglected. However, the rotational dependence of turbulent viscosity in our simulations is in good agreement with theoretical prediction. Our results are in qualitative agreement with global convection simulations in that a NSSL obtained near the equator.

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. 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. 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 vs. non-axisymmetric, and poloidal vs. toroidal components, and compare them to the Rossby numbers of the stars. We present five new ZDI-maps of young, solar-type stars from Dec 2017. Of special interest is the case of V1358 Ori, that has 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 more axisymmetric field, and possibly 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 similar levels as in earlier observations. We find evidence for a magnetic polarity reversal having occurred in V1358 Ori. \c{hi}1 Ori could be an interesting target for future observations, with a possible 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 each other out, and are not resolved in our maps.

The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid., Comment: 7 pages, submitted to the Journal for Open Source Software (JOSS)

One interpretation of the activity and magnetism of late-type stars is that these both intensify with decreasing Rossby number up to a saturation level, suggesting that stellar dynamos depend on both rotation and convective turbulence. Some studies have claimed, however, that rotation alone suffices to parametrise this scaling adequately. Here, we tackle the question of the relevance of turbulence to stellar dynamos by including evolved, post main sequence stars in the analysis of the rotation-activity relation. These stars rotate very slowly compared with main sequence stars, but exhibit similar activity levels. We show that the two evolutionary stages fall together in the rotation-activity diagram and form a single sequence in the unsaturated regime in relation only to Rossby numbers derived from stellar models, confirming earlier preliminary results that relied on a more simplistic parametrisation of the convective turnover time. This mirrors recent results of fully convective M dwarfs, which likewise fall on the same rotation-activity sequence as partially convective solar-type stars. Our results demonstrate that turbulence plays a crucial role in driving stellar dynamos and suggest that there is a common turbulence-related dynamo mechanism explaining the magnetic activity of all late-type stars., Comment: This is the initial submission version of the paper "Common dynamo scaling in slowly rotating young and evolved stars" (Lehtinen et al. 2020). The final version accepted in Nature Astronomy can be found at https://www.nature.com/articles/s41550-020-1039-x

In this study we present a compressible test-field method (CTFM) for computing $\alpha$ effect and turbulent magnetic diffusivity tensors, as well as those relevant for mean ponderomotive force and mass source, applied to the full MHD equations. We describe the theoretical background of the method, and compare it to the quasi-kinematic test-field method, and to the previously studied variant working in simplified MHD (SMHD). We present several test cases using velocity and magnetic fields of the Roberts geometry, and also compare with the imposed-field method. We show that, for moderate imposed field strengths, the nonlinear CTFM (nCTFM) gives results in agreement with the imposed-field method. Comparison of different flavors of the nCTFM in the shear dynamo case also agree up to equipartition field strengths. Some deviations between the CTFM and SMHD variants exist. As a relevant physical application, we study non-helically forced shear flows, which exhibit large-scale dynamo action, and present a re-analysis of low Reynolds number, moderate shear systems, where we previously neglected the pressure gradient in the momentum equation, and found no coherent shear-current effect. Another key difference is that in the earlier study we used magnetic forcing to mimic small-scale dynamo action, while here it is self-consistently driven by purely kinetic forcing. The kinematic CTFM with general validity forms the core of our analysis. We still find no coherent shear-current effect, but do recover strong large-scale dynamo action that, according to our analysis, is driven through the incoherent effects., Comment: 23 pages, 11 figures, submitted to the Astrophysical Journal

(abridged) Context: The ratio of kinematic viscosity to thermal diffusivity, the Prandtl number, is much smaller than unity in stellar convection zones. Aims: To study the statistics of convective flows and energy transport as functions of the Prandtl number. Methods: Three-dimensional numerical simulations convection in Cartesian geometry are used. The convection zone (CZ) is embedded between two stably stratified layers. Statistics and transport properties of up- and downflows are studied separately. Results: The rms velocity increases with decreasing Prandtl number. At the same time the filling factor of downflows decreases and leads to stronger downflows at lower Prandtl numbers, and to a strong dependence of overshooting on the Prandtl number. Velocity power spectra do not show marked changes as a function of Prandtl number. At the highest Reynolds numbers the velocity power spectra are compatible with the Bolgiano-Obukhov $k^{-11/5}$ scaling. The horizontally averaged convected energy flux ($\overline{F}_{\rm conv}$) is independent of the Prandtl number within the CZ. However, the upflows (downflows) are the dominant contribution to the convected flux at low (high) Prandtl number. These results are similar to those from Rayleigh-Ben\'ard convection in the low Prandtl number regime where convection is vigorously turbulent but inefficient at transporting energy. Conclusions: The current results indicate a strong dependence of convective overshooting and energy flux on the Prandtl number. Numerical simulations of astrophysical convection often use Prandtl number of unity. The current results suggest that this can lead to misleading results and that the astrophysically relevant low Prandtl number regime is qualitatively different from the parameters regimes explored in typical simulations., Comment: 15 pages, 20 figures, submitted to Astron. Astrophys

The role of turbulent effects for dynamos in the Sun and stars continues to be debated. Mean-field (MF) theory provides a broadly used framework to connect these effects to fundamental magnetohydrodynamics. While inaccessible observationally, turbulent effects can be directly studied using global convective dynamo (GCD) simulations. We measure the turbulent effects in terms of turbulent transport coefficients, based on the MF framework, from an exemplary GCD simulation using the test-field method. These coefficients are then used as an input into an MF model. We find a good agreement between the MF and GCD solutions, which validates our theoretical approach. This agreement requires all turbulent effects to be included, even those which have been regarded as unimportant so far. Our results suggest that simple dynamo models, as are commonly used in the solar and stellar community, relying on very few, precisely fine-tuned turbulent effects, may not be representative of the full dynamics of dynamos in global convective simulations and astronomical objects., Comment: 7 pages with 3 figures, plus 3 pages and 2 figures in the appendix, published in The Astrophysical Journal Letters

In this study we present a dynamical agent-based model to investigate the interplay between the socio-economy of and SEIRS-type epidemic spreading over a geographical area, divided to smaller area districts and further to smallest area cells. The model treats the populations of cells and authorities of districts as agents, such that the former can reduce their economic activity and the latter can recommend economic activity reduction both with the overall goal to slow down the epidemic spreading. The agents make decisions with the aim of attaining as high socio-economic standings as possible relative to other agents of the same type by evaluating their standings based on the local and regional infection rates, compliance to the authorities' regulations, regional drops in economic activity, and efforts to mitigate the spread of epidemic. We find that the willingness of population to comply with authorities' recommendations has the most drastic effect on the epidemic spreading: periodic waves spread almost unimpeded in non-compliant populations, while in compliant ones the spread is minimal with chaotic spreading pattern and significantly lower infection rates. Health and economic concerns of agents turn out to have lesser roles, the former increasing their efforts and the latter decreasing them., Comment: 13 pages, 6 figures. Submitted to Physica A: Statistical Mechanics and Its Applications

Obtaining observational constraints on the role of turbulent effects for the solar dynamo is a difficult, yet crucial, task. Without such knowledge, the full picture of the operation mechanism of the solar dynamo cannot be formed. The magnetic helicity spectrum provides important information about the $\alpha$ effect. Here we demonstrate a formalism in spherical geometry to infer magnetic helicity spectra directly from observations of the magnetic field, taking into account the sign change of magnetic helicity across the Sun's equator. Using an angular correlation function of the magnetic field, we develop a method to infer spectra for magnetic energy and helicity. The retrieval of the latter relies on a fundamental definition of helicity in terms of linkage of magnetic flux. We apply the two-scale approach, previously used in Cartesian geometry, to spherical geometry for systems where a sign reversal of helicity is expected across the equator at both small and large scales. We test the method by applying it to an analytical model of a fully helical field, and to magneto-hydrodynamic simulations of a turbulent dynamo. The helicity spectra computed from the vector potential available in the models are in excellent agreement to the spectra computed solely from the magnetic field using our method. In a next test, we use our method to obtain the helicity spectrum from a synoptic magnetic field map corresponding to a Carrington rotation. We observe clear signs of a bihelical spectrum of magnetic helicity. Our formalism makes it possible to infer magnetic helicity in spherical geometry, without the necessity of computing the magnetic vector potential. This has the advantage of being gauge invariant. It has many applications in solar and stellar observations, but can also be used to analyze global magnetoconvection models of stars and compare them with observations., Comment: 10 pages, 6 figures, 3 appendices, submitted to Astronomy and Astrophysics

Modern compute nodes in high-performance computing provide a tremendous level of parallelism and processing power. However, as arithmetic performance has been observed to increase at a faster rate relative to memory and network bandwidths, optimizing data movement has become critical for achieving strong scaling in many communication-heavy applications. This performance gap has been further accentuated with the introduction of graphics processing units, which can provide by multiple factors higher throughput in data-parallel tasks than central processing units. In this work, we explore the computational aspects of iterative stencil loops and implement a generic communication scheme using CUDA-aware MPI, which we use to accelerate magnetohydrodynamics simulations based on high-order finite differences and third-order Runge-Kutta integration. We put particular focus on improving intra-node locality of workloads. Our GPU implementation scales strongly from one to $64$ devices at $50\%$--$87\%$ of the expected efficiency based on a theoretical performance model. Compared with a multi-core CPU solver, our implementation exhibits $20$--$60\times$ speedup and $9$--$12\times$ improved energy efficiency in compute-bound benchmarks on $16$ nodes., Comment: 15 pages, 15 figures. Updated with the accepted manuscript. More extensive tests added and wording clarified in several places. Please refer to the published article for the most polished version

The possibility to detect circumbinary planets and to study stellar magnetic fields through eclipsing time variations (ETVs) in binary stars has sparked an increase of interest in this area of research. We revisit the connection between stellar magnetic fields and the gravitational quadrupole moment $Q_{xx}$ and compare different dynamo-generated ETV models with our simulations. We present magnetohydrodynamical simulations of solar mass stars with rotation periods of 8.3, 1.2, and 0.8 days and perform a detailed analysis of the magnetic and quadrupole moment using spherical harmonic decomposition. The extrema of $Q_{xx}$ are associated with changes in the magnetic field structure. This is evident in the simulation with a rotation period of 1.2 days. Its magnetic field has a more complex behavior than in the other models, as the large-scale nonaxisymmetric field dominates throughout the simulation and the axisymmetric component is predominantly hemispheric. This triggers variations in the density field that follow the magnetic field asymmetry with respect to the equator, affecting the $zz$ component of the inertia tensor, and thus modulating $Q_{xx}$. The magnetic field of the two other runs are less variable in time and more symmetric with respect to the equator, such that the variations in the density are weaker, and therefore only small variations in $Q_{xx}$ are seen. If interpreted via the classical Applegate mechanism (tidal locking), the quadrupole moment variations obtained in the current simulations are about two orders of magnitude below those deduced from observations of post-common-envelope binaries. However, if no tidal locking is assumed, our results are compatible with the observed ETVs., Comment: 15 pages, 16 figures. Accepted for publication in A&A

To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. We present here the first test-field measurements from our higher-resolution runs with improved heat conduction description. They indicate significant changes in the profiles of the most crucial inductive effect related to solar and stellar dynamo mechanisms. Higher resolution runs, currently undertaken, will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way., Comment: Originally published in "High Performance Computing Science and Engineering, Garching/Munich, 2020". Editors: P. Bastina, D. Kranzm\"uller, H. Br\"euchle, M. Brehm. ISBN: 978-3-9816675-4-7

Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale (LSD) and small-scale dynamos (SSD) is not fully understood. We performed a systematic study of turbulent dynamos driven by isotropic forcing in isothermal MHD with magnetic Prandtl number of unity, focusing on the exponential growth stage. Both helical and non-helical forcing was employed to separate the effects of LSD and SSD in a periodic domain. Reynolds numbers (Rm) up to $\approx 250$ were examined and multiple resolutions used for convergence checks. We ran our simulations with the Astaroth code, designed to accelerate 3D stencil computations on graphics processing units (GPUs) and to employ multiple GPUs with peer-to-peer communication. We observed a speedup of $\approx 35$ in single-node performance compared to the widely used multi-CPU MHD solver Pencil Code. We estimated the growth rates both from the averaged magnetic fields and their power spectra. At low Rm, LSD growth dominates, but at high Rm SSD appears to dominate in both helically and non-helically forced cases. Pure SSD growth rates follow a logarithmic scaling as a function of Rm. Probability density functions of the magnetic field from the growth stage exhibit SSD behaviour in helically forced cases even at intermediate Rm. We estimated mean-field turbulence transport coefficients using closures like the second-order correlation approximation (SOCA). They yield growth rates similar to the directly measured ones and provide evidence of $\alpha$ quenching. Our results are consistent with the SSD inhibiting the growth of the LSD at moderate Rm, while the dynamo growth is enhanced at higher Rm., Comment: 22 pages, 23 figures, 2 tables, Accepted for publication in the Astrophysical Journal

(abridged) Context: Main-sequence late-type stars with masses less than $0.35 M_\odot$ are fully convective. Aims: The goal is to study convection, differential rotation, and dynamos as functions of rotation in fully convective stars. Methods: Three-dimensional hydrodynamic and magnetohydrodynamic numerical simulations with a star-in-a-box model, where a spherical star is immersed inside of a Cartesian cube, are used. The model corresponds to a $0.2M_\odot$ M5 dwarf. Rotation periods ($P_{\rm rot}$) between 4.3 and 430 days are explored. Results: The slowly rotating model with $P_{\rm rot}=430$ days produces anti-solar differential rotation with a slow equator and fast poles, along with predominantly axisymmetric quasi-steady large-scale magnetic fields. For intermediate rotation ($P_{\rm rot}=144$ and $43$ days) differential rotation is solar-like (fast equator, slow poles) and large-scale magnetic fields are mostly axisymmetric and either quasi-stationary or cyclic. The latter occurs in a similar parameter regime as in other numerical studies in spherical shells, and the cycle period is similar to observed cycles in fully convective stars with comparable $P_{\rm rot}$. In the rapid rotation regime the differential rotation is weak and the large-scale magnetic fields are increasingly non-axisymmetric with a dominating $m=1$ mode. This large-scale non-axisymmetric field also exhibits azimuthal dynamo waves. Conclusions: The results of the star-in-a-box models agree with simulations of partially convective late-type stars in spherical shells in that the transitions in differential rotation and dynamo regimes occur at similar rotational regimes in terms of the Coriolis (inverse Rossby) number. This similarity between partially and fully convective stars suggests that the processes generating differential rotation and large-scale magnetism are insensitive to the geometry of the star., Comment: 17 pages, 11 figures, submitted to Astron. Astrophys, revised as per referee report

Global and semi-global convective dynamo simulations of solar-like stars are known to show a transition from an anti-solar (fast poles, slow equator) to solar-like (fast equator, slow poles) differential rotation (DR) for increasing rotation rate. The dynamo solutions in the latter regime can exhibit regular cyclic modes, whereas in the former one, only stationary or temporally irregular solutions have been obtained so far. In this paper we present a semi-global dynamo simulation in the transition region, exhibiting two coexisting dynamo modes, a cyclic and a stationary one, both being dynamically significant. We seek to understand how such a dynamo is driven by analyzing the large-scale flow properties (DR and meridional circulation) together with the turbulent transport coefficients obtained with the test-field method. Neither an $\alpha\Omega$ dynamo wave nor an advection-dominated dynamo are able to explain the cycle period and the propagation direction of the mean magnetic field. Furthermore, we find that the $\alpha$ effect is comparable or even larger than the $\Omega$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $\alpha^2\Omega$ or $\alpha^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping., Comment: 13 pages, 7 figures. Published in ApJ

(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${\rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${\rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${\rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${\rm Re}$ whereas for small ${\rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${\rm Pm}_{\rm t}$ found from the simulations ($\approx0.9\ldots0.95$) at large ${\rm Re}$ and scale separation ratio is higher than any of the analytic predictions., Comment: 13 pages, 12 figures, accepted in Astron. Astrophys

The magnetic activity of late-type stars is correlated with their rotation rates. Up to a certain limit, stars with smaller Rossby numbers, defined as the rotation period divided by the convective turnover time, have higher activity. A more detailed look at this rotation-activity relation reveals that, rather than being a simple power law relation, the activity scaling has a shallower slope for the low-Rossby stars than for the high-Rossby ones. We find that, for the chromospheric CaII H&K activity, this scaling relation is well modelled by a broken two-piece power law. Furthermore, the knee-point of the relation coincides with the axisymmetry to non-axisymmetry transition seen in both the spot activity and surface magnetic field configuration of active stars. We interpret this knee-point as a dynamo transition between dominating axi- and non-axisymmetric dynamo regimes with a different dependence on rotation and discuss this hypothesis in the light of current numerical dynamo models., Comment: 10 pages, 5 figures, 2 tables

The possibility of explaining shear flow dynamos by a magnetic shear-current (MSC) effect is examined via numerical simulations. Our primary diagnostics is the determination of the turbulent magnetic diffusivity tensor $\boldsymbol{\eta}$. In our setup, a negative sign of its component $\eta_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficients from systems with magnetic background turbulence, we present an extension of the test-field method (TFM), applicable to our setup where the pressure gradient is dropped from the momentum equation: the nonlinear TFM (NLTFM). Our momentum equation is related to Burgers' equation and the resulting flows are referred to as magnetized burgulence. We use both stochastic kinetic and magnetic forcings to mimic cases without and with simultaneous small-scale dynamo action (SSD). When we force only kinetically, negative $\eta_{yx}$ are obtained with exponential growth in both the radial and azimuthal mean magnetic field components. Using isotropic magnetic forcing, the field growth is no longer exponential, while NLTFM yields positive $\eta_{yx}$. By employing an alternative forcing from which wavevectors having small components are removed, the exponential growth is recovered, but the NLTFM results do not change significantly. Analyzing the dynamo excitation conditions for the coherent SC and incoherent $\alpha$ and SC effects shows that the incoherent effects are the main drivers of the dynamo in the majority of cases. We find no evidence for MSC-effect-driven dynamos in our simulations., Comment: 21 pages, 11 figures, published in Astrophysical Journal

Context. Results from global magnetoconvection simulations of solar-like stars are at odds with observations in many respects: They show a surplus of energy in the kinetic power spectrum at large scales, anti-solar differential rotation profiles, with accelerated poles and a slow equator, for the solar rotation rate, and a transition from axi- to non-axisymmetric dynamos at a much lower rotation rate than what is observed. Even though the simulations reproduce the observed active longitudes in fast rotators, their motion in the rotational frame (the so-called azimuthal dynamo wave, ADW) is retrograde, in contrast to the prevalent prograde motion in observations. Aims. We study the effect of a more realistic treatment of heat conductivity in alleviating the discrepancies between observations and simulations. Methods. We use physically-motivated heat conduction, by applying Kramers opacity law, on a semi-global spherical setup describing convective envelopes of solar-like stars, instead of a prescribed heat conduction profile from mixing-length arguments. Results. We find that some aspects of the results now better correspond to observations: The axi- to non-axisymmetric transition point is shifted towards higher rotation rates. We also find a change in the propagation direction of ADWs so that also prograde waves are now found. The transition from anti-solar to solar-like rotation profile, however, is also shifted towards higher rotation rates, leaving the models into an even more unrealistic regime. Conclusions. Although a Kramers-based heat conduction does not help in reproducing the solar rotation profile, it does help in the faster rotation regime, where the dynamo solutions now match better with observations., Comment: 10 pages, 6 figures, 1 appendix. Submitted to A&A

Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other stellar parameters. We also study the differences in the shape of the solar cycle, depending on what activity indicator is used. The observations are also compared to simulated activity cycles. We use the chromospheric Ca II H&K data from the Mount Wilson Observatory HK Project. From this data set we identify 47 individual cycles from 18 stars. We use the statistical skewness of a cycle as a measure of its asymmetry, and compare this to other stellar parameters. A similar analysis has been done to magnetic cycles extracted from direct numerical magnetohydrodynamic simulations of solar-type convection zones. The shape of the solar cycle (fast rise and slower decline) is common in other stars as well, although the Sun has particularly asymmetric cycles. Cycle-to-cycle variations are large, but the average shape of a cycle is still fairly well represented by a sinusoid. We find only slight correlations between the cycle asymmetry and other stellar parameters. There are large differences in the shape of the solar cycle, depending on what activity indicator is used. In the simulated cycles, there is a difference in the symmetry of global simulations covering the full longitudinal range, hence capable of exciting non-axisymmetric dynamo modes, versus wedge simulations covering a partial extent in longitude, where only axisymmetric modes are possible. The former produce preferentially positive skewness, while the latter a negative one., Comment: Accepted to A&A; 12 pages, 12 figures

The alpha effect is believed to play a key role in the generation of the solar magnetic field. A fundamental test for its significance in the solar dynamo is to look for magnetic helicity of opposite signs in the two hemispheres, and at small and large scales. However, measuring magnetic helicity is compromised by the inability to fully infer the magnetic field vector from observations of solar spectra, caused by what is known as the "pi ambiguity" of spectropolarimetric observations. We decompose linear polarisation into parity-even and parity-odd E and B polarisations, which are not affected by the "pi ambiguity". Furthermore, we study whether the correlations of spatial Fourier spectra of B and parity-even quantities such as E or temperature T are a robust proxy for magnetic helicity of solar magnetic fields. We analyse polarisation measurements of active regions observed by the Helioseismic and Magnetic Imager on board the Solar Dynamics observatory. Theory predicts the magnetic helicity of active regions to have, statistically, opposite signs in the two hemispheres. We then compute the parity-odd E B and T B correlations, and test for systematic preference of their sign based on the hemisphere of the active regions. We find that: (i) E B and T B correlations are a reliable proxy for magnetic helicity, when computed from linear polarisation measurements away from spectral line cores, and (ii) E polarisation reverses its sign close to the line core. Our analysis reveals Faraday rotation to not have a significant influence on the computed parity-odd correlations. The EB decomposition of linear polarisation appears to be a good proxy for magnetic helicity independent of the "pi ambiguity". This allows us to routinely infer magnetic helicity directly from polarisation measurements., Comment: 10 pages, 7 figures, 2 appendices with 2 figures, published in Astronomy & Astrophysics

Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the $f$-mode, which emerges from the density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localized bipolar magnetic structures in the middle of the domain. Some notable findings are: (i) compared to the purely hydrodynamic case, the strength of the $f$-mode is significantly larger at high horizontal wavenumbers $k$, but the fanning is weaker for the localized subsurface magnetic field concentrations investigated here than the periodic ones studied earlier; (ii) when the strength of the magnetic field is enhanced at a fixed depth below the surface, the fanning of the $f$-mode in the $k\omega$ diagram increases proportionally in such a way that the normalized $f$-mode strengths remain nearly the same in different such cases; (iii) the unstable Bloch modes reported previously in case of harmonically varying magnetic fields are now completely absent when more realistic localized magnetic field concentrations are imposed beneath the surface, thus suggesting that the Bloch modes are unlikely to be supported during most phases of the solar cycle; (iv) the $f$-mode strength appears to depend also on the depth of magnetic field concentrations such that it shows a relative decrement when the maximum of the magnetic field is moved to a deeper layer. We argue that detections of $f$-mode perturbations such as those being explored here could be effective tracers of solar magnetic fields below the photosphere before these are directly detectable as visible manifestations in terms of active regions or sunspots., Comment: 15 pages, 7 figures, revised version submitted to GAFD special issue on "Physics and Algorithms of the Pencil Code"

We test the sensitivity of hydrodynamic and magnetohydrodynamic turbulent convection simulations with respect to Mach number, thermal and magnetic boundary conditions, and the centrifugal force. We find that varying the luminosity, which also controls the Mach number, has only a minor effect on the large-scale dynamics. A similar conclusion can also be drawn from the comparison of two formulations of the lower magnetic boundary condition with either vanishing electric field or current density. The centrifugal force has an effect on the solutions, but only if its magnitude with respect to acceleration due to gravity is by two orders of magnitude greater than in the Sun. Finally, we find that the parameterisation of the photospheric physics, either by an explicit cooling term or enhanced radiative diffusion, is more important than the thermal boundary condition. In particular, runs with cooling tend to lead to more anisotropic convection and stronger deviations from the Taylor-Proudman state. In summary, the fully compressible approach taken here with the Pencil Code is found to be valid, while still allowing the disparate timescales to be taken into account., Comment: 26 pages, 14 figures, accepted to GAFD Special issue on 'Physics and Algorithms of the Pencil Code'

Using a recently developed two-scale formalism to determine the magnetic helicity spectrum (Brandenburg et al. 2017), we analyze synoptic vector magnetograms built with data from the Vector Spectromagnetograph (VSM) instrument on the \emph{Synoptic Optical Long-term Investigations of the Sun} (SOLIS) telescope during January 2010-July 2016. In contrast to an earlier study using only three Carrington rotations, our analysis includes 74 synoptic Carrington rotation maps. We recover here bihelical spectra at different phases of solar cycle~24, where the net magnetic helicity in the majority of the data is consistent with a large-scale dynamo with helical turbulence operating in the Sun. More than $20\%$ of the analyzed maps, however, show violations of the expected sign rule., Comment: 10 pages, 10 figures, submitted to the Astrophysical Journal

We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We find that the convective heat transport is strongly concentrated toward the equatorial and polar regions in the cases without a substantial radiative layer below the convection zone. The presence of a stable layer below the convection zone significantly reduces the anisotropy of radial enthalpy transport. Furthermore, the dynamo solutions are sensitive to subtle changes in the convection zone structure. We find that the kinetic helicity changes sign in the deeper parts of the convection zone at high latitudes in all runs. This region expands progressively toward the equator in runs with a thicker stably stratified layer., Comment: 32 pages, 18 figures, accepted to GAFD Special issue on 'Recent Developments in Natural Dynamos'

Small-scale dynamo action is often held responsible for the generation of quiet-Sun magnetic fields. We aim to determine the excitation conditions and saturation level of small-scale dynamos in non-rotating turbulent convection at low magnetic Prandtl numbers. We use high resolution direct numerical simulations of weakly stratified turbulent convection. We find that the critical magnetic Reynolds number for dynamo excitation increases as the magnetic Prandtl number is decreased, which might suggest that small-scale dynamo action is not automatically evident in bodies with small magnetic Prandtl numbers as the Sun. As a function of the magnetic Reynolds number (${\rm Rm}$), the growth rate of the dynamo is consistent with an ${\rm Rm}^{1/2}$ scaling. No evidence for a logarithmic increase of the growth rate with ${\rm Rm}$ is found., Comment: 6 pages, 5 figures, submitted to Astron. Nachr

The influence of fluctuating conductivity on the coefficients known from the mean-field electrodynamics is considered. If the conductivity fluctuations are assumed as uncorrelated with the turbulent velocity field then only the effective magnetic diffusivity of the fluid is reduced and the decay time of a large-scale magnetic field is increased. If the fluctuations of conductivity and flow are correlated in a certain direction then an additional diamagnetic pumping effect results transporting magnetic field in opposite direction to the resistivity flux vector $\langle \eta'\vec{u}'\rangle$. Even for homogeneous turbulence fields in the presence of rotation an alpha effect appears. With the characteristic values of the outer core of the Earth or the solar convection zone, however, the dynamo number of the alpha effect never reaches supercritical values to operate as an $\alpha^2$-dynamo., Comment: 11 pages, 4 figures

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

Previous studies have related surface temperature maps, obtained with the Doppler imaging (DI) technique, of LQ Hya with long-term photometry. We compare surface magnetic field maps, obtained with the Zeeman Doppler imaging (ZDI) technique, with contemporaneous photometry, with the aim of quantifying the star's magnetic cycle characteristics. We inverted Stokes IV spectropolarimetry into magnetic field and surface brightness maps using a tomographic inversion code that models high signal-to-noise ratio mean line profiles produced by the least squares deconvolution (LSD) technique. The magnetic field and surface brightness maps reveal similar patterns to previous DI and ZDI studies: non-axisymmetric polar magnetic field structure, void of fields at mid-latitudes, and a complex structure in the equatorial regions. There is a weak but clear tendency of the polar structures to be linked with strong radial field and the equatorial ones with the azimuthal. We find a polarity reversal in the radial field between 2016 and 2017 coincident with an activity minimum seen in the long-term photometry. The inverted field strengths cannot easily be related with the observed spottedness, but we find that they are partially connected with the retrieved field complexity. Comparing to global magnetoconvection models for rapidly rotating young Suns, this field topology and dominance of the poloidal field component could be explained by a turbulent dynamo, where differential rotation does not play a major role (so called alpha^2 Omega or alpha^2 dynamos), and axi- and non-axisymmetric modes are excited simultaneously. The complex equatorial magnetic field structure could arise from the twisted (helical) wreaths often seen in these simulations, while the polar feature would be connected to the mostly poloidal non-axisymmetric component having a smooth spatial structure., Comment: 13 pages, 8 figures, 7 tables

Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes., Comment: 16 pages, 24 figures, accepted for publication with MNRAS

Context. LQ Hya is one of the most frequently studied young solar analogue stars. Recently, it has been observed to show intriguing behaviour by analysing long-term photometry: during 2003--2009, a coherent spot structure migrating in the rotational frame has been reported by various authors, but since that time the star has entered a chaotic state where coherent structures seem to have disappeared and rapid phase jumps of the photometric minima occur irregularly over time. Aims. LQ Hya is one of the stars included in the SOFIN/FIES long-term monitoring campaign extending over 25 years. Here we publish new temperature maps for the star during 2006--2017, covering the chaotic state of the star. Methods. We use a Doppler imaging technique to derive surface temperature maps from high-resolution spectra. Results. From the mean temperatures of the Doppler maps we see a weak but systematic increase in the surface temperature of the star. This is consistent with the simultaneously increasing photometric magnitude. During nearly all observing seasons we see a high-latitude spot structure which is clearly nonaxisymmetric. The phase behaviour of this structure is very chaotic but agrees reasonably well with the photometry. Equatorial spots are also frequently seen, but many of them we interpret to be artefacts due to the poor to moderate phase coverage. Conclusions. Even during the chaotic phase of the star, the spot topology has remained very similar to the higher activity epochs with more coherent and long-lived spot structures: we see high-latitude and equatorial spot activity, the mid latitude range still being most often void of spots. We interpret the erratic jumps and drifts in phase of the photometric minima to be caused by changes in the high-latitude spot structure rather than the equatorial spots., Comment: 9 pages, 5 figures

The $\Lambda$ effect describes a rotation-induced non-diffusive contribution to the Reynolds stress. It is commonly held responsible for maintaining the observed differential rotation of the Sun and other late-type stars. Here the sensitivity of the $\Lambda$ effect to small-scale magnetic fields and compressibility is studied by means of forced turbulence simulations either with anisotropic forcing in fully periodic cubes or in density-stratified domains with isotropic forcing. Effects of small-scale magnetic fields are studied in cases where the magnetic fields are self-consistently generated by a small-scale dynamo. The results show that small-scale magnetic fields lead to a quenching of the $\Lambda$ effect which is milder than in cases where also a large-scale field is present. The effect of compressibility on the $\Lambda$ effect is negligible in the range of Mach numbers from 0.015 to 0.8. Density stratification induces a marked anisotropy in the turbulence and a vertical $\Lambda$ effect if the forcing scale is roughly two times larger than the density scale height., Comment: 8 pages, 7 figures, published in Astron. Nachr

(abridged) Context: Convective motions overshooting to regions that are formally convectively stable cause extended mixing. Aims: To determine the scaling of overshooting depth ($d_{\rm os}$) at the base of the convection zone as a function of imposed energy flux ($\mathscr{F}_{\rm n}$) and to estimate the extent of overshooting at the base of the solar convection zone. Methods: Three-dimensional Cartesian simulations of compressible non-rotating convection with unstable and stable layers are used. The simulations use either a fixed heat conduction profile or a temperature and density dependent formulation based on Kramers opacity law. The simulations cover a range of almost four orders of magnitude in the imposed flux. Results: A smooth heat conduction profile (either fixed or through Kramers opacity law) leads to a relatively shallow power law with $d_{\rm os}\propto \mathscr{F}_{\rm n}^{0.08}$ for low $\mathscr{F}_{\rm n}$. A fixed step-profile of the heat conductivity at the bottom of the convection zone leads to a somewhat steeper dependency with $d_{\rm os}\propto \mathscr{F}_{\rm n}^{0.12}$. Experiments with and without subgrid-scale entropy diffusion revealed a strong dependence on the effective Prandtl number which is likely to explain the steep power laws as a function of $\mathscr{F}_{\rm n}$ reported in the literature. Furthermore, changing the heat conductivity artificially below the convection zone is shown to lead to substantial underestimation of overshooting depth. Conclusions: Extrapolating from the results obtained with smooth heat conductivity profiles suggest that the overshooting depth for the solar flux is of the order of $0.2$ pressure scale heights at the base of the convection zone which is two to four times higher than estimates from helioseismology., Comment: 15 pages, 13 figures, published in Astron. Astrophys

Context. Studying the spots of late-type stars is crucial for distinguishing between the various proposed dynamo mechanisms believed to be the main cause of starspot activity. For this research it is important to collect observation time series that are long enough to unravel both long- and short-term spot evolution. Doppler imaging is a very efficient method for studying spots of stars that cannot be angularly resolved. Aims. High-resolution spectral observations during 1994-2017 are analysed in order to reveal long- and short-term changes in the spot activity of the FK Comae-type subgiant HD 199178. Methods. Most of the observations were collected with the Nordic Optical Telescope. The Doppler imaging temperature maps were calculated using an inversion technique based on Tikhonov regularisation and utilising multiple spectral lines. Results. We present a unique series of 41 temperature maps spanning more than 23 years. All reliable images show a large cool spot region centred near the visible rotation pole. Some lower latitude cool features are also recovered, although the reliability of these is questionable. There is an expected anti-correlation between the mean surface temperature and the spot coverage. Using the Doppler images, we construct the equivalent of a solar butterfly diagram for HD 199178. Conclusions. HD 199178 clearly has a long-term large and cool spot structure at the rotational pole. This spot structure dominated the spot activity during the years 1994-2017. The size and position of the structure has evolved with time, with a gradual increase during the last years. The lack of lower latitude features prevents the determination of a possible differential rotation., Comment: 11 pages, 7 figures, Astronomy and Astrophysics, parallel version of published article

We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m=1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what is called the inactive branch in observations, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio with increasing rotation rate. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value., Comment: 19 pages, 17 figures, accepted for publication in A&A

Context: Convectively-driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional nonlinear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, it is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity ("grand minima"-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field., Comment: 16 pages, 16 figures, to appear in Astronomy & Astrophysics. The material in the Appendix could form the basis for a possible benchmarking exercise; we would encourage anyone who might be interested in participating in such an exercise to contact the authors

We focus on implementing and optimizing a sixth-order finite-difference solver for simulating compressible fluids on a GPU using third-order Runge-Kutta integration. Since graphics processing units perform well in data-parallel tasks, this makes them an attractive platform for fluid simulation. However, high-order stencil computation is memory-intensive with respect to both main memory and the caches of the GPU. We present two approaches for simulating compressible fluids using 55-point and 19-point stencils. We seek to reduce the requirements for memory bandwidth and cache size in our methods by using cache blocking and decomposing a latency-bound kernel into several bandwidth-bound kernels. Our fastest implementation is bandwidth-bound and integrates $343$ million grid points per second on a Tesla K40t GPU, achieving a $3.6 \times$ speedup over a comparable hydrodynamics solver benchmarked on two Intel Xeon E5-2690v3 processors. Our alternative GPU implementation is latency-bound and achieves the rate of $168$ million updates per second., Comment: 14 pages, 7 figures