Weak magnetic fields have recently been detected in a number of A-type stars, including Vega and Sirius. At the same time, space photometry observations of A- and late B-type stars from Kepler and TESS have highlighted the existence of rotational modulation of surface features akin to stellar spots. Here we explore the possibility that surface magnetic spots might be caused by the presence of small envelope convective layers at or just below the stellar surface, caused by recombination of H and He. Using 1D stellar evolution calculations and assuming an equipartition dynamo, we make simple estimates of field strength at the photosphere. For most models the largest effects are caused by a convective layer driven by second helium ionization. While it is difficult to predict the geometry of the magnetic field, we conclude that the majority of intermediate-mass stars should have dynamo-generated magnetic fields of order a few gauss at the surface. These magnetic fields can appear at the surface as bright spots, and cause photometric variability via rotational modulation, which could also be wide-spread in A-stars. The amplitude of surface magnetic fields and their associated photometric variability is expected to decrease with increasing stellar mass and surface temperature, so that magnetic spots and their observational effects should be much harder to detect in late B-type stars.

We review the current state of knowledge of magnetic fields inside stars, concentrating on recent developments concerning magnetic fields in stably stratified (zones of) stars, leaving out convective dynamo theories and observations of convective envelopes. We include the observational properties of A, B and O-type main-sequence stars, which have radiative envelopes, and the fossil field model which is normally invoked to explain the strong fields sometimes seen in these stars. Observations seem to show that Ap-type stable fields are excluded in stars with convective envelopes. Most stars contain both radiative and convective zones, and there are potentially important effects arising from the interaction of magnetic fields at the boundaries between them. Related to this, we discuss whether the Sun could harbour a magnetic field in its core. Recent developments regarding the various convective and radiative layers near the surfaces of early-type stars and their observational effects are examined. We look at possible dynamo mechanisms that run on differential rotation rather than convection. Finally we turn to neutron stars with a discussion of the possible origins for their magnetic fields.

Type I X-ray bursts are produced by thermonuclear runaways that develop on accreting neutron stars. Once one location ignites, the flame propagates across the surface of the star. Flame propagation is fundamental in order to understand burst properties like rise time and burst oscillations. Previous work quantified the effects of rotation on the front, showing that the flame propagates as a deflagration and that the front strongly resembles a hurricane. However the effect of magnetic fields was not investigated, despite the fact that magnetic fields strong enough to have an effect on the propagating flame are expected to be present on many bursters. In this paper we show how the coupling between fluid layers introduced by an initially vertical magnetic field plays a decisive role in determining the character of the fronts that are responsible for the Type I bursts. In particular, on a star spinning at 450 Hz (typical among the bursters) we test seed magnetic fields of $10^{7} - 10^{10}$ G and find that for the medium fields the magnetic stresses that develop during the burst can speed up the velocity of the burning front, bringing the simulated burst rise time close to the observed values. By contrast, in a magnetic slow rotator like IGR J17480--2446, spinning at 11 Hz, a seed field $\gtrsim 10^9$ G is required to allow localized ignition and the magnetic field plays an integral role in generating the burst oscillations observed during the bursts.

The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B stars were collected with the instruments ESPaDOnS at the Canada-France-Hawaii Telescope, Narval at the T\'elescope Bernard Lyot, and HARPSpol at the European Southern Observatory La Silla 3.6m telescope, making MiMeS by far the largest systematic investigation of massive star magnetism ever undertaken. In this paper, the first in a series reporting the general results of the survey, we introduce the scientific motivation and goals, describe the sample of targets, review the instrumentation and observational techniques used, explain the exposure time calculation designed to provide sensitivity to surface dipole fields above approximately 100 G, discuss the polarimetric performance, stability and uncertainty of the instrumentation, and summarize the previous and forthcoming publications.

From pulsar scintillations we infer the presence of sheet-like structures in the ISM; it has been suggested that these are current sheets. Current sheets probably play an important role in heating the solar corona, and there is evidence for their presence in the solar wind. Such magnetic discontinuities have been found in numerical simulations with particular boundary conditions, as well as in simulations using an incompressible equation of state. Here, I investigate their formation under more general circumstances by means of topological considerations as well as numerical simulations of the relaxation of an arbitrary smoothly-varying magnetic field. The simulations are performed with a variety of parameters and boundary conditions: in low, high and of-order-unity plasma-$\beta$ regimes, with periodic and fixed boundaries, with and without a friction force, at various resolutions and with various diffusivities. Current sheets form, over a dynamical timescale, under {\it all} conditions explored. At higher resolution they are thinner, and there is a greater number of weaker current sheets. The magnetic field eventually relaxes into a smooth minimum energy state, the energy of which depends on the magnetic helicity, as well as on the nature of the boundaries.

