The magnetospheres around neutron stars should be very particular because of their strong magnetic field and rapid rotation. A study of the pulsar magnetospheres is of a crucial importance since it is the key issue to understand how energy inflow to the exterior is produced. In this paper, I discuss magnetohydrodynamic processes in the pulsar magnetosphere. I consider in detail the properties of magnetohydrodynamic waves that can exist in the magnetosphere and their instabilities.These instabilities lead to formation of magnetic structures and can be responsible for a short-term variability ofthe pulsar emission.

A new magnetic field-generation mechanism based on the Nernst effect is considered in hot massive stars. This mechanism can operate in the upper atmospheres of O and B stars where departures from the LTE form a region with the inverse temperature gradient.

We consider diffusion caused by a combined influence of the electric current and the Hall effect, and argue that such diffusion can form inhomogeneities of the chemical composition in plasma. The considered mechanism can be responsible for a formation of element spots in laboratory and astrophysical plasmas. This current-driven diffusion can be accompanied by propagation of a particular type waves in which the impurity number density oscillate alone. These compositional waves exist if the magnetic pressure in plasma is much greater than the gas pressure,

The stars of the middle main sequence often have relatively quiescent outer layers and spot-like chemical structures may develope in their atmospheres. Recent observations show that abundance peculiarities can change as stars evolve on the main sequence and the timescale of these changes lies in a wide range from million years to months. These observations imply that, perhaps, our understanding of diffusion processes at work in magnetic stars is incomplete and a more detailed analysis of these processes is required. In the present paper, we consider diffusion caused by a combined influence of the electric current and the Hall effect.Such diffusion has a number of very particular properties and, generally, can change the surface chemistry of stars in combination with other diffusion processes. For instance, current-driven diffusion is accompanied by a propagation of the special type of waves in which only the impurity number density oscillates. Propagation of such waves changes the shape and size of spots as well as chemical abundances within them. The period of compositional waves depends on the parameters of plasma (magnetic field, electric current, temperature, etc.) and can be different for the waves of different elements. Compositional waves exist in the regions where the magnetic pressure is greater than the gas pressure. These waves can be the reason of variations with different timescales in the abundance peculiarities of magnetic stars.

The stars of the middle main sequence often have spot-like chemical structures at their surfaces. We consider the diffusion process caused by electric currents that can lead to the formation of such chemical spots. Diffusion is considered using the partial momentum equations derived by the Chapman-Enskog method. We argue that diffusion caused by electric currents can substantially change the surface chemistry of stars and form spotted chemical structures even in a relatively weak magnetic field. The considered mechanism can be responsible for a formation of element spots in HgMn- and Ap-stars.

We consider combined rotational, magnetic, and thermal evolution of the neutron star during the accretion phase in a binary system. A rapid accretion-driven decay of the magnetic field decreases substantially the efficiency of angular momentum transfer. As a result, the neutron star cannot spin up to periods shorter than some limiting value even if accretion is very long and accretion rate is high. The proposed mechanism can explain a discrepancy between the shortest detected period and minimal possible spin period of neutron stars.

The recent discovery of long-period, low magnetic field pulsars in low-mass X-ray binaries (LMXBs) represents a challenge for the standard evolutionary scenario. These pulsars have a magnetic field strength comparable to that of millisecond pulsars ($\sim 10^8 - 10^9$ G), but their period is at least an order of magnitude longer. We discuss the origin of this new class of pulsars within the standard picture of LMXBs formation and apply our results to the case of IGR J17480-2446. The magnetothermal evolution of the binary system is studied numerically by taking into account the effect of different accretion rates during the Roche-lobe overflow in the framework of the minimal cooling scenario. We show that, in addition to standard millisecond pulsars, long-period low magnetic field pulsars should also be expected as a possible outcome of the binary evolution, depending on the strength of the accretion rate during the Roche-lobe overflow. In particular, we argue that IGR J17480-2446 belongs to this new class of objects.

