I review our current knowledge of magnetic fields in stars more massive than around 1.5 M⊙, in particular their nature and origin. This includes the strong magnetic fields found in a subset of the population and the fossil field theory invoked to explain them; the subgauss fields detected in Vega and Sirius and their possible origin; and what we can infer about magnetic activity in massive stars and how it might be linked to subsurface convection.

We test the stability of a magnetic equilibrium configuration using numerical simulations and semi-analytical tools. The tested configuration is, as described by Duez & Mathis (2010), the lowest energy state for a given helicity in a stellar radiation zone. We show using 3D magneto-hydrodynamic (MHD) simulations that the present configuration is stable with respect to all submitted perturbations, that would lead to the development of kink-type instabilities in the case of purely poloidal or toroidal fields, both well known to be unstable. We also discuss, using semi-analytic work, the stabilizing influence of one component on the other and show that the found configuration actually lies in the stability domain predicted by a linear analysis of resonant modes.

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred to as the “iron convection zone” (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

We conduct 3D magneto-hydrodynamic (MHD) simulations in order to test the stability of the magnetic equilibrium configuration described by Duez & Mathis (2010). This analytically-derived configuration describes the lowest energy state for a given helicity in a stellar radiation zone. The necessity of taking into account the non force-free property of the large-scale, global field is here emphasized. We then show that this configuration is stable. It therefore provides a useful model to initialize the magnetic topology in upcoming MHD simulations and stellar evolution codes taking into account magneto-rotational transport processes.

The understanding of fossil fields origin, topology, and stability is one of the corner stones of the stellar magnetism theory. On one hand, since they survive on secular time scales, they may modify the structure and the evolution of their host stars. On the other hand, they must have a complex stable structure since it has been demonstrated that the simplest purely poloidal or toroidal fields are unstable on dynamical time scales. In this context, the only stable stellar configurations found today are those resulting from numerical simulations by Braithwaite and collaborators who studied the evolution of an initial stochastic magnetic field, which relaxes with a selective decay of magnetic helicity and energy, on mixed stable configurations (poloidal and toroidal) that seem to be in equilibrium and then diffuse. In this talk, we report the semi-analytical investigation of such an equilibrium field in the axisymmetric case. We use variational methods, which describe selective decay of magnetic helicity and energy during MHD relaxation, and we identify a supplementary invariant due to the stable stratification of stellar radiation zones. This leads to states that generalize force-free Taylor's relaxation states studied in plasma laboratory experiments that become non force-free in the stellar case. Moreover, astrophysical applications are presented and the stability of obtained configurations is studied.