Showing 1–50 of 51 results for author: Käpylä, M J

Zeeman-Doppler imaging of five young solar-type stars

Authors: T. Willamo, J. J. Lehtinen, T. Hackman, M. J. Käpylä, O. Kochukhov, S. V. Jeffers, H. Korhonen, S. C. Marsden

Abstract: 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… ▽ More 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. △ Less

Submitted 10 May, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Journal ref: A&A 659, A71 (2022)

Compressible test-field method and its application to shear dynamos

Authors: Maarit J. Käpylä, Matthias Rheinhardt, Axel Brandenburg

Abstract: In this study we present a compressible test-field method (CTFM) for computing $α$ 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… ▽ More In this study we present a compressible test-field method (CTFM) for computing $α$ 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. △ Less

Submitted 3 February, 2022; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 23 pages, 11 figures, submitted to the Astrophysical Journal

Report number: NORDITA-2021-040

Investigating global convective dynamos with mean-field models: full spectrum of turbulent effects required

Authors: Jörn Warnecke, Matthias Rheinhardt, Mariangela Viviani, Frederick Gent, Simo Tuomisto, Maarit J. Käpylä

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: 7 pages with 3 figures, plus 3 pages and 2 figures in the appendix, published in The Astrophysical Journal Letters

Report number: NORDITA-2021-026

Journal ref: 2021 ApJL 919 L13

Inferring magnetic helicity spectrum in spherical domains: the method and example applications

Authors: A. P. Prabhu, N. K. Singh, M. J. Käpylä, A. Lagg

Abstract: 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 $α$ effect. Here we demonstrate a formalism in spherical geometry to infer magnetic helicity spec… ▽ More 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 $α$ 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. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 10 pages, 6 figures, 3 appendices, submitted to Astronomy and Astrophysics

Report number: NORDITA-2021-027

Journal ref: A&A 654, A3 (2021)

Hunting down the cause of solar magnetism

Authors: M. Viviani, A. Prabhu, J. Warnecke, L. Duarte, J. Pekkilä, M. Rheinhardt, M. J. Käpylä

Abstract: 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… ▽ More 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. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Comments: Originally published in "High Performance Computing Science and Engineering, Garching/Munich, 2020". Editors: P. Bastina, D. Kranzmüller, H. Brëuchle, M. Brehm. ISBN: 978-3-9816675-4-7

Interaction of large- and small-scale dynamos in isotropic turbulent flows from GPU-accelerated simulations

Authors: Miikka S. Väisälä, Johannes Pekkilä, Maarit J. Käpylä, Matthias Rheinhardt, Hsien Shang, Ruben Krasnopolsky

Abstract: 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… ▽ More 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 $α$ 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. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 22 pages, 23 figures, 2 tables, Accepted for publication in the Astrophysical Journal

Report number: NORDITA-2020-067

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Authors: A. Barekat, M. J. Käpylä, P. J. Käpylä, E. P. Gilson, H. Ji

Abstract: 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 transpo… ▽ More 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. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Report number: NORDITA-2020-130

Small-Scale Dynamo in Supernova-Driven Interstellar Turbulence

Authors: Frederick A. Gent, Mordecai-Mark Mac Low, Maarit J. Kapyla, Nishant K. Singh

Abstract: Magnetic fields grow quickly even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium of galaxies. Many studies have focused on idealized turbulent driving of the SSD. Here we simulate more realistic supernova-driven turbulence to determine whether it can drive an SSD. Magnetic field growth occurring in our models appears inconsistent with si… ▽ More Magnetic fields grow quickly even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium of galaxies. Many studies have focused on idealized turbulent driving of the SSD. Here we simulate more realistic supernova-driven turbulence to determine whether it can drive an SSD. Magnetic field growth occurring in our models appears inconsistent with simple tangling of magnetic fields, but consistent with SSD action, reproducing and confirming models by Balsara et al. (2004) that did not include physical resistivity $η$. We vary $η$, as well as the numerical resolution and supernova rate, $\dotσ$, to delineate the regime in which an SSD occurs. For a given $\dotσ$ we find convergence for SSD growth rate with resolution of a parsec. For $\dotσ\simeq\dotσ_{\rm sn}$, with $\dotσ_{\rm sn}$ the solar neighbourhood rate, the critical resistivity below which an SSD occurs is $0.005>η_{\rm crit}>0.001\,\rm kpc^{-1}\,\rm km s^{-1}$, and this increases with the supernova rate. Across the modelled range of 0.5--4 pc resolution we find that for $η<η_{\rm crit}$, the SSD saturates at about 5% of kinetic energy equipartition, independent of growth rate. In the range $0.2\dotσ_{\rm sn}\leq \dotσ\leq8\dotσ_{\rm sn}$ growth rate increases with $\dotσ$. SSDs in the supernova-driven interstellar medium commonly exhibit erratic growth. △ Less

Submitted 12 March, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

