A Thorough Investigation of Distance and Age of the Pulsar Wind Nebula 3C58

Abstract: A growing number of researchers present evidence that the pulsar wind nebula 3C58 is much older than predicted by its proposed connection to the historical supernova of A.D. 1181. There is also a great diversity of arguments. The strongest of these arguments rely heavily on the assumed distance of 3.2 kpc determined with HI absorption measurements. This publication aims at determining a more accurate distance for 3C58 and re-evaluating the arguments for a larger age. I have re-visited the distance determination of 3C58 based on new HI data from the Canadian Galactic Plane Survey and our recent improvements in the knowledge of the rotation curve of the outer Milky Way Galaxy. I have also used newly determined distances to objects in the neighbourhood, which are based on direct measurements by trigonometric parallax. I have derived a new more reliable distance estimate of 2 kpc for 3C58. This makes the connection between the pulsar wind nebula and the historical event from A.D. 1181 once again much more viable.

\section{Introduction}

Of all the supernova remnants (SNRs) that are linked to historically observed
supernova events the most disputed one is probably the connection between
the pulsar wind nebula (PWN) 3C\,58 and the supernova explosion observed
in A.D.~1181 by Chinese and Japanese astronomers. \citet{step71} and
\citet{step02} claim that there is a
high probability for a connection between the Guest Star from A.D.~1181 and the
supernova explosion that created 3C\,58. The length of the visibility of the
Guest Star and the similarity in the description in independent sources
indicate a supernova explosion. The lack of any other supernova remnant
candidate in the area strongly supports the identification of 3C\,58 with
this historical supernova explosion. Although it is difficult to escape
this argument, strong evidence against such a young age for
3C\,58 has been mounting with a great diversity of support from expansion
studies, theoretical
modeling of the evolutionary path, and comparisons of the pulsar's and the
pulsar wind nebula's characteristics to other SNRs and PWNe of known age
\citep[for a list of the main arguments see e.g.][Table 3]{fese08}.

Many of the strong arguments against a young age for 3C\,58 rely heavily
on the assumed distance of 3.2~kpc \citep{robe93}. This distance was
determined kinematically from H~{\sc{i}} absorption measurements by
comparing the resulting systemic velocity for 3C\,58 with a flat
rotation curve for the Outer Galaxy. For
Perseus arm objects in particular this can lead to a significant overestimate of the
distance. A spiral shock in the Perseus arm is ``pushing''
objects, including H~{\sc{ii}} regions, SNRs, and PWNe towards the Galactic centre,
giving them -- from our point of view -- a
higher negative radial velocity, which makes them appear to
be farther away than they actually are \citep{robe72}. Examples of this
effect on distance estimates of SNRs and PWNe can be found, e.g.
in \citet{koth02,koth03} and \citet{fost04} and most recently in \citet{koth12}. SNRs and PWNe
in particular rely heavily on kinematic distance determinations, since
their distances cannot easily be determined
directly or by relating them to embedded stars as in the case for
H~{\sc{ii}} regions.

In this paper I discuss a new distance estimate for 3C\,58.
Based on this newly determined distance I will re-evaluate
evidence presented in the literature against the historical connection between
3C\,58 and the Guest Star of A.D. 1181. I will show that the new distance
changes some of the PWN's characteristics quite dramatically. This leads
once more to a higher probability for its historical connection. I will
describe the H~{\sc{i}} observations of the Canadian Galactic Plane Survey in
Section 2. H~{\sc{i}} emission around 3C\,58, its newly determined H~{\sc{i}}
absorption characteristics, and the resulting systemic velocity are discussed
in Section 3 and a thorough evaluation of the distance is presented in
Section 4. In section 5 I will re-visit all arguments contra and pro a young
age, develop a model in Section 6 for the evolution of 3C\,58 that
explains all observations, and give a short summary in Section 7.
Preliminary results of this project were published in \citet{koth10}.

\section{Observations}

\begin{figure*}[!ht]
\begin{center}
  \includegraphics[bb = 100 260 475 625,width=0.80\textwidth,clip]{f1.eps}
\end{center}
 \caption{{\bf Top}: H~{\sc{i}} absorption profiles of 3C\,58. The relative absorption
     of the radio continuum signal at 1420~MHz is displayed as a function of radial
     velocity $v_{LSR}$. Left: Absorption profiles of 3 different positions are
     displayed as indicated in the bottom right panel. Right: Profile
     calculated using all pixels on 3C\,58 with $T_b \ge 100$\,K. Each pixel was
     weighted by its intensity.
     {\bf Bottom}: Left: H~{\sc{i}} emission profiles used to determine the
     absorption profiles in the top left. These were averaged
     over the positions marked in the right panel. Right: Continuum image of
     3C\,58 taken from the CGPS. The "on" and "off" positions used to calculate the
     H~{\sc{i}} absorption profiles are indicated by solid and dashed circles,
     respectively.}
   \label{abs}
\end{figure*}

The data for this new H~{\sc{i}} absorption study towards 3C\,58
were obtained with the Synthesis Telescope at the Dominion
Radio Astrophysical Observatory \citep[DRAO,][]{land00} as part of the Canadian
Galactic Plane Survey \citep[CGPS,][]{tayl03}. Single antenna data were
incorporated into the interferometer maps after suitable filtering in the
Fourier domain. This assures the accurate representation of atomic hydrogen
emission from the largest
structures down to the resolution limit of about 1\arcmin. The low spatial
frequency data were drawn from the Low-Resolution DRAO Survey (LRDS) of the CGPS
region, which was observed with the 26-m radio telescope at DRAO
\citep{higg00}. The resolution in the final data sets varies slightly across the
images as about $1\arcmin \times 1\arcmin {\rm cosec(Declination)}$. At the
centre of 3C\,58 the resolution of the CGPS H~{\sc{i}} data is $58\arcsec \times 65\arcsec$
at an angle of $75\degr$ for the major axis (counter-clockwise from the
Galactic longitude axis). The RMS noise is about 3\,K $T_B$ in each velocity
channel of width 0.82446\,km~s$^{-1}$. The frequency band used for the CGPS
H~{\sc{i}} observations is 1\,MHz wide, which gives a velocity coverage of about
210\,km\,\,s$^{-1}$ at a velocity resolution of 1.3\,km~s$^{-1}$. In this area
of the Milky Way Galaxy the CGPS H~{\sc{i}} data are centred at a radial
velocity of $-60$\,km~s$^{-1}$ relative to the local standard of rest.

\section{Results}

\subsection{H~{\sc{i}} Absorption Towards 3C\,58}

The first attempt to derive an H~{\sc{i}} absorption distance to 3C\,58 was carried
out by \citet{will73}, who found absorption at negative velocities down to
$-95$\,km\,\,s$^{-1}$, which
was translated to a distance of 8.2~kpc. However, those observations were obtained
with a single-dish telescope and the ``off position'' to determine the background
H~{\sc{i}} emission for the absorption profile
was taken $9\arcmin$ away from the centre of 3C\,58. \citet{gree82} who
did the first absorption study towards 3C\,58 with an interferometer
argue that most of the absorption components at high negative velocities in the work
of \citet{will73} are artificial due to the highly structured interstellar medium. They
found no absorption beyond $-34$~km~s$^{-1}$ and determined a systemic velocity of
about $-37$~km~s$^{-1}$, which they translated with a Schmidt rotation model
\citep{schm65} for our Galaxy to a distance of 2.6~kpc.
The latest distance
determination for 3C\,58 by H~{\sc{i}} absorption measurements resulted in a
systemic velocity of
$\sim -38$\,km~s$^{-1}$ and a Perseus arm location \citep{robe93}. This was translated
with a flat rotation model for our Galaxy
and the IAU supported values for the Sun's Galacto-centric distance of
$R_\odot = 8.5$~kpc and the Sun's circular motion around the Galactic centre of
$v_\odot = 220$\,km~s$^{-1}$ to a distance of 3.2~kpc.

\begin{figure*}[!ht]
\begin{center}
 \includegraphics[bb = 35 125 555 480,width=0.85\textwidth,clip]{f2.eps}
\end{center}
\caption{H~{\sc{i}} channel maps of 3C\,58 taken at the peaks in the
absorption spectra at +1 (top left), $-$7 (top centre), $-$34 (bottom left), and
$-$41~km\,s$^{-1}$ (bottom centre). The right panels show 10 channels averaged over
the inter-arm velocities between the Local and Perseus arm (top) and 10 channels
averaged over the velocities just outside the last absorption feature.
The continuum emission of 3C\,58 is indicated
by the white contours. The black dashed lines represent the H~{\sc{i}} contours at the
levels indicated by the labels of the colour bars.}
\label{abs2}
\end{figure*}

