Only eleven O-type stars have been confirmed to possess large-scale organized magnetic fields. The presence of a -600G longitudinal magnetic field in the O9.7V star HD54879 with a lower limit of the dipole strength of ~2kG, was discovered a few years ago in the framework of the ESO large program "B-fields in OB stars". Our FORS2 spectropolarimetric observations from 2017 October 4 to 2018 February 21 reveal the presence of short and long-term spectral variability and a gradual magnetic field decrease from about -300G down to about -90G. Surprisingly, we discover on the night of 2018 February 17 a sudden, short-term increase of the magnetic field strength and measure a longitudinal magnetic field of -833G. The inspection of the FORS2 spectrum acquired during the observed magnetic field increase indicates a very strong change in spectral appearance with a significantly lower photospheric temperature and a decrease of the radial velocity by several 10km/s. Different scenarios are discussed in an attempt to interpret our observations. The FORS2 radial velocity measurements indicate that HD54879 is a member of a long-period binary.

Although the occurrence of solar irradiance variations induced by magnetic surface features (e.g., sunspots, faculae, magnetic network) is generally accepted, the existence of intrinsic luminosity changes due to the internal magnetic fields is still controversial. This additional contribution is expected to be accompanied by radius variations, and to be potentially significant for the climate of the Earth. We aim to constrain theoretically the radius and luminosity variations of the Sun that are due to the effect of the variable magnetic fields in its interior associated with the dynamo cycle. We have extended a one-dimensional stellar evolution code to include several effects of the magnetic fields on the interior structure. We investigate different magnetic configurations, based on both observational constraints and on the output of state-of-the-art mean field dynamo models. We explore both step-like and simply periodic time dependences of the magnetic field peak strength. We find that the luminosity and radius variations are in anti-phase and in phase, respectively, with the magnetic field strength. For peak magnetic field strengths of the order of tens of kilogauss, luminosity variations ranging between 10^{-6} and 10^{-3} (in modulus) and radius variations between 10^{-6} and 10^{-5} are obtained. Modest but significant radius variations (up to 10^{-5} in relative terms) are obtained for magnetic fields of realistic strength and geometry, providing a potentially observable signature of the intrinsic variations. Establishing their existence in addition to the accepted surface effects would have very important implications for the understanding of solar-induced long-term trends on climate.

Motivated by the empirical finding that the known hemispheric rules for the current helicity at the solar surface are not strict, the excitation of small-scale current helicity by the influence of a large-scale helical magnetic background fields on nonrotating magnetoconvection is demonstrated. It is shown within a quasilinear analytic theory of driven turbulence and by nonlinear simulations of magnetoconvection that the resulting small-scale current helicity has the same sign as the large-scale current helicity while the ratio of both pseudo-scalars is of the order of the magnetic Reynolds number of the convection. The same models do not provide finite values of the small-scale kinetic helicity. On the other hand, a turbulence-induced electromotive force is produced including the diamagnetic pumping term as well as the eddy diffusivity but no alpha effect. It is thus argued that the relations by Pouquet & Patterson (1978) and Keinigs (1983) for the simultaneous existence of small-scale current helicity and alpha effect do not hold for the considered model of nonrotating magnetoconvection. Calculations for various values of the magnetic Prandtl number demonstrate that for the considered diffusivities the current helicity grows for growing magnetic Reynolds number which is not true for the velocity of the diamagnetic pumping -- in agreement with the results of the quasilinear analytical approximation.

The empirical finding with data from the Solar Dynamics Observatory (SDO) of positive (negative) horizontal Reynolds stress at the northern (southern) hemisphere for solar giant cells (Hathaway et al. 2013) is discussed for its consequences for the theory of the solar/stellar differential rotation. Solving the nonlinear Reynolds equation for the angular velocity under neglect of the meridional circulation we show that the horizontal Reynolds stress of the northern hemisphere is always negative at the surface but it is positive in the bulk of the solar convection zone by the action of the Lambda effect. The Lambda effect, which describes the angular momentum transport of rigidly rotating anisotropic turbulence and which avoids a rigid-body rotation of the convection zones, is in horizontal direction of cubic power in the rotation rate $\Omega$ and it is always equatorwards directed. Theories without Lambda effect which may also provide the observed solar rotation law only by the action of a meridional circulation lead to a horizontal Reynolds stress with the opposite sign as observed.

The quasilinear mean-field theory for driven MHD turbulence leads to the result that the observed cross helicity <u \cdot b> may directly yield the magnetic eddy diffusivity \eta_{T} of the quiet Sun. In order to model the cross helicity at the solar surface, magnetoconvection under the presence of a vertical large-scale magnetic field is simulated with the nonlinear MHD code NIRVANA. The very robust result of the calculations is that < u_z b_z >\simeq 2 < u\cdot b> independent of the applied magnetic field amplitude. The correlation coefficient for the cross helicity is about 10%. Of similar robustness is the finding that the rms value of the magnetic perturbations exceeds the mean-field amplitude (only) by a factor of five. The characteristic helicity speed u_{\eta} as the ratio of the eddy diffusivity and the density scale height for an isothermal sound velocity of 6.6 km/s proves to be 1 km/s for weak fields. This value well coincides with empirical results obtained from the data of the HINODE satellite and the Swedish 1-m Solar Telescope (SST) providing the cross helicity component < u_z b_z >. Both simulations and observations thus lead to a numerical value of \eta_{T} \simeq 10^12 cm^2 /s as characteristic for the surface of the quiet Sun.