We present 3D MHD simulations of purely toroidal and mixed poloidal-toroidal magnetic field configurations to study the behavior of the Tayler instability. For the first time the simultaneous action of rotation and magnetic diffusion are taken into account and the effects of a poloidal field on the dynamic evolution of unstable toroidal magnetic fields is included. In the absence of diffusion, fast rotation (rotation rate compared to Alfv\'en frequency) is able to suppress the instability when the rotation and magnetic axes are aligned and when the radial field strength gradient p < 1.5. When diffusion is included, this system turns unstable for diffusion dominated and marginally diffusive dominated regions. If the magnetic and rotation axes are perpendicular to each other the stabilizing effect induced by the Coriolis force is scale dependent and decreases with increasing wavenumber. In toroidal fields with radial field gradients bigger than p > 1.5, rapid rotation does not suppress the instability but instead introduces a damping factor to the growth rate in agreement with the analytic predictions. For the mixed poloidal-toroidal fields we find an unstable axisymmetric mode, not predicted analytically, right at the stability threshold for the non-axisymmetric modes; it has been argued that an axisymmetric mode is necessary for the closure of the Tayler-Spruit dynamo loop.

In a previous study on thermonuclear (type I) nursts on accreting neutron stars we addressed and demonstrated the importance of the effects of rotation, through the Coriolis force, on the propagation of the burning flame. However, that study only analysed cases of longitudinal propagation, where the Coriolis force coefficient $2\Omega\cos\theta$ was constant. In this paper, we study the effects of rotation on propagation in the meridional (latitudinal) direction, where the Coriolis force changes from its maximum at the poles to zero at the equator. We find that the zero Coriolis force at the equator, while affecting the structure of the flame, does not prevent its propagation from one hemisphere to another. We also observe structural differences between the flame propagating towards the equator and that propagating towards the pole, the second being faster. In the light of the recent discovery of the low spin frequency of burster IGR~J17480-2446 rotating at 11 Hz (for which Coriolis effects should be negligible) we also extend our simulations to slow rotation.

In stably stratified stars, numerical magneto-hydrodynamics simulations have shown that arbitrary initial magnetic fields evolve into stable equilibrium configurations, usually containing nearly axisymmetric, linked poloidal and toroidal fields that stabilize each other. In this work, we test the hypothesis that stable stratification is a requirement for the existence of such stable equilibria. For this purpose, we follow numerically the evolution of magnetic fields in barotropic (and thus neutrally stable) stars, starting from two different types of initial conditions, namely random disordered magnetic fields, as well as linked poloidal-toroidal configurations resembling the previously found equilibria. With many trials, we always find a decay of the magnetic field over a few Alfv\'en times, never a stable equilibrium. This strongly suggests that there are no stable equilibria in barotropic stars, thus clearly invalidating the assumption of barotropic equations of state often imposed on the search of magnetic equilibria. It also supports the hypothesis that, as dissipative processes erode the stable stratification, they might destabilize previously stable magnetic field configurations, leading to their decay.

The nature of the mechanisms apparently driving X-rays from intermediate mass stars lacking strong convection zones or massive winds remains poorly understood, and the possible role of hidden, lower mass close companions is still unclear. A 20ks Chandra HRC-I observation of HR 4796A, an 8 Myr old main sequence A0 star devoid of close stellar companions, has been used to search for a signature or remnant of magnetic activity from the Herbig Ae phase. X-rays were not detected and the X-ray luminosity upper limit was L_X =< 1.3x10^27 erg/s. The result is discussed in the context of various scenarios for generating magnetic activity, including rotational shear and subsurface convection. A dynamo driven by natal differential rotation is unlikely to produce observable X rays, chiefly because of the difficulty in getting the dissipated energy up to the surface of the star. A subsurface convection layer produced by the ionisation of helium could host a dynamo that should be effective throughout the main-sequence but can only produce X-ray luminosities of order 10^25 erg/s. This luminosity lies only moderately below the current detection limit for Vega. Our study supports the idea that X-ray production in Herbig Ae/Be stars is linked largely to the accretion process rather than the properties of the underlying star, and that early A stars generally decline in X-ray luminosity at least 100,000 fold in only a few million years.

I review our current knowledge of magnetic fields in stars more massive than around 1.5 Msun, 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.

It is now believed that magnetohydrodynamic equilibria can exist in stably stratified stars due to the seminal works of Braithwaite & Spruit (2004) and Braithwaite & Nordlund (2006). What is still not known is whether magnetohydrodynamic equilibria can exist in a barotropic star, in which stable stratification is not present. It has been conjectured by Reisenegger (2009) that there will likely not exist any magnetohydrodynamical equilibria in barotropic stars. We aim to test this claim by presenting preliminary MHD simulations of barotropic stars using the three dimensional stagger code of Nordlund & Galsgaard (1995).

In gravitationally stratified fluids, length scales are normally much greater in the horizontal direction than in the vertical one. When modelling these fluids it can be advantageous to use the hydrostatic approximation, which filters out vertically propagating sound waves and thus allows a greater timestep. We briefly review this approximation, which is commonplace in atmospheric physics, and compare it to other approximations used in astrophysics such as Boussinesq and anelastic, finding that it should be the best approximation to use in context such as radiative stellar zones, compact objects, stellar or planetary atmospheres and other contexts. We describe a finite-difference numerical scheme which uses this approximation, which includes magnetic fields.

We present the first vertically resolved hydrodynamic simulations of a laterally propagating, deflagrating flame in the thin helium ocean of a rotating accreting neutron star. We use a new hydrodynamics solver tailored to deal with the large discrepancy in horizontal and vertical length scales typical of neutron star oceans, and which filters out sound waves that would otherwise limit our timesteps. We find that the flame moves horizontally with velocities of order $10^5$ cm s$^{-1}$, crossing the ocean in few seconds, broadly consistent with the rise times of Type I X-ray bursts. We address the open question of what drives flame propagation, and find that heat is transported from burning to unburnt fuel by a combination of top-to-bottom conduction and mixing driven by a baroclinic instability. The speed of the flame propagation is therefore a sensitive function of the ocean conductivity and spin: we explore this dependence for an astrophysically relevant range of parameters and find that in general flame propagation is faster for slower rotation and higher conductivity.

Strong magnetic fields play a crucial role in the removal of angular momentum from collapsing clouds and protostellar discs and are necessary for the formation of disc winds as well as jets from the inner disc and indeed, strong large-scale poloidal magnetic fields are observed in protostellar discs at all radii down to $\sim 10 R_\odot$. Nevertheless, by the time the star is visible virtually all of the original magnetic flux has vanished. I explore mechanisms for removing this flux during the formation of the protostar once it is magnetically disconnected from the parent cloud, looking at both radiative and convective protostars. This includes a numerical investigation of buoyant magnetic field removal from convective stars. It is found that if the star goes through a fully convective phase all remaining flux can easily be removed from the protostar, essentially on an Alfv\'en timescale. If on the other hand the protostar has no fully convective phase then some flux can be retained, the quantity depending on the net magnetic helicity, which is probably quite small. Only some fraction of this flux is visible at the stellar surface. I also look at how the same mechanisms could prevent flux from accreting onto the star at all, meaning that mass would only accrete as fast as it is able to slip past the flux.

Weak magnetic fields have recently been detected in Vega and Sirius. Here, we explore the possibility that these fields are the remnants of some field inherited or created during or shortly after star formation and, unlike true fossil fields, are still evolving as we observe them. The timescale of this evolution is given in terms of the Alfven timescale and the rotation frequency by tau_evol ~ tau_A^2 Omega, which would be comparable to the age of the star. It is shown that it is likely that all intermediate- and high-mass stars contain fields of at least the order of the strength found so far in Vega and Sirius. Faster rotators are expected to have stronger magnetic fields. Stars may experience an increase in field strength during their early main-sequence, but for most of their lives field strength will decrease slowly. The length scale of the magnetic structure on the surface may be small in very young stars but should quickly increase to at least very approximately a fifth of the stellar radius.

A remarkable fraction of supernovae (SNe) have no obvious host galaxy. Two possible explanations are that (i) the host galaxy is simply not detected within the sensitivity of the available data or that (ii) the progenitor is a hypervelocity star that has escaped its parent galaxy. We use the Type IIb SN 2009Z as a prototype of case (i), an example of how a very faint (here Low Surface Brightness; LSB) galaxy can be discovered via the observation of a seemingly host-less SN. By identifying and studying LSB galaxies that host SNe related to the death of massive stars, we can place constraints on the stellar population and environment of LSB galaxies, which at present are poorly understood. From an HI spectrum, a redshift of z = 0.02513+-0.00001 and an HI mass of (2.96+-0.12) 10^9 M_sun are computed. This redshift is consistent with that obtained from optical emission lines of SN 2009Z. Furthermore, a gas mass fraction of f_g = 0.87+-0.04 is obtained, one of the highest fractions ever measured. The host galaxy shows signs of recently enhanced star formation activity with a far-UV derived extinction-corrected Star Formation Rate (SFR) of 0.44+-0.34 M_sun/yr. Based on the B-band luminosity we estimate an extinction-corrected metallicity following the calibration by Pilyugin (2001) of 12 + log(O/H) = 8.24+-0.70. The presence of a Type IIb SN in an LSB galaxy suggests, contrary to popular belief, that massive stars can be formed in this type of galaxies. Furthermore, our results imply that LSB galaxies undergo phases of small, local burst activity intermittent with longer phases of inactivity.

Hot luminous stars show a variety of phenomena in their photospheres and winds which still lack clear physical explanation. Among these phenomena are photospheric turbulence, line profile variability (LPV), non-thermal emission, non-radial pulsations, discrete absorption components (DACs) and wind clumping. Cantiello et al. (2009) argued that a convection zone close to the stellar surface could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity associated with iron-group elements and is referred to as the "iron convection zone" (FeCZ). Assuming dynamo action producing magnetic fields at equipartition in the FeCZ, we investigate the occurrence of subsurface magnetism in OB stars. Then we study the surface emergence of these magnetic fields and discuss possible observational signatures of magnetic spots. Simple estimates are made using the subsurface properties of massive stars, as calculated in 1D stellar evolution models. We find that magnetic fields of sufficient amplitude to affect the wind could emerge at the surface via magnetic buoyancy. While at this stage it is difficult to predict the geometry of these features, we show that magnetic spots of size comparable to the local pressure scale height can manifest themselves as hot, bright spots. Localized magnetic fields could be widespread in those early type stars that have subsurface convection. This type of surface magnetism could be responsible for photometric variability and play a role in X-ray emission and wind clumping.

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

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 existence of stable magnetic configurations in white dwarfs, neutron stars and various non-convective stellar {regions} is now well recognized. It has recently been shown numerically that various families of equilibria, including axisymmetric mixed poloidal-toroidal configurations, are stable. Here we test the stability of an analytically-derived non force-free magnetic equilibrium, using three-dimensional magnetohydrodynamic simulations: the mixed configuration is compared with the dynamical evolution of its purely poloidal and purely toroidal components, both known to be unstable. The mixed equilibrium shows no sign of instability under white noise perturbations. {This configuration therefore provides a good description of magnetic equilibrium topology inside non-convective stellar objects and will be useful to initialize magneto-rotational transport in stellar evolution codes.

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.

Observations of clusters of galaxies show ubiquitous presence of X-ray cavities, presumably blown by the AGN jets. We consider magnetic field structures of these cavities. Stability requires that they contain both toroidal and poloidal magnetic fields, while realistic configurations should have vanishing magnetic field on the boundary. For axisymmetric configurations embedded in unmagnetized plasma, the continuity of poloidal and toroidal magnetic field components on the surface of the bubble then requires solving the elliptical Grad-Shafranov equation with both Dirichlet and Neumann boundary conditions. This leads to a double eigenvalue problem, relating the pressure gradients and the toroidal magnetic field to the radius of the bubble. We have found fully analytical stable solutions. This result is confirmed by numerical simulation. We present synthetic X-ray images and synchrotron emission profiles and evaluate the rotation measure for radiation traversing the bubble.

X-ray images of galaxy clusters often display underdense bubbles which are apparently inflated by AGN outflow. I consider the evolution of the magnetic field inside such a bubble, using a mixture of analytic and numerical methods. It is found that the field relaxes into an equilibrium filling the entire volume of the bubble. The timescale on which this happens depends critically on the magnetisation and helicity of the outflow as well as on properties of the surrounding ICM. If the outflow is strongly magnetised, the magnetic field undergoes reconnection on a short timescale, magnetic energy being converted into heat whilst the characteristic length scale of the field rises; this process stops when a global equilibrium is reached. The strength of the equilibrium field is determined by the magnetic helicity injected into the bubble by the AGN: if the outflow has a consistent net flux and consequently a large helicity then a global equilibrium will be reached on a short timescale, whereas a low-helicity outflow results in no global equilibrium being reached and at the time of observation reconnection will be ongoing. However, localised flux-tube equilibria will form. If, on the other hand, the outflow is very weakly magnetised, no reconnection occurs and the magnetic field inside the bubble remains small-scale and passive. These results have implications for the internal composition of the bubbles, their interaction with ICM -- in particular to explain how bubbles could move a large distance through the ICM without breaking up -- as well as for the cooling flow problem in general. In addition, reconnection sites in a bubble could be a convenient source of energetic particles, circumventing the problem of synchrotron emitters having a shorter lifetime than the age of the bubble they inhabit.

We review the interaction in intermediate and high mass stars between their evolution and magnetic and chemical properties. We describe the theory of Ap-star `fossil' fields, before touching on the expected secular diffusive processes which give rise to evolution of the field. We then present recent results from a spectropolarimetric survey of Herbig Ae/Be stars, showing that magnetic fields of the kind seen on the main-sequence already exist during the pre-main sequence phase, in agreement with fossil field theory, and that the origin of the slow rotation of Ap/Bp stars also lies early in the pre-main sequence evolution; we also present results confirming a lack of stars with fields below a few hundred gauss. We then seek which macroscopic motions compete with atomic diffusion in determining the surface abundances of AmFm stars. While turbulent transport and mass loss, in competition with atomic diffusion, are both able to explain observed surface abundances, the interior abundance distribution is different enough to potentially lead to a test using asterosismology. Finally we review progress on the turbulence-driving and mixing processes in stellar radiative zones.

A detailed study spanning approximately a year has been conducted on the Type Ib supernova 2007Y. Imaging was obtained from X-ray to radio wavelengths, and a comprehensive set of multi-band (w2m2w1u'g'r'i'UBVYJHKs) light curves and optical spectroscopy is presented. A virtually complete bolometric light curve is derived, from which we infer a (56)Ni-mass of 0.06 M_sun. The early spectrum strongly resembles SN 2005bf and exhibits high-velocity features of CaII and H_alpha; during late epochs the spectrum shows evidence of a ejecta-wind interaction. Nebular emission lines have similar widths and exhibit profiles that indicate a lack of major asymmetry in the ejecta. Late phase spectra are modeled with a non-LTE code, from which we find (56)Ni, O and total-ejecta masses (excluding He) to be 0.06, 0.2 and 0.42 M_sun, respectively, below 4,500 km/s. The (56)Ni mass confirms results obtained from the bolometric light curve. The oxygen abundance suggests the progenitor was most likely a ~3.3 M_sun He core star that evolved from a zero-age-main-sequence mass of 10-13 M_sun. The explosion energy is determined to be ~10^50 erg, and the mass-loss rate of the progenitor is constrained from X-ray and radio observations to be <~10^-6 M_sun/yr. SN 2007Y is among the least energetic normal Type Ib supernovae ever studied.

In this third paper in a series on stable magnetic equilibria in stars, I look at the stability of axisymmetric field configurations and in particular at the relative strengths of the toroidal and poloidal components. Both toroidal and poloidal fields are unstable on their own, and stability is achieved by adding the two together in some ratio. I use Tayler's (1973) stability conditions for toroidal fields and other analytic tools to predict the range of stable ratios and then check these predictions by running numerical simulations. It is found that while the poloidal field can account for no more than approximately 80% of the total energy, it can account for a very small fraction of the energy, i.e. that the toroidal field can be -- and is likely to be -- significantly stronger than the poloidal. Furthermore, the weaker the field, the weaker the poloidal component can be in relation to the toroidal. The implications of this result are discussed in various contexts such as the emission of gravitational waves by neutron stars, free precession, and a `hidden' energy source for magnetars.

In previous work stable approximately axisymmetric equilibrium configurations for magnetic stars were found by numerical simulation. Here I investigate the conditions under which more complex, non-axisymmetric configurations can form. I present numerical simulations of the formation of stable equilibria from turbulent initial conditions and demonstrate the existence of non-axisymmetric equilibria consisting of twisted flux tubes lying horizontally below the surface of the star, meandering around the star in random patterns. Whether such a non-axisymmetric equilibrium or a simple axisymmetric equilibrium forms depends on the radial profile of the strength of the initial magnetic field. The results could explain observations of non-dipolar fields on stars such as the B0.2 main-sequence star tau-Sco or the pulsar 1E 1207.4-5209. The secular evolution of these equilibria due to Ohmic and buoyancy processes is also examined.

An arbitrary initial magnetic field in an A star evolves into a stable equilibrium. Simulations are presented of the formation of non-axisymmetric equilibria consisting of twisted flux tubes meandering under the surface of the star, and analytic arguments are given relating the stability and form of these equilibria. These results may help to explain observations of Ap stars with very non-dipolar fields. This work is also applicable to other essentially non-convective stars such as white dwarfs and neutron stars.

The stability of large-scale magnetic fields in rotating stars is explored, using 3D numerical hydrodynamics to follow the evolution of an initial poloidal field. It is found that the field is subject to an instability, located initially on the magnetic equator, whereby the gas is displaced in a direction parallel to the magnetic axis. If the magnetic axis is parallel to the rotation axis, the rotation does not affect the initial linear growth of the instability, but does restrict the growth of the instability outside of the equatorial zone. The magnetic energy decays on a timescale which is a function of the Alfv\'en crossing time and the rotation speed, but short compared to any evolutionary timescale. No evidence is found for a possible stable end state to evolve from an initial axisymmetric poloidal field. The field of an oblique rotator is similarly unstable, in both cases regardless of the rotation speed.

We present numerical models of hydromagnetic instabilities under the conditions prevailing in a stably stratified, non-convective stellar interior, and compare them with previous results of analytic work on instabilities in purely toroidal fields. We confirm that an $m=1$ mode (`kink') is the dominant instability in a toroidal field in which the field strength is proportional to distance from the axis, such as the field formed by the winding up of a weak field by differential rotation. We measure the growth rate of the instability as a function of field strength and rotation rate $\Omega$, and investigate the effects of a stabilising thermal stratification as well as magnetic and thermal diffusion on the stability. Where comparison is computationally feasible, the results agree with analytic predictions.

We investigate the 50-year old hypothesis that the magnetic fields of the Ap stars are stable equilibria that have survived in these stars since their formation. With numerical simulations we find that stable magnetic field configurations indeed appear to exist under the conditions in the radiative interior of a star. Confirming a hypothesis by Prendergast (1956), the configurations have roughly equal poloidal and toroidal field strengths. We find that tori of such twisted fields can form as remnants of the decay of an unstable random initial field. In agreement with observations, the appearance at the surface is an approximate dipole with smaller contributions from higher multipoles, and the surface field strength can increase with the age of the star. The results of this paper were summarised by Braithwaite & Spruit (2004).

We use numerical MHD to look at the stability of a possible poloidal field in neutron stars (Flowers & Ruderman 1977), and follow its unstable evolution, which leads to the complete decay of the field. We then model a neutron star after the formation of a solid crust of high conductivity. As the initial magnetic field we use the stable `twisted torus' field which was the result of our earlier work (Braithwaite & Nordlund 2005), since this field is likely to exist in the interior of the star at the time of crust formation. We follow the evolution of the field under the influence of diffusion, and find that large stresses build up in the crust, which will presumably lead to cracking. We put this forward as a model for outbursts in soft gamma repeaters.

We present numerical simulations of a self-sustaining magnetic field in a differentially rotating non-convective stellar interior. A weak initial field is wound up by the differential rotation; the resulting azimuthal field becomes unstable and produces a new meridional field component, which is then wound up anew, thus completing the `dynamo loop'. This effect is observed both with and without a stable stratification. A self-sustained field is actually obtained more easily in the presence of a stable stratification. The results confirm the analytical expectations of the role of Tayler instability.

Some main-sequence stars of spectral type A are observed to have a strong (300-30000 G), static, large-scale magnetic field, of a chiefly dipolar shape -- the `Ap stars' (for example Alioth, the fifth star in the Big Dipper). Following their discovery fifty years ago, it was speculated that these fields a remnants of the star's formation, a `fossil' field. Alternatively, they could be generated by a dynamo process in the star's convective core. The dynamo hypothesis has difficulty explaining high field strengths and the lack of a correlation with rotation. The weakness of the fossil-field theory has been the absence of known field configurations stable enough to survive in a star over its lifetime. We demonstrate here the formation of stable magnetic field configurations, with properties agreeing with those observed, by evolution from arbitrary, unstable initial fields. The results are applicable equally to Ap stars, magnetic white dwarfs and some highly magnetised neutron stars known as magnetars. This establishes fossil fields as the natural, unifying explanation for the magnetism of all these stars.