Magnetohydrodynamic (MHD) instabilities can play an important role in the dynamics of the pulsar magnetosphere and can be responsible for the formation of various structures. We consider the instability caused by a gradient of the magnetic pressure which can occur in a non-neutral magnetospheric plasma. We arguethat the pulsar magnetospheres are always unstable. The unstable disturbances have a form of filaments directed along the magnetic field lines with plasma motions being almost parallel (or anti-parallel) to the magnetic field. The growth rate of instability is high and can reach a fraction of $ck$, where $k$ is the wavevector of unstable disturbances. The instability can be responsible for fluctuations of plasma and the short-term variability of pulsar emission.

The thermomagnetic evolution of the young neutron star in Cassiopea A is studied by considering fast neutrino emission processes. In particular, we consider neutron star models obtained from the equation of state computed in the framework of the Brueckner-Bethe-Goldstone many-body theory and variational methods, and models obtained with the Akmal-Pandharipande-Ravenhall equation of state. It is shown that it is possible to explain a fast cooling regime as the one observed in the neutron star in Cassiopea A if the Joule heating produced by dissipation of the small-scale magnetic field in the crust is taken into account. We thus argue that it is difficult to put severe constraints on the superfluid gap if the Joule heating is considered.

The study of the magnetic field in stellar radiation zones is an important topic in modern astrophysics because the magnetic field can play an important role in several transport phenomena such as mixing and angular momentum transport. We consider the influence of rotation on stability of a predominantly toroidal magnetic field in the radiation zone. We find that the effect of rotation on the stability depends on the magnetic configuration of the basic state. If the toroidal field increases sufficiently rapidly with the spherical radius, the instability cannot be suppressed entirely even by a very fast rotation although the strength of the instability can be significantly reduced. On the other hand, if the field increases slowly enough with the radius or decreases, the instability has a threshold and can be completely suppressed in rapidly rotating stars. We find that in the regions where the instability is entirely suppressed a particular type of magnetohydrodynamic waves may exist which are marginally stable.

The magnetic field in stellar radiation zones can play an important role in phenomena such as mixing, angular momentum transport, etc. We study the effect of rotation on the stability of a predominantly toroidal magnetic field in the radiation zone. In particular we considered the stability in spherical geometry by means of a linear analysis in the Boussinesq approximation. It is found that the effect of rotation on the stability depends on a magnetic configuration. If the toroidal field increases with the spherical radius, the instability cannot be suppressed entirely even by a very fast rotation. Rotation can only decrease the growth rate of instability. If the field decreases with the radius, the instability has a threshold and can be completey suppressed.

Magnetohydrodynamic (MHD) instabilities can play an important role in the structure and dynamics of the pulsar magnetosphere. We consider the instabilitycaused by differential rotation that is suggested by many theoretical models. Stability is considered by means of a linear analysis within the frane of the force-free MHD. We argue that differentially rotating magnetospheres are unstable for any particular geometry of the magnetic field and rotation law. The characteristic growth time of instability is of the order of the rotation period. The instability can lead to fluctuations of the emission and enhancement of diffusion in the magnetosphere.

We review the evidence for and against the possibility that a strong enough poloidal field stabilizes an axisymmetric magnetostatic field configuration. We show that there does exist a class of resonant MHD waves which produce instability for any value of the ratio of poloidal and toroidal field strength. We argue that recent investigations of the stability of mixed poloidal and toroidal field configurations based on 3-d numerical simulations, can miss this instability because of the very large azimuthal wave numbers involved and its resonant character.

MHD waves can be responsible for plasma fluctuations and short-term variations of the pulsar emission. We consider the properties of plane and cylindrical waves that can exist in the force-free magnetosphere. Waves are considered by means of a linear analysis of the force-free MHD equations. We argue that these particular types of waves can exist in the magnetosphere of pulsars. These waves are closely related to the Alfven waves of the standard magnetohydrodynamics but are modified by the force-free condition and non-zero charge density. We derive the dispersion relation for magnetospheric waves and show that the wave periods are likely within the range $\sim 10^{-2}-10^(-4}$ s depending on the magnetospheric parameters.