We study the variation of surface differential rotation and meridional flow along the lower part of the zero age main sequence (ZAMS). We first compute a sequence of stellar models with masses from 0.3 to 1.5 solar masses. We then construct mean field models of their outer convection zones and compute differential rotation and meridional flows by solving the Reynolds equation with transport coefficients from the second order correlation approximation. For a fixed rotation period of 2.5 d we find a strong dependence of the surface differential rotation on the effective temperature with weak surface shear for M dwarfs and very large values for F stars. The increase with effective temperature is modest below 6000 K but very steep above 6000 K. The meridional flow shows a similar variation with temperature but the increase with temperature is not quite so steep. Both the surface rotation and the meridional circulation are solar-type over the entire temperature range. We also study the dependence of differential rotation and meridional flow on the rotation period for masses. from 0.3 to 1.1 solar masses. The variation of the differential rotation with period is weak except for very rapid rotation. The meridional flow shows a systematic increase of the flow speed with the rotation rate. Numerical experiments in which either the $\Lambda$ effect is dropped in the Reynolds stress or the baroclinic term in the equation of motion is cancelled show that for effective temperatures below 6000 K the Reynolds stress is the dominant driver of differential rotation.

We present models of temperature distribution in the crust of a neutron star in the presence of a strong toroidal component superposed to the poloidal component of the magnetic field. The presence of such a toroidal field hinders heat flow toward the surface in a large part of the crust. As a result, the neutron star surface presents two warm regions surrounded by extended cold regions and has a thermal luminosity much lower than in the case the magnetic field is purely poloidal. We apply these models to calculate the thermal evolution of such neutron stars and show that the lowered photon luminosity naturally extends their life-time as detectable thermal X-ray sources.

The differential rotation of the F5V-IV star Procyon A is computed for a class of models which are consistent with recent astrometric and asteroseismic data. The rotation pattern is determined by solving the Reynolds equation for motion, including the convective energy transport, where the latter is anisotropic owing to the Coriolis force action which produces a horizontal temperature gradient between the poles and the equator. All the models show a decrease of the rotation rate with increasing radius and solar-type isorotation surfaces with the equator rotating faster than the poles, the horizontal rotational shear being much smaller for models with a less extended convective envelope. The meridional flow circulation can be either clockwise or counter-clockwise, and in some cases a double latitudinal cell appears. The rotational splittings are calculated for low degree $p$-modes with $l=1, m=1$ and $l=2, m=1,2$, and it is shown that, for modes with $m=1$, the stronger is the horizontal differential rotation shear the weaker the effect on the average rotational splitting expected, whilst the opposite happens for the mode with $m=2$. On the basis of the present study, a resolution of $10 \rm nHz$ in individual oscillation frequencies seems to be necessary to test the different dynamical behaviour of the proposed models, that appears barely achievable even in the forthcoming space missions. However, the average over several observed splittings could produce the required accuracy.

We continue the study of the effects of a strong magnetic field on the temperature distribution in the crust of a magnetized neutron star (NS) and its impact on the observable surface temperature. Extending the approach initiated in Geppert et al.(2004), we consider more complex and, hence, more realistic, magnetic field structures but still restrict ourselves to axisymmetric configurations. We put special emphasis on the heat blanketing effect of a toroidal field component. We show that asymmetric temperature distributions can occur and a crustal field consisting of dipolar poloidal and toroidal components will cause one polar spot to be larger than the opposing one. These two warm regions can be separated by an extended cold equatorial belt. We present an internal magnetic field structure which can explain both the X-ray and optical spectra of the isolated NS RXJ 1856-3754. We investigate the effects of the resulting surface temperature profiles on the observable lightcurve which an isolated thermally emitting NS would produce for different field geometries. The lightcurves will be both qualitatively (deviations from sinusoidal shape) and quantitatively (larger pulsed fraction for the same observational geometry) different from those of a NS with an isothermal crust. This opens the possibility to determine the interna} magnetic field strengths and structures in NSs by modeling their X-ray lightcurves and spectra. The striking similarities of our model calculations with the observed spectra and pulse profiles of isolated thermally emitting NSs is an indication for the existence of strong magnetic field components maintained by crustal currents.

