N. Olspert

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

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

Multiperiodicity, modulations and flip-flops in variable star light curves - II. Analysis of II Pegasus photometry during 1979–2010

Aims. According to previously published Doppler images of the magnetically active primary giant component of the RS CVn binary II Peg, the surface of the star was dominated by one single active longitude that was clearly drifting in the rotational frame of the binary system during 1994-2002; later imaging for 2004–2010, however, showed decreased and chaotic spot activity, with no signs of the drift pattern. Here we set out to investigate from a more extensive photometric dataset whether this drift is a persistent phenomenon, in which case it could be caused either by an azimuthal dynamo wave or be an indication that the binary system’s orbital synchronization is still incomplete. On a differentially rotating stellar surface, spot structures preferentially on a certain latitude band could also cause such a drift, the disruption of which could arise from the change of the preferred spot latitude. Methods. We analyzed the datasets using the carrier fit (CF) method, which is especially suitable for analyzing time series in which a fast clocking frequency (such as the rotation of the star) is modulated with a slower process (such as the stellar activity cycle). Results. We combined all collected photometric data into one single data set and analyzed it with the CF method. We confirm the previously published results that the spot activity has been dominated by one primary spotted region almost through the entire data set and also confirm a persistent, nearly linear drift. Disruptions of the linear trend and complicated phase behavior are also seen, but the period analysis reveals a rather stable periodicity with Pspot = 6. d 71054 ± 0. 00005. After removing the linear trend from the data, we identified several abrupt phase jumps, three of which are analyzed in more detail with the CF method. These phase jumps closely resemble what is called a flip-flop event, but the new spot configurations do not persist for longer than a few months in most cases. Conclusions. There is some evidence that the regular drift without phase jumps is related to the high state, while the complex phase behavior and disrupted drift pattern are related to the low state of magnetic activity. The most natural explanation of the drift is weak anti-solar (pole rotating faster than the equator) differential rotation with a coefficient k ≈ 0.002 combined with the preferred latitude of the spot structure.

Convection-driven spherical shell dynamos at varying Prandtl numbers

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

Multiperiodicity, modulations and flip-flops in variable star light curves - I. Carrier fit method

5-151

Extended Subadiabatic Layer in Simulations of Overshooting Convection

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

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

Context. The light curves of variable stars are commonly described using simple trigonometric models, that make use of the assumption that the model parameters are constant in time. This assumption, however, is often violated, and consequently, time series models with components that vary slowly in time are of great interest. Aims. In this paper we introduce a class of data analysis and visualization methods which can be applied in many different contexts of variable star research, for example spotted stars, variables showing the Blazhko effect, and the spin-down of rapid rotators. The methods proposed are of explorative type, and can be of significant aid when performing a more thorough data analysis and interpretation with a more conventional method. Methods. Our methods are based on a straightforward decomposition of the input time series into a fast “clocking” periodicity and smooth modulating curves. The fast frequency, referred to as the carrier frequency, can be obtained from earlier observations (for instance in the case of photometric data the period can be obtained from independently measured radial velocities), postulated using some simple physical principles (Keplerian rotation laws in accretion disks), or estimated from the data as a certain mean frequency. The smooth modulating curves are described by trigonometric polynomials or splines. The data approximation procedures are based on standard computational packages implementing simple or constrained least-squares fit -type algorithms. Results. Using both artificially generated data sets and observed data, we demonstrate the utility of the proposed methods. Our interest is mainly focused on cases where multiperiodicity, trends or abrupt changes take place in the variable star light curves. Conclusions. The presented examples show that the proposed methods significantly enrich the traditional toolbox for variable star researchers. Applications of the methods to solve problems of astrophysical interest will be presented in the next papers of the series.

Estimating activity cycles with probabilistic methods: I. Bayesian generalised Lomb-Scargle periodogram with trend

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

Method of frequency dependent correlations: investigating the variability of total solar irradiance

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.

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

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 noise model 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 different from constant noise model depends on the density of the data sampling as well as on the true noise model of the process.

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

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.