I used the H~{\sc{i}} data of the CGPS to conduct a new H~{\sc{i}} absorption study
towards 3C\,58 to investigate not only the integrated absorption spectrum, but
also any possible significant changes over the PWN: if absorbing clouds are
small then the absorption by those clouds may be washed out in an integrated
profile. This could cause the loss of important information about the
velocity range in which 3C\,58 is absorbed. In Fig.~\ref{abs} I show three absorption
profiles towards individual positions on 3C\,58. To calculate the displayed
spectra I subtracted the averaged emission profiles as indicated on the map
in the lower right panel of Fig.~\ref{abs}. For the integrated absorption
profile I averaged over all pixels on 3C\,58 that display a surface brightness in excess
of 100\,K~$T_b$ in the 1420\,MHz continuum maps of the CGPS. Each pixel was
weighted by its brightness above the background. For the background
emission profile the average over all locations marked in Fig.~\ref{abs} was
calculated. The resulting profiles are shown in Fig.~\ref{abs}.
H~{\sc{i}} channel maps
taken at the major peaks of the integrated absorption spectrum are displayed in
Fig.~\ref{abs2}. These peaks can be found at $+1$ and $-7$\,km~s$^{-1}$ which
represents local gas, at $-34$\,km~s$^{-1}$ representing the Perseus spiral arm
and for a small area in the right part of 3C\,58 at $-41$\,km~s$^{-1}$.
In addition I added two maps averaged over 10 velocity
channels in the
inter-arm area between the Local and the Perseus arm and 10 velocity channels
just beyond the last visible absorption component at $-41$\,km\,\,s$^{-1}$.

The integrated absorption spectrum is almost identical to the one published
by \citet{robe93}. Every peak and wiggle is perfectly reproduced, which is
not surprising since 3C\,58 is a very bright source. At the local velocities and
in the inter-arm area between the Local arm and Perseus arm the absorption looks quite
constant over the PWN. However, I found that the weak peak, seen as a faint wiggle
at around $-41$~km~s$^{-1}$ is actually a real absorption feature. It is not
seen over the entire PWN, but only on the right hand side at position 3
(see Figs.~\ref{abs} and \ref{abs2}).
This is likely the reason why this feature was not obvious in previously
published integrated absorption spectra. The Perseus arm absorption component
around $-34$\,km\,\,s$^{-1}$ seems to be confined to the left hand side of the
PWN and is absent from the lower right. The reason for this is probably the structured
interstellar medium. In the local medium clouds have a bigger angular size than
in the Perseus arm simply because they are much closer. Hence the absorbing clouds in
the Perseus arm show structure on smaller angular scale.

\begin{figure*}[!ht]
\begin{center}
 \includegraphics[bb = 25 25 535 740,width=0.85\textwidth,clip]{f3new.eps}
\end{center}
\caption{H~{\sc{i}} channel maps of the area around 3C\,58. Each channel is
1.0~km\,s$^{-1}$ wide and the centre velocity is indicated. The location
of 3C\,58 is marked by the white asterisk. The black ellipse indicates the
approximate extend of the bipolar bubble discovered by \citet{wall94}.}
\label{hiemi}
\end{figure*}

\subsection{H~{\sc{i}} Emission Around 3C\,58}

\citet{wall94} discussed the possible association of 3C\,58 with a large
bi-polar interstellar bubble which they detected in low resolution H~{\sc{i}}
data.
An inspection of the CGPS H~{\sc{i}} data sets reveals that this double-shell
structure is actually a relatively smooth object (see Fig.~\ref{hiemi}).
This
bubble is visible from about $-33$ to $\bf -39$\,km\,s$^{-1}$ and its
structure is
essentially unchanging over that entire velocity range. No caps are obvious in
emission, neither at more nor at less negative velocities. The systemic velocity
of this structure is about $-36$\,km\,s$^{-1}$. It seems to appear at
about $-32$ to $-33$\,km\,s$^{-1}$ without any dynamical connection to all the
other filaments that are visible in this velocity range and disappears the same
way at about $-39$ to $-40$\,km\,s$^{-1}$.
The fact that this bubble seems to be stationary through many velocity channels
and does not show any cap moving towards or away from us does not contradict
its interpretation as an expanding bubble. \citet{cazz05} simulated H~{\sc{i}}
observations of expanding stellar wind bubbles and found that this is
not unexpected.
Stationary rings or shells and weak or undetectable caps are the result of velocity
dispersion. A very nice example for an actual observation of such a
structure is the stellar wind bubble around the SNR 3C434.1, where the receding
cap is seen absorbing warm background H~{\sc{i}} emission \citep{fost04}.

3C\,58 can be found projected onto the top left area of this bubble. If 3C\,58
was located inside this structure, the last absorption feature at about
$-41$\,km\,s$^{-1}$ must be produced by the cap of this bubble which is
expanding towards us. After correction for the off-centre position of 3C\,58 this would result
in an expansion velocity of 6 to 7\,km~s$^{-1}$ for the bubble.

\subsection{The Systemic Velocity of 3C\,58}

\citet{robe93} found a systemic velocity of about
$v_{LSR} = -38$\,km\,s$^{-1}$ (LSR = Local Standard of Rest), which they translated
to a distance of 3.2\,kpc, in the Perseus arm.
This systemic velocity was confirmed with the discovery of the large bubble
by \citet{wall94}. The CGPS H~{\sc{i}} absorption profile (Fig.~\ref{abs}) shows an
additional absorption feature at about $-41$\,km\,s$^{-1}$ and a systemic
velocity of about $-36$\,km~s$^{-1}$ for the bubble. There is no absorption
at velocities beyond $-41$\,km~s$^{-1}$.
In Fig.~\ref{abs2} the map averaged
around $-25$\,km~s$^{-1}$, representative of the inter-arm between the Local arm and
the Perseus arm shows significant absorption. However, the map averaged over velocities
beyond $-41$\,km~s$^{-1}$ is free of absorption even though this velocity range
should represent a location inside a spiral arm with denser gas than
the inter-arm around $-25$\,km~s$^{-1}$.
This confirms the Perseus arm location of 3C\,58. 3C\,58 is most likely
located inside the bubble shown in Fig.~\ref{hiemi} and the last absorption feature,
which is only partly visible in the absorption profile of 3C\,58 is a part of that
bubble moving
towards us. The emission profiles show bright emission peaks for all major absorption
components indicating that those are major parts of the ISM. The only exception
is the absorption feature at $-41$\,km~s$^{-1}$, which has no obvious
related emission component. There is some enhanced H~{\sc{i}} emission
around $-41$\,km~s$^{-1}$ in off-positions 2 and 3 (see Fig.~\ref{abs}),
but only the absorption profile in position 3 shows absorption. This
indicates that this absorption feature must be caused by a small
localized cloud, which
is not part of the grand kinematic scheme of the Galaxy.

This would, however, be typical for the broken cap of an expanding bubble,
which are very difficult to
detect due to velocity dispersion \citep{cazz05}. However, such a cap could
reveal
itself in absorption in particular with a bright background source like 3C\,58.
\citet{koth02b} show many examples of expanding H~{\sc{i}} shells around
compact H~{\sc{ii}} regions which are revealed by H~{\sc{i}} absorption of the
H~{\sc{ii}} region's radio continuum emission by the approaching shell.

Based on these results I will adopt a systemic velocity of about
$-36$\,km~s$^{-1}$ for 3C\,58, which is consistent with earlier estimates.
This is independent of whether 3C\,58 is located inside the bubble or
not, because the absorption feature at $-41$\,km~s$^{-1}$ represents the
absorption by a small cloud that would not be taken into account for
the determination of a systemic velocity. And the last absorption
feature, related to a major structure of the ISM, seen in bright
supernova remnants is typically very close to the actual
systemic velocity of these objects \citep[e.g.][]{koth03}.

\section{The Distance to 3C\,58}

The method of \citet{fost06} was used to determine a new distance-velocity relation in
the direction of 3C\,58 (Fig.~\ref{dvplot}). This is a kinematic-based distance
tool, which accounts for not only circular motion of the Galaxy as in a
flat rotation curve model, but also for radial motions produced by spiral shocks.
To determine the distance to 3C\,58, \ion{H}{i} profiles were obtained from a
$10\degr \times 10\degr$ CGPS mosaic. This mosaic was first processed by a new
cloud filtering technique described in \citet{fost10}, which produces a
datacube of the smooth intercloud medium by removing unwanted small clouds that may not follow
large-scale Galactic dynamics.

\begin{figure}[!ht]
\begin{center}
 \includegraphics[bb = 55 45 535 500,width=0.48\textwidth,clip]{f4.eps}
