Periodic changes in a thermal soft X-ray flux of a rotating neutron star indicate a non-uniform distribution of the surface temperature. A possible cause of this phenomenon is a suppression of the heat flux across the magnetic field lines in a crust and an envelope of magnetized neutron stars. In this paper we study three-dimensional effects, associated with non-axisymmetric magnetic fields in neutron stars. We calculate the surface temperature distribution by solving numerically a three dimensional heat transfer equation in a magnetized neutron star crust. We adopt an anisotropic (tensorial) electron thermal conductivity coefficient, which is derived as an analytical solution of the Boltzmann equation with a Chapman-Enskog method. To calculate the surface temperature distribution, we construct a local one-dimensional plane-parallel model ("Ts-Tb"-relationship) of a magnetized neutron star envelope. We then use it as an outer boundary condition for the three-dimensional problem in the crust to find the self-consistent solution. To study possible observational manifestations from anisotropic temperature distributions we calculate light curves with a composite black-body model. Our calculations show, that a non-axisymmetric magnetic field distribution can lead to the irregular non-sinusoidal shape of a pulse profile as well as in some cases a significant amplification of pulsations of the thermal flux in comparison to the pure-dipolar magnetic field configurations.

A solution to the Boltzmann equation is obtained for a magnetized plasma with strongly degenerate nonrelativistic electrons and nondegenerate nuclei. The components of the diffusion, thermal diffusion and diffusion thermoeffect tensors in a nonquantizing magnetic field are calculated in the Lorentz approximation without allowance for electron-electron collisions, which is asymptotically accurate for plasma with strongly degenerate electrons. Asymptotically accurate analytical expressions for the electron diffusion, thermal diffusion and diffusion thermoeffect tensors in the presence of a magnetic field are obtained for the first time. The expressions reveal a considerably more complicated dependence on magnetic field than analogous dependences derived in the previous publications on this subject.

The solution of Boltzmann equation for plasma in magnetic field, with arbitrarily degenerate electrons and non-degenerate nuclei, is obtained by Chapman-Enskog method. Functions, generalizing Sonin polynomials are used for obtaining an approximate solution. Fully ionized plasma is considered. The tensor of the heat conductivity coefficients in non-quantized magnetic field is calculated. For non-degenerate and strongly degenerate plasma the asymptotic analytic formulas are obtained, which are compared with results of previous authors. The Lorentz approximation, with neglecting of electron-electron encounters, is asymptotically exact for strongly degenerate plasma. We obtain, for the first time, in three polynomial approximation, with account of electron-electron collisions, analytical expressions for the heat conductivity tensor for non-degenerate electrons, in presence of a magnetic field. Account of the third polynomial improved substantially the precision of results. In two polynomial approximation our solution coincides with the published results. For strongly degenerate electrons we obtain, for the first time, an asymptotically exact analytical solution for the heat conductivity tensor in presence of a magnetic field. This solution has considerably more complicated dependence on the magnetic field than those in previous publications, and gives several times smaller relative value of a thermal conductivity across the magnetic field at $\omega\tau\gtrsim 0.8$.

The coefficients that determine the electron heat transfer and diffusion in the crust of neutron stars are calculated on the basis of a solution of the Boltzmann equation with allowance for degeneracy.