Authors: A. Brandenburg, A. Johansen, P. A. Bourdin, W. Dobler, W. Lyra, M. Rheinhardt, S. Bingert, N. E. L. Haugen, A. Mee, F. Gent, N. Babkovskaia, C. -C. Yang, T. Heinemann, B. Dintrans, D. Mitra, S. Candelaresi, J. Warnecke, P. J. Käpylä, A. Schreiber, P. Chatterjee, M. J. Käpylä, X. -Y. Li, J. Krüger, J. R. Aarnes, G. R. Sarson , et al. (12 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 7 pages, submitted to the Journal for Open Source Software (JOSS)

Report number: NORDITA-2020-087

Journal ref: Journal of Open Source Software 6, 2807 (2021)

A Knee-Point in the Rotation-Activity Scaling of Late-type Stars with a Connection to Dynamo Transitions

Authors: Jyri J. Lehtinen, Maarit J. Käpylä, Nigul Olspert, Federico Spada

Abstract: 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 th… ▽ More 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. △ Less

Submitted 15 December, 2021; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 10 pages, 5 figures, 2 tables

Report number: NORDITA-2020-052

Journal ref: ApJ 910 110 (2021)

On the existence of shear-current effects in magnetized burgulence

Authors: Maarit J. Käpylä, Javier Álvarez Vizoso, Matthias Rheinhardt, Axel Brandenburg, Nishant K. Singh

Abstract: 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η$. In our setup, a negative sign of its component $η_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficie… ▽ More 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η$. In our setup, a negative sign of its component $η_{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 $η_{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 $η_{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 $α$ 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. △ Less

Submitted 12 January, 2021; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 21 pages, 11 figures, published in Astrophysical Journal

Report number: NORDITA-2020-052

Journal ref: Astrophys. J. 905, 179 (2020)

Physically motivated heat conduction treatment in simulations of solar-like stars: effects on dynamo transitions

Authors: M. Viviani, M. J. Käpylä

Abstract: 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 r… ▽ More 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. △ Less

Submitted 11 June, 2020; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 10 pages, 6 figures, 1 appendix. Submitted to A&A

Report number: NORDITA-2020-057

Journal ref: A&A 645, A141 (2021)

Common dynamo scaling in slowly rotating young and evolved stars

Authors: Jyri J. Lehtinen, Federico Spada, Maarit J. Käpylä, Nigul Olspert, Petri J. Käpylä

Abstract: 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… ▽ More 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. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 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

Shapes of stellar activity cycles

Authors: T. Willamo, T. Hackman, J. J. Lehtinen, M. J. Käpylä, N. Olspert, M. Viviani, J. Warnecke

Abstract: 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 ste… ▽ More 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. △ Less

Submitted 29 April, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: Accepted to A&A; 12 pages, 12 figures

Journal ref: A&A 638, A69 (2020)

Helicity proxies from linear polarisation of solar active regions

Authors: A. Prabhu, A. Brandenburg, M. J. Käpylä, A. Lagg

Abstract: 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 sp… ▽ More 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. △ Less

Submitted 6 December, 2020; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: 10 pages, 7 figures, 2 appendices with 2 figures, published in Astronomy & Astrophysics

Report number: NORDITA-2020-009

Journal ref: A&A 641, A46 (2020)

Rotational dependence of turbulent transport coefficients in global convective dynamo simulations of solar-like stars

Authors: J. Warnecke, M. J. Käpylä

Abstract: 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… ▽ More 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 $α$ tensor increases for moderate rotation rates with Co$^{0.5}$ and levels off for rapid rotation. This behavior agrees with the kinetic $α$, if one considers the decrease of the convective scale with increasing rotation. The $α$ tensor becomes highly anisotropic for Co$\gtrsim 1$. Furthermore, $α_{rr}$ dominates for moderate rotation (1<Co<10), and $α_{φφ}$ for rapid rotation (Co $\gtrsim 10$). The turbulent pumping effect is dominating the meridional transport of the magnetic field. Taking all dynamo effects into account, we find three distinct regimes. For slow rotation, the $α$ and Rädler effects are dominating in the presence of anti-solar differential rotation. For moderate rotation, $α$ and $Ω$ effects are dominant, indicative of $αΩ$ or $α^2Ω$ dynamos in operation, producing equatorward-migrating dynamo waves with a qualitatively solar-like rotation profile. For rapid rotation, an $α^2$ mechanism, with an influence from the Rädler effect appears to be the most probable driver of the dynamo. Our study reveals the presence of a large variety of dynamo effects beyond the classical $αΩ$ mechanism. △ Less

Submitted 22 July, 2020; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: 16 pages with 13 figures, plus 3 pages with 4 figures in the appendix, accepted for publication in A&A

Journal ref: A&A 642, A66 (2020)

Topological changes in the magnetic field of LQ Hya during an activity minimum

Authors: Jyri J. Lehtinen, Maarit J. Käpylä, Thomas Hackman, Oleg Kochukhov, Teemu Willamo, Stephen C. Marsden, Sandra V. Jeffers, Gregory W. Henry, Lauri Jetsu

Abstract: 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… ▽ More 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. △ Less

Submitted 15 December, 2021; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: 13 pages, 8 figures, 7 tables

Journal ref: A&A 660, A141 (2022)

Spot evolution on LQ Hya during 2006--2017: temperature maps based on SOFIN and FIES data

Authors: Elizabeth M. Cole-Kodikara, Maarit J. Käpylä, Jyri J. Lehtinen, Thomas Hackman, Ilya V. Ilyin, Nikolai Piskunov, Oleg Kochukhov

Abstract: 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 disappeare… ▽ More 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. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 9 pages, 5 figures

Journal ref: A&A 629, A120 (2019)

Stellar dynamos in the transition regime: multiple dynamo modes and anti-solar differential rotation

Authors: M. Viviani, M. J. Käpylä, J. Warnecke, P. J. Käpylä, M. Rheinhardt

Abstract: 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… ▽ More 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 $αΩ$ 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 $α$ effect is comparable or even larger than the $Ω$ effect in generating the toroidal magnetic field, and therefore, the dynamo seems to be $α^2Ω$ or $α^2$ type. We further find that the effective large-scale flows are significantly altered by turbulent pumping. △ Less

Submitted 22 November, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 13 pages, 7 figures. Published in ApJ

Journal ref: The Astrophysical Journal, Volume 886, Issue 1, article id. 21, 10 pp. (2019)

Turbulent viscosity and effective magnetic Prandtl number from simulations of isotropically forced turbulence

Authors: Petri J. Käpylä, Matthias Rheinhardt, Axel Brandenburg, Maarit J. Käpylä

Abstract: (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 turbule… ▽ More (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. △ Less

Submitted 17 February, 2020; v1 submitted 3 January, 2019; originally announced January 2019.

Comments: 13 pages, 12 figures, accepted in Astron. Astrophys

Report number: NORDITA-2018-126

Journal ref: A&A 636, A93 (2020)

Starspot activity of HD 199178. Doppler images from 1994--2017

Authors: Thomas Hackman, Iya Ilyin, Jyri J. Lehtinen, Oleg Kochukhov, Maarit J. Käpylä, Nikolai Piskunov, Teemu Willamo

Abstract: 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 canno… ▽ More 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. △ Less

Submitted 7 June, 2019; v1 submitted 5 December, 2018; originally announced December 2018.

Comments: 11 pages, 7 figures, Astronomy and Astrophysics, parallel version of published article

Journal ref: A&A 625, A79 (2019)

Long-term spot monitoring of the young solar analogue V889 Her

Authors: T. Willamo, T. Hackman, J. J. Lehtinen, M. J. Käpylä, I. Ilyin, G. W. Henry, L. Jetsu, O. Kochukhov, N. Piskunov

Abstract: Starspots are important manifestations of stellar magnetic activity. By studying their behaviour in young solar analogues, we can unravel the properties of their magnetic cycles. This gives crucial information of the underlying dynamo process. Comparisons with the solar cycle enable us to infer knowledge about how the solar dynamo has evolved during the Sun's lifetime. Here we study the correlatio… ▽ More Starspots are important manifestations of stellar magnetic activity. By studying their behaviour in young solar analogues, we can unravel the properties of their magnetic cycles. This gives crucial information of the underlying dynamo process. Comparisons with the solar cycle enable us to infer knowledge about how the solar dynamo has evolved during the Sun's lifetime. Here we study the correlation between photometric brightness variations, spottedness and mean temperature in V889 Her, a young solar analogue. Our data covers 18 years of spectroscopic and 25 years of photometric observations. We use Doppler imaging to derive temperature maps from high-resolution spectra. We use the Continuous Period Search method to retrieve mean V-magnitudes from photometric data. Our Doppler imaging maps show a persistent polar spot structure varying in strength. This structure is centered slightly off the rotational pole. The mean temperature derived from the maps shows an overall decreasing trend, as does the photometric mean brightness, until it reaches its minimum around 2017. The filling factor of cool spots, however, shows only a weak tendency to anti-correlate with the decreasing mean brightness. We interpret V889 Her to have entered into a grand maximum in its activity. The clear relation between the mean temperature of the Doppler imaging surface maps and the mean magnitude supports the reliability of the Doppler images. The lack of correlation between the mean magnitude and the spottedness may indicate that bright features in the Doppler images are real. △ Less

Submitted 27 February, 2019; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: 9 pages, 7 figures, 4 tables

Journal ref: A&A 622, A170 (2019)

$f$-mode strengthening from a localized bipolar subsurface magnetic field

Authors: Nishant K. Singh, Harsha Raichur, Maarit J. Käpylä, Matthias Rheinhardt, Axel Brandenburg, Petri J. Käpylä

Abstract: 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. So… ▽ More 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ω$ 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. △ Less

Submitted 12 April, 2019; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: 15 pages, 7 figures, revised version submitted to GAFD special issue on "Physics and Algorithms of the Pencil Code"

Report number: NORDITA-2018-084

Journal ref: Geophys. Astro. Fluid 2020, 114, 196-212 (2020)

Sensitivity to luminosity, centrifugal force, and boundary conditions in spherical shell convection

Authors: Petri J. Käpylä, Frederick A. Gent, Nigul Olspert, Maarit J. Käpylä, Axel Brandenburg

Abstract: 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… ▽ More 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. △ Less

Submitted 2 January, 2019; v1 submitted 24 July, 2018; originally announced July 2018.

Comments: 26 pages, 14 figures, accepted to GAFD Special issue on 'Physics and Algorithms of the Pencil Code'

Report number: NORDITA-2018-061

Journal ref: Geophys. Astrophys. Fluid Dyn. 114, 8-34 (2020)

Modelling supernova driven turbulence

Authors: Frederick A. Gent, Mordecai-Mark Mac Low, Maarit J. Käpylä, Graeme R. Sarson, James F. Hollins

Abstract: High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to h… ▽ More High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and a novel generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the one-dimensional (1D) Riemann shock tube (Sod, J. Comput. Phys. 1978, 27), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the three-dimensional (3D) numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and compare these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi, McKee and Bertschinger (Astrophys. J. 1988, 334) results incorporating cooling and heating processes. We show that our new timestep algorithm leads to linear rather than quadratic resolution dependence as the strength of the artificial viscosity varies, because of the corresponding change in the strength of interzone gradients. △ Less

Submitted 18 June, 2019; v1 submitted 5 June, 2018; originally announced June 2018.

Comments: 25 pages, 20 figures, Pencil Code special edition

MSC Class: 35; 85; 68

Journal ref: GAFD, 114, 77-105, 2020

Bihelical spectrum of solar magnetic helicity and its evolution

Authors: Nishant K. Singh, Maarit J. Käpylä, Axel Brandenburg, Petri J. Käpylä, Andreas Lagg, Ilpo Virtanen

Abstract: 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 Carringt… ▽ More 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. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: 10 pages, 10 figures, submitted to the Astrophysical Journal

Report number: NORDITA-2018-031

Journal ref: Astrophys. J. 863, 182 (2018)

Effects of a subadiabatic layer on convection and dynamos in spherical wedge simulations

Authors: Petri J. Käpylä, Mariangela Viviani, Maarit J. Käpylä, Axel Brandenburg, Federico Spada

Abstract: 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 fi… ▽ More 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. △ Less

Submitted 2 January, 2019; v1 submitted 15 March, 2018; originally announced March 2018.

Comments: 32 pages, 18 figures, accepted to GAFD Special issue on 'Recent Developments in Natural Dynamos'

Report number: NORDITA 2018-022

Journal ref: Geophys. Astrophys. Fluid Dyn. 113, 149-183 (2020)

Small-scale dynamos in simulations of stratified turbulent convection

Authors: Petri J. Käpylä, Maarit J. Käpylä, Axel Brandenburg

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: 6 pages, 5 figures, submitted to Astron. Nachr

Report number: NORDITA-2018-011

Journal ref: Astron. Nachr. 339, 127-133 (2018)

Estimating activity cycles with probabilistic methods II. The Mount Wilson Ca H&K data

Authors: N. Olspert, J. Lehtinen, M. J. Käpylä, J. Pelt, A. Grigorievskiy

Abstract: Debate over the existence of branches in the stellar activity-rotation diagrams continues. Application of modern time series analysis tools to study the mean cycle periods in chromospheric activity index is lacking. We develop such models, based on Gaussian processes, for one-dimensional time series and apply it to the extended Mount Wilson Ca H&K sample. Our main aim is to study how the previousl… ▽ More Debate over the existence of branches in the stellar activity-rotation diagrams continues. Application of modern time series analysis tools to study the mean cycle periods in chromospheric activity index is lacking. We develop such models, based on Gaussian processes, for one-dimensional time series and apply it to the extended Mount Wilson Ca H&K sample. Our main aim is to study how the previously commonly used assumption of strict harmonicity of the stellar cycles as well as handling of the linear trends affects the results. We introduce three methods of different complexity, starting with the simple Bayesian harmonic model and followed by Gaussian Process models with periodic and quasi-periodic covariance functions. We confirm the existence of two populations in the activity-period diagram. We find only one significant trend in the inactive population, namely that the cycle periods get shorter with increasing rotation. This is in contrast with earlier studies, that postulate the existence of trends in both of the populations. In terms of rotation to cycle period ratio, our data is consistent with only two activity branches such that the active branch merges together with the transitional one. The retrieved stellar cycles are uniformly distributed over the R'HK activity index, indicating that the operation of stellar large-scale dynamos carries smoothly over the Vaughan-Preston gap. At around the solar activity index, however, indications of a disruption in the cyclic dynamo action are seen. Our study shows that stellar cycle estimates depend significantly on the model applied. Such model-dependent aspects include the improper treatment of linear trends, while the assumption of strict harmonicity can result in the appearance of double cyclicities that seem more likely to be explained by the quasi-periodicity of the cycles. △ Less

Submitted 12 July, 2018; v1 submitted 21 December, 2017; originally announced December 2017.

Estimating activity cycles with probabilistic methods I. Bayesian Generalised Lomb-Scargle Periodogram with Trend

Authors: N. Olspert, J. Pelt, M. J. Käpylä, J. Lehtinen

Abstract: Period estimation is one of the central topics in astronomical time series analysis, where data is often unevenly sampled. Especially challenging are studies of stellar magnetic cycles, as there the periods looked for are of the order of the same length than the datasets themselves. The datasets often contain trends, the origin of which is either a real long-term cycle or an instrumental effect, b… ▽ More Period estimation is one of the central topics in astronomical time series analysis, where data is often unevenly sampled. Especially challenging are studies of stellar magnetic cycles, as there the periods looked for are of the order of the same length than the datasets themselves. The datasets often contain trends, the origin of which is either a real long-term cycle or an instrumental effect, but these effects cannot be reliably separated, while they can lead to erroneous period determinations if not properly handled. In this study we aim at developing a method that can handle the trends properly, and by performing extensive set of testing, we show that this is the optimal procedure when contrasted with methods that do not include the trend directly to the model. The effect of the form of the noise (whether constant or heteroscedastic) on the results is also investigated. We introduce a Bayesian Generalised Lomb-Scargle Periodogram with Trend (BGLST), which is a probabilistic linear regression model using Gaussian priors for the coefficients and uniform prior for the frequency parameter. We show, using synthetic data, that when there is no prior information on whether and to what extent the true model of the data contains a linear trend, the introduced BGLST method is preferable to the methods which either detrend the data or leave the data untrended before fitting the periodic model. Whether to use noise with different than constant variance in the model depends on the density of the data sampling as well as on the true noise type of the process. △ Less

Submitted 13 March, 2018; v1 submitted 21 December, 2017; originally announced December 2017.

Journal ref: A&A 615, A111 (2018)

Transition from axi- to nonaxisymmetric dynamo modes in spherical convection models of solar-like stars

Authors: Mariangela Viviani, Jörn Warnecke, Maarit J. Käpylä, Petri J. Käpylä, Nigul Olspert, Elizabeth M. Cole-Kodikara, Jyri J. Lehtinen, Axel Brandenburg

Abstract: 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… ▽ More 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. △ Less

Submitted 7 May, 2018; v1 submitted 27 October, 2017; originally announced October 2017.

Comments: 19 pages, 17 figures, accepted for publication in A&A

Report number: NORDITA-2017-113

Journal ref: A&A 616, A160 (2018)

Large-scale dynamos in rapidly rotating plane layer convection

Authors: P. J. Bushby, P. J. Käpylä, Y. Masada, A. Brandenburg, B. Favier, C. Guervilly, M. J. Käpylä

Abstract: 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 up… ▽ More 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. △ Less

Submitted 15 January, 2018; v1 submitted 9 October, 2017; originally announced October 2017.

Comments: 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

Report number: NORDITA-2017-072

Journal ref: A&A 612, A97 (2018)

Long-term variations of turbulent transport coefficients in a solar-like convective dynamo simulation

Authors: Frederick A. Gent, Maarit J. Käpylä, Jörn Warnecke

Abstract: The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dip… ▽ More The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dipolar and quadrupolar modes. We compute the turbulent transport coefficients, describing the effects of the turbulent velocity field on the mean magnetic field, using the test-field method. The test-field analysis provides a plausible explanation of the missing cycle in terms of the reduction of $α_{φφ}$ in advance of the reduced surface activity, and enhanced downward turbulent pumping during the event to confine some of the magnetic field at the bottom of the convection zone, where local maximum of magnetic energy is observed during the event. At the same time, however, a quenching of the turbulent magnetic diffusivities is observed, albeit differently distributed in depth compared to the other transport coefficients. Therefore, dedicated mean-field modelling is required for verification. △ Less

Submitted 6 October, 2017; v1 submitted 1 September, 2017; originally announced September 2017.

Comments: 11 pages, 12 figures, accepted by AN for 14th Potsdam Thinkshop

Methods for compressible fluid simulation on GPUs using high-order finite differences

Authors: Johannes Pekkilä, Miikka S. Väisälä, Maarit J. Käpylä, Petri J. Käpylä, Omer Anjum

Abstract: 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 cac… ▽ More 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. △ Less

Submitted 27 July, 2017; originally announced July 2017.

Comments: 14 pages, 7 figures

Journal ref: Computer Physics Communications, Volume 217, August 2017, Pages 11-22

The supernova-regulated ISM. III. Generation of vorticity, helicity and mean flows

Authors: Maarit J. Käpylä, Frederick A. Gent, Miikka S. Väisälä, Graeme R. Sarson

Abstract: The forcing of interstellar turbulence, driven mainly by supernova explosions, is irrotational in nature, but the development of significant amounts of vorticity and helicity, accompanied by large-scale dynamo action, has been reported. Several earlier investigations examined vorticity production in simpler systems; here all the relevant processes can be considered simultaneously. We also investig… ▽ More The forcing of interstellar turbulence, driven mainly by supernova explosions, is irrotational in nature, but the development of significant amounts of vorticity and helicity, accompanied by large-scale dynamo action, has been reported. Several earlier investigations examined vorticity production in simpler systems; here all the relevant processes can be considered simultaneously. We also investigate the mechanisms for the generation of net helicity and large-scale flow in the system. We use a three-dimensional, stratified, rotating and shearing local simulation domain of the size 1x1x2 kpc$^3$, forced with SN explosions occurring at the rate typical of the solar neighbourhood in the Milky Way. In addition to the nominal simulation run with realistic Milky Way parameters, we vary the rotation and shear rates, but keep the absolute value of their ratio fixed. Reversing the sign of shear vs. rotation allows us to separate the rotation- and shear-generated contributions. As in earlier studies, we find the generation of significant amounts of vorticity, with on average 65% of the kinetic energy being in the rotational modes. The vorticity production can be related to the baroclinicity of the flow, especially in the regions of hot, dilute clustered supernova bubbles. In these regions, the vortex stretching acts as a sink of vorticity. The net helicities produced by rotation and shear are of opposite signs for physically motivated rotation laws, with the solar neighbourhood parameters resulting in the near cancellation of the total net helicity. We also find the excitation of oscillatory mean flows, the strength and oscillation period of which depend on the Coriolis and shear parameters; we interpret these as signatures of the anisotropic kinetic (AKA) effect. We use the method of moments to fit for the turbulent transport coeffcients, and find $α_{\rm AKA}$ values of the order 3-5 km/s. △ Less

Submitted 19 October, 2017; v1 submitted 24 May, 2017; originally announced May 2017.

Comments: 19 pages, 12 figures, Astronomy and Astrophysics, accepted for publication

Journal ref: A&A 611, A15 (2018)

The supernova-regulated ISM. IV. A comparison of simulated polarization with Planck observations

Authors: Miikka S. Väisälä, Frederick A. Gent, Mika Juvela, Maarit J. Käpylä

Abstract: The efforts for comparing dust polarization measurements with synthetic observations from MHD models have previously concentrated on the scale of molecular clouds. Here we extend the model comparisons to kiloparsec scales, taking into account hot shocked gas generated by supernovae, and a non-uniform dynamo-generated magnetic field at both large and small scales down to 4 pc spatial resolution. Ra… ▽ More The efforts for comparing dust polarization measurements with synthetic observations from MHD models have previously concentrated on the scale of molecular clouds. Here we extend the model comparisons to kiloparsec scales, taking into account hot shocked gas generated by supernovae, and a non-uniform dynamo-generated magnetic field at both large and small scales down to 4 pc spatial resolution. Radiative transfer calculations are used to model dust emission and polarization on the top of MHD simulations. We compute synthetic maps of column density $N_H$, polarization fraction $p$, and polarization angle dispersion $S$, and study their dependencies on the important properties of the MHD simulations. These include the large-scale magnetic field and its orientation, the small-scale magnetic field, and supernova-driven shocks. Similar filament-like structures of $S$ as seen in the Planck all-sky maps are visible in our synthetic results, although the smallest scale structures are absent from our maps. Supernova-driven shock fronts and $S$ do not show significant correlation. Instead, $S$ can clearly be attributed to the distribution of the small-scale magnetic field. We also find that the large-scale magnetic field influences the polarization properties, such that, for a given strength of magnetic fluctuation, a strong plane-of-the-sky mean field weakens the observed $S$, while strengthening $p$. The anticorrelation of $p$ and $S$, and decreasing $p$ as a function of $N_H$ are consistent across all synthetic observations. The magnetic fluctuations follow an exponential distribution, rather than Gaussian, characteristic of flows with intermittent repetitive shocks. The observed polarization properties and column densities are sensitive to the line of sight distance over which the emission is integrated. △ Less

Submitted 13 June, 2018; v1 submitted 21 March, 2017; originally announced March 2017.

Comments: 21 pages, 18 figures, 2 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 614, A101 (2018)

Extended subadiabatic layer in simulations of overshooting convection

Authors: Petri J. Käpylä, Matthias Rheinhardt, Axel Brandenburg, Rainer Arlt, Maarit J. Käpylä, Andreas Lagg, Nigul Olspert, Jörn Warnecke

Abstract: We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kr… ▽ More We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ. △ Less

Submitted 23 August, 2017; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: 6 pages, 3 figures, published in ApJL

Report number: NORDITA-2017-25

Journal ref: Astrophys. J. Lett. 845, L23 (2017)

Convection-driven spherical shell dynamos at varying Prandtl numbers

Authors: P. J. Käpylä, M. J. Käpylä, N. Olspert, J. Warnecke, A. Brandenburg

Abstract: (abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at v… ▽ More (abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at various thermal and magnetic Prandtl numbers (from 0.25 to 2 and 5, respectively). Results: We find that the rotation profiles for high thermal diffusivity show a monotonically increasing angular velocity from the bottom of the convection zone to the top and from the poles toward the equator. For sufficiently rapid rotation, a region of negative radial shear develops at mid-latitudes as the thermal diffusivity is decreased. This coincides with a change in the dynamo mode from poleward propagating activity belts to equatorward propagating ones. Furthermore, the cyclic solutions disappear at the highest magnetic Reynolds numbers. The total magnetic energy increases with the magnetic Reynolds number in the range studied here ($5-151$), but the energies of the mean magnetic fields level off at high magnetic Reynolds numbers. The differential rotation is strongly affected by the magnetic fields and almost vanishes at the highest magnetic Reynolds numbers. In some of our most turbulent cases we find that two regimes are possible where either differential rotation is strong and mean magnetic fields relatively weak or vice versa. Conclusions: Our simulations indicate a strong non-linear feedback of magnetic fields on differential rotation, leading to qualitative changes in the behaviors of large-scale dynamos at high magnetic Reynolds numbers. Furthermore, we do not find indications of the simulations approaching an asymptotic regime where the results would be independent of diffusion coefficients. △ Less

Submitted 12 September, 2016; v1 submitted 19 May, 2016; originally announced May 2016.

Comments: 18 pages, 10 figures, submitted to Astron. Astrophys

Report number: NORDITA-2016-38

Journal ref: A&A 599, A4 (2017)

Robustness of oscillatory $α^2$ dynamos in spherical wedges

Authors: E. Cole, A. Brandenburg, P. J. Käpylä, M. J. Käpylä

Abstract: We study the connection between spherical wedge and full spherical shell geometries using simple mean-field $α^2$ dynamos. We solve the equations for a one-dimensional time-dependent mean-field dynamo to examine the effects of varying the polar angle $θ_0$ between the latitudinal boundaries and the poles in spherical coordinates. We investigate the effects of turbulent magnetic diffusivity and… ▽ More We study the connection between spherical wedge and full spherical shell geometries using simple mean-field $α^2$ dynamos. We solve the equations for a one-dimensional time-dependent mean-field dynamo to examine the effects of varying the polar angle $θ_0$ between the latitudinal boundaries and the poles in spherical coordinates. We investigate the effects of turbulent magnetic diffusivity and $α$ effect profiles as well as different latitudinal boundary conditions to isolate parameter regimes where oscillatory solutions are found. Finally, we add shear along with a damping term mimicking radial gradients to study the resulting dynamo regimes. We find that the commonly used perfect conductor boundary condition leads to oscillatory $α^2$ dynamo solutions only if the wedge boundary is at least one degree away from the poles. Other boundary conditions always produce stationary solutions. By varying the profile of the turbulent magnetic diffusivity alone, oscillatory solutions are achieved with models extending to the poles, but the magnetic field is strongly concentrated near the poles and the oscillation period is very long. By introducing radial shear and a damping term mimicking radial gradients, we again see oscillatory dynamos, and the direction of drift follows the Parker--Yoshimura rule. Oscillatory solutions in the weak shear regime are found only in the wedge case with $θ_0 = 1^\circ$ and perfect conductor boundaries. A reduced $α$ effect near the poles with a turbulent diffusivity concentrated toward the equator yields oscillatory dynamos with equatorward migration and reproduces best the solutions in spherical wedges. △ Less

Submitted 20 January, 2016; originally announced January 2016.

Comments: 9 pages, 10 figures, submitted to A&A

Report number: NORDITA-2016-4

Journal ref: A&A 593, A134 (2016)

Turbulent transport coefficients in spherical wedge dynamo simulations of solar-like stars

Authors: Jörn Warnecke, Matthias Rheinhardt, S. Tuomisto, Petri J. Käpylä, Maarit J. Käpylä, Axel Brandenburg

Abstract: We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($φ$) averaged fields to determine t… ▽ More We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($φ$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $α$ tensor. This method has previously been used either in simulations in Cartesian coordinates or in the geodynamo context and is applied here for the first time to fully compressible simulations of solar-like dynamos. We find that the $φφ$-component of the $α$ tensor does not follow the profile expected from that of kinetic helicity. The turbulent pumping velocities significantly alter the effective mean flows acting on the magnetic field and therefore challenge the flux transport dynamo concept. All coefficients are significantly affected by dynamically important magnetic fields. Quenching as well as enhancement are being observed. This leads to a modulation of the coefficients with the activity cycle. The temporal variations are found to be comparable to the time-averaged values and seem to be responsible for a nonlinear feedback on the magnetic field generation. Furthermore, we quantify the validity of the Parker-Yoshimura rule for the equatorward propagation of the mean magnetic field in the present case. △ Less

Submitted 16 January, 2018; v1 submitted 14 January, 2016; originally announced January 2016.

Comments: 21 pages, 21 figures, published

Report number: NORDITA-2016-3

Journal ref: A&A 609, A51 (2018)

Magnetic flux concentrations from turbulent stratified convection

Authors: P. J. Käpylä, A. Brandenburg, N. Kleeorin, M. J. Käpylä, I. Rogachevskii

Abstract: (abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal unifo… ▽ More (abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation. △ Less

Submitted 11 November, 2015; originally announced November 2015.

Comments: 13 pages, 13 figures, submitted to Astronomy & Astrophysics, associated animations can be found at http://research.ics.aalto.fi/cmdaa/group-Movies.shtml

Report number: NORDITA-2015-129

Journal ref: A&A 588, A150 (2016)

Zeeman-Doppler imaging of active young solar type stars

Authors: Thomas Hackman, Jyri Lehtinen, Lisa Rosén, Oleg Kochukhov, Maarit J. Käpylä

Abstract: By studying young magnetically active late-type stars, i.e. analogues to the young Sun, one can draw conclusions on the evolution of the solar dynamo. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric spots in three young late-type stars. High-resolution spectropolarimetry of the targets were obtained with the HARPSpol i… ▽ More By studying young magnetically active late-type stars, i.e. analogues to the young Sun, one can draw conclusions on the evolution of the solar dynamo. We determine the topology of the surface magnetic field and study the relation between the magnetic field and cool photospheric spots in three young late-type stars. High-resolution spectropolarimetry of the targets were obtained with the HARPSpol instrument mounted at the ESO 3.6 m telescope. The signal-to-noise ratio of the Stokes IV measurements were boosted by combining the signal from a large number of spectroscopic absorption lines through the least squares deconvolution technique. Surface brightness and magnetic field maps were calculated using the Zeeman-Doppler imaging technique. All the three targets show clear signs of both magnetic fields and cool spots. Only one of the targets, namely V1358 Ori, shows evidence of the dominance of non-axisymmetric modes. In two of the targets, the poloidal field is significantly stronger than the toroidal one, indicative of an $α^2$-type of a dynamo, in which convective turbulence effects dominate over the weak differential rotation. In two of the cases there is a slight anti-correlation between the cool spots and the strength of the radial magnetic field. However, even in these cases the correlation is much weaker than in the case of sunspots. The weak correlation between the measured radial magnetic field and cool spots may indicate a more complex magnetic field structure in the spots or spot groups involving mixed magnetic polarities. Comparison with a previously published magnetic field map shows that on one of the stars, HD 29615, the underlying magnetic field has changed its polarity between 2009 and 2013. △ Less

Submitted 1 March, 2016; v1 submitted 8 September, 2015; originally announced September 2015.

Comments: 7 pages, 6 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 587, A28 (2016)

Multiple dynamo modes as a mechanism for long-term solar activity variations

Authors: Maarit J. Käpylä, Petri J. Käpylä, Nigul Olspert, Axel Brandenburg, Jörn Warnecke, Bidya B. Karak, Jaan Pelt

Abstract: Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years,… ▽ More Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years, during which epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. The DNS employed is a semi-global (wedge) magnetoconvection model. For data analysis we use Ensemble Empirical Mode Decomposition (EEMD) and phase dispersion ($D^2$) methods. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like. This mode is accompanied by a higher frequency mode near the surface and a low-frequency mode in the bottom of the convection zone. The overall behavior of the dynamo solution is very complex exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum. We interpret the overall irregular behavior to be due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the differential rotation profile via the suppression of the Reynolds stresses. △ Less

Submitted 22 January, 2016; v1 submitted 20 July, 2015; originally announced July 2015.

Comments: 24 pages, 16 figures, Astronomy and Astrophysics, accepted for publication

Report number: NORDITA-2015-90

Journal ref: A&A 589, A56 (2016)

Doppler Imaging of LQ Hya for 1998-2002

Authors: E. M. Cole, T. Hackman, M. J. Käpylä, I. Ilyin, O. Kochukhov, N. Piskunov

Abstract: We study the spot distribution on the surface of LQ~Hya during the observing seasons October 1998 -- November 2002. We look for persistent active longitudes, trends in the level of spot activity and compare to photometric data, specifically to the derived time epochs of the lightcurve minima. We apply the Doppler Imaging technique on photospheric spectral lines using an inversion code to retrieve… ▽ More We study the spot distribution on the surface of LQ~Hya during the observing seasons October 1998 -- November 2002. We look for persistent active longitudes, trends in the level of spot activity and compare to photometric data, specifically to the derived time epochs of the lightcurve minima. We apply the Doppler Imaging technique on photospheric spectral lines using an inversion code to retrieve images of the surface temperature. We present new temperature maps using multiple spectral lines for a total of 7 seasons. We calculate spot coverage fractions from each map, and as a result we find a general trend that is interpreted to be an indication of a spot cycle. There is a minimum during the observing season of March 1999. After this the activity increases until November 2000, followed by a general decrease in activity again. We find no evidence for active longitudes persisting over multiple observing seasons. The spot activity appears to be concentrated to two latitude regions. The high latitude spots are particularly strong when the spot coverage is at a maximum. Using the currently accepted rotation period, we find spot structures to show a trend in the phase-time plot, indicative of a need for a longer period. We conclude that the long-term activity of LQ~Hya is more chaotic than that of some magnetically active binary stars analyzed with similar methods, but still with clear indications of an activity cycle. △ Less

Submitted 14 April, 2015; originally announced April 2015.

Comments: 9 pages, 10 figures

Journal ref: A&A 581, A69 (2015)

Influence of a coronal envelope as a free boundary to global convective dynamo simulations

Authors: Jörn Warnecke, Petri J. Käpylä, Maarit J. Käpylä, Axel Brandenburg

Abstract: We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above… ▽ More We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above. The interface represents a free surface. We compare our model to models that have no coronal envelope. The presence of a coronal envelope is found to modify the Reynolds stress and the $Λ$ effect resulting in a weaker and non-cylindrical differential rotation. This is related to the reduced latitudinal temperature variations that are caused by and dependent on the angular velocity. Some simulations develop a near-surface shear layer that we can relate to a sign change in the meridional Reynolds stress term in the thermal wind balance equation. Furthermore, the presence of a free surface changes the magnetic field evolution since the toroidal field is concentrated closer to the surface. In all simulations, however, the migration direction of the mean magnetic field can be explained by the Parker--Yoshimura rule, which is consistent with earlier findings. A realistic treatment of the upper boundary in spherical dynamo simulations is crucial for the dynamics of the flow and magnetic field evolution. △ Less

Submitted 19 December, 2016; v1 submitted 17 March, 2015; originally announced March 2015.

Comments: 16 pages, 13 figures, Astronomy and Astrophysics, published

Report number: NORDITA-2015-020

Journal ref: A&A 596, A115 (2016)

Multiperiodicity, modulations and flip-flops in variable star light curves III. Carrier fit analysis of LQ Hya photometry for 1982-2014

Authors: Nigul Olspert, Maarit J. Käpylä, Jaan Pelt, Elizabeth M. Cole, Thomas Hackman, Jyri Lehtinen, Gregory W. Henry

Abstract: We study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the perio… ▽ More We study LQ Hya photometry for 1982-2014 with the carrier fit (CF) -method and compare our results to earlier photometric analysis and recent Doppler imaging maps. We first utilize different types of statistical methods to estimate various candidates for the carrier period for the CF method. Secondly, a global fit to the whole data set and local fits to shorter segments are computed with the period that is found to be the optimal one. The harmonic least-squares analysis of all the available data reveals a short period close to 1.6 days as a limiting value for a set of significant frequencies. We interpret this as the rotation period of the spots near the equatorial region. In addition, the distribution of the significant periods is found to be bimodal, hinting of a longer-term modulating period, which we set out to study with a two-harmonic CF model. Weak modulation signal is, indeed retrieved, with a period of roughly 6.9 years. The phase dispersion analysis gives a clear symmetric minimum for coherence times lower than and around 100 days. We interpret this as the mean rotation period of the spots (1.60514 days), and this value is chosen to be used as the carrier period for the CF analysis. With the CF method we seek for any systematic trends in the spot distribution in the global time frame, and locally look for abrupt phase changes earlier reported in rapidly rotating objects. During 2005-2008 the global CF reveals a coherent structure rotating with a period of 1.6037 days, while during most other times the spot distribution appears rather random in phase. The evolution of the spot distribution of the object is found to be very chaotic, with no clear signs of an azimuthal dynamo wave that would persist over longer time scales, although the short-lived coherent structures observed occasionally do not rotate with the same speed as the mean spot distribution. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: 14 pages, 10 figures, submitted to Astronomy and Astrophysics

Journal ref: A&A 577, A120 (2015)

On the cause of solar-like equatorward migration in global convective dynamo simulations

Authors: Jörn Warnecke, Petri J. Käpylä, Maarit J. Käpylä, Axel Brandenburg

Abstract: We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can be explained by an $αΩ$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due… ▽ More We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can be explained by an $αΩ$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $α$ effect in the northern (southern) hemisphere and a negative radial gradient of $Ω$ outside the inner tangent cylinder of these models. This idea is supported by a strong correlation between negative radial shear and toroidal field strength in the region of equatorward propagation. △ Less

Submitted 28 November, 2014; v1 submitted 10 September, 2014; originally announced September 2014.

Comments: 6 pages, 5 figures, published in ApJ Letters

Report number: NORDITA-2014-103

Journal ref: Astrophys. J. Lett. 796, L12 (2014)

Magnetically controlled stellar differential rotation near the transition from solar to anti-solar profiles

Authors: Bidya Binay Karak, Petri J. Kapyla, Maarit J. Kapyla, Axel Brandenburg, N. Olspert, J. Pelt

Abstract: Late-type stars rotate differentially owing to anisotropic turbulence in their outer convection zones. The rotation is called solar-like (SL) when the equator rotates fastest and anti-solar (AS) otherwise. Hydrodynamic simulations show a transition from SL to AS rotation as the influence of rotation on convection is reduced, but the opposite transition occurs at a different point in the parameter… ▽ More Late-type stars rotate differentially owing to anisotropic turbulence in their outer convection zones. The rotation is called solar-like (SL) when the equator rotates fastest and anti-solar (AS) otherwise. Hydrodynamic simulations show a transition from SL to AS rotation as the influence of rotation on convection is reduced, but the opposite transition occurs at a different point in the parameter space. The system is bistable, i.e., SL and AS rotation profiles can both be stable. We study the effect of a dynamo-generated magnetic field on the large-scale flows, particularly on the possibility of bistable behavior of differential rotation. We solve the hydromagnetic equations numerically in a rotating spherical shell for a set of different radiative conductivities controlling the relative importance of convection. In agreement with earlier findings, our models display SL rotation profiles when the rotational influence on convection is strong and a transition to AS when the rotational influence decreases. We find that dynamo-generated magnetic fields help to produce SL differential rotation compared to the hydrodynamic simulations. We do not observe any bistable states of differential rotation. In the AS cases we get coherent single-cell meridional circulation, whereas in SL cases we get multi-cellular patterns. In both cases, we obtain poleward circulation near the surface with a magnitude close to that observed in the Sun. Moreover, both differential rotation and meridional circulation have significant magnetic cycle-related variations that are similar in strength to those of the Sun. Purely hydrodynamic simulations of differential rotation and meridional circulation are shown to be of limited relevance as magnetic fields, self-consistently generated by dynamo action, significantly affect the flows. △ Less

Submitted 21 January, 2015; v1 submitted 3 July, 2014; originally announced July 2014.

Comments: Accepted for publication in A&A

Report number: NORDITA-2014-84

Journal ref: A&A 576, A26 (2015)

Quenching and anisotropy of hydromagnetic turbulent transport

Authors: Bidya Binay Karak, Matthias Rheinhardt, Axel Brandenburg, Petri J. Kapyla, Maarit J. Kapyla

Abstract: Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $α$ effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Ro… ▽ More Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $α$ effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Roberts flow, isotropically forced non-helical turbulence, and rotating thermal convection. We see significant quenching only when the mean magnetic field is larger than the equipartition value of the turbulence. Expressing the magnetic field in terms of the equipartition value of the {\it quenched} flows, we obtain for the quenching exponents of the turbulent magnetic diffusivity about 1.3, 1.1, and 1.3 for Roberts flow, forced turbulence, and convection, respectively. However, when the magnetic field is expressed in terms of the equipartition value of the unquenched flows these quenching exponents become about 4, 1.5, and 2.3, respectively. For the $α$ effect, the exponent is about 1.3 for the Roberts flow and 2 for convection in the first case, but 4 and 3, respectively, in the second. In convection, the quenching of turbulent pumping follows the same power law as turbulent diffusion, while for the coefficient describing the $\bf Ω\times \bf J$ effect nearly the same quenching exponent is obtained as for $α$. For forced turbulence, turbulent diffusion proportional to the second derivative along the mean magnetic field is quenched much less, especially for larger values of the magnetic Reynolds number. However, we find that in corresponding axisymmetric mean-field dynamos with dominant toroidal field the quenched diffusion coefficients are the same for the poloidal and toroidal field constituents. △ Less

Submitted 20 August, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Comments: Accepted for publication in ApJ

Report number: NORDITA-2014-72

Journal ref: ApJ, 795, 16 (2014)

Confirmation of bistable stellar differential rotation profiles

Authors: P. J. Käpylä, M. J. Käpylä, A. Brandenburg

Abstract: (abridged) Context: Solar-like differential rotation is characterized by a rapidly rotating equator and slower poles. However, theoretical models and numerical simulations can result in a slower equator and faster poles when the rotation is slow. Aims: We study the critical rotational influence under which differential rotation flips from solar-like (fast equator, slow poles) to an anti-solar one… ▽ More (abridged) Context: Solar-like differential rotation is characterized by a rapidly rotating equator and slower poles. However, theoretical models and numerical simulations can result in a slower equator and faster poles when the rotation is slow. Aims: We study the critical rotational influence under which differential rotation flips from solar-like (fast equator, slow poles) to an anti-solar one (slow equator, fast poles). We estimate the non-diffusive ($Λ$ effect) and diffusive (turbulent viscosity) contributions to the Reynolds stress. Methods: We present the results of three-dimensional numerical simulations of mildly turbulent convection in spherical wedge geometry. Here we apply a fully compressible setup which would suffer from a prohibitive time step constraint if the real solar luminosity was used. We regulate the convective velocities by varying the amount of heat transported by thermal conduction, turbulent diffusion, and resolved convection. Results: Increasing the efficiency of resolved convection leads to a reduction of the rotational influence on the flow and a sharp transition from solar-like to anti-solar differential rotation for Coriolis numbers around 1.3. We confirm the recent finding of a large-scale flow bistability: contrasted with running the models from an initial condition with unprescribed differential rotation, the initialization of the model with certain kind of rotation profile sustains the solution over a wider parameter range. Conclusions: Our results may have implications for real stars that start their lives as rapid rotators implying solar-like rotation in the early main-sequence evolution. As they slow down, they might be able to retain solar-like rotation for lower Coriolis numbers before switching to anti-solar rotation. This could partially explain the puzzling findings of anti-solar rotation profiles for models in the solar parameter regime. △ Less

Submitted 16 October, 2014; v1 submitted 13 January, 2014; originally announced January 2014.

Comments: 11 pages, 10 figures, published in Astronomy & Astrophysics

Report number: NORDITA-2014-2

Journal ref: Astron. Astrophys. 570, A43 (2014)