\end{center}
\caption{Plot of the distance-velocity relation in the direction of 3C\,58
determined with the method by \citet{fost06} (solid line). The flat rotation
curve with $R_\odot = 8.5$~kpc and $v_\odot = 220$\,km~s$^{-1}$ used by
\citet{robe93} is indicated
by the dashed line. The dotted lines mark the velocity range observed for
the bubble around 3C\,58.}
\label{dvplot}
\end{figure}

To produce the distance-velocity relation displayed in Fig.~\ref{dvplot}, \ion{H}{i}
profiles along ten separate lines of sight (LOS) around 3C\,58 and a location in the
midplane at the same longitude were fitted independently.
For the systemic velocity of 3C\,58, between $-35$ and $-36$\,km\,s$^{-1}$, the only possible location
is at the velocity discontinuity produced by the Perseus arm spiral shock at 2\,kpc.
The distance to the spiral shock varies only slightly between the different
lines of sight, resulting in a very small uncertainty. This leads to a distance of
$2.0\pm 0.3$~kpc for 3C\,58.

There is a second independent method to determine a distance to 3C\,58:
relate the PWN to a nearby object or nearby objects which themselves have an independent or a more
reliable distance estimate. The W\,3/4/5 H~{\sc{ii}} region
complex and the related SNR HB\,3 are just a few degrees away from 3C\,58. Both, HB\,3
and W\,3/4/5, have very similar systemic velocities \citep[e.g.][]{rout91}, indicating that they are all
located in Perseus spiral arm. The distances to HB\,3, W\,5, and the nearby \ion{H}{ii} region Sh\,2-190
were also determined by
\citet{fost06} with the same method used here for 3C\,58. For all three objects the distance was found to
be $2.3\pm 0.5$~kpc. This is consistent with our current result. The descrepancy in
the results and errors is likely due to the application of the cloud filter in our study, which
significantly reduces the uncertainties caused by small-scale clouds.

For W\,3 the distance was independently determined with trigonometric
parallax to related masers to be $1.95\pm0.04$\,kpc and $2.04\pm0.07$\,kpc
by \citet{xu06} and \citet{hach06}, respectively. This very nicely confirms the 2\,kpc
I estimated for 3C\,58 from the systemic velocity of its host stellar wind bubble.

I believe there is compelling evidence for a location of 3C\,58 in the Perseus arm
spiral shock and from now on I will assume that 3C\,58 is located at a
distance $d = 2\pm 0.3$~kpc.

\subsection{The Progenitor of 3C\,58}

Interstellar material is compressed in the spiral shock, forms molecules, and later
stars. The shock provides the perfect environment for the formation of massive stars that later
explode as a supernova like the progenitor star of 3C\,58. It takes a long time
for those stars -- some $10^8$ years -- to migrate
to a position beyond the Perseus arm spiral shock
\citep[e.g.][for similar calculations]{arvi09,koth12}. At a distance of 2\,kpc the average radius
of the bubble would be about 35\,pc. With an average expansion velocity of
6.5\,km\,s$^{-1}$ this results in a dynamic age of $3.2\times 10^6$\,yr
for a pure stellar wind bubble produced by a single star \citep{weav77}.
Considering the age of this bubble the location of 3C\,58, about $30\arcmin$ away
from the geometric centre, requires a velocity of about 5 to 6\,km\,s$^{-1}$ perpendicular
to the line of sight for the progenitor star, which is certainly not unusual.

From the age of the stellar wind bubble we can determine approximately the type of the progenitor
star assuming that it alone was responsible for the bubble's formation. A lifetime of
$3.2\times 10^6$\,yr would point to a late O-type star with an initial mass between
20 and 30~M$_\odot$.

\section{The Age of 3C\,58}

There is a large number of arguments against and only very few for a connection
between the PWN 3C\,58 and the historical supernova explosion from A.D.~1181. A quite
comprehensive list can be found in \citet[Table~3]{fese08}. The discussion
in the literature of the age of 3C\,58 has been based on an assumed distance of
3.2~kpc. I will discuss
and re-evaluate all of the arguments, pro and contra a young age based on the
newly determined distance of 2~kpc. The results of the discussion below
are summarized in Table~\ref{tab:age}. In the end I will try to create
a compelling picture that explains the results of all observations and
theoretical models.

\subsection{The Guest Star from A.D. 1181}

The strongest arguments for the connection between 3C\,58 and the Guest Star
of A.D.~1181 remain those by \citet{step71} and \citet{step02}.
The length of the visibility of the Guest Star and the similarity of the
description in the four independent Chinese and Japanese sources
indicate that this Guest star was the result of a supernova explosion.
The lack of any
other supernova remnant candidate for the event other than 3C\,58 strongly
supports the identification of 3C\,58 with the supernova
remnant, in particular, since such a recent event should leave an
easily detectable remnant behind. The connection of 3C\,58 with the Guest
Star of A.D. 1181 implies an age of about 830~yr.

\subsubsection{Supernova Peak Brightness}

According to \citet{step02} the supernova had at its peak a visual apparent brightness
of about 0 magnitudes, maybe somewhat brighter. This would translate to an absolute brightness
between -13.4 and -14.3~mag, using the range of possible values for visual extinction
$A_v$ \citep[1.9 to 2.8,]{shib08} and the newly determined distance of 2~kpc. This is a
very low brightness
even if we assume an additional uncertainty of $\pm 1$~mag.

However, other low luminosity
supernova explosions have been observed before.
\citet{spir09} have reviewed the literature about underluminous type II-P supernovae
and reach the conclusion that those are a quite homogenous group of supernovae with
low peak luminosities, red colours, slow moving ejecta, and small ejected $^{56}$Ni mass.
However, the presence of the stellar wind bubble around 3C\,58 indicates that the progenitor
star produced a strong stellar wind implying a supernova of type Ib or Ic, or at least
a type IIL or IIb supernova since the progenitor must have lost a significant amount of material
through that wind, unless the bubble is unrelated to 3C\,58, in which case it could
have been a member of the underluminous type II-P group.

There have been observations of supernovae of type IIL and IIb which showed a low
visual peak brightness. A good example is SN1987a, which was a faint type IIb supernova
explosion with an absolute visual brightness of -15.5\,mag at its peak
\citep[][and references therein]{arne89}. The mass of the progenitor
star of SN1987a would be in approximate agreement with the progenitor mass proposed for
3C\,58 through the age of the stellar wind bubble \citep[][and references therein]{arne89}.
However, even SN1987a was still more than 1~mag
brighter at its peak than the predicted maximum brightness for 3C\,58. And to the best
of my knowledge there has not been a core-collapse supernova, other than a type IIP,
observed with
a visual peak brightness comparable to the value predicted for 3C\,58. The low absolute
visual brightness at peak is still problematic.
It is possible that one of the main reasons
for this low brightness
is a low explosion energy, and that may explain why we still do not see the radio remnant
of the supernova explosion. Other possible reasons are small ejected $^{56}$Ni masses
\citep{spir09} and a small radius of the progenitor star, which would lead to a longer
period with large diffusion time scales. During this period the expanding star will
lose heat by adiabatic expansion. The latter reason was proposed for the low
luminosity of SN1987a \citep{arne89}.

\subsection{Spin Down Age, Neutron Star Cooling, and Other
Pulsar-Related Arguments}

\subsubsection{Spin Down Age}
\citet{murr02} found a pulsar spin down age, also called the characteristic
age, of $\tau = 5380$~yr for J\,0205+6449, the pulsar inside 3C\,58
($\tau = \frac{P}{2\,\dot{P}}$, $P$: pulsar period).
However, a large discrepancy between the characteristic age and the real age
is not unusual for young pulsars, and \citet{murr02} point out that their result
is not in disagreement with the connection to the supernova of
A.D. 1181. The characteristic age is only a good estimate of the
real age if the intrinsic period $P_0$ is significantly shorter than the current
period. There are many other examples of pulsars
with a characteristic age significantly larger than their actual age. One
extreme example
would be the SNR G11.2$-$0.3 which is linked to a historically observed supernova
explosion from 386~A.D. \citep{step02} and hosts the pulsar J1811$-$1925. This
pulsar is a little over 1600~yr old and has a characteristic age of 23,300 yr
\citep{tori99}.
The radio pulsar wind nebula in
G11.2$-$0.3 is
much fainter than 3C\,58 \citep{koth01}, which probably can be explained by
the much
lower energy loss rate resulting from a weaker magnetic field at birth.

\subsubsection{Neutron Star Cooling}
\citet{slan02} fitted neutron star cooling models to CHANDRA X-ray
observations of J\,0205+6449 and calculated strong upper limits for the
thermal emission originating from the surface of the neutron star in
3C\,58. These fall far below predictions for standard neutron star
cooling. However, this does not necessarily imply a larger age for the pulsar in
3C\,58 but it implies the presence of exotic cooling processes no
matter whether we assume an age of 830~yr or any other age proposed for
3C\,58 in the literature. The minimum age for 3C\,58 if cooling were to follow
standard models is about
25,000~yr \citep[see Fig.~4 in][]{slan02}, far beyond any age that has
been proposed for 3C\,58. Exotic cooling processes must be in operation
in this pulsar, and this result does not favour any of the proposed ages.