Understanding the stability of the magnetic field in radiation zones is of crucial importance for various processes in stellar interior like mixing, circulation and angular momentum transport. The stability properties of a star containing a prominent toroidal field in a radiation zone is investigated by means of a linear stability analysis in the Boussinesq approximation taking into account the effect of thermal conductivity. The growth rate of the instability is explicitly calculated and the effects of stable stratification and heat transport are discussed in detail. It is argued that the stabilizing influence of gravity can never entirely suppress the instability caused by electric currents in radiation zones although the stable stratification can significantly decrease the growth rate of instability

We briefly review the turbulent mean-field dynamo action in protoneutron stars that are subject to convective and neutron finger instabilities during the early evolutionary phase. By solving the mean-field induction equation with the simplest model of $\alpha$-quenching we estimate the strength of the generated magnetic field. If the initial period of the protoneutron star is short, then the generated large-scale field is very strong ($> 3 \times 10^{13}$G) and exceeds the small-scale field at the neutron star surface, while if the rotation is moderate, then the pulsars are formed with more or less standard dipole fields ($< 3 \times 10^{13}$G) but with surface small-scale magnetic fields stronger than the dipole field. If rotation is very slow, then the mean-field dynamo does not operate, and the neutron star has no global field.

Magnetohydrodynamic instabilities can be responsible for the formation of structures with various scales in astrophysical jets. We consider the stability properties of jets containing both the azimuthal and axial field of subthermal strength. A magnetic field with complex topology in jets is suggested by theoretical models and is consistent with recent observations. Stability is discussed by means of a linear analysis of the ideal magnetohydrodynamic equations. We argue that in azimuthal and axial magnetic fields the jet is always unstable to non-axisymmetric perturbations. Stabilization does not occur even if the strengths of these field components are comparable. If the axial field is weaker than the azimuthal one, instability occurs for perturbations with any azimuthal wave number $m$, and the growth rate reaches a saturation value for low values of $m$. If the axial field is stronger than the toroidal one, the instability shows for perturbations with relatively high $m$.

A linear stability analysis of ionized discs with a temperature gradient and an external axial magnetic field is presented. It is shown that both hydromagnetic and thermomagnetic effects can lead to the amplification of waves and make discs unstable. The conditions under which the instabilities grow are found and the characteristic growth rate is calculated. The regimes at which both the thermomagnetic and magnetorotational instabilities can operate are discussed.

Magnetorotational instability (MRI) has been suggested to lead a rapid growth of the magnetic field in core collapse supernovae and produce departures from spherical syymmetry that can be important in determining the explosion mechanism. We address the problem of stability in differentially rotating magnetized proto-neutron stars at the beginning of their evolution. Criteria for MRI in proto-neutron stars are derived without simplying assumptions about a weak magnetic field and are substantially different from the standard condition. If the magnetic field is strong, MRI can occur only in the neighbourhood of the region where the spherical radial component of the magnetic field vanishes. The growth rate of MRI is relatively low except for perturbations with very small scales which usually are not detected in numerical simulations. We find that MRI in proto-neutron stars grows more slowly than than the double diffusive instability analogous the Goldreich-Schubert-Fricke instability in ordinary stars.

The MHD instabilities can generate complex field topologies even if the initial field configuration is a very simple one. We consider the stability properties of magnetic configurations containing a toroidal and an axial field. In this paper, we concentrate mainly on the behavior of non-axisymmetric perturbations in axisymmetric magnetic configurations. The stability is treated by a linear analysis of ideal MHD equations.In the presence of an axial field, it is shown that the instability can occur for a wide range of the azimuthal wavenumber $m$, and its growth rate increases with increasing $m$. At given $m$, the growth rate is at its maximum for perturbations with the axial wave-vector that makes the Alfv\'en frequency approximately vanishing. We argue that the instability of magnetic configurations in the ideal MHD can typically be dominated by perturbations with very short azimuthal and axial wavelengths.

