The heating of the real polar cap surface of radio pulsars by the bombardment of ultra-relativistic charges is studied. The real polar cap is a significantly smaller area within or close by the conventional polar cap which is encircled by the last open field lines of the dipolar field $\vec{B}_d$. It is surrounded by those field lines of the small scale local surface field $\vec{B}_s$ that join the last open field lines of $\vec{B}_d$ in a height of $\sim 10^5$ cm above the cap. As the ratio of radii of the conventional and real polar cap $R_{dip}/R_{pc}\sim 10$, flux conservation requires $B_s/B_d\sim 100$. For rotational periods $P\sim 0.5$ s, $B_s\sim 10^{14}$ G creates a strong electric potential gap that forms the inner accelerating region (IAR) in which charges gain kinetic energies $\sim 3\times 10^{14}$ eV. This sets an upper limit for the energy that back flowing charges can release as heat in the surface layers of the real polar cap. Within the IAR, which is flown through with a dense stream of extremely energetic charges, no stable atmosphere of hydrogen can survive. Therefore, we consider the polar cap as a solidified "naked" surface consisting of fully ionized iron ions. We discuss the physical situation at the real polar cap, calculate its surface temperatures $T_s$ as functions of $B_s$ and $P$, and compare the results with X-ray observations of radio pulsars.

Pulsars are neutron stars, stellar corpses left over after supernova explosions of stars about ten times as massive as our Sun, with densities comparable to the atomic nucleus', spinning with periods from few milliseconds up to few seconds, and endowed with magnetic fields thousands billion times stronger than the Earth's, where particles are accelerated to the relativistic regime producing electromagnetic radiation across the entire spectrum. Although there is a general consensus on the fact that pulsars' radio emission is coherent in nature, whereas the emission from the optical to high-energy $\gamma$-rays is due to incoherent processes, it has not been established yet at which wavelengths the transition occurs, a key information for all emission models of pulsar magnetospheres. Of course, to address this issue covering the spectral region between high-frequency radio waves and the mid-IR is crucial. We used the Atacama Large Millimetre Array (ALMA) to observe the Vela pulsar (PSR\, B0833$-$45), one of the very few observed in radio and from the mid-infrared up to the very high-energy $\gamma$-rays. We detected Vela at frequencies of 97.5, 145, 233, 343.5 GHz and found that its energy density spectrum follows a power-law of spectral index $\alpha = -0.93 \pm 0.16$. The ALMA spectrum corresponds to very high brightness temperatures - from $10^{17}$ to $10^{15}$ K - suggesting that a coherent radiative process still contributes to the mm/sub-mm emission. This is, therefore, the first indication of coherent emission other than radio originating in pulsars. At the same time, we identified an extended structure, preliminarily detected in ground-based observations. We support its interpretation as a counter-jet protruding from the pulsar.

The magnetic and thermal evolution of neutron stars is a very complex process with many nonlinear interactions. For a decent understanding of neutron star physics, these evolutions cannot be considered isolated. A brief overview is presented, which describes the main magnetothermal interactions that determine the fate of both isolated neutron stars and accreting ones. Special attention is devoted to the interplay of thermal and magnetic evolution at the polar cap of radio pulsars. There, a strong meridional temperature gradient is maintained over the lifetime of radio pulsars. It may be strong enough to drive thermoelectric magnetic field creation which perpetuate a toroidal magnetic field around the polar cap rim. Such a local field component may amplify and curve the poloidal surface field at the cap, forming a strong and small scale magnetic field as required for the radio emission of pulsars

We constrain the X-ray properties of the nearby $(360\,{\rm pc})$, old ($5\,{\rm Myr}$) pulsar B1133+16 with $\sim 100\,{\rm ks}$ effective exposure time by {\it XMM-Newton}. The observed pulsar flux in the 0.2-3 keV energy range is $\sim 10^{-14} \, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$, which results in the recording of $\sim 600$ source counts with the EPIC pn and MOS detectors. The X-ray radiation is dominated by nonthermal radiation and is well described by both a single power-law model (PL) and a sum of blackbody and power-law emission (BB+PL). The BB+PL model results in a spectral photon index $\Gamma=2.4^{+0.4}_{-0.3}$ and a nonthermal flux in the 0.2-3 keV energy range of $(7\pm 2) \times 10^{-15}\, {\rm erg \, cm}^{-2} \, {\rm s}^{-1}$. The thermal emission is consistent with the blackbody emission from a small hot spot with a radius of $R_{\rm pc} \approx 14^{+7}_{-5} \, {\rm m}$ and a temperature of $T_{\rm s} = 2.9^{+0.6}_{-0.4} \, {\rm MK}$. Assuming that the hot spot corresponds to the polar cap of the pulsar, we can use the magnetic flux conservation law to estimate the magnetic field at the surface $B_{\rm s} \approx 3.9 \times 10^{14} \, {\rm G}$. The observations are in good agreement with the predictions of the partially screened gap model, which assumes the existence of small-scale surface magnetic field structures in the polar cap region.