We develop a model based on 3D mean-field MHD for the generation of large scale magnetic fields in fully convective objects like low-mass stars, brown dwarfs and possibly gaseous planets. The dynamo process is of alpha^2 type and thus differs from the shell-dynamo at work in more massive stars.The alpha^2 dynamo becomes supercritical for Rossby numbers Ro\la 10. It generates a large-scale, non-axisymmetric, steady field that is symmetric with respect to the equatorial plane. Saturation of the alpha^2-generated field at the equipartition yields strengths of several kGauss, in agreement with observations of active M dwarfs, and provides a qualitative explanation for the observed activity saturation in late M stars. For brown dwarfs with a conductive core, as occurs at the center of the most massive and oldest of these objects, we have also studied an alpha^2-Omega dynamo, i.e. the effect of differential rotation. In this case the field is predominantly toroidal, axisymmetric and oscillatory, like the solar field. The topology of the field in the fully convective objects exhibits a high order multipole character that differs from the aligned dipole field generated by the alpha-Omega dynamo. The strong reduction of the dipolar component due to the field non-axisymmetry should considerably reduce the Alfven radius and thus the efficiency of magnetic braking, providing an appealing explanation for the decreasing angular momentum loss rate observed in low-mass stars and brown dwarfs. This may have also implications for cataclysmic variables below the period gap. In spite of this large-scale field, the decreasing conductivity in the dominantly neutral atmosphere of these objects prevents the current generation necessary to support a chromosphere and thus activity. (Abridged)

We compute the differential rotation of main sequence stars of the spectral types F, G, K, and M by solving the equation of motion and the equation of convective heat transport in a mean-field formulation. For each spectral type the rotation rate is varied to study the dependence of the surface shear on this parameter. The resulting rotation patterns are all solar-type. The horizontal shear turns out to depend strongly on the effective temperature and only weakly on the rotation rate. The meridional flow depends more strongly on the rotation rate and has different directions in the cases of very slow and very fast rotation, respectively.

The differential rotation of a 1.2 $M_\odot$ zero age MS star (spectral type F8) is computed and the results are compared with those from a similar model of the Sun. The rotation pattern is determined by solving the Reynolds equation including the convective energy transport. The latter is anisotropic due to the Coriolis force causing a horizontal temperature gradient of ~ 7 K between the poles and the equator. Comparison of the transport mechanisms of angular momentum (the eddy viscosity, the $\Lambda$-effect and the meridional flow) shows that for the F star the $\Lambda$-effect is the most powerful transporter for rotation periods of 7 d or less. In the limit of very fast rotation the $\Lambda$-effect is balanced by the meridional flow alone and the rotation is nearly rigid. The rotation pattern found for the F star is very similar to the solar rotation law, but the horizontal shear is about twice the solar value. As a function of the rotation period, the total equator-pole difference of the angular velocity has a (slight) maximum at a period of 7 d and (slowly) vanishes in both the limiting cases of very fast and very slow rotation. A comparison of the solar models with those for the F-type star shows a much stronger dependence of the differential surface rotation on the stellar luminosity rather than on the rotation rate.

The three components of the heat-flux vector $F =\rho C_p <u^\prime T^\prime> are numerically computed for a stratified rotating turbulent convection using the NIRVANA code in a flat box. The latitudinal component $F_\theta$ proves to be negative (positive) in the northern (southern) hemisphere so that the heat always flows towards the poles. As a surprise, the radial heat-flux $F_r$ peaks at the equator rather than at the poles (Taylor numbers O(10^6)). The same behavior is observed for the radial turbulence intensity $<u_r^{\prime 2}>$ which for \emph{free} turbulence is also believed to peak at the poles (see Eq. (19) below). As we can show, however, the consequences of this unexpected result (also obtained by Kaepylae, Korpi and Tuominen 2004) for the theory of differential rotation are small as mainly the $F_\theta$ is responsible to solve the `Taylor number puzzle'. In all our simulations the azimuthal component $F_\phi$ proves to be negative so that the rotating turbulence produces an westwards directed azimuthal heat-flux which should be observable. Fluctuations with higher temperature are expected to be anticorrelated with their own angular velocity fluctuations. We find this rotation-induced result as understandable as the $F_\phi$ is closely related to the radial \Lambda-effect which is known to be also negative in stratified and rapidly rotating convection zones.

We investigate the influence of different magnetic field configurations on the temperature distribution in neutron star crusts. We consider axisymmetric dipolar fields which are either restricted to the stellar crust, ``crustal fields'', or allowed to penetrate the core, ``core fields''. By integrating the two-dimensional heat transport equation in the crust, taking into account the classical (Larmor) anisotropy of the heat conductivity, we obtain the crustal temperature distribution, assuming an isothermal core. Including quantum magnetic field effects in the envelope as a boundary condition, we deduce the corresponding surface temperature distributions. We find that core fields result in practically isothermal crusts unless the surface field strength is well above $10^{15}$ G while for crustal fields with surface strength above a few times $10^{12}$ G significant deviations from isothermality occur at core temperatures inferior or equal to $10^8$ K. At the stellar surface, the cold equatorial region produced by the quantum suppression of heat transport perpendicular to the field in the envelope, present for both core and crustal fields, is significantly extended by the classical suppression at higher densities in the case of crustal fields. This can result, for crustal fields, in two small warm polar regions which will have observational consequences: the neutron star has a small effective thermally emitting area and the X-ray pulse profiles are expected to have a distinctively different shape compared to the case of a neutron star with a core field. These features, when compared with X-ray data on thermal emission of young cooling neutron stars, will open a way to provide observational evidence in favor, or against, the two radically different configurations of crustal or core magnetic fields.