Cool weakly ionized gaseous rotating disk form the basis for many models in astrophysics objects. Instabilities against perturbations in such disks play an important role in the theory of the formation of stars and planets. Traditionally, axisymmetric magnetohydrodynamic (MHD) and recently Hall-MHD instabilities have been thoroughly studied as providers of an efficient mechanism for radial transfer of angular momentum, and of density radial stratification. In the current work, the Hall instability against axisymmetric perturbations in incompressible rotating fluid in external poloidal and toroidal magnetic field is considered.

We study the efficiency of Joule heating in the crustal layers of young neutron stars. It is shown that dissipation of the magnetic field is highly inhomogeneous in the crust with much faster dissipation in relatively low density layers. In young neutron stars, the rate of Joule heating in the crust can exceed the standard luminosity of non-magnetic star and can even be comparable to the luminosity of magnetars. The results of calculations are compared with the available observational data. We argue that the crustal field model can well account for the data on the surface temperature and magnetic field of young neutron stars.

Observations of protostellar disks indicate the presence of the magnetic field of thermal (or superthermal) strength. In such a strong magnetic field, many MHD instabilities responsible for turbulent transport of the angular momentum are suppressed. We consider the shear-driven instability that can occur in protostellar disks even if the field is superthermal. This instability is caused by the combined influence of shear and compressibility in a magnetized gas and can be an efficient mechanism to generate turbulence in disks. The typical growth time is of the order of several rotation periods.

Stability properties of magnetic-field configurations containing the toroidal and axial field are considered. The stability is treated by making use of linear analysis. It is shown that the conditions required for the onset of instability are essentially different from those given by the necessary condition $d (s B_{\phi})/ds > 0$, where $s$ is the cylindrical radius. The growth rate of instability is calculated for a wide range of the parameters. We argue that the instability can operate in two different regimes depending on the strength of the axial field and the profile of the toroidal field.

We evaluate the generation of magnetosonic waves in differentially rotating magnetized plasma. Differential rotation leads to an increase of the azimuthal field by winding up the poloidal field lines into the toroidal field lines. An amplification of weak seed perturbations is considered in this time-dependent background state. It is shown that seed perturbations can be amplified by several orders of magnitude in a differentially rotating flow. The only necessary condition for this amplification is the presence of a non-vanishing component of the magnetic field in the direction of the angular velocity gradient.

We study the stability of a compressible magnetic plane Couette flow and show that compressibility profoundly alters the stability properties if the magnetic field has a component perpendicular to the direction of flow. The necessary condition of a newly found instability can be satisfied in a wide variety of flows in laboratory and astrophysical conditions. The instability can operate even in a very strong magnetic field which entirely suppresses other MHD instabilities. The growth time of this instability can be rather short and reach $\sim 10$ shear timescales.

We consider stability of differentially rotating non-magnetic proto-neutron stars. When neutrino transport is efficient, the star can be subject to a diffusive instability that can occur even in the convectively stable region. The instability arises on the time-scale comparable to the time-scale of thermal diffusion. Hydrodynamic motions driven by the instability can lead to anisotropy in the neutrino flux since the instability is suppressed near the equator and rotation axis.