A precondition for the radio emission of pulsars is the existence of strong, small-scale magnetic field structures (`magnetic spots') in the polar cap region. Their creation can proceed via crustal Hall drift out of two qualitatively and quantitatively different initial magnetic field configurations: a field confined completely to the crust and another which penetrates the whole star. The aim of this study is to explore whether these magnetic structures in the crust can deform the star sufficiently to make it an observable source of gravitational waves. We model the evolution of these field configurations, which can develop, within $\sim 10^4$ -- $10^5$ yr, magnetic spots with local surface field strengths $\sim 10^{14}$ G maintained over $\gtrsim 10^6$ yr. Deformations caused by the magnetic forces are calculated. We show that, under favourable initial conditions, a star undergoing crustal Hall drift can have ellipticity $\epsilon\sim 10^{-6}$, even with sub-magnetar polar field strengths, after $\sim 10^5$ yr. A pulsar rotating at $\sim 10^2$ Hz with such $\epsilon$ is a promising gravitational-wave source candidate. Since such large deformations can be caused only by a particular magnetic field configuration that penetrates the whole star and whose maximum magnetic energy is concentrated in the outer core region, gravitational wave emission observed from radio pulsars can thus inform us about the internal field structures of young neutron stars.

According to the partially screened gap scenario, an efficient electron-positron pair creation, a general precondition of radio-pulsar activity, relies on the existence of magnetic spots, i.e., local concentrations of strong and small scale magnetic field structures at the surface of neutron stars. They have a strong impact on the surface temperature, which is potentially observable. Here we reinforce the idea that such magnetic spots can be formed by extracting magnetic energy from the toroidal field that resides in deep crustal layers, via Hall drift. We study and discuss the magneto-thermal evolution of qualitatively different neutron star models and initial magnetic field configurations that lead to the creation of magnetic spots. We find that magnetic spots can be created on a timescale of $10^4$ years with magnetic field strengths $\gtrsim 5\times 10^{13}$ G, provided almost the whole magnetic energy is stored in its toroidal component, and that the conductivity in the inner crust is not too large. The lifetime of the magnetic spots is at least $\sim$one million of years, being longer if the initial field permeates both core and crust.

Models of pulsar radio emission that are based on an inner accelerating region require the existence of very strong and small scale surface magnetic field structures at or near the canonical polar cap. The aim of this paper is to identify a mechanism that creates such field structures and maintains them over a pulsar's lifetime. The likely physical process that can create the required 'magnetic spots' is the Hall drift occurring in the crust of a neutron star. It is demonstrated, that the Hall drift can produce small scale strong surface magnetic field anomalies (spots) on timescales of $10^4$ years by means of non-linear interaction between poloidal and toroidal components of the subsurface magnetic field. These anomalies are characterized by strengths of about $10^{14}$ G and curvature radii of field lines of about $10^6$ cm, both of which are fundamental for generation of observable radio emission.

(Abridged) The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index, which is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. By performing axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs, we find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n~3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n<3); in a second stage, the evolution is governed by almost pure Ohmic field decay, and a braking index n>3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of high absolute value and both positive and negative signs. Models with strong (1e14 G) multipolar or toroidal components, even with a weak (~1e12 G) dipolar field are consistent with the observed trend of the timing properties.

The radio pulsar models based on the existence of an inner accelerating gap located above the polar cap rely on the existence of a small scale, strong surface magnetic field $B_s$. This field exceeds the dipolar field $B_d$, responsible for the braking of the pulsar rotation, by at least one order of magnitude. Neither magnetospheric currents nor small scale field components generated during neutron star's birth can provide such field structures in old pulsars. While the former are too weak to create $B_s \gtrsim 5\times 10^{13}$G$\;\gg B_d$, the ohmic decay time of the latter is much shorter than $10^6$ years. We suggest that a large amount of magnetic energy is stored in a toroidal field component that is confined in deeper layers of the crust, where the ohmic decay time exceeds $10^7$ years. This toroidal field may be created by various processes acting early in a neutron star's life. The Hall drift is a non-linear mechanism that, due to the coupling between different components and scales, may be able to create the demanded strong, small scale, magnetic spots. Taking into account both realistic crustal microphysics and a minimal cooling scenario, we show that, in axial symmetry, these field structures are created on a Hall time scale of $10^3$-$10^4$ years. These magnetic spots can be long-lived, thereby fulfilling the pre-conditions for the appearance of the radio pulsar activity. Such magnetic structures created by the Hall drift are not static, and dynamical variations on the Hall time scale are expected in the polar cap region.

The non-linear Hall term present in the induction equation in the electron-magneto-hydrodynamics limit is responsible for the Hall drift of the magnetic field and, in some cases, for the so-called Hall instability. We investigate whether or not the growth rates and eigenfunctions found in the linear analysis are consistent with the results of non-linear numerical simulations. Following the linear analysis of Rheinhardt & Geppert, we study the same cases for which the Hall instability was predicted by solving the non-linear Hall induction equation using a two-dimensional conservative and divergence-free finite difference scheme that overcomes intrinsic difficulties of pseudo-spectral methods and can describe situations with arbitrarily high magnetic Reynolds numbers. We show that unstable modes can grow to the level of the background field without being overwhelmed by the Hall cascade, and cause a complete rearrangement of the field geometry. We confirm both the growth rates and eigenfunctions found in the linearized analysis and hence the instability. In the non-linear regime, after the unstable modes grow to the background level, the naturally selected modes become stable and oscillatory. Later on, the evolution tends to select the modes with the longest possible wavelengths, but this process occurs on the magnetic diffusion timescale. We confirm the existence of the Hall instability. We argue against using the misleading terminology that associates the non-linear Hall term with a turbulent Hall cascade, since small-scale structures are not created everywhere. The field evolves instead in a Burgers-like manner, forms local structures with strong gradients which become shocks in the zero resistivity limit, and Hall waves are launched and propagated through the entire domain.

We study the mutual influence of thermal and magnetic evolution in a neutron star's crust in axial symmetry. Taking into account realistic microphysical inputs, we find the heat released by Joule effect consistent with the circulation of currents in the crust, and we incorporate its effects in 2D cooling calculations. We solve the induction equation numerically using a hybrid method (spectral in angles, but a finite--differences scheme in the radial direction), coupled to the thermal diffusion equation. We present the first long term 2D simulations of the coupled magneto-thermal evolution of neutron stars. This substantially improves previous works in which a very crude approximation in at least one of the parts (thermal or magnetic diffusion) has been adopted. Our results show that the feedback between Joule heating and magnetic diffusion is strong, resulting in a faster dissipation of the stronger fields during the first million years of a NS's life. As a consequence, all neutron stars born with fields larger than a critical value (about 5 10^13 G) reach similar field strengths (approximately 2-3 10^{13} G) at late times. Irrespectively of the initial magnetic field strength, after $10^6$ years the temperature becomes so low that the magnetic diffusion timescale becomes longer than the typical ages of radio--pulsars, thus resulting in apparently no dissipation of the field in old NS. We also confirm the strong correlation between the magnetic field and the surface temperature of relatively young NSs discussed in preliminary works. The effective temperature of models with strong internal toroidal components are systematically higher than those of models with purely poloidal fields, due to the additional energy reservoir stored in the toroidal field that is gradually released as the field dissipates.

We report the X-ray observations of two radio pulsars with drifting subpulses: B0834 + 06 and B0826 - 34 using \xmm\. PSR B0834 + 06 was detected with a total of 70 counts from the three EPIC instruments over 50 ks exposure time. Its spectrum was best described as that of a blackbody (BB) with temperature $T_s=(2.0^{+2.0}_{-0.9}) \times 10^6$ K and bolometric luminosity of $L_b=(8.6^{+14.2}_{-4.4}) \times 10^{28}$ erg s$^{-1}$. As it is typical in pulsars with BB thermal components in their X-ray spectra, the hot spot surface area is much smaller than that of the canonical polar cap, implying a non-dipolar surface magnetic field much stronger than the dipolar component derived from the pulsar spin-down (in this case about 50 times smaller and stronger, respectively). The second pulsar PSR B0826 - 34 was not detected over 50 ks exposure time, giving an upper limit for the bolometric luminosity $L_b \leq 1.4 \times 10^{29}$ erg s$^{-1}$. We use these data as well as the radio emission data concerned with drifting subpulses to test the Partially Screened Gap (PSG) model of the inner accelerator in pulsars.

We study the non--linear evolution of magnetic fields in neutron star crusts with special attention to the influence of the Hall drift. Our goal is to understand the conditions for fast dissipation due to the Hall term in the induction equation. We study the interplay of Ohmic dissipation and Hall drift in order to find a timescale for the overall crustal field decay. We solve numerically the Hall induction equation by means of a hybrid method (spectral in angles but finite differences in the radial coordinate). The microphysical input consists of the most modern available crustal equation of state, composition and electrical conductivities. We present the first long term simulations of the non--linear magnetic field evolution in realistic neutron star crusts with a stratified electron number density and temperature dependent conductivity. We show that Hall drift influenced Ohmic dissipation takes place in neutron star crusts on a timescale of 1 Myr. When the initial magnetic field has magnetar strength, the fast Hall drift results in an initial rapid dissipation stage that lasts 10-50 kyr. The interplay of the Hall drift with the temporal variation and spatial gradient of conductivity tends to favor the displacement of toroidal fields toward the inner crust, where stable configurations can last for 1 Myr. We show that the thermally emitting isolated neutron stars, as the Magnificent Seven, are very likely descendants of neutron stars born as magnetars.

We show the existence of a strong trend between neutron star surface temperature and the dipolar component of the magnetic field extending through three orders of field magnitude, a range that includes magnetars, radio-quiet isolated neutron stars, and many ordinary radio pulsars. We suggest that this trend can be explained by the decay of currents in the crust over a time scale of few Myr. We estimate the minimum temperature that a NS with a given magnetic field can reach in this interpretation.

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

During the life of isolated neutron stars (NSs) their magnetic field passes through a variety of evolutionary phases. Depending on its strength and structure and on the physical state of the NS (e.g. cooling, rotation), the field looks qualitatively and quantitatively different after each of these phases. Three of them, the phase of MHD instabilities immediately after NS's birth, the phase of fallback which may take place hours to months after NS's birth, and the phase when strong temperature gradients may drive thermoelectric instabilities, are concentrated in a period lasting from the end of the proto--NS phase until 100, perhaps 1000 years, when the NS has become almost isothermal. The further evolution of the magnetic field proceeds in general inconspicuous since the star is in isolation. However, as soon as the product of Larmor frequency and electron relaxation time, the so-called magnetization parameter, locally and/or temporally considerably exceeds unity, phases, also unstable ones, of dramatic changes of the field structure and magnitude can appear. An overview is given about that field evolution phases, the outcome of which makes a qualitative decision regarding the further evolution of the magnetic field and its host NS.

(abridged) We report on XMM-Newton observations of SN1987A in the Large Magellanic Cloud. The large collecting area telescopes together with the European Photon Imaging Cameras (EPIC) provide X-ray spectra with unprecedented statistical quality and make it possible to investigate the spectral evolution during the brightening observed since the discovery in X-rays. High resolution spectra from the Reflection Grating Spectrometers yield a complementary view and allow us to perform more detailed investigations of prominent emission lines. The X-ray spectra were modeled with two-temperature emission components from a hot plasma in collisional ionization equilibrium and in non-equilibrium (NEI). We find a temperature for the equilibrium component of 0.24+/-0.02 keV in January 2000 and April 2001 which increased to 0.30+/-0.02 keV in May 2003 and also an indication for a temperature increase in the hot NEI component from \~2 keV to ~3 keV. Emission line ratios inferred from the RGS spectra suggest temperatures as low as 100 eV and an increase in the ionization state of oxygen and neon consistent with the observed temperature increases. The fast readout of the EPIC-pn instrument yields X-ray fluxes free of CCD pile-up effects which we used to normalize pile-up corrections for the published Chandra fluxes. The corrected X-ray light curve of SN1987A in the 0.5-2.0 keV energy band is best represented by a linear increase up to about day 4000 after the explosion and an exponential rise afterwards until the last published Chandra observation on day 6716. Modeling the light curve by emission from the inner ring which is approximated by a circular torus a central density n_H = 1.15x10^4 cm^-3 is found. In this model the forward shock has just passed the center of the torus.

We study the stability/establishment of dipolar magnetostatic equilibrium configurations in new--born neutron stars (NSs) in dependence on the rotational velocity $\Omega$ and on the initial angle $\alpha$ between rotation and magnetic axis. The NS is modeled as a sphere of a highly magnetized ($B \sim 10^{15}$G) incompressible fluid of uniform density which rotates rigidly. For the initial dipolar background magnetic field, which defines the magnetic axis, two different configurations are assumed. We solve the 3D non--linear MHD equations by use of a spectral code. The problem in dimensionless form is completely defined by the initial field strength (for a fixed field geometry), the magnetic Prandtl number $\Pm$, and the normalized rotation rate. The evolution of the magnetic and velocity fields is considered for initial magnetic field strengths characterized by the ratio of ohmic diffusion and initial \Alf{} travel times $\ttOhm/\ttAO \approx 1000$, for $\Pm = 0.1, 1, 10$, and the ratio of rotation period and initial \Alf{} travel time, $P/\ttAO = 0.012, 0.12, 1.2, 12$. We find hints for the existence of a unique stable dipolar magnetostatic configuration for any specific $\alpha$, independent of the initial field geometry. Comparing NSs possessing the same field structure at the end of their proto--NS phase, it turns out that sufficiently fast rotating NSs ($P\la6 $ms) with $\alpha \la 45^0$ retain their magnetar field, while the others lose almost all of their initial magnetic energy by transferring it into magnetic and kinetic energy of relatively small--scaled fields and continue their life as radio pulsars with a dipolar surface field of $10^{12...13}$G.

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

The cooling of a compact star depends very sensitively on the state of dense matter at supranuclear densities, which essentially controls the neutrino emission, as well as on the structure of the stellar outer layers which control the photon emission. Open issues concern the hyperon population, the presence of meson condensates, superfluidity and superconductivity, and the transition of confined hadronic matter to quark matter. This paper describes these issues and presents cooling calculations based on a broad collection of equations of state for neutron star matter and strange matter. These results are tested against the body of observed cooling data.

We study convective motions taken from hydrodynamic simulations of rotating proto--neutron stars (PNSs) with respect to their ability to excite a dynamo instability which may be responsible for the giant neutron star magnetic fields. Since it is impossible to simulate the magnetic field evolution employing the actual magnetic Reynolds numbers ($\Rm$) resulting from the hydrodynamic simulations, (smallest) critical $\Rm$s and the corresponding field geometries are derived on the kinematic level by rescaling the velocity amplitudes. It turns out that the actual values of $\Rm$ are by many orders of magnitude larger than the critical values found. A dynamo might therefore start to act vigorously very soon after the onset of convection. But as in general dynamo growth rates are non--monotonous functions of $\Rm$ the later fate of the magnetic field is uncertain. Hence, no reliable statements on the existence and efficiency of PNS dynamos can be drawn without considering the interplay of magnetic field and convection from the beginning. Likewise, in so far as convection inside the PNS is regarded to be essential in re--launching the supernova explosion, a revision of its role in this respect could turn out to be necessary.

In the paper we comment on (R\"udiger & Shalybkov, Phys. Rev. E. 69, 016303 (2004) (RS)), the instability of the Taylor--Couette flow interacting with a homogeneous background field subject to Hall effect is studied. We correct a falsely generalizing interpretation of results presented there which could be taken to disprove the existence of the Hall--drift induced magnetic instability described in Rheinhardt and Geppert, Phys. Rev. Lett. 88, 101103. It is shown that in contrast to what is suggested by RS, no additional shear flow is necessary to enable such an instability with a non--potential magnetic background field, whereas for a curl--free one it is. In the latter case, the instabilities found in RS in situations where neither a hydrodynamic nor a magneto--rotational instability exists are demonstrated to be most likely magnetic instead of magnetohydrodynamic. Further, some minor inaccuracies are clarified.

We report the results of a deep near-infrared survey with VLT/ISAAC of the environment of the supernova remnant N49 in the Large Magellanic Cloud, which contains the soft gamma-ray repeater SGR 0526-66. Two of the four confirmed SGRs are potentially associated with compact stellar clusters. We thus searched for a similar association of SGR 0526-66, and imaged a young stellar cluster at a projected distance of ~30 pc from the SGR. This constitutes the third cluster-SGR link, and lends support to scenarios in which SGR progenitors originate in young, dusty clusters. If confirmed, the cluster-SGR association constrains the age and thus the initial mass of SGR progenitors.

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

In former papers we showed that during the decay of a neutron star's magnetic field under the influence of the Hall--drift, an unstable rise of small--scale field structures at the expense of the large--scale background field may happen. This linear stability analysis was based on the assumption of a uniform density throughout the neutron star crust, whereas in reality the density and all transport coefficients vary by many orders of magnitude. Here, we extend the investigation of the Hall--drift induced instability by considering realistic profiles of density and chemical composition, as well as background fields with more justified radial profiles. Two neutron star models are considered differing primarily in the assumption on the core matter equation of state. For their cooling history and radial profiles of density and composition we use known results to infer the conductivity profiles. These were fed into linear calculations of a dipolar field decay starting from various initial configurations. At different stages of the decay, snapshots of the magnetic fields at the equator were taken to yield background field profiles for the stability analysis. The main result is that the Hall instability may really occur in neutron star crusts. Characteristic growth times are in the order of \lesssim 10^4 ... 10^6 yrs depending on cooling age and background field strength. The influence of the equation of state and of the initial field configuration is discussed.

There is growing evidence, based on both X-ray and radio observations of isolated neutron stars, that besides the large--scale (dipolar) magnetic field, which determines the pulsar spin--down behaviour, small--scale poloidal field components are present, which have surface strengths one to two orders of magnitude larger than the dipolar component. We argue in this paper that the Hall--effect can be an efficient process in producing such small--scale field structures just above the neutron star surface. It is shown that due to a Hall--drift induced instability, poloidal magnetic field structures can be generated from strong subsurface toroidal fields, which are the result of either a dynamo or a thermoelectric instability acting at early times of a neutron star's life. The geometrical structure of these small--scale surface anomalies of the magnetic field resembles that of some types of ``star--spots''. The magnetic field strength and the length--scales are comparable with values that can be derived from various observations.

The classical vacuum gap model of Ruderman & Sutherland, in which spark-associated subbeams of subpulse emission circulate around the magnetic axis due to the EB drift, provides a natural and plausible physical mechanism of the subpulse drift phenomenon. Recent progress in the analysis of drifting subpulses in pulsars has provided a strong support to this model by revealing a number of subbeams circulating around the magnetic axis in a manner compatible with theoretical predictions. However, a more detailed analysis revealed that the circulation speed in a pure vacuum gap is too high when compared with observations. Moreover, some pulsars demonstrate significant time variations of the drift rate, including a change of the apparent drift direction, which is obviously inconsistent with the EB drift scenario in a pure vacuum gap. We resolved these discrepancies by considering a partial flow of iron ions from the positively charged polar cap, coexisting with the production of outflowing electron-positron plasmas. By fitting the observationally deduced drift-rates to the theoretical values, we managed to estimate polar cap surface temperatures in a number of pulsars. The estimated surface temperatures correspond to a small charge depletion of the order of a few percent of the corotational charge density. We also argue that if the thermionic electron outflow from the surface of a negatively charged polar cap is slightly below the Goldreich-Julian density, then the resulting small charge depletion will have similar consequences as in the case of the ions outflow. We thus believe that the sparking discharge of a partially shielded acceleration potential drop occurs in all pulsars, with both positively (``pulsars'') and negatively (``anti-pulsars'') charged polar caps.

There exists both theoretical and observational evidence that the magnetic field decay in neutron stars may proceed in a pronounced non--linear way during a certain episode of the neutron star's life. In the presence of a strong magnetic field the Hall--drift dominates the field evolution in the crust and/or the superfluid core of neutron stars. Analysing observations of $P$ and $\dot{P}$ for sufficiently old isolated pulsars we gain strong hints for a significantly non--linear magnetic field decay. Under certain conditions with respect to the geometry and strength of a large--scale magnetic background field an instability is shown to occur which rapidly raises small--scale magnetic field modes. Their growth rates increase with the background field strength and may reach $\sim 10^4$ times the ohmic decay rate. Consequences for the rotational and thermal evolution as well as for the cracking of the crust of neutron stars are discussed.

In the presence of a strong magnetic field and under conditions as realized in the crust and the superfluid core of neutron stars the Hall--drift dominates the field evolution. We show by a linear analysis that for a sufficiently strong large--scale background field depending at least quadratically on position in a plane conducting slab an instability occurs which rapidly generates small--scale fields. Their growth rates depend on the choice of the boundary conditions, increase with the background field strength and may reach $10^3$ times the ohmic decay rate. The effect of that instability on the rotational and thermal evolution of neutron stars is discussed.

We consider the expulsion of proton fluxoids along neutron vortices from the superfluid/superconductive core of neutron star with weak ($B<10^{10}$ G) magnetic field. The velocity of fluxoids is calculated from the balance of buoyancy, drag and crustal forces. We show, that the proton fluxoids can leave the superfluid core sliding {\it along} the neutron vortices on a timescale of about $10^7$ years. An alternative possibility is that fluxoids are aligned with the vortices on the same timescale. As the result, non--aligned surface magnetic fields of millisecond pulsars can be sustained for $\ga 10^9$ years only in case of a comparable dissipation timescale of the currents in the neutron star crust. This defines upper limits of the impurity concentration in the neutron star crust: $Q \la 0.1$ if a stiff equation of state determines the density profile.

We apply the model of flux expulsion from the superfluid and superconductive core of a neutron star, developed by Konenkov & Geppert (2000), both to neutron star models based on different equations of state and to different initial magnetic field structures. When initially the core and the surface magnetic field are of the same order of magnitude, the rate of flux expulsion from the core is almost independent of the equation of state, and the evolution of the surface field decouples from the core field evolution with increasing stiffness. When the surface field is initially much stronger than the core field, the magnetic and rotational evolution resembles to those of a neutron star with a purely crustal field configuration; the only difference is the occurence of a residual field. In case of an initially submerged field significant differences from the standard evolution occur only during the early period of neutron star's life, until the field has been rediffused to the surface. The reminder of the episode of submergence is a correlation of the residual field strength with the submergence depth of the initial field. We discuss the effect of the rediffusion of the magnetic field on to the difference between the real and the active age of young pulsars and on their braking indices. Finally, we estimate the shear stresses built up by the moving fluxoids at the crust--core interface and show that preferentially in neutron stars with a soft equation of state these stresses may cause crust cracking.

We explore, using an exact cooling code, the thermal evolution of a neutron star undergoing episodes of intense accretion, alternated by long periods of quiescence (e.g. Soft X-Ray Transients). We find that the soft component of the quiescent luminosity of Aql X-1, 4U 1608-522 and of SAX J1808.4-3658 can be understood as thermal emission from a cooling neutron star with negligible neutrino emission. In the case of Cen X-4 strong neutrino emission from the inner core is necessary to explain the observation: this may indicate that the neutron star of Cen X-4 is heavier than 1.4 Msun. This study opens the possibility of using the quiescent emission of Soft X-Ray Transients as a tool for probing the core superfluidity in relation to the mass of the neutron star.

We explore the thermal evolution of a neutron star undergoing episodes of intense accretion, separated by long periods of quiescence. By using an exact cooling code we follow in detail the flow of heat in the star due to the time-dependent accretion-induced heating from pycno-nuclear reactions in the stellar crust, to the surface photon emission, and the neutrino cooling. These models allow us to study the neutron stars of the Soft X-Ray Transients. In agreement with Brown, Bildsten and Rutledge (1998) we conclude that the soft component of the quiescent luminosity of Aql X-1, 4U 1608-522, and of the recently discovered SAX J1808.4, can be understood as thermal emission from a cooling neutron star with negligible neutrino emission. However, we show that in the case of Cen X-4, despite its long recurrence time, strong neutrino emission from the neutron star inner core is necessary to understand the observed low ratio of quiescent to outburst luminosity. This result implies that the neutron star in Cen X-4 is heavier than the one in the other systems and the pairing critical temperature in its center must be low enough (well below 10^9 K) to avoid a strong suppression of the neutrino emission.

We investigate the thermal, magnetic and rotational evolution of isolated neutron stars assuming that the dipolar magnetic field is confined to the crust. Our treatment, for the first time, uses a fully general relativistic formalism not only for the thermal but also for the magnetic part, and includes partial general relativistic effects in the rotational part. Due to the fact that the combined evolution depends crucially upon the compactness of the star, three different equations of state have been employed in the calculations. In the absence of general relativistic effects, while upon increasing compactness a decrease of the crust thickness takes place leading into an accelerating field decay, the inclusion of general relativistic effects intend to ``decelerate this acceleration''. As a consequence we find that within the crustal field hypothesis, a given equation of state is compatible with the observed periods $P$ and period derivative $\dot P$ provided the initial field strength and current location as well as the magnitude of the impurity content are constrained appropriately. Finally, we access the flexibility of the soft, medium and stiff classes of equations of state as candidates in describing the state of the matter in the neutron star interiors. The comparison of our model calculations with observations, together with the consideration of independent information about neutron star evolution, suggests that a not too soft equation of state describes neutron star interiors and its cooling proceeds along the `standard' scenario.

An analysis of the role of general relativistic effects on the decay of neutron star's magnetic field is presented. At first, a generalized induction equation on an arbitrary static background geometry has been derived and, secondly, by a combination of analytical and numerical techniques, a comparison of the time scales for the decay of an initial dipole magnetic field in flat and curved spacetime is discussed. For the case of very simple neutron star models, rotation not accounted for and in the absence of cooling effects, we find that the inclusion of general relativistic effects result, on the average, in an enlargement of the decay time of the field in comparison to the flat spacetime case. Via numerical techniques we show that, the enlargement factor depends upon the dimensionless compactness ratio ${\epsilon}={2GM \over c^{2}R}$, and for ${\epsilon}$ in the range $(0.3~,~0.5)$, corresponding to compactness ratio of realistic neutron star models, this factor is between 1.2 to 1.3. The present analysis shows that general relativistic effects on the magnetic field decay ought to be examined more carefully than hitherto. A brief discussion of our findings on the impact of neutron stars physics is also presented.

We confront theoretical models for the rotational, magnetic, and thermal evolution of an ultra-magnetized neutron star, or magnetar, with available data on the Anomalous X-ray Pulsars (AXPs). We argue that, if the AXPs are interpreted as magnetars, their clustering of spin periods between 6 and 12 seconds, observed at present in this class of objects, their period derivatives, their thermal X-ray luminosities, and the association of two of them with young supernova remnants, can only be understood globally if the magnetic field in magnetars decays significantly on a time scale of the order of 10^4 years.

The paper deals with a simple spherical mean-field dynamo model of $\alpha^2$-type in the case of high electrical conductivity. A spherically symmetric distribution of turbulent motions is assumed inside a spherical fluid body surrounded by free space, and some complete form of the corresponding mean electromotive force is taken into account. For a turbulence lacking reflectional symmetry finite growth rates of magnetic fields prove to be possible even in the limit of perfect conductivity, that is, the model corresponds to a fast dynamo. In accordance with the theorem by Bondi and Gold the dynamo-generated field is completely confined in the fluid body.

We investigate the statistical evolution of magnetic neutron stars recycled in Low Mass Binary (LMB) systems, simulating synthetic populations. Irrespective to the details of the physical models, we find to be significant the fraction of neutron stars spinning close to their mass shedding limit relative to the millisecond pulsar population. The simulated neutron stars show a tail in their period distribution at periods shorter than 1.558 ms, the minimum detected so far. Crustal magnetic field decay models predict also the existence of massive rapidly spinning neutron stars with very low magnetic moment.

Soft Gamma-ray Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) are interpreted as young highly magnetized neutron stars (NSs). Their X-ray luminosity in quiescence, exceeding 10^{35} erg s^{-1} cannot be explained as due to cooling of a highly magnetized NS, but requires as an extra heat source the decay of its magnetic field (MF). We study numerically the coupled evolution of the MF, temperature and spin period under the assumption that the currents maintaining the field are confined in the crust of the star. The decay of the field depends on the field strength itself (Hall-drift), on the temperature and injects heat into the star, but is controlled by neutrino emission. Finally we consider the spin down from magnetic dipole braking with this decaying field to track the long term evolution. We find reasonable initial conditions for the MF strength and structure to explain their current observational values both of their rotational period, its time derivative and the X-ray luminosity of AXPs and SGRs.the X-ray luminosity of AXPs and SGRs.