\subsubsection{Glitches}

\citet{livi09} found two large spin up glitches in a long-term timing analysis
of the pulsar in 3C\,58. They point out that such glitches are more
typical for pulsars with characteristic ages between 5 and 10~kyr, but suggest that
glitches could be related to the surface temperature of the neutron star
rather than the actual age, and we have just noted (Section 5.2.2) that
the pulsar in 3C\,58 has a very unusually low surface temperature.

\subsubsection{Pulsar Proper Motion}

\citet{gott07} found a circular thermal X-ray shell surrounding the pulsar
inside 3C\,58 from {\it XMM-Newton} observations. However, the pulsar is offset
by $27''\pm 5''$ from the centroid of this shell. If 3C\,58 is the result of the 1181
supernova explosion the pulsar requires a velocity of 500~km\,s$^{-1}$ (300~km\,s$^{-1}$)
perpendicular to the line of sight to reach its current position assuming that it was
born at the very centre of this spherical shell and that its distance
is 3.2~kpc (2.0~kpc). \citet{gott07} argue that if the pulsar in 3C\,58 had the same
transverse velocity as the Crab pulsar ($\sim 70$~km\,s$^{-1}$) it must be much older
than the event from A.D. 1181.

\citet{hobb05} have made a statistical analysis of pulsar proper motions; the
proper motion of pulsar J\,0205+6449 inside 3C\,58 is close to the average
of their sample.
The average for normal pulsars is about 246~km\,s$^{-1}$, for pulsars in SNRs
227~km\,s$^{-1}$, and for pulsars with a characteristic age below $3\times 10^{6}$~yr
310~km\,s$^{-1}$. Hence, the pulsars's predicted proper motion at a distance
of 2~kpc supports its
connection with the supernova of A.D. 1181.

\subsubsection{Other Pulsar Characteristics}

\citet{shea08} say that the low optical efficiency of the emission and the
low radio to X-ray luminosity ratio observed
from the pulsar inside 3C\,58
indicate a low age, when compared with other pulsars
detected at optical wavelength. \citet{abdo09} point out that the ratio
between X-ray and $\gamma$-ray luminosity of that pulsar is more similar
to the middle
aged Vela pulsar than the young Crab pulsar and B1509-58. \citet{kuip10} say
that the spectral properties of PSR J0205+6449 inside 3C\,58
confirm that we are dealing with a young pulsar, and suggest a real
age between those of Vela and PSR B1509-58, favouring its
characteristic age of 5.4 kyr over the connection to SN 1181. However,
these emission characteristics cannot be taken as a means to determine
the age of the pulsar, quite simply because we have no information about
the pulsar's emission characteristics from the time it was born.

\subsection{PWN Evolution and Energetics}

\citet{chev04} discussed evolutionary models of the expansion of the 3C\,58
pulsar wind
nebula into the expanding ejecta of a supernova explosion. He assumed a
power law density distribution for the ejecta into which the PWN is
expanding and that the ejecta is swept up into a thin spherically symmetric
shell. He finds for the ejecta mass:

\begin{equation}
\label{ejecta}
M_{ej}~[{\rm M}_\odot] = 3.5 (7.0)\,\dot{E}_{38}^{0.508}\,E_{51}^{0.492}\,R_{p}^{-2}\,t_3^{2.508}
= 2.0(4.0)\,R_{p}^{-2}\,t_3^{2.508}
\end{equation}

and for the swept-up mass:

\begin{equation}
\label{sweptup}
M_{sw}~[{\rm M}_\odot] = 0.17(0.34)\,\dot{E}_{38}\,R_{p}^{-2}\,t_3^3
= 0.054(0.108)\,R_{p}^{-2}\,t_3^3.
\end{equation}

Here, $\dot{E}_{38}$ is the energy loss rate of the pulsar in
$10^{38}$~erg\,s$^{-1}$, $E_{51}$ the explosion energy in $10^{51}$~erg, $R_{p}$
the PWN radius in pc, and $t_3$ the age of the PWN in $10^3$~yr.

For
the explosion energy I assume the canonical value of $10^{51}$~erg and
for $\dot{E}_{38}$ the average energy loss rate of the pulsar since birth
of $3.2\,10^{37}$~erg\,s$^{-1}$, assuming an age of 830~yr. The lower and
upper limits (in brackets) represent maximum and minimum coupling between
the pulsar wind bubble and the swept up material, respectively \citep{chev04}.
In \citet{chev04} the theoretical results for swept up and ejecta mass and
internal energy are compared with those from observations. All of these results rely
heavily on the assumed distance.

In the following subsections I apply these equations to various aspects of the 3C\,58
situation, and I emphasize the major changes in the interpretation that flow from adopting
the distance 2~kpc as opposed to 3.2~kpc.

\subsubsection{Ejecta Mass}

With equation~\ref{ejecta}, a distance of 3.2 kpc, and an age of 830~yr the ejecta mass
results in 0.1 - 0.2~M$_\odot$ for 3C\,58. This changes
to $M_{ej} = $0.3 - 0.6~M$_\odot$ for a distance of 2~kpc. This is still a
very low value for a core collapse supernova explosion \citep{chev04}.
If the explosion energy is much lower
than the canonical value of $10^{51}$\,erg, which is a possible
explanation for the low peak brightness of the supernova explosion (see section
5.1.1), the ejecta mass would be again much lower.

With deep optical observations \citet{fese08} revealed two distinct kinematic
populations of expanding optical knots within the boundary of the pulsar wind
nebula, the slower population is interpreted as
circumstellar mass-loss material from the 3C\,58 progenitor and the higher
velocity population as expanding supernova ejecta forming a thick shell. This
indicates that the density distribution the PWN
is expanding into is not well represented by a power law and the ejecta is not
entirely swept up into a thin shell. In addition \citet{fese08} reported that the
distribution of the ejecta in their observations is not spherically symmetric but
exhibits a strong bipolar expansion pattern. Therefore the theoretical estimate
of the ejecta mass becomes very uncertain.

\subsubsection{Swept Up Mass}

\citet{bocc01} determined from their X-ray observations a swept-up mass
for the expanding PWN of $M_{sw} = 0.1 d_{3.2}^{2.5}$\,M$_\odot$, assuming a radius $R$ of
$2.5\arcmin$, which results in 0.1\,M$_\odot$ for a distance $d$ of 3.2\,kpc.

We can calculate the ratio between the mass of  observed and theoretical
(Equation~\ref{sweptup}) swept up material to be: $18(9)\,d_{3.2}^5\,t_{830}^{-3}$.
Hence, for an age of 830~yr and the distance of 3.2~kpc the observed swept up mass
is about 10-20 times higher than the theoretical one. Even with the
relatively high uncertainty of this kind of calculation the large discrepancy
suggests a much larger age for 3C\,58. For a distance of 2\,kpc this ratio decreases
to 0.9 - 1.7 so that the two results for the swept-up mass now agree within the
uncertainties.
According to \citet{chev04} the theoretical estimate of the swept up mass
does not rely strongly on the distribution of the ejecta. This is easily
understandable since in principle the equation calculates the mass that could have been
swept up based on the available energy. To extrapolate the total amount of
ejecta from this of course relies heavily on the assumed distribution of
the ejecta material and makes the theoretical estimate of the ejecta mass
very uncertain.

The strong dependence of the theoretically calculated swept-up mass on the
age now actually pretty much negates an age significantly larger than
830 yr. For an age of 1500~yr the ratio between observed and theoretically calculated
swept up mass decreases to 0.1 - 0.2, which despite the uncertainties in theoretical
estimates of that kind makes an age significantly higher than 830~yr very unlikely.

\subsubsection{Energy}