We study the stability of compressible cylindrical differentially rotating flow in the presence of the magnetic field, and show that compressibility alters qualitatively the stability properties of flows. Apart from the well-known magnetorotational instability that can occur even in incompressible flow, there exist a new instability caused by compressibility. The necessary condition of the newly found instability can easily be satisfied in various flows in laboratory and astrophysical conditions and reads $B_{s} B_{\phi} \Omega' \neq 0$ where $B_{s}$ and $B_{\phi}$ are the radial and azimuthal magnetic fields, $\Omega' =d \Omega/ds$ with $s$ being the cylindrical radius. Contrary to the magnetorotational instability that occurs only if $\Omega$ decreases with $s$, the newly found instability operates at any sign of $\Omega'$. The considered instability can arise even in a very strong magnetic field that suppresses the magnetorotational instability.

We study the stability of a compressible differentially rotating flows in the presence of the magnetic field, and we show that the compressibility profoundly alters the previous results for a magnetized incompressible flow. The necessary condition of newly found instability can be easily satisfied in various flows in laboratory and astrophysical conditions and reads $B_{s} B_{\phi} \Omega' \neq 0$ where $B_{s}$ and $B_{\phi}$ are the radial and azimuthal components of the magnetic field, $\Omega' = d \Omega/ds$ with $s$ being the cylindrical radius. Contrary to the well-known magnetorotational instability that occurs only if $\Omega$ decreases with $s$, the instability considered in this paper may occur at any sign of $\Omega'$. The instability can operate even in a very strong magnetic field which entirely suppresses the standard magnetorotational instability. The growth time of instability can be as short as few rotation periods.

We consider dynamo action under the combined influence of turbulence and large-scale shear in sheared jets. Shear can stretch turbulent magnetic field lines in such a way that even turbulent motions showing mirror symmetry become suitable for generation of a large-scale magnetic field. We derive the integral induction equation governing the behaviour of the mean field in jets. The main result is that sheared jets may generate a large-scale magnetic field if shear is sufficiently strong. The generated mean field is mainly concentrated in a magnetic sheath surrounding the central region of a jet, and it exhibits sign reversals in the direction of the jet axis. Typically, the magnetic field in a sheath is dominated by the component along the jet that can reach equipartition with the kinetic energy of particles, The field in the central region of jets has a more disordered structure.

We discuss the mean-field dynamo action in protoneutron stars that are subject to instabilities during the early evolutionary phase. The mean field is generated in the neutron-finger unstable region where the Rossby number is $\sim 1$ and mean-field dynamo is efficient. Depending on the rotation rate, the mean-field dynamo can lead to the formation of three different types of pulsars. If the initial period of the protoneutron star is short, then the generated large-scale field is very strong ($> 3 \times 10^{13}$G) and exceeds the small-scale field at the neutron star surface. If rotation is moderate, then the pulsars are formed with more or less standard dipole fields ($< 3 \times 10^{13}$G) but with surface small-scale magnetic fields stronger than the dipole field. If rotation is very slow, then the mean-field dynamo does not operate, and the neutron star has no global field. Nevertheless, strong small-scale fields are generated in such pulsars, and they can manifest themselves as objects with very low spin-down rate but with a strong magnetic field inferred from the spectral features.

We study the secular instability of magnetized differentially rotating radiative zones taking account of viscosity and magnetic and thermal diffusivities. The considered instability generalizes the well-known Goldreich-Schubert-Fricke instability for the case of a sufficiently strong magnetic field. In magnetized stars, instability can lead to a formation of non-spherical unstable zones where weak turbulence mixes the material between the surface and interiors. Such unstable zones can manifest themselves by a non-spherical distribution of abundance anormalies on the stellar surface.

We study the relation between the mean effective surface temperature T_s and the internal temperature T_b for magnetic neutron stars, assuming that the magnetic field near the surface has a presumably small-scale structure. The heavy-element (iron) and light-element (accreted) heat-blanketing envelopes are considered, and the results are compared with the case of a dipole magnetic field. We argue that the difference in the T_b(T_s)-relation for different magnetic configurations is always much smaller than a possible difference caused by variations of the chemical composition in the envelope.