We consider the expulsion of the magnetic field from the super-conducting core of a neutron star and its subsequent decay in the crust. Particular attention is paid to a strong feedback of the distortion of magnetic field lines in the crust on the expulsion of the flux from the core. This causes a considerable delay of the core flux expulsion if the initial field strength is larger than 10^{11} G. It is shown that the hypothesis on the magnetic field expulsion induced by the neutron star spin-down is adequate only for a relatively weak initial magnetic field $B \approx 10^{11}$ G. The expulsion time-scale depends not only on the conductivity of the crust, but also on the initial magnetic field strength itself. Our model of the field evolution naturally explains the existence of the residual magnetic field of neutron stars. Its strength is correlated with the impurity concentration in neutron star crusts and anti-correlated with the initial field strengths.

We investigate the statistical evolution of magnetic neutron stars, recycled in binary systems, simulating synthetic populations. To bracket uncertainties, we consider a soft (FP) and a stiff (PS) equation of state (EoS) for nuclear matter and explore the hypothesis that the magnetic field is confined in the stellar crust. We follow the magneto-rotational evolution within a simple recycling scenario, including the possibility of magnetospheric propeller. We find the presence of a tail in the period distribution of the synthetic populations at periods shorter than 1.558 ms, the minimum detected so far. For the soft EoS the recycling gives rise to a spin distribution which is increasing monotonically toward short periods and a clear ``barrier'' forms at the minimum period for the onset of mass shedding. For the stiff EoS the distribution is flatter displaying a broad maximum about 2-4 ms. The estimated fraction of neutron stars spinning close to their shedding limit over the millisecond pulsar population is found to be significant. Crustal magnetic field decay models predict the existence of massive (M>1.4 M_sun) rapidly spinning neutron stars with very low magnetic moment.

The evolution of neutron stars in close binary systems with a low-mass companion is considered assuming the magnetic field to be confined within the solid crust. We adopt the standard scenario of the evolution in a close binary system in accordance with which the neutron star passes throughout four evolutionary phases ("isolated pulsar" -- "propeller" -- accretion from the wind of a companion -- accretion due to Roche-lobe overflow). Calculations have been performed for a great variety of parameters characterizing the properties both of the neutron star and low-mass companion. We find that neutron stars with more or less standard magnetic field and spin period being processed in low-mass binaries can evolve to low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long as $10^{10}$ yr, the neutron star can maintain a relatively strong magnetic field to the end of the accretion phase. The considered model can well account for the origin of millisecond pulsars.