-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Differences in the solar cycle variability of simple and complex active regions during 1996-2018
Authors:
Shabnam Nikbakhsh,
Eija Tanskanen,
Maarit Käpylä,
Thomas Hackman
Abstract:
Aims. Our aim is to examine the solar cycle variability of magnetically simple and complex active region. Methods. We studied simple ($α$ and $β$) and complex ($βγ$ and $βγδ$) active regions based on the Mount Wilson magnetic classification by applying our newly developed daily approach. We analyzed the daily number of the simple active regions (SARs) and compared that to the abundance of the comp…
▽ More
Aims. Our aim is to examine the solar cycle variability of magnetically simple and complex active region. Methods. We studied simple ($α$ and $β$) and complex ($βγ$ and $βγδ$) active regions based on the Mount Wilson magnetic classification by applying our newly developed daily approach. We analyzed the daily number of the simple active regions (SARs) and compared that to the abundance of the complex active regions (CARs) over the entire solar cycle 23 and cycle 24 until December 2018. Results. We show that CARs evolve differently over the solar cycle from SARs. The time evolution of SARs and CARs on different hemispheres also shows differences, even though on average their latitudinal distributions are shown to be similar. The time evolution of SARs closely follows that of the sunspot number, and their maximum abundance was observed to occur during the early maximum phase, while that of the CARs was seen roughly two years later. We furthermore found that the peak of CARs was reached before the latitudinal width of the activity band starts to decease. Conclusions. Our results suggest that the active region formation process is a competition between the large-scale dynamo (LSD) and the small-scale dynamo (SSD) near the surface, the former varying cyclically and the latter being independent of the solar cycle. During solar maximum, LSD is dominant, giving a preference to SARs, while during the declining phase the relative role of SSD increases. Therefore, a preference for CARs is seen due to the influence of the SSD on the emerging flux.
△ Less
Submitted 6 August, 2019;
originally announced August 2019.
-
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.
-
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.
-
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.
-
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.
-
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.
-
$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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Method of frequency dependent correlations: investigating the variability of total solar irradiance
Authors:
Jaan Pelt,
Maarit Käpylä,
Nigul Olspert
Abstract:
This paper contributes to the field of modeling and hindcasting of the total solar irradiance (TSI) based on different proxy data that extend further back in time than the TSI that is measured from satellites.
We introduce a simple method to analyze persistent frequency-dependent correlations (FDCs) between the time series and use these correlations to hindcast missing historical TSI values. We…
▽ More
This paper contributes to the field of modeling and hindcasting of the total solar irradiance (TSI) based on different proxy data that extend further back in time than the TSI that is measured from satellites.
We introduce a simple method to analyze persistent frequency-dependent correlations (FDCs) between the time series and use these correlations to hindcast missing historical TSI values. We try to avoid arbitrary choices of the free parameters of the model by computing them using an optimization procedure. The method can be regarded as a general tool for pairs of data sets, where correlating and anticorrelating components can be separated into non-overlapping regions in frequency domain.
Our method is based on low-pass and band-pass filtering with a Gaussian transfer function combined with de-trending and computation of envelope curves.
We find a major controversy between the historical proxies and satellite-measured targets: a large variance is detected between the low-frequency parts of targets, while the low-frequency proxy behavior of different measurement series is consistent with high precision. We also show that even though the rotational signal is not strongly manifested in the targets and proxies, it becomes clearly visible in FDC spectrum.
The application of the new method to solar data allows us to obtain important insights into the different TSI modeling procedures and their capabilities for hindcasting based on the directly observed time intervals.
△ Less
Submitted 22 December, 2016;
originally announced December 2016.
-
Method for estimating cycle lengths from multidimensional time series: Test cases and application to a massive "in silico" dataset
Authors:
Nigul Olspert,
Maarit Käpylä,
Jaan Pelt
Abstract:
Many real world systems exhibit cyclic behavior that is, for example, due to the nearly harmonic oscillations being perturbed by the strong fluctuations present in the regime of significant non-linearities. For the investigation of such sys- tems special techniques relaxing the assumption to periodicity are required. In this paper, we present the generalization of one of such techniques, namely th…
▽ More
Many real world systems exhibit cyclic behavior that is, for example, due to the nearly harmonic oscillations being perturbed by the strong fluctuations present in the regime of significant non-linearities. For the investigation of such sys- tems special techniques relaxing the assumption to periodicity are required. In this paper, we present the generalization of one of such techniques, namely the D2 phase dispersion statistic, to multidimensional datasets, especially suited for the analysis of the outputs from three-dimensional numerical simulations of the full magnetohydrodynamic equations. We present the motivation and need for the usage of such a method with simple test cases, and present an application to a solar-like semi-global numerical dynamo simulation covering nearly 150 magnetic cycles.
△ Less
Submitted 6 December, 2016;
originally announced December 2016.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
The Maunder minimum (1645--1715) was indeed a Grand minimum: A reassessment of multiple datasets
Authors:
Ilya G. Usoskin,
Rainer Arlt,
Eleanna Asvestari,
Ed Hawkins,
Maarit Käpylä,
Gennady A. Kovaltsov,
Natalie Krivova,
Michael Lockwood,
Kalevi Mursula,
Jezebel O'Reilly,
Matthew Owens,
Chris J. Scott,
Dmitry D. Sokoloff,
Sami K. Solanki,
Willie Soon,
José M. Vaquero
Abstract:
Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity durin…
▽ More
Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum.
Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum.
Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets.
Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level.
△ Less
Submitted 18 July, 2015;
originally announced July 2015.
-
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.
-
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.
-
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.