We consider the stability properties of the ocean of accreting magnetic neutron stars. It turns out that the ocean is always unstable due to the combined influence of the temperature and chemical composition gradients along the surface and of the Hall effect. Both the oscillatory and non-oscillatory modes can be unstable in accreting stars. The oscillatory instability grows on a short timescale of the order os 0.1-10 s depending on the lengthscale of a surface inhomogeneity and the wavelength of perturbations. The instability of non-oscillatory modes is typically much slower and can develop on a timescale of hours or days. Instability generates a weak turbulence that can be responsible for mixing between the surface and deep ocean layers and for spreading the accreted material over the stellar surface. Spectral features of heavy elements can be detected in the atmospheres of accreting stars due to mixing, and these features should be different in neutron stars with stable and unstable burning. Motions caused by instability can also be the reason for slow variations in the luminosity.

We have investigated the turbulent mean-field dynamo action in protoneutron stars that are subject to convective and neutron finger instabilities during the early evolutionary phase. While the first one develops mostly in the inner regions of the star, the second one is favored in the outer regions, where the Rossby number is much smaller and a mean-field dynamo action is more efficient. By solving the mean-field induction equation we have computed the critical spin period below which no dynamo action is possible and found it to be $\sim 1$ s for a wide range of stellar models and for both axisymmetric and non-axisymmetric magnetic fields. Because this critical period is substantially longer than the characteristic spin period of very young pulsars, we expect that a mean-field dynamo will be effective for most protoneutron stars. The saturation dipole field estimated by making use of the model of ``global'' quenching fits well the pulsar magnetic fields inferred from the spin-down data. Apart from the large scale magnetic field, our model predicts also a generation of small scale fields which are typically stronger than the poloidal field and can survive during the lifetime of pulsars. Extremely rapidly rotating protoneutron stars ($P \sim 1$ ms) may have the dipole field $\sim (3-6) \times 10^{14}$ G.

We study the conditions for convective instability in rotating, non-magnetic proto--neutron stars. The criteria that determine stability of nascent neutron stars are analogous to the Solberg--Hoiland conditions but including the presence of lepton gradients. Our results show that, for standard angular velocity profiles, convectively unstable modes with wave-vectors parallel to the rotation axis are suppressed by a stable angular momentum profile, while unstable modes with wave-vectors perpendicular to the axis remain unaltered. Since the wave-vector is perpendicular to the velocity perturbation, the directional selection of the unstable modes may result in fluid motions along the direction of the rotation axis. This occurs in rigidly rotating stars as well as in the inner core of differentially rotating stars. Our results provide a natural source of asymmetry for proto--neutron stars with the only requirement that angular velocities be of the order of the convective characteristic frequency.

We consider generation and evolution of small-scale magnetic fields in neutron stars. These fields can be generated by small-scale turbulent dynamo action soon after the collapse when the proto-neutron star is subject to convective and neutron finger instabilities. After instabilities stop, small-scale fields should be frozen into the crust that forms initially at high density about 10^{14} g/cm^{3} and then spreads to the surface. Because of high crustal conductivity, magnetic fields with the lengthscale about 1-3 km can survive in the crust as long as 10-100 Myr and form a sunspot-like structure at the surface of radiopulsars.

We have investigated the turbulent mean-field dynamo action in protoneutron stars that are subject to convective and neutron finger instabilities. While the first one develops mostly in the inner regions of the star, the second one is favoured in the outer regions, where the Rossby number is much smaller and a mean-field dynamo action is more efficient. By solving the mean-field induction equation we have computed the critical spin period below which no dynamo action is possible and found it to be ~1 s for a wide range of stellar models. Because this critical period is substantially longer than the characteristic spin period of very young pulsars, we expect that a mean-field dynamo will be effective for most protoneutron stars.