Another argument of \citet{chev04} against an age of 830\,yr is that the
minimum energy $E_{min}$ required to produce the observed synchrotron emission is
higher than the total energy released by the pulsar.
$E_{min}$ is about $10^{48}$\,erg assuming equipartition
between the energy in the magnetic field and the relativistic particles
with $E_{min} \sim d^\frac{17}{7}$ \citep{chev04}. The total energy $E_{tot}$ released
by the pulsar into the nebula, however, can be approximated by
$E_{tot} \approx \dot{E} t = 7\times 10^{47}$ \citep{chev04}
assuming $t = 830$\,yr. $E_{tot}$ is distance independent. We can calculate
this value a little more accurately by integrating the equation for
the evolution of the energy loss rate, $\dot{E}$, of a
pulsar over its lifetime. $\dot{E}$ is defined by \citep[e.g.][]{paci73}:
\begin{equation}
\label{edot}
\dot{E}(t) = \frac{\dot{E_0}}{(1+\frac{t}{\tau_0})^\beta}.
\end{equation}
Here
$\beta$ is related to the nature of the braking torque ($\beta =
\frac{n + 1}{n-1}$, and $n$ is the braking index of the pulsar),
$t$ is the real
age of the pulsar, $\dot{E_0}$ is its energy loss rate
at the time it was born, and $\tau_0$ is its intrinsic characteristic age.
Integrating
equation~\ref{edot} over the lifetime
of 3C\,58 assuming an age of
830~\,yr and a dipolar magnetic field results in
$E_{tot}=8.35\times 10^{47}$\,erg. This
is still significantly lower than $E_{min}$. However, for a distance of
2\,kpc $E_{min}$
decreases to $3.2\times 10^{47}$\,erg, which is now much lower than the total
energy released by
the pulsar, which remains unchanged.

\subsubsection{MHD simulations}

\citet{swal03} performed MHD simulations of a magnetized pulsar wind
interacting with either a uniform interstellar medium or freely expanding
ejecta. For 3C\,58 he found that the observed radio expansion
velocity does not agree with the radius and the age
implied by the connection with the supernova of A.D.~1181. \citet{swal03}
suggest that there are two ways out of this. One would be a much larger age. The
second possibility is that 3C\,58 is being compressed by a reverse shock.
According to \citet{swal03} this
would not only explain the decrease in expansion velocity, but also the radio
flux density increase of 3C\,58 \citep{gree87}, a very low cooling break, and it would decrease
the wind termination shock. The ratio between the wind termination shock
and the PWN radius is remarkably small in 3C\,58 \citep{frai93}.
The question remains, why do we not detect a radio bright shell
corresponding to the blast-wave of the supernova if indeed the PWN is compressed
by a reverse shock? And why don't we detect thermal X-ray emission from
the entire ejecta mass?

\subsection{Expansion Studies}

\begin{figure}[ht]
\begin{center}
 \includegraphics[bb = 55 45 530 480,width=0.48\textwidth,clip]{f5.eps}
\end{center}
\caption{Plot of the optical knots found by \citet{fese08} in a
radius-velocity diagram. The radius is the distance of the knots from
the central pulsar projected to the plane of the sky, assuming a distance
of 2~kpc for 3C\,58. The observed absolute radial velocity is shown on the
y-axis. The dashed line represents the outer edge of this distribution assuming
no deceleration since birth and a connection to SN\,1181. The dotted line represents
the outer edge of this distribution for Sedov expansion \citep{sedo59} and a
connection to SN\,1181. The solid line represents
the best fit outer edge of this distribution.}
\label{exp}
\end{figure}

A supernova explosion is a one-time event.
For any expansion study of a supernova remnant observed features related to this event
could have been decelerated since the time of the explosion but not
accelerated.
This may not be true for synchrotron emitting
filaments related to a pulsar wind nebula, because those have a
continuous source of energy in the central pulsar.
Hence, the interpretation of expansion
studies of synchrotron filaments inside a pulsar wind nebula
is not straightforward.

A comparison of the two major expansion studies in optical
\citep[$\Rightarrow$ age $t \approx 3000$\,yr]{fese08} and radio
\citep[$\Rightarrow t \approx 7000$\,yr]{biet06} already implies
significant deceleration for the synchrotron filaments relative
to the optically observed knots unless the energy released by the
pulsar significantly accelerated all the supernova ejecta. In this
context -- given the age approximations of the two expansion studies --
significantly means that more than 80\,\% of the kinetic energy
of the optical knots must have been provided by the pulsar.
A comparison of the typical explosion energy of $10^{51}$\,erg
released in a supernova with the approximate energy released by the
pulsar since birth of $8.35\times 10^{47}$\,erg (see above) negates
that possibility entirely. Hence, we can assume that the low expansion
velocity observed for the radio synchrotron filaments is the result
of significant deceleration.

Most of the optical filaments are created by material accelerated
through the supernova explosion. Because those filaments could have
been decelerated but not accelerated' a simple averaging of the
expansion velocity of all the filaments would not necessarily lead
to a good age estimate unless the scatter is entirely produced
by uncertainties in the observations or systematic errors, which is
certainly not the case in the careful study of \citet{fese08}. The fastest
filament should be taken, since it presumably shows the lowest
deceleration. Or better in this case would be the ``youngest'' filament
which is the filament that indicates the lowest age by simply dividing
its distance from the pulsar by its velocity. There is quite
a wide spread in velocity in the optical study by \citet{fese08}, even
among those knots that are, according to that study, related to the
explosion and not to the circumstellar material ejected by the progenitor star's
wind. This again
indicates that a lot of the emitting material must have been
decelerated. Therefore the assumption of insignificant deceleration for the
optical knots is invalid as well.
This weakens the case for a large age based on the optical and
radio expansion studies.

In Fig.~\ref{exp} I display a diagram of the optical knots constructed from the
data by
\citet{fese08}. The absolute radial velocity is plotted as a function of radius
from the central pulsar. This radius was calculated assuming that 3C\,58 is 2~kpc
away from the sun. Since the optical knots can have been decelerated but
not accelerated any valid expansion law for the PWN has to encompass all
of these knots.

In Fig.~\ref{exp} I also display several possible expansion laws. For
all a spherical geometry has been adopted. Assuming
a connection to the supernova from 1181~A.D. the resulting best fitting functions
to encompass all expanding knots do not seem to fit very well. Free expansion
leads to a velocity-radius ratio which is too high and for the Sedov solution this
ratio is too low. The best fit expansion law indicated by the solid line represents
either free expansion for an age of about 1600~yr or Sedov expansion for an age of
about 630~yr. Sedov expansion would result in about 0.012 for the ratio of explosion energy
$E_0$ ($10^{51}$~erg) and the ambient density $n_0$ (cm$^{-3}$), which implies either
a very high ambient density of a very low explosion energy.

Free expansion gives only an upper limit for the age in any case since we neglect any kind
of deceleration and the Sedov age can only be a lower limit, since the ejecta mass is
neglected. Therefore the age of 3C\,58 must be between 630 and 1600~yr. We must be in the
transition between free and Sedov expansion where neither law is sufficiently accurate to
describe the expansion.
The large spread in radial velocities for the optical knots indicates a lower rather than
a larger age, since it indicates that a lot of the ejecta must have been
decelerated already. But in principle the
results support any age between 630 and 1600~yr. More accurate calculations are
unfortunately not possible since we lack information about the total mass
of the supernova ejecta and the explosion energy.

\subsection{Radio Flux Evolution}

A very strong argument for a young age is the fact that the radio flux density
of 3C\,58 at 408~MHz increased between 1967 and 1986 at a rate of
$0.32\pm0.13$~\%~yr$^{-1}$ \citep{gree87}.
This is not an insignificant amount considering that the pulsar has been pumping
energy into the nebula since its birth and at a decreasing rate. The radio
synchrotron emission of 3C\,58 depends on the energy
released by the pulsar into the nebula. Increasing radio flux implies that
the energy in magnetic field and relativistic particles
is still being generated by the pulsar faster than it is being radiated away
by the nebula.

Although it is not easy to directly link the radio synchrotron emission to the
total amount of energy in the nebula by
a simple equation, we can find an estimate for the age by
assuming equipartition between relativistic particles and magnetic field. In this
case we find that the flux density $S$ depends on the minimum energy $E_{min}$:
$S \sim E_{min}^\frac{7}{4}$.

\begin{figure}[ht]
\begin{center}
 \includegraphics[bb = 55 50 540 485,width=0.48\textwidth,clip]{f6.eps}
\end{center}
\caption{Plot of the rate of radio flux density increase of 3C\,58
over the last few hundred years predicted for different ages for the PWN.
Negligible energy loss over its lifetime has been assumed. The flux
density increase observed by \citet{gree87} is indicated.}
\label{eev}
\end{figure}

A flux density increase of $0.32 \pm 0.13$~\%~yr$^{-1}$ implies an
energy increase between about 0.11 and 0.26~\%~yr$^{-1}$.
With equation~\ref{edot} this results in an age range between 350 and
810~yr \citep[for the time range discussed in][]{gree87},
assuming no energy loss and a dipolar magnetic field for the
pulsar. Constant or
increasing energy loss would decrease that age range. In Fig.~\ref{eev} I
displayed the rate of radio flux density increase of 3C\,58
over the last few hundred years predicted for different possible ages
for 3C\,58 assuming negligible energy loss. Only a current age of 830~yr would
be located within the error bars observed by \citet{gree87}.
However, even for 830~yr the presence of an additional energy source is more
likely. This additional energy source could be the reverse shock produced by
the interaction of the supernova shock-wave with circumstellar material,
which would be required according to the MHD simulations of \citet{swal03}
(see also Section 5.3.4) to explain the observations. For any
larger age the presence of a significant reverse shock would be essential. This
would imply the remnant has already entered Sedov phase, which implies significant
deceleration, which would indicate that the ages derived from expansion studies
are significant overestimates. A connection to the supernova of 1181~A.D. is
the only possibility to make all the observations consistent.