The stability properties of inviscid protostellar disks are examined taking into account the Hall effect and buoyancy. Depending on the parameters, different types of instabilities can exist in different regions of disks. In a very low ionized region, the instability associated with baroclinic effects of buoyancy is likely most efficient. The Hall-driven shear instability can lead to destabilization of regions with a higher ionization. The magnetorotational instability modified by buoyancy can only be a destabilizing factor in regions with strong magnetic field or a relatively high conductivity ($\sigma B^2 / \rho > a_e \Omega$, with $a_e$ the magnetization parameter of electrons).

We consider the transverse structure and stability properties of relativistic jets formed in the course of the collapse of a massive progenitor. Our numerical simulations show the presence of a strong shear in the bulk velocity of such jets. This shear can be responsible for a very rapid shear--driven instability that arises for any velocity profile. This conclusion has been confirmed both by numerical simulations and theoretical analysis. The instability leads to rapid fluctuations of the main hydrodynamical parameters (density, pressure, Lorentz factor, etc.). However, the perturbations of the density are effectively decoupled from those of the pressure because the beam of the jet is radiation--dominated. The characteristic growth time of instability is much shorter than the life time of the jet and, therefore, may lead to a complete turbulent beam. In the course of the non-linear evolution, these fluctuations may yield to internal shocks which can be randomly distributed in the jet. In the case that internal shocks in a ultrarelativistic outflow are responsible for the observed phenomenology of gamma-ray bursts, the proposed instability can well account for the short-term variability of gamma-ray light curves down to milliseconds.

The stability properties of newly born neutron stars, or proto--neutron stars, are considered. We take into account dissipative processes, such as neutrino transport and viscosity, in the presence of a magnetic field. In order to find the regions of the star subject to different sorts of instability, we derive the general instability criteria and apply it to evolutionary models of PNSs. The influence of the magnetic field on instabilities is analyzed and the critical magnetic field stabilizing the star is obtained. In the light of our results, we estimate of the maximum poloidal magnetic field that might be present in young pulsars or magnetars.

The linear hydrodynamic stability of proto-neutron stars (PNSs) is considered taking into account dissipative processes such as neutrino transport and viscosity. We obtain the general instability criteria which differ essentially from the well-known Ledoux criterion used in previous studies. We apply the criteria to evolutive models of PNSs that, in general, can be subject to the various known regimes such as neutron fingers and convective instabilities. Our results indicate that the fingers instability arises in a more extended region of the stellar volume and lasts a longer time than expected.

We consider Joule heating caused by dissipation of the magnetic field in the neutron star crust. This mechanism may be efficient in maintaining a relatively high surface temperature in very old neutron stars. Calculations of the thermal evolution show that, at the late evolutionary stage ($t \geq 10$ Myr), the luminosity of the neutron star is approximately equal to the energy released due to the field dissipation and is practically independent of the atmosphere models. At this stage, the surface temperature can be of the order of $3 \times 10^{4} - 10^{5}$K. Joule heating can maintain this high temperature during extremely long time ($\geq 100$ Myr), comparable with the decay time of the magnetic field.

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.

We consider the magnetic and spin evolution of isolated neutron stars assuming that the magnetic field is initially confined to the crust. The evolution of the crustal field is determined by the conductive properties of the crust which, in its turn, depend on the thermal history of the neutron star. Due to this fact, a study of the magnetic field decay may be a powerful diagnostic of the properties of matter in the core where the density is above the nuclear density. We treat the evolution of neutron stars for different possible equations of state and cooling scenarios (standard and, so called, accelerated cooling). The spin evolution is strongly influenced by the behaviour of the magnetic field. Assuming that the spin-down rate of the neutron star is determined by the magnetodipole radiation, we calculated the evolutionary tracks of isolated pulsars in the $B - \tau$ and $B-P$ planes, where $B$ and $P$ are the magnetic field and period, respectively, and $\tau$ is the spin-down age. The calculated tracks are compared with observational data on the magnetic field and period of pulsars. This comparison allows to infer the most suitable equations of state of nuclear matter and cooling model and to determine the range of parameters of the original magnetic configurations of pulsars.