\begin{table}
\caption{\label{tab:age}Arguments Pro and Contra 830~yr. $\surd$ indicates that this
age is supported and $\times$ that it is not supported by the results of that method.}
\centerline{
\begin{tabular}{lcccc}
\hline \hline
Method & 830~yr & 1500~yr & $\ge 2500$~yr & Section  \\ \hline
A.D. 1181 Gueststar & $\surd$ & $\times$ & $\times$ & 5.1  \\
Peak Brigthness of SN\,1181 & $\times$ & $\surd$ & $\surd$ & 5.1.1  \\
$\tau = 5380$~yr & $\surd$ & $\surd$ & $\surd$ & 5.2.1  \\
NS Cooling & $\surd$ & $\surd$ & $\surd$ & 5.2.2  \\
Glitches & $\surd$ & $\surd$ & $\surd$ & 5.2.1  \\
Proper Motion & $\surd$ & $\surd$ & $\surd$ & 5.2.3  \\
Ejecta Mass & $\surd$ & $\surd$ &  $\surd$ &5.3.1  \\
Swept Up Mass & $\surd$ & $\times$ & $\times$ & 5.3.2  \\
Energetics & $\surd$ & $\surd$ & $\surd$ & 5.3.3  \\
MHD Simulations & $\surd$ & $\surd$ & $\surd$ & 5.3.4  \\
Radio Expansion & $\surd$ & $\surd$ & $\surd$ & 5.4  \\
Optical Expansion & $\surd$ & $\surd$ & $\times$ & 5.4  \\
Radio Flux Evolution & $\surd$ & $\times$ & $\times$ & 5.5  \\ \hline
\end{tabular}}
\end{table}

\section{The Story 3C\,58, a Proposal}

The progenitor star of 3C\,58 must have been a massive star, which was of late O-type during
its Main Sequence phase. Its strong stellar wind carved out the large stellar wind bubble
we find in H~{\sc{i}}. The dynamic age of that bubble of about $3.2\times 10^6$~yr
indicates a mass between 20 and 30~M$_\odot$ for the progenitor star when it was formed,
assuming the bubble was blown by that single star.
In its late evolutionary phase this star must have been a red giant, supergiant star or
even a blue super giant similar to the progenitor of SN1987a,
as indicated by the production of a massive low velocity
wind which is implied by the slowly expanding optical knots found by \citet{fese08}.
When the star exploded it left the neutron star that powers 3C\,58 behind and the
shock-wave of the explosion expanded into the very low density interior of the stellar
wind bubble which was partly filled with high density clouds expelled by the progenitor
star in its wind before it exploded. Consequently, the early expansion of the supernova shock wave
has
to be treated as expanding into a cloudy instead of a homogeneous medium. This explains the
not limb-brightened partly centrally peaked thermal X-ray emission which was observed by
\citet[their Fig.~8, right panel]{gott07} and the lack of non-thermal radio emission
related to the supernova remnant (not the PWN) \citep[e.g.][]{whit91}. The low brightness
of the supernova explosion at the peak of its light curve can be explained by a low
explosion energy, small mass of ejected $^{56}$Ni, a small progenitor radius, or various
other situations.  The interaction
with the circumstellar material is sapping the shock wave of its energy and generates
reverse shocks moving back into the interior. These are decelerating the ejecta, creating
the large velocity spread inside a thick shell. This reverse shock is also decelerating and
compressing the PWN thereby increasing its energy. This
explains the continuous significant increase of radio flux from the PWN.

Although the new distance estimate explains many of the unusual characteristics found for
3C\,58 a few problems still remain:
\begin{itemize}
\item Why did the supernova explosion that created 3C\,58 have such a low peak brightness?
\item Where is the radio remnant of the supernova explosion?
\item Why is the extrapolated mass of the ejecta so low?
\item Why does the neutron star display such a low surface temperature?
\end{itemize}

\section{Summary}

I have derived a new more reliable distance of 2\,kpc to the PWN 3C\,58, by means of
H~{\sc{i}} absorption in combination with a newly determined distance-velocity relation
and by relating this PWN to a nearby H~{\sc{ii}} region and SNR complex of well known
distance. This new distance changes many characteristics of this PWN quite
dramatically. The new results point once again to the historical
connection of 3C\,58, because only an age of 830~yr makes most of the observations
and theoretical calculations consistent.

\begin{acknowledgements}
I would like to thank Tyler Foster for providing me with a distance-velocity
diagram in the direction of 3C\,58 using his modeling technique of the Galactic hydrogen
distribution. I would also like to thank Tom Landecker for careful reading of this
manuscript. The Dominion Radio Astrophysical Observatory is a National
Facility
operated by the National Research Council of Canada. The Canadian Galactic
Plane Survey is a Canadian project with international partners, and is
supported by the Natural Sciences and Engineering Research Council
(NSERC).
\end{acknowledgements}

\bibliographystyle{aa}

Astrophysics of Galaxies (astro-ph.GA); High Energy Astrophysical Phenomena (astro-ph.HE)


Magnetic Fields in Stars: Origin and Impact

Abstract: Various types of magnetic fields occur in stars: small scale fields, large scale fields, and internal toroidal fields. While the latter may be ubiquitous in stars due to differential rotation, small scale fields (spots) may be associated with envelop convection in all low and high mass stars. The stable large scale fields found in only about 10 per cent of intermediate mass and massive stars may be understood as a consequence of dynamical binary interaction, e.g., the merging of two stars in a binary. We relate these ideas to magnetic fields in white dwarfs and neutron stars, and to their role in core-collapse and thermonuclear supernova explosions.

\section{Introduction}
Magnetic fields play a vital role in all stages of stellar evolution.
This is already true during star formation. The magnetic
support in collapsing molecular cloud cores is fundamentally affecting the
fragmentation process (Price \& Bate 2007). Later-on, during the accretion process,
magnetic fields provide the required viscosity to bring in mass and
to remove surplus angular momentum (Donati et al. 2007).

In this paper, we investigate the role of magnetic fields in stars once they are born.
In order to do so, we distinguish various types of magnetic fields.
First, it is useful to distinguish stable from dynamo fields.
As stable fields we consider those which have a decay time of the order of
the stellar life time or more, and which therefore do not need a dynamo
action to continuously replenish them. Braithwaite \& Spruit (2004)
showed that combined toroidal-poloidal magnetic fields can survive in the radiative envelopes
of stars for a long time, which they suggested to exist in magnetic A\,stars and in magnetic
white dwarfs.

Other magnetic field geometries have so far been found unstable, e.g., such fields are
expected to decay on their Alfv\'en time scale (e.g., Tayler 1973). However,
inherently unstable field configurations may be present in stars over long time scales,
if a dynamo process is continuously regenerating the field (Brandenburg \& Subramanian).
It may be expected that this regeneration process leads to some time variability.
The prime example may be the Solar magnetic field, which is produced by a so called
$\alpha\Omega$-dynamo, where the B-field is generated by an interplay between the differential
rotation, which winds up poloidal field and generates toroidal field, and the $\alpha$-effect,
which generates a poloidal field from a toroidal one (R\"udiger et al. 2013).

While only stable and dynamo fields are long-lived and thus accessible to observations,
there is some evidence for intermittent fields playing a role as well (Langer 2012).
In particular during dynamical stellar merger events, which are suspected to lead to stable fields
in the merger product (see below), the fields during the merger event itself are
thought to be significantly stronger than thereafter. This strong intermittent component
may be responsible for a removal of a large fraction of the angular momentum during the merging process.

From the observational perspective, it is also useful to distinguish various types of fields.
For once, there are large scale fields, i.e. fields where the length scale over which
local field maxima occur at the stellar surface is comparable to the size of the star itself.
A classical example is a dipole field, i.e. a field which has only two points of maximum field
strength at the stellar surface, which are located at different sides of the star.
Dipole fields which have their magnetic axis inclined to the axis of rotation are in fact common
amongst intermediate mass and high mass stars, although somewhat more complicated but still
large scale field geometries occur as well.

In contrast to the large scale fields are small scale fields, for which the length scale of significant
field variation is small compared to the stellar radius. An example for this are the
Solar sunspots. The Sun also shows that stellar magnetic fields can have various components,
as the small scale sunspots with field strengths of the order of 1000\,G, and the
global Solar dipole field with a strengths of about 1\,G.

Finally, from the observational perspective, we want to distinguish toroidal magnetic fields
as a third type, since toroidal fields are essentially hidden from direct observations.
Still, they may strongly influence the evolution of stars, and may thus produce indirect evidence
of their existence.

\section{Toroidal fields: ubiquitous?}

\begin{figure}[b]
\begin{center}
 \includegraphics[width=3.4in,angle=-90]{norbert.langer_fig1.ps}
 \caption{
Average core specific angular momentum versus initial mass for
the low and intermediate mass models of Suijs et al. (2008)
and for the 15$\mso$ model of Heger, Woosley \& Spruit (2005)
evolving from the zero age main sequence to their end stages (full drawn lines).
The upper line corresponds
to the initial models. Filled triangles mark the final models
of the non-magnetic sequences, and filled squares  the final models
of the magnetic sequences. The dashed horizontal line indicates the
spectroscopic upper limit on the white dwarf spins obtained by Berger
et al. (2005). Star symbols represent asteroseismic measurements from
ZZ Ceti stars, where smaller symbols correspond to less certain measurements.
The green hatched area is populated by magnetic white dwarfs. The three black open
pentagons correspond to the youngest galactic neutron stars,
while the green pentagon is thought to roughly correspond
to magnetars. See Suijs et al. (2008) for details.
}
\end{center}
\end{figure}

Spruit (2002) has suggested that a dynamo process can operate in differentially rotating
radiative stellar envelopes. The main component of the produced magnetic field is toroidal,
which is thought to counteract the differential rotation by producing a torque which
transports angular momentum against the angular momentum gradient.
While the model of Spruit has been criticized (Zahn, Brun \& Mathis 2007),
the main effect has been confirmed in simplified MHD models (Braithwaite 2006).

While the toroidal fields are not directly observable, there are currently two lines of
observational evidence in their support. First, the nearly rigid rotation in the Sun
beneath the Solar convection zone has been reproduced by Eggenberger et al. (2005) relying on the
Spruit mechanism. Second, Mosser et al. (2012) found through oscillation measurements in
a large sample of pulsating giants that the red giant cores rotate much slower than
expected when only non-magnetic angular momentum transport processes are taken into account.
While they claim that their results agree with the slow observed spins of white dwarfs,
the evolutionary models of Suijs et al. (2007) which include angular momentum transport
through the Spruit mechanism predict indeed white dwarf and neutron star spin periods
which are close to the observed values (Fig.~1). While it can not be excluded that
non-magnetic transport processes like gravity waves (Talon \& Charbonnel 2008) could
also reproduce the observational constraints, the results quoted above may speak in favor of the
Spruit mechanism.

The only prerequisite of the Spruit mechanism is differential rotation.
In case the Spruit mechanism works as expected, we may than conclude that toroidal
magnetic fields are present in {\em all} stars,
perhaps with the exception of stars with strong internal large scale fields, which
perhaps rotate as rigid bodies,
and of fully convective stars --- where the Spruit mechanism may be overpowered by the
predominance of convection.

\section{Low mass stars}

We define low mass stars as such stars which have convective envelopes during core
hydrogen burning. The Sun is a low mass star. Magnetic activity in low mass stars
is investigated since may decades, and many papers in these proceedings give the
status of the current research. We therefore restrict ourselves here to address the magnetic
fields in low mass stars just for comparison to those in more massive stars.

The present conclusion is that {\em all} low mass stars show magnetic fields. I.e., the presence
of a convective envelope is sufficient to develop a field. While rotation is a necessary ingredient
to the $\alpha\Omega$-dynamo, and faster rotators tend to show higher magnetic activity,
even slow rotators as our Sun possess an appreciable magnetic field. It is clear that
these fields lead to an efficient angular momentum loss over the lifetime of these stars,
to the extent that the their spin rate is a function of their age.

\section{Intermediate mass stars}

Intermediate mass stars prove mostly to be non-magnetic. Only about 10\% of the
core hydrogen burning stars in this mass range
show a strong large scale field of more than a few hundred Gauss, while the remaining
90\% appear to have fields which are weaker than about one Gauss.
While the fields in the magnetic fraction of intermediate mass
main sequence stars appear, partly, to have a rather complex morphology
(Donati et al. 2006), their structure appears to be simple compared to the fields of low mass stars
(Donati \& Landstreet 2009).

The magnetic intermediate mass stars are generally slow rotators. I.e., while the field strengths
of low mass stars are larger for faster rotators, the situation is almost the reverse for
intermediate mass stars. While their cores are convective and could produce a magnetic
field deep down (Brun et al. 2005), it appears
unlikely that this field is transported deeply into the radiative envelope
or even to the stellar surface (Charbonneau \& MacGregor 2001, MacGregor \& Cassinelli 2003,
MacDonald \& Mullan 2004).
Since also their envelops are radiative, it appears not to be
possible to explain their magnetic fields through a dynamo process.

As mentioned above, Braithwaite \& Spruit (2004) found stable magnetic field configurations
to be able to exist in these stars. However, while the suggested field geometries are
quite compatible to those observed, this does not allow any conclusion on the
origin of these fields. We return to this question on Sect.~6.

\section{Massive stars}

Evidence is accumulating that massive main sequence stars show both, the small scale
fields produced by convective envelop dynamos, and the large scale stable fields
just as they are observed in the intermediate mass stars.

\subsection{Small-scale fields}

\begin{figure}[b]
\begin{center}
 \includegraphics[width=3.4in, angle=-00]{norbert.langer_fig2.ps}
 \caption{
Evolution of the radial extent of the subsurface helium and iron convective regions (hatched)
as function of time, from the zero age main sequence to roughly the end of core hydrogen burning,
for a 60  $~{M}_\odot$ star (Cantiello et al. 2009). The top of the plot represents the stellar surface. Only the upper
4$~{R}_\odot$ of the star are shown in the plot, while the stellar radius itself increases
during the evolution. The star has a metallicity of Z=0.02, and its effective temperature decreases
from 48 000 K to 18 000 K during the main sequence phase.
}
   \label{fig1}
\end{center}
\end{figure}

Cantiello et al. (2009) pointed out that massive main sequence stars have convective envelops
which may be capable to produce observable magnetic fields. While the convection zones occur
beneath the stellar surface due to opacity peaks produced by iron and helium recombination,
their distance to the surface is so small that magnetic flux tubes can buoyantly float to the surface
in a short time. Their spatial extent is a significant fraction of the stellar radius
(Fig.~2).

There is multiple observational evidence for the existence of these sub-surface convection zones.
First, Cantiello et al. (2009) showed that the predicted dependence of the kinematic signature
of these zones at the stellar surface on stellar mass, surface temperature and metallicity
agrees with the observations of micro-turbulence as determined from spectroscopic measurements of
a large number of O~and early B~stars in the Galaxy and the Magellanic Clouds. Secondly,
stochastically excited pulsations have been measured in several massive main sequence stars
(cf., Belkacem et al. 2009). And thirdly, the velocity field induced by the sub-surface convection
may lead to a clumping of the hot star winds very near to their surface, in agreement with
observations (Cantiello et al. 2009,  Sundqvist \& Owocki 2013).

Finally, and most relevant in the present context, magnetic spots at the surface
of massive main sequence stars as predicted by Cantiello et al. (2009) may provide an
explanation of the discrete absorption components (DACs) in their UV spectra lines.
The DACs phenomenon appears to be best explained by assuming a disturbance of the
radiation driven wind of the star due to hot spots at its surface (Cranmer \& Owocki
1996). As DACs are ubiquitous in O stars (Howarth \& Prinja 1989), the tentative implication is
that so are their small scale magnetic fields.

\subsection{Large-scale fields}

\begin{figure}[b]
\begin{center}
 \includegraphics[width=3.4in,angle=-90]{norbert.langer_fig3.ps}
 \caption{
Hunter diagram for LMC early B~type stars from the VLT-FLAMES Survey of Massive Stars
(symbols), showing projected rotational velocity against their nitrogen surface
abundance.
Single stars are plotted as circles, radial velocity variables as triangles.
A population synthesis simulation based on single star evolution models with rotational
mixing (Brott et al. 2011a) is shown as a density plot in the background (Brott et al. 2011b).
The color coding corresponds to the number of predicted stars per pixel.
The cross in the lower right corner shows the typical error on the observations.
}
\end{center}
\end{figure}

In recent years, the evidence has been growing that concerning large scale stable fields,
the massive stars behave essentially as their intermediate mass counterparts.
The magnetic fraction of OB stars has been determined to be of the order of 10\%
(Grunhut \& Wade 2012),
and the magnetic topologies are similar to those found in intermediate mass stars,
with a predominance of highly inclined magnetic dipoles or low-order multipoles.

Indirect evidence for a magnetic fraction of massive stars comes from the Hunter diagram
of LMC early B-stars (Fig. 3), which shows that 15\% of them are nitrogen-rich slow rotators
(Hunter et al. 2008, Brott et al. 2011b). While this group of stars is not reproduced from
models of rotating single stars (Brott et al. 2011a), a magnetic field would allow to
explain their slow rotation. Morel et al. (2008) identified a similar population in
our Galaxy, and showed that a large fraction of these objects does indeed show a magnetic field.
Further evidence is provided by Dufton et al. (2013), who showed that the velocity distribution
of LMC early B-type stars is bimodal, with $\sim$20\% of them rotating with values below
100\,km/s, very reminiscent of the situation in the A-type stars (Zorec \& Royer 2012).

\section{The origin of stable large-scale fields}

Two ideas are pursued to understand the origin of the stable, large-scale magnetic fields
in intermediate and massive stars. One is that the field inside the main sequence stars
is a relic of the interstellar field present in the molecular cloud at the time when
it formed the stars (e.g., Mestel 2001, Moss 2001).
While it may appear difficult to understand how the memory of the pristine B-field is preserved
when matter is funneled through the MRI-driven accretion disc, this {\em fossil field}
hypothesis has the basic problem to explain why it is $\sim$10\% of the stars that obtain
a large scale field in this way.

Alternatively, it has been postulated that the magnetic field as been acquired by the
star in an earlier phase of its evolution, i.e., not during its formation.
Confusingly, fields according to this hypothesis are also sometimes called {\em fossil}.
Since the event with the largest appeal for field generation is strong close binary
interaction, we speak here of {\em binary-induced fields}.
The general idea that magnetic fields are generated through stellar mergers is
supported by the dearth of close companions
to magnetic main sequence stars at intermediate and high mass (Carrier et al. 2002).

Several people have suggested that strong binary interaction, in particular stellar merger,
can result in the generation of strong, stable magnetic fields. Ferrario et al. (2009)
and Tutukov \& Fedorova (2010)
suggested to explain the magnetic intermediate mass and massive stars through
pre-main sequence mergers, while Tout et al. (2008) argued that the fields in
magnetic white dwarfs may originate from white dwarf mergers.

However, there is evidence that a large fraction of the observed magnetic intermediate
mass and massive stars are remnants of mergers between two main sequence stars.
The latest determination of O~star main sequence binary parameter distributions by
Sana et al. (2012) implies that $8^{+9}_{-4}$\% of all Galactic O~stars are
indeed the product of a merger between two main sequence stars (de Mink et al. 2014).
Furthermore, Glebbeek et al. (2013) predict a nitrogen enrichment in main sequence
merger products which is well compatible with the values found in the magnetic
early B-stars analyzed by Morel et al. (2008; cf. Sect.~5.2).

While indeed also a significant fraction of the intermediate mass Herbig stars
is found to be magnetic (Hubrig et al. 2004, Wade et al. 2007),
those might also be merger products on the pre-main sequence,
possibly induced by circumstellar tides as proposed by Krontreff et al. (2012).
The ratio of nitrogen-rich to nitrogen-normal magnetic massive main sequence stars
may thus give an indication of the ratio of pre-main sequence to main sequence mergers.

A stellar merging process is a most drastic event which induces strong differential rotation
on a timescale close to the dynamical timescale. Binary evolution may produce such a situation
also in some cases where the final merger is avoided in the end, i.e. during common envelope evolution
which leads to an {\em almost-merger}. Perhaps Plaskett's star, which is a very massive
close binary just past its rapid mass transfer phase, and of which the mass gainer has been
found to host a strong magnetic field (Grunhut et al. 2013), and the Polars,
Cataclysmic Variables with magnetic white dwarf companions (Tout et al. 2008),
belong to this category.

It appears in any case unlikely that thermal timescale mass transfer can lead to
strong magnetic fields in the mass gainer, neither accretion during star formation, which occurs
on a similar time scale. The point is that the products are not observed to be magnetic.
Mass transfer in massive binaries is known to lead to Be stars, which may often
be seen as single stars because their companion exploded in a supernova explosion
(de Mink et al. 2013). While perhaps not all Be stars are binary products,
{\em none} of the many analyzed Be stars has been found to be magnetic by Grunhut \& Wade (2012).
Analogously, if accretion during star formation would induce strong stable fields,
all stars should possess such fields, which is obviously not the case.

\section{Evolution}

While in the sections above, we were mostly assessing main sequence star, one may wonder
how the stable magnetic fields survive during the post main sequence evolution
of intermediate mass and massive stars. Theoretical ideas for this are scarce, since
after the main sequence, ever changing parts of the post main sequence stars become
convective. In particular, most of these stars may evolve into red giants and supergiants,
which have convective cores and deep convective envelopes and thus leave little room
for a stable magnetic field. On the other hand, there is no stage where the whole star
would become convective, and it may thus be possible that the stable field of the main
sequence stars is preserved throughout their post-main sequence evolution deep inside the star.

It has indeed been postulated that the magnetic white dwarfs --- again: about 10\% of all
white dwarfs --- are the remnants of magnetic intermediate mass
main sequence stars, as the magnetic fluxes of both are quite comparable.
In the light of the previous section, it may be plausible
to assume that the magnetic white dwarfs have perhaps two components, one evolving from
magnetic main sequence stars, and the other from white dwarf mergers. The latter channel
may lead to an on average larger mass of magnetic white dwarfs, while the former might
stand out by a slower-than-average spin. If magnetic white dwarfs are really merger
products, it may be interesting to note that they will only play a minor role in
Type~Ia supernovae, if any at all. However, magnetic fields could form in double
degenerate mergers, which are thought to provide one channel towards Type~Ia explosions.

While the fraction of neutron stars which have extreme fields, the magnetars, is not
well established, it appears again compatible with an order of magnitude of $\sim 10\%$.
Also the flux freezing argument could apply. If massive main sequence stars
had a B-field of $\sim 10^4$\,G in their core --- which appears plausible as they have surface
fields of up to $\sim 10^3$\,G ---, the resulting B-field in the neutron stars would be of the order of
$10^{14}\,$G, which is two orders of magnitude larger than typical neutron star magnetic field strengths
(Ferrario \& Wickramasinghe 2006).
In contrast to the scenario by Duncan \& Thompson (1992), where the magnetar field forms from an
extremely rapidly rotating collapsing iron core, magnetars as successors of magnetic main sequence stars
would form slowly rotating neutron stars.
This appears not only to be more compatible with the young supernova remnants
surrounding some magnetars (Vink 2008), but also with the dearth of progenitors which can
produce rapidly rotating iron cores in a high metallicity environment as our Galaxy (Yoon et al. 2006),
and would argue against magnetar-powered supernovae (Woosley 2010).
As also proposed by Duncan \& Thompson (1992), the proto-neutron stars with ordinary
spin rates ( corresponding to $j\simeq 10^{14}\,$cm$^2$/s; see Fig.~1) may
well produce the $10^{12}\,$G fields found in most neutron stars.

\section{Outlook}

While all low mass stars appear to have small scale, dynamo-produced surface fields,
the stronger the faster they rotate, only a fraction
of $\sim 10$\% of the intermediate mass stars possess strong B-field, which are large scale
and occur mostly in slowly rotating single stars.
Both types of fields may be combined in the massive stars, the small scale ones in all of the,
the large scale one again in a fraction of about 10\%.
In addition, all stars may contain internal toroidal magnetic fields induced by differential rotation,
which couples their core and envelope spins.

The hypothesis for the formation of the large scale fields
which is consistent with all currently known constraints is that of strong binary
interaction, preferentially via stellar mergers. In contrast to the fossil field hypothesis,
it makes several clear predictions.
Due to the large observational efforts currently underway,
we can expect that this topic will be settled within the next years.

The question of the influence and survival of the large scale fields during the post-main
sequence evolution of intermediate mass and massive stars appears more difficult to answer.
Since red supergiants will likely not allow to assess this question observationally due to
their deep convective envelopes --- which may produce its own field through a dynamo process ---,
it may be interesting to focus on blue supergiants and
Wolf-Rayet stars. If descendants of magnetic main sequence stars evolve into long-lived
blue supergiants, a fair fraction of them might show surface magnetic fields, although
considerably weaker ones if the magnetic flux is conserved. Also some Wolf-Rayet stars may
have magnetic main sequence stars as precursors, but a field detection in these objects
appears difficult due to their strong winds (de la Chevroti\`ere et al. 2013).

Whether the large scale fields survive even until the formation of the compact remnant
remains an open question, although there may be more arguments in favor of this idea
than against it (cf. Sect. 7). However, it remains a challenge to produce solid theoretical predictions
about the survival of the field during the post-main sequence evolution,
as well as to identify observational strategies which would allow to settle the case.

Solar and Stellar Astrophysics (astro